Mesoscopic model for DNA G-quadruplex unfolding

[EN] Genomes contain rare guanine-rich sequences capable of assembling into four-stranded helical structures, termed G-quadruplexes, with potential roles in gene regulation and chromosome stability. Their mechanical unfolding has only been reported to date by all-atom simulations, which cannot dissect the major physical interactions responsible for their cohesion. Here, we propose a mesoscopic model to describe both the mechanical and thermal stability of DNA G-quadruplexes, where each nucleotide of the structure, as well as each central cation located at the inner channel, is mapped onto a single bead. In this framework we are able to simulate loading rates similar to the experimental ones, which are not reachable in simulations with atomistic resolution. In this regard, we present single-molecule force-induced unfolding experiments by a high-resolution optical tweezers on a DNA telomeric sequence capable of adopting a G-quadruplex conformation. Fitting the parameters of the model to the experiments we find a correct prediction of the rupture-force kinetics and a good agreement with previous near equilibrium measurements. Since G-quadruplex unfolding dynamics is halfway in complexity between secondary nucleic acids and tertiary protein structures, our model entails a nanoscale paradigm for non-equilibrium processes in the cell.Work supported by the Spanish Ministry of Economy and Competitiveness (MINECO), grant No. FIS2014-55867, co-financed by FEDER funds. We also thank the support of the Aragon Government and Fondo Social Europeo to FENOL group. Work in J.R.A.-G. laboratory was supported by a grant from MINECO, No. MAT2015-71806-R).Bergues-Pupo, A.; Gutiérrez, I.; Arias-Gonzalez, JR.; Falo, F.; Fiasconaro, A. (2017). Mesoscopic model for DNA G-quadruplex unfolding. Scientific Reports. 7:1-13. https://doi.org/10.1038/s41598-017-10849-2S1137Arias-Gonzalez, J. R. Single-molecule portrait of DNA and RNA double helices. Integr. Biol. 6, 904 (2014).Burge, S., Parkinson, G. N., Hazel, P., Todd, A. K. & Neidle, S. Quadruplex DNA: sequence, topology and structure. Nucleic Acids Res. 34, 5402 (2006).Lam, E. Y., Beraldi, D., Tannahill, D. & Balasubramanian, S. G-quadruplex structures are stable and detectable in human genomic DNA. Nat. Commun. 4, 1796 (2013).Siddiqui-Jain, A., Grand, C. L., Bearss, D. J. & Hurley, L. H. Direct evidence for a G-quadruplex in a promoter region and its targeting with a small molecule to repress c-MYC transcription. Proc. Natl. Acad. Sci. USA 99, 11593 (2002).Endoh, T. & Sugimoto, N. Mechanical insights into ribosomal progression overcoming RNA G-quadruplex from periodical translation suppression in cells. Sci. Rep. 6, 1 (2016).Hänsel-Hertsch, R., Di Antonio, M. & Balasubramanian, S. DNA G-quadruplexes in the human genome: detection, functions and therapeutic potential. Nat. Rev. Mol. Cell Biol. 18, 279 (2017).de Messieres, M., Chang, J. C., Brawn-Cinani, B. & La Porta, A. Single-molecule study of G-quadruplex disruption using dynamic force spectroscopy. Phys. Rev. Lett. 109, 058101 (2012).Koirala, D. et al. A single-molecule platform for investigation of interactions between G-quadruplexes and small-molecule ligands. Nat. Chem. 3, 782 (2011).Long, X. et al. Mechanical unfolding of human telomere G-quadruplex DNA probed by integrated fluorescence and magnetic tweezers spectroscopy. Nucleic Acids Res. 41, 2746 (2013).Ghimire, C. et al. Direct Quantification of Loop Interaction and pi-pi Stacking for G-Quadruplex Stability at the Submolecular Level. J. Am. Chem. Soc. 136, 15544 (2014).Garavís, M. et al. Mechanical Unfolding of Long Human Telomeric RNA (TERRA). Chem. Commun. 49, 6397 (2013).Fonseca Guerra, C., Zijlstra, H., Paragi, G. & Bickelhaupt, F. M. Telomere structure and stability: covalency in hydrogen bonds, not resonance assistance, causes cooperativity in guanine quartets. Chemistry-A European Journal 17, 12612 (2011).Yurenko, Y. P., Novotn, J., Sklen, V. & Marek, R. Exploring non-covalent interactions in guanine-and xanthine-based model DNA quadruplex structures: a comprehensive quantum chemical approach. Phys. Chem. Chem. Phys. 16, 2072 (2014).Poudel, L. et al. Implication of the solvent effect, metal ions and topology in the electronic structure and hydrogen bonding of human telomeric G-quadruplex DNA. Phys. Chem. Chem. Phys. 18, 21573 (2016).Li, M. H., Luo, Q., Xue, X. G. & Li, Z. S. Toward a full structural characterization of G-quadruplex DNA in aqueous solution: Molecular dynamics simulations of four G-quadruplex molecules. J. Mol. Struct-Theochem. 952, 96 (2010).Islam, B. et al. Conformational dynamics of the human propeller telomeric DNA quadruplex on a microsecond time scale. Nucleic Acids Res. 41, 2723 (2013).Stadlbauer, P., Krepl, M., Cheatham, T. E., Koca, J. & Sponer, J. Structural dynamics of possible late-stage intermediates in folding of quadruplex DNA studied by molecular simulations. Nucleic Acids Res. 41, 7128 (2013).Li, H., Cao, E. & Gisler, T. Force-induced unfolding of human telomeric G-quadruplex: a steered molecular dynamics simulation study. Biochem. Bioph. Res. Co. 379, 70 (2009).Yang, C., Jang, S. & Pak, Y. Multiple stepwise pattern for potential of mean force in unfolding the thrombin binding aptamer in complex with Sr2+. J. Chem. Phys. 135, 225104 (2011).Bergues-Pupo, A. E., Arias-Gonzalez, J. R., Morón, M. C., Fiasconaro, A. & Falo, F. Role of the central cations in the mechanical unfolding of DNA and RNA G-quadruplexes. Nucleic Acids Res. 43, 7638 (2015).Linak, M. C., Tourdot, R. & Dorfman, K. D. Moving beyond Watson-Crick models of coarse grained DNA dynamics. J. Chem Phys. 135, 205102 (2011).Rebi, M., Mocci, F., Laaksonen, A. & Ulin, J. Multiscale simulations of human telomeric G-quadruplex DNA. J. Phys. Chem. B 119, 105 (2014).Stadlbauer, P. et al. Coarse-Grained Simulations Complemented by Atomistic Molecular Dynamics Provide New Insights into Folding and Unfolding of Human Telomeric G-Quadruplexes. J. Chem. Theory Comput. 12, 6077 (2016).Parkinson, G. N., Lee, M. P. & Neidle, S. Crystal structure of parallel quadruplexes from human telomeric DNA. Nature 417, 876 (2002).Bhattacharya, D., Arachchilageand, G. M. & Basu, S. Metal Cations in G-Quadruplex Folding and Stability. Frontiers in Chemistry 4, 38 (2016).de Lorenzo, S., Ribezzi-Crivellari, M., Arias-Gonzalez, J. R., Smith, S. B. & Ritort, F. A Temperature-Jump Optical Trap for Single-Molecule Manipulation. Biophys. J. 108, 2854 (2015).Smith, S. B., Cui, Y. & Bustamante, C. Optical-trap force transducer that operates by direct measurement of light momentum. Methods Enzymol. 361, 134 (2003).Mergny, J. L., Phan, A. T. & Lacroix, L. Following G-quartet formation by UV-spectroscopy. FEBS letters 435, 74 (1998).Torrie, G. M. & Valleau, J. P. Nonphysical sampling distributions in Monte Carlo free-energy estimation: umbrella sampling. J. Comput. Phys. 23, 187 (1977).Kumar, S., Bouzida, D., Swendsen, R. H., Kollman, P. A. & Rosenberg, J. M. The weighted histogram analysis method for free-energy calculations on biomolecules I. The method. J. Comput. Chem. 13, 1011 (1992).Evans, E. & Ritchie, K. Dynamic strength of molecular adhesion bonds. Biophys. J. 72, 1541 (1997).Dudko, O. K., Hummer, G. & Szabo, A. Intrinsic rates and activation free energies from single-molecule pulling experiments. Phys. Rev. Lett. 96, 108101 (2006).Friddle, R. W., Noy, A. & De Yoreo, J. J. Interpreting the widespread nonlinear force spectra of intermolecular bonds. Proc. Natl. Acad. Sci. 109, 13573 (2012)

The katG mRNA of Mycobacterium tuberculosis and Mycobacterium smegmatis is processed at its 5' end and is stabilized by both a polypurine sequence and translation initiation

<p>Abstract</p> <p>Background</p> <p>In <it>Mycobacterium tuberculosis </it>and in <it>Mycobacterium smegmatis </it>the <it>furA</it>-<it>katG </it>loci, encoding the FurA regulatory protein and the KatG catalase-peroxidase, are highly conserved. In <it>M. tuberculosis furA-katG </it>constitute a single operon, whereas in <it>M. smegmatis </it>a single mRNA covering both genes could not be found. In both species, specific 5' ends have been identified: the first one, located upstream of the <it>furA </it>gene, corresponds to transcription initiation from the <it>furA </it>promoter; the second one is the <it>katG </it>mRNA 5' end, located in the terminal part of <it>furA</it>.</p> <p>Results</p> <p>In this work we demonstrate by in vitro transcription and by RNA polymerase Chromatin immunoprecipitation that no promoter is present in the <it>M. smegmatis </it>region covering the latter 5' end, suggesting that it is produced by specific processing of longer transcripts. Several DNA fragments of <it>M. tuberculosis </it>and <it>M. smegmatis </it>were inserted in a plasmid between the <it>sigA </it>promoter and the <it>lacZ </it>reporter gene, and expression of the reporter gene was measured. A polypurine sequence, located four bp upstream of the <it>katG </it>translation start codon, increased beta-galactosidase activity and stabilized the <it>lacZ </it>transcript. Mutagenesis of this sequence led to destabilization of the mRNA. Analysis of constructs, in which the polypurine sequence of <it>M. smegmatis </it>was followed by an increasing number of <it>katG </it>codons, demonstrated that mRNA stability requires translation of at least 20 amino acids. In order to define the requirements for the 5' processing of the <it>katG </it>transcript, we created several mutations in this region and analyzed the 5' ends of the transcripts: the distance from the polypurine sequence does not seem to influence the processing, neither the sequence around the cutting point. Only mutations which create a double stranded region around the processing site prevented RNA processing.</p> <p>Conclusion</p> <p>This is the first reported case in mycobacteria, in which both a polypurine sequence and translation initiation are shown to contribute to mRNA stability. The <it>furA-katG </it>mRNA is transcribed from the <it>furA </it>promoter and immediately processed; this processing is prevented by a double stranded RNA at the cutting site, suggesting that the endoribonuclease responsible for the cleavage cuts single stranded RNA.</p

Characterization of a Novel Interaction between Bcl-2 Members Diva and Harakiri

Interactions within proteins of the Bcl-2 family are key in the regulation of apoptosis. The death-inducing members control apoptotic mechanisms partly by antagonizing the prosurvival proteins through heterodimer formation. Structural and biophysical studies on these complexes are providing important clues to understand their function. To help improve our knowledge on protein-protein interactions within the Bcl-2 family we have studied the binding between two of its members: mouse Diva and human Harakiri. Diva has been shown to perform both prosurvival and killing activity. In contrast, Harakiri induces cell death by interacting with antiapoptotic Bcl-2 members. Here we show using ELISA and NMR that Diva and Harakiri can interact in vitro. Combining the NMR data with the previously reported three-dimensional structure of Diva we find that Harakiri binds to a specific region in Diva. This interacting surface is equivalent to the known binding area of prosurvival Bcl-2 members from the reported structures of the complexes, suggesting that Diva could function at the structural level similarly to the antiapoptotic proteins of the Bcl-2 family. We illustrate this result by building a structural model of the heterodimer using molecular docking and the NMR data as restraints. Moreover, combining circular dichroism and NMR we also show that Harakiri is largely unstructured with residual (13%) α-helical conformation. This result agrees with intrinsic disorder previously observed in other Bcl-2 members. In addition, Harakiri constructs of different length were studied to identify the region critical for the interaction. Differential affinity for Diva of these constructs suggests that the amino acid sequence flanking the interacting region could play an important role in binding

Measurement of Lifetime and Decay-Width Difference in B0s -> J/psi phi Decays

We measure the mean lifetime, tau=2/(Gamma_L+Gamma_H), and the width difference, DeltaGamma=Gamma_L-Gamma_H, of the light and heavy mass eigenstates of the B0s meson, B0sL and B0sH, in B0s -> J/psi phi decays using 1.7 fb^-1 of data collected with the CDF II detector at the Fermilab Tevatron ppbar collider. Assuming CP conservation, a good approximation for the B0s system in the Standard Model, we obtain DeltaGamma = 0.076^+0.059_-0.063 (stat.) +- 0.006 (syst.) ps^-1 and tau = 1.52 +- 0.04 (stat.) +- 0.02 (syst.) ps, the most precise measurements to date. Our constraints on the weak phase and DeltaGamma are consistent with CP conservation. Dedicated to the memory of our dear friend and colleague, Michael P. Schmid

Limits on Anomalous Triple Gauge Couplings in ppbar Collisions at sqrt{s}=1.96 TeV

We present a search for anomalous triple gauge couplings (ATGC) in WW and WZ boson production. The boson pairs are produced in ppbar collisions at sqrt{s}=1.96 TeV, and the data sample corresponds to 350 pb-1 of integrated luminosity collected with the CDF II detector at the Fermilab Tevatron. In this search one W decays to leptons, and the other boson (W or Z) decays hadronically. Combining with a previously published CDF measurement of Wgamma boson production yields ATGC limits of -0.18 < lambda < 0.17 and -0.46 < Delta kappa < 0.39 at the 95% confidence level, using a cut-off scale Lambda=1.5 TeV.Comment: 7 pages, 3 figures. Submitted to Phys. Rev.

Search for Pair Production of Scalar Top Quarks Decaying to a tau Lepton and a b Quark in ppbar Collisions at sqrt{s}=1.96 TeV

We search for pair production of supersymmetric top quarks (~t_1), followed by R-parity violating decay ~t_1 -> tau b with a branching ratio beta, using 322 pb^-1 of ppbar collisions at sqrt{s}=1.96 TeV collected by the CDF II detector at Fermilab. Two candidate events pass our final selection criteria, consistent with the standard model expectation. We set upper limits on the cross section sigma(~t_1 ~tbar_1)*beta^2 as a function of the stop mass m(~t_1). Assuming beta=1, we set a 95% confidence level limit m(~t_1)>153 GeV/c^2. The limits are also applicable to the case of a third generation scalar leptoquark (LQ_3) decaying LQ_3 -> tau b.Comment: 7 pages, 2 eps figure

Studying the Underlying Event in Drell-Yan and High Transverse Momentum Jet Production at the Tevatron

We study the underlying event in proton-antiproton collisions by examining the behavior of charged particles (transverse momentum pT > 0.5 GeV/c, pseudorapidity |\eta| < 1) produced in association with large transverse momentum jets (~2.2 fb-1) or with Drell-Yan lepton-pairs (~2.7 fb-1) in the Z-boson mass region (70 < M(pair) < 110 GeV/c2) as measured by CDF at 1.96 TeV center-of-mass energy. We use the direction of the lepton-pair (in Drell-Yan production) or the leading jet (in high-pT jet production) in each event to define three regions of \eta-\phi space; toward, away, and transverse, where \phi is the azimuthal scattering angle. For Drell-Yan production (excluding the leptons) both the toward and transverse regions are very sensitive to the underlying event. In high-pT jet production the transverse region is very sensitive to the underlying event and is separated into a MAX and MIN transverse region, which helps separate the hard component (initial and final-state radiation) from the beam-beam remnant and multiple parton interaction components of the scattering. The data are corrected to the particle level to remove detector effects and are then compared with several QCD Monte-Carlo models. The goal of this analysis is to provide data that can be used to test and improve the QCD Monte-Carlo models of the underlying event that are used to simulate hadron-hadron collisions.Comment: Submitted to Phys.Rev.

Measurement of the W+W−W^+W^- Production Cross Section and Search for Anomalous WWγWW\gamma and WWZWWZ Couplings in ppˉp \bar p Collisions at s=1.96\sqrt{s} = 1.96 TeV

This Letter describes the current most precise measurement of the $W$ boson pair production cross section and most sensitive test of anomalous $WW\gamma$ and $WWZ$ couplings in $p \bar p$ collisions at a center-of-mass energy of 1.96 TeV. The $WW$ candidates are reconstructed from decays containing two charged leptons and two neutrinos, where the charged leptons are either electrons or muons. Using data collected by the CDF II detector from 3.6 fb$^{-1}$ of integrated luminosity, a total of 654 candidate events are observed with an expected background contribution of $320 \pm 47$ events. The measured total cross section is $\sigma (p \bar p \to W^+ W^- + X) = 12.1 \pm 0.9 \textrm{(stat)} ^{+1.6}_{-1.4} \textrm{(syst)}$ pb, which is in good agreement with the standard model prediction. The same data sample is used to place constraints on anomalous $WW\gamma$ and $WWZ$ couplings.Comment: submitted to Phys. Rev. Let

Observation of associated near-side and away-side long-range correlations in √sNN=5.02 TeV proton-lead collisions with the ATLAS detector

Two-particle correlations in relative azimuthal angle (Δϕ) and pseudorapidity (Δη) are measured in √sNN=5.02  TeV p+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using approximately 1  μb-1 of data as a function of transverse momentum (pT) and the transverse energy (ΣETPb) summed over 3.1<η<4.9 in the direction of the Pb beam. The correlation function, constructed from charged particles, exhibits a long-range (2<|Δη|<5) “near-side” (Δϕ∼0) correlation that grows rapidly with increasing ΣETPb. A long-range “away-side” (Δϕ∼π) correlation, obtained by subtracting the expected contributions from recoiling dijets and other sources estimated using events with small ΣETPb, is found to match the near-side correlation in magnitude, shape (in Δη and Δϕ) and ΣETPb dependence. The resultant Δϕ correlation is approximately symmetric about π/2, and is consistent with a dominant cos⁡2Δϕ modulation for all ΣETPb ranges and particle pT

Measurement of the Forward-Backward Asymmetry in the B -> K(*) mu+ mu- Decay and First Observation of the Bs -> phi mu+ mu- Decay

We reconstruct the rare decays $B^+ \to K^+\mu^+\mu^-$, $B^0 \to K^{*}(892)^0\mu^+\mu^-$, and $B^0_s \to \phi(1020)\mu^+\mu^-$ in a data sample corresponding to $4.4 {\rm fb^{-1}}$ collected in $p\bar{p}$ collisions at $\sqrt{s}=1.96 {\rm TeV}$ by the CDF II detector at the Fermilab Tevatron Collider. Using $121 \pm 16$ $B^+ \to K^+\mu^+\mu^-$ and $101 \pm 12$ $B^0 \to K^{*0}\mu^+\mu^-$ decays we report the branching ratios. In addition, we report the measurement of the differential branching ratio and the muon forward-backward asymmetry in the $B^+$ and $B^0$ decay modes, and the $K^{*0}$ longitudinal polarization in the $B^0$ decay mode with respect to the squared dimuon mass. These are consistent with the theoretical prediction from the standard model, and most recent determinations from other experiments and of comparable accuracy. We also report the first observation of the $B^0_s \to \phi\mu^+\mu^- decay and measure its branching ratio${\mathcal{B}}(B^0_s \to \phi\mu^+\mu^-) = [1.44 \pm 0.33 \pm 0.46] \times 10^{-6}$using$27 \pm 6$signal events. This is currently the most rare$B^0_s$ decay observed.Comment: 7 pages, 2 figures, 3 tables. Submitted to Phys. Rev. Let