Building pangenome graphs.

Pangenome graphs can represent all variation between multiple reference genomes, but current approaches to build them exclude complex sequences or are based upon a single reference. In response, we developed the PanGenome Graph Builder, a pipeline for constructing pangenome graphs without bias or exclusion. The PanGenome Graph Builder uses all-to-all alignments to build a variation graph in which we can identify variation, measure conservation, detect recombination events and infer phylogenetic relationships.

BIMSA: accelerating long sequence alignment using processing-in-memory.

Motivation: Recent advances in sequencing technologies have stressed the critical role of sequence analysis algorithms and tools in genomics and healthcare research. In particular, sequence alignment is a fundamental building block in many sequence analysis pipelines and is frequently a performance bottleneck both in terms of execution time and memory usage. Classical sequence alignment algorithms are based on dynamic programming and often require quadratic time and memory with respect to the sequence length.

Effects of removal of the axial methionine heme ligand on the binding of S. cerevisiae iso-1 cytochrome c to cardiolipin.

The cleavage of the axial S(Met) - Fe bond in cytochrome c (cytc) upon binding to cardiolipin (CL), a glycerophospholipid of the inner mitochondrial membrane, is one of the key molecular changes that impart cytc with (lipo)peroxidase activity essential to its pro-apoptotic function. In this work, UV - VIS, CD, MCD and fluorescence spectroscopies were used to address the role of the Fe - M80 bond in controlling the cytc-CL interaction, by studying the binding of the Met80Ala (M80A) variant of S. cerevisiae iso-1 cytc (ycc) to CL liposomes in comparison with the wt protein [Paradisi et al.

WFA-GPU: gap-affine pairwise read-alignment using GPUs.

Motivation: Advances in genomics and sequencing technologies demand faster and more scalable analysis methods that can process longer sequences with higher accuracy. However, classical pairwise alignment methods, based on dynamic programming (DP), impose impractical computational requirements to align long and noisy sequences like those produced by PacBio and Nanopore technologies. The recently proposed wavefront alignment (WFA) algorithm paves the way for more efficient alignment tools, improving time and memory complexity over previous methods.

A draft human pangenome reference.

Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals. These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels.

Building pangenome graphs.

Pangenome graphs can represent all variation between multiple reference genomes, but current approaches to build them exclude complex sequences or are based upon a single reference. In response, we developed the PanGenome Graph Builder (PGGB), a pipeline for constructing pangenome graphs without bias or exclusion. PGGB uses all-to-all alignments to build a variation graph in which we can identify variation, measure conservation, detect recombination events, and infer phylogenetic relationships.

Optimal gap-affine alignment in O(s) space.

Motivation: Pairwise sequence alignment remains a fundamental problem in computational biology and bioinformatics. Recent advances in genomics and sequencing technologies demand faster and scalable algorithms that can cope with the ever-increasing sequence lengths. Classical pairwise alignment algorithms based on dynamic programming are strongly limited by quadratic requirements in time and memory.

Thermodynamics and Kinetics of Electron Transfer of Electrode-Immobilized Small Laccase from .

The thermodynamic and kinetic properties for heterogeneous electron transfer (ET) were measured for the electrode-immobilized small laccase (SLAC) from subjected to different electrostatic and covalent protein-electrode linkages, using cyclic voltammetry. Once immobilized electrostatically onto a gold electrode using mixed carboxyl- and hydroxy-terminated alkane-thiolate SAMs or covalently exploiting the same SAM subjected to N-hydroxysuccinimide+1-Ethyl-3-(3-dimethylaminopropyl)carbodiimide (NHS-EDC) chemistry, the SLAC-electrode electron flow occurs through the T1 center. The E°' values (from +0.

Assessing the Functional and Structural Stability of the Met80Ala Mutant of Cytochrome in Dimethylsulfoxide.

The Met80Ala variant of yeast cytochrome is known to possess electrocatalytic properties that are absent in the wild type form and that make it a promising candidate for biocatalysis and biosensing. The versatility of an enzyme is enhanced by the stability in mixed aqueous/organic solvents that would allow poorly water-soluble substrates to be targeted. In this work, we have evaluated the effect of dimethylsulfoxide (DMSO) on the functionality of the Met80Ala cytochrome mutant, by investigating the thermodynamics and kinetics of electron transfer in mixed water/DMSO solutions up to 50% DMSO /.

Hydrogen peroxide induces heme degradation and protein aggregation in human neuroglobin: roles of the disulfide bridge and hydrogen-bonding in the distal heme cavity.

In the present study, human neuroglobin (hNgb) was found to undergo H O -induced breakdown of the heme center at a much slower rate than other globins, namely in the timescale of hours against minutes. We investigated how the rate of the process is affected by the Cys46/Cys55 disulfide bond and the network of non-covalent interactions in the distal heme side involving Tyr44, Lys67, the His64 heme iron axial ligand and the heme propionate-7. The rate is increased by the Tyr44 to Ala and Phe mutations; however the rate is lowered by Lys67 to Ala swapping.

AnchorWave: Sensitive alignment of genomes with high sequence diversity, extensive structural polymorphism, and whole-genome duplication.

Millions of species are currently being sequenced, and their genomes are being compared. Many of them have more complex genomes than model systems and raise novel challenges for genome alignment. Widely used local alignment strategies often produce limited or incongruous results when applied to genomes with dispersed repeats, long indels, and highly diverse sequences.

Activity and substrate specificity of lytic polysaccharide monooxygenases: An ATR FTIR-based sensitive assay tested on a novel species from Pseudomonas putida.

Pseudomonas putida W619 is a soil Gram-negative bacterium commonly used in environmental studies thanks to its ability in degrading many aromatic compounds. Its genome contains several putative carbohydrate-active enzymes such as glycoside hydrolases and lytic polysaccharide monooxygenases (PMOs). In this study, we have heterologously produced in Escherichia coli and characterized a new enzyme belonging to the AA10 family, named PpAA10 (Uniprot: B1J2U9), which contains a chitin-binding type-4 module and showed activity toward β-chitin.

Phosphodiesterase (PDE) 5 inhibitors sildenafil, tadalafil and vardenafil impact cAMP-specific PDE8 isoforms-linked second messengers and steroid production in a mouse Leydig tumor cell line.

Type 5 phosphodiesterase (PDE5) blockade by inhibitors (PDE5i) results in intracellular cyclic guanosine monophosphate (cGMP) increase and smooth muscle relaxation and are used for the treatment of men erectile dysfunction. Although they have high specificity for PDE5, these inhibitors are suspected to cross-interact also with cyclic adenosine monophosphate (cAMP)-specific PDEs, inducing the intracellular accumulation of this cyclic nucleotide and related testosterone increase, positively impacting male reproductive parameters. However, the link between the use of PDE5i and the activation of cAMP-mediated steroidogenesis is still unclear.

Pseudoperoxidase activity, conformational stability, and aggregation propensity of the His98Tyr myoglobin variant: implications for the onset of myoglobinopathy.

The autosomal dominant striated muscle disease myoglobinopathy is due to the single point mutation His98Tyr in human myoglobin (MB), the heme protein responsible for binding, storage, and controlled release of O in striated muscle. In order to understand the molecular basis of this disease, a comprehensive biochemical and biophysical study on wt MB and the variant H98Y has been performed. Although only small differences exist between the active site architectures of the two proteins, the mutant (a) exhibits an increased reactivity toward hydrogen peroxide, (b) exhibits a higher tendency to form high-molecular-weight aggregates, and (c) is more prone to heme bleaching, possibly as a consequence of the observed H O -induced formation of the Tyr98 radical close to the metal center.

How to Turn an Electron Transfer Protein into a Redox Enzyme for Biosensing.

Cytochrome is a small globular protein whose main physiological role is to shuttle electrons within the mitochondrial electron transport chain. This protein has been widely investigated, especially as a paradigmatic system for understanding the fundamental aspects of biological electron transfer and protein folding. Nevertheless, cytochrome can also be endowed with a non-native catalytic activity and be immobilized on an electrode surface for the development of third generation biosensors.

Fast gap-affine pairwise alignment using the wavefront algorithm.

Motivation: Pairwise alignment of sequences is a fundamental method in modern molecular biology, implemented within multiple bioinformatics tools and libraries. Current advances in sequencing technologies press for the development of faster pairwise alignment algorithms that can scale with increasing read lengths and production yields.

Results: In this article, we present the wavefront alignment algorithm (WFA), an exact gap-affine algorithm that takes advantage of homologous regions between the sequences to accelerate the alignment process.

Adsorbing surface strongly influences the pseudoperoxidase and nitrite reductase activity of electrode-bound yeast cytochrome c. The effect of hydrophobic immobilization.

The Met80Ala and Met80Ala/Tyr67Ala variants of S. cerevisiae iso-1 cytochrome c (ycc) and their adducts with cardiolipin immobilized onto a gold electrode coated with a hydrophobic self-assembled monolayer (SAM) of decane-1-thiol were studied through cyclic voltammetry and surface-enhanced resonance Raman spectroscopy (SERRS). The electroactive species - containing a six-coordinate His/His axially ligated heme and a five-coordinate His/- heme stable in the oxidized and reduced state, respectively - and the pseudoperoxidase activity match those found previously for the wt species and are only slightly affected by CL binding.

Binding of S. cerevisiae iso-1 cytochrome c and its surface lysine-to-alanine variants to cardiolipin: charge effects and the role of the lipid to protein ratio.

The interaction of cytochrome c with cardiolipin (CL) is a critical step in the initial stages of apoptosis and is mediated by a positively charged region on the protein surface comprising several lysine residues (site A). Here, the interaction of wt S. cerevisiae cytochrome c (ycc) and its K72A/K73A, K72A/K79A, K73A/K79A and K72A/K73A/K79A variants with CL was studied through UV-Vis and MCD spectroscopies at pH 7 and molecular dynamics (MD) simulations, to clarify the role of the mutated lysines.

Enamel peptides reveal the sex of the Late Antique 'Lovers of Modena'.

Recent work has disclosed the critical role played by enamel peptides in sex classification of old skeletal remains. In particular, protein AMELY (amelogenin isoform Y) is present in the enamel dental tissue of male individuals only, while AMELX (isoform X) can be found in both sexes. AMELY can be easily detected by LC-MS/MS in the ion extracted chromatograms of the SMIRPPY peptide (monoisotopic [M + 2 H] mass = 440.

Trapping at the Solid-Gas Interface: Selective Adsorption of Naphthalene by Montmorillonite Intercalated with a Fe(III)-Phenanthroline Complex.

In this study, stable hybrid materials (Mt-Fe(III)Phen), made by the μ-oxo Fe(III)-phenanthroline complex [(OH)(Phen)FeOFe(Phen)(OH)] (Fe(III)Phen) intercalated in different amounts into montmorillonite (Mt), were used as a trap for immobilizing gaseous benzene and naphthalene and their mono chloro-derivatives at 25 and 50 °C. The entrapping process was studied through elemental analysis, magic angle spinning NMR spectroscopy, thermal analysis, and evolved gas mass spectrometry. Naphthalene and 1-chloronaphthalene were found to be immobilized in large amount at both temperatures.

Probing the Effect of Sildenafil on Progesterone and Testosterone Production by an Intracellular FRET/BRET Combined Approach.

Forster resonance energy transfer (FRET)-based biosensors have been recently applied to the study of biological pathways. In this study, a new biosensor was validated for the first time in live HEK293 and steroidogenic MLTC-1 cell lines for studying the effect of the PDE5 inhibitor on the hCG/LH-induced steroidogenic pathway. The sensor improves FRET between a donor (D), the fluorescein-like diarsenical probe that can covalently bind a tetracysteine motif fused to the PDE5 catalytic domain, and an acceptor (A), the rhodamine probe conjugated to the pseudosubstrate cGMPS.

gemBS: high throughput processing for DNA methylation data from bisulfite sequencing.

Motivation: DNA methylation is essential for normal embryogenesis and development in mammals and can be captured at single base pair resolution by whole genome bisulfite sequencing (WGBS). Current available analysis tools are becoming rapidly outdated as they lack sensible functionality and efficiency to handle large amounts of data now commonly created.

Results: We developed gemBS, a fast high-throughput bioinformatics pipeline specifically designed for large scale BS-Seq analysis that combines a high performance BS-mapper (GEM3) and a variant caller specifically for BS-Seq data (BScall).

Catalytic Mechanism of Fungal Lytic Polysaccharide Monooxygenases Investigated by First-Principles Calculations.

Lytic polysaccharide monooxygenases (LPMOs) are Cu-containing enzymes that facilitate the degradation of recalcitrant polysaccharides by the oxidative cleavage of glycosidic bonds. They are gaining rapidly increasing attention as key players in biomass conversion, especially for the production of second-generation biofuels. Elucidation of the detailed mechanism of the LPMO reaction is a major step toward the assessment and optimization of LPMO efficacy in industrial biotechnology, paving the way to utilization of sustainable fuel sources.

Fluorometric detection of protein-ligand engagement: The case of phosphodiesterase5.

Phosphodiesterases (PDEs) regulate the intracellular levels of cAMP and cGMP. The great clinical success of the PDE5 inhibitors, Sildenafil (Viagra), Vardenafil (Levitra) and Tadalafil (Cialis) has led to an increasing interest for this class of enzymes. Recent studies have shown a correlation between tumor growth and PDE5 overexpression, making PDE5-selective inhibitors promising candidates for cancer treatment.

The influence of the Cys46/Cys55 disulfide bond on the redox and spectroscopic properties of human neuroglobin.

Neuroglobin is a monomeric globin containing a six-coordinate heme b, expressed in the nervous system, which exerts an important neuroprotective role. In the human protein (hNgb), Cys46 and Cys55 form an intramolecular disulfide bond under oxidizing conditions, whose cleavage induces a helix-to-strand rearrangement of the CD loop that strengthens the bond between the heme iron and the distal histidine. Hence, it is conceivable that the intramolecular disulfide bridge modulates the functionality of human neuroglobin by controlling exogenous ligand binding.