Large introns in relation to alternative splicing and gene evolution: a case study of Drosophila bruno-3

Background: Alternative splicing (AS) of maturing mRNA can generate structurally and functionally distinct transcripts from the same gene. Recent bioinformatic analyses of available genome databases inferred a positive correlation between intron length and AS. To study the interplay between intron length and AS empirically and in more detail, we analyzed the diversity of alternatively spliced transcripts (ASTs) in the Drosophila RNA-binding Bruno-3 (Bru-3) gene. This gene was known to encode thirteen exons separated by introns of diverse sizes, ranging from 71 to 41,973 nucleotides in D. melanogaster. Although Bru-3's structure is expected to be conducive to AS, only two ASTs of this gene were previously described. Results: Cloning of RT-PCR products of the entire ORF from four species representing three diverged Drosophila lineages provided an evolutionary perspective, high sensitivity, and long-range contiguity of splice choices currently unattainable by high-throughput methods. Consequently, we identified three new exons, a new exon fragment and thirty-three previously unknown ASTs of Bru-3. All exon-skipping events in the gene were mapped to the exons surrounded by introns of at least 800 nucleotides, whereas exons split by introns of less than 250 nucleotides were always spliced contiguously in mRNA. Cases of exon loss and creation during Bru-3 evolution in Drosophila were also localized within large introns. Notably, we identified a true de novo exon gain: exon 8 was created along the lineage of the obscura group from intronic sequence between cryptic splice sites conserved among all Drosophila species surveyed. Exon 8 was included in mature mRNA by the species representing all the major branches of the obscura group. To our knowledge, the origin of exon 8 is the first documented case of exonization of intronic sequence outside vertebrates. Conclusion: We found that large introns can promote AS via exon-skipping and exon turnover during evolution likely due to frequent errors in their removal from maturing mRNA. Large introns could be a reservoir of genetic diversity, because they have a greater number of mutable sites than short introns. Taken together, gene structure can constrain and/or promote gene evolution

Interleaved Parton Showers and Tuning Prospects

General-purpose Monte Carlo event generators have become important tools in particle physics, allowing the simulation of exclusive hadronic final states. In this article we examine the Pythia 8 generator, in particular focusing on its parton-shower algorithms. Some relevant new additions to the code are introduced, that should allow for a better description of data. We also implement and compare with 2 to 3 real-emission QCD matrix elements, to check how well the shower algorithm fills the phase space away from the soft and collinear regions. A tuning of the generator to Tevatron data is performed for two PDF sets and the impact of first new LHC data is examined

Exclusive Leptoproduction of rho^0 Mesons from Hydrogen at Intermediate Virtual Photon Energies

Measurements of the cross section for exclusive virtual-photoproduction of rho^0 mesons from hydrogen are reported. The data were collected by the HERMES experiment using 27.5 GeV positrons incident on a hydrogen gas target in the HERA storage ring. The invariant mass W of the photon-nucleon system ranges from 4.0 to 6.0 GeV, while the negative squared four-momentum Q^2 of the virtual photon varies from 0.7 to 5.0 GeV^2. The present data together with most of the previous data at W > 4 GeV are well described by a model that infers the W-dependence of the cross section from the dependence on the Bjorken scaling variable x of the unpolarized structure function for deep-inelastic scattering. In addition, a model calculation based on Off-Forward Parton Distributions gives a fairly good account of the longitudinal component of the rho^0 production cross section for Q^2 > 2 GeV^2.Comment: 10 pages, 6 embedded figures, LaTeX for SVJour(epj) document class. Revisions: curves added to Fig. 1, several clarifications added to tex

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal

Jet energy measurement with the ATLAS detector in proton-proton collisions at root s=7 TeV

The jet energy scale and its systematic uncertainty are determined for jets measured with the ATLAS detector at the LHC in proton-proton collision data at a centre-of-mass energy of √s = 7TeV corresponding to an integrated luminosity of 38 pb-1. Jets are reconstructed with the anti-kt algorithm with distance parameters R=0. 4 or R=0. 6. Jet energy and angle corrections are determined from Monte Carlo simulations to calibrate jets with transverse momenta pT≥20 GeV and pseudorapidities {pipe}η{pipe}<4. 5. The jet energy systematic uncertainty is estimated using the single isolated hadron response measured in situ and in test-beams, exploiting the transverse momentum balance between central and forward jets in events with dijet topologies and studying systematic variations in Monte Carlo simulations. The jet energy uncertainty is less than 2. 5 % in the central calorimeter region ({pipe}η{pipe}<0. 8) for jets with 60≤pT<800 GeV, and is maximally 14 % for pT<30 GeV in the most forward region 3. 2≤{pipe}η{pipe}<4. 5. The jet energy is validated for jet transverse momenta up to 1 TeV to the level of a few percent using several in situ techniques by comparing a well-known reference such as the recoiling photon pT, the sum of the transverse momenta of tracks associated to the jet, or a system of low-pT jets recoiling against a high-pT jet. More sophisticated jet calibration schemes are presented based on calorimeter cell energy density weighting or hadronic properties of jets, aiming for an improved jet energy resolution and a reduced flavour dependence of the jet response. The systematic uncertainty of the jet energy determined from a combination of in situ techniques is consistent with the one derived from single hadron response measurements over a wide kinematic range. The nominal corrections and uncertainties are derived for isolated jets in an inclusive sample of high-pT jets. Special cases such as event topologies with close-by jets, or selections of samples with an enhanced content of jets originating from light quarks, heavy quarks or gluons are also discussed and the corresponding uncertainties are determined. © 2013 CERN for the benefit of the ATLAS collaboration

Hypothesis for the evolution of three-helix Chl a/b and Chl a/c light-harvesting antenna proteins from two-helix and four-helix ancestors

The nuclear-encoded Chl a/b and Chl a/c antenna proteins of photosynthetic eukaryotes are part of an extended family of proteins that also includes the early light-induced proteins (ELIPs) and the 22 kDa intrinsic protein of PS II (encoded by psb S gene). All members of this family have three transmembrane helices except for the psb S protein, which has four. The amino acid sequences of these proteins are compared and related to the three-dimensional structure of pea LHC II Type I (Kühlbrandt and Wang, Nature 350: 130–134, 1991). The similarity of psb S to the three-helix members of the family suggests that the latter arose from a four-helix ancestor that lost its C-terminal helix by deletion. Strong internal similarity between the two halves of the psb S protein suggests that it in turn arose as the result of the duplication of a gene encoding a two-helix protein. Since psb S is reported to be present in at least one cyanobacterium, the ancestral four-helix protein may have been present prior to the endosymbiotic event or events that gave rise to the photosynthetic eukaryotes. The Chl a/b and Chl a/c antenna proteins, and the immunologically-related proteins in the rhodophytes may have had a common ancestor which was present in the early photosynthetic eukaryotes, and predated their division into rhodophyte, chromophyte and chlorophyte lineages. The LHC I-LHC II divergence probably occurred before the separation of higher plants from chlorophyte algae and euglenophytes, and the different Types of LHC I and LHC II proteins arose prior to the separation of angiosperms and gymnosperms.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/43538/1/11120_2004_Article_BF00029382.pd

The MDM2-p53 pathway is involved in preconditioning-induced neuronal tolerance to ischemia.

Brain preconditioning (PC) refers to a state of transient tolerance against a lethal insult that can be evoked by a prior mild event. It is thought that PC may induce different pathways responsible for neuroprotection, which may involve the attenuation of cell damage pathways, including the apoptotic cell death. In this context, p53 is a stress sensor that accumulates during brain ischemia leading to neuronal death. The murine double minute 2 gene (MDM2), a p53-specific E3 ubiquitin ligase, is the main cellular antagonist of p53, mediating its degradation by the proteasome. Here, we study the role of MDM2-p53 pathway on PC-induced neuroprotection both in cultured neurons (in vitro) and rat brain (in vivo). Our results show that PC increased neuronal MDM2 protein levels, which prevented ischemiainduced p53 stabilization and neuronal death. Indeed, PC attenuated ischemia-induced activation of the p53/PUMA/caspase-3 signaling pathway. Pharmacological inhibition of MDM2-p53 interaction in neurons abrogated PC-induced neuroprotection against ischemia. Finally, the relevance of the MDM2-p53 pathway was confirmed in rat brain using a PC model in vivo. These findings demonstrate the key role of the MDM2-p53 pathway in PC-induced neuroprotection against a subsequent ischemic insult and poses MDM2 as an essential target in ischemic tolerance

Active Nuclear Receptors Exhibit Highly Correlated AF-2 Domain Motions

Nuclear receptor ligand binding domains (LBDs) convert ligand binding events into changes in gene expression by recruiting transcriptional coregulators to a conserved activation function-2 (AF-2) surface. While most nuclear receptor LBDs form homo- or heterodimers, the human nuclear receptor pregnane X receptor (PXR) forms a unique and essential homodimer and is proposed to assemble into a functional heterotetramer with the retinoid X receptor (RXR). How the homodimer interface, which is located 30 Å from the AF-2, would affect function at this critical surface has remained unclear. By using 20- to 30-ns molecular dynamics simulations on PXR in various oligomerization states, we observed a remarkably high degree of correlated motion in the PXR–RXR heterotetramer, most notably in the four helices that create the AF-2 domain. The function of such correlation may be to create “active-capable” receptor complexes that are ready to bind to transcriptional coactivators. Indeed, we found in additional simulations that active-capable receptor complexes involving other orphan or steroid nuclear receptors also exhibit highly correlated AF-2 domain motions. We further propose a mechanism for the transmission of long-range motions through the nuclear receptor LBD to the AF-2 surface. Taken together, our findings indicate that long-range motions within the LBD scaffold are critical to nuclear receptor function by promoting a mobile AF-2 state ready to bind coactivators

Characterizing Structural Transitions Using Localized Free Energy Landscape Analysis

Structural changes in molecules are frequently observed during biological processes like replication, transcription and translation. These structural changes can usually be traced to specific distortions in the backbones of the macromolecules involved. Quantitative energetic characterization of such distortions can greatly advance the atomic-level understanding of the dynamic character of these biological processes.Molecular dynamics simulations combined with a variation of the Weighted Histogram Analysis Method for potential of mean force determination are applied to characterize localized structural changes for the test case of cytosine (underlined) base flipping in a GTCAGCGCATGG DNA duplex. Free energy landscapes for backbone torsion and sugar pucker degrees of freedom in the DNA are used to understand their behavior in response to the base flipping perturbation. By simplifying the base flipping structural change into a two-state model, a free energy difference of upto 14 kcal/mol can be attributed to the flipped state relative to the stacked Watson-Crick base paired state. This two-state classification allows precise evaluation of the effect of base flipping on local backbone degrees of freedom.The calculated free energy landscapes of individual backbone and sugar degrees of freedom expectedly show the greatest change in the vicinity of the flipping base itself, but specific delocalized effects can be discerned upto four nucleotide positions away in both 5' and 3' directions. Free energy landscape analysis thus provides a quantitative method to pinpoint the determinants of structural change on the atomic scale and also delineate the extent of propagation of the perturbation along the molecule. In addition to nucleic acids, this methodology is anticipated to be useful for studying conformational changes in all macromolecules, including carbohydrates, lipids, and proteins

A Novel Side-Chain Orientation Dependent Potential Derived from Random-Walk Reference State for Protein Fold Selection and Structure Prediction

An accurate potential function is essential to attack protein folding and structure prediction problems. The key to developing efficient knowledge-based potential functions is to design reference states that can appropriately counteract generic interactions. The reference states of many knowledge-based distance-dependent atomic potential functions were derived from non-interacting particles such as ideal gas, however, which ignored the inherent sequence connectivity and entropic elasticity of proteins.We developed a new pair-wise distance-dependent, atomic statistical potential function (RW), using an ideal random-walk chain as reference state, which was optimized on CASP models and then benchmarked on nine structural decoy sets. Second, we incorporated a new side-chain orientation-dependent energy term into RW (RWplus) and found that the side-chain packing orientation specificity can further improve the decoy recognition ability of the statistical potential.RW and RWplus demonstrate a significantly better ability than the best performing pair-wise distance-dependent atomic potential functions in both native and near-native model selections. It has higher energy-RMSD and energy-TM-score correlations compared with other potentials of the same type in real-life structure assembly decoys. When benchmarked with a comprehensive list of publicly available potentials, RW and RWplus shows comparable performance to the state-of-the-art scoring functions, including those combining terms from multiple resources. These data demonstrate the usefulness of random-walk chain as reference states which correctly account for sequence connectivity and entropic elasticity of proteins. It shows potential usefulness in structure recognition and protein folding simulations. The RW and RWplus potentials, as well as the newly generated I-TASSER decoys, are freely available in http://zhanglab.ccmb.med.umich.edu/RW