Map conformational landscapes of intrinsically disordered proteins with polymer physics quantities

Disordered proteins are conformationally flexible proteins that are biologically important and have been implicated in devastating diseases such as Alzheimer's disease and cancer. Unlike stably folded structured proteins, disordered proteins sample a range of different conformations that needs to be accounted for. Here, we treat disordered proteins as polymer chains, and compute a dimensionless quantity called instantaneous shape ratio (Rs), as Rs = Ree2/Rg2, where Ree is end-to-end distance and Rg is radius of gyration. Extended protein conformations tend to have high Ree compared with Rg, and thus have high Rs values, whereas compact conformations have smaller Rs values. We use a scatter plot of Rs (representing shape) against Rg (representing size) as a simple map of conformational landscapes. We first examine the conformational landscape of simple polymer models such as Random Walk, Self-Avoiding Walk, and Gaussian Walk (GW), and we notice that all protein/polymer maps lie within the boundaries of the GW map. We thus use the GW map as a reference and, to assess conformational diversity, we compute the fraction of the GW conformations (fC) covered by each protein/polymer. Disordered proteins all have high fC scores, consistent with their disordered nature. Each disordered protein accesses a different region of the reference map, revealing differences in their conformational ensembles. We additionally examine the conformational maps of the nonviral gene delivery vector polyethyleneimine at various protonation states, and find that they resemble disordered proteins, with coverage of the reference map decreasing with increasing protonation state, indicating decreasing conformational diversity. We propose that our method of combining Rs and Rg in a scatter plot generates a simple, meaningful map of the conformational landscape of a disordered protein, which in turn can be used to assess conformational diversity of disordered proteins.

Gene network based analysis identifies a coexpression module involved in regulating plasma lipids with high-fat diet response

Plasma lipids are modulated by gene variants and many environmental factors, including diet-associated weight gain. However, understanding how these factors jointly interact to influence molecular networks that regulate plasma lipid levels is limited. Here, we took advantage of the BXD recombinant inbred family of mice to query weight gain as an environmental stressor on plasma lipids. Coexpression networks were examined in both nonobese and obese livers, and a network was identified that specifically responded to the obesogenic diet. This obesity-associated module was significantly associated with plasma lipid levels and enriched with genes known to have functions related to inflammation and lipid homeostasis. We identified key drivers of the module, including Cidec, Cidea, Pparg, Cd36, and Apoa4. The Pparg emerged as a potential master regulator of the module as it can directly target 19 of the top 30 hub genes. Importantly, activation of this module is causally linked to lipid metabolism in humans, as illustrated by correlation analysis and inverse-variance weighed Mendelian randomization. Our findings provide novel insights into gene-by-environment interactions for plasma lipid metabolism that may ultimately contribute to new biomarkers, better diagnostics, and improved approaches to prevent or treat dyslipidemia in patients.

Exploring Structures and Dynamics of Protamine Molecules through Molecular Dynamics Simulations

Protamines are arginine-rich proteins that condense DNA in sperm. Despite their importance in reproduction, information on protamine structure is scarce. We, therefore, used molecular dynamics to examine the structures of salmon, bull P1, and human P1 protamines. The sizes and shapes of each protamine varied widely, indicating that they were disordered with structures covering a broad conformational landscape, from hairpin loop structures to extended coils. Despite their general disorder, the protamines did form secondary structures, including helices and hairpin loops. In eutherians, hairpins may promote disulfide bonding that facilitates protamine-DNA condensation, but the specifics of this bonding is not well established. We examined inter-residue distances in the simulations to predict residue pairs likely to form intramolecular bonds, leading to the identification of bonding pairs consistent with previous results in bull and human. These results support a model for eutherian protamine structures where a highly charged center is surrounded by disulfide-bond-stabilized loops.

Gene-by-environment modulation of lifespan and weight gain in the murine BXD family

How lifespan and body weight vary as a function of diet and genetic differences is not well understood. Here we quantify the impact of differences in diet on lifespan in a genetically diverse family of female mice, split into matched isogenic cohorts fed a low-fat chow diet (CD, n = 663) or a high-fat diet (HFD, n = 685). We further generate key metabolic data in a parallel cohort euthanized at four time points. HFD feeding shortens lifespan by 12%: equivalent to a decade in humans. Initial body weight and early weight gains account for longevity differences of roughly 4-6 days per gram. At 500 days, animals on a HFD typically gain four times as much weight as control, but variation in weight gain does not correlate with lifespan. Classic serum metabolites, often regarded as health biomarkers, are not necessarily strong predictors of longevity. Our data indicate that responses to a HFD are substantially modulated by gene-by-environment interactions, highlighting the importance of genetic variation in making accurate individualized dietary recommendations.

Functional Analysis of Genetic Variants and Somatic Mutations Impacting MicroRNA-Target Recognition: Bioinformatics Resources

MicroRNAs are small noncoding RNA molecules with great importance in regulating a large number of diverse biological processes in health and disease. MicroRNAs can bind to both coding and noncoding RNAs and regulate their stability and expression. Genetic variants and somatic mutations may alter microRNA sequences and their target sites and therefore impact microRNA-target recognition. Aberrant microRNA-target interactions have been associated with many diseases. In recent years, computational resources have been developed for retrieving, annotating, and analyzing the impact of mutations on microRNA-target recognition. In this chapter, we provide an overview on the computational analysis of mutations impacting microRNA target recognition, followed by a detailed tutorial on how to use three major Web-based bioinformatics resources: PolymiRTS ( http://compbio.uthsc.edu/miRSNP ), a database of genetic variants impacting microRNA target recognition; SomamiR ( http://compbio.uthsc.edu/SomamiR ), a database of somatic mutations affecting the interactions between microRNAs and their targets in mRNAs and noncoding RNAs; and miR2GO ( http://compbio.uthsc.edu/miR2GO ), a computational tool for knowledge-based functional analysis of genetic variants and somatic mutations in microRNA seed regions.

Structural Comparisons of PEI/DNA and PEI/siRNA Complexes Revealed with Molecular Dynamics Simulations

Polyplexes composed of polyethyleneimine (PEI) and DNA or siRNA have attracted great attention for their use in gene therapy. Although many physicochemical characteristics of these polyplexes remain unknown, PEI/DNA complexes have been repeatedly shown to be more stable than their PEI/siRNA counterparts. Here, we examine potential causes for this difference using atomistic molecular dynamics simulations of complexation between linear PEI and DNA or siRNA duplexes containing the same number of bases. The two types of polyplexes are stabilized by similar interactions, as PEI amines primarily interact with nucleic acid phosphate groups but also occasionally interact with groove atoms of both nucleic acids. However, the number of interactions in PEI/DNA complexes is greater than in comparable PEI/siRNA complexes, with interactions between protonated PEI amines and DNA being particularly enhanced. These results indicate that structural differences between DNA and siRNA may play a role in the increased stability of PEI/DNA complexes. In addition, we investigate the binding of PEI chains to polyplexes that have a net positive charge. The binding of PEI to these overcharged complexes involves interactions between PEI and areas on the nucleic acid surface that have maintained a negative electrostatic potential and is facilitated by the release of water from the nucleic acid.

Interactions of complex polymers with nanoporous substrate

With the advance of polymer synthesis, polymers that possess unique architectures such as stars or cyclic chains, and unique chemical composition distributions such as block copolymers or statistical copolymers have become frequently encountered. Characterization of these complex polymer systems drives the development of interactive chromatography where the adsorption of polymers on the porous substrate in chromatography columns is finely tuned. Liquid Chromatography at the Critical Condition (LCCC) in particular makes use of the existence of the Critical Adsorption Point (CAP) of polymers on solid surfaces and has been successfully applied to characterization of complex polymer systems. Interpretation and understanding of chromatography behaviour of complex polymers in interactive chromatography motivates theoretical/computational studies on the CAP of polymers and partitioning of these complex polymers near the CAP. This review article covers the theoretical questions encountered in chromatographic studies of complex polymers.

Photosensitizer-loaded gold nanorods for near infrared photodynamic and photothermal cancer therapy

Despite the advancement of photodynamic therapy and photothermal therapy, the ability to form compact nanocomplex for combined photodynamic and photothermal cancer therapy under a single near infrared irradiation remains limited. In this work, we prepared an integrated sub-100 nm nanosystem for simultaneous near infrared photodynamic and photothermal cancer therapy. The nanosystem was formed by adsorption of silicon 2,3-naphthalocyanine dihydroxide onto gold nanorod followed by covalent stabilization with alkylthiol linked polyethylene glycol. The effects of alkylthiol chain length on drug loading, release and cell killing efficacy were examined using 6-mercaptohexanoic acid, 11-mercaptoundecanoic acid and 16-mercaptohexadecanoic acid. We found that the loading efficiency of silicon 2,3-naphthalocyanine dihydroxide increased and the release rate decreased with the increase of the alkylthiol chain length. We demonstrated that the combined near infrared photodynamic and photothermal therapy using the silicon 2,3-naphthalocyanine dihydroxide-loaded gold nanorods exhibit superior efficacy in cancer cell destruction as compared to photodynamic therapy and photothermal therapy alone. The nanocomplex stabilized with 16-mercaptohexadecanoic acid linked polyethylene glycol provided highest cell killing efficiency as compared to those stabilized with the other two stabilizers under low drug dose. This new nanosystem has potential to completely eradicate tumors via noninvasive phototherapy, preventing tumor reoccurrence and metastasis.

Pre-transplant predictors of one yr weight gain after kidney transplantation

Clinically useful predictors of weight gain could be used to reduce the epidemic of post-kidney transplant obesity and resulting co-morbidities. The purpose of this study was to identify predictors of weight gain at 12 months following kidney transplant in a cohort of 96 recipients. Demographic, clinical, and environmental data were obtained at transplant and 12 months. Descriptive, correlational, and Bayesian network analysis were used to identify predictors. For the 52 (55.9%) recipients who gained weight, the average amount gained was 9.18 ± 6.59 kg. From the 15 baseline factors that met inclusion criteria, Bayesian network modeling identified four baseline predictors for weight gain: younger age, higher carbohydrate consumption, higher trunk fat percentage, and higher perception of mental health quality of life. Three are modifiable through either pre- or immediate post-transplant clinical intervention programs.

Interplay of Coil-Globule Transition and Surface Adsorption of a Lattice HP Protein Model

An end-grafted hydrophobic-polar (HP) model protein chain with alternating H and P monomers is studied to examine interactions between the critical adsorption transition due to surface attraction and the collapse transition due to pairwise attractive H–H interactions. We find that the critical adsorption phenomenon can always be observed; however, the critical adsorption temperature T_CAP is influenced by the attractive H–H interactions in some cases. When the collapse temperature T_c is lower than T_CAP, the critical adsorption of the HP chain is similar to that of a homopolymer without intrachain attractions and T_CAP remains unchanged, whereas the collapse transition is suppressed by the adsorption. In contrast, for cases where T_c is close to or higher than T_CAP, T_CAP of the HP chain is increased, indicating that a collapsed chain is more easily adsorbed on the surface. The strength of the H–H attraction also influences the statistical size and shape of the polymer, with strong H–H attractions resulting in adsorbed and collapsed chains adopting two-dimensional, circular conformations.

Comparison of monovalent and divalent ion distributions around a DNA duplex with molecular dynamics simulation and a Poisson-Boltzmann approach

The ion atmosphere created by monovalent (Na(+) ) or divalent (Mg(2+) ) cations surrounding a B-form DNA duplex were examined using atomistic molecular dynamics (MD) simulations and the nonlinear Poisson-Boltzmann (PB) equation. The ion distributions predicted by the two methods were compared using plots of radial and two-dimensional cation concentrations and by calculating the total number of cations and net solution charge surrounding the DNA. Na(+) ion distributions near the DNA were more diffuse in PB calculations than in corresponding MD simulations, with PB calculations predicting lower concentrations near DNA groove sites and phosphate groups and a higher concentration in the region between these locations. Other than this difference, the Na(+) distributions generated by the two methods largely agreed, as both predicted similar locations of high Na(+) concentration and nearly identical values of the number of cations and the net solution charge at all distances from the DNA. In contrast, there was greater disagreement between the two methods for Mg(2+) cation concentration profiles, as both the locations and magnitudes of peaks in Mg(2+) concentration were different. Despite experimental and simulation observations that Mg(2+) typically maintains its first solvation shell when interacting with nucleic acids, modeling Mg(2+) as an unsolvated ion during PB calculations improved the agreement of the Mg(2+) ion atmosphere predicted by the two methods and allowed for values of the number of bound ions and net solution charge surrounding the DNA from PB calculations that approached the values observed in MD simulations.

PolymiRTS Database 3.0: linking polymorphisms in microRNAs and their target sites with human diseases and biological pathways

Polymorphisms in microRNAs (miRNAs) and their target sites (PolymiRTS) are known to disrupt miRNA function, leading to the development of disease and variation in physiological and behavioral phenotypes. Here, we describe recent updates to the PolymiRTS database (http://compbio.uthsc.edu/miRSNP), an integrated platform for analyzing the functional impact of genetic polymorphisms in miRNA seed regions and miRNA target sites. Recent advances in genomic technologies have made it possible to identify miRNA–mRNA binding sites from direct mapping experiments such as CLASH (cross linking, ligation and sequencing of hybrids). We have integrated data from CLASH experiments in the PolymiRTS database to provide more complete and accurate miRNA–mRNA interactions. Other significant new features include (i) small insertions and deletions in miRNA seed regions and miRNA target sites, (ii) TargetScan context + score differences for assessing the impact of polymorphic miRNA–mRNA interactions and (iii) biological pathways. The browse and search pages of PolymiRTS allow users to explore the relations between the PolymiRTSs and gene expression traits, physiological and behavioral phenotypes, human diseases and biological pathways.

Bayesian Network Webserver: a comprehensive tool for biological network modeling

SUMMARY

The Bayesian Network Webserver (BNW) is a platform for comprehensive network modeling of systems genetics and other biological datasets. It allows users to quickly and seamlessly upload a dataset, learn the structure of the network model that best explains the data and use the model to understand relationships between network variables. Many datasets, including those used to create genetic network models, contain both discrete (e.g. genotype) and continuous (e.g. gene expression traits) variables, and BNW allows for modeling hybrid datasets. Users of BNW can incorporate prior knowledge during structure learning through an easy-to-use structural constraint interface. After structure learning, users are immediately presented with an interactive network model, which can be used to make testable hypotheses about network relationships.

AVAILABILITY AND IMPLEMENTATION

BNW, including a downloadable structure learning package, is available at http://compbio.uthsc.edu/BNW. (The BNW interface for adding structural constraints uses HTML5 features that are not supported by current version of Internet Explorer. We suggest using other browsers (e.g. Google Chrome or Mozilla Firefox) when accessing BNW).

CONTACT

ycui2@uthsc.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CTCFBSDB 2.0: a database for CTCF-binding sites and genome organization

CTCF is a highly conserved transcriptional regulator protein that performs diverse functions such as regulating gene expression and organizing the 3D structure of the genome. Here, we describe recent updates to a database of CTCF-binding sites, CTCFBSDB (http://insulatordb.uthsc.edu/), which now contains almost 15 million CTCF-binding sequences in 10 species. Since the original publication of the database, studies of the 3D structure of the genome, such as those provided by Hi-C experiments, have suggested that CTCF plays an important role in mediating intra- and inter-chromosomal interactions. To reflect this important progress, we have integrated CTCF-binding sites with genomic topological domains defined using Hi-C data. Additionally, the updated database includes new features enabled by new CTCF-binding site data, including binding site occupancy and the ability to visualize overlapping CTCF-binding sites determined in separate experiments.

SomamiR: a database for somatic mutations impacting microRNA function in cancer

Whole-genome sequencing of cancers has begun to identify thousands of somatic mutations that distinguish the genomes of normal tissues from cancers. While many germline mutations within microRNAs (miRNAs) and their targets have been shown to alter miRNA function in cancers and have been associated with cancer risk, the impact of somatic mutations on miRNA function has received relatively little attention. Here, we have created the SomamiR database (http://compbio.uthsc.edu/SomamiR/) to provide a comprehensive resource that integrates several types of data for use in investigating the impact of somatic and germline mutations on miRNA function in cancer. The database contains somatic mutations that may create or disrupt miRNA target sites and integrates these somatic mutations with germline mutations within the same target sites, genome-wide and candidate gene association studies of cancer and functional annotations that link genes containing mutations with cancer. Additionally, the database contains a collection of germline and somatic mutations in miRNAs and their targets that have been experimentally shown to impact miRNA function and have been associated with cancer.

Integrative analysis of somatic mutations altering microRNA targeting in cancer genomes

Determining the functional impact of somatic mutations is crucial to understanding tumorigenesis and metastasis. Recent sequences of several cancers have provided comprehensive lists of somatic mutations across entire genomes, enabling investigation of the functional impact of somatic mutations in non-coding regions. Here, we study somatic mutations in 3'UTRs of genes that have been identified in four cancers and computationally predict how they may alter miRNA targeting, potentially resulting in dysregulation of the expression of the genes harboring these mutations. We find that somatic mutations create or disrupt putative miRNA target sites in the 3'UTRs of many genes, including several genes, such as MITF, EPHA3, TAL1, SCG3, and GSDMA, which have been previously associated with cancer. We also integrate the somatic mutations with germline mutations and results of association studies. Specifically, we identify putative miRNA target sites in the 3'UTRs of BMPR1B, KLK3, and SPRY4 that are disrupted by both somatic and germline mutations and, also, are in linkage disequilibrium blocks with high scoring markers from cancer association studies. The somatic mutation in BMPR1B is located in a target site of miR-125b; germline mutations in this target site have previously been both shown to disrupt regulation of BMPR1B by miR-125b and linked with cancer.

Systematic analysis of microRNA targeting impacted by small insertions and deletions in human genome

MicroRNAs (miRNAs) are small noncoding RNA that play an important role in posttranscriptional regulation of mRNA. Genetic variations in miRNAs or their target sites have been shown to alter miRNA function and have been associated with risk for several diseases. Previous studies have focused on the most abundant type of genetic variations, single nucleotide polymorphisms (SNPs) that affect miRNA-mRNA interactions. Here, we systematically identified small insertions and deletions (indels) in miRNAs and their target sites, and investigated the effects of indels on miRNA targeting. We studied the distribution of indels in miRNAs and their target sites and found that indels in mature miRNAs, experimentally supported miRNA target sites and PAR-CLIP footprints have significantly lower density compared to flanking regions. We identified over 20 indels in the seed regions of miRNAs, which may disrupt the interactions between these miRNAs and their target genes. We also identified hundreds of indels that alter experimentally supported miRNA target sites. We mapped these genes to human disease pathways to identify indels that affect miRNA targeting in these pathways. We also used the results of genome-wide association studies (GWAS) to identify potential links between miRNA-related indels and diseases.

Treatment- and population-dependent activity patterns of behavioral and expression QTLs

Genetic control of gene expression and higher-order phenotypes is almost invariably dependent on environment and experimental conditions. We use two families of recombinant inbred strains of mice (LXS and BXD) to study treatment- and genotype-dependent control of hippocampal gene expression and behavioral phenotypes. We analyzed responses to all combinations of two experimental perturbations, ethanol and restraint stress, in both families, allowing for comparisons across 8 combinations of treatment and population. We introduce the concept of QTL activity patterns to characterize how associations between genomic loci and traits vary across treatments. We identified several significant behavioral QTLs and many expression QTLs (eQTLs). The behavioral QTLs are highly dependent on treatment and population. We classified eQTLs into three groups: cis-eQTLs (expression variation that maps to within 5 Mb of the cognate gene), syntenic trans-eQTLs (the gene and the QTL are on the same chromosome but not within 5 Mb), and non-syntenic trans-eQTLs (the gene and the QTL are on different chromosomes). We found that most non-syntenic trans-eQTLs were treatment-specific whereas both classes of syntenic eQTLs were more conserved across treatments. We also found there was a correlation between regions along the genome enriched for eQTLs and SNPs that were conserved across the LXS and BXD families. Genes with eQTLs that co-localized with the behavioral QTLs and displayed similar QTL activity patterns were identified as potential candidate genes associated with the phenotypes, yielding identification of novel genes as well as genes that have been previously associated with responses to ethanol.

PolymiRTS Database 2.0: linking polymorphisms in microRNA target sites with human diseases and complex traits

The polymorphism in microRNA target site (PolymiRTS) database aims to identify single-nucleotide polymorphisms (SNPs) that affect miRNA targeting in human and mouse. These polymorphisms can disrupt the regulation of gene expression by miRNAs and are candidate genetic variants responsible for transcriptional and phenotypic variation. The database is therefore organized to provide links between SNPs in miRNA target sites, cis-acting expression quantitative trait loci (eQTLs), and the results of genome-wide association studies (GWAS) of human diseases. Here, we describe new features that have been integrated in the PolymiRTS database, including: (i) polymiRTSs in genes associated with human diseases and traits in GWAS, (ii) polymorphisms in target sites that have been supported by a variety of experimental methods and (iii) polymorphisms in miRNA seed regions. A large number of newly identified microRNAs and SNPs, recently published mouse phenotypes, and human and mouse eQTLs have also been integrated into the database. The PolymiRTS database is available at http://compbio.uthsc.edu/miRSNP/.

Understanding the protonation behavior of linear polyethylenimine in solutions through Monte Carlo simulations

The success of polyethyleneimine (PEI) as a nonviral-based gene delivery vector has been attributed to its proton buffering capacity. Despite the great interest in PEI for its use in nonviral-based gene delivery, the protonation behavior of PEI in solution is not well understood. Earlier experimental studies have reported inconsistent values of the protonation state of PEI. In this work, we report our investigation of the protonation behavior of a realistic linear PEI (lPEI) with computational approaches. Reported experimental pK(a) values of several diamine compounds are first examined. A screened Coulombic interaction with a distance dependence dielectric is shown to reproduce the shifted pK(a) values of the model diamine compounds. Then atomistic molecular dynamic simulations of lPEI chain with 20 repeating units are performed and the results are used to provide parameters for a coarse-grained polyamine model. The screened Coulombic interaction is then incorporated in the coarse-grained lPEI chain and computational titrations are performed. The obtained computational titration curves of lPEI in solutions were found to be in best agreement with experimental results by Smits et al., but the computational titration curves have too strong of a dependence on salt concentration compared to the experimental results by Smits et al. Disregarding the discrepancy in the salt dependence, our computational titrations reveal that approximately 55% of the lPEI amine groups are protonated under physiological conditions in solution with a nearly alternating arrangement of protonated and nonprotonated amines. Titrations of lPEI in the presence of a polyanion are also performed to determine how the charge state of lPEI could be affected by complexation with DNA in gene therapy preparations. While the presence of the polyanion increases the degree of protonation of the PEI, many of PEI amines remain unprotonated under physiological conditions, providing evidence that PEI complexed with DNA could still have proton buffering capacity. Potential sources of error that have resulted in the inconsistency of previously reported protonation states of PEI were also discussed.