Bacillus thuringiensis (Bt) is a bacterium capable of producing Cry toxins, which are recognized for their bio-controlling actions against insects. However, a few Bt strains encode proteins lacking insecticidal activity but showing cytotoxic activity against different cancer cell lines and low or no cytotoxicity toward normal human cells. A subset of Cry anticancer proteins, termed parasporins (PSs), has recently arisen as a potential alternative for cancer treatment. However, the molecular receptors that allow the binding of PSs to cells and their cytotoxic mechanisms of action have not been well established. Nonetheless, their selective cytotoxic activity against different types of cancer cell lines places PSs as a promising alternative treatment modality. In this review, we provide an overview of the classification, structures, mechanisms of action, and insights obtained from genetic modification approaches for PS proteins.

Objective. To identify DNA methylation patterns of heavy smokers in oral rinse samples. Methods. Genome-wide DNA methylation data was imported from Gene Expression Omnibus GSE70977 using the GEOquery package. Two independent sets were analyzed: (a) 71 epigenomes of cancer-free subjects (heavy smokers n = 37 vs. non-smokers n = 31); for concordance assessment (b) 139 oral-cancer patients’ epigenomes (heavy smokers n = 92 vs. non-smokers n = 47). Differential DNA methylation for CpG positions and at the regional level was determined using Limma and DMRcate Bioconductor packages. The linear model included sex, age, and alcohol consumption. The statistical threshold was set to p < 0.05. Functional gene prioritization analysis was performed for gene-targeted analysis. Results. In individuals without cancer and heavy smokers, the FAM184B gene was found with two CpG positions differentially hypermethylated (p = 0.012 after FDR adjustment), in a region of 48 bp with an absolute methylation difference >10% between groups (p = 1.76 × 10⁻⁸). In the analysis corresponding to oral-cancer patients, we found AHRR differentially hypomethylated cancer patients, but also in subjects without oral cancer in the targeted analyses. Remarkably, ADAMTS2 was found differentially hypermethylated in heavy smokers without a diagnosis of cancer in two consecutive probes cg05575921 (p = 3.13 × 10⁻⁷) and cg10208897 (p = 1.36 × 10⁻⁵). Conclusions: Differentially methylated AHRR, ADAMTS2, and FAM184B genes are biomarker candidates in oral rinse samples.

The Cry toxins, or δ-endotoxins, are a diverse group of proteins produced by Bacillus thuringiensis. While DNA secondary structures are biologically relevant, it is unknown if such structures are formed in regions encoding conserved domains of Cry toxins under shuffling conditions. We analyzed 5 holotypes that encode Cry toxins and that grouped into 4 clusters according to their phylogenetic closeness. The mean number of DNA secondary structures that formed and the mean Gibbs free energy $$ \left(\overline{\varDelta G}\right) $$ were determined by an in silico analysis using different experimental DNA shuffling scenarios. In terms of spontaneity, shuffling efficiency was directly proportional to the formation of secondary structures but inversely proportional to ∆G. The results showed a shared thermodynamic pattern for each cluster and relationships among sequences that are phylogenetically close at the protein level. The regions of the cry11Aa, Ba and Bb genes that encode domain I showed more spontaneity and thus a greater tendency to form secondary structures (<∆G). In the region of domain III; this tendency was lower (>∆G) in the cry11Ba and Bb genes. Proteins that are phylogenetically closer to Cry11Ba and Cry11Bb, such as Cry2Aa and Cry18Aa, maintained the same thermodynamic pattern. More distant proteins, such as Cry1Aa, Cry1Ab, Cry30Aa and Cry30Ca, featured different thermodynamic patterns in their DNA. These results suggest the presence of thermodynamic variations associated to the formation of secondary structures and an evolutionary relationship with regions that encode highly conserved domains in Cry proteins. The findings of this study may have a role in the in silico design of cry gene assembly by DNA shuffling techniques.

Parasporin 2 has cytotoxic effects against numerous colon cancer cell lines, making it a viable alternative to traditional treatments. However, its mechanism of action and receptors remain unknown. In this study, site-directed mutagenesis was used to obtain PS2Aa1 mutants with variation in domain I at positions 256 and 257. Variants 015, 002, 3-3, 3-35, and 3-45 presented G256A, G256E, G257A, G257V, and G257E substitutions, respectively. Cytotoxicity tests were performed for the cell viability of cell lines SW480, SW620, and CaCo-2. Mutants 3-3, 3-35, and 3-45 efficiently killed the cell lines. It was found that the activated forms of caspase-3 and PARP were in higher abundance as well as increased production of γH2AX when 3-35 was used to treat CaCo-2 and SW480. To assess possible membrane-binding receptors involved in the interaction, an APN receptor blocking assay showed reduced activity of some parasporins. Hence, we performed molecular docking and molecular dynamics simulations to analyze the stability of possible interactions and identify the residues that could be involved in the protein–protein interaction of PS2Aa1 and APN. We found that residues 256 and 257 facilitate the interaction. Parasporin 3-35 is promising because it has higher cytotoxicity than PS2Aa1.

Directed evolution methods mimic in vitro Darwinian evolution, inducing random mutations and selective pressure in genes to obtain proteins with enhanced characteristics. These techniques are developed using trial-and-error testing at an experimental level with a high degree of uncertainty. Therefore, in silico modeling of directed evolution is required to support experimental assays. Several in silico approaches have reproduced directed evolution, using statistical, thermodynamic, and kinetic models in an attempt to recreate experimental conditions. Likewise, optimization techniques using heuristic models have been used to understand and find the best scenarios of directed evolution. Our study uses an in silico model named HeurIstics DirecteD EvolutioN, which is based on a genetic algorithm designed to generate chimeric libraries from 2 parental genes, cry11Aa and cry11Ba, of Bacillus thuringiensis. These genes encode crystal-shaped δ-endotoxins with 3 conserved domains. Cry11 toxins are of biotechnological interest because they have shown to be effective as biopesticides for disease-spreading vectors. With our heuristic model, we considered experimental parameters such as DNA fragmentation length, number of generations or simulation cycles, and mutation rate, to get characteristics of Cry11 chimeric libraries such as percentage of population identity, truncation of variants obtained from the presence of internal stop codons, percentage of thermodynamic diversity, and stability of variants. Our study allowed us to focus on experimental conditions that may be useful for the design of in vitro and in silico experiments of directed evolution with Cry toxins of 3 conserved domains. Furthermore, we obtained in silico libraries of Cry11 variants, in which structural characteristics of wild Cry families were observed in a review of a sample of in silico sequences. We consider that future studies could use our in silico libraries and heuristic computational models, as the one suggested here, to support in vitro experiments of directed evolution.

Polymerase Chain Reaction (PCR) based techniques for DNA methylation techniques includes MS-HRM technique. Methylation Sensitive High-Resolution Melting (MS-HRM) primer-design requires a set of necessary recommendations for such DNA methylation assessment. However, there were not any available software that allows an automatic design of this kind primers. We present Softepigen, the first complete MS-HRM primer design software. Softepigen allows to search for primers in a genomic region following Wojdacz's recommendations and targets primer binding regions with high linguistic complexity sequences that increase the specificity of the converted sequence of the human genome. We performed in-silico PCR analysis through BiSearch ePCR tool to validate the specificity of the of the primers designed using Softepigen. Softepigen for MS-HRM performance in our genomic regions of interest show satisfactory specificity measurements, and we implemented it for freely available use in web-based interface in www.soft-epigen.com.