milRNApredictor: Genome-free prediction of fungi milRNAs by incorporating k-mer scheme and distance-dependent pair potential.

MicroRNA-like small RNAs (milRNAs) with length of 21-22 nucleotides are a type of small non-coding RNAs that are firstly found in Neurospora crassa in 2010. Identifying milRNAs of species without genomic information is a difficult problem. Here, knowledge-based energy features are developed to identify milRNAs by tactfully incorporating k-mer scheme and distance-dependent pair potential. Compared with k-mer scheme, features developed here can alleviate the inherent curse of dimensionality in k-scheme once k becomes large. In addition, milRNApredictor built on novel features performs comparably to k-mer scheme, and achieves sensitivity of 74.21%, and specificity of 75.72% based on 10-fold cross-validation. Furthermore, for novel miRNA prediction, there exists high overlap of results from milRNApredictor and state-of-the-art mirnovo. However, milRNApredictor is simpler to use with reduced requirements of input data and dependencies. Taken together, milRNApredictor can be used to de novo identify fungi milRNAs and other very short small RNAs of non-model organisms.

ANDIS: an atomic angle- and distance-dependent statistical potential for protein structure quality assessment.

BACKGROUND: The knowledge-based statistical potential has been widely used in protein structure modeling and model quality assessment. They are commonly evaluated based on their abilities of native recognition as well as decoy discrimination. However, these two aspects are found to be mutually exclusive in many statistical potentials. RESULTS: We developed an atomic ANgle- and DIStance-dependent (ANDIS) statistical potential for protein structure quality assessment with distance cutoff being a tunable parameter. When distance cutoff is ≤9.0 Å, "effective atomic interaction" is employed to enhance the ability of native recognition. For a distance cutoff of ≥10 Å, the distance-dependent atom-pair potential with random-walk reference state is combined to strengthen the ability of decoy discrimination. Benchmark tests on 632 structural decoy sets from diverse sources demonstrate that ANDIS outperforms other state-of-the-art potentials in both native recognition and decoy discrimination. CONCLUSIONS: Distance cutoff is a crucial parameter for distance-dependent statistical potentials. A lower distance cutoff is better for native recognition, while a higher one is favorable for decoy discrimination. The ANDIS potential is freely available as a standalone application at http://qbp.hzau.edu.cn/ANDIS/ .

Protein structure prediction.

Predicting 3D structure of protein from its amino acid sequence is one of the most important unsolved problems in biophysics and computational biology. This paper attempts to give a comprehensive introduction of the most recent effort and progress on protein structure prediction. Following the general flowchart of structure prediction, related concepts and methods are presented and discussed. Moreover, brief introductions are made to several widely-used prediction methods and the community-wide critical assessment of protein structure prediction (CASP) experiments.

Diverse effects of distance cutoff and residue interval on the performance of distance-dependent atom-pair potential in protein structure prediction.

BACKGROUND: As one of the most successful knowledge-based energy functions, the distance-dependent atom-pair potential is widely used in all aspects of protein structure prediction, including conformational search, model refinement, and model assessment. During the last two decades, great efforts have been made to improve the reference state of the potential, while other factors that also strongly affect the performance of the potential have been relatively less investigated. RESULTS: Based on different distance cutoffs (from 5 to 22 Å) and residue intervals (from 0 to 15) as well as six different reference states, we constructed a series of distance-dependent atom-pair potentials and tested them on several groups of structural decoy sets collected from diverse sources. A comprehensive investigation has been performed to clarify the effects of distance cutoff and residue interval on the potential's performance. Our results provide a new perspective as well as a practical guidance for optimizing distance-dependent statistical potentials. CONCLUSIONS: The optimal distance cutoff and residue interval are highly related with the reference state that the potential is based on, the measurements of the potential's performance, and the decoy sets that the potential is applied to. The performance of distance-dependent statistical potential can be significantly improved when the best statistical parameters for the specific application environment are adopted.

3DRobot: automated generation of diverse and well-packed protein structure decoys.

MOTIVATION: Computationally generated non-native protein structure conformations (or decoys) are often used for designing protein folding simulation methods and force fields. However, almost all the decoy sets currently used in literature suffer from uneven root mean square deviation (RMSD) distribution with bias to non-protein like hydrogen-bonding and compactness patterns. Meanwhile, most protein decoy sets are pre-calculated and there is a lack of methods for automated generation of high-quality decoys for any target proteins. RESULTS: We developed a new algorithm, 3DRobot, to create protein structure decoys by free fragment assembly with enhanced hydrogen-bonding and compactness interactions. The method was benchmarked with three widely used decoy sets from ab initio folding and comparative modeling simulations. The decoys generated by 3DRobot are shown to have significantly enhanced diversity and evenness with a continuous distribution in the RMSD space. The new energy terms introduced in 3DRobot improve the hydrogen-bonding network and compactness of decoys, which eliminates the possibility of native structure recognition by trivial potentials. Algorithms that can automatically create such diverse and well-packed non-native conformations from any protein structure should have a broad impact on the development of advanced protein force field and folding simulation methods. AVAILIABLITY AND IMPLEMENTATION: http://zhanglab.ccmb.med.umich.edu/3DRobot/ CONTACT: jiay@phy.ccnu.edu.cn; zhng@umich.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

What is the best reference state for designing statistical atomic potentials in protein structure prediction?

Many statistical potentials were developed in last two decades for protein folding and protein structure recognition. The major difference of these potentials is on the selection of reference states to offset sampling bias. However, since these potentials used different databases and parameter cutoffs, it is difficult to judge what the best reference states are by examining the original programs. In this study, we aim to address this issue and evaluate the reference states by a unified database and programming environment. We constructed distance-specific atomic potentials using six widely-used reference states based on 1022 high-resolution protein structures, which are applied to rank modeling in six sets of structure decoys. The reference state on random-walk chain outperforms others in three decoy sets while those using ideal-gas, quasi-chemical approximation and averaging sample stand out in one set separately. Nevertheless, the performance of the potentials relies on the origin of decoy generations and no reference state can clearly outperform others in all decoy sets. Further analysis reveals that the statistical potentials have a contradiction between the universality and pertinence, and optimal reference states should be extracted based on specific application environments and decoy spaces.

Impact of bounded noise on the formation and instability of spiral wave in a 2D Lattice of neurons.

Spiral waves in the neocortex may provide a spatial framework to organize cortical oscillations, thus help signal communication. However, noise influences spiral wave. Many previous theoretical studies about noise mainly focus on unbounded Gaussian noise, which contradicts that a real physical quantity is always bounded. Furthermore, non-Gaussian noise is also important for dynamical behaviors of excitable media. Nevertheless, there are no results concerning the effect of bounded noise on spiral wave till now. Based on Hodgkin-Huxley neuron model subjected to bounded noise with the form of Asin[ωt + σW(t)], the influences of bounded noise on the formation and instability of spiral wave in a two-dimensional (2D) square lattice of neurons are investigated in detail by separately adjusting the intensity σ, amplitude A, and frequency f of bounded noise. It is found that the increased intensity σ can facilitate the formation of spiral wave while the increased amplitude A tends to destroy spiral wave. Furthermore, frequency of bounded noise has the effect of facilitation or inhibition on pattern synchronization. Interestingly, for the appropriate intensity, amplitude and frequency can separately induce resonance-like phenomenon.

Impact of Bounded Noise and Rewiring on the Formation and Instability of Spiral Waves in a Small-World Network of Hodgkin-Huxley Neurons.

Spiral waves are observed in the chemical, physical and biological systems, and the emergence of spiral waves in cardiac tissue is linked to some diseases such as heart ventricular fibrillation and epilepsy; thus it has importance in theoretical studies and potential medical applications. Noise is inevitable in neuronal systems and can change the electrical activities of neuron in different ways. Many previous theoretical studies about the impacts of noise on spiral waves focus an unbounded Gaussian noise and even colored noise. In this paper, the impacts of bounded noise and rewiring of network on the formation and instability of spiral waves are discussed in small-world (SW) network of Hodgkin-Huxley (HH) neurons through numerical simulations, and possible statistical analysis will be carried out. Firstly, we present SW network of HH neurons subjected to bounded noise. Then, it is numerically demonstrated that bounded noise with proper intensity σ, amplitude A, or frequency f can facilitate the formation of spiral waves when rewiring probability p is below certain thresholds. In other words, bounded noise-induced resonant behavior can occur in the SW network of neurons. In addition, rewiring probability p always impairs spiral waves, while spiral waves are confirmed to be robust for small p, thus shortcut-induced phase transition of spiral wave with the increase of p is induced. Furthermore, statistical factors of synchronization are calculated to discern the phase transition of spatial pattern, and it is confirmed that larger factor of synchronization is approached with increasing of rewiring probability p, and the stability of spiral wave is destroyed.

Dynamics and robustness of the cardiac progenitor cell induced pluripotent stem cell network during cell phenotypes transition.

Robustness is a fundamental characteristic of biological systems since all living systems need to adapt to internal or external perturbations, unpredictable environments, stochastic events and unreliable components, and so on. A long-term challenge in systems biology is to reveal the origin of robustness underlying molecular regulator network. In this study, a simple Boolean model is used to investigate the global dynamic properties and robustness of cardiac progenitor cell (CPC) induced pluripotent stem cell network that governs reprogramming and directed differentiation process. It is demonstrated that two major attractors correspond to source and target cell phenotypes, respectively, and two dominating attracting trajectories characterise the biological pathways between two major cell phenotypes. In particular, the experimentally observed transition between different cell phenotypes can be reproduced and explained theoretically. Furthermore, the robustness of major attractors and trajectories is largely maintained with respect to small perturbations to the network. Taken together, the CPC-induced pluripotent stem cell network is extremely robustly designed for their functions.

plantMirP: an efficient computational program for the prediction of plant pre-miRNA by incorporating knowledge-based energy features.

MicroRNAs are a predominant type of small non-coding RNAs approximately 21 nucleotides in length that play an essential role at the post-transcriptional level by either RNA degradation, translational repression or both through an RNA-induced silencing complex. Identification of these molecules can aid the dissecting of their regulatory functions. The secondary structures of plant pre-miRNAs are much more complex than those of animal pre-miRNAs. In contrast to prediction tools for animal pre-miRNAs, much less effort has been contributed to plant pre-miRNAs. In this study, a set of novel knowledge-based energy features that has very high discriminatory power is proposed and incorporated with the existing features for specifically distinguishing the hairpins of real/pseudo plant pre-miRNAs. A promising performance area under a receiver operating characteristic curve of 0.9444 indicates that 5 knowledge-based energy features have very high discriminatory power. The 10-fold cross-validation result demonstrates that plantMirP with full features has a promising sensitivity of 92.61% and a specificity of 98.88%. Based on various different datasets, it was found that plantMirP has a higher prediction performance by comparison with miPlantPreMat, PlantMiRNAPred, triplet-SVM, and microPred. Meanwhile, plantMirP can greatly balance sensitivity and specificity for real/pseudo plant pre-miRNAs. Taken together, we developed a promising SVM-based program, plantMirP, for predicting plant pre-miRNAs by incorporating knowledge-based energy features. This study shows it to be a valuable tool for miRNA-related studies.

qiRNApredictor: A Novel Computational Program for the Prediction of qiRNAs in Neurospora crassa.

Recently, a new type of small interfering RNAs (qiRNAs) of typically 20~21 nucleotides was found in Neurospora crassa and rice and has been shown to regulate gene silencing in the DNA damage response. Identification of qiRNAs is fundamental for dissecting regulatory functions and molecular mechanisms. In contrast to other expensive and time-consuming experimental methods, the computational prediction of qiRNAs is a conveniently rapid method for gaining valuable information for a subsequent experimental verification. However, no tool existed to date for the prediction of qiRNAs. To this purpose, we developed the novel qiRNA prediction software package qiRNApredictor. This software demonstrates a promising sensitivity of 93.55% and a specificity of 71.61% from the leave-one-out validation. These studies might be beneficial for further experimental investigation. Furthermore, the local package of qiRNApredictor was implemented and made freely available to the academic community at Supplementary material.

Noise Decomposition Principle in a Coherent Feed-Forward Transcriptional Regulatory Loop.

Coherent feed-forward loops exist extensively in realistic biological regulatory systems, and are common signaling motifs. Here, we study the characteristics and the propagation mechanism of the output noise in a coherent feed-forward transcriptional regulatory loop that can be divided into a main road and branch. Using the linear noise approximation, we derive analytical formulae for the total noise of the full loop, the noise of the branch, and the noise of the main road, which are verified by the Gillespie algorithm. Importantly, we find that (i) compared with the branch motif or the main road motif, the full motif can effectively attenuate the output noise level; (ii) there is a transition point of system state such that the noise of the main road is dominated when the underlying system is below this point, whereas the noise of the branch is dominated when the system is beyond the point. The entire analysis reveals the mechanism of how the noise is generated and propagated in a simple yet representative signaling module.