Thermostabilizing mechanisms of canonical single amino acid substitutions at a GH1 ß-glucosidase probed by multiple MD and computational approaches.

ß-glucosidases play a pivotal role in second-generation biofuel (2G-biofuel) production. For this application, thermostable enzymes are essential due to the denaturing conditions on the bioreactors. Random amino acid substitutions have originated new thermostable ß-glucosidases, but without a clear understanding of their molecular mechanisms. Here, we probe by different molecular dynamics simulation approaches with distinct force fields and submitting the results to various computational analyses, the molecular bases of the thermostabilization of the Paenibacillus polymyxa GH1 ß-glucosidase by two-point mutations E96K (TR1) and M416I (TR2). Equilibrium molecular dynamic simulations (eMD) at different temperatures, principal component analysis (PCA), virtual docking, metadynamics (MetaDy), accelerated molecular dynamics (aMD), Poisson-Boltzmann surface analysis, grid inhomogeneous solvation theory and colony method estimation of conformational entropy allow to converge to the idea that the stabilization carried by both substitutions depend on different contributions of three classic mechanisms: (i) electrostatic surface stabilization; (ii) efficient isolation of the hydrophobic core from the solvent, with energetic advantages at the solvation cap; (iii) higher distribution of the protein dynamics at the mobile active site loops than at the protein core, with functional and entropic advantages. Mechanisms i and ii predominate for TR1, while in TR2, mechanism iii is dominant. Loop A integrity and loops A, C, D, and E dynamics play critical roles in such mechanisms. Comparison of the dynamic and topological changes observed between the thermostable mutants and the wildtype protein with amino acid co-evolutive networks and thermostabilizing hotspots from the literature allow inferring that the mechanisms here recovered can be related to the thermostability obtained by different substitutions along the whole family GH1. We hope the results and insights discussed here can be helpful for future rational approaches to the engineering of optimized ß-glucosidases for 2G-biofuel production for industry, biotechnology, and science.

A Computational Method to Propose Mutations in Enzymes Based on Structural Signature Variation (SSV).

With the use of genetic engineering, modified and sometimes more efficient enzymes can be created for different purposes, including industrial applications. However, building modified enzymes depends on several in vitro experiments, which may result in the process being expensive and time-consuming. Therefore, computational approaches could reduce costs and accelerate the discovery of new technological products. In this study, we present a method, called structural signature variation (SSV), to propose mutations for improving enzymes' activity. SSV uses the structural signature variation between target enzymes and template enzymes (obtained from the literature) to determine if randomly suggested mutations may provide some benefit for an enzyme, such as improvement of catalytic activity, half-life, and thermostability, or resistance to inhibition. To evaluate SSV, we carried out a case study that suggested mutations in ß-glucosidases: Essential enzymes used in biofuel production that suffer inhibition by their product. We collected 27 mutations described in the literature, and manually classified them as beneficial or not. SSV was able to classify the mutations with values of 0.89 and 0.92 for precision and specificity, respectively. Then, we used SSV to propose mutations for Bgl1B, a low-performance ß-glucosidase. We detected 15 mutations that could be beneficial. Three of these mutations (H228C, H228T, and H228V) have been related in the literature to the mechanism of glucose tolerance and stimulation in GH1 ß-glucosidase. Hence, SSV was capable of detecting promising mutations, already validated by in vitro experiments, that improved the inhibition resistance of a ß-glucosidase and, consequently, its catalytic activity. SSV might be useful for the engineering of enzymes used in biofuel production or other industrial applications.

Molecular Dynamics Gives New Insights into the Glucose Tolerance and Inhibition Mechanisms on ß-Glucosidases.

ß-Glucosidases are enzymes with high importance for many industrial processes, catalyzing the last and limiting step of the conversion of lignocellulosic material into fermentable sugars for biofuel production. However, ß-glucosidases are inhibited by high concentrations of the product (glucose), which limits the biofuel production on an industrial scale. For this reason, the structural mechanisms of tolerance to product inhibition have been the target of several studies. In this study, we performed in silico experiments, such as molecular dynamics (MD) simulations, free energy landscape (FEL) estimate, Poisson-Boltzmann surface area (PBSA), and grid inhomogeneous solvation theory (GIST) seeking a better understanding of the glucose tolerance and inhibition mechanisms of a representative GH1 ß-glucosidase and a GH3 one. Our results suggest that the hydrophobic residues Y180, W350, and F349, as well the polar one D238 act in a mechanism for glucose releasing, herein called "slingshot mechanism", dependent also on an allosteric channel (AC). In addition, water activity modulation and the protein loop motions suggest that GH1 ß-Glucosidases present an active site more adapted to glucose withdrawal than GH3, in consonance with the GH1s lower product inhibition. The results presented here provide directions on the understanding of the molecular mechanisms governing inhibition and tolerance to the product in ß-glucosidases and can be useful for the rational design of optimized enzymes for industrial interests.

Whole-genome optical mapping reveals a mis-assembly between two rRNA operons of Corynebacterium pseudotuberculosis strain 1002.

BACKGROUND: Studies have detected mis-assemblies in genomes of the species Corynebacterium pseudotuberculosis. These new discover have been possible due to the evolution of the Next-Generation Sequencing platforms, which have provided sequencing with accuracy and reduced costs. In addition, the improving of techniques for construction of high accuracy genomic maps, for example, Whole-genome mapping (WGM) (OpGen Inc), have allow high-resolution assembly that can detect large rearrangements. RESULTS: In this work, we present the resequencing of Corynebacterium pseudotuberculosis strain 1002 (Cp1002). Cp1002 was the first strain of this species sequenced in Brazil, and its genome has been used as model for several studies in silico of caseous lymphadenitis disease. The sequencing was performed using the platform Ion PGM and fragment library (200 bp kit). A restriction map was constructed, using the technique of WGM with the enzyme KpnI. After the new assembly process, using WGM as scaffolder, we detected a large inversion with size bigger than one-half of genome. A specific analysis using BLAST and NR database shows that the inversion occurs between two homology RNA ribosomal regions. CONCLUSION: In conclusion, the results showed by WGM could be used to detect mismatches in assemblies, providing genomic maps with high resolution and allow assemblies with more accuracy and completeness. The new assembly of C. pseudotuberculosis was deposited in GenBank under the accession no. CP012837.

Interaction of DisBa01 peptide from Bothrops alternatus venom with BRAF melanoma receptors: Modeling and molecular docking.

Metastatic melanoma is highly aggressive and challenging, often leading to a grim prognosis. Its progression is swift, especially when mutations like BRAFV600E continuously activate pathways vital for cell growth and survival. Although several treatments target this mutation, resistance typically emerges over time. In recent decades, research has underscored the potential of snake venoms and peptides as bioactive substances for innovative drugs, including anti-coagulants, anti-microbial, and anti-cancer agents. Leveraging this knowledge, we propose employing a bioinformatics simulation approach to: a) Predict how well a peptide (DisBa01) from Bothrops alternatus snake venom binds to the melanoma receptor BRAFV600E via Molecular Docking. b) Identify the specific peptide binding sites on receptors and analyze their proximity to active receptor sites. c) Evaluate the behavior of resulting complexes through molecular dynamics simulations. d) Assess whether this peptide qualifies as a candidate for anti-melanoma therapy. Our findings reveal that DisBa01 enhances stability in the BRAFV600E melanoma receptor structure by binding to its RGD motif, an interaction absent in the BRAF WT model. Consequently, both docking and molecular dynamics simulations suggest that DisBa01 shows promise as a BRAFV600E inhibitor.

Conformational flexibility correlates with glucose tolerance for point mutations in ß-glucosidases - a computational study.

ß-glucosidases (EC 3.2.1.21) have been described as essential to second-generation biofuel production. They act in the last step of the lignocellulosic saccharification, cleaving the ß - 1,4 glycosidic bonds in cellobiose to produce two molecules of glucose. However, ß-glucosidases have been described as strongly inhibited by glucose, causing an increment of cellobiose concentration. Also, cellobiose is an inhibitor of other enzymes used in this process, such as exoglucanases and endoglucanases. Hence, the engineering of thermostable and glucose-tolerant ß-glucosidases has been targeted by many studies. In this study, we performed high sampling accelerated molecular dynamics for a wild glucose-tolerant GH1 ß-glucosidase (Bgl1A), a wild non-tolerant (Bgl1B), and a set of glucose-tolerant Bgl1B's mutants: V302F, N301Q/V302F, F172I, V227M, G246S, T299S, and H228T. Our results suggest that point mutations promissory to induce glucose tolerance trend to enhance the mobility of the flexible loops around the active site. Mutations affected B and C loops regions, and an αß-hairpin motif between them. Conformational clusters and free energy landscape profiles suggest that the mobility acquired by mutants allows a higher closure of the substrate channel. This closure is compatible with a higher impedance for glucose entrance and stimulus of its withdrawal. Based on mutants' structural analyses, we inferred that both the direct stereochemical effect on the glucose path and the changes in the mobility affect glucose tolerance. We hope these results be useful for the rational design of glucose-tolerant and industrially promising enzymes.Communicated by Ramaswamy H. Sarma.

Complete genome sequence of Streptococcus agalactiae strain GBS85147 serotype of type Ia isolated from human oropharynx.

Streptococcus agalactiae, also referred to as Group B Streptococcus, is a frequent resident of the rectovaginal tract in humans, and a major cause of neonatal infection. The pathogen can also infect adults with underlying disease, particularly the elderly and immunocompromised ones. In addition, S. agalactiae is a known fish pathogen, which compromises food safety and represents a zoonotic hazard. This study provides valuable structural, functional and evolutionary genomic information of a human S. agalactiae serotype Ia (ST-103) GBS85147 strain isolated from the oropharynx of an adult patient from Rio de Janeiro, thereby representing the first human isolate in Brazil. We used the Ion Torrent PGM platform with the 200 bp fragment library sequencing kit. The sequencing generated 578,082,183 bp, distributed among 2,973,022 reads, resulting in an approximately 246-fold mean coverage depth and was assembled using the Mira Assembler v3.9.18. The S. agalactiae strain GBS85147 comprises of a circular chromosome with a final genome length of 1,996,151 bp containing 1,915 protein-coding genes, 18 rRNA, 63 tRNA, 2 pseudogenes and a G + C content of 35.48 %.

Genome Sequence of Corynebacterium ulcerans Strain FRC11.

Here, we present the genome sequence of Corynebacterium ulcerans strain FRC11. The genome includes one circular chromosome of 2,442,826 bp (53.35% G+C content), and 2,210 genes were predicted, 2,146 of which are putative protein-coding genes, with 12 rRNAs and 51 tRNAs; 1 pseudogene was also identified.