基于贝叶斯网络的DNA序列剪接位点预测

采用基于贝叶斯网络的建模方法,预测真核生物DNA序列中的剪接位点．分别建立了供体位点和受体位点模型,并根据两种位点的生物学特性,对模型的拓扑结构和上下游节点的选择进行了优化．通过贝叶斯网络的最大似然学习算法求出模型参数后,利用10分组交互验证方法对测试数据进行剪接位点预测。结果显示,受体位点的平均预测准确率为92．5％,伪受体位点的平均预测准确率为94．0％,供体位点的平均预测准确率为92．3％,伪供体位点的平均预测准确率为93．5％,整体效果要好于基于使用独立和条件概率矩阵、以及隐Markov模型的预测方法．表明利用贝叶斯网络对剪接位点建模是预测剪接位点的一种有效手段．     

A statistical framework for genomic data fusion

MOTIVATION: During the past decade, the new focus on genomics has highlighted a particular challenge: to integrate the different views of the genome that are provided by various types of experimental data. RESULTS: This paper describes a computational framework for integrating and drawing inferences from a collection of genome-wide measurements. Each dataset is represented via a kernel function, which defines generalized similarity relationships between pairs of entities, such as genes or proteins. The kernel representation is both flexible and efficient, and can be applied to many different types of data. Furthermore, kernel functions derived from different types of data can be combined in a straightforward fashion. Recent advances in the theory of kernel methods have provided efficient algorithms to perform such combinations in a way that minimizes a statistical loss function. These methods exploit semidefinite programming techniques to reduce the problem of finding optimizing kernel combinations to a convex optimization problem. Computational experiments performed using yeast genome-wide datasets, including amino acid sequences, hydropathy profiles, gene expression data and known protein-protein interactions, demonstrate the utility of this approach. A statistical learning algorithm trained from all of these data to recognize particular classes of proteins--membrane proteins and ribosomal proteins--performs significantly better than the same algorithm trained on any single type of data. AVAILABILITY: Supplementary data at http://noble.gs.washington.edu/proj/sdp-svm     

Recombination Losses Above and Below the Transport Percolation Threshold in Bulk Heterojunction Organic Solar Cells

Achieving the highest power conversion efficiencies in bulk heterojunction organic solar cells requires a morphology that delivers electron and hole percolation pathways for optimized transport, plus sufficient donor:acceptor contact area for near unity charge transfer state formation. This is a significant structural challenge, particularly in semiconducting polymer:fullerene systems. This balancing act in the model high efficiency PTB7:PC70BM blend is studied by tuning the donor:acceptor ratio, with a view to understanding the recombination loss mechanisms above and below the fullerene transport percolation threshold. The internal quantum efficiency is found to be strongly correlated to the slower carrier mobility in agreement with other recent studies. Furthermore, second‐order recombination losses dominate the shape of the current density–voltage curve in efficient blend combinations, where the fullerene phase is percolated. However, below the charge transport percolation threshold, there is an electric‐field dependence of first‐order losses, which includes electric‐field‐dependent photogeneration. In the intermediate regime, the fill factor appears to be limited by both first‐ and second‐order losses. These findings provide additional basic understanding of the interplay between the bulk heterojunction morphology and the order of recombination in organic solar cells. They also shed light on the limitations of widely used transport models below the percolation threshold.     

Charge Transport in Pure and Mixed Phases in Organic Solar Cells

In organic solar cells continuous donor and acceptor networks are considered necessary for charge extraction, whereas discontinuous neat phases and molecularly mixed donor–acceptor phases are generally regarded as detrimental. However, the impact of different levels of domain continuity, purity, and donor–acceptor mixing on charge transport remains only semiquantitatively described. Here, cosublimed donor–acceptor mixtures, where the distance between the donor sites is varied in a controlled manner from homogeneously diluted donor sites to a continuous donor network are studied. Using transient measurements, spanning from sub‐picoseconds to microseconds photogenerated charge motion is measured in complete photovoltaic devices, to show that even highly diluted donor sites (5.7%–10% molar) in a buckminsterfullerene matrix enable hole transport. Hopping between isolated donor sites can occur by long‐range hole tunneling through several buckminsterfullerene molecules, over distances of up to ≈4 nm. Hence, these results question the relevance of “pristine” phases and whether a continuous interpenetrating donor–acceptor network is the ideal morphology for charge transport.     

Interactions between pairs of DNA-binding dyes: results and implications of chromosome analysis

A number of DNA-binding dyes, with spectral properties making them suitable as components of energy donor-acceptor pairs, are described. If such pairs are used to stain metaphase chromosomes, and if the energy acceptor (e.g., actinomycin D or methyl green) has a binding specificity opposite to the binding or fluorescence specificity of the donor (e.g., 33258 Hoechst, quinacrine or chromomycin A3), contrast in donor fluorescence can be enhanced, leading to patterns selectively highlighting standard or reverse chromosome bands or particular polymorphic regions. Such results presumably reflect chromosomal regions enriched in 10-20 base pair clusters to which the donor binds and fluoresces but to which the acceptor cannot bind. For other pairs, involving counterstains such as netropsin or echinomycin, which are not suitable as energy acceptors, specific changes observed in polymorphic region fluorescence are most likely due to binding competition between dyes. Dye pairs producing contrast by either method can be used to differentiate between homologous chromosomes or to facilitate detection of specific chromosomal rearrangements. Preliminary data indicate that contrast enhancement generated in fixed metaphase chromosomes spread on microscopic slides can also be observed in suspensions of unfixed metaphase chromosomes, reinforcing the expectation that the methodology described will be of use in flow cytometry.     

Donor functionalized quinoline based organic sensitizers for dye sensitized solar cell (DSSC) applications: DFT and TD-DFT investigations

The influence of different donor groups in quinoline based novel sensitizers for dye sensitized solar cell (DSSC) applications is analyzed by using density functional theory (DFT) and time dependent density functional theory (TD-DFT). Quinoline and donor functionalized quinoline based novel organic sensitizers have been designed with different π-spacers for DSSC applications. The ground state molecular structure of novel organic sensitizers is fully optimized by DFT calculation in both gas and chloroform phases. Electronic absorption characteristics are predicted by the TD-DFT calculation in both gas and chloroform phases. The polarizable continuum model is used for solvent phase optimization. The net electron transfer from the donor to acceptor is calculated from natural bond orbital (NBO) analysis. The injection energy and dye regeneration energy values are also calculated. Different donor groups are substituted in quinoline, and these substituted quinoline donors are used as the donor group. Cyanovinyl and thiophene groups act as π-spacers and cyanoacrylic acid acts as an acceptor. DFT and TD-DFT studies of the quinoline and donor functionalized quinoline sensitizers show that the coumarin based and N-hexyltetrahydroquinoline donors are more efficient for DSSC application.     

On the Molecular Origin of Charge Separation at the Donor–Acceptor Interface

Fullerene‐based acceptors have dominated organic solar cells for almost two decades. It is only within the last few years that alternative acceptors rival their dominance, introducing much more flexibility in the optoelectronic properties of these material blends. However, a fundamental physical understanding of the processes that drive charge separation at organic heterojunctions is still missing, but urgently needed to direct further material improvements. Here a combined experimental and theoretical approach is used to understand the intimate mechanisms by which molecular structure contributes to exciton dissociation, charge separation, and charge recombination at the donor–acceptor (D–A) interface. Model systems comprised of polythiophene‐based donor and rylene diimide‐based acceptor polymers are used and a detailed density functional theory (DFT) investigation is performed. The results point to the roles that geometric deformations and direct‐contact intermolecular polarization play in establishing a driving force (energy gradient) for the optoelectronic processes taking place at the interface. A substantial impact for this driving force is found to stem from polymer deformations at the interface, a finding that can clearly lead to new design approaches in the development of the next generation of conjugated polymers and small molecules.     

Pyrimidone-based series of glucokinase activators with alternative donor–acceptor motif

Glucokinase activators are a class of experimental agents under investigation as a therapy for Type 2 diabetes mellitus. An X-ray crystal structure of a modestly potent agent revealed the potential to substitute the common heterocyclic amide donor–acceptor motif for a pyridone moiety. We have successfully demonstrated that both pyridone and pyrimidone heterocycles can be used as a potent donor–acceptor substituent. Several sub-micromolar analogs that possess the desired partial activator profile were synthesized and characterized. Unfortunately, the most potent activators suffered from sub-optimal pharmacokinetic properties. Nonetheless, these donor–acceptor motifs may find utility in other glucokinase activator series or beyond.     

RNA canonical and non-canonical base pairing types: a recognition method and complete repertoire

The problem of systematic and objective identification of canonical and non-canonical base pairs in RNA three-dimensional (3D) structures was studied. A probabilistic approach was applied, and an algorithm and its implementation in a computer program that detects and analyzes all the base pairs contained in RNA 3D structures were developed. The algorithm objectively distinguishes among canonical and non-canonical base pairing types formed by three, two and one hydrogen bonds (H-bonds), as well as those containing bifurcated and C-H...X H-bonds. The nodes of a bipartite graph are used to encode the donor and acceptor atoms of a 3D structure. The capacities of the edges correspond to probabilities computed from the geometry of the donor and acceptor groups to form H-bonds. The maximum flow from donors to acceptors irectly identifies base pairs and their types. A complete repertoire of base pairing types was built from the detected H-bonds of all X-ray crystal structures of a resolution of 3.0 Å or better, including the large and small ribosomal subunits. The base pairing types are labeled using an extension of the nomenclature recently introduced by Leontis and Westhof. The probabilistic method was implemented in MC-Annotate, an RNA structure analysis computer program used to determine the base pairing parameters of the 3D modeling system MC-Sym.     

Prediction of human mRNA donor and acceptor sites from the DNA sequence

Artificial neural networks have been applied to the prediction of splice site location in human pre-mRNA. A joint prediction scheme where prediction of transition regions between introns and exons regulates a cutoff level for splice site assignment was able to predict splice site locations with confidence levels far better than previously reported in the literature. The problem of predicting donor and acceptor sites in human genes is hampered by the presence of numerous amounts of false positives: here, the distribution of these false splice sites is examined and linked to a possible scenario for the splicing mechanism in vivo. When the presented method detects 95% of the true donor and acceptor sites, it makes less than 0.1% false donor site assignments and less than 0.4% false acceptor site assignments. For the large data set used in this study, this means that on average there are one and a half false donor sites per true donor site and six false acceptor sites per true acceptor site. With the joint assignment method, more than a fifth of the true donor sites and around one fourth of the true acceptor sites could be detected without accompaniment of any false positive predictions. Highly confident splice sites could not be isolated with a widely used weight matrix method or by separate splice site networks. A complementary relation between the confidence levels of the coding/non-coding and the separate splice site networks was observed, with many weak splice sites having sharp transitions in the coding/non-coding signal and many stronger splice sites having more ill-defined transitions between coding and non-coding.     

Dynamic map of protein interactions in the Escherichia coli chemotaxis pathway

Protein–protein interactions play key roles in virtually all cellular processes, often forming complex regulatory networks. A powerful tool to study interactions in vivo is fluorescence resonance energy transfer (FRET), which is based on the distance‐dependent energy transfer from an excited donor to an acceptor fluorophore. Here, we used FRET to systematically map all protein interactions in the chemotaxis signaling pathway in Escherichia coli, one of the most studied models of signal transduction, and to determine stimulation‐induced changes in the pathway. Our FRET analysis identified 19 positive FRET pairs out of the 28 possible protein combinations, with 9 pairs being responsive to chemotactic stimulation. Six stimulation‐dependent and five stimulation‐independent interactions were direct, whereas other interactions were apparently mediated by scaffolding proteins. Characterization of stimulation‐induced responses revealed an additional regulation through activity dependence of interactions involving the adaptation enzyme CheB, and showed complex rearrangement of chemosensory receptors. Our study illustrates how FRET can be efficiently employed to study dynamic protein networks in vivo.     

Sequentially Deposited versus Conventional Nonfullerene Organic Solar Cells: Interfacial Trap States,Vertical Stratification,and Exciton Dissociation

Bulk heterojunction (BHJ) nonfullerene organic solar cells prepared from sequentially deposited donor and acceptor layers (sq‐BHJ) have recently been shown to be highly efficient, environmentally friendly, and compatible with large area and roll‐to‐roll fabrication. However, the related photophysics at donor‐acceptor interface and the vertical heterogeneity of donor‐acceptor distribution, critical for exciton dissociation and device performance, have been largely unexplored. Herein, steady‐state and time‐resolved optical and electrical techniques are employed to characterize the interfacial trap states. Correlating with the luminescent efficiency of interfacial states and its nonradiative recombination, interfacial trap states are characterized to be about 40% more populated in the sq‐BHJ devices than the as‐cast BHJ (c‐BHJ), which probably limits the device voltage output. Cross‐sectional energy‐dispersive X‐ray spectroscopy and ultraviolet photoemission spectroscopy depth profiling directly visualize the donor–acceptor vertical stratification with a precision of 1–2 nm. From the proposed “needle” model, the high exciton dissociation efficiency is rationalized. This study highlights the promise of sequential deposition to fabricate efficient solar cells, and points toward improving the voltage output and overall device performance via eliminating interfacial trap states.     

基于RNA-Seq数据识别果蝇剪接位点和可变剪接事件

完整基因结构的预测是当前生命科学研究的一个重要基础课题,其中一个关键环节是剪接位点和各种可变剪接事件的精确识别.基于转录组测序（RNA-seq）数据,识别剪接位点和可变剪接事件是近几年随着新一代测序技术发展起来的新技术策略和方法.本工作基于黑腹果蝇睾丸RNA-seq数据,使用TopHat软件成功识别出39718个果蝇剪接位点,其中有10584个新剪接位点.同时,基于剪接位点的不同组合,针对各类型可变剪接特征开发出计算识别算法,成功识别了8477个可变剪接事件（其中新识别的可变剪接事件3922个）,包括可变供体位点、可变受体位点、内含子保留和外显子缺失4种类型.RT-PCR实验验证了2个果蝇基因上新识别的可变剪接事件,发现了全新的剪接异构体.进一步表明,RNA-seq数据可有效应用于识别剪接位点和可变剪接事件,为深入揭示剪接机制及可变剪接生物学功能提供新思路和新手段.     

On the possibility of long-wavelength long-lifetime high-quantum-yield luminophores

We describe an approach to creating a new class of luminophores which display both long wavelength emissions exceeding 600 nm and long lifetimes. These luminophores are based on resonance energy transfer (RET) from a long lifetime donor to a short lifetime but long wavelength acceptor. We demonstrated the possibility of obtaining these desirable spectral properties using donors and acceptors noncovalently bound to DNA. The donor was a ruthenium (Ru) metal-ligand complex in which one of the diimine ligands intercalated into double-helix DNA. The acceptors were either nile blue, TOTO-3, or TO-PRO-3. Upon binding of the acceptor to donor-labeled DNA, we found that the acceptor quantum yield was remarkably enhanced so that the wavelength-integrated intensities of the donor and acceptor bound to DNA were many-fold greater than the intensity of the donor and acceptor alone when separately bound to DNA. The origin of this effect is efficient energy transfer from the donor. Under these conditions the effective overall quantum yield approaches that of the acceptor. Importantly, the increased quantum yield can be obtained while maintaining usefully long apparent acceptor lifetimes of 30 to 80 ns. The effect of an increased quantum yield from a low quantum yield donor may find use in assays to detect macromolecular binding interactions. These results suggest the synthesis of covalently linked donor-acceptor pairs with the desirable spectral properties of long wavelength emission, high quantum yield, and moderately long lifetimes for gated detection.     

Structural changes of yellow Cameleon domains observed by quantitative FRET analysis and polarized fluorescence correlation spectroscopy

Förster resonance energy transfer (FRET) is a widely used method for monitoring interactions between or within biological macromolecules conjugated with suitable donor-acceptor pairs. Donor fluorescence lifetimes in absence and presence of acceptor molecules are often measured for the observation of FRET. However, these lifetimes may originate from interacting and noninteracting molecules, which hampers quantitative interpretation of FRET data. We describe a methodology for the detection of FRET that monitors the rise time of acceptor fluorescence on donor excitation thereby detecting only those molecules undergoing FRET. The large advantage of this method, as compared to donor fluorescence quenching method used more commonly, is that the transfer rate of FRET can be determined accurately even in cases where the FRET efficiencies approach 100% yielding highly quenched donor fluorescence. Subsequently, the relative orientation between donor and acceptor chromophores is obtained from time-dependent fluorescence anisotropy measurements carried out under identical conditions of donor excitation and acceptor detection. The FRET based calcium sensor Yellow Cameleon 3.60 (YC3.60) was used because it changes its conformation on calcium binding, thereby increasing the FRET efficiency. After mapping distances and orientation angles between the FRET moieties in YC3.60, cartoon models of this FRET sensor with and without calcium could be created. Independent support for these representations came from experiments where the hydrodynamic properties of YC3.60 under ensemble and single-molecule conditions on selective excitation of the acceptor were determined. From rotational diffusion times as found by fluorescence correlation spectroscopy and consistently by fluorescence anisotropy decay analysis it could be concluded that the open structure (without calcium) is flexible as opposed to the rather rigid closed conformation. The combination of two independent methods gives consistent results and presents a rapid and specific methodology to analyze structural and dynamical changes in a protein on ligand binding.     

Engineering resonance energy transfer for advanced immunoassays: the case of celiac disease

Celiac disease (CD) is an immune-mediated disorder affecting genetically predisposed subjects. It is caused by the ingestion of wheat gluten and related prolamins. A final diagnosis for this disease can be obtained by examination of jejunal biopsies. Nevertheless, different analytical approaches have been established to detect the presence of anti-tissue transglutaminase antibodies that represent a serological hallmark of the disease. In this work, we explored a new method for the diagnosis of CD based on the detection of serum anti-transglutaminase antibodies by resonance energy transfer (RET) between donor molecules and acceptor molecules. In particular, we labeled the liver transglutaminase (tTG) enzyme from guinea pig and the rabbit anti-tTG antibodies with a couple of fluorescence probes that are able to make RET if they are located within with Förster distance. We labeled tTG with the fluorescence probe DyLight 594 as donor and the anti-tTG antibodies with the fluorescence probe DyLight 649 as acceptor. However, due to the large size of the formed complex (tTG/anti-tTG), and consequently to the low efficiency energy transfer process between the donor–acceptor molecules, we explored a new experimental approach that allows us to extend the utilizable range of RET between donor:acceptor pairs by using one single molecule as donor and multiple molecules as energy acceptors, instead of using a single acceptor molecule as usually occurs in RET experiments. The obtained results clearly show that the use of one donor and multiacceptor strategy enables for a simple and rapid detection of serum anti-transglutaminase antibodies. In addition, our results point out that it is possible to consider this approach as a new method for a wide variety of analytical assays.     

Hologram quantitative structure-activity relationships for a series of farnesoid X receptor activators

The farnesoid X receptor (FXR) is an attractive drug target for the development of novel therapeutic agents for the treatment of dyslipidemia and cholestasis. Hologram quantitative structure-activity relationship (HQSAR) studies were conducted on a series of potent FXR activators originated from natural product-like libraries. A training set containing 82 compounds served to establish the models. The best HQSAR model was generated using atoms, bonds, connections, chirality, and donor and acceptor as fragment distinction and fragment size default (4-7) with six components. The model was used to predict the potency of 20 test set compounds that were not included in the training set, and the predicted values were in good agreement with the experimental results. The final HQSAR model and the information obtained from HQSAR 2D contribution maps should be useful for the design of novel FXR ligands having improved potency.     

A new semiempirical codon substitution model based on principal component analysis of mammalian sequences

Codon substitution models have traditionally been parametric Markov models, but recently, empirical and semiempirical models also have been proposed. Parametric codon models are typically based on 61×61 rate matrices that are derived from a small number of parameters. These parameters are rooted in experience and theoretical considerations and generally show good performance but are still relatively arbitrary. We have previously used principal component analysis (PCA) on data obtained from mammalian sequence alignments to empirically identify the most relevant parameters for codon substitution models, thereby confirming some commonly used parameters but also suggesting new ones. Here, we present a new semiempirical codon substitution model that is directly based on those PCA results. The substitution rate matrix is constructed from linear combinations of the first few (the most important) principal components with the coefficients being free model parameters. Thus, the model is not only based on empirical rates but also uses the empirically determined most relevant parameters for a codon model to adjust to the particularities of individual data sets. In comparisons against established parametric and semiempirical models, the new model consistently achieves the highest likelihood values when applied to sequences of vertebrates, which include the taxonomic class where the model was trained on.     

Prediction of the lowest charge-transfer excited-state energy at the donor–acceptor interface in a condensed phase using ground-state DFT calculations with generalized Kohn–Sham functionals

Knowledge of the charge (electron) transfer process at the donor–acceptor interface is required to understand the working mechanisms of different organic photovoltaic materials. Investigating the lowest charge-transfer state in organic donor–acceptor solar cells is important as it allows the destruction/formation of excitons and polarons to be studied, and is directly related to the open circuit voltage. By performing low-cost and feasible calculations of ground-state electronic structures using the Mulliken rule as well as the optimally tuned range-separated hybrid (OTRSH) density functional and a regular long-range corrected functional, the lowest charge-transfer (CT) state energies of a series of dimers containing functionalized anthracene (the donor) and tetracyanoethylene (the acceptor) were obtained. The jumping distances of excited electrons during CT were calculated. The polarizable continuum model was applied to account for the effects of the solvent methylene chloride (CH₂Cl₂) on the lowest CT state energies obtained from gas-phase calculations. The calculated lowest CT state energies of the dimers were close to the corresponding experimental results, with a root mean square deviation (RMSD) of 0.22 eV.