Evaluating Hydrogen Bonds and Base Stacking of Single, Tandem and Terminal GU Mismatches in RNA with a Mesoscopic Model.

Guanine-Uracil (GU) mismatches are crucial to the stability of the RNA double helix and need to be considered in RNA folding algorithms for numerous biotechnological applications. Yet despite its importance, many aspects of GU base pairs are still poorly understood. There is also a lack of parametrization which prevents it to be considered in mesoscopic models. Here, we adapted the mesoscopic Peyrard-Bishop model to deal with context-dependent hydrogen bonds of GU mismatches and calculated the model parameters related to hydrogen bonding and base stacking from available experimental melting temperatures. The context-dependence causes a proliferation of parameters which made the problem computationally very demanding. We were able to overcome this problem by systematically regrouping the parameters during the minimization procedure. Our results not only provide the much needed parametrization but also answer several questions about the general properties of GU base pairs, as they can be associated straightforwardly to hydrogen bonding and base stacking. In particular, we found a very small Morse potential for tandem 5'-GU-3', which confirms a single hydrogen bond for this configuration, answering a long-standing question over conflicting experimental findings. Terminal GU base pairs are known to increase the duplex stability, but it is not clear why. Our results suggest that the increased terminal stability is mostly due to stronger hydrogen bonding.

DNA terminal base pairs have weaker hydrogen bonds especially for AT under low salt concentration.

DNA base pairs are known to open more easily at the helix terminal, a process usually called end fraying, the details of which are still poorly understood. Here, we present a mesoscopic model calculation based on available experimental data where we consider separately the terminal base pairs of a DNA duplex. Our results show an important reduction of hydrogen bond strength for terminal cytosine-guanine (CG) base pairs which is uniform over the whole range of salt concentrations, while for AT base pairs, we obtain a nearly 1/3 reduction but only at low salt concentrations. At higher salt concentrations, terminal adenine-thymine (AT) pair has almost the same hydrogen bond strength than interior bases. The calculated terminal stacking interaction parameters display some peculiarly contrasting behavior. While there is mostly no perceptible difference to internal stacking, for some cases, we observe an unusually strong dependence with salt concentration which does not appear follow any pattern or trend.

Effects of airway surface liquid height on the kinetics of extracellular nucleotides in airway epithelia.

Experimental techniques aimed at measuring the concentration of signaling molecules in the airway surface liquid (ASL) often require an unrealistically large ASL volume to facilitate sampling. This experimental limitation, prompted by the difficulty of pipetting liquid from a very shallow layer (~15 µm), leads to dilution and the under-prediction of physiologic concentrations of signaling molecules that are vital to the regulation of mucociliary clearance. Here, we use a computational model to describe the effect of liquid height on the kinetics of extracellular nucleotides in the airway surface liquid coating respiratory epithelia. The model consists of a reaction-diffusion equation with boundary conditions that represent the enzymatic reactions occurring on the epithelial surface. The simulations reproduce successfully the kinetics of extracellular ATP following hypotonic challenge for ASL volumes ranging from 25 µl to 500 µl in a 12-mm diameter cell culture. The model reveals that [ATP] and [ADO] reach 1200 nM and 2200 nM at the epithelial surface, respectively, while their volumetric averages remain less than 200 nM at all times in experiments with a large ASL volume (500 µl). These findings imply that activation of P2Y2 and A2B receptors is robust after hypotonic challenge, in contrast to what could be concluded based on experimental measurements of volumetric concentrations in large ASL volumes. Finally, given the central role that ATP and ADO play in regulating mucociliary clearance, we investigated which enzymes, when inhibited, provide the greatest increase in ATP and ADO concentrations. Our findings suggest that inhibition of NTPDase1/highTNAP would cause the greatest increase in [ATP] after hypotonic challenge, while inhibition of the transporter CNT3 would provide the greatest increase in [ADO].

Salt dependent mesoscopic model for RNA at multiple strand concentrations.

Mesoscopic models can be used for the description of the thermodynamic properties of RNA duplexes. With the use of experimental melting temperatures, its parametrization can provide important insights into its hydrogen bonds and stacking interactions as has been done for high sodium concentrations. However, the RNA parametrization for lower salt concentrations is still missing due to the limited amount of published melting temperature data. While the Peyrard-Bishop (PB) parametrization was found to be largely independent of strand concentrations, it requires that all temperatures are provided at the same strand concentrations. Here we adapted the PB model to handle multiple strand concentrations and in this way we were able to make use of an experimental set of temperatures to model the hydrogen bond and stacking interactions at low and intermediate sodium concentrations. For the parametrizations we make a distinction between terminal and internal base pairs, and the resulting potentials were qualitatively similar as we obtained previously for DNA. The main difference from DNA parameters, was the Morse potentials at low sodium concentrations for terminal r(AU) which is stronger than d(AT), suggesting higher hydrogen bond strength.

Spectrum of mutational signatures in T-cell lymphoma reveals a key role for UV radiation in cutaneous T-cell lymphoma.

T-cell non-Hodgkin's lymphomas develop following transformation of tissue resident T-cells. We performed a meta-analysis of whole exome sequencing data from 403 patients with eight subtypes of T-cell non-Hodgkin's lymphoma to identify mutational signatures and associated recurrent gene mutations. Signature 1, indicative of age-related deamination, was prevalent across all T-cell lymphomas, reflecting the derivation of these malignancies from memory T-cells. Adult T-cell leukemia-lymphoma was specifically associated with signature 17, which was found to correlate with the IRF4 K59R mutation that is exclusive to Adult T-cell leukemia-lymphoma. Signature 7, implicating UV exposure was uniquely identified in cutaneous T-cell lymphoma (CTCL), contributing 52% of the mutational burden in mycosis fungoides and 23% in Sezary syndrome. Importantly this UV signature was observed in CD4 + T-cells isolated from the blood of Sezary syndrome patients suggesting extensive re-circulation of these T-cells through skin and blood. Analysis of non-Hodgkin's T-cell lymphoma cases submitted to the national 100,000 WGS project confirmed that signature 7 was only identified in CTCL strongly implicating UV radiation in the pathogenesis of cutaneous T-cell lymphoma.