A fast and flexible approach to oligonucleotide probe design for genomes and gene families

MOTIVATION: With hundreds of completely sequenced microbial genomes available, and advancements in DNA microarray technology, the detection of genes in microbial communities consisting of hundreds of thousands of sequences may be possible. The existing strategies developed for DNA probe design, geared toward identifying specific sequences, are not suitable due to the lack of coverage, flexibility and efficiency necessary for applications in metagenomics. METHODS: ProDesign is a tool developed for the selection of oligonucleotide probes to detect members of gene families present in environmental samples. Gene family-specific probe sequences are generated based on specific and shared words, which are found with the spaced seed hashing algorithm. To detect more sequences, those sharing some common words are re-clustered into new families, then probes specific for the new families are generated. RESULTS: The program is very flexible in that it can be used for designing probes for detecting many genes families simultaneously and specifically in one or more genomes. Neither the length nor the melting temperature of the probes needs to be predefined. We have found that ProDesign provides more flexibility, coverage and speed than other software programs used in the selection of probes for genomic and gene family arrays. AVAILABILITY: ProDesign is licensed free of charge to academic users. ProDesign and Supplementary Material can be obtained by contacting the authors. A web server for ProDesign is available at http://www.uhnresearch.ca/labs/tillier/ProDesign/ProDesign.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

FiberSim: A flexible open-source model of myofilament-level contraction

FiberSim is a flexible open-source model of myofilament-level contraction. The code uses a spatially explicit technique, meaning that it tracks the position and status of each contractile molecule within the lattice framework. This allows the model to simulate some of the mechanical effects modulated by myosin-binding protein C, as well as the dose dependence of myotropes and the effects of varying isoform expression levels. This paper provides a short introduction to FiberSim and presents simulations of tension-pCa curves with and without regulation of thick and thin filament activation by myosin-binding protein C. A myotrope dose-dependent response as well as slack/re-stretch maneuvers to assess rates of tension recovery are also presented. The software was designed to be flexible (the user can define their own model and/or protocol) and computationally efficient (simulations can be performed on a regular laptop). We hope that other investigators will use FiberSim to explore myofilament level mechanisms and to accelerate research focusing on the contractile properties of sarcomeres.     

A short, novel, and cheaper procedure for oligonucleotide synthesis using automated solid phase synthesizer

Dimethylthiuram disulfide (DTD) has been developed as an efficient thiolation reagent during automated synthesis of oligonucleotides using phosphoramidite chemistry. Simultaneous thiolation and capping was accomplished by mixing DTD with capping solution B, which saved 20% of solvent consumption and compressed the four-step synthesis cycle to three. Large-scale (1 mmol) synthesis of phosphorothioate oligonucleotides has been demonstrated with excellent yield and purity.     

The GlycanBuilder: a fast, intuitive and flexible software tool for building and displaying glycan structures

Background  

Carbohydrates play a critical role in human diseases and their potential utility as biomarkers for pathological conditions is a major driver for characterization of the glycome. However, the additional complexity of glycans compared to proteins and nucleic acids has slowed the advancement of glycomics in comparison to genomics and proteomics. The branched nature of carbohydrates, the great diversity of their constituents and the numerous alternative symbolic notations, make the input and display of glycans not as straightforward as for example the amino-acid sequence of a protein. Every glycoinformatic tool providing a user interface would benefit from a fast, intuitive, appealing mechanism for input and output of glycan structures in a computer readable format.     

R-Gada: a fast and flexible pipeline for copy number analysis in association studies

Background  

Genome-wide association studies (GWAS) using Copy Number Variation (CNV) are becoming a central focus of genetic research. CNVs have successfully provided target genome regions for some disease conditions where simple genetic variation (i.e., SNPs) has previously failed to provide a clear association.     

In situ oligonucleotide synthesis on poly(dimethylsiloxane): a flexible substrate for microarray fabrication

In this paper, we demonstrate in situ synthesis of oligonucleotide probes on poly(dimethylsiloxane) (PDMS) microchannels through use of conventional phosphoramidite chemistry. PDMS polymer was moulded into a series of microchannels using standard soft lithography (micro-moulding), with dimensions <100 μm. The surface of the PDMS was derivatized by exposure to ultraviolet/ozone followed by vapour phase deposition of glycidoxypropyltrimethoxysilane and reaction with poly(ethylene glycol) spacer, resulting in a reactive surface for oligonucleotide coupling. High, reproducible yields were achieved for both 6mer and 21mer probes as assessed by hybridization to fluorescent oligonucleotides. Oligonucleotide surface density was comparable with that obtained on glass substrates. These results suggest PDMS as a stable and flexible alternative to glass as a suitable substrate in the fabrication and synthesis of DNA microarrays.     

A guide to Open-JIP,a low-cost open-source chlorophyll fluorometer

Photosynthesis Research - Chlorophyll a fluorescence is the most widely used method to study photosynthesis and plant stress. While several commercial fluorometers are available, there is a need...     

OpenAWSEM with Open3SPN2: A fast,flexible, and accessible framework for large-scale coarse-grained biomolecular simulations

We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.     

Construction of linear functional expression elements with DNA and RNA hybrid primers: a flexible and fast method for proteomics

Two methods of constructing linear functional expression elements (LFEE) using hybrid DNA and RNA primers in DNA amplification for rapid gene expression are described. In both methods, it is not necessary to have additional transformation or bacterial propagation. The promoter, open reading frame (ORF) and terminator are amplified using Pfu or Taq DNA polymerase. Three elements containing DNA or RNA overhang are covalently ligated by T4 DNA ligase. The recombinant molecule is amplified with element-specific primers. The LFEE can be generated by both methods in a few hours and can be expressed in mammalian cells.     

Minimus: a fast,lightweight genome assembler

Background  

Genome assemblers have grown very large and complex in response to the need for algorithms to handle the challenges of large whole-genome sequencing projects. Many of the most common uses of assemblers, however, are best served by a simpler type of assembler that requires fewer software components, uses less memory, and is far easier to install and run.     

Technical note: GAMMORA,a free,open-source,and validated GATE-based model for Monte-Carlo simulations of the Varian TrueBeam

PurposeMonte Carlo (MC) is the reference computation method for medical physics. In radiotherapy, MC computations are necessary for some issues (such as assessing figures of merit, double checks, and dose conversions). A tool based on GATE is proposed to easily create full MC simulations of the Varian TrueBeam STx.MethodsGAMMORA is a package that contains photon phase spaces as a pre-trained generative adversarial network (GAN) and the TrueBeam’s full geometry. It allows users to easily create MC simulations for simple or complex radiotherapy plans such as VMAT. To validate the model, the characteristics of generated photons are first compared to those provided by Varian (IAEA format). Simulated data are also compared to measurements in water and heterogeneous media. Simulations of 8 SBRT plans are compared to measurements (in a phantom). Two examples of applications (a second check and interplay effect assessment) are presented.ResultsThe simulated photons generated by the GAN have the same characteristics (energy, position, and direction) as the IAEA data. Computed dose distributions of simple cases (in water) and complex plans delivered in a phantom are compared to measurements, and the Gamma index (3%/3mm) was always superior to 98%. The feasibility of both clinical applications is shown.ConclusionsThis model is now shared as a free and open-source tool that generates radiotherapy MC simulations. It has been validated and used for five years. Several applications can be envisaged for research and clinical purposes.     

The SpikerBox: a low cost, open-source bioamplifier for increasing public participation in neuroscience inquiry

Although people are generally interested in how the brain functions, neuroscience education for the public is hampered by a lack of low cost and engaging teaching materials. To address this, we developed an open-source tool, the SpikerBox, which is appropriate for use in middle/high school educational programs and by amateurs. This device can be used in easy experiments in which students insert sewing pins into the leg of a cockroach, or other invertebrate, to amplify and listen to the electrical activity of neurons. With the cockroach leg preparation, students can hear and see (using a smartphone oscilloscope app we have developed) the dramatic changes in activity caused by touching the mechanosensitive barbs. Students can also experiment with other manipulations such as temperature, drugs, and microstimulation that affect the neural activity. We include teaching guides and other resources in the supplemental materials. These hands-on lessons with the SpikerBox have proven to be effective in teaching basic neuroscience.     

Affordable biocomputing for everyone: using the Internet,freeware and open-source software

How to build your own complete working biocomputing platform with nothing more than a desktop computer and an Internet connection.     

Using an aryl phenanthroimidazole moiety as a conjugated flexible intercalator to improve the hybridization efficiency of a triplex-forming oligonucleotide

When inserting 2-phenyl or 2-naphth-1-yl-phenanthroimidazole intercalators (X and Y, respectively) as bulges into triplex-forming oligonucleotides, both intercalators show extraordinary high thermal stability of the corresponding Hoogsteen-type triplexes and Hoogsteen-type parallel duplexes with high discrimination to Hoogsteen mismatches. Molecular modeling shows that the phenyl or the naphthyl ring stacks with the nucleobases in the TFO, while the phenanthroimidazol moiety stacks with the base pairs of the dsDNA. DNA-strands containing the intercalator X show higher thermal triplex stability than DNA-strands containing the intercalator Y. The difference can be explained by a lower degree of planarity of the intercalator in the case of naphthyl. It was also observed that triplex stability was considerably reduced when the intercalators X or Y was replaced by 2-(naphthlen-1-yl)imidazole. This confirms intercalation as the important factor for triplex stabilization and it rules out an alternative complexation of protonated imidazole with two phosphate groups. The intercalating nucleic acid monomers X and Y were obtained via a condensation reaction of 9,10-phenanthrenequinone (4) with (S)-4-(2-(2,2-dimethyl-1,3-dioxolan-4-yl)ethoxy)benzaldehyde (3a) or (S)-4-(2-(2,2-dimethyl-1,3-dioxolan-4-yl)ethoxy)-1-naphthaldehyde (3b), respectively, in the presence of acetic acid and ammonium acetate. The required monomers for DNA synthesis using amidite chemistry were obtained by standard deprotection of the hydroxy groups followed by 4,4'-dimethoxytritylation and phosphitylation.     

The spliceosome: a flexible, reversible macromolecular machine

With more than a hundred individual RNA and protein parts and a highly dynamic assembly and disassembly pathway, the spliceosome is arguably the most complicated macromolecular machine in the eukaryotic cell. This complexity has made kinetic and mechanistic analysis of splicing incredibly challenging. Yet, recent technological advances are now providing tools for understanding this process in much greater detail. Ranging from genome-wide analyses of splicing and creation of an orthogonal spliceosome in vivo, to purification of active spliceosomes and observation of single molecules in vitro, such new experimental approaches are yielding significant insight into the inner workings of this remarkable machine. These experiments are rewriting the textbooks, with a new picture emerging of a dynamic, malleable machine heavily influenced by the identity of its pre-mRNA substrate.