-
Autoregressive Enzyme Function Prediction with Multi-scale Multi-modality Fusion
Authors:
Dingyi Rong,
Wenzhuo Zheng,
Bozitao Zhong,
Zhouhan Lin,
Liang Hong,
Ning Liu
Abstract:
Accurate prediction of enzyme function is crucial for elucidating biological mechanisms and driving innovation across various sectors. Existing deep learning methods tend to rely solely on either sequence data or structural data and predict the EC number as a whole, neglecting the intrinsic hierarchical structure of EC numbers. To address these limitations, we introduce MAPred, a novel multi-modal…
▽ More
Accurate prediction of enzyme function is crucial for elucidating biological mechanisms and driving innovation across various sectors. Existing deep learning methods tend to rely solely on either sequence data or structural data and predict the EC number as a whole, neglecting the intrinsic hierarchical structure of EC numbers. To address these limitations, we introduce MAPred, a novel multi-modality and multi-scale model designed to autoregressively predict the EC number of proteins. MAPred integrates both the primary amino acid sequence and the 3D tokens of proteins, employing a dual-pathway approach to capture comprehensive protein characteristics and essential local functional sites. Additionally, MAPred utilizes an autoregressive prediction network to sequentially predict the digits of the EC number, leveraging the hierarchical organization of EC classifications. Evaluations on benchmark datasets, including New-392, Price, and New-815, demonstrate that our method outperforms existing models, marking a significant advance in the reliability and granularity of protein function prediction within bioinformatics.
△ Less
Submitted 11 August, 2024;
originally announced August 2024.
-
Geneverse: A collection of Open-source Multimodal Large Language Models for Genomic and Proteomic Research
Authors:
Tianyu Liu,
Yijia Xiao,
Xiao Luo,
Hua Xu,
W. Jim Zheng,
Hongyu Zhao
Abstract:
The applications of large language models (LLMs) are promising for biomedical and healthcare research. Despite the availability of open-source LLMs trained using a wide range of biomedical data, current research on the applications of LLMs to genomics and proteomics is still limited. To fill this gap, we propose a collection of finetuned LLMs and multimodal LLMs (MLLMs), known as Geneverse, for th…
▽ More
The applications of large language models (LLMs) are promising for biomedical and healthcare research. Despite the availability of open-source LLMs trained using a wide range of biomedical data, current research on the applications of LLMs to genomics and proteomics is still limited. To fill this gap, we propose a collection of finetuned LLMs and multimodal LLMs (MLLMs), known as Geneverse, for three novel tasks in genomic and proteomic research. The models in Geneverse are trained and evaluated based on domain-specific datasets, and we use advanced parameter-efficient finetuning techniques to achieve the model adaptation for tasks including the generation of descriptions for gene functions, protein function inference from its structure, and marker gene selection from spatial transcriptomic data. We demonstrate that adapted LLMs and MLLMs perform well for these tasks and may outperform closed-source large-scale models based on our evaluations focusing on both truthfulness and structural correctness. All of the training strategies and base models we used are freely accessible.
△ Less
Submitted 21 June, 2024;
originally announced June 2024.
-
Seeing through the Brain: Image Reconstruction of Visual Perception from Human Brain Signals
Authors:
Yu-Ting Lan,
Kan Ren,
Yansen Wang,
Wei-Long Zheng,
Dongsheng Li,
Bao-Liang Lu,
Lili Qiu
Abstract:
Seeing is believing, however, the underlying mechanism of how human visual perceptions are intertwined with our cognitions is still a mystery. Thanks to the recent advances in both neuroscience and artificial intelligence, we have been able to record the visually evoked brain activities and mimic the visual perception ability through computational approaches. In this paper, we pay attention to vis…
▽ More
Seeing is believing, however, the underlying mechanism of how human visual perceptions are intertwined with our cognitions is still a mystery. Thanks to the recent advances in both neuroscience and artificial intelligence, we have been able to record the visually evoked brain activities and mimic the visual perception ability through computational approaches. In this paper, we pay attention to visual stimuli reconstruction by reconstructing the observed images based on portably accessible brain signals, i.e., electroencephalography (EEG) data. Since EEG signals are dynamic in the time-series format and are notorious to be noisy, processing and extracting useful information requires more dedicated efforts; In this paper, we propose a comprehensive pipeline, named NeuroImagen, for reconstructing visual stimuli images from EEG signals. Specifically, we incorporate a novel multi-level perceptual information decoding to draw multi-grained outputs from the given EEG data. A latent diffusion model will then leverage the extracted information to reconstruct the high-resolution visual stimuli images. The experimental results have illustrated the effectiveness of image reconstruction and superior quantitative performance of our proposed method.
△ Less
Submitted 16 August, 2023; v1 submitted 27 July, 2023;
originally announced August 2023.
-
A Deep Learning Approach to Predicting Ventilator Parameters for Mechanically Ventilated Septic Patients
Authors:
Zhijun Zeng,
Zhen Hou,
Ting Li,
Lei Deng,
Jianguo Hou,
Xinran Huang,
Jun Li,
Meirou Sun,
Yunhan Wang,
Qiyu Wu,
Wenhao Zheng,
Hua Jiang,
Qi Wang
Abstract:
We develop a deep learning approach to predicting a set of ventilator parameters for a mechanically ventilated septic patient using a long and short term memory (LSTM) recurrent neural network (RNN) model. We focus on short-term predictions of a set of ventilator parameters for the septic patient in emergency intensive care unit (EICU). The short-term predictability of the model provides attending…
▽ More
We develop a deep learning approach to predicting a set of ventilator parameters for a mechanically ventilated septic patient using a long and short term memory (LSTM) recurrent neural network (RNN) model. We focus on short-term predictions of a set of ventilator parameters for the septic patient in emergency intensive care unit (EICU). The short-term predictability of the model provides attending physicians with early warnings to make timely adjustment to the treatment of the patient in the EICU. The patient specific deep learning model can be trained on any given critically ill patient, making it an intelligent aide for physicians to use in emergent medical situations.
△ Less
Submitted 20 February, 2022;
originally announced February 2022.
-
Mining of high throughput screening database reveals AP-1 and autophagy pathways as potential targets for COVID-19 therapeutics
Authors:
Hu Zhu,
Catherine Z. Chen,
Srilatha Sakamuru,
Anton Simeonov,
Mathew D. Hall,
Menghang Xia,
Wei Zheng,
Ruili Huang
Abstract:
The recent global pandemic of Coronavirus Disease 2019 (COVID-19) caused by the new coronavirus SARS-CoV-2 presents an urgent need for new therapeutic candidates. Many efforts have been devoted to screening existing drug libraries with the hope to repurpose approved drugs as potential treatments for COVID-19. However, the antiviral mechanisms of action for the drugs found active in these phenotypi…
▽ More
The recent global pandemic of Coronavirus Disease 2019 (COVID-19) caused by the new coronavirus SARS-CoV-2 presents an urgent need for new therapeutic candidates. Many efforts have been devoted to screening existing drug libraries with the hope to repurpose approved drugs as potential treatments for COVID-19. However, the antiviral mechanisms of action for the drugs found active in these phenotypic screens are largely unknown. To deconvolute the viral targets for more effective anti-COVID-19 drug development, we mined our in-house database of approved drug screens against 994 assays and compared their activity profiles with the drug activity profile in a cytopathic effect (CPE) assay of SARS-CoV-2. We found that the autophagy and AP-1 signaling pathway activity profiles are significantly correlated with the anti-SARS-CoV-2 activity profile. In addition, a class of neurology/psychiatry drugs was found significantly enriched with anti-SARS-CoV-2 activity. Taken together, these results have provided new insights into SARS-CoV-2 infection and potential targets for COVID-19 therapeutics.
△ Less
Submitted 23 July, 2020;
originally announced July 2020.
-
Total Variation Regularization for Compartmental Epidemic Models with Time-Varying Dynamics
Authors:
Wenjie Zheng
Abstract:
Compartmental epidemic models are among the most popular ones in epidemiology. For the parameters (e.g., the transmission rate) characterizing these models, the majority of researchers simplify them as constants, while some others manage to detect their continuous variations. In this paper, we aim at capturing, on the other hand, discontinuous variations, which better describe the impact of many n…
▽ More
Compartmental epidemic models are among the most popular ones in epidemiology. For the parameters (e.g., the transmission rate) characterizing these models, the majority of researchers simplify them as constants, while some others manage to detect their continuous variations. In this paper, we aim at capturing, on the other hand, discontinuous variations, which better describe the impact of many noteworthy events, such as city lockdowns, the opening of field hospitals, and the mutation of the virus, whose effect should be instant. To achieve this, we balance the model's likelihood by total variation, which regulates the temporal variations of the model parameters. To infer these parameters, instead of using Monte Carlo methods, we design a novel yet straightforward optimization algorithm, dubbed Iterated Nelder--Mead, which repeatedly applies the Nelder--Mead algorithm. Experiments conducted on the simulated data demonstrate that our approach can reproduce these discontinuities and precisely depict the epidemics.
△ Less
Submitted 29 May, 2020; v1 submitted 1 April, 2020;
originally announced April 2020.
-
Protein structure and sequence re-analysis of 2019-nCoV genome does not indicate snakes as its intermediate host or the unique similarity between its spike protein insertions and HIV-1
Authors:
Chengxin Zhang,
Wei Zheng,
Xiaoqiang Huang,
Eric W. Bell,
Xiaogen Zhou,
Yang Zhang
Abstract:
As the infection of 2019-nCoV coronavirus is quickly developing into a global pneumonia epidemic, careful analysis of its transmission and cellular mechanisms is sorely needed. In this report, we re-analyzed the computational approaches and findings presented in two recent manuscripts by Ji et al. (https://doi.org/10.1002/jmv.25682) and by Pradhan et al. (https://doi.org/10.1101/2020.01.30.927871)…
▽ More
As the infection of 2019-nCoV coronavirus is quickly developing into a global pneumonia epidemic, careful analysis of its transmission and cellular mechanisms is sorely needed. In this report, we re-analyzed the computational approaches and findings presented in two recent manuscripts by Ji et al. (https://doi.org/10.1002/jmv.25682) and by Pradhan et al. (https://doi.org/10.1101/2020.01.30.927871), which concluded that snakes are the intermediate hosts of 2019-nCoV and that the 2019-nCoV spike protein insertions shared a unique similarity to HIV-1. Results from our re-implementation of the analyses, built on larger-scale datasets using state-of-the-art bioinformatics methods and databases, do not support the conclusions proposed by these manuscripts. Based on our analyses and existing data of coronaviruses, we concluded that the intermediate hosts of 2019-nCoV are more likely to be mammals and birds than snakes, and that the "novel insertions" observed in the spike protein are naturally evolved from bat coronaviruses.
△ Less
Submitted 8 February, 2020;
originally announced February 2020.
-
$hv$-Block Cross Validation is not a BIBD: a Note on the Paper by Jeff Racine (2000)
Authors:
Wenjie Zheng
Abstract:
This note corrects a mistake in the paper "consistent cross-validatory model-selection for dependent data: $hv$-block cross-validation" by Racine (2000). In his paper, he implied that the therein proposed $hv$-block cross-validation is consistent in the sense of Shao (1993). To get this intuition, he relied on the speculation that $hv$-block is a balanced incomplete block design (BIBD). This note…
▽ More
This note corrects a mistake in the paper "consistent cross-validatory model-selection for dependent data: $hv$-block cross-validation" by Racine (2000). In his paper, he implied that the therein proposed $hv$-block cross-validation is consistent in the sense of Shao (1993). To get this intuition, he relied on the speculation that $hv$-block is a balanced incomplete block design (BIBD). This note demonstrates that this is not the case, and thus the theoretical consistency of $hv$-block remains an open question. In addition, I also provide a Python program counting the number of occurrences of each sample and each pair of samples.
△ Less
Submitted 20 October, 2019;
originally announced October 2019.
-
A Random Interaction Forest for Prioritizing Predictive Biomarkers
Authors:
Zhen Zeng,
Yuefeng Lu,
Judong Shen,
Wei Zheng,
Peter Shaw,
Mary Beth Dorr
Abstract:
Precision medicine is becoming a focus in medical research recently, as its implementation brings values to all stakeholders in the healthcare system. Various statistical methodologies have been developed tackling problems in different aspects of this field, e.g., assessing treatment heterogeneity, identifying patient subgroups, or building treatment decision models. However, there is a lack of ne…
▽ More
Precision medicine is becoming a focus in medical research recently, as its implementation brings values to all stakeholders in the healthcare system. Various statistical methodologies have been developed tackling problems in different aspects of this field, e.g., assessing treatment heterogeneity, identifying patient subgroups, or building treatment decision models. However, there is a lack of new tools devoted to selecting and prioritizing predictive biomarkers. We propose a novel tree-based ensemble method, random interaction forest (RIF), to generate predictive importance scores and prioritize candidate biomarkers for constructing refined treatment decision models. RIF was evaluated by comparing with the conventional random forest and univariable regression methods and showed favorable properties under various simulation scenarios. We applied the proposed RIF method to a biomarker dataset from two phase III clinical trials of bezlotoxumab on $\textit{Clostridium difficile}$ infection recurrence and obtained biologically meaningful results.
△ Less
Submitted 3 October, 2019;
originally announced October 2019.
-
A Novel Bi-hemispheric Discrepancy Model for EEG Emotion Recognition
Authors:
Yang Li,
Wenming Zheng,
Lei Wang,
Yuan Zong,
Lei Qi,
Zhen Cui,
Tong Zhang,
Tengfei Song
Abstract:
The neuroscience study has revealed the discrepancy of emotion expression between left and right hemispheres of human brain. Inspired by this study, in this paper, we propose a novel bi-hemispheric discrepancy model (BiHDM) to learn the asymmetric differences between two hemispheres for electroencephalograph (EEG) emotion recognition. Concretely, we first employ four directed recurrent neural netw…
▽ More
The neuroscience study has revealed the discrepancy of emotion expression between left and right hemispheres of human brain. Inspired by this study, in this paper, we propose a novel bi-hemispheric discrepancy model (BiHDM) to learn the asymmetric differences between two hemispheres for electroencephalograph (EEG) emotion recognition. Concretely, we first employ four directed recurrent neural networks (RNNs) based on two spatial orientations to traverse electrode signals on two separate brain regions, which enables the model to obtain the deep representations of all the EEG electrodes' signals while keeping the intrinsic spatial dependence. Then we design a pairwise subnetwork to capture the discrepancy information between two hemispheres and extract higher-level features for final classification. Besides, in order to reduce the domain shift between training and testing data, we use a domain discriminator that adversarially induces the overall feature learning module to generate emotion-related but domain-invariant feature, which can further promote EEG emotion recognition. We conduct experiments on three public EEG emotional datasets, and the experiments show that the new state-of-the-art results can be achieved.
△ Less
Submitted 10 May, 2019;
originally announced June 2019.
-
Predicting protein inter-residue contacts using composite likelihood maximization and deep learning
Authors:
Haicang Zhang,
Qi Zhang,
Fusong Ju,
Jianwei Zhu,
Shiwei Sun,
Yujuan Gao,
Ziwei Xie,
Minghua Deng,
Shiwei Sun,
Wei-Mou Zheng,
Dongbo Bu
Abstract:
Accurate prediction of inter-residue contacts of a protein is important to calcu- lating its tertiary structure. Analysis of co-evolutionary events among residues has been proved effective to inferring inter-residue contacts. The Markov ran- dom field (MRF) technique, although being widely used for contact prediction, suffers from the following dilemma: the actual likelihood function of MRF is acc…
▽ More
Accurate prediction of inter-residue contacts of a protein is important to calcu- lating its tertiary structure. Analysis of co-evolutionary events among residues has been proved effective to inferring inter-residue contacts. The Markov ran- dom field (MRF) technique, although being widely used for contact prediction, suffers from the following dilemma: the actual likelihood function of MRF is accurate but time-consuming to calculate, in contrast, approximations to the actual likelihood, say pseudo-likelihood, are efficient to calculate but inaccu- rate. Thus, how to achieve both accuracy and efficiency simultaneously remains a challenge. In this study, we present such an approach (called clmDCA) for contact prediction. Unlike plmDCA using pseudo-likelihood, i.e., the product of conditional probability of individual residues, our approach uses composite- likelihood, i.e., the product of conditional probability of all residue pairs. Com- posite likelihood has been theoretically proved as a better approximation to the actual likelihood function than pseudo-likelihood. Meanwhile, composite likelihood is still efficient to maximize, thus ensuring the efficiency of clmDCA. We present comprehensive experiments on popular benchmark datasets, includ- ing PSICOV dataset and CASP-11 dataset, to show that: i) clmDCA alone outperforms the existing MRF-based approaches in prediction accuracy. ii) When equipped with deep learning technique for refinement, the prediction ac- curacy of clmDCA was further significantly improved, suggesting the suitability of clmDCA for subsequent refinement procedure. We further present successful application of the predicted contacts to accurately build tertiary structures for proteins in the PSICOV dataset.
Accessibility: The software clmDCA and a server are publicly accessible through http://protein.ict.ac.cn/clmDCA/.
△ Less
Submitted 31 August, 2018;
originally announced September 2018.
-
Improving protein threading accuracy via combining local and global potential using TreeCRF model
Authors:
Haicang Zhang,
Mingfu Shao,
Chao Wang,
Jianwei Zhu,
Wei-Mou Zheng,
Dongbo Bu
Abstract:
Protein structure prediction remains to be an open problem in bioinformatics. There are two main categories of methods for protein structure prediction: Free Modeling (FM) and Template Based Modeling (TBM). Protein threading, belonging to the category of template based modeling, identifies the most likely fold with the target by making a sequence-structure alignment between target protein and temp…
▽ More
Protein structure prediction remains to be an open problem in bioinformatics. There are two main categories of methods for protein structure prediction: Free Modeling (FM) and Template Based Modeling (TBM). Protein threading, belonging to the category of template based modeling, identifies the most likely fold with the target by making a sequence-structure alignment between target protein and template protein. Though protein threading has been shown to more be successful for protein structure prediction, it performs poorly for remote homology detection.
△ Less
Submitted 11 September, 2015;
originally announced September 2015.
-
TOPO: Improving remote homologue recognition via identifying common protein structure framework
Authors:
Jianwei Zhu,
Haicang Zhang,
Chao Wang,
Bin Ling,
Wei-Mou Zheng,
Dongbo Bu
Abstract:
Protein structure prediction remains a challenge in the field of computational biology. Traditional protein structure prediction approaches include template-based modelling (say, homology modelling, and threading), and ab initio. A threading algorithm takes a query protein sequence as input, recognizes the most likely fold, and finally reports the alignments of the query sequence to structure-know…
▽ More
Protein structure prediction remains a challenge in the field of computational biology. Traditional protein structure prediction approaches include template-based modelling (say, homology modelling, and threading), and ab initio. A threading algorithm takes a query protein sequence as input, recognizes the most likely fold, and finally reports the alignments of the query sequence to structure-known templates as output. The existing threading approaches mainly utilizes the information of protein sequence profile, solvent accessibility, contact probability, etc., and correctly recognize folds for some proteins. However, the existing threading approaches show poorly performance for remote homology proteins. How to improve the fold recognition for remote homology proteins remains to be a difficult task for protein structure prediction.
△ Less
Submitted 12 July, 2015;
originally announced July 2015.
-
Learning to Fold Proteins Using Energy Landscape Theory
Authors:
N. P. Schafer,
B. L. Kim,
W. Zheng,
P. G. Wolynes
Abstract:
This review is a tutorial for scientists interested in the problem of protein structure prediction, particularly those interested in using coarse-grained molecular dynamics models that are optimized using lessons learned from the energy landscape theory of protein folding. We also present a review of the results of the AMH/AMC/AMW/AWSEM family of coarse-grained molecular dynamics protein folding m…
▽ More
This review is a tutorial for scientists interested in the problem of protein structure prediction, particularly those interested in using coarse-grained molecular dynamics models that are optimized using lessons learned from the energy landscape theory of protein folding. We also present a review of the results of the AMH/AMC/AMW/AWSEM family of coarse-grained molecular dynamics protein folding models to illustrate the points covered in the first part of the article. Accurate coarse-grained structure prediction models can be used to investigate a wide range of conceptual and mechanistic issues outside of protein structure prediction; specifically, the paper concludes by reviewing how AWSEM has in recent years been able to elucidate questions related to the unusual kinetic behavior of artificially designed proteins, multidomain protein misfolding, and the initial stages of protein aggregation.
△ Less
Submitted 24 December, 2013;
originally announced December 2013.
-
Optimizing weights of protein energy function to improve ab initio protein structure prediction
Authors:
Chao Wang,
Yi Wei,
Juntao Liu,
Haicang Zhang,
Bin Ling,
Shuai Cheng Li,
Wei-Mou Zheng,
Dongbo Bu
Abstract:
Predicting protein 3D structure from amino acid sequence remains as a challenge in the field of computational biology. If protein structure homologues are not found, one has to construct 3D structural conformations from the very beginning by the so-called ab initio approach, using some empirical energy functions. A successful algorithm in this category, Rosetta, creates an ensemble of decoy confor…
▽ More
Predicting protein 3D structure from amino acid sequence remains as a challenge in the field of computational biology. If protein structure homologues are not found, one has to construct 3D structural conformations from the very beginning by the so-called ab initio approach, using some empirical energy functions. A successful algorithm in this category, Rosetta, creates an ensemble of decoy conformations by assembling selected best short fragments of known protein structures and then recognizes the native state as the highly populated one with a very low energy. Typically, an energy function is a combination of a variety of terms characterizing different structural features, say hydrophobic interactions, van der Waals force, hydrogen bonding, etc. It is critical for an energy function to be capable to distinguish native-like conformations from non-native ones and to drive most initial conformations assembled from fragments to a native-like one in a conformation search process. In this paper we propose a linear programming algorithm to optimize weighting of a total of 14 energy terms used in Rosetta. We reverse the Monte Carlo process of Rosetta to approach native-like conformations to a process generating from the native state an ensemble of initial conformations most relevant to the native state. Intuitively, an ideal weighting scheme would result in a large "basin of attraction" of the native structure, which leads to an objective function for the linear programming. We have examined the proposal on several benchmark proteins, and the experimental results suggest that the optimized weights enlarge the attraction basin of the native state and improve the quality of the predicted native states as well. In addition, a comparison of optimal weighting schema for proteins of different classes indicates that in different protein classes energy terms may have different effects.
△ Less
Submitted 10 January, 2014; v1 submitted 25 December, 2013;
originally announced December 2013.
-
Genome-wide scan of 29,141 African Americans finds no evidence of selection since admixture
Authors:
Gaurav Bhatia,
Arti Tandon,
Melinda C. Aldrich,
Christine B. Ambrosone,
Christopher Amos,
Elisa V. Bandera,
Sonja I. Berndt,
Leslie Bernstein,
William J. Blot,
Cathryn H. Bock,
Neil Caporaso,
Graham Casey,
Sandra L. Deming,
W. Ryan Diver,
Susan M. Gapstur,
Elizabeth M. Gillanders,
Curtis C. Harris,
Brian E. Henderson,
Sue A. Ingles,
William Isaacs,
Esther M. John,
Rick A. Kittles,
Emma Larkin,
Lorna H. McNeill,
Robert C. Millikan
, et al. (22 additional authors not shown)
Abstract:
We scanned through the genomes of 29,141 African Americans, searching for loci where the average proportion of African ancestry deviates significantly from the genome-wide average. We failed to find any genome-wide significant deviations, and conclude that any selection in African Americans since admixture is sufficiently weak that it falls below the threshold of our power to detect it using a lar…
▽ More
We scanned through the genomes of 29,141 African Americans, searching for loci where the average proportion of African ancestry deviates significantly from the genome-wide average. We failed to find any genome-wide significant deviations, and conclude that any selection in African Americans since admixture is sufficiently weak that it falls below the threshold of our power to detect it using a large sample size. These results stand in contrast to the findings of a recent study of selection in African Americans. That study, which had 15 times fewer samples, reported six loci with significant deviations. We show that the discrepancy is likely due to insufficient correction for multiple hypothesis testing in the previous study. The same study reported 14 loci that showed greater population differentiation between African Americans and Nigerian Yoruba than would be expected in the absence of natural selection. Four such loci were previously shown to be genome-wide significant and likely to be affected by selection, but we show that most of the 10 additional loci are likely to be false positives. Additionally, the most parsimonious explanation for the loci that have significant evidence of unusual differentiation in frequency between Nigerians and Africans Americans is selection in Africa prior to their forced migration to the Americas.
△ Less
Submitted 10 December, 2013;
originally announced December 2013.
-
Looking for packing units of the protein structure
Authors:
Wei-Mou Zheng,
Hui Zeng,
Dong-Bo Bu,
Ming-Fu Shao,
Ke-Song Liu,
Chao Wang
Abstract:
Lattice-model simulations and experiments of some small proteins suggest that folding is essentially controlled by a few conserved contacts. Residues of these conserved contacts form the minimum set of native contacts needed to ensure foldability. Keeping such conserved specific contacts in mind, we examine contacts made by two secondary structure elements of different helices or sheets and look f…
▽ More
Lattice-model simulations and experiments of some small proteins suggest that folding is essentially controlled by a few conserved contacts. Residues of these conserved contacts form the minimum set of native contacts needed to ensure foldability. Keeping such conserved specific contacts in mind, we examine contacts made by two secondary structure elements of different helices or sheets and look for possible 'packing units' of the protein structure. Two short backbone fragments of width five centred at the C? atoms in contact is called an H-form, which serves as a candidate for the packing units. The structural alignment of protein family members or even across families indicates that there are conservative H-forms which are similar both in their sequences and local geometry, and consistent with the structural alignment. Carrying strong sequence signals, such packing units would provide 3D constraints as a complement of the potential functions for the structure prediction.
△ Less
Submitted 13 September, 2011;
originally announced September 2011.
-
Evidence Of Protein Collective Motions On The Picosecond Time Scale
Authors:
Yunfen He,
Jing-Yin Chen,
Joseph R. Knab,
Wenjun Zheng,
Andrea G. Markelz
Abstract:
We investigate the presence of structural collective motions on a picosecond time scale for the heme protein, cytochrome c, as a function of oxidation and hydration, using terahertz (THz) time-domain spectroscopy and molecular dynamics simulations. The THz response dramatically increases with oxidation, with the largest increase for lowest hydrations and highest frequencies. For both oxidation sta…
▽ More
We investigate the presence of structural collective motions on a picosecond time scale for the heme protein, cytochrome c, as a function of oxidation and hydration, using terahertz (THz) time-domain spectroscopy and molecular dynamics simulations. The THz response dramatically increases with oxidation, with the largest increase for lowest hydrations and highest frequencies. For both oxidation states the THz response rapidly increases with hydration saturating above ~25% (g H2O/g protein). Quasi-harmonic vibrational modes and dipole-dipole correlation functions are calculated from molecular dynamics trajectories. The collective mode density of states alone reproduces the measured hydration dependence providing strong evidence of the existence of these motions. The large oxidation dependence is reproduced only by the dipole-dipole correlation function, indicating the contrast arises from diffusive motions consistent with structural changes occurring in the vicinity of a buried internal water molecule.
△ Less
Submitted 7 June, 2010;
originally announced June 2010.
-
CLePAPS: Fast Pair Alignment of Protein Structures Based on Conformational Letters
Authors:
Sheng Wang,
Wei-Mou Zheng
Abstract:
Fast, efficient and reliable algorithms for pairwise alignment of protein structures are in ever increasing demand for analyzing the rapidly growing data of protein structures. CLePAPS is a tool developed for this purpose. It distinguishes itself from other existing algorithms by the use of conformational letters, which are discretized states of 3D segmental structural states. A letter correspon…
▽ More
Fast, efficient and reliable algorithms for pairwise alignment of protein structures are in ever increasing demand for analyzing the rapidly growing data of protein structures. CLePAPS is a tool developed for this purpose. It distinguishes itself from other existing algorithms by the use of conformational letters, which are discretized states of 3D segmental structural states. A letter corresponds to a cluster of combinations of the three angles formed by C_alpha pseudobonds of four contiguous residues. A substitution matrix called CLESUM is available to measure similarity between any two such letters. CLePAPS regards an aligned fragment pair (AFP) as an ungapped string pair with a high sum of pairwise CLESUM scores. Using CLESUM scores as the similarity measure, CLePAPS searches for AFPs by simple string comparison. The transformation which best superimposes a highly similar AFP can be used to superimpose the structure pairs under comparison. A highly scored AFP which is consistent with several other AFPs determines an initial alignment. CLePAPS then joins consistent AFPs guided by their similarity scores to extend the alignment by several `zoom-in' iteration steps. A follow-up refinement produces the final alignment. CLePAPS does not implement dynamic programming. The utility of CLePAPS is tested on various protein structure pairs.
△ Less
Submitted 31 October, 2007;
originally announced October 2007.
-
Dynamic Studies of Scaffold-dependent Mating Pathway in Yeast
Authors:
Danying Shao,
Wen Zheng,
Wenjun Qiu,
Qi Ouyang,
Chao Tang
Abstract:
The mating pathway in \emph{Saccharomyces cerevisiae} is one of the best understood signal transduction pathways in eukaryotes. It transmits the mating signal from plasma membrane into the nucleus through the G-protein coupled receptor and the mitogen-activated protein kinase (MAPK) cascade. According to the current understandings of the mating pathway, we construct a system of ordinary differen…
▽ More
The mating pathway in \emph{Saccharomyces cerevisiae} is one of the best understood signal transduction pathways in eukaryotes. It transmits the mating signal from plasma membrane into the nucleus through the G-protein coupled receptor and the mitogen-activated protein kinase (MAPK) cascade. According to the current understandings of the mating pathway, we construct a system of ordinary differential equations to describe the process. Our model is consistent with a wide range of experiments, indicating that it captures some main characteristics of the signal transduction along the pathway. Investigation with the model reveals that the shuttling of the scaffold protein and the dephosphorylation of kinases involved in the MAPK cascade cooperate to regulate the response upon pheromone induction and to help preserving the fidelity of the mating signaling. We explored factors affecting the dose-response curves of this pathway and found that both negative feedback and concentrations of the proteins involved in the MAPK cascade play crucial role. Contrary to some other MAPK systems where signaling sensitivity is being amplified successively along the cascade, here the mating signal is transmitted through the cascade in an almost linear fashion.
△ Less
Submitted 11 September, 2006;
originally announced September 2006.
-
A protein structural alphabet and its substitution matrix CLESUM
Authors:
Wei-Mou Zheng,
Xin Liu
Abstract:
By using a mixture model for the density distribution of the three pseudobond angles formed by $C_α$ atoms of four consecutive residues, the local structural states are discretized as 17 conformational letters of a protein structural alphabet. This coarse-graining procedure converts a 3D structure to a 1D code sequence. A substitution matrix between these letters is constructed based on the stru…
▽ More
By using a mixture model for the density distribution of the three pseudobond angles formed by $C_α$ atoms of four consecutive residues, the local structural states are discretized as 17 conformational letters of a protein structural alphabet. This coarse-graining procedure converts a 3D structure to a 1D code sequence. A substitution matrix between these letters is constructed based on the structural alignments of the FSSP database.
△ Less
Submitted 27 December, 2004;
originally announced December 2004.
-
Amino acid substitution matrices for protein conformation identification
Authors:
Xin Liu,
Wei-Mou Zheng
Abstract:
Methods for alignment of protein sequences typically measure similarity by using substitution matrix with scores for all possible exchanges of one amino acid with another. Although widely used, the matrices derived from homologous sequence segments, such as Dayhoff's PAM matrices and Henikoff's BLOSUM matrices, are not specific for protein conformation identification. Using a different approach,…
▽ More
Methods for alignment of protein sequences typically measure similarity by using substitution matrix with scores for all possible exchanges of one amino acid with another. Although widely used, the matrices derived from homologous sequence segments, such as Dayhoff's PAM matrices and Henikoff's BLOSUM matrices, are not specific for protein conformation identification. Using a different approach, we got many amino acid segment blocks. For each of them, the protein secondary structure is identical. Based on these blocks, we have derived new amino acid substitution matrices. The application of these matrices led to marked improvements in conformation segment search and homologues detection in twilight zone.
△ Less
Submitted 15 June, 2004;
originally announced June 2004.
-
Protein secondary structure prediction by combining hidden Markov models and sliding window scores
Authors:
Wei-Mou Zheng
Abstract:
Instead of conformation states of single residues, refined conformation states of quintuplets are proposed to reflect conformation correlation. Simple hidden Markov models combining with sliding window scores are used for predicting secondary structure of a protein from its amino acid sequence. Since the length of protein conformation segments varies in a narrow range, we ignore the duration eff…
▽ More
Instead of conformation states of single residues, refined conformation states of quintuplets are proposed to reflect conformation correlation. Simple hidden Markov models combining with sliding window scores are used for predicting secondary structure of a protein from its amino acid sequence. Since the length of protein conformation segments varies in a narrow range, we ignore the duration effect of the length distribution. The window scores for residues are a window version of the Chou-Fasman propensities estimated under an approximation of conditional independency. Different window widths are examined, and the optimal width is found to be 17. A high accuracy about 70% is achieved.
△ Less
Submitted 20 October, 2003;
originally announced October 2003.
-
Protein secondary structure prediction based on quintuplets
Authors:
Wei-Mou Zheng
Abstract:
Simple hidden Markov models are proposed for predicting secondary structure of a protein from its amino acid sequence. Since the length of protein conformation segments varies in a narrow range, we ignore the duration effect of length distribution, and focus on inclusion of short range correlations of residues and of conformation states in the models. Conformation-independent and -dependent amin…
▽ More
Simple hidden Markov models are proposed for predicting secondary structure of a protein from its amino acid sequence. Since the length of protein conformation segments varies in a narrow range, we ignore the duration effect of length distribution, and focus on inclusion of short range correlations of residues and of conformation states in the models. Conformation-independent and -dependent amino acid coarse-graining schemes are designed for the models by means of proper mutual information. We compare models of different level of complexity, and establish a practical model with a high prediction accuracy.
△ Less
Submitted 15 July, 2003;
originally announced July 2003.
-
Prediction of protein secondary structure based on residue pairs
Authors:
Xin Liu,
Li-mei Zhang,
Wei-Mou Zheng
Abstract:
The GOR program for predicting protein secondary structure is extended to include triple correlation. A score system for a residue pair to be at certain conformation state is derived from the conditional weight matrix describing amino acid frequencies at each position of a window flanking the pair under the condition for the pair to be at the fixed state. A program using this score system to pre…
▽ More
The GOR program for predicting protein secondary structure is extended to include triple correlation. A score system for a residue pair to be at certain conformation state is derived from the conditional weight matrix describing amino acid frequencies at each position of a window flanking the pair under the condition for the pair to be at the fixed state. A program using this score system to predict protein secondary structure is established. After training the model with a learning set created from PDB\_SELECT, the program is tested with two test sets. As a method using single sequence for predicting secondary structures, the approach achieves a high accuracy near 70%.
△ Less
Submitted 15 June, 2004; v1 submitted 18 December, 2002;
originally announced December 2002.
-
Distances and classification of amino acids for different protein secondary structures
Authors:
Xin Liu,
Li-mei Zhang,
Shan Guan,
Wei-Mou Zheng
Abstract:
Window profiles of amino acids in protein sequences are taken as a description of the amino acid environment. The relative entropy or Kullback-Leibler distance derived from profiles is used as a measure of dissimilarity for comparison of amino acids and secondary structure conformations. Distance matrices of amino acid pairs at different conformations are obtained, which display a non-negligible…
▽ More
Window profiles of amino acids in protein sequences are taken as a description of the amino acid environment. The relative entropy or Kullback-Leibler distance derived from profiles is used as a measure of dissimilarity for comparison of amino acids and secondary structure conformations. Distance matrices of amino acid pairs at different conformations are obtained, which display a non-negligible dependence of amino acid similarity on conformations. Based on the conformation specific distances clustering analysis for amino acids is conducted.
△ Less
Submitted 16 June, 2003; v1 submitted 7 November, 2002;
originally announced November 2002.
-
Simplified amino acid alphabets based on deviation of conditional probability from random background
Authors:
Xin Liu,
Di Liu,
Ji Qi,
Wei-Mou Zheng
Abstract:
The primitive data for deducing the Miyazawa-Jernigan contact energy or BLOSUM score matrix consists of pair frequency counts. Each amino acid corresponds to a conditional probability distribution. Based on the deviation of such conditional probability from random background, a scheme for reduction of amino acid alphabet is proposed. It is observed that evident discrepancy exists between reduced…
▽ More
The primitive data for deducing the Miyazawa-Jernigan contact energy or BLOSUM score matrix consists of pair frequency counts. Each amino acid corresponds to a conditional probability distribution. Based on the deviation of such conditional probability from random background, a scheme for reduction of amino acid alphabet is proposed. It is observed that evident discrepancy exists between reduced alphabets obtained from raw data of the Miyazawa-Jernigan's and BLOSUM's residue pair counts. Taking homologous sequence database SCOP40 as a test set, we detect homology with the obtained coarse-grained substitution matrices. It is verified that the reduced alphabets obtained well preserve information contained in the original 20-letter alphabet.
△ Less
Submitted 6 November, 2002;
originally announced November 2002.
-
Entropic Approach for Reduction of Amino Acid Alphabets
Authors:
Wei-Mou Zheng
Abstract:
The primitive data for deducing the Miyazawa-Jernigan contact energy or BLOSUM score metrix are the pair frequency counts. Each amino acid corresponds to a distribution. Taking the Kullback-Leibler distance of two probability distributions as resemblance coefficient and relating cluster to mixed population, we perform cluster analysis of amino acids based on the frequecy counts data. Furthermore…
▽ More
The primitive data for deducing the Miyazawa-Jernigan contact energy or BLOSUM score metrix are the pair frequency counts. Each amino acid corresponds to a distribution. Taking the Kullback-Leibler distance of two probability distributions as resemblance coefficient and relating cluster to mixed population, we perform cluster analysis of amino acids based on the frequecy counts data. Furthermore, Ward's clustering is also obtained by adopting the average score as an objective function. An ordinal cophenetic is introduced to compare results from different clustering methods.
△ Less
Submitted 22 June, 2001;
originally announced June 2001.
-
Self-Organizing Approach for Finding Borders of DNA Coding Regions
Authors:
Fang Wu,
Wei-Mou Zheng
Abstract:
A self-organizing approach is proposed for gene finding based on the model of codon usage for coding regions and positional preference for noncoding regions. The symmetry between the direct and reverse coding regions is adopted for reducing the number of parameters. Without requiring prior training, parameters are estimated by iteration. By employing the window sliding technique and likelihood r…
▽ More
A self-organizing approach is proposed for gene finding based on the model of codon usage for coding regions and positional preference for noncoding regions. The symmetry between the direct and reverse coding regions is adopted for reducing the number of parameters. Without requiring prior training, parameters are estimated by iteration. By employing the window sliding technique and likelihood ratio, a very accurate segmentation is obtained.
△ Less
Submitted 16 February, 2001;
originally announced February 2001.
