Examining DNA breathing with pyDNA-EPBD.

MOTIVATION: The two strands of the DNA double helix locally and spontaneously separate and recombine in living cells due to the inherent thermal DNA motion. This dynamics results in transient openings in the double helix and is referred to as "DNA breathing" or "DNA bubbles." The propensity to form local transient openings is important in a wide range of biological processes, such as transcription, replication, and transcription factors binding. However, the modeling and computer simulation of these phenomena, have remained a challenge due to the complex interplay of numerous factors, such as, temperature, salt content, DNA sequence, hydrogen bonding, base stacking, and others. RESULTS: We present pyDNA-EPBD, a parallel software implementation of the Extended Peyrard-Bishop-Dauxois (EPBD) nonlinear DNA model that allows us to describe some features of DNA dynamics in detail. The pyDNA-EPBD generates genomic scale profiles of average base-pair openings, base flipping probability, DNA bubble probability, and calculations of the characteristically dynamic length indicating the number of base pairs statistically significantly affected by a single point mutation using the Markov Chain Monte Carlo algorithm. AVAILABILITY AND IMPLEMENTATION: pyDNA-EPBD is supported across most operating systems and is freely available at https://github.com/lanl/pyDNA_EPBD. Extensive documentation can be found at https://lanl.github.io/pyDNA_EPBD/.

A rare functional noncoding variant at the GWAS-implicated MIR137/MIR2682 locus might confer risk to schizophrenia and bipolar disorder.

Schizophrenia (SZ) genome-wide association studies (GWASs) have identified common risk variants in >100 susceptibility loci; however, the contribution of rare variants at these loci remains largely unexplored. One of the strongly associated loci spans MIR137 (miR137) and MIR2682 (miR2682), two microRNA genes important for neuronal function. We sequenced ∼6.9 kb MIR137/MIR2682 and upstream regulatory sequences in 2,610 SZ cases and 2,611 controls of European ancestry. We identified 133 rare variants with minor allele frequency (MAF) <0.5%. The rare variant burden in promoters and enhancers, but not insulators, was associated with SZ (p = 0.021 for MAF < 0.5%, p = 0.003 for MAF < 0.1%). A rare enhancer SNP, 1:g.98515539A>T, presented exclusively in 11 SZ cases (nominal p = 4.8 × 10(-4)). We further identified its risk allele T in 2 of 2,434 additional SZ cases, 11 of 4,339 bipolar (BP) cases, and 3 of 3,572 SZ/BP study controls and 1,688 population controls; yielding combined p values of 0.0007, 0.0013, and 0.0001 for SZ, BP, and SZ/BP, respectively. The risk allele T of 1:g.98515539A>T reduced enhancer activity of its flanking sequence by >50% in human neuroblastoma cells, predicting lower expression of MIR137/MIR2682. Both empirical and computational analyses showed weaker transcription factor (YY1) binding by the risk allele. Chromatin conformation capture (3C) assay further indicated that 1:g.98515539A>T influenced MIR137/MIR2682, but not the nearby DPYD or LOC729987. Our results suggest that rare noncoding risk variants are associated with SZ and BP at MIR137/MIR2682 locus, with risk alleles decreasing MIR137/MIR2682 expression.

CD39 and CD161 modulate Th17 responses in Crohn's disease.

CD39 (ENTPD1) is expressed by subsets of pathogenic human CD4(+) T cells, such as Th17 cells. These Th17 cells are considered important in intestinal inflammation, such as seen in Crohn's disease (CD). Recently, CD161 (NKR-P1A) was shown to be a phenotypic marker of human Th17 cells. In this study, we report that coexpression of CD161 and CD39 not only identifies these cells but also promotes Th17 generation. We note that human CD4(+)CD39(+)CD161(+) T cells can be induced under stimulatory conditions that promote Th17 in vitro. Furthermore, CD4(+)CD39(+)CD161(+) cells purified from blood and intestinal tissues, from both healthy controls and patients with CD, are of the Th17 phenotype and exhibit proinflammatory functions. CD39 is coexpressed with CD161, and this association augments acid sphingomyelinase (ASM) activity upon stimulation of CD4(+) T cells. These pathways regulate mammalian target of rapamycin and STAT3 signaling to drive the Th17 phenotype. Inhibition of ASM activity by pharmacological blockers or knockdown of ASM abrogates STAT3 signaling, thereby limiting IL-17 production in CD4(+) T cells obtained from both controls and patients with active CD. Increased levels of CD39(+)CD161(+) CD4(+) T cells in blood or lamina propria are noted in patients with CD, and levels directly correlate with clinical disease activity. Hence, coexpression of CD39 and CD161 by CD4(+) T cells might serve as a biomarker to monitor Th17 responsiveness. Collectively, CD39 and CD161 modulate human Th17 responses in CD through alterations in purinergic nucleotide-mediated responses and ASM catalytic bioactivity, respectively.

Lysophosphatidic acid generation by pulmonary NKT cell ENPP-2/autotaxin exacerbates hyperoxic lung injury.

Hyperoxia is still broadly used in clinical practice in order to assure organ oxygenation in critically ill patients, albeit known toxic effects. In this present study, we hypothesize that lysophosphatidic acid (LPA) mediates NKT cell activation in a mouse model of hyperoxic lung injury. In vitro, pulmonary NKT cells were exposed to hyperoxia for 72 h, and the induction of the ectonucleotide pyrophosphatase/phosphodiesterase 2 (ENPP-2) was examined and production of lysophosphatidic acid (LPA) was measured. In vivo, animals were exposed to 100 % oxygen for 72 h and lungs and serum were harvested. Pulmonary NKT cells were then incubated with the LPA antagonist Brp-LPA. Animals received BrP-LPA prior to oxygen exposure. Autotaxin (ATX, ENPP-2) was significantly up-regulated on pulmonary NKT cells after hyperoxia (p < 0.01) in vitro. LPA levels were increased in supernatants of hyperoxia-exposed pulmonary NKT cells. LPA levels were significantly reduced by incubating NKT cells with LPA-BrP during oxygen exposure (p < 0,05) in vitro. Hyperoxia-exposed animals showed significantly increased serum levels of LPA (p ≤ 0,05) as well as increased pulmonary NKT cell numbers in vivo. BrP-LPA injection significantly improved survival as well as significantly decreased lung injury and lowered pulmonary NKT cell numbers. We conclude that NKT cell-induced hyperoxic lung injury is mediated by pro-inflammatory LPA generation, at least in part, secondary to ENPP-2 up-regulation on pulmonary NKT cells. Being a potent LPA antagonist, BrP-LPA prevents hyperoxia-induced lung injury in vitro and in vivo.

Soluble ecto-5'-nucleotidase (5'-NT), alkaline phosphatase, and adenosine deaminase (ADA1) activities in neonatal blood favor elevated extracellular adenosine.

Extracellular adenosine, a key regulator of physiology and immune cell function that is found at elevated levels in neonatal blood, is generated by phosphohydrolysis of adenine nucleotides released from cells and catabolized by deamination to inosine. Generation of adenosine monophosphate (AMP) in blood is driven by cell-associated enzymes, whereas conversion of AMP to adenosine is largely mediated by soluble enzymes. The identities of the enzymes responsible for these activities in whole blood of neonates have been defined in this study and contrasted to adult blood. We demonstrate that soluble 5'-nucleotidase (5'-NT) and alkaline phosphatase (AP) mediate conversion of AMP to adenosine, whereas soluble adenosine deaminase (ADA) catabolizes adenosine to inosine. Newborn blood plasma demonstrates substantially higher adenosine-generating 5'-NT and AP activity and lower adenosine-metabolizing ADA activity than adult plasma. In addition to a role in soluble purine metabolism, abundant AP expressed on the surface of circulating neonatal neutrophils is the dominant AMPase on these cells. Plasma samples from infant observational cohorts reveal a relative plasma ADA deficiency at birth, followed by a gradual maturation of plasma ADA through infancy. The robust adenosine-generating capacity of neonates appears functionally relevant because supplementation with AMP inhibited whereas selective pharmacologic inhibition of 5'-NT enhanced Toll-like receptor-mediated TNF-α production in neonatal whole blood. Overall, we have characterized previously unrecognized age-dependent expression patterns of plasma purine-metabolizing enzymes that result in elevated plasma concentrations of anti-inflammatory adenosine in newborns. Targeted manipulation of purine-metabolizing enzymes may benefit this vulnerable population.

Disordered purinergic signaling and abnormal cellular metabolism are associated with development of liver cancer in Cd39/ENTPD1 null mice.

UNLABELLED: Liver cancer is associated with chronic inflammation, which is linked to immune dysregulation, disordered metabolism, and aberrant cell proliferation. Nucleoside triphosphate diphosphohydrolase-1; (CD39/ENTPD1) is an ectonucleotidase that regulates extracellular nucleotide/nucleoside concentrations by scavenging nucleotides to ultimately generate adenosine. These properties inhibit antitumor immune responses and promote angiogenesis, being permissive for the growth of transplanted tumors. Here we show that Cd39 deletion promotes development of both induced and spontaneous autochthonous liver cancer in mice. Loss of Cd39 results in higher concentrations of extracellular nucleotides, which stimulate proliferation of hepatocytes, abrogate autophagy, and disrupt glycolytic metabolism. Constitutive activation of Ras-mitogen-activated protein kinase (MAPK) and mammalian target of rapamycin (mTOR)-S6K1 pathways occurs in both quiescent Cd39 null hepatocytes in vitro and liver tissues in vivo. Exogenous adenosine 5'-triphosphate (ATP) boosts these signaling pathways, whereas rapamycin inhibits such aberrant responses in hepatocytes. CONCLUSION: Deletion of Cd39 and resulting changes in disordered purinergic signaling perturb hepatocellular metabolic/proliferative responses, paradoxically resulting in malignant transformation. These findings might impact adjunctive therapies for cancer. Our studies indicate that the biology of autochthonous and transplanted tumors is quite distinct.

Binding of nucleoid-associated protein fis to DNA is regulated by DNA breathing dynamics.

Physicochemical properties of DNA, such as shape, affect protein-DNA recognition. However, the properties of DNA that are most relevant for predicting the binding sites of particular transcription factors (TFs) or classes of TFs have yet to be fully understood. Here, using a model that accurately captures the melting behavior and breathing dynamics (spontaneous local openings of the double helix) of double-stranded DNA, we simulated the dynamics of known binding sites of the TF and nucleoid-associated protein Fis in Escherichia coli. Our study involves simulations of breathing dynamics, analysis of large published in vitro and genomic datasets, and targeted experimental tests of our predictions. Our simulation results and available in vitro binding data indicate a strong correlation between DNA breathing dynamics and Fis binding. Indeed, we can define an average DNA breathing profile that is characteristic of Fis binding sites. This profile is significantly enriched among the identified in vivo E. coli Fis binding sites. To test our understanding of how Fis binding is influenced by DNA breathing dynamics, we designed base-pair substitutions, mismatch, and methylation modifications of DNA regions that are known to interact (or not interact) with Fis. The goal in each case was to make the local DNA breathing dynamics either closer to or farther from the breathing profile characteristic of a strong Fis binding site. For the modified DNA segments, we found that Fis-DNA binding, as assessed by gel-shift assay, changed in accordance with our expectations. We conclude that Fis binding is associated with DNA breathing dynamics, which in turn may be regulated by various nucleotide modifications.

DNA breathing dynamics distinguish binding from nonbinding consensus sites for transcription factor YY1 in cells.

The genome-wide mapping of the major gene expression regulators, the transcription factors (TFs) and their DNA binding sites, is of great importance for describing cellular behavior and phenotypic diversity. Presently, the methods for prediction of genomic TF binding produce a large number of false positives, most likely due to insufficient description of the physiochemical mechanisms of protein-DNA binding. Growing evidence suggests that, in the cell, the double-stranded DNA (dsDNA) is subject to local transient strands separations (breathing) that contribute to genomic functions. By using site-specific chromatin immunopecipitations, gel shifts, BIOBASE data, and our model that accurately describes the melting behavior and breathing dynamics of dsDNA we report a specific DNA breathing profile found at YY1 binding sites in cells. We find that the genomic flanking sequence variations and SNPs, may exert long-range effects on DNA dynamics and predetermine YY1 binding. The ubiquitous TF YY1 has a fundamental role in essential biological processes by activating, initiating or repressing transcription depending upon the sequence context it binds. We anticipate that consensus binding sequences together with the related DNA dynamics profile may significantly improve the accuracy of genomic TF binding sites and TF binding-related functional SNPs.

Advancing Transcription Factor Binding Site Prediction Using DNA Breathing Dynamics and Sequence Transformers via Cross Attention.

Understanding the impact of genomic variants on transcription factor binding and gene regulation remains a key area of research, with implications for unraveling the complex mechanisms underlying various functional effects. Our study delves into the role of DNA's biophysical properties, including thermodynamic stability, shape, and flexibility in transcription factor (TF) binding. We developed a multi-modal deep learning model integrating these properties with DNA sequence data. Trained on ChIP-Seq (chromatin immunoprecipitation sequencing) data in vivo involving 690 TF-DNA binding events in human genome, our model significantly improves prediction performance in over 660 binding events, with up to 9.6% increase in AUROC metric compared to the baseline model when using no DNA biophysical properties explicitly. Further, we expanded our analysis to in vitro high-throughput Systematic Evolution of Ligands by Exponential enrichment (SELEX) and Protein Binding Microarray (PBM) datasets, comparing our model with established frameworks. The inclusion of DNA breathing features consistently improved TF binding predictions across different cell lines in these datasets. Notably, for complex ChIP-Seq datasets, integrating DNABERT2 with a cross-attention mechanism provided greater predictive capabilities and insights into the mechanisms of disease-related non-coding variants found in genome-wide association studies. This work highlights the importance of DNA biophysical characteristics in TF binding and the effectiveness of multi-modal deep learning models in gene regulation studies.

Essential roles of Raf/extracellular signal-regulated kinase/mitogen-activated protein kinase pathway, YY1, and Ca2+ influx in growth arrest of human vascular smooth muscle cells by bilirubin.

The biological effects of bilirubin, still poorly understood, are concentration-dependent ranging from cell protection to toxicity. Here we present data that at high nontoxic physiological concentrations, bilirubin inhibits growth of proliferating human coronary artery smooth muscle cells by three events. It impairs the activation of Raf/ERK/MAPK pathway and the cellular Raf and cyclin D1 content that results in retinoblastoma protein hypophosphorylation on amino acids S608 and S780. These events impede the release of YY1 to the nuclei and its availability to regulate the expression of genes and to support cellular proliferation. Moreover, altered calcium influx and calpain II protease activation leads to proteolytical degradation of transcription factor YY1. We conclude that in the serum-stimulated human vascular smooth muscle primary cell cultures, bilirubin favors growth arrest, and we propose that this activity is regulated by its interaction with the Raf/ERK/MAPK pathway, effect on cyclin D1 and Raf content, altered retinoblastoma protein profile of hypophosphorylation, calcium influx, and YY1 proteolysis. We propose that these activities together culminate in diminished 5 S and 45 S ribosomal RNA synthesis and cell growth arrest. The observations provide important mechanistic insight into the molecular mechanisms underlying the transition of human vascular smooth muscle cells from proliferative to contractile phenotype and the role of bilirubin in this transition.

Adenosine generation catalyzed by CD39 and CD73 expressed on regulatory T cells mediates immune suppression.

The study of T regulatory cells (T reg cells) has been limited by the lack of specific surface markers and an inability to define mechanisms of suppression. We show that the expression of CD39/ENTPD1 in concert with CD73/ecto-5'-nucleotidase distinguishes CD4(+)/CD25(+)/Foxp3(+) T reg cells from other T cells. These ectoenzymes generate pericellular adenosine from extracellular nucleotides. The coordinated expression of CD39/CD73 on T reg cells and the adenosine A2A receptor on activated T effector cells generates immunosuppressive loops, indicating roles in the inhibitory function of T reg cells. Consequently, T reg cells from Cd39-null mice show impaired suppressive properties in vitro and fail to block allograft rejection in vivo. We conclude that CD39 and CD73 are surface markers of T reg cells that impart a specific biochemical signature characterized by adenosine generation that has functional relevance for cellular immunoregulation.

Complement factors are secreted in human follicular fluid by granulosa cells and are possible oocyte maturation factors.

AIMS: In this study, we identify components of the complement system present in human follicular fluid that affect oocyte development and maturation. MATERIAL AND METHODS: Using bottom-up liquid chromatography/mass spectrometry/mass spectrometry, we identified complement factors as consistently present in human follicular fluid from 15 different subjects. RESULTS: According to our gene-chip data, these complement factors are actively produced by granulosa cells. CONCLUSIONS: By applying the computational Ingenuity Pathway Analysis software and database we have identified complement pathways that play a role in oocyte maturation and follicular development.

Impaired cardiac glycolysis and glycogen depletion are linked to poor myocardial outcomes in juvenile male swine with metabolic syndrome and ischemia.

Obesity continues to rise in the juveniles and obese children are more likely to develop metabolic syndrome (MetS) and related cardiovascular disease. Unfortunately, effective prevention and long-term treatment options remain limited. We determined the juvenile cardiac response to MetS in a swine model. Juvenile male swine were fed either an obesogenic diet, to induce MetS, or a lean diet, as a control (LD). Myocardial ischemia was induced with surgically placed ameroid constrictor on the left circumflex artery. Physiological data were recorded and at 22 weeks of age the animals underwent a terminal harvest procedure and myocardial tissue was extracted for total metabolic and proteomic LC/MS-MS, RNA-seq analysis, and data underwent nonnegative matrix factorization for metabolic signatures. Significantly altered in MetS versus. LD were the glycolysis-related metabolites and enzymes. In MetS compared with LD Glycogen synthase 1 (GYS1)-glycogen phosphorylases (PYGM/PYGL) expression disbalance resulted in a loss of myocardial glycogen. Our findings are consistent with the concept that transcriptionally driven myocardial changes in glycogen and glucose metabolism-related enzymes lead to a deficiency of their metabolite products in MetS. This abnormal energy metabolism provides insight into the pathogenesis of the juvenile heart in MetS. This study reveals that MetS and ischemia diminishes ATP availability in the myocardium via altering the glucose-G6P-pyruvate axis at the level of metabolites and gene expression of related enzymes. The observed severe glycogen depletion in MetS coincides with disbalance in expression of GYS1 and both PYGM and PYGL. This altered energy substrate metabolism is a potential target of pharmacological agents for improving juvenile myocardial function in MetS and ischemia.

Examining DNA Breathing with pyDNA-EPBD.

Motivation: The two strands of the DNA double helix locally and spontaneously separate and recombine in living cells due to the inherent thermal DNA motion.This dynamics results in transient openings in the double helix and is referred to as "DNA breathing" or "DNA bubbles." The propensity to form local transient openings is important in a wide range of biological processes, such as transcription, replication, and transcription factors binding. However, the modeling and computer simulation of these phenomena, have remained a challenge due to the complex interplay of numerous factors, such as, temperature, salt content, DNA sequence, hydrogen bonding, base stacking, and others. Results: We present pyDNA-EPBD, a parallel software implementation of the Extended Peyrard-Bishop- Dauxois (EPBD) nonlinear DNA model that allows us to describe some features of DNA dynamics in detail. The pyDNA-EPBD generates genomic scale profiles of average base-pair openings, base flipping probability, DNA bubble probability, and calculations of the characteristically dynamic length indicating the number of base pairs statistically significantly affected by a single point mutation using the Markov Chain Monte Carlo (MCMC) algorithm.

Canagliflozin improves coronary microvascular vasodilation and increases absolute blood flow to the myocardium independent of angiogenesis.

OBJECTIVES: Sodium-glucose cotransporter-2 inhibitor, canagliflozin, improves myocardial perfusion to ischemic territory without accompanying changes in vascular density. We aimed to (1) characterize effects on angiogenic pathways, (2) use multiomics to identify gene expression and metabolite profiles relevant to regulation of myocardial blood flow, and (3) investigate drug effect on coronary microvascular reactivity. METHODS: A nondiabetic swine model of chronic myocardial ischemia and nondiabetic rat model were used to study functional and molecular effects of canagliflozin on myocardium and in vitro microvascular reactivity. RESULTS: Canagliflozin resulted in increased coronary microvascular vasodilation and decreased vasoconstriction (P < .05) without changes in microvascular density (P > .3). Expression of the angiogenic modulator, endostatin, increased (P = .008), along with its precursor, collagen 18 (P < .001), and factors that increase its production, including cathepsin L (P = .004). Endostatin and collagen 18 levels trended toward an inverse correlation with blood flow to ischemic territory at rest. Proangiogenic fibroblast growth factor receptor was increased (P = .03) and matrix metalloproteinase-9 was decreased (P < .001) with canagliflozin treatment. Proangiogenic vascular endothelial growth factor A (P = .13), Tie-2 (P = .10), and Ras (P = .18) were not significantly altered. Gene expression related to the cardiac renin-angiotensin system was significantly decreased. CONCLUSIONS: In chronic myocardial ischemia, canagliflozin increased absolute blood flow to the myocardium without robustly increasing vascular density or proangiogenic signaling. Canagliflozin resulted in altered coronary microvascular reactivity to favor vasodilation, likely through direct effect on vascular smooth muscle. Downregulation of cardiac renin-angiotensin system demonstrated local regulation of perfusion. VIDEO ABSTRACT.

Effects of canagliflozin on myocardial microvascular density, oxidative stress, and proteomic profile☆.

Introduction: Sodium-glucose cotransporter-2 inhibitors (SGLT2i) are cardioprotective, and canagliflozin (CANA), an SGLT2i, has been shown to improve perfusion, AMPK signaling, and oxidative stress in chronically ischemic myocardium. The aim of this study is to determine the effects of CANA in nonischemic myocardium on coronary collateralization, oxidative stress, and other molecular pathways determined by proteomic profiling. Methods: Yorkshire swine underwent placement of an ameroid constrictor to the left circumflex artery. Two weeks later, pigs received no drug (CON, n = 8) or 300 mg CANA daily (n = 8). Treatment continued for five weeks, followed by tissue harvest of nonischemic myocardium. Results: CANA was associated with decreased capillary density (p = 0.05) compared to CON, without changes in arteriolar density. Reduced capillary density did not correlate with reduced perfusion. Oxidative stress was reduced with CANA (22 % decrease). In the CANA group, there was a trend towards increased p-eNOS and eNOS, without a change in p-eNOS/eNOS ratio, p-Akt, Akt, and p-Akt/Akt ratio. There was no change in p-ERK1/2, but a decrease in total ERK1/2 and increase in p-ERK1/2/ERK1/2 ratio. There were no changes in expression of p-AMPK, AMPK, with a trend towards increased ratio of p-AMPK/AMPK. Proteomics analysis identified 2819 common proteins, of which 120 were upregulated and 425 were downregulated with CANA. Pathway analysis demonstrated wide regulation of metabolic proteins. Conclusions: The effects of CANA on myocardial perfusion and AMPK signaling in chronically ischemic myocardium are not found in nonischemic territory, despite attenuation of oxidative stress. Metabolic proteins are widely regulated in nonischemic myocardium with CANA.

DNA dynamics play a role as a basal transcription factor in the positioning and regulation of gene transcription initiation.

We assess the role of DNA breathing dynamics as a determinant of promoter strength and transcription start site (TSS) location. We compare DNA Langevin dynamic profiles of representative gene promoters, calculated with the extended non-linear PBD model of DNA with experimental data on transcription factor binding and transcriptional activity. Our results demonstrate that DNA dynamic activity at the TSS can be suppressed by mutations that do not affect basal transcription factor binding-DNA contacts. We use this effect to establish the separate contributions of transcription factor binding and DNA dynamics to transcriptional activity. Our results argue against a purely 'transcription factor-centric' view of transcription initiation, suggesting that both DNA dynamics and transcription factor binding are necessary conditions for transcription initiation.

Carbon monoxide suppresses arteriosclerotic lesions associated with chronic graft rejection and with balloon injury.

Carbon monoxide (CO), one of the products of heme oxygenase action on heme, prevents arteriosclerotic lesions that occur following aorta transplantation; pre-exposure to 250 parts per million of CO for 1 hour before injury suppresses stenosis after carotid balloon injury in rats as well as in mice. The protective effect of CO is associated with a profound inhibition of graft leukocyte infiltration/activation as well as with inhibition of smooth muscle cell proliferation. The anti-proliferative effect of CO in vitro requires the activation of guanylate cyclase, the generation of cGMP, the activation of p38 mitogen-activated protein kinases and the expression of the cell cycle inhibitor p21Cip1. These findings demonstrate a protective role for CO in vascular injury and support its use as a therapeutic agent.

A nonlinear dynamic model of DNA with a sequence-dependent stacking term.

No simple model exists that accurately describes the melting behavior and breathing dynamics of double-stranded DNA as a function of nucleotide sequence. This is especially true for homogenous and periodic DNA sequences, which exhibit large deviations in melting temperature from predictions made by additive thermodynamic contributions. Currently, no method exists for analysis of the DNA breathing dynamics of repeats and of highly G/C- or A/T-rich regions, even though such sequences are widespread in vertebrate genomes. Here, we extend the nonlinear Peyrard-Bishop-Dauxois (PBD) model of DNA to include a sequence-dependent stacking term, resulting in a model that can accurately describe the melting behavior of homogenous and periodic sequences. We collect melting data for several DNA oligos, and apply Monte Carlo simulations to establish force constants for the 10 dinucleotide steps (CG, CA, GC, AT, AG, AA, AC, TA, GG, TC). The experiments and numerical simulations confirm that the GG/CC dinucleotide stacking is remarkably unstable, compared with the stacking in GC/CG and CG/GC dinucleotide steps. The extended PBD model will facilitate thermodynamic and dynamic simulations of important genomic regions such as CpG islands and disease-related repeats.

The cardiac molecular setting of metabolic syndrome in pigs reveals disease susceptibility and suggests mechanisms that exacerbate COVID-19 outcomes in patients.

Although metabolic syndrome (MetS) is linked to an elevated risk of cardiovascular disease (CVD), the cardiac-specific risk mechanism is unknown. Obesity, hypertension, and diabetes (all MetS components) are the most common form of CVD and represent risk factors for worse COVID-19 outcomes compared to their non MetS peers. Here, we use obese Yorkshire pigs as a highly relevant animal model of human MetS, where pigs develop the hallmarks of human MetS and reproducibly mimics the myocardial pathophysiology in patients. Myocardium-specific mass spectroscopy-derived metabolomics, proteomics, and transcriptomics enabled the identity and quality of proteins and metabolites to be investigated in the myocardium to greater depth. Myocardium-specific deregulation of pro-inflammatory markers, propensity for arterial thrombosis, and platelet aggregation was revealed by computational analysis of differentially enriched pathways between MetS and control animals. While key components of the complement pathway and the immune response to viruses are under expressed, key N6-methyladenosin RNA methylation enzymes are largely overexpressed in MetS. Blood tests do not capture the entirety of metabolic changes that the myocardium undergoes, making this analysis of greater value than blood component analysis alone. Our findings create data associations to further characterize the MetS myocardium and disease vulnerability, emphasize the need for a multimodal therapeutic approach, and suggests a mechanism for observed worse outcomes in MetS patients with COVID-19 comorbidity.