1,418 research outputs found
Origin of Native Driving Force in Protein Folding
We derive an expression with four adjustable parameters that reproduces well
the 20x20 Miyazawa-Jernigan potential matrix extracted from known protein
structures. The numerical values of the parameters can be approximately
computed from the surface tension of water, water-screened dipole interactions
between residues and water and among residues, and average exposures of
residues in folded proteins.Comment: LaTeX file, Postscript file; 4 pages, 1 figure (mij.eps), 2 table
An Analytical Approach to the Protein Designability Problem
We present an analytical method for determining the designability of protein
structures. We apply our method to the case of two-dimensional lattice
structures, and give a systematic solution for the spectrum of any structure.
Using this spectrum, the designability of a structure can be estimated. We
outline a heirarchy of structures, from most to least designable, and show that
this heirarchy depends on the potential that is used.Comment: 16 pages 4 figure
Geometrically Reduced Number of Protein Ground State Candidates
Geometrical properties of protein ground states are studied using an
algebraic approach. It is shown that independent from inter-monomer
interactions, the collection of ground state candidates for any folded protein
is unexpectedly small: For the case of a two-parameter Hydrophobic-Polar
lattice model for -mers, the number of these candidates grows only as .
Moreover, the space of the interaction parameters of the model breaks up into
well-defined domains, each corresponding to one ground state candidate, which
are separated by sharp boundaries. In addition, by exact enumeration, we show
there are some sequences which have one absolute unique native state. These
absolute ground states have perfect stability against change of inter-monomer
interaction potential.Comment: 9 page, 4 ps figures are include
Design of Force Fields from Data at Finite Temperature
We investigate the problem of how to obtain the force field between atoms of
an experimentally determined structure. We show how this problem can be
efficiently solved, even at finite temperature, where the position of the atoms
differs substantially from the ground state. We apply our method to systems
modeling proteins and demonstrate that the correct potentials can be recovered
even in the presence of thermal noise.Comment: 10 pages, 1 postcript figure, Late
Freezing Transition of Random Heteropolymers Consisting of an Arbitrary Set of Monomers
Mean field replica theory is employed to analyze the freezing transition of
random heteropolymers comprised of an arbitrary number () of types of
monomers. Our formalism assumes that interactions are short range and
heterogeneity comes only from pairwise interactions, which are defined by an
arbitrary matrix. We show that, in general, there exists a
freezing transition from a random globule, in which the thermodynamic
equilibrium is comprised of an essentially infinite number polymer
conformations, to a frozen globule, in which equilibrium ensemble is dominated
by one or very few conformations. We also examine some special cases of
interaction matrices to analyze the relationship between the freezing
transition and the nature of interactions involved.Comment: 30 pages, 1 postscript figur
Is Heteropolymer Freezing Well Described by the Random Energy Model?
It is widely held that the Random Energy Model (REM) describes the freezing
transition of a variety of types of heteropolymers. We demonstrate that the
hallmark property of REM, statistical independence of the energies of states
over disorder, is violated in different ways for models commonly employed in
heteropolymer freezing studies. The implications for proteins are also
discussed.Comment: 4 pages, 3 eps figures To appear in Physical Review Letters, May 199
Modeling study on the validity of a possibly simplified representation of proteins
The folding characteristics of sequences reduced with a possibly simplified
representation of five types of residues are shown to be similar to their
original ones with the natural set of residues (20 types or 20 letters). The
reduced sequences have a good foldability and fold to the same native structure
of their optimized original ones. A large ground state gap for the native
structure shows the thermodynamic stability of the reduced sequences. The
general validity of such a five-letter reduction is further studied via the
correlation between the reduced sequences and the original ones. As a
comparison, a reduction with two letters is found not to reproduce the native
structure of the original sequences due to its homopolymeric features.Comment: 6 pages with 4 figure
Selective Constraints on Amino Acids Estimated by a Mechanistic Codon Substitution Model with Multiple Nucleotide Changes
Empirical substitution matrices represent the average tendencies of
substitutions over various protein families by sacrificing gene-level
resolution. We develop a codon-based model, in which mutational tendencies of
codon, a genetic code, and the strength of selective constraints against amino
acid replacements can be tailored to a given gene. First, selective constraints
averaged over proteins are estimated by maximizing the likelihood of each 1-PAM
matrix of empirical amino acid (JTT, WAG, and LG) and codon (KHG) substitution
matrices. Then, selective constraints specific to given proteins are
approximated as a linear function of those estimated from the empirical
substitution matrices.
Akaike information criterion (AIC) values indicate that a model allowing
multiple nucleotide changes fits the empirical substitution matrices
significantly better. Also, the ML estimates of transition-transversion bias
obtained from these empirical matrices are not so large as previously
estimated. The selective constraints are characteristic of proteins rather than
species. However, their relative strengths among amino acid pairs can be
approximated not to depend very much on protein families but amino acid pairs,
because the present model, in which selective constraints are approximated to
be a linear function of those estimated from the JTT/WAG/LG/KHG matrices, can
provide a good fit to other empirical substitution matrices including cpREV for
chloroplast proteins and mtREV for vertebrate mitochondrial proteins.
The present codon-based model with the ML estimates of selective constraints
and with adjustable mutation rates of nucleotide would be useful as a simple
substitution model in ML and Bayesian inferences of molecular phylogenetic
trees, and enables us to obtain biologically meaningful information at both
nucleotide and amino acid levels from codon and protein sequences.Comment: Table 9 in this article includes corrections for errata in the Table
9 published in 10.1371/journal.pone.0017244. Supporting information is
attached at the end of the article, and a computer-readable dataset of the ML
estimates of selective constraints is available from
10.1371/journal.pone.001724
- …