baerhunter: an R package for the discovery and analysis of expressed non-coding regions in bacterial RNA-seq data.

SUMMARY: Standard bioinformatics pipelines for the analysis of bacterial transcriptomic data commonly ignore non-coding but functional elements e.g. small RNAs, long antisense RNAs or untranslated regions (UTRs) of mRNA transcripts. The root of this problem is the use of incomplete genome annotation files. Here, we present baerhunter, a coverage-based method implemented in R, that automates the discovery of expressed non-coding RNAs and UTRs from RNA-seq reads mapped to a reference genome. The core algorithm is part of a pipeline that facilitates downstream analysis of both coding and non-coding features. The method is simple, easy to extend and customize and, in limited tests with simulated and real data, compares favourably against the currently most popular alternative. AVAILABILITY AND IMPLEMENTATION: The baerhunter R package is available from: https://github.com/irilenia/baerhunter. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Molecular dynamics simulations of the interaction of wild type and mutant human CYP2J2 with polyunsaturated fatty acids.

OBJECTIVES: The data presented here is part of a study that was aimed at characterizing the molecular mechanisms of polyunsaturated fatty acid metabolism by CYP2J2, the main cytochrome P450 enzyme active in the human cardiovasculature. This part comprises the molecular dynamics simulations of the binding of three eicosanoid substrates to wild type and mutant forms of the enzyme. These simulations were carried out with the aim of dissecting the importance of individual residues in the active site and the roles they might play in dictating the binding and catalytic specificity exhibited by CYP2J2. DATA DESCRIPTION: The data comprise: (a) a new homology model of CYP2J2, (b) a number of predicted low-energy complexes of CYP2J2 with arachidonic acid, docosahexaenoic acid and eicosapentaenoic acid, produced with molecular docking and (c) a series of molecular dynamics simulations of the wild type and four mutants interacting with arachidonic acid as well as simulations of the wild type interacting with the two other eicosanoid ligands. The simulations may be helpful in identifying the determinants of substrate specificity of this enzyme and in unraveling the role of individual mutations on its function. They may also help guide the generation of mutants with altered substrate preferences.

Computational modelling of the binding of arachidonic acid to the human monooxygenase CYP2J2.

An experimentally determined structure for human CYP2J2-a member of the cytochrome P450 family with significant and diverse roles across a number of tissues-does not yet exist. Our understanding of how CYP2J2 accommodates its cognate substrates and how it might be inhibited by other ligands thus relies on our ability to computationally predict such interactions using modelling techniques. In this study we present a computational investigation of the binding of arachidonic acid (AA) to CYP2J2 using homology modelling, induced fit docking (IFD) and molecular dynamics (MD) simulations. Our study reveals a catalytically competent binding mode for AA that is distinct from a recently published study that followed a different computational pipeline. Our proposed binding mode for AA is supported by crystal structures of complexes of related enzymes to inhibitors, and evolutionary conservation of a residue whose role appears essential for placing AA in the right site for catalysis. Graphical Abstract Arachidonic acid docked in the active site of CYP2J2 assumes a catalytically competent binding mode stabilised by hydrogen bonds to Arg117.

On the molecular discrimination between adenine and guanine by proteins.

The molecular recognition and discrimination of adenine and guanine ligand moieties in complexes with proteins have been studied using empirical observations on carefully selected crystal structures. The distribution of protein folds that bind these purines has been found to differ significantly from that across the whole PDB, but the most populated architectures and folds are also the most common in three genomes from the three different domains of life. The protein environments around the two nucleic acid bases were significantly different, in terms of the propensities of amino acid residues to be in the binding site, as well as their propensities to form hydrogen bonds to the bases. Plots of the distribution of protein atoms around the two purines clearly show different clustering of hydrogen bond donors and acceptors opposite complimentary acceptors and donors in the rings, with hydrophobic areas below and above the rings. However, the clustering pattern is fuzzy, reflecting the variety of ways that proteins have evolved to recognise the same molecular moiety. Furthermore, an analysis of the conservation of residues in the protein chains binding guanine shows that residues in contact with the base are in general better conserved than the rest of the chain.