Journal of Bioinformatics and Computational Biology最新文献

Motivation: The synthesis of proteins with novel desired properties is challenging but sought after by the industry and academia. The dominating approach is based on trial-and-error inducing point mutations, assisted by structural information or predictive models built with paired data that are difficult to collect. This study proposes a sequence-based unpaired-sample of novel protein inventor (SUNI) to build ThermalProGAN for generating thermally stable proteins based on sequence information.

Results: The ThermalProGAN can strongly mutate the input sequence with a median number of 32 residues. A known normal protein, 1RG0, was used to generate a thermally stable form by mutating 51 residues. After superimposing the two structures, high similarity is shown, indicating that the basic function would be conserved. Eighty four molecular dynamics simulation results of 1RG0 and the COVID-19 vaccine candidates with a total simulation time of 840[Formula: see text]ns indicate that the thermal stability increased.

Conclusion: This proof of concept demonstrated that transfer of a desired protein property from one set of proteins is feasible. Availability and implementation: The source code of ThermalProGAN can be freely accessed at https://github.com/markliou/ThermalProGAN/ with an MIT license. The website is https://thermalprogan.markliou.tw:433. Supplementary information: Supplementary data are available on Github.

Some prediction methods use probability to rank their predictions, while some other prediction methods do not rank their predictions and instead use [Formula: see text]-values to support their predictions. This disparity renders direct cross-comparison of these two kinds of methods difficult. In particular, approaches such as the Bayes Factor upper Bound (BFB) for [Formula: see text]-value conversion may not make correct assumptions for this kind of cross-comparisons. Here, using a well-established case study on renal cancer proteomics and in the context of missing protein prediction, we demonstrate how to compare these two kinds of prediction methods using two different strategies. The first strategy is based on false discovery rate (FDR) estimation, which does not make the same naïve assumptions as BFB conversions. The second strategy is a powerful approach which we colloquially call "home ground testing". Both strategies perform better than BFB conversions. Thus, we recommend comparing prediction methods by standardization to a common performance benchmark such as a global FDR. And where this is not possible, we recommend reciprocal "home ground testing".

Lung adenocarcinoma (LUAD) seriously threatens human health and generally results from dysfunction of relevant module molecules, which dynamically change with time and conditions, rather than that of an individual molecule. In this study, a novel network construction algorithm for identifying early warning network signals (IEWNS) is proposed for improving the performance of LUAD early diagnosis. To this end, we theoretically derived a dynamic criterion, namely, the relationship of variation (RV), to construct dynamic networks. RV infers correlation [Formula: see text] statistics to measure dynamic changes in molecular relationships during the process of disease development. Based on the dynamic networks constructed by IEWNS, network warning signals used to represent the occurrence of LUAD deterioration can be defined without human intervention. IEWNS was employed to perform a comprehensive analysis of gene expression profiles of LUAD from The Cancer Genome Atlas (TCGA) database and the Gene Expression Omnibus (GEO) database. The experimental results suggest that the potential biomarkers selected by IEWNS can facilitate a better understanding of pathogenetic mechanisms and help to achieve effective early diagnosis of LUAD. In conclusion, IEWNS provides novel insight into the initiation and progression of LUAD and helps to define prospective biomarkers for assessing disease deterioration.

Gel zymography quantifies the activity of certain enzymes in tumor processes. These enzymes are widely used in medical diagnosis. In order to analyze them, experts classify the zymography spots into various classes according to their tonalities. This classification is done by visual analysis, which is what makes it a subjective process. This work proposes a methodology to carry out this classifications with a process that involves an unsupervised learning algorithm in the images, denoted as the GI algorithm. With the experiments shown in this paper, this methodology could constitute a tool that bioinformatics scientists can trust to perform the desired classification since it is a quantitative indicator to order the enzymatic activity of the spots in a zymography.

In this study, we propose Feedback-AVPGAN, a system that aims to computationally generate novel antiviral peptides (AVPs). This system relies on the key premise of the Generative Adversarial Network (GAN) model and the Feedback method. GAN, a generative modeling approach that uses deep learning methods, comprises a generator and a discriminator. The generator is used to generate peptides; the generated proteins are fed to the discriminator to distinguish between the AVPs and non-AVPs. The original GAN design uses actual data to train the discriminator. However, not many AVPs have been experimentally obtained. To solve this problem, we used the Feedback method to allow the discriminator to learn from the existing as well as generated synthetic data. We implemented this method using a classifier module that classifies each peptide sequence generated by the GAN generator as AVP or non-AVP. The classifier uses the transformer network and achieves high classification accuracy. This mechanism enables the efficient generation of peptides with a high probability of exhibiting antiviral activity. Using the Feedback method, we evaluated various algorithms and their performance. Moreover, we modeled the structure of the generated peptides using AlphaFold2 and determined the peptides having similar physicochemical properties and structures to those of known AVPs, although with different sequences.

We present hybrid system-based gene regulatory network models for lambda, HK022, and Mu bacteriophages together with dynamics analysis of the modeled networks. The proposed lambda phage model LPH2 is based on an earlier work and incorporates more recent biological assumptions about the underlying gene regulatory mechanism, HK022, and Mu phage models are new. All three models provide accurate representations of experimentally observed lytic and lysogenic behavioral cycles. Importantly, the models also imply that lysis and lysogeny are the only stable behaviors that can occur in the modeled networks. In addition, the models allow to derive switching conditions that irrevocably lead to either lytic or lysogenic behavioral cycle as well as constraints that are required for their biological feasibility. For LPH2 model the feasibility constraints place two mutually independent requirements on comparative order of cro and cI protein binding site affinities. However, HK022 model, while broadly similar, does not require any of these constraints. Biologically very different lysis-lysogeny switching mechanism of Mu phage is also accurately reproduced by its model. In general the results show that hybrid system model (HSM) hybrid system framework can be successfully applied to modeling small ([Formula: see text] gene) regulatory networks and used for comprehensive analysis of model dynamics and stable behavior regions.

The profound impact of in silico studies for a fast-paced drug discovery pipeline is undeniable for pharmaceutical community. The rational design of novel drug candidates necessitates considering optimization of their different aspects prior to synthesis and biological evaluations. The affinity prediction of small ligands to target of interest for rank-ordering the potential ligands is one of the most routinely used steps in the context of virtual screening. So, the end-point methods were employed for binding free energy estimation focusing on evaluating simulation time effect. Then, a set of human aldose reductase inhibitors were selected for molecular dynamics (MD)-based binding free energy calculations. A total of 100[Formula: see text]ns MD simulation time was conducted for the ligand-receptor complexes followed by prediction of binding free energies using MM/PB(GB)SA and LIE approaches under different simulation time. The results revealed that a maximum of 30[Formula: see text]ns simulation time is sufficient for determination of binding affinities inferred from steady trend of squared correlation values (R²) between experimental and predicted [Formula: see text]G as a function of MD simulation time. In conclusion, the MM/PB(GB)SA algorithms performed well in terms of binding affinity prediction compared to LIE approach. The results provide new insights for large-scale applications of such predictions in an affordable computational cost.