Modeling and simulation analysis of propyl-thiouracil (PTU), an anti-thyroid drug on thyroid peroxidase (TPO), thyroid stimulating hormone receptor (TSHR), and sodium iodide (NIS) symporter based on systems biology approach

Modeling and simulation of atrazine biodegradation in bacteria and its effect in other living systems

Atrazine is the most commonly used herbicide worldwide in the agricultural system. The increased environmental concentration of the atrazine showed the toxic effects on the non-target living species. Biodegradation of the atrazine is possible with the bacterial systems. The present study investigated biodegradation potential of atrazine degrading bacteria and the impact of atrazine on environmental systems. Model of atrazine fate in ecological systems constructed using the cell designer. The used model further analyzed and simulated to know the biochemistry and physiology of the atrazine in different cellular networks. Topological analysis of the atrazine degradation confirmed the 289 nodes and 300 edges. Our results showed that the overall biomagnification of the atrazine in the different environmental systems. Atrazine is showing toxic effects on humans and plants, whereas degraded by the bacterial systems. To date, no one has analyzed the complete degradation and poisonous effects of the atrazine in the environment. Therefore, this study is useful for overall system biology based modeling and simulation analysis of atrazine in living systems.Communicated by Ramaswamy H. Sarma.

Insights into zinc-sensing metalloregulator 'Zur' deciphering mechanism of zinc transportation in Bacillus spp. by modeling, simulation and molecular docking

To comprehend the molecular mechanism of zinc transportation by bacteria tends to be a very complicated and time-consuming method. To date, fragmented and scanty studies are available about the mechanism of zinc transportation at molecular level. So, the present study scrutinizes in silico pathways of zinc fractions transportation, specifically in Bacillus spp. stimulating dynamic performance of zinc. For this, the constructed model reveals Zur to be the prime regulatory transport protein maintaining bacterial survivability at fluctuation in zinc concentrations, thereby attaining zinc homeostasis. Topology for hub nodes displays appropriate evidence of the molecular basis of bacterial zinc imports and exports. Further, the molecular docking reveals interaction of Zur protein with the zinc ligands (ZnCO₃ and ZnSO₄). By validation of binding affinity, binding energy and docking score via Autodock Vina and X-Score, the ZnSO₄ compound was found to possess excellent stability in the active pocket site of Zur, stating Zur-ZnSO₄ complex to be the most potential. Owing to which, the Zur-ZnSO₄ complex was selected and subjected to molecular dynamics simulation, revealing RMSD, RG, RMSF, SASA and interaction energy for 20 ns trajectory period. Henceforth,the study provides novel insight into revealing the unrecognized Zur protein pathway, assisting zinc transportation, besides retaining best interaction with ZnSO₄ ligand. This is the first system biology where molecular docking and molecular dynamics simulation-based investigation decipher the role of Zur transport protein system and interaction of its amino acids with zinc ligands in a simpler and economical form via in silico techniques.Communicated by Ramaswamy H. Sarma.

Biomolecules of 2-Thiouracil, 4-Thiouracil and 2,4-Dithiouracil: A DFT Study of the Hydration, Molecular Docking and Effect in DNA:RNAMicrohelixes

The molecular structure of 2-thiouracil, 4-thiouracil and 2,4-dithiouracil was analyzed under the effect of the first and second hydration shell by using the B3LYP density functional (DFT) method, and the results were compared to those obtained for the uracil molecule. A slight difference in the water distribution appears in these molecules. On the hydration of these molecules several trends in bond lengths and atomic charges were established. The ring in uracil molecule appears easier to be deformed and adapted to different environments as compared to that when it is thio-substituted. Molecular docking calculations of 2-thiouracil against three different pathogens: Bacillus subtilis, Escherichia coli and Candida albicans were carried out. Docking calculations of 2,4-dithiouracil ligand with various targeted proteins were also performed. Different DNA: RNA hybrid microhelixes with uridine, 2-thiouridine, 4-thiouridine and 2,4-dithiouridine nucleosides were optimized in a simple model with three nucleotide base pairs. Two main types of microhelixes were analyzed in detail depending on the intramolecular H-bond of the 2'-OH group. The weaker Watson-Crick (WC) base pair formed with thio-substituted uracil than with unsubstituted ones slightly deforms the helical and backbone parameters, especially with 2,4-dithiouridine. However, the thio-substitution significantly increases the dipole moment of the A-type microhelixes, as well as the rise and propeller twist parameters.

Modelling of the methyl halide biodegradation in bacteria and its effect on environmental systems

Methyl halide group of pesticides are being used widely in past decades as fumigant but due to their hazardous effect, these pesticides are not sold directly. They are volatile and gaseous in nature and may easily come in the contact of trophosphere and stratosphere. In troposphere, they are harmful to the living beings; nevertheless, in stratosphere they react with ozone and degrade the ozone layers. In this study, we have investigated the in-silico pathways of methyl halide and its toxic effect on living systems like pest, humans and environment. Till date, limited studies provide the understanding of degradation of methyl halide and its effect on the environment. This leads to availability of scanty information for overall bio-magnifications of methyl halides at molecular and cellular level. The model developed in the present study explains how a volatile toxic compound not only affects living systems on earth but also on environmental layers. Hub nodes were also evaluated by investigating the developed model topologically. Methyl transferase system is identified as promising enzyme in response to degradation of methyl halides.

Modeling of the jasmonate signaling pathway in Arabidopsis thaliana with respect to pathophysiology of Alternaria blight in Brassica

The productivity of Oilseed Brassica, one of the economically important crops of India, is seriously affected by the disease, Alternaria blight. The disease is mainly caused by two major necrotrophic fungi, Alternaria brassicae and Alternaria brassicicola which are responsible for significant yield losses. Till date, no resistant source is available against Alternaria blight, hence plant breeding methods can not be used to develop disease resistant varieties. Jasmonate mediated signalling pathway, which is known to play crucial role during defense response against necrotrophs, could be strengthened in Brassica plants to combat the disease. Since scanty information is available in Brassica-Alternaria pathosystems at molecular level therefore, in the present study efforts have been made to model jasmonic acid pathway in Arabidopsis thaliana to simulate the dynamic behaviour of molecular species in the model. Besides, the developed model was also analyzed topologically for investigation of the hubs node. COI1 is identified as one of the promising candidate genes in response to Alternaria and other linked components of plant defense mechanisms against the pathogens. The findings from present study are therefore informative for understanding the molecular basis of pathophysiology and rational management of Alternaria blight for securing food and nutritional security.

ZBIT Bioinformatics Toolbox: A Web-Platform for Systems Biology and Expression Data Analysis

Bioinformatics analysis has become an integral part of research in biology. However, installation and use of scientific software can be difficult and often requires technical expert knowledge. Reasons are dependencies on certain operating systems or required third-party libraries, missing graphical user interfaces and documentation, or nonstandard input and output formats. In order to make bioinformatics software easily accessible to researchers, we here present a web-based platform. The Center for Bioinformatics Tuebingen (ZBIT) Bioinformatics Toolbox provides web-based access to a collection of bioinformatics tools developed for systems biology, protein sequence annotation, and expression data analysis. Currently, the collection encompasses software for conversion and processing of community standards SBML and BioPAX, transcription factor analysis, and analysis of microarray data from transcriptomics and proteomics studies. All tools are hosted on a customized Galaxy instance and run on a dedicated computation cluster. Users only need a web browser and an active internet connection in order to benefit from this service. The web platform is designed to facilitate the usage of the bioinformatics tools for researchers without advanced technical background. Users can combine tools for complex analyses or use predefined, customizable workflows. All results are stored persistently and reproducible. For each tool, we provide documentation, tutorials, and example data to maximize usability. The ZBIT Bioinformatics Toolbox is freely available at https://webservices.cs.uni-tuebingen.de/.

SBMLsqueezer 2: context-sensitive creation of kinetic equations in biochemical networks

BACKGROUND

The size and complexity of published biochemical network reconstructions are steadily increasing, expanding the potential scale of derived computational models. However, the construction of large biochemical network models is a laborious and error-prone task. Automated methods have simplified the network reconstruction process, but building kinetic models for these systems is still a manually intensive task. Appropriate kinetic equations, based upon reaction rate laws, must be constructed and parameterized for each reaction. The complex test-and-evaluation cycles that can be involved during kinetic model construction would thus benefit from automated methods for rate law assignment.

RESULTS

We present a high-throughput algorithm to automatically suggest and create suitable rate laws based upon reaction type according to several criteria. The criteria for choices made by the algorithm can be influenced in order to assign the desired type of rate law to each reaction. This algorithm is implemented in the software package SBMLsqueezer 2. In addition, this program contains an integrated connection to the kinetics database SABIO-RK to obtain experimentally-derived rate laws when desired.

CONCLUSIONS

The described approach fills a heretofore absent niche in workflows for large-scale biochemical kinetic model construction. In several applications the algorithm has already been demonstrated to be useful and scalable. SBMLsqueezer is platform independent and can be used as a stand-alone package, as an integrated plugin, or through a web interface, enabling flexible solutions and use-case scenarios.

A comprehensive metabolic modeling of thyroid pathway in relation to thyroid pathophysiology and therapeutics

The thyroid pathway represents a complex interaction of different glands for thyroid hormone synthesis. Thyrotropin releasing hormone is synthesized in the hypothalamus and regulates thyrotropin stimulating hormone gene expression in the pituitary gland. In order to understand the complexity of the thyroid pathways, and using experimental data retrieved from the biomedical literature (e.g., NCBI, HuGE Navigator, Protein Data Bank, and KEGG), we constructed a metabolic map of the thyroid hormone pathway at a molecular level and analyzed it topologically. A total of five hub nodes were predicted in regards to the transcription thyroid receptor (TR), cAMP response element-binding protein (CREB), signal transducer and activator of transcription 3 (STAT 3), nuclear factor kappa-light-chain-enhancer of activated B cells (NFkB), and activator protein 1 (AP-1) as being potentially important in study of thyroid disorders and as novel putative therapeutic drug targets. Notably, the thyroid receptor is a highly connected hub node and currently used as a therapeutic target in hypothyroidism. Our analysis represents the first comprehensive description of the thyroid pathway, which pertains to understanding the function of the protein and gene interaction networks. The findings from this study are therefore informative for pathophysiology and rational therapeutics of thyroid disorders.