Portosystemic shunt placement reveals blood signatures for the development of hepatic encephalopathy through mass spectrometry.

Elective transjugular intrahepatic portosystemic shunt (TIPS) placement can worsen cognitive dysfunction in hepatic encephalopathy (HE) patients due to toxins, including possible microbial metabolites, entering the systemic circulation. We conducted untargeted metabolomics on a prospective cohort of 22 patients with cirrhosis undergoing elective TIPS placement and followed them up to one year post TIPS for HE development. Here we suggest that pre-existing intrahepatic shunting predicts HE severity post-TIPS. Bile acid levels decrease in the peripheral vein post-TIPS, and the abundances of three specific conjugated di- and tri-hydroxylated bile acids are inversely correlated with HE grade. Bilirubins and glycerophosphocholines undergo chemical modifications pre- to post-TIPS and based on HE grade. Our results suggest that TIPS-induced metabolome changes can impact HE development, and that pre-existing intrahepatic shunting could be used to predict HE severity post-TIPS.

DNA binding specificity of all four Saccharomyces cerevisiae forkhead transcription factors.

Quantifying the nucleotide preferences of DNA binding proteins is essential to understanding how transcription factors (TFs) interact with their targets in the genome. High-throughput in vitro binding assays have been used to identify the inherent DNA binding preferences of TFs in a controlled environment isolated from confounding factors such as genome accessibility, DNA methylation, and TF binding cooperativity. Unfortunately, many of the most common approaches for measuring binding preferences are not sensitive enough for the study of moderate-to-low affinity binding sites, and are unable to detect small-scale differences between closely related homologs. The Forkhead box (FOX) family of TFs is known to play a crucial role in regulating a variety of key processes from proliferation and development to tumor suppression and aging. By using the high-sequencing depth SELEX-seq approach to study all four FOX homologs in Saccharomyces cerevisiae, we have been able to precisely quantify the contribution and importance of nucleotide positions all along an extended binding site. Essential to this process was the alignment of our SELEX-seq reads to a set of candidate core sequences determined using a recently developed tool for the alignment of enriched k-mers and a newly developed approach for the reprioritization of candidate cores.

Pre- and Post-Portosystemic Shunt Placement Metabolomics Reveal Molecular Signatures for the Development of Hepatic Encephalopathy.

Hepatic encephalopathy (HE) is a common complication of advanced liver disease causing brain dysfunction. This is likely due to the accumulation of unfiltered toxins within the bloodstream. A known risk factor for developing or worsening HE is the placement of a transjugular intrahepatic portosystemic shunt (TIPS), which connects the pre-hepatic and post-hepatic circulation allowing some blood to bypass the dysfunctional liver and decreases portal hypertension. To better understand the pathophysiology of post-TIPS HE, we conducted a multi-center prospective cohort study employing metabolomic analyses on hepatic vein and peripheral vein blood samples from participants with cirrhosis undergoing elective TIPS placement, measuring chemical modifications and changes in concentrations of metabolites resulting from TIPS placement. In doing so, we identified numerous alterations in metabolites, including bile acids, glycerophosphocholines, and bilirubins possibly implicated in the development and severity of HE.

Diet and feeding pattern modulate diurnal dynamics of the ileal microbiome and transcriptome.

Compositional oscillations of the gut microbiome are essential for normal peripheral circadian rhythms, both of which are disrupted in diet-induced obesity (DIO). Although time-restricted feeding (TRF) maintains circadian synchrony and protects against DIO, its impact on the dynamics of the cecal gut microbiome is modest. Thus, other regions of the gut, particularly the ileum, the nexus for incretin and bile acid signaling, may play an important role in entraining peripheral circadian rhythms. We demonstrate the effect of diet and feeding rhythms on the ileal microbiome composition and transcriptome in mice. The dynamic rhythms of ileal microbiome composition and transcriptome are dampened in DIO. TRF partially restores diurnal rhythms of the ileal microbiome and transcriptome, increases GLP-1 release, and alters the ileal bile acid pool and farnesoid X receptor (FXR) signaling, which could explain how TRF exerts its metabolic benefits. Finally, we provide a web resource for exploration of ileal microbiome and transcriptome circadian data.

Intermittent Hypoxia and Hypercapnia Alter Diurnal Rhythms of Luminal Gut Microbiome and Metabolome.

Obstructive sleep apnea (OSA), characterized by intermittent hypoxia and hypercapnia (IHC), affects the composition of the gut microbiome and metabolome. The gut microbiome has diurnal oscillations that play a crucial role in regulating circadian and overall metabolic homeostasis. Thus, we hypothesized that IHC adversely alters the gut luminal dynamics of key microbial families and metabolites. The objective of this study was to determine the diurnal dynamics of the fecal microbiome and metabolome of Apoe-/- mice after a week of IHC exposure. Individually housed, 10-week-old Apoe-/- mice on an atherogenic diet were split into two groups. One group was exposed to daily IHC conditions for 10 h (Zeitgeber time 2 [ZT2] to ZT12), while the other was maintained in room air. Six days after the initiation of the IHC conditions, fecal samples were collected every 4 h for 24 h (6 time points). We performed 16S rRNA gene amplicon sequencing and untargeted liquid chromatography-mass spectrometry (LC-MS) to assess changes in the microbiome and metabolome. IHC induced global changes in the cyclical dynamics of the gut microbiome and metabolome. Ruminococcaceae, Lachnospiraceae, S24-7, and Verrucomicrobiaceae had the greatest shifts in their diurnal oscillations. In the metabolome, bile acids, glycerolipids (phosphocholines and phosphoethanolamines), and acylcarnitines were greatly affected. Multi-omic analysis of these results demonstrated that Ruminococcaceae and tauro-ß-muricholic acid (TßMCA) cooccur and are associated with IHC conditions and that Coriobacteriaceae and chenodeoxycholic acid (CDCA) cooccur and are associated with control conditions. IHC significantly change the diurnal dynamics of the fecal microbiome and metabolome, increasing members and metabolites that are proinflammatory and proatherogenic while decreasing protective ones. IMPORTANCE People with obstructive sleep apnea are at a higher risk of high blood pressure, type 2 diabetes, cardiac arrhythmias, stroke, and sudden cardiac death. We wanted to understand whether the gut microbiome changes induced by obstructive sleep apnea could potentially explain some of these medical problems. By collecting stool from a mouse model of this disease at multiple time points during the day, we studied how obstructive sleep apnea changed the day-night patterns of microbes and metabolites of the gut. Since the oscillations of the gut microbiome play a crucial role in regulating metabolism, changes in these oscillations can explain why these patients can develop so many metabolic problems. We found changes in microbial families and metabolites that regulate many metabolic pathways contributing to the increased risk for heart disease seen in patients with obstructive sleep apnea.

The RNA Polymerase α Subunit Recognizes the DNA Shape of the Upstream Promoter Element.

We demonstrate here that the α subunit C-terminal domain of Escherichia coli RNA polymerase (αCTD) recognizes the upstream promoter (UP) DNA element via its characteristic minor groove shape and electrostatic potential. In two compositionally distinct crystallized assemblies, a pair of αCTD subunits bind in tandem to the UP element consensus A-tract that is 6 bp in length (A6-tract), each with their arginine 265 guanidinium group inserted into the minor groove. The A6-tract minor groove is significantly narrowed in these crystal structures, as well as in computationally predicted structures of free and bound DNA duplexes derived by Monte Carlo and molecular dynamics simulations, respectively. The negative electrostatic potential of free A6-tract DNA is substantially enhanced compared to that of generic DNA. Shortening the A-tract by 1 bp is shown to "knock out" binding of the second αCTD through widening of the minor groove. Furthermore, in computationally derived structures with arginine 265 mutated to alanine in either αCTD, either with or without the "knockout" DNA mutation, contact with the DNA is perturbed, highlighting the importance of arginine 265 in achieving αCTD-DNA binding. These results demonstrate that the importance of the DNA shape in sequence-dependent recognition of DNA by RNA polymerase is comparable to that of certain transcription factors.

Landscape of DNA binding signatures of myocyte enhancer factor-2B reveals a unique interplay of base and shape readout.

Myocyte enhancer factor-2B (MEF2B) has the unique capability of binding to its DNA target sites with a degenerate motif, while still functioning as a gene-specific transcriptional regulator. Identifying its DNA targets is crucial given regulatory roles exerted by members of the MEF2 family and MEF2B's involvement in B-cell lymphoma. Analyzing structural data and SELEX-seq experimental results, we deduced the DNA sequence and shape determinants of MEF2B target sites on a high-throughput basis in vitro for wild-type and mutant proteins. Quantitative modeling of MEF2B binding affinities and computational simulations exposed the DNA readout mechanisms of MEF2B. The resulting binding signature of MEF2B revealed distinct intricacies of DNA recognition compared to other transcription factors. MEF2B uses base readout at its half-sites combined with shape readout at the center of its degenerate motif, where A-tract polarity dictates nuances of binding. The predominant role of shape readout at the center of the core motif, with most contacts formed in the minor groove, differs from previously observed protein-DNA readout modes. MEF2B, therefore, represents a unique protein for studies of the role of DNA shape in achieving binding specificity. MEF2B-DNA recognition mechanisms are likely representative for other members of the MEF2 family.

Crystal Structures of Ternary Complexes of MEF2 and NKX2-5 Bound to DNA Reveal a Disease Related Protein-Protein Interaction Interface.

MEF2 and NKX2-5 transcription factors interact with each other in cardiogenesis and are necessary for normal heart formation. Despite evidence suggesting that these two transcription factors function synergistically and possibly through direct physical interactions, molecular mechanisms by which they interact are not clear. Here we determined the crystal structures of ternary complexes of MEF2 and NKX2-5 bound to myocardin enhancer DNA in two crystal forms. These crystal structures are the first example of human MADS-box/homeobox ternary complex structures involved in cardiogenesis. Our structures reveal two possible modes of interactions between MEF2 and NKX2-5: MEF2 and NKX bind to adjacent DNA sites to recognize DNA in cis; and MEF2 and NKX bind to different DNA strands to interact with each other in trans via a conserved protein-protein interface observed in both crystal forms. Disease-related mutations are mapped to the observed protein-protein interface. Our structural studies provide a starting point to understand and further study the molecular mechanisms of the interactions between MEF2 and NKX2.5 and their roles in cardiogenesis.

Structure of the Forkhead Domain of FOXA2 Bound to a Complete DNA Consensus Site.

FOXA2, a member of the forkhead family of transcription factors, plays essential roles in liver development and bile acid homeostasis. In this study, we report a 2.8 Å co-crystal structure of the FOXA2 DNA-binding domain (FOXA2-DBD) bound to a DNA duplex containing a forkhead consensus binding site (GTAAACA). The FOXA2-DBD adopts the canonical winged-helix fold, with helix H3 and wing 1 regions mainly mediating the DNA recognition. Although the wing 2 region was not defined in the structure, isothermal titration calorimetry assays suggested that this region was required for optimal DNA binding. Structure comparison with the FOXA3-DBD bound to DNA revealed more major groove contacts and fewer minor groove contacts in the FOXA2 structure than in the FOXA3 structure. Structure comparison with the FOXO1-DBD bound to DNA showed that different forkhead proteins could induce different DNA conformations upon binding to identical DNA sequences. Our findings provide the structural basis for FOXA2 protein binding to a consensus forkhead site and elucidate how members of the forkhead protein family bind different DNA sites.

Unraveling determinants of transcription factor binding outside the core binding site.

Binding of transcription factors (TFs) to regulatory sequences is a pivotal step in the control of gene expression. Despite many advances in the characterization of sequence motifs recognized by TFs, our ability to quantitatively predict TF binding to different regulatory sequences is still limited. Here, we present a novel experimental assay termed BunDLE-seq that provides quantitative measurements of TF binding to thousands of fully designed sequences of 200 bp in length within a single experiment. Applying this binding assay to two yeast TFs, we demonstrate that sequences outside the core TF binding site profoundly affect TF binding. We show that TF-specific models based on the sequence or DNA shape of the regions flanking the core binding site are highly predictive of the measured differential TF binding. We further characterize the dependence of TF binding, accounting for measurements of single and co-occurring binding events, on the number and location of binding sites and on the TF concentration. Finally, by coupling our in vitro TF binding measurements, and another application of our method probing nucleosome formation, to in vivo expression measurements carried out with the same template sequences serving as promoters, we offer insights into mechanisms that may determine the different expression outcomes observed. Our assay thus paves the way to a more comprehensive understanding of TF binding to regulatory sequences and allows the characterization of TF binding determinants within and outside of core binding sites.

Evolving insights on how cytosine methylation affects protein-DNA binding.

Many anecdotal observations exist of a regulatory effect of DNA methylation on gene expression. However, in general, the underlying mechanisms of this effect are poorly understood. In this review, we summarize what is currently known about how this important, but mysterious, epigenetic mark impacts cellular functions. Cytosine methylation can abrogate or enhance interactions with DNA-binding proteins, or it may have no effect, depending on the context. Despite being only a small chemical change, the addition of a methyl group to cytosine can affect base readout via hydrophobic contacts in the major groove and shape readout via electrostatic contacts in the minor groove. We discuss the recent discovery that CpG methylation increases DNase I cleavage at adjacent positions by an order of magnitude through altering the local 3D DNA shape and the possible implications of this structural insight for understanding the methylation sensitivity of transcription factors (TFs). Additionally, 5-methylcytosines change the stability of nucleosomes and, thus, affect the local chromatin structure and access of TFs to genomic DNA. Given these complexities, it seems unlikely that the influence of DNA methylation on protein-DNA binding can be captured in a small set of general rules. Hence, data-driven approaches may be essential to gain a better understanding of these mechanisms.

Absence of a simple code: how transcription factors read the genome.

Transcription factors (TFs) influence cell fate by interpreting the regulatory DNA within a genome. TFs recognize DNA in a specific manner; the mechanisms underlying this specificity have been identified for many TFs based on 3D structures of protein-DNA complexes. More recently, structural views have been complemented with data from high-throughput in vitro and in vivo explorations of the DNA-binding preferences of many TFs. Together, these approaches have greatly expanded our understanding of TF-DNA interactions. However, the mechanisms by which TFs select in vivo binding sites and alter gene expression remain unclear. Recent work has highlighted the many variables that influence TF-DNA binding, while demonstrating that a biophysical understanding of these many factors will be central to understanding TF function.

Conformations of p53 response elements in solution deduced using site-directed spin labeling and Monte Carlo sampling.

The tumor suppressor protein p53 regulates numerous signaling pathways by specifically recognizing diverse p53 response elements (REs). Understanding the mechanisms of p53-DNA interaction requires structural information on p53 REs. However, such information is limited as a 3D structure of any RE in the unbound form is not available yet. Here, site-directed spin labeling was used to probe the solution structures of REs involved in p53 regulation of the p21 and Bax genes. Multiple nanometer distances in the p21-RE and BAX-RE, measured using a nucleotide-independent nitroxide probe and double-electron-electron-resonance spectroscopy, were used to derive molecular models of unbound REs from pools of all-atom structures generated by Monte-Carlo simulations, thus enabling analyses to reveal sequence-dependent DNA shape features of unbound REs in solution. The data revealed distinct RE conformational changes on binding to the p53 core domain, and support the hypothesis that sequence-dependent properties encoded in REs are exploited by p53 to achieve the energetically most favorable mode of deformation, consequently enhancing binding specificity. This work reveals mechanisms of p53-DNA recognition, and establishes a new experimental/computational approach for studying DNA shape in solution that has far-reaching implications for studying protein-DNA interactions.

Structure of p53 binding to the BAX response element reveals DNA unwinding and compression to accommodate base-pair insertion.

The p53 core domain binds to response elements (REs) that contain two continuous half-sites as a cooperative tetramer, but how p53 recognizes discontinuous REs is not well understood. Here we describe the crystal structure of the p53 core domain bound to a naturally occurring RE located at the promoter of the Bcl-2-associated X protein (BAX) gene, which contains a one base-pair insertion between the two half-sites. Surprisingly, p53 forms a tetramer on the BAX-RE that is nearly identical to what has been reported on other REs with a 0-bp spacer. Each p53 dimer of the tetramer binds in register to a half-site and maintains the same protein-DNA interactions as previously observed, and the two dimers retain all the protein-protein contacts without undergoing rotation or translation. To accommodate the additional base pair, the DNA is deformed and partially disordered around the spacer region, resulting in an apparent unwinding and compression, such that the interactions between the dimers are maintained. Furthermore, DNA deformation within the p53-bound BAX-RE is confirmed in solution by site-directed spin labeling measurements. Our results provide a structural insight into the mechanism by which p53 binds to discontinuous sites with one base-pair spacer.

DNAshape: a method for the high-throughput prediction of DNA structural features on a genomic scale.

We present a method and web server for predicting DNA structural features in a high-throughput (HT) manner for massive sequence data. This approach provides the framework for the integration of DNA sequence and shape analyses in genome-wide studies. The HT methodology uses a sliding-window approach to mine DNA structural information obtained from Monte Carlo simulations. It requires only nucleotide sequence as input and instantly predicts multiple structural features of DNA (minor groove width, roll, propeller twist and helix twist). The results of rigorous validations of the HT predictions based on DNA structures solved by X-ray crystallography and NMR spectroscopy, hydroxyl radical cleavage data, statistical analysis and cross-validation, and molecular dynamics simulations provide strong confidence in this approach. The DNAshape web server is freely available at http://rohslab.cmb.usc.edu/DNAshape/.

Probing DNA shape and methylation state on a genomic scale with DNase I.

DNA binding proteins find their cognate sequences within genomic DNA through recognition of specific chemical and structural features. Here we demonstrate that high-resolution DNase I cleavage profiles can provide detailed information about the shape and chemical modification status of genomic DNA. Analyzing millions of DNA backbone hydrolysis events on naked genomic DNA, we show that the intrinsic rate of cleavage by DNase I closely tracks the width of the minor groove. Integration of these DNase I cleavage data with bisulfite sequencing data for the same cell type's genome reveals that cleavage directly adjacent to cytosine-phosphate-guanine (CpG) dinucleotides is enhanced at least eightfold by cytosine methylation. This phenomenon we show to be attributable to methylation-induced narrowing of the minor groove. Furthermore, we demonstrate that it enables simultaneous mapping of DNase I hypersensitivity and regional DNA methylation levels using dense in vivo cleavage data. Taken together, our results suggest a general mechanism by which CpG methylation can modulate protein-DNA interaction strength via the remodeling of DNA shape.