A protein evolution model with independent sites that reproduces site-specific amino acid distributions from the Protein Data Bank

被引:36
作者
Bastolla, Ugo [1 ]
Porto, Markus
Roman, H. Eduardo
Vendruscolo, Michele
机构
[1] Univ Autonoma Madrid, CSIC, Ctr Biol Mol Severo Ochoa, E-28049 Madrid, Spain
[2] Tech Univ Darmstadt, Inst Festkorperphys, D-64289 Darmstadt, Germany
[3] Univ Milano Bicocca, Dipartimento Fis, I-20126 Milan, Italy
[4] Univ Cambridge, Dept Chem, Cambridge CB2 1EW, England
关键词
D O I
10.1186/1471-2148-6-43
中图分类号
Q [生物科学];
学科分类号
07 ; 0710 ; 09 ;
摘要
Background: Since thermodynamic stability is a global property of proteins that has to be conserved during evolution, the selective pressure at a given site of a protein sequence depends on the amino acids present at other sites. However, models of molecular evolution that aim at reconstructing the evolutionary history of macromolecules become computationally intractable if such correlations between sites are explicitly taken into account. Results: We introduce an evolutionary model with sites evolving independently under a global constraint on the conservation of structural stability. This model consists of a selection process, which depends on two hydrophobicity parameters that can be computed from protein sequences without any fit, and a mutation process for which we consider various models. It reproduces quantitatively the results of Structurally Constrained Neutral (SCN) simulations of protein evolution in which the stability of the native state is explicitly computed and conserved. We then compare the predicted site-specific amino acid distributions with those sampled from the Protein Data Bank (PDB). The parameters of the mutation model, whose number varies between zero and five, are fitted from the data. The mean correlation coefficient between predicted and observed site-specific amino acid distributions is larger than < r > = 0.70 for a mutation model with no free parameters and no genetic code. In contrast, considering only the mutation process with no selection yields a mean correlation coefficient of < r > = 0.56 with three fitted parameters. The mutation model that best fits the data takes into account increased mutation rate at CpG dinucleotides, yielding < r > = 0.90 with five parameters. Conclusion: The effective selection process that we propose reproduces well amino acid distributions as observed in the protein sequences in the PDB. Its simplicity makes it very promising for likelihood calculations in phylogenetic studies. Interestingly, in this approach the mutation process influences the effective selection process, i. e. selection and mutation must be entangled in order to obtain effectively independent sites. This interdependence between mutation and selection reflects the deep influence that mutation has on the evolutionary process: The bias in the mutation influences the thermodynamic properties of the evolving proteins, in agreement with comparative studies of bacterial proteomes, and it also influences the rate of accepted mutations.
引用
收藏
页数:21
相关论文
共 92 条
[1]   FREE-ENERGY LANDSCAPE FOR PROTEIN-FOLDING KINETICS - INTERMEDIATES, TRAPS, AND MULTIPLE PATHWAYS IN THEORY AND LATTICE MODEL SIMULATIONS [J].
ABKEVICH, VI ;
GUTIN, AM ;
SHAKHNOVICH, EI .
JOURNAL OF CHEMICAL PHYSICS, 1994, 101 (07) :6052-6062
[2]   Roles of molecular chaperones in cytoplasmic protein folding [J].
Agashe, VR ;
Hartl, FU .
SEMINARS IN CELL & DEVELOPMENTAL BIOLOGY, 2000, 11 (01) :15-25
[3]  
[Anonymous], 1998, Science
[4]  
[Anonymous], 1991, The Causes of Molecular Evolution
[5]   Neutral networks in protein space: a computational study based on knowledge-based potentials of mean force [J].
Babajide, A ;
Hofacker, IL ;
Sippl, MJ ;
Stadler, PF .
FOLDING & DESIGN, 1997, 2 (05) :261-269
[6]   Neutral evolution of model proteins: Diffusion in sequence space and overdispersion [J].
Bastolla, U ;
Roman, HE ;
Vendruscolo, M .
JOURNAL OF THEORETICAL BIOLOGY, 1999, 200 (01) :49-64
[7]   Looking at structure, stability, and evolution of proteins through the principal eigenvector of contact matrices and hydrophobicity profiles [J].
Bastolla, U ;
Porto, M ;
Roman, HE ;
Vendrusolo, M .
GENE, 2005, 347 (02) :219-230
[8]   Principal eigenvector of contact matrices and hydrophobicity profiles in proteins [J].
Bastolla, U ;
Porto, M ;
Roman, HE ;
Vendruscolo, M .
PROTEINS-STRUCTURE FUNCTION AND BIOINFORMATICS, 2005, 58 (01) :22-30
[9]   Stability constraints and protein evolution: the role of chain length, composition and disulfide bonds [J].
Bastolla, U ;
Demetrius, L .
PROTEIN ENGINEERING DESIGN & SELECTION, 2005, 18 (09) :405-415
[10]   Genomic determinants of protein folding thermodynamics in prokaryotic organisms [J].
Bastolla, U ;
Moya, A ;
Viguera, E ;
van Ham, RCHJ .
JOURNAL OF MOLECULAR BIOLOGY, 2004, 343 (05) :1451-1466