Search PubMed⌕ Search

Biomedical subjects

A Kolinski

Publications and source records attributed to A Kolinski.

53 records · Page 3Linked to original sources
PubMed1994

Prediction of the folding pathways and structure of the GCN4 leucine zipper.

A hierarchical approach is described for the prediction of the three-dimensional structure and folding pathway of the GCN4 leucine zipper. Dimer assembly is simulated by Monte Carlo dynamics. The resulting lowest energy structures undergo cooperative rearrangement of their hydrophobic core leading to side-chain fixation. The coarse-grained structures are further refined using a molecular dynamics annealing protocol. This produces full atom models with a backbone root-mean-square deviation from the crystal structure of 0.81 A. Thus, we demonstrate the predictive ability of our approach to yield high resolution structures of small coiled coils from their sequence.

Amino Acid Sequence↗
PubMed1994

Monte Carlo simulations of protein folding. I. Lattice model and interaction scheme.

A new hierarchical method for the simulation of the protein folding process and the de novo prediction of protein three-dimensional structure is proposed. The reduced representation of the protein alpha-carbon backbone employs lattice discretizations of increasing geometrical resolution and a single ball representation of side chain rotamers. In particular, coarser and finer lattice backbone descriptions are used. The coarser (finer) lattice represents C alpha traces of native proteins with an accuracy of 1.0 (0.7) A rms. Folding is simulated by means of very fast Monte Carlo lattice dynamics. The potential of mean force, predominantly of statistical origin, contains several novel terms that facilitate the cooperative assembly of secondary structure elements and the cooperative packing of the side chains. Particular contributions to the interaction scheme are discussed in detail. In the accompanying paper (Kolinski, A., Skolnick, J. Monte Carlo simulation of protein folding. II. Application to protein A, ROP, and crambin. Proteins 18:353-366, 1994), the method is applied to three small globular proteins.

Amino Acid Sequence↗
PubMed1994

Monte Carlo simulations of protein folding. II. Application to protein A, ROP, and crambin.

The hierarchy of lattice Monte Carlo models described in the accompanying paper (Kolinski, A., Skolnick, J. Monte Carlo simulations of protein folding. I. Lattice model and interaction scheme. Proteins 18:338-352, 1994) is applied to the simulation of protein folding and the prediction of 3-dimensional structure. Using sequence information alone, three proteins have been successfully folded: the B domain of staphylococcal protein A, a 120 residue, monomeric version of ROP dimer, and crambin. Starting from a random expanded conformation, the model proteins fold along relatively well-defined folding pathways. These involve a collection of early intermediates, which are followed by the final (and rate-determining) transition from compact intermediates closely resembling the molten globule state to the native-like state. The predicted structures are rather unique, with native-like packing of the side chains. The accuracy of the predicted native conformations is better than those obtained in previous folding simulations. The best (but by no means atypical) folds of protein A have a coordinate rms of 2.25 A from the native C alpha trace, and the best coordinate rms from crambin is 3.18 A. For ROP monomer, the lowest coordinate rms from equivalent C alpha s of ROP dimer is 3.65 A. Thus, for two simple helical proteins and a small alpha/beta protein, the ability to predict protein structure from sequence has been demonstrated.

Amino Acid Sequence↗
PubMed1993

A method for predicting protein structure from sequence.

BACKGROUND: The ability to predict the native conformation of a globular protein from its amino-acid sequence is an important unsolved problem of molecular biology. We have previously reported a method in which reduced representations of proteins are folded on a lattice by Monte Carlo simulation, using statistically-derived potentials. When applied to sequences designed to fold into four-helix bundles, this method generated predicted conformations closely resembling the real ones. RESULTS: We now report a hierarchical approach to protein-structure prediction, in which two cycles of the above-mentioned lattice method (the second on a finer lattice) are followed by a full-atom molecular dynamics simulation. The end product of the simulations is thus a full-atom representation of the predicted structure. The application of this procedure to the 60 residue, B domain of staphylococcal protein A predicts a three-helix bundle with a backbone root mean square (rms) deviation of 2.25-3 A from the experimentally determined structure. Further application to a designed, 120 residue monomeric protein, mROP, based on the dimeric ROP protein of Escherichia coli, predicts a left turning, four-helix bundle native state. Although the ultimate assessment of the quality of this prediction awaits the experimental determination of the mROP structure, a comparison of this structure with the set of equivalent residues in the ROP dime- crystal structure indicates that they have a rms deviation of approximately 3.6-4.2 A. CONCLUSION: Thus, for a set of helical proteins that have simple native topologies, the native folds of the proteins can be predicted with reasonable accuracy from their sequences alone. Our approach suggest a direction for future work addressing the protein-folding problem.

Journal Article↗
PubMed1993

De novo and inverse folding predictions of protein structure and dynamics.

In the last two years, the use of simplified models has facilitated major progress in the globular protein folding problem, viz., the prediction of the three-dimensional (3D) structure of a globular protein from its amino acid sequence. A number of groups have addressed the inverse folding problem where one examines the compatibility of a given sequence with a given (and already determined) structure. A comparison of extant inverse protein-folding algorithms is presented, and methodologies for identifying sequences likely to adopt identical folding topologies, even when they lack sequence homology, are described. Extension to produce structural templates or fingerprints from idealized structures is discussed, and for eight-membered beta-barrel proteins, it is shown that idealized fingerprints constructed from simple topology diagrams can correctly identify sequences having the appropriate topology. Furthermore, this inverse folding algorithm is generalized to predict elements of supersecondary structure including beta-hairpins, helical hairpins and alpha/beta/alpha fragments. Then, we describe a very high coordination number lattice model that can predict the 3D structure of a number of globular proteins de novo; i.e. using just the amino acid sequence. Applications to sequences designed by DeGrado and co-workers [Biophys. J., 61 (1992) A265] predict folding intermediates, native states and relative stabilities in accord with experiment. The methodology has also been applied to the four-helix bundle designed by Richardson and co-workers [Science, 249 (1990) 884] and a redesigned monomeric version of a naturally occurring four-helix dimer, rop. Based on comparison to the rop dimer, the simulations predict conformations with rms values of 3-4 A from native. Furthermore, the de novo algorithms can assess the stability of the folds predicted from the inverse algorithm, while the inverse folding algorithms can assess the quality of the de novo models. Thus, the synergism of the de novo and inverse folding algorithm approaches provides a set of complementary tools that will facilitate further progress on the protein-folding problem.

Algorithms↗
PubMed1993

Regularities in interaction patterns of globular proteins.

The description of protein structure in the language of side chain contact maps is shown to offer many advantages over more traditional approaches. Because it focuses on side chain interactions, it aids in the discovery, study and classification of similarities between interactions defining particular protein folds and offers new insights into the rules of protein structure. For example, there is a small number of characteristic patterns of interactions between protein supersecondary structural fragments, which can be seen in various non-related proteins. Furthermore, the overlap of the side chain contact maps of two proteins provides a new measure of protein structure similarity. As shown in several examples, alignments based on contact map overlaps are a powerful alternative to other structure-based alignments.

Computer Simulation↗
PubMed1992

Topology fingerprint approach to the inverse protein folding problem.

We describe the most general solution to date of the problem of matching globular protein sequences to the appropriate three-dimensional structures. The screening template, against which sequences are tested, is provided by a protein "structural fingerprint" library based on the contact map and the buried/exposed pattern of residues. Then, a lattice Monte Carlo algorithm validates or dismisses the stability of the proposed fold. Examples of known structural similarities between proteins having weakly or unrelated sequences such as the globins and phycocyanins, the eight-member alpha/beta fold of triose phosphate isomerase and even a close structural equivalence between azurin and immunoglobulins are found.

Algorithms↗
PubMed1992

Simulations of the folding pathway of triose phosphate isomerase-type alpha/beta barrel proteins.

Simulations of the folding pathways of two large alpha/beta proteins, the alpha subunit of tryptophan synthase and triose phosphate isomerase, are reported using the knight's walk lattice model of globular proteins and Monte Carlo dynamics. Starting from randomly generated unfolded states and with no assumptions regarding the nature of the folding intermediates, for the tryptophan synthase subunit these simulations predict, in agreement with experiment, the existence and location of a stable equilibrium intermediate comprised of six beta strands on the amino terminus of the molecule. For the case of triose phosphate isomerase, the simulations predict that both amino- and carboxyl-terminal intermediates should be observed. In a significant modification of previous lattice models, this model includes a full heavy atom side chain description and is capable of representing native conformations at the level of 2.5- to 3-A rms deviation for the C alpha positions, as compared to the crystal structure. With a well-balanced compromise between accuracy of the protein description and the computer requirements necessary to perform simulations spanning biologically significant amounts of time, the lattice model described here brings the possibility of studying important biological processes to present-day computers.

Animals↗
PubMed1991

Dynamic Monte Carlo simulations of a new lattice model of globular protein folding, structure and dynamics.

A long-standing problem of molecular biology is the prediction of globular protein tertiary structure from the primary sequence. In the context of a new, 24-nearest-neighbor lattice model of proteins that includes both alpha and beta-carbon atoms, the requirements for folding to a unique four-member beta-barrel, four-helix bundles and a model alpha/beta-bundle have been explored. A number of distinct situations are examined, but the common requirements for the formation of a unique native conformation are tertiary interactions plus the presence of relatively small (but not irrelevant) intrinsic turn preferences that select out the native conformer from a manifold of compact states. When side-chains are explicitly included, there are many conformations having the same or a slightly greater number of side-chain contacts as in the native conformation, and it is the local intrinsic turn preferences that produce the conformational selectivity on collapse. The local preference for helix or beta-sheet secondary structure may be at odds with the secondary structure ultimately found in the native conformation. The requisite intrinsic turn populations are about 0.3% for beta-proteins, 2% for mixed alpha/beta-proteins and 6% for helix bundles. In addition, an idealized model of an allosteric conformational transition has been examined. Folding occurs predominantly by a sequential on-site assembly mechanism with folding initiating either at a turn or from an isolated helix or beta-strand (where appropriate). For helical and beta-protein models, similar folding pathways were obtained in diamond lattice simulations, using an entirely different set of local Monte Carlo moves. This argues strongly that the results are universal; that is, they are independent of lattice, protein model or the particular realization of Monte Carlo dynamics. Overall, these simulations demonstrate that the folding of all known protein motifs can be achieved in the context of a single class of lattice models that includes realistic backbone structures and idealized side-chains.

Algorithms↗
PubMed1990

Dynamic Monte Carlo simulations of globular protein folding/unfolding pathways. I. Six-member, Greek key beta-barrel proteins.

In the context of a simplified diamond lattice model of a six-member, Greek key beta-barrel protein that is closely related in topology to plastocyanin, the nature of the folding and unfolding pathways have been investigated using dynamic Monte Carlo techniques. The mechanism of Greek key assembly is best described as punctuated "on site construction". Folding typically starts at or near a beta-turn, and then the beta-strands sequentially form by using existing folded structure as a scaffold onto which subsequent tertiary structure assembles. On average, beta-strands tend to zip up from one tight bend to the next. After the four-member, beta-barrel assembles, there is a long pause as the random coil portion of the chain containing the long loop thrahes about trying to find the native state. Thus, there is an entropic barrier that must be surmounted. However, while a given piece of the protein may be folding, another section may be unfolding. A competition therefore exists to assemble a fairly stable intermediate before it dissolves. Folding may initiate at any of the tight turns, but the turn closer to the N terminus seems to be preferred due to well-known excluded volume effects. When the protein first starts to fold, there are a multiplicity of folding pathways, but the number of options is reduced as the system gets closer to the native state. In the early stages, the excluded volume effect exerted by the already assembled protein helps subsequent assembly. Then, near the native conformation, the folded parts reduce the accessible conformational space available to the remaining unfolded sections. Unfolding essentially occurs in reverse. Employing a simple statistical mechanical theory, the configurational free energy along the reaction co-ordinate for this model has been constructed. The free energy surface, in agreement with the simulations, provides the following predictions. The transition state is quite near the native state, and consists of five of the six beta-strands being fully assembled, with the remaining long loop plus sixth beta-strand in place, but only partially assembled. It is separated from the beta-barrel intermediate by a free energy barrier of mainly entropic origin and from the native state by a barrier that is primarily energetic in origin. The latter feature is in agreement with the "Cardboard Box" model described by Goldenberg and Creighton but, unlike their model, the transition state is not a high-energy distorted form of the native state.(ABSTRACT TRUNCATED AT 250 WORDS)

Algorithms↗
PubMed1989

Monte Carlo studies on equilibrium globular protein folding. II. Beta-barrel globular protein models.

In the context of dynamic Monte Carlo simulations on a model protein confined to a tetrahedral lattice, the interplay of protein size and tertiary structure, and the requirements for an all-or-none transition to a unique native state, are investigated. Small model proteins having a primary sequence consisting of a central bend neutral region flanked by two tails having an alternating hydrophobic/hydrophilic pattern of residues are seen to undergo a continuous transition to a beta-hairpin collapsed state. On increasing the length of the tails, the beta-hairpin structural motif is found to be in equilibrium with a four-member beta-barrel. Further increase of the tail length results in the shift of the structural equilibrium to the four-member beta-barrel. The random coil to beta-barrel transition is of an all-or-none character, but while the central turn is always the desired native bend, the location of the turns involving the two external strands is variable. That is, beta-barrels having the external stands that are two residues out of register are also observed in the transition region. Introduction into the primary sequence of two additional regions that are at the very least neutral toward turn formation produces an all-or-none transition to the unique, native, four-member beta-barrel. Various factors that can augment the stability of the native conformation are explored. Overall, these folding simulations strongly indicate that the general rules of globular protein folding are rather robust--namely, one requires a general pattern of hydrophobic/hydrophilic residues that allow the protein to have a well-defined interior and exterior and the presence of regions in the amino acid sequence that at the very least are locally indifferent to turn formation. Since no site-specific interactions between hydrophobic and hydrophilic residues are required to produce a unique four-member beta-barrel, these simulations strongly suggest that site specificity is involved in structural fine-tuning.

Algorithms↗
PubMed1989

Dynamic Monte Carlo study of the folding of a six-stranded Greek key globular protein.

To help elucidate the general rules of equilibrium globular protein folding, dynamic Monte Carlo simulations of a model beta-barrel globular protein having the six-stranded Greek key motif characteristic of real globular proteins were undertaken. The model protein possesses a typical beta-barrel amino acid sequence; however, all residues of a given type (e.g. hydrophobic residues) are identical. Even in the absence of site-specific interactions, starting from a high-temperature denatured state, these models undergo an all-or-none transition to a structurally unique six-stranded beta-barrel. These simulations suggest that the general rules of globular protein folding are rather robust in that the overall tertiary structure is determined by the general pattern of hydrophobic, hydrophilic, and turn-type residues, with site-specific interactions mainly involved in structural fine tuning of a given topology. Finally, these studies suggest that loops may play an important role in producing a unique native state. Depending on the stability of the native conformation of the long loop in the Greek key, the conformational transition can be described by a two-state, three-state, or even larger number of multiple equilibrium states model.

Models, Molecular↗
PubMed1989

Computer simulations of globular protein folding and tertiary structure.

In summary, although a large number of disparate techniques have been applied to predict the tertiary structure of globular proteins from their amino acid sequence, the solution is not yet at hand. Methodologies for predicting the conformation of constrained, small protein fragments appear to be successful. As the size of the system increases, the level of detail of the treatment decreases; approaches that employ very detailed potentials appear to be limited to about 30-40 residues. Although this is a major advance, methods that reduce the effective number of degrees of freedom are clearly required. Lattice representations coupled to highly efficient Monte Carlo procedures appear to be one such approach. Thus, although a number of theoretical advances in the computer simulation of globular protein structure have been made, much work remains to be done before the globular protein folding problem is solved.

Computer Simulation↗
PubMed1988

Monte Carlo simulations of the folding of beta-barrel globular proteins.

With the use of dynamic Monte Carlo simulations, the necessary conditions for the collapse from a random-coil denatured state to a structurally unique four-member beta-barrel native state of a model globular protein have been investigated. These systems are free to roam through all of configuration space--both native and nonnative interactions are allowed. The relative importance of hydrophobic and hydrophilic interactions and the presence or absence of statistical bend-forming regions for the formation of a unique native state are examined, and the conditions necessary for a denatured-to-native (and vice versa) conformational transition that is thermodynamically all-or-none and which always results in collapse to the same, four-member beta-barrel are explored. These conditions are found to be a general pattern of hydrophobic/hydrophilic residues that allows the native state to differentiate the interior from the exterior of the protein and the presence of regions that are, at the very least, neutral toward turn formation. The former set of interactions seems to define the mean length of the beta-stretch, and the latter set serves to lock the native state into the lowest free energy state, the native conformation. These folding simulations strongly suggest that the general rules of protein folding are rather robust and that site-specific tertiary interactions are only involved in structural fine tuning. The conditions required for the formation of a structurally unique native state from a manifold of collapsed conformations that are originally quite close in energy is highly suggestive of a mechanism of protein evolution by means of random mutations. The implications of these folding studies for such a mechanism are qualitatively explored.

Chemical Phenomena↗
PubMed1986

Monte Carlo simulations on an equilibrium globular protein folding model.

Monte Carlo simulations were performed on a diamond lattice, globular protein model in which the trans conformational state is energetically favored over the gauche states (thereby perhaps favoring a beta-sheet secondary structure) and in which nonspecific nonbonded nearest-neighbor attractive interactions are allowed. If the attractive interactions are sufficiently weak that the molecule possesses a relatively high fraction of trans states in the denatured state, then on collapse, a beta-barrel tertiary structure, highly reminiscent of the "native" structure seen in beta-proteins, spontaneously forms. If, however, the attractive interactions are dominant, a coil-to-random globule collapse transition is observed. The roles of short-, medium-, and long-range interactions and topological constraints in determining the observed tertiary structure are addressed, and the implications and limitations of the simulations for the equilibrium folding process in renal globular proteins are explored.

Models, Structural↗