Organizational Unit:
Center for the Study of Systems Biology

Research Organization Registry ID
Description
Previous Names
Parent Organization
Parent Organization
Organizational Unit
Includes Organization(s)
ArchiveSpace Name Record

Publication Search Results

Now showing 1 - 6 of 6
Thumbnail Image
Item

Unfolding of globular proteins: monte carlo dynamics of a realistic reduced model

2003-11 , Kolinski, Andrzej , Klein, Piotr , Romiszowski, Piotr , Skolnick, Jeffrey

Reduced lattice models of proteins and Monte Carlo dynamics were used to simulate the initial stages of the unfolding of several proteins of various structural types, and the results were compared to experiment. The models semiquantitatively reproduce the approximate order of events of unfolding as well as subtle mutation effects and effects resulting from differences in sequences of similar folds. The short-time mobility of particular residues, observed in simulations, correlates with the crystallographic temperature factor. The main factor controlling unfolding is the native state topology, with sequence playing a less important role. The correlation with various experiments, especially for sequence-specific effects, strongly suggests that properly designed reduced models of proteins can be used for qualitative studies (or prediction) of protein unfolding pathways

Thumbnail Image
Item

Numerical study of the entropy loss of dimerization and the folding thermodynamics of the GCN4 leucine zipper

2002-11 , Viñals, Jorge , Kolinski, Andrzej , Skolnick, Jeffrey

A lattice-based model of a protein and the Monte Carlo simulation method are used to calculate the entropy loss of dimerization of the GCN4 leucine zipper. In the representation used, a protein is a sequence of interaction centers arranged on a cubic lattice, with effective interaction potentials that are both of physical and statistical nature. The Monte Carlo simulation method is then used to sample the partition functions of both the monomer and dimer forms as a function of temperature. A method is described to estimate the entropy loss upon dimerization, a quantity that enters the free energy difference between monomer and dimer, and the corresponding dimerization reaction constant. As expected, but contrary to previous numerical studies, we find that the entropy loss of dimerization is a strong function of energy (or temperature), except in the limit of large energies in which the motion of the two dimer chains becomes largely uncorrelated. At the monomer-dimer transition temperature we find that the entropy loss of dimerization is approximately five times smaller than the value that would result from ideal gas statistics, a result that is qualitatively consistent with a recent experimental determination of the entropy loss of dimerization of a synthetic peptide that also forms a two-stranded -helical coiled coil.

Thumbnail Image
Item

TOUCHSTONE II: a new approach to ab initio protein structure prediction

2003-08 , Zhang, Yang , Kolinski, Andrzej , Skolnick, Jeffrey

We have developed a new combined approach for ab initio protein structure prediction. The protein conformation is described as a lattice chain connecting Ca atoms, with attached Cb atoms and side-chain centers of mass. The model force field includes various short-range and long-range knowledge-based potentials derived from a statistical analysis of the regularities of protein structures. The combination of these energy terms is optimized through the maximization of correlation for 30 3 60,000 decoys between the root mean square deviation (RMSD) to native and energies, as well as the energy gap between native and the decoy ensemble. To accelerate the conformational search, a newly developed parallel hyperbolic sampling algorithm with a composite movement set is used in the Monte Carlo simulation processes. We exploit this strategy to successfully fold 41/100 small proteins (36 ; 120 residues) with predicted structures having a RMSD from native below 6.5 A˚ in the top five cluster centroids. To fold larger-size proteins as well as to improve the folding yield of small proteins, we incorporate into the basic force field side-chain contact predictions from our threading program PROSPECTOR where homologous proteins were excluded from the data base. With these threading-based restraints, the program can fold 83/125 test proteins (36 ; 174 residues) with structures having a RMSD to native below 6.5 A˚ in the top five cluster centroids. This shows the significant improvement of folding by using predicted tertiary restraints, especially when the accuracy of side-chain contact prediction is [20%. For native fold selection, we introduce quantities dependent on the cluster density and the combination of energy and free energy, which show a higher discriminative power to select the native structure than the previously used cluster energy or cluster size, and which can be used in native structure identification in blind simulations. These procedures are readily automated and are being implemented on a genomic scale.

Thumbnail Image
Item

Computational studies of protein folding

2001-09 , Skolnick, Jeffrey , Kolinski, Andrzej

The authors describe the state of the art in the field of protein structure prediction. They also introduce Prospector, a newly developed, iterative threading algorithm for protein structure prediction that can also be applied to ab initio protein folding, and discuss the promising results of its large-scale application.

Thumbnail Image
Item

A minimal physically realistic protein-like lattice model: designing an energy landscape that ensures all-or-none folding to a unique native state

2003-03 , Pokarowski, Piotr , Kolinski, Andrzej , Skolnick, Jeffrey

A simple protein model restricted to the face-centered cubic lattice has been studied. The model interaction scheme includes attractive interactions between hydrophobic (H) residues, repulsive interactions between hydrophobic and polar (P) residues, and orientation-dependent P-P interactions. Additionally, there is a potential that favors extended b-type conformations. A sequence has been designed that adopts a native structure, consisting of an antiparallel, six-member Greekkey b-barrel with protein-like structural degeneracy. It has been shown that the proposed model is a minimal one, i.e., all the above listed types of interactions are necessary for cooperative (all-or-none) type folding to the native state. Simulations were performed via the Replica Exchange Monte Carlo method and the numerical data analyzed via a multihistogram method.

Thumbnail Image
Item

Comparison of three Monte Carlo conformational search strategies for a proteinlike homopolymer model: Folding thermodynamics and identification of low-energy structures

2000-09-22 , Gront, Dominik , Kolinski, Andrzej , Skolnick, Jeffrey

Entropy sampling Monte Carlo, the replica method, and the classical Metropolis scheme were applied in numerical studies of the collapse transition in a simple face-centered cubic lattice polymer. The force field of the model consists of pairwise, contact-type, long-range interactions and a short-range potential based on the β -sheet definition assumed in the model. The ability to find the lowest energy conformation by various Monte Carlo methods and the computational cost associated with each was examined. It is shown that all of the methods generally provide the same picture of the collapse transition. However, the most complete thermodynamic description of the transition derives from the results of entropy sampling Monte Carlo simulations, but this is the most time-consuming method. The replica method is shown to be the most effective and efficient in searching for the lowest energy conformation. The possible consequences of these findings for the development of simulation strategies for the folding of model proteins are discussed briefly.