Refined template selection and combination algorithm significantly improves template-based modeling accuracy

Computationally Decoding NudF Residues To Enhance the Yield of the DXP Pathway

Terpenoids form a large pool of highly diverse organic compounds possessing several economically important properties, including nutritional, aromatic, and pharmacological properties. The 1-deoxy-d-xylulose 5-phosphate (DXP) pathway's end enzyme, nuclear distribution protein (NudF), interacting with isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP), is critical for the synthesis of isoprenol/prenol/downstream compounds. The enzyme is yet to be thoroughly investigated to increase the overall yield of terpenoids in the Bacillus subtilis, which is widely used in industry and is generally regarded as a safe (GRAS) bacterium. The study aims to analyze the evolutionary conservation across the active site for mapping the key residues for mutagenesis studies. The 37-sequence data set, extracted from 103 Bacillus subtilis entries, shows a high phylogenetic divergence, and only six one-motif sequences ASB92783.1, ASB69297.1, ASB56714.1, AOR97677.1, AOL97023.1, and OAZ71765.1 show a monophyly relationship, unlike a complete polyphyly relationship between the other 31 three-motif sequences. Furthermore, only 47 of 179 residues of the representative sequence CUB50584.1 are observed to be significantly conserved. Docking analysis suggests a preferential bias of adenosine diphosphate (ADP)-ribose pyrophosphatase toward IPP, and a nearly threefold energetic difference is observed between IPP and DMAPP. The loops are hereby shown to play a regulatory role in guiding the promiscuity of NudF toward a specific ligand. Computational saturation mutagenesis of the seven hotspot residues identifies two key positions LYS78 and PHE116, orderly encoded within loop1 and loop7, majorly interacting with the ligands DMAPP and IPP, and their mutants K78I/K78L and PHE116D/PHE116E are found to stabilize the overall conformation. Molecular dynamics analysis shows that the IPP complex is significantly more stable than the DMAPP complex, and the NudF structure is very unstable. Besides showing a promiscuous binding of NudF with ligands, the analysis suggests its rate-limiting nature. The study would allow us to customize the metabolic load toward the synthesis of any of the downstream molecules. The findings would pave the way for the development of catalytically improved NudF mutants for the large-scale production of specific terpenoids with significant nutraceutical or commercial value.

Computational and synthetic biology approaches for the biosynthesis of antiviral and anticancer terpenoids from Bacillus subtilis

Recent advancements in medicinal research have identified several antiviral and anticancer terpenoids that are usually deployed as a source of flavor, fragrances and pharmaceuticals. Under the current COVID-19 pandemic conditions, natural therapeutics with least side effects are the need of the hour to save the patients, especially, which are pre-affected with other medical complications. Although, plants are the major sources of terpenoids; however, for the environmental concerns, the global interest has shifted to the biocatalytic production of molecules from microbial sources. The gram-positive bacterium Bacillus subtilis is a suitable host in this regard due to its GRAS (generally regarded as safe) status, ease in genetic manipulations and wide industrial acceptability. The B. subtilis synthesizes its terpenoid molecules from 1-deoxy-d-xylulose-5-phosphate (DXP) pathway, a common route in almost all microbial strains. Here, we summarize the computational and synthetic biology approaches to improve the production of terpenoid-based therapeutics from B. subtilis by utilizing DXP pathway. We focus on the in-silico approaches for screening the functionally improved enzyme-variants of the two crucial enzymes namely, the DXP synthase (DXS) and farnesyl pyrophosphate synthase (FPPS). The approaches for engineering the active sites are subsequently explained. It will be helpful to construct the functionally improved enzymes for the high-yield production of terpenoid-based anticancer and antiviral metabolites, which would help to reduce the cost and improve the availability of such therapeutics for the humankind.

Engineering a multi epitope vaccine against SARS-CoV-2 by exploiting its non structural and structural proteins

SARS-CoV-2, the causative agent behind the ongoing pandemic exhibits an enhanced potential for infection when compared to its related family members- the SARS-CoV and MERS-CoV; which have caused similar disease outbreaks in the past. The severity of the global health burden, increasing mortality rate and the emergent economic crisis urgently demands the development of next generation vaccines. Amongst such emergent next generation vaccines are the multi-epitope subunit vaccines, which hold promise in combating deadly pathogens. In this study we have exploited immunoinformatics applications to delineate a vaccine candidate possessing multiple B and T cells epitopes by utilizing the SARS-CoV-2 non structural and structural proteins. The antigenicity potential, safety, structural stability and the production feasibility of the designed construct was evaluated computationally. Furthermore, due to the known role of human TLR-3 immune receptor in viral sensing, which facilitates host cells activation for an immune response, the vaccine construct was examined for its binding efficiency using molecular docking and molecular dynamics simulation studies, which resulted in strong and stable interactions. Finally, the immune simulation studies suggested an effective immune response on vaccine administration. Overall, the immunoinformatics analysis advocates that the proposed vaccine candidate is safe and immunogenic and therefore can be pushed as a lead for in vitro and in vivo investigations.Communicated by Ramaswamy H. Sarma.

Probabilistic divergence of a template-based modelling methodology from the ideal protocol

Protein structural information is essential for the detailed mapping of a functional protein network. For a higher modelling accuracy and quicker implementation, template-based algorithms have been extensively deployed and redefined. The methods only assess the predicted structure against its native state/template and do not estimate the accuracy for each modelling step. A divergence measure is therefore postulated to estimate the modelling accuracy against its theoretical optimal benchmark. By freezing the domain boundaries, the divergence measures are predicted for the most crucial steps of a modelling algorithm. To precisely refine the score using weighting constants, big data analysis could further be deployed.

Excavating the functionally crucial active-site residues of the DXS protein of Bacillus subtilis by exploring its closest homologues

Background

To achieve a high yield of terpenoid-based therapeutics, 1-deoxy-d-xylulose-5-phosphate (DXP) pathway has been significantly exploited for the production of downstream enzymes. The DXP synthase (DXS) enzyme, the initiator of this pathway, is pivotal for the convergence of carbon flux, and is computationally studied well for the industrially utilized generally regarded as safe (GRAS) bacterium Bacillus subtilis to decode its vital regions for aiding the construction of a functionally improved mutant library.

Results

For the 546 sequence dataset of DXS sequences, a representative set of 108 sequences is created, and it shows a significant evolutionary divergence across different species clubbed into 37 clades, whereas three clades are observed for the 76 sequence dataset of Bacillus subtilis. The DXS enzyme, sharing a statistically significant homology to transketolase, is shown to be evolutionarily too distant. By the mutual information-based co-evolutionary network and hotspot analysis, the most crucial loci within the active site are deciphered. The 650-residue representative structure displays a complete conservation of 114 loci, and only two co-evolving residues ASP154 and ILE371 are found to be the conserved ones. Lastly, P318D is predicted to be the top-ranked mutation causing the increase in the thermodynamic stability of 6OUW.

Conclusion

The study excavates the vital functional, phylogenetic, and conserved residues across the active site of the DXS protein, the key rate-limiting controller of the entire pathway. It would aid to computationally understand the evolutionary landscape of this industrially useful enzyme and would allow us to widen its substrate repertoire to increase the enzymatic yield of unnatural molecules for in vivo and in vitro applications.

Decoding the vital segments in human ATP-dependent RNA helicase

An analysis of the ATP-dependent RNA helicase using known functionally close analogs helps disclose the structural and functional information of the enzyme. The enzyme plays several interlinked biological functions and there is an urgent need to interpret its key active-site residues to infer function and establish role. The human protein q96c10.1 is annotated using tools such as interpro, go and cdd. The physicochemical properties are estimated using the tool protparam. We describe the enzyme protein model developed using modeller to identify active site residues. We used consurf to estimate the structural conservation and is evolutionary relationship is inferred using known close sequence homologs. The active site is predicted using castp and its topological flexibility is estimated through cabs-flex. The protein is annotated as a hydrolase using available data and ddx58 is found as its top-ranked interacting protein partner. We show that about 124 residues are found to be highly conserved among 259 homologs, clustered in 7 clades with the active-site showing low sequence conservation. It is further shown that only 9 loci among the 42 active-site residues are conserved with limited structural fluctuation from the wild type structure. Thus, we document various useful information linked to function, sequence similarity and phylogeny of the enzyme for annotation as potential helicase as designated by uniprot. Data shows limited degree of conserved sequence segments with topological flexibility unlike in other subfamily members of the protein.