ALBERT

Description: Motivation: Loops in proteins are often involved in biochemical functions. Their irregularity and flexibility make experimental structure determination and computational modeling challenging. Most current loop modeling methods focus on modeling single loops. In protein structure prediction, multiple loops often need to be modeled simultaneously. As interactions among loops in spatial proximity can be rather complex, sampling the conformations of multiple interacting loops is a challenging task. Results: In this study, we report a new method called m ulti-loop Di stance-guided S equential chain- Gro wth Monte Carlo ( M -D i SG ro ) for prediction of the conformations of multiple interacting loops in proteins. Our method achieves an average RMSD of 1.93 Å for lowest energy conformations of 36 pairs of interacting protein loops with the total length ranging from 12 to 24 residues. We further constructed a data set containing proteins with 2, 3 and 4 interacting loops. For the most challenging target proteins with four loops, the average RMSD of the lowest energy conformations is 2.35 Å. Our method is also tested for predicting multiple loops in β-barrel membrane proteins. For outer-membrane protein G, the lowest energy conformation has a RMSD of 2.62 Å for the three extracellular interacting loops with a total length of 34 residues (12, 12 and 10 residues in each loop). Availability and implementation : The software is freely available at: tanto.bioe.uic.edu/m-DiSGro. Contact: jinfeng@stat.fsu.edu or jliang@uic.edu Supplementary information : Supplementary data are available at Bioinformatics online.

Description: Viruses with single-stranded (ss) DNA genomes infect hosts in all three domains of life and include many medically, ecologically, and economically important pathogens. Recently, a new group of ssDNA viruses with chimeric genomes has been discovered through viral metagenomics. These chimeric viruses combine capsid protein genes and replicative protein genes that, respectively, appear to have been inherited from viruses with positive-strand RNA genomes, such as tombusviruses, and ssDNA genomes, such as circoviruses, nanoviruses or geminiviruses. Here, we describe the genome sequence of a new representative of this virus group and reveal an additional layer of chimerism among ssDNA viruses. We show that not only do these viruses encompass genes for capsid proteins and replicative proteins that have distinct evolutionary histories, but also the replicative genes themselves are chimeras of functional domains inherited from viruses of different families. Our results underscore the importance of horizontal gene transfer in the evolution of ssDNA viruses and the role of genetic recombination in the emergence of novel virus groups.

Description: A number of studies have been conducted that link mitochondrial dysfunction (MD) to Huntington's disease (HD); however, contradicting results had resulted in a lack of a clear mechanism that links expression of mutant Huntingtin protein and MD. Mouse homozygous (HM) and heterozygous (HT) mutant striatal cells with two or one allele encoding for a mutant huntingtin protein with 111 polyGln repeats showed a significant impairment of the mitochondrial disulfide relay system (MDRS). This system (consisting of two proteins, Gfer and Mia40) is involved in the mitochondrial import of Cys-rich proteins. The Gfer-to-Mia40 ratio was significantly altered in HM cells compared with controls, along with the expression of mitochondrial proteins considered substrates of the MDRS. In progenitors and differentiated neuron-like HM cells, impairment of MDRS were accompanied by deficient oxidative phosphorylation, Complex I, IV and V activities, decreased mtDNA copy number and transcripts, accumulation of mtDNA deletions and changes in mitochondrial morphology, consistent with other MDRS-deficient biological models, thus providing a framework for the energy deficits observed in this HD model. The majority (〉90%) of the mitochondrial outcomes exhibited a gene–dose dependency with the expression of mutant Htt. Finally, decreases in the mtDNA copy number, along with the accumulation of mtDNA deletions, provide a mechanism for the progressive neurodegeneration observed in HD patients.

Description: Prader–Willi syndrome (PWS), a genetic disorder of obesity, intellectual disability and sleep abnormalities, is caused by loss of non-coding RNAs on paternal chromosome 15q11-q13. The imprinted minimal PWS locus encompasses a long non-coding RNA (lncRNA) transcript processed into multiple SNORD116 small nucleolar RNAs and the spliced exons of the host gene, 116HG . However, both the molecular function and the disease relevance of the spliced lncRNA 116HG are unknown. Here, we show that 116HG forms a subnuclear RNA cloud that co-purifies with the transcriptional activator RBBP5 and active metabolic genes, remains tethered to the site of its transcription and increases in size in post-natal neurons and during sleep. Snord116del mice lacking 116HG exhibited increased energy expenditure corresponding to the dysregulation of diurnally expressed Mtor and circadian genes Clock , Cry1 and Per2 . These combined genomic and metabolic analyses demonstrate that 116HG regulates the diurnal energy expenditure of the brain. These novel molecular insights into the energy imbalance in PWS should lead to improved therapies and understanding of lncRNA roles in complex neurodevelopmental and metabolic disorders.

Description: Systemic lupus erythematosus (SLE) is a prototypic autoimmune disease that affects mainly females. What role the X chromosome plays in the disease has always been an intriguing question. In this study, we examined the genetic variants on the X chromosome through meta-analysis of two genome-wide association studies (GWAS) on SLE on Chinese Han populations. Prominent association signals from the meta-analysis were replicated in 4 additional Asian cohorts, with a total of 5373 cases and 9166 matched controls. We identified a novel variant in PRPS2 on Xp22.3 as associated with SLE with genome-wide significance (rs7062536, OR = 0.84, P = 1.00E–08). Association of the L1CAM-MECP2 region with SLE was reported previously. In this study, we identified independent contributors in this region in NAA10 (rs2071128, OR = 0.81, P = 2.19E–13) and TMEM187 (rs17422, OR = 0.75, P = 1.47E–15), in addition to replicating the association from IRAK1-MECP2 region (rs1059702, OR = 0.71, P = 2.40E–18) in Asian cohorts. The X-linked susceptibility variants showed higher effect size in males than that in females, similar to results from a genome-wide survey of associated SNPs on the autosomes. These results suggest that susceptibility genes identified on the X chromosome, while contributing to disease predisposition, might not contribute significantly to the female predominance of this prototype autoimmune disease.

Description: Oligodontia, which is the congenital absence of six or more permanent teeth, excluding the third molars, may contribute to masticatory dysfunction, speech alteration, aesthetic problems and malocclusion. Msh homeobox 1 (MSX1) was the first gene identified as causing non-syndromic oligodontia. In this study, we identified a novel heterozygous non-stop mutation (c.910_911dupTA, p.*304Tyrext*48) in MSX1 in a Chinese family with autosomal dominant non-syndromic oligodontia. This novel mutation substitutes the stop codon with a tyrosine residue, potentially adding 48 amino acids to the C-terminus of MSX1. Further in vitro study found that mutant MSX1 could be expressed but had lost its ability to enter the nucleus. This is the first report indicating that a non-stop mutation in MSX1 is responsible for oligodontia. This study broadens the mutation spectrum for MSX1 and provides a new way to clarify the mechanism of MSX1 in tooth agenesis.

Description: : Systematic studies of drug repositioning require the integration of multi-level drug data, including basic chemical information (such as SMILES), drug targets, target-related signaling pathways, clinical trial information and Food and Drug Administration (FDA)-approval information, to predict new potential indications of existing drugs. Currently available databases, however, lack query support for multi-level drug information and thus are not designed to support drug repositioning studies. DrugMap Central (DMC), an online tool, is developed to help fill the gap. DMC enables the users to integrate, query, visualize, interrogate, and download multi-level data of known drugs or compounds quickly for drug repositioning studies all within one system. Availability : DMC is accessible at http://r2d2drug.org/DMC.aspx . Contact : STWong@tmhs.org

Description: We aimed to understand the relation of photosynthetic rate ( A ) with g s and electron transport rate (ETR) in species of great taxonomic range and light adaptation capability during photosynthetic light induction. We studied three woody species ( Alnus formosana , Ardisia crenata and Ardisia cornudentata ) and four fern species ( Pyrrosia lingus , Asplenium antiquum , Diplazium donianum and Archangiopteris somai ) with different light adaptation capabilities. Pot-grown materials received 100 and/or 10% sunlight according to their light adaptation capabilities. At least 4 months after light acclimation, CO 2 and H 2 O exchange and chlorophyll fluorescence were measured simultaneously by equipment in the laboratory. In plants adapted or acclimated to low light, dark-adapted leaves exposed to 500 or 2000 µmol m –2  s –1 photosynthetic photon flux (PPF) for 30 min showed low gross photosynthetic rate ( P g ) and short time required to reach 90% of maximum P g (). At the initiation of illumination, two broad-leaved understory shrubs and the four ferns, especially ferns adapted to heavy shade, showed higher stomatal conductance ( g s ) than pioneer tree species; materials with higher g s had short at both 500 and 2000 µmol m –2  s –1 PPF. With 500 or 2000 µmol m –2  s –1 PPF, the g s for the three woody species increased from 2 to 30 min after the start of illumination, but little change in the g s of the four ferns. Thus, P g and g s were not correlated for all material measured at the same PPF and induction time. However, P g was positively correlated with ETR, even though CO 2 assimilation may be influenced by stomatal, biochemical and photoinhibitory limitations. In addition, was closely related to time required to reach 90% maximal ETR for all materials and with two levels of PPF combined. Thus, ETR is a good indicator for estimating the light induction of photosynthetic rate of species, across a wide taxonomic range and light adaptation and acclimation capability.

Description: Motivation: Currently there are no curative anticancer drugs, and drug resistance is often acquired after drug treatment. One of the reasons is that cancers are complex diseases, regulated by multiple signaling pathways and cross talks among the pathways. It is expected that drug combinations can reduce drug resistance and improve patients’ outcomes. In clinical practice, the ideal and feasible drug combinations are combinations of existing Food and Drug Administration-approved drugs or bioactive compounds that are already used on patients or have entered clinical trials and passed safety tests. These drug combinations could directly be used on patients with less concern of toxic effects. However, there is so far no effective computational approach to search effective drug combinations from the enormous number of possibilities. Results: In this study, we propose a novel systematic computational tool DrugComboRanker to prioritize synergistic drug combinations and uncover their mechanisms of action. We first build a drug functional network based on their genomic profiles, and partition the network into numerous drug network communities by using a Bayesian non-negative matrix factorization approach. As drugs within overlapping community share common mechanisms of action, we next uncover potential targets of drugs by applying a recommendation system on drug communities. We meanwhile build disease-specific signaling networks based on patients’ genomic profiles and interactome data. We then identify drug combinations by searching drugs whose targets are enriched in the complementary signaling modules of the disease signaling network. The novel method was evaluated on lung adenocarcinoma and endocrine receptor positive breast cancer, and compared with other drug combination approaches. These case studies discovered a set of effective drug combinations top ranked in our prediction list, and mapped the drug targets on the disease signaling network to highlight the mechanisms of action of the drug combinations. Availability and implementation: The program is available on request. Contact: stwong@tmhs.org

Description: Systemic lupus erythematosus (SLE) has a complex etiology and is affected by both genetic and environmental factors. Although more than 40 loci have shown robust association with SLE, the details of these loci, such as the independent contributors and the genes involved, are still unclear. In this study, we performed meta-analysis of two existing genome-wide association studies (GWASs) on Chinese Han populations from Hong Kong and Anhui, China, and followed the findings by further replication on three additional Chinese and Thailand cohorts with a total of 4254 cases and 6262 controls matched geographically and ethnically. We discovered multiple susceptibility variants for SLE in the 11q23.3 region, including variants in/near PHLDB1 (rs11603023, P_ combined = 1.25E–08, OR = 1.20), DDX6 (rs638893, P_ combined = 5.19E–07, OR = 1.22) and CXCR5 (rs10892301, P_ combined = 2.51E–08, OR = 0.85). Genetic contributions from the newly identified variants were all independent of SNP rs4639966, whose association was reported from the previous GWAS. In addition, the three newly identified variants all showed independent association with the disease through modeling by both stepwise and conditional logistic regression. The presence of multiple independent variants in this region emphasizes its role in SLE susceptibility, and also hints the possibility that distinct biological mechanisms might be involved in the disease involving this genomic region.