Computational approach to the protein-folding problem

Physics of proteins

Globular proteins are a key component of the network of life. Over many decades much experimental data on proteins have been gathered, yet theoretical progress has been somewhat limited. We show that the results accumulated over the years inexorably lead to a unified framework for understanding proteins. The framework predicts the existence of a fixed menu of folds determined by geometry, clarifies the role of the amino acid sequence in selecting the native-state structure from this menu, and explains the propensity for amyloid formation. The experimental data and the new approach reveal an astonishing simplicity underlying the protein problem.

A crumpled surface having transverse attractive interactions as a simplified model with biological significance

Using extensive analogical simulations with square sheets of paper we investigate the influence of short-range transverse attractive interactions on the packing properties of a crumpled surface. These interactions are due to transverse connections or local bridges associated with a given number of binding sites localized on the two-dimensional surface and distributed in several patterns in the three-dimensional physical space. Geometrical relations and critical exponents describing the statistical properties of the crumpled surface are obtained as a function of the strength of the attractive interactions. Our model suggests how the presence of short-range interactions as, e.g. van der Waals forces can be important for the geometric plasticity of biological molecules, which in turn is important for biological function. The relevance of our results to the study of molecular conformation of proteins and membranes is discussed, and a comparison is also made between the behavior of the crumpled surface studied here and other important non-equilibrium fractal structures.

Unified perspective on proteins: a physics approach

We study a physical system which, while devoid of the complexity one usually associates with proteins, nevertheless displays a remarkable array of proteinlike properties. The constructive hypothesis that this striking resemblance is not accidental not only leads to a unified framework for understanding protein folding, amyloid formation, and protein interactions but also has implications for natural selection.

Recognition of coarse-grained protein tertiary structure

A model of the protein backbone is considered in which each residue is characterized by the location of its C(alpha) atom and one of a discrete set of conformal (phi, psi) states. We investigate the key differences between a description that offers a locally precise fit to known backbone structures and one that provides a globally accurate fit to protein structures. Using a statistical scoring scheme and threading, a protein's local best-fit conformation is highly recognizable, but its global structure cannot be directly determined from an amino acid sequence. The incorporation of information about the conformal states of neighboring residues along the chain allows one to accurately translate the local structure into a global structure. We present a two-step algorithm, which recognizes up to 95% of the tested protein native-state structures to within a 2.5 A root mean square deviation.

Solvation effects and driving forces for protein thermodynamic and kinetic cooperativity: how adequate is native-centric topological modeling?

What energetic and solvation effects underlie the remarkable two-state thermodynamics and folding/unfolding kinetics of small single-domain proteins? To address this question, we investigate the folding and unfolding of a hierarchy of continuum Langevin dynamics models of chymotrypsin inhibitor 2. We find that residue-based additive Gō-like contact energies, although native-centric, are by themselves insufficient for protein-like calorimetric two-state cooperativity. Further native biases by local conformational preferences are necessary for protein-like thermodynamics. Kinetically, however, even models with both contact and local native-centric energies do not produce simple two-state chevron plots. Thus a model protein's thermodynamic cooperativity is not sufficient for simple two-state kinetics. The models tested appear to have increasing internal friction with increasing native stability, leading to chevron rollovers that typify kinetics that are commonly referred to as non-two-state. The free energy profiles of these models are found to be sensitive to the choice of native contacts and the presumed spatial ranges of the contact interactions. Motivated by explicit-water considerations, we explore recent treatments of solvent granularity that incorporate desolvation free energy barriers into effective implicit-solvent intraprotein interactions. This additional feature reduces both folding and unfolding rates vis-à-vis that of the corresponding models without desolvation barriers, but the kinetics remain non-two-state. Taken together, our observations suggest that interaction mechanisms more intricate than simple Gō-like constructs and pairwise additive solvation-like contributions are needed to rationalize some of the most basic generic protein properties. Therefore, as experimental constraints on protein chain models, requiring a consistent account of protein-like thermodynamic and kinetic cooperativity can be more stringent and productive for some applications than simply requiring a model heteropolymer to fold to a target structure.

A Mutation in the gene-encoding bacteriophage T7 DNA polymerase that renders the phage temperature-sensitive

Gene 5 of bacteriophage T7 encodes a DNA polymerase essential for phage replication. A single point mutation in gene 5 confers temperature sensitivity for phage growth. The mutation results in an alanine to valine substitution at residue 73 in the exonuclease domain. Upon infection of Escherichia coli by the temperature-sensitive phage at 42 degrees C, there is no detectable T7 DNA synthesis in vivo. DNA polymerase activity in these phage-infected cell extracts is undetectable at assay temperatures of 30 degrees C or 42 degrees C. Upon infection at 30 degrees C, both DNA synthesis in vivo and DNA polymerase activity in cell extracts assayed at 30 degrees C or 42 degrees C approach levels observed using wild-type T7 phage. The amount of soluble gene 5 protein produced at 42 degrees C is comparable to that produced at 30 degrees C, indicating that the temperature-sensitive phenotype is not due to reduced expression, stability, or solubility. Thus the polymerase induced at elevated temperatures by the temperature-sensitive phage is functionally inactive. Consistent with this observation, biochemical properties and heat inactivation profiles of the genetically altered enzyme over-produced at 30 degrees C closely resemble that of wild-type T7 DNA polymerase. It is likely that the polymerase produced at elevated temperatures is a misfolded intermediate in its folding pathway.

Role of native-state topology in the stabilization of intracellular antibodies

The role played by the geometric position of each amino acid in the folding process of the immunoglobulin (Ig) variable domain is identified and measured through molecular dynamics simulations of models based on the topology of its native state. This measure allows identifying the parts of the protein that, for geometrical reasons, when mutated, would result in relevant protein stability changes. Simulations were performed without considering the covalent disulfide bond present in most of the Ig domains. The results are in good agreement with site-directed mutagenesis experiments on the folding of intracellular antibodies in which the disulfide bond does not form. We also found agreement with data on amino acid conservation in the Ig variable domain sequences. This indicates a new way for a rational approach to the design of intracellular antibodies more resistant to the suppression of the disulfide bond that occurs in the cytoplasm.

Modeling two-state cooperativity in protein folding

A protein model with the pairwise interaction energies varying as the local environment changes, i.e., including some kind of collective effect between the contacts, is proposed. Lattice Monte Carlo simulations on the thermodynamical characteristics and free energy profile show a well-defined two-state behavior and cooperativity of folding for such a model. As a comparison, related simulations for the usual G&o model, where the interaction energies are independent of the local conformation, are also made. Our results indicate that the evolution of interactions during the folding process plays an important role in the two-state cooperativity in protein folding.