Invariant TAD Boundaries Constrain Cell-Type-Specific Looping Interactions between Promoters and Distal Elements around the CFTR Locus

Three-dimensional genome structure plays an important role in gene regulation. Globally, chromosomes are organized into active and inactive compartments while, at the gene level, looping interactions connect promoters to regulatory elements. Topologically associating domains (TADs), typically several hundred kilobases in size, form an intermediate level of organization. Major questions include how TADs are formed and how they are related to looping interactions between genes and regulatory elements. Here we performed a focused 5C analysis of a 2.8 Mb chromosome 7 region surrounding CFTR in a panel of cell types. We find that the same TAD boundaries are present in all cell types, indicating that TADs represent a universal chromosome architecture. Furthermore, we find that these TAD boundaries are present irrespective of the expression and looping of genes located between them. In contrast, looping interactions between promoters and regulatory elements are cell-type specific and occur mostly within TADs. This is exemplified by the CFTR promoter that in different cell types interacts with distinct sets of distal cell-type-specific regulatory elements that are all located within the same TAD. Finally, we find that long-range associations between loci located in different TADs are also detected, but these display much lower interaction frequencies than looping interactions within TADs. Interestingly, interactions between TADs are also highly cell-type-specific and often involve loci clustered around TAD boundaries. These data point to key roles of invariant TAD boundaries in constraining as well as mediating cell-type-specific long-range interactions and gene regulation.     

Guanidine-HCl Dependent Structural Unfolding of M-Crystallin: Fluctuating Native State Like Topologies and Intermolecular Association

Numerous experimental techniques and computational studies, proposed in recent times, have revolutionized the understanding of protein-folding paradigm. The complete understanding of protein folding and intermediates are of medical relevance, as the aggregation of misfolding proteins underlies various diseases, including some neurodegenerative disorders. Here, we describe the unfolding of M-crystallin, a βγ-crystallin homologue protein from archaea, from its native state to its denatured state using multidimensional NMR and other biophysical techniques. The protein, which was earlier characterized to be a predominantly β-sheet protein in its native state, shows different structural propensities (α and β), under different denaturing conditions. In 2 M GdmCl, the protein starts showing two distinct sets of peaks, with one arising from a partially unfolded state and the other from a completely folded state. The native secondary structural elements start disappearing as the denaturant concentration approaches 4 M. Subsequently, the protein is completely unfolded when the denaturant concentration is 6 M. The ¹⁵N relaxation data (T₁/T₂), heteronuclear ¹H-¹⁵N Overhauser effects (nOes), NOESY data, and other biophysical data taken together indicate that the protein shows a consistent, gradual change in its structural and motional preferences with increasing GdmCl concentration.     

Quantifying Missing Heritability at Known GWAS Loci

Recent work has shown that much of the missing heritability of complex traits can be resolved by estimates of heritability explained by all genotyped SNPs. However, it is currently unknown how much heritability is missing due to poor tagging or additional causal variants at known GWAS loci. Here, we use variance components to quantify the heritability explained by all SNPs at known GWAS loci in nine diseases from WTCCC1 and WTCCC2. After accounting for expectation, we observed all SNPs at known GWAS loci to explain more heritability than GWAS-associated SNPs on average (). For some diseases, this increase was individually significant: for Multiple Sclerosis (MS) () and for Crohn''s Disease (CD) (); all analyses of autoimmune diseases excluded the well-studied MHC region. Additionally, we found that GWAS loci from other related traits also explained significant heritability. The union of all autoimmune disease loci explained more MS heritability than known MS SNPs () and more CD heritability than known CD SNPs (), with an analogous increase for all autoimmune diseases analyzed. We also observed significant increases in an analysis of Rheumatoid Arthritis (RA) samples typed on ImmunoChip, with more heritability from all SNPs at GWAS loci () and more heritability from all autoimmune disease loci () compared to known RA SNPs (including those identified in this cohort). Our methods adjust for LD between SNPs, which can bias standard estimates of heritability from SNPs even if all causal variants are typed. By comparing adjusted estimates, we hypothesize that the genome-wide distribution of causal variants is enriched for low-frequency alleles, but that causal variants at known GWAS loci are skewed towards common alleles. These findings have important ramifications for fine-mapping study design and our understanding of complex disease architecture.     

Mapping DNA damage‐dependent genetic interactions in yeast via party mating and barcode fusion genetics

Condition‐dependent genetic interactions can reveal functional relationships between genes that are not evident under standard culture conditions. State‐of‐the‐art yeast genetic interaction mapping, which relies on robotic manipulation of arrays of double‐mutant strains, does not scale readily to multi‐condition studies. Here, we describe barcode fusion genetics to map genetic interactions (BFG‐GI), by which double‐mutant strains generated via en masse “party” mating can also be monitored en masse for growth to detect genetic interactions. By using site‐specific recombination to fuse two DNA barcodes, each representing a specific gene deletion, BFG‐GI enables multiplexed quantitative tracking of double mutants via next‐generation sequencing. We applied BFG‐GI to a matrix of DNA repair genes under nine different conditions, including methyl methanesulfonate (MMS), 4‐nitroquinoline 1‐oxide (4NQO), bleomycin, zeocin, and three other DNA‐damaging environments. BFG‐GI recapitulated known genetic interactions and yielded new condition‐dependent genetic interactions. We validated and further explored a subnetwork of condition‐dependent genetic interactions involving MAG1, SLX4, and genes encoding the Shu complex, and inferred that loss of the Shu complex leads to an increase in the activation of the checkpoint protein kinase Rad53.     

Energetic localization of saxitoxin in its channel binding site

Saxitoxin (STX) selectively blocks the voltage-gated sodium channel at the outer vestibule lined by P-loops of the four domains. Neosaxitoxin has an additional -OH group at the N1 position of the 1,2,3 guanidinium (N1-OH) that interacts with domains I and IV of the Na(+) channel. Determination of a second toxin interaction with the channel would fix the location of STX. Gonyautoxin 2,3 and Gonyautoxin 1,4 are C-11 sulfated derivatives of saxitoxin and neosaxitoxin, respectively. We used these variants to constrain the STX docking orientation by energetically localizing the C-11 sulfate in the outer vestibule. Interactions between the C-11 sulfate and each of the four domains of the channel were determined by a systematic approach to mutant cycle analysis in which all known carboxyl groups important for site 1 toxin binding were neutralized, allowing energetic triangulation of the toxin sulfate and overcoming some limitations of mutant cycles. Toxin IC(50)s were measured by two-electrode voltage clamp from Xenopus oocytes injected with the channel mRNA. Three unique types of analysis based on the coupling results localized the C-11 sulfate between domains III and IV. Combined with our previous report, the data establish the orientation of STX in the outer vestibule and confirm the clockwise arrangement of the channel domains.     

Molecular Dynamic Simulations Reveal the Structural Determinants of Fatty Acid Binding to Oxy-Myoglobin

The mechanism(s) by which fatty acids are sequestered and transported in muscle have not been fully elucidated. A potential key player in this process is the protein myoglobin (Mb). Indeed, there is a catalogue of empirical evidence supporting direct interaction of globins with fatty acid metabolites; however, the binding pocket and regulation of the interaction remains to be established. In this study, we employed a computational strategy to elucidate the structural determinants of fatty acids (palmitic & oleic acid) binding to Mb. Sequence analysis and docking simulations with a horse (Equus caballus) structural Mb reference reveals a fatty acid-binding site in the hydrophobic cleft near the heme region in Mb. Both palmitic acid and oleic acid attain a “U” shaped structure similar to their conformation in pockets of other fatty acid-binding proteins. Specifically, we found that the carboxyl head group of palmitic acid coordinates with the amino group of Lys45, whereas the carboxyl group of oleic acid coordinates with both the amino groups of Lys45 and Lys63. The alkyl tails of both fatty acids are supported by surrounding hydrophobic residues Leu29, Leu32, Phe33, Phe43, Phe46, Val67, Val68 and Ile107. In the saturated palmitic acid, the hydrophobic tail moves freely and occasionally penetrates deeper inside the hydrophobic cleft, making additional contacts with Val28, Leu69, Leu72 and Ile111. Our simulations reveal a dynamic and stable binding pocket in which the oxygen molecule and heme group in Mb are required for additional hydrophobic interactions. Taken together, these findings support a mechanism in which Mb acts as a muscle transporter for fatty acid when it is in the oxygenated state and releases fatty acid when Mb converts to deoxygenated state.