Identification of a reliable sacral-sparing examination to assess the ASIA impairment scale in patients with traumatic spinal cord injury

OBJECTIVES

We evaluated the time course of the American Spinal Cord Injury Association (ASIA) impairment scale (AIS) for up to three months in participants within 72 h after traumatic spinal cord injury (TSCI) with complete paralysis. We aimed to determine the most useful sacral-sparing examination (deep anal pressure [DAP], voluntary anal contraction [VAC], S4-5 light touch [LT], or pin prick [PP] sensation) in determining AIS grades.

DESIGN

Retrospective cohort study.

SETTING

Spinal Injuries Center, Fukuoka, Japan.

PARTICIPANTS

Among 668 TSCI participants registered in the Japan Single Center study for Spinal Cord Injury Database (JSSCI-DB) between January 2012 and May 2020, we extracted the data of 80 patients with AIS grade A within 72 h after injury and neurological level of injury (NLI) at T12 or higher.

INTERVENTIONS

None.

OUTCOME MEASURES

The sacral-sparing examination at the time of the change to incomplete paralysis was compared to the AIS determination using a standard algorithm and with each assessment including the VAC, DAP, S4-5LT, and S4-5PP examinations at the time of AIS functional change. Agreement among assessments was evaluated using weighted kappa coefficients. The relationship was evaluated using Spearman's rank correlation coefficients.

RESULTS

Fifteen participants (18.8%) improved to incomplete paralysis (AIS B to D) within three months after injury. The single assessment among the sacral-sparing examinations with the highest agreement and strongest correlation with AIS determination was the S4-5LT examination (k = 0.89, P < 0.01, r = 0.84, P < 0.01).

CONCLUSIONS

The S4-5LT examination is key in determining complete or incomplete paralysis due to its high discriminatory power.

Stack-DHUpred: Advancing the accuracy of dihydrouridine modification sites detection via stacking approach

Dihydrouridine (DHU, D) is one of the most abundant post-transcriptional uridine modifications found in tRNA, mRNA, and snoRNA, closely associated with disease pathogenesis and various biological processes in eukaryotes. Identifying D sites is important for understanding the modification mechanisms and/or epigenetic regulation. However, biological experiments for detecting D sites are time-consuming and expensive. Given these challenges, computational methods have been developed for accurately identifying the D sites in genome-wide datasets. However, existing methods have some limitations, and their prediction performance needs to be improved. In this work, we have developed a new computational predictor for accurately identifying D sites called Stack-DHUpred. Briefly, we trained 66 baseline models or single-feature models by connecting six machine learning classifiers with eleven different feature encoding methods and stacked different baseline models to build stacked ensemble learning models. Subsequently, the optimal combination of the baseline models was identified for the construction of the final stacked model. Remarkably, the Stack-DHUpred outperformed the existing predictors on our new independent dataset, indicating that the stacking approach significantly improved the prediction performance. We have made Stack-DHUpred available to the public through a web server (http://kurata35.bio.kyutech.ac.jp/Stack-DHUpred) and a standalone program (https://github.com/kuratahiroyuki/Stack-DHUpred). We believe that Stack-DHUpred will be a valuable tool for accelerating the discovery of D modifications and understanding their role in post-transcriptional regulation.

Elucidating Key Characteristics of PFAS Binding to Human Peroxisome Proliferator-Activated Receptor Alpha: An Explainable Machine Learning Approach

Per- and polyfluoroalkyl substances (PFAS) are widely employed anthropogenic fluorinated chemicals known to disrupt hepatic lipid metabolism by binding to human peroxisome proliferator-activated receptor alpha (PPARα). Therefore, screening for PFAS that bind to PPARα is of critical importance. Machine learning approaches are promising techniques for rapid screening of PFAS. However, traditional machine learning approaches lack interpretability, posing challenges in investigating the relationship between molecular descriptors and PPARα binding. In this study, we aimed to develop a novel, explainable machine learning approach to rapidly screen for PFAS that bind to PPARα. We calculated the PPARα-PFAS binding score and 206 molecular descriptors for PFAS. Through systematic and objective selection of important molecular descriptors, we developed a machine learning model with good predictive performance using only three descriptors. The molecular size (b_single) and electrostatic properties (BCUT_PEOE_3 and PEOE_VSA_PPOS) are important for PPARα-PFAS binding. Alternative PFAS are considered safer than their legacy predecessors. However, we found that alternative PFAS with many carbon atoms and ether groups exhibited a higher affinity for PPARα. Therefore, confirming the toxicity of these alternative PFAS compounds with such characteristics through biological experiments is important.

Automatic Generation of SBML Kinetic Models from Natural Language Texts Using GPT

Kinetic modeling is an essential tool in systems biology research, enabling the quantitative analysis of biological systems and predicting their behavior. However, the development of kinetic models is a complex and time-consuming process. In this article, we propose a novel approach called KinModGPT, which generates kinetic models directly from natural language text. KinModGPT employs GPT as a natural language interpreter and Tellurium as an SBML generator. We demonstrate the effectiveness of KinModGPT in creating SBML kinetic models from complex natural language descriptions of biochemical reactions. KinModGPT successfully generates valid SBML models from a range of natural language model descriptions of metabolic pathways, protein-protein interaction networks, and heat shock response. This article demonstrates the potential of KinModGPT in kinetic modeling automation.

Decomposition of Unpolarized Fluorescence Spectrum of Uniaxially Oriented 1,3,5-Triphenylbenzene Microcrystals Into Polarized Fluorescence Spectra

Luminescence from solids such as crystals and aggregates is of growing academic and industrial interest. In this study, we report decomposition of the unpolarized fluorescence spectrum of uniaxially oriented 1,3,5-triphenylbenzene (TPB) microcrystals into four polarized spectra measured with polarizer (V: vertical and H: horizontal) and analyser (V: vertical and H: horizontal), where V and H indicate perpendicular and parallel to the layer of TPB molecules in the crystal, respectively. Resolved spectra were interpreted in terms of the molecular and excimer like (J- and H-dimer) emissions. The origin of the excimer like emissions was discussed in relation to the molecular packing in the crystal. It was shown that polarized crystal fluorescence can provide insight into the excitation/emission process in the crystal. Although preliminary, this study demonstrates the potential of polarized fluorescence to elucidate the luminescent mechanism.

CNN6mA: Interpretable neural network model based on position-specific CNN and cross-interactive network for 6mA site prediction

N6-methyladenine (6mA) plays a critical role in various epigenetic processing including DNA replication, DNA repair, silencing, transcription, and diseases such as cancer. To understand such epigenetic mechanisms, 6 mA has been detected by high-throughput technologies on a genome-wide scale at single-base resolution, together with conventional methods such as immunoprecipitation, mass spectrometry and capillary electrophoresis, but these experimental approaches are time-consuming and laborious. To complement these problems, we have developed a CNN-based 6 mA site predictor, named CNN6mA, which proposed two new architectures: a position-specific 1-D convolutional layer and a cross-interactive network. In the position-specific 1-D convolutional layer, position-specific filters with different window sizes were applied to an inquiry sequence instead of sharing the same filters over all positions in order to extract the position-specific features at different levels. The cross-interactive network explored the relationships between all the nucleotide patterns within the inquiry sequence. Consequently, CNN6mA outperformed the existing state-of-the-art models in many species and created the contribution score vector that intelligibly interpret the prediction mechanism. The source codes and web application in CNN6mA are freely accessible at https://github.com/kuratahiroyuki/CNN6mA.git and http://kurata35.bio.kyutech.ac.jp/CNN6mA/, respectively.

MLAGO: machine learning-aided global optimization for Michaelis constant estimation of kinetic modeling

BACKGROUND

Kinetic modeling is a powerful tool for understanding the dynamic behavior of biochemical systems. For kinetic modeling, determination of a number of kinetic parameters, such as the Michaelis constant (K_m), is necessary, and global optimization algorithms have long been used for parameter estimation. However, the conventional global optimization approach has three problems: (i) It is computationally demanding. (ii) It often yields unrealistic parameter values because it simply seeks a better model fitting to experimentally observed behaviors. (iii) It has difficulty in identifying a unique solution because multiple parameter sets can allow a kinetic model to fit experimental data equally well (the non-identifiability problem).

RESULTS

To solve these problems, we propose the Machine Learning-Aided Global Optimization (MLAGO) method for K_m estimation of kinetic modeling. First, we use a machine learning-based K_m predictor based only on three factors: EC number, KEGG Compound ID, and Organism ID, then conduct a constrained global optimization-based parameter estimation by using the machine learning-predicted K_m values as the reference values. The machine learning model achieved relatively good prediction scores: RMSE = 0.795 and R² = 0.536, making the subsequent global optimization easy and practical. The MLAGO approach reduced the error between simulation and experimental data while keeping K_m values close to the machine learning-predicted values. As a result, the MLAGO approach successfully estimated K_m values with less computational cost than the conventional method. Moreover, the MLAGO approach uniquely estimated K_m values, which were close to the measured values.

CONCLUSIONS

MLAGO overcomes the major problems in parameter estimation, accelerates kinetic modeling, and thus ultimately leads to better understanding of complex cellular systems. The web application for our machine learning-based K_m predictor is accessible at https://sites.google.com/view/kazuhiro-maeda/software-tools-web-apps , which helps modelers perform MLAGO on their own parameter estimation tasks.

ICAN: Interpretable cross-attention network for identifying drug and target protein interactions

Drug-target protein interaction (DTI) identification is fundamental for drug discovery and drug repositioning, because therapeutic drugs act on disease-causing proteins. However, the DTI identification process often requires expensive and time-consuming tasks, including biological experiments involving large numbers of candidate compounds. Thus, a variety of computation approaches have been developed. Of the many approaches available, chemo-genomics feature-based methods have attracted considerable attention. These methods compute the feature descriptors of drugs and proteins as the input data to train machine and deep learning models to enable accurate prediction of unknown DTIs. In addition, attention-based learning methods have been proposed to identify and interpret DTI mechanisms. However, improvements are needed for enhancing prediction performance and DTI mechanism elucidation. To address these problems, we developed an attention-based method designated the interpretable cross-attention network (ICAN), which predicts DTIs using the Simplified Molecular Input Line Entry System of drugs and amino acid sequences of target proteins. We optimized the attention mechanism architecture by exploring the cross-attention or self-attention, attention layer depth, and selection of the context matrixes from the attention mechanism. We found that a plain attention mechanism that decodes drug-related protein context features without any protein-related drug context features effectively achieved high performance. The ICAN outperformed state-of-the-art methods in several metrics on the DAVIS dataset and first revealed with statistical significance that some weighted sites in the cross-attention weight matrix represent experimental binding sites, thus demonstrating the high interpretability of the results. The program is freely available at https://github.com/kuratahiroyuki/ICAN.

iACVP: markedly enhanced identification of anti-coronavirus peptides using a dataset-specific word2vec model

The COVID-19 pandemic caused several million deaths worldwide. Development of anti-coronavirus drugs is thus urgent. Unlike conventional non-peptide drugs, antiviral peptide drugs are highly specific, easy to synthesize and modify, and not highly susceptible to drug resistance. To reduce the time and expense involved in screening thousands of peptides and assaying their antiviral activity, computational predictors for identifying anti-coronavirus peptides (ACVPs) are needed. However, few experimentally verified ACVP samples are available, even though a relatively large number of antiviral peptides (AVPs) have been discovered. In this study, we attempted to predict ACVPs using an AVP dataset and a small collection of ACVPs. Using conventional features, a binary profile and a word-embedding word2vec (W2V), we systematically explored five different machine learning methods: Transformer, Convolutional Neural Network, bidirectional Long Short-Term Memory, Random Forest (RF) and Support Vector Machine. Via exhaustive searches, we found that the RF classifier with W2V consistently achieved better performance on different datasets. The two main controlling factors were: (i) the dataset-specific W2V dictionary was generated from the training and independent test datasets instead of the widely used general UniProt proteome and (ii) a systematic search was conducted and determined the optimal k-mer value in W2V, which provides greater discrimination between positive and negative samples. Therefore, our proposed method, named iACVP, consistently provides better prediction performance compared with existing state-of-the-art methods. To assist experimentalists in identifying putative ACVPs, we implemented our model as a web server accessible via the following link: http://kurata35.bio.kyutech.ac.jp/iACVP.

Deepm5C: A deep learning-based hybrid framework for identifying human RNA N5-methylcytosine sites using a stacking strategy

As one of the most prevalent post-transcriptional epigenetic modifications, N5-methylcytosine (m5C), plays an essential role in various cellular processes and disease pathogenesis. Therefore, it is important accurately identify m5C modifications in order to gain a deeper understanding of cellular processes and other possible functional mechanisms. Although a few computational methods have been proposed, their respective models have been developed using small training datasets. Hence, their practical application is quite limited in genome-wide detection. To overcome the existing limitations, we propose Deepm5C, a bioinformatics method to identify RNA m5C sites in the throughout human genome. To develop Deepm5C, we constructed a novel benchmarking dataset and investigated a mixture of three conventional feature encoding algorithms and a feature derived from word embedding approaches. Afterwards, four variants of deep learning classifiers and four commonly used conventional classifiers were employed and trained with the four encodings, ultimately obtaining 32 baseline models. A stacking strategy is effectively utilized by integrating the predicted output of the optimal baseline models and trained with a 1-D convolutional neural network. As a result, the Deepm5C predictor achieved excellent performance during cross-validation with a Matthews correlation coefficient and accuracy of 0.697 and 0.855, respectively. The corresponding metrics during the independent test were 0.691 and 0.852, respectively. Overall, Deepm5C achieved a more accurate and stable performance than the baseline models and significantly outperformed the existing predictors, demonstrating the effectiveness of our proposed hybrid framework. Furthermore, Deepm5C is expected to assist community-wide efforts in identifying putative m5Cs and formulate the novel testable biological hypothesis.

BERT6mA: prediction of DNA N6-methyladenine site using deep learning-based approaches

N6-methyladenine (6mA) is associated with important roles in DNA replication, DNA repair, transcription, regulation of gene expression. Several experimental methods were used to identify DNA modifications. However, these experimental methods are costly and time-consuming. To detect the 6mA and complement these shortcomings of experimental methods, we proposed a novel, deep leaning approach called BERT6mA. To compare the BERT6mA with other deep learning approaches, we used the benchmark datasets including 11 species. The BERT6mA presented the highest AUCs in eight species in independent tests. Furthermore, BERT6mA showed higher and comparable performance with the state-of-the-art models while the BERT6mA showed poor performances in a few species with a small sample size. To overcome this issue, pretraining and fine-tuning between two species were applied to the BERT6mA. The pretrained and fine-tuned models on specific species presented higher performances than other models even for the species with a small sample size. In addition to the prediction, we analyzed the attention weights generated by BERT6mA to reveal how the BERT6mA model extracts critical features responsible for the 6mA prediction. To facilitate biological sciences, the BERT6mA online web server and its source codes are freely accessible at https://github.com/kuratahiroyuki/BERT6mA.git, respectively.

Late onset of coronary ostial stenosis following surgical aortic valve replacement in a patient with anomalous origin of the right coronary artery

Surgical aortic valve replacement (SAVR) in patients with anomalous origination of a coronary artery from the opposite sinus is associated with risk for myocardial ischemia during the perioperative period. [1] However, iatrogenic coronary ostial stenosis (ICOS) generally occurs within the first 6 months after SAVR. We present an unusual case of a 74-year-old man with anomalous origination of the right coronary artery from the left coronary sinus, who developed effort angina due to ICOS 19 months following SAVR and ascending aorta replacement. Angiography and computed tomography were utilized to perform a comparison before and after the procedure. From the results, it was evident that the flattened mild stenosis preoperatively was caused by anomalous origination of a coronary artery from the opposite sinus and progressed to severe stenosis by ICOS after the procedure. The patient was successfully treated with percutaneous coronary intervention. <Learning objective: Surgical aortic valve replacement in patients with anomalous origination of a coronary artery from the opposite sinus may lead to coronary ischemic events due to iatrogenic coronary ostial stenosis in the long term. Percutaneous coronary intervention is a reasonable treatment option for the patients with iatrogenic coronary ostial stenosis after surgical aortic valve replacement.>.

LSTM-PHV: prediction of human-virus protein-protein interactions by LSTM with word2vec

Viral infection involves a large number of protein-protein interactions (PPIs) between human and virus. The PPIs range from the initial binding of viral coat proteins to host membrane receptors to the hijacking of host transcription machinery. However, few interspecies PPIs have been identified, because experimental methods including mass spectrometry are time-consuming and expensive, and molecular dynamic simulation is limited only to the proteins whose 3D structures are solved. Sequence-based machine learning methods are expected to overcome these problems. We have first developed the LSTM model with word2vec to predict PPIs between human and virus, named LSTM-PHV, by using amino acid sequences alone. The LSTM-PHV effectively learnt the training data with a highly imbalanced ratio of positive to negative samples and achieved AUCs of 0.976 and 0.973 and accuracies of 0.984 and 0.985 on the training and independent datasets, respectively. In predicting PPIs between human and unknown or new virus, the LSTM-PHV learned greatly outperformed the existing state-of-the-art PPI predictors. Interestingly, learning of only sequence contexts as words is sufficient for PPI prediction. Use of uniform manifold approximation and projection demonstrated that the LSTM-PHV clearly distinguished the positive PPI samples from the negative ones. We presented the LSTM-PHV online web server and support data that are freely available at http://kurata35.bio.kyutech.ac.jp/LSTM-PHV.

Recent development of bioinformatics tools for microRNA target prediction

MicroRNAs (miRNAs) are central players that regulate the post-transcriptional processes of gene expression. Binding of miRNAs to target mRNAs can repress their translation by inducing the degradation or by inhibiting the translation of the target mRNAs. High-throughput experimental approaches for miRNA target identification are costly and time-consuming, depending on various factors. It is vitally important to develop the bioinformatics methods for accurately predicting miRNA targets. With the increase of RNA sequences in the post-genomic era, bioinformatics methods are being developed for miRNA studies specially for miRNA target prediction. This review summarizes the current development of state-of-the-art bioinformatics tools for miRNA target prediction, points out the progress and limitations of the available miRNA databases, and their working principles. Finally, we discuss the caveat and perspectives of the next-generation algorithms for the prediction of miRNA targets.

Fabrication and Characterisation of Organic EL Devices in the Presence of Cyclodextrin as an Interlayer

This study examined glass-based organic electroluminescence in the presence of a cyclodextrin polymer as an interlayer. Glass-based organic electroluminescence was achieved by the deposition of five layers of N,N'-Bis(3-methylphenyl)N,N'-bis(phenyl)-benzidine, cyclodextrin polymer (CDP), tris-(8-hydroxyquinolinato) aluminium LiF and Al on an indium tin oxide-coated glass substrate. The glass-based OEL exhibited green emission owing to the fluorescence of tris-(8-hydroxyquinolinato) aluminium. The highest luminance was 19,620 cd m-2. Moreover, the glass-based organic electroluminescence device showed green emission at 6 V in the curved state because of the inhibited aggregation of the cyclodextrin polymer. All organic molecules are insulating, but except CDP, they are standard molecules in conventional organic electroluminescence devices. In this device, the CDP layer contained pores that could allow conventional organic molecules to enter the pores and affect the organic electroluminescence interface. In particular, self-association was suppressed, efficiency was improved, and light emission was observed without the need for a high voltage. Overall, the glass-based organic electroluminescence device using CDP is an environmentally friendly device with a range of potential energy saving applications.

NeuroPred-FRL: an interpretable prediction model for identifying neuropeptide using feature representation learning

Neuropeptides (NPs) are the most versatile neurotransmitters in the immune systems that regulate various central anxious hormones. An efficient and effective bioinformatics tool for rapid and accurate large-scale identification of NPs is critical in immunoinformatics, which is indispensable for basic research and drug development. Although a few NP prediction tools have been developed, it is mandatory to improve their NPs' prediction performances. In this study, we have developed a machine learning-based meta-predictor called NeuroPred-FRL by employing the feature representation learning approach. First, we generated 66 optimal baseline models by employing 11 different encodings, six different classifiers and a two-step feature selection approach. The predicted probability scores of NPs based on the 66 baseline models were combined to be deemed as the input feature vector. Second, in order to enhance the feature representation ability, we applied the two-step feature selection approach to optimize the 66-D probability feature vector and then inputted the optimal one into a random forest classifier for the final meta-model (NeuroPred-FRL) construction. Benchmarking experiments based on both cross-validation and independent tests indicate that the NeuroPred-FRL achieves a superior prediction performance of NPs compared with the other state-of-the-art predictors. We believe that the proposed NeuroPred-FRL can serve as a powerful tool for large-scale identification of NPs, facilitating the characterization of their functional mechanisms and expediting their applications in clinical therapy. Moreover, we interpreted some model mechanisms of NeuroPred-FRL by leveraging the robust SHapley Additive exPlanation algorithm.

PredNTS: Improved and Robust Prediction of Nitrotyrosine Sites by Integrating Multiple Sequence Features

Nitrotyrosine, which is generated by numerous reactive nitrogen species, is a type of protein post-translational modification. Identification of site-specific nitration modification on tyrosine is a prerequisite to understanding the molecular function of nitrated proteins. Thanks to the progress of machine learning, computational prediction can play a vital role before the biological experimentation. Herein, we developed a computational predictor PredNTS by integrating multiple sequence features including K-mer, composition of k-spaced amino acid pairs (CKSAAP), AAindex, and binary encoding schemes. The important features were selected by the recursive feature elimination approach using a random forest classifier. Finally, we linearly combined the successive random forest (RF) probability scores generated by the different, single encoding-employing RF models. The resultant PredNTS predictor achieved an area under a curve (AUC) of 0.910 using five-fold cross validation. It outperformed the existing predictors on a comprehensive and independent dataset. Furthermore, we investigated several machine learning algorithms to demonstrate the superiority of the employed RF algorithm. The PredNTS is a useful computational resource for the prediction of nitrotyrosine sites. The web-application with the curated datasets of the PredNTS is publicly available.

PUP-Fuse: Prediction of Protein Pupylation Sites by Integrating Multiple Sequence Representations

Pupylation is a type of reversible post-translational modification of proteins, which plays a key role in the cellular function of microbial organisms. Several proteomics methods have been developed for the prediction and analysis of pupylated proteins and pupylation sites. However, the traditional experimental methods are laborious and time-consuming. Hence, computational algorithms are highly needed that can predict potential pupylation sites using sequence features. In this research, a new prediction model, PUP-Fuse, has been developed for pupylation site prediction by integrating multiple sequence representations. Meanwhile, we explored the five types of feature encoding approaches and three machine learning (ML) algorithms. In the final model, we integrated the successive ML scores using a linear regression model. The PUP-Fuse achieved a Mathew correlation value of 0.768 by a 10-fold cross-validation test. It also outperformed existing predictors in an independent test. The web server of the PUP-Fuse with curated datasets is freely available.

Critical evaluation of web-based DNA N6-methyladenine site prediction tools

Methylation of DNA N6-methyladenosine (6mA) is a type of epigenetic modification that plays pivotal roles in various biological processes. The accurate genome-wide identification of 6mA is a challenging task that leads to understanding the biological functions. For the last 5 years, a number of bioinformatics approaches and tools for 6mA site prediction have been established, and some of them are easily accessible as web application. Nevertheless, the accurate genome-wide identification of 6mA is still one of the challenging works that lead to understanding the biological functions. Especially in practical applications, these tools have implemented diverse encoding schemes, machine learning algorithms and feature selection methods, whereas few systematic performance comparisons of 6mA site predictors have been reported. In this review, 11 publicly available 6mA predictors evaluated with seven different species-specific datasets (Arabidopsis thaliana, Tolypocladium, Diospyros lotus, Saccharomyces cerevisiae, Drosophila melanogaster, Caenorhabditis elegans and Escherichia coli). Of those, few species are close homologs, and the remaining datasets are distant sequences. Our independent, validation tests demonstrated that Meta-i6mA and MM-6mAPred models for A. thaliana, Tolypocladium, S. cerevisiae and D. melanogaster achieved excellent overall performance when compared with their counterparts. However, none of the existing methods were suitable for E. coli, C. elegans and D. lotus. A feasibility of the existing predictors is also discussed for the seven species. Our evaluation provides useful guidelines for the development of 6mA site predictors and helps biologists selecting suitable prediction tools.

IRC-Fuse: improved and robust prediction of redox-sensitive cysteine by fusing of multiple feature representations

Redox-sensitive cysteine (RSC) thiol contributes to many biological processes. The identification of RSC plays an important role in clarifying some mechanisms of redox-sensitive factors; nonetheless, experimental investigation of RSCs is expensive and time-consuming. The computational approaches that quickly and accurately identify candidate RSCs using the sequence information are urgently needed. Herein, an improved and robust computational predictor named IRC-Fuse was developed to identify the RSC by fusing of multiple feature representations. To enhance the performance of our model, we integrated the probability scores evaluated by the random forest models implementing different encoding schemes. Cross-validation results exhibited that the IRC-Fuse achieved accuracy and AUC of 0.741 and 0.807, respectively. The IRC-Fuse outperformed exiting methods with improvement of 10% and 13% on accuracy and MCC, respectively, over independent test data. Comparative analysis suggested that the IRC-Fuse was more effective and promising than the existing predictors. For the convenience of experimental scientists, the IRC-Fuse online web server was implemented and publicly accessible at http://kurata14.bio.kyutech.ac.jp/IRC-Fuse/ .

iLBE for Computational Identification of Linear B-cell Epitopes by Integrating Sequence and Evolutionary Features

Linear B-cell epitopes are critically important for immunological applications, such as vaccine design, immunodiagnostic test, and antibody production, as well as disease diagnosis and therapy. The accurate identification of linear B-cell epitopes remains challenging despite several decades of research. In this work, we have developed a novel predictor, Identification of Linear B-cell Epitope (iLBE), by integrating evolutionary and sequence-based features. The successive feature vectors were optimized by a Wilcoxon-rank sum test. Then the random forest (RF) algorithm using the optimal consecutive feature vectors was applied to predict linear B-cell epitopes. We combined the RF scores by the logistic regression to enhance the prediction accuracy. iLBE yielded an area under curve score of 0.809 on the training dataset and outperformed other prediction models on a comprehensive independent dataset. iLBE is a powerful computational tool to identify the linear B-cell epitopes and would help to develop penetrating diagnostic tests. A web application with curated datasets for iLBE is freely accessible at http://kurata14.bio.kyutech.ac.jp/iLBE/.

Meta-i6mA: an interspecies predictor for identifying DNA N6-methyladenine sites of plant genomes by exploiting informative features in an integrative machine-learning framework

DNA N6-methyladenine (6mA) represents important epigenetic modifications, which are responsible for various cellular processes. The accurate identification of 6mA sites is one of the challenging tasks in genome analysis, which leads to an understanding of their biological functions. To date, several species-specific machine learning (ML)-based models have been proposed, but majority of them did not test their model to other species. Hence, their practical application to other plant species is quite limited. In this study, we explored 10 different feature encoding schemes, with the goal of capturing key characteristics around 6mA sites. We selected five feature encoding schemes based on physicochemical and position-specific information that possesses high discriminative capability. The resultant feature sets were inputted to six commonly used ML methods (random forest, support vector machine, extremely randomized tree, logistic regression, naïve Bayes and AdaBoost). The Rosaceae genome was employed to train the above classifiers, which generated 30 baseline models. To integrate their individual strength, Meta-i6mA was proposed that combined the baseline models using the meta-predictor approach. In extensive independent test, Meta-i6mA showed high Matthews correlation coefficient values of 0.918, 0.827 and 0.635 on Rosaceae, rice and Arabidopsis thaliana, respectively and outperformed the existing predictors. We anticipate that the Meta-i6mA can be applied across different plant species. Furthermore, we developed an online user-friendly web server, which is available at http://kurata14.bio.kyutech.ac.jp/Meta-i6mA/.

Evolution of Sequence-based Bioinformatics Tools for Protein-protein Interaction Prediction

Protein-protein interactions (PPIs) are the physical connections between two or more proteins via electrostatic forces or hydrophobic effects. Identification of the PPIs is pivotal, which contributes to many biological processes including protein function, disease incidence, and therapy design. The experimental identification of PPIs via high-throughput technology is time-consuming and expensive. Bioinformatics approaches are expected to solve such restrictions. In this review, our main goal is to provide an inclusive view of the existing sequence-based computational prediction of PPIs. Initially, we briefly introduce the currently available PPI databases and then review the state-of-the-art bioinformatics approaches, working principles, and their performances. Finally, we discuss the caveats and future perspective of the next generation algorithms for the prediction of PPIs.

Asymmetric Amination of meso-Epoxide with Vegetable Powder as a Low-Toxicity Catalyst

This paper describes the scope and limitation of substrates subjected to asymmetric amination with epoxides catalyzed by a soluble soybean polysaccharide (Soyafibe S-DN), which we recently discovered from the reaction of 1,2-epoxycyclohexane with cyclopropylamine. Various meso-epoxides reacted with various amines afforded the corresponding products with good enantiomeric selectivity. Since it was found that pectin was found to have a catalytic ability after screening commercially available polysaccharides, we studied 33 different vegetable powders having pectic substances, and we found that many vegetable powders showed catalytic ability. These results should guide in using vegetable components as low-toxic catalysts for the production of pharmaceuticals.

Coupling protocol of interlocked feedback oscillators in circadian clocks

Circadian rhythms (approx. 24 h) show the robustness of key oscillatory features such as phase, period and amplitude against external and internal variations. The robustness of Drosophila circadian clocks can be generated by interlocked transcriptional-translational feedback loops, where two negative feedback loops are coupled through mutual activations. The mechanisms by which such coupling protocols have survived out of many possible protocols remain to be revealed. To address this question, we investigated two distinct coupling protocols: activator-coupled oscillators (ACO) and repressor-coupled oscillators (RCO). We focused on the two coupling parameters: coupling dissociation constant and coupling time-delay. Interestingly, the ACO was able to produce anti-phase or morning-evening cycles, whereas the RCO produced in-phase ones. Deterministic and stochastic analyses demonstrated that the anti-phase ACO provided greater fluctuations in amplitude not only with respect to changes in coupling parameters but also to random parameter perturbations than the in-phase RCO. Moreover, the ACO deteriorated the entrainability to the day-night master clock, whereas the RCO produced high entrainability. Considering that the real, interlocked feedback loops have evolved as the ACO, instead of the RCO, we first proposed a hypothesis that the morning-evening or anti-phase cycle is more essential for Drosophila than achieving robustness and entrainability.

i6mA-Fuse: improved and robust prediction of DNA 6 mA sites in the Rosaceae genome by fusing multiple feature representation

DNA N⁶-methyladenine (6 mA) is one of the most vital epigenetic modifications and involved in controlling the various gene expression levels. With the avalanche of DNA sequences generated in numerous databases, the accurate identification of 6 mA plays an essential role for understanding molecular mechanisms. Because the experimental approaches are time-consuming and costly, it is desirable to develop a computation model for rapidly and accurately identifying 6 mA. To the best of our knowledge, we first proposed a computational model named i6mA-Fuse to predict 6 mA sites from the Rosaceae genomes, especially in Rosa chinensis and Fragaria vesca. We implemented the five encoding schemes, i.e., mononucleotide binary, dinucleotide binary, k-space spectral nucleotide, k-mer, and electron-ion interaction pseudo potential compositions, to build the five, single-encoding random forest (RF) models. The i6mA-Fuse uses a linear regression model to combine the predicted probability scores of the five, single encoding-based RF models. The resultant species-specific i6mA-Fuse achieved remarkably high performances with AUCs of 0.982 and 0.978 and with MCCs of 0.869 and 0.858 on the independent datasets of Rosa chinensis and Fragaria vesca, respectively. In the F. vesca-specific i6mA-Fuse, the MBE and EIIP contributed to 75% and 25% of the total prediction; in the R. chinensis-specific i6mA-Fuse, Kmer, MBE, and EIIP contribute to 15%, 65%, and 20% of the total prediction. To assist high-throughput prediction for DNA 6 mA identification, the i6mA-Fuse is publicly accessible at https://kurata14.bio.kyutech.ac.jp/i6mA-Fuse/.

Computer-Aided Rational Design of Efficient NADPH Production System by Escherichia coli pgi Mutant Using a Mixture of Glucose and Xylose

Lignocellulosic biomass can be hydrolyzed into two major sugars of glucose and xylose, and thus the strategy for the efficient consumption of both sugars is highly desirable. NADPH is the essential molecule for the production of industrially important value-added chemicals, and thus its availability is quite important. Escherichia coli mutant lacking the pgi gene encoding phosphoglucose isomerase (Pgi) has been preferentially used to overproduce the NADPH. However, there exists a disadvantage that the cell growth rate becomes low for the mutant grown on glucose. This limits the efficient NADPH production, and therefore, it is quite important to investigate how addition of different carbon source such as xylose (other than glucose) effectively improves the NADPH production. In this study, we have developed a kinetic model to propose an efficient NADPH production system using E. coli pgi-knockout mutant with a mixture of glucose and xylose. The proposed system adds xylose to glucose medium to recover the suppressed growth of the pgi mutant, and determines the xylose content to maximize the NADPH productivity. Finally, we have designed a mevalonate (MVA) production system by implementing ArcA overexpression into the pgi-knockout mutant using a mixture of glucose and xylose. In addition to NADPH overproduction, the accumulation of acetyl-CoA (AcCoA) is necessary for the efficient MVA production. In the present study, therefore, we considered to overexpress ArcA, where ArcA overexpression suppresses the TCA cycle, causing the overflow of AcCoA, a precursor of MVA. We predicted the xylose content that maximizes the MVA production. This approach demonstrates the possibility of a great progress in the computer-aided rational design of the microbial cell factories for useful metabolite production.

Platypnea-Orthodeoxia Syndrome Associated with Spontaneously Ruptured Chordae Tendineae of Tricuspid Valve

•

POS is a rare clinical condition.

•

Symptoms are upright posture dyspnea and arterial desaturation.

•

Supine and sitting position oxygen saturation levels suggest right-to-left shunt.

•

First report of uncommon spontaneous tricuspid valve chordal rupture–induced POS.

Recent Development of Machine Learning Methods in Microbial Phosphorylation Sites

A variety of protein post-translational modifications has been identified that control many cellular functions. Phosphorylation studies in mycobacterial organisms have shown critical importance in diverse biological processes, such as intercellular communication and cell division. Recent technical advances in high-precision mass spectrometry have determined a large number of microbial phosphorylated proteins and phosphorylation sites throughout the proteome analysis. Identification of phosphorylated proteins with specific modified residues through experimentation is often labor-intensive, costly and time-consuming. All these limitations could be overcome through the application of machine learning (ML) approaches. However, only a limited number of computational phosphorylation site prediction tools have been developed so far. This work aims to present a complete survey of the existing ML-predictors for microbial phosphorylation. We cover a variety of important aspects for developing a successful predictor, including operating ML algorithms, feature selection methods, window size, and software utility. Initially, we review the currently available phosphorylation site databases of the microbiome, the state-of-the-art ML approaches, working principles, and their performances. Lastly, we discuss the limitations and future directions of the computational ML methods for the prediction of phosphorylation.

Improvement of the memory function of a mutual repression network in a stochastic environment by negative autoregulation

BACKGROUND

Cellular memory is a ubiquitous function of biological systems. By generating a sustained response to a transient inductive stimulus, often due to bistability, memory is central to the robust control of many important biological processes. However, our understanding of the origins of cellular memory remains incomplete. Stochastic fluctuations that are inherent to most biological systems have been shown to hamper memory function. Yet, how stochasticity changes the behavior of genetic circuits is generally not clear from a deterministic analysis of the network alone. Here, we apply deterministic rate equations, stochastic simulations, and theoretical analyses of Fokker-Planck equations to investigate how intrinsic noise affects the memory function in a mutual repression network.

RESULTS

We find that the addition of negative autoregulation improves the persistence of memory in a small gene regulatory network by reducing stochastic fluctuations. Our theoretical analyses reveal that this improved memory function stems from an increased stability of the steady states of the system. Moreover, we show how the tuning of critical network parameters can further enhance memory.

CONCLUSIONS

Our work illuminates the power of stochastic and theoretical approaches to understanding biological circuits, and the importance of considering stochasticity when designing synthetic circuits with memory function.

i4mC-ROSE, a bioinformatics tool for the identification of DNA N4-methylcytosine sites in the Rosaceae genome

One of the most important epigenetic modifications is N4-methylcytosine, which regulates many biological processes including DNA replication and chromosome stability. Identification of N4-methylcytosine sites is pivotal to understand specific biological functions. Herein, we developed the first bioinformatics tool called i4mC-ROSE for identifying N4-methylcytosine sites in the genomes of Fragaria vesca and Rosa chinensis in the Rosaceae, which utilizes a random forest classifier with six encoding methods that cover various aspects of DNA sequence information. The i4mC-ROSE predictor achieves area under the curve scores of 0.883 and 0.889 for the two genomes during cross-validation. Moreover, the i4mC-ROSE outperforms other classifiers tested in this study when objectively evaluated on the independent datasets. The proposed i4mC-ROSE tool can serve users' demand for the prediction of 4mC sites in the Rosaceae genome. The i4mC-ROSE predictor and utilized datasets are publicly accessible at http://kurata14.bio.kyutech.ac.jp/i4mC-ROSE/.

Self-replenishment cycles generate a threshold response

Many metabolic cycles, including the tricarboxylic acid cycle, glyoxylate cycle, Calvin cycle, urea cycle, coenzyme recycling, and substrate cycles, are well known to catabolize and anabolize different metabolites for efficient energy and mass conversion. In terms of stoichiometric structure, this study explicitly identifies two types of metabolic cycles. One is the well-known, elementary cycle that converts multiple substrates into different products and recycles one of the products as a substrate, where the recycled substrate is supplied from the outside to run the cycle. The other is the self-replenishment cycle that merges multiple substrates into two or multiple identical products and reuses one of the products as a substrate. The substrates are autonomously supplied within the cycle. This study first defines the self-replenishment cycles that many scientists have overlooked despite its functional importance. Theoretical analysis has revealed the design principle of the self-replenishment cycle that presents a threshold response without any bistability nor cooperativity. To verify the principle, three detailed kinetic models of self-replenishment cycles embedded in an E. coli metabolic system were simulated. They presented the threshold response or digital switch-like function that steeply shift metabolic status.

P315In-vivo and vitro mitochondrial transfer from adipose-derived mesenchymal stem cell to ischemic cardiomyocyte

Background

Although transplantation of human Adipose-derived Mesenchymal stem cell (hADSC) shows efficacy in the treatment of ischemic cardiomyopathy, its therapeutic mechanisms have not been fully elucidated. It has been already reported that mitochondria transfer to recipient cells have impact on resistance to injury and tissue regeneration, however this phenomenon has not been elucidated in the damaged heart. Therefore, we hypothesized that ADSC transfer own mitochondria to cardiomyocytes in-vivo and in-vitro under ischemic condition, resulting in the functional recovery of cardiomyocyte.

Method and result

Transplantation of hADSC (group A) to the heart surface or sham operation (group C) was performed in rats that were subjected to LAD ligation 2 weeks prior to the treatment (n=10 each). The number of transplant cell was 1x106/body. Three days after transplantation, transferred hADSCs' mitochondria were observed in recipient cardiomyocytes histologically (Figure). Quantitative PCR analysis revealed that mitochondrial genome of recipient myocytes increased over time. The cardiac function assessed with echocardiography was significantly better in group A. Furthermore, live-imaging of hADSC transplantation revealed the suspected transfer of mitochondria to beating heart.

In-vitro, the co-culture of rat cardiomyocytes (rCM) and hADSC was observed with time-lapse photography and demonstrated mitochondrial transfer under the hypoxic condition. The measuring the oxygen consumption rate (OCR) of these cells showed that OCR of rCM was reinforced by co-culture with hADSC conspicuously.

Figure 1

Conclusion

Mitochondrial transfer from hADSC to rCM was suggested in-vivo and in-vitro ischemic condition and suspected to be related to functional recovery of ischemic cardiomyocyte.

Efficient computational model for identification of antitubercular peptides by integrating amino acid patterns and properties

Tuberculosis (TB) is a leading killer caused by Mycobacterium tuberculosis. Recently, anti-TB peptides have provided an alternative approach to combat antibiotic tolerance. We have developed an effective computational predictor, identification of antitubercular peptides (iAntiTB), by the integration of multiple feature vectors deriving from the amino acid sequences via random forest (RF) and support vector machine (SVM) classifiers. The iAntiTB combines the RF and SVM scores via linear regression to enhance the prediction accuracy. To make a robust and accurate predictor, we prepared the two datasets with different types of negative samples. The iAntiTB achieved area under the ROC curve values of 0.896 and 0.946 on the training datasets of the first and second datasets, respectively. The iAntiTB outperformed the other existing predictors.

Computational identification of microbial phosphorylation sites by the enhanced characteristics of sequence information

Protein phosphorylation on serine (S) and threonine (T) has emerged as a key device in the control of many biological processes. Recently phosphorylation in microbial organisms has attracted much attention for its critical roles in various cellular processes such as cell growth and cell division. Here a novel machine learning predictor, MPSite (Microbial Phosphorylation Site predictor), was developed to identify microbial phosphorylation sites using the enhanced characteristics of sequence features. The final feature vectors optimized via a Wilcoxon rank sum test. A random forest classifier was then trained using the optimum features to build the predictor. Benchmarking investigation using the 5-fold cross-validation and independent datasets test showed that the MPSite is able to achieve robust performance on the S- and T-phosphorylation site prediction. It also outperformed other existing methods on the comprehensive independent datasets. We anticipate that the MPSite is a powerful tool for proteome-wide prediction of microbial phosphorylation sites and facilitates hypothesis-driven functional interrogation of phosphorylation proteins. A web application with the curated datasets is freely available at http://kurata14.bio.kyutech.ac.jp/MPSite/.

Ranking network mechanisms by how they fit diverse experiments and deciding on E. coli's ammonium transport and assimilation network

The complex ammonium transport and assimilation network of E. coli involves the ammonium transporter AmtB, the regulatory proteins GlnK and GlnB, and the central N-assimilating enzymes together with their highly complex interactions. The engineering and modelling of such a complex network seem impossible because functioning depends critically on a gamut of data known at patchy accuracy. We developed a way out of this predicament, which employs: (i) a constrained optimization-based technology for the simultaneous fitting of models to heterogeneous experimental data sets gathered through diverse experimental set-ups, (ii) a 'rubber band method' to deal with different degrees of uncertainty, both in experimentally determined or estimated parameter values and in measured transient or steady-state variables (training data sets), (iii) integration of human expertise to decide on accuracies of both parameters and variables, (iv) massive computation employing a fast algorithm and a supercomputer, (v) an objective way of quantifying the plausibility of models, which makes it possible to decide which model is the best and how much better that model is than the others. We applied the new technology to the ammonium transport and assimilation network, integrating recent and older data of various accuracies, from different expert laboratories. The kinetic model objectively ranked best, has E. coli's AmtB as an active transporter of ammonia to be assimilated with GlnK minimizing the futile cycling that is an inevitable consequence of intracellular ammonium accumulation. It is 130 times better than a model with facilitated passive transport of ammonia.

PreAIP: Computational Prediction of Anti-inflammatory Peptides by Integrating Multiple Complementary Features

Numerous inflammatory diseases and autoimmune disorders by therapeutic peptides have received substantial consideration; however, the exploration of anti-inflammatory peptides via biological experiments is often a time-consuming and expensive task. The development of novel in silico predictors is desired to classify potential anti-inflammatory peptides prior to in vitro investigation. Herein, an accurate predictor, called PreAIP (Predictor of Anti-Inflammatory Peptides) was developed by integrating multiple complementary features. We systematically investigated different types of features including primary sequence, evolutionary and structural information through a random forest classifier. The final PreAIP model achieved an AUC value of 0.833 in the training dataset via 10-fold cross-validation test, which was better than that of existing models. Moreover, we assessed the performance of the PreAIP with an AUC value of 0.840 on a test dataset to demonstrate that the proposed method outperformed the two existing methods. These results indicated that the PreAIP is an accurate predictor for identifying AIPs and contributes to the development of AIPs therapeutics and biomedical research. The curated datasets and the PreAIP are freely available at http://kurata14.bio.kyutech.ac.jp/PreAIP/.

Large-Scale Assessment of Bioinformatics Tools for Lysine Succinylation Sites

Lysine succinylation is a form of posttranslational modification of the proteins that play an essential functional role in every aspect of cell metabolism in both prokaryotes and eukaryotes. Aside from experimental identification of succinylation sites, there has been an intense effort geared towards the development of sequence-based prediction through machine learning, due to its promising and essential properties of being highly accurate, robust and cost-effective. In spite of these advantages, there are several problems that are in need of attention in the design and development of succinylation site predictors. Notwithstanding of many studies on the employment of machine learning approaches, few articles have examined this bioinformatics field in a systematic manner. Thus, we review the advancements regarding the current state-of-the-art prediction models, datasets, and online resources and illustrate the challenges and limitations to present a useful guideline for developing powerful succinylation site prediction tools.

Prediction of S-nitrosylation sites by integrating support vector machines and random forest

Cysteine S-nitrosylation is a type of reversible post-translational modification of proteins, which controls diverse biological processes.

GPSuc: Global Prediction of Generic and Species-specific Succinylation Sites by aggregating multiple sequence features

Lysine succinylation is one of the dominant post-translational modification of the protein that contributes to many biological processes including cell cycle, growth and signal transduction pathways. Identification of succinylation sites is an important step for understanding the function of proteins. The complicated sequence patterns of protein succinylation revealed by proteomic studies highlight the necessity of developing effective species-specific in silico strategies for global prediction succinylation sites. Here we have developed the generic and nine species-specific succinylation site classifiers through aggregating multiple complementary features. We optimized the consecutive features using the Wilcoxon-rank feature selection scheme. The final feature vectors were trained by a random forest (RF) classifier. With an integration of RF scores via logistic regression, the resulting predictor termed GPSuc achieved better performance than other existing generic and species-specific succinylation site predictors. To reveal the mechanism of succinylation and assist hypothesis-driven experimental design, our predictor serves as a valuable resource. To provide a promising performance in large-scale datasets, a web application was developed at http://kurata14.bio.kyutech.ac.jp/GPSuc/.

Long negative feedback loop enhances period tunability of biological oscillators

Oscillatory phenomena play a major role in organisms. In some biological oscillations such as cell cycles and heartbeats, the period can be tuned without significant changes in the amplitude. This property is called (period) tunability, one of the prominent features of biological oscillations. However, how biological oscillators produce tunable oscillations remains largely unexplored. We tackle this question using computational experiments. It has been reported that positive-plus-negative feedback oscillators produce tunable oscillations through the hysteresis-based mechanism. First, in this study, we confirmed that positive-plus-negative feedback oscillators generate tunable oscillations. Second, we found that tunability is positively correlated with the dynamic range of oscillations. Third, we showed that long negative feedback oscillators without any additional positive feedback loops can produce tunable oscillations. Finally, we computationally demonstrated that by lengthening the negative feedback loop, the Repressilator, known as a non-tunable synthetic gene oscillator, can be converted into a tunable oscillator. This work provides synthetic biologists with clues to design tunable gene oscillators.

A novel isolation method for cancer prognostic factors via the p53 pathway by a combination of in vitro and in silico analyses

Identifying new therapeutic target genes affecting the survival of patients with cancer is crucial for the development of new cancer therapies. Here, we developed a novel technology combining in vitro short hairpin RNA (shRNA) library screening and in silico analysis of the tumor transcriptome to identify prognostic factors via the p53 tumor-suppressor pathway. For initial screening, we screened 5,000 genes through selection of shRNAs in p53 wild-type tumor cells that altered sensitivity to the p53 activator actinomycin D (ActD) to identify p53 regulatory genes; shRNAs targeting 322 genes were obtained. Among these 322 genes, seven were prognostic factor candidates whose high expression increased ActD sensitivity while prolonging the survival period in patients with the p53 wild-type genotype. Conversely, we identified 33 genes as prognostic factor candidates among ActD-resistant genes related to a shortened survival period only in p53 wild-type tumors. These 40 genes had biological functions such as apoptosis, drug response, cell cycle checkpoint, and cell proliferation. The 40 genes selected by this method contained many known genes related to the p53 pathway and prognosis in patients with cancer. In summary, we developed an efficient screening method to identify p53-dependent prognostic factors with in vitro experimental data and database analysis.

Modeling and simulation of the redox regulation of the metabolism in Escherichia coli at different oxygen concentrations

BACKGROUND

Microbial production of biofuels and biochemicals from renewable feedstocks has received considerable recent attention from environmental protection and energy production perspectives. Many biofuels and biochemicals are produced by fermentation under oxygen-limited conditions following initiation of aerobic cultivation to enhance the cell growth rate. Thus, it is of significant interest to investigate the effect of dissolved oxygen concentration on redox regulation in Escherichia coli, a particularly popular cellular factory due to its high growth rate and well-characterized physiology. For this, the systems biology approach such as modeling is powerful for the analysis of the metabolism and for the design of microbial cellular factories.

RESULTS

Here, we developed a kinetic model that describes the dynamics of fermentation by taking into account transcription factors such as ArcA/B and Fnr, respiratory chain reactions and fermentative pathways, and catabolite regulation. The hallmark of the kinetic model is its ability to predict the dynamics of metabolism at different dissolved oxygen levels and facilitate the rational design of cultivation methods. The kinetic model was verified based on the experimental data for a wild-type E. coli strain. The model reasonably predicted the metabolic characteristics and molecular mechanisms of fnr and arcA gene-knockout mutants. Moreover, an aerobic-microaerobic dual-phase cultivation method for lactate production in a pfl-knockout mutant exhibited promising yield and productivity.

CONCLUSIONS

It is quite important to understand metabolic regulation mechanisms from both scientific and engineering points of view. In particular, redox regulation in response to oxygen limitation is critically important in the practical production of biofuel and biochemical compounds. The developed model can thus be used as a platform for designing microbial factories to produce a variety of biofuels and biochemicals.