Evolutionary computer programming of protein folding and structure predictions

In order to understand the mechanism of protein folding and to assist the rational de-novo design of fast-folding, non-aggregating and stable artificial enzymes it is very helpful to be able to simulate protein folding reactions and to predict the structures of proteins and other biomacromolecules. Here, we use a method of computer programming called "evolutionary computer programming" in which a program evolves depending on the evolutionary pressure exerted on the program. In the case of the presented application of this method on a computer program for folding simulations, the evolutionary pressure exerted was towards faster finding deep minima in the energy landscape of protein folding. Already after 20 evolution steps, the evolved program was able to find deep minima in the energy landscape more than 10 times faster than the original program prior to the evolution process.     

Tracing Primordial Protein Evolution through Structurally Guided Stepwise Segment Elongation

The understanding of how primordial proteins emerged has been a fundamental and longstanding issue in biology and biochemistry. For a better understanding of primordial protein evolution, we synthesized an artificial protein on the basis of an evolutionary hypothesis, segment-based elongation starting from an autonomously foldable short peptide. A 10-residue protein, chignolin, the smallest foldable polypeptide ever reported, was used as a structural support to facilitate higher structural organization and gain-of-function in the development of an artificial protein. Repetitive cycles of segment elongation and subsequent phage display selection successfully produced a 25-residue protein, termed AF.2A1, with nanomolar affinity against the Fc region of immunoglobulin G. AF.2A1 shows exquisite molecular recognition ability such that it can distinguish conformational differences of the same molecule. The structure determined by NMR measurements demonstrated that AF.2A1 forms a globular protein-like conformation with the chignolin-derived β-hairpin and a tryptophan-mediated hydrophobic core. Using sequence analysis and a mutation study, we discovered that the structural organization and gain-of-function emerged from the vicinity of the chignolin segment, revealing that the structural support served as the core in both structural and functional development. Here, we propose an evolutionary model for primordial proteins in which a foldable segment serves as the evolving core to facilitate structural and functional evolution. This study provides insights into primordial protein evolution and also presents a novel methodology for designing small sized proteins useful for industrial and pharmaceutical applications.     

Meta-structure correlation in protein space unveils different selection rules for folded and intrinsically disordered proteins

The number of existing protein sequences spans a very small fraction of sequence space. Natural proteins have overcome a strong negative selective pressure to avoid the formation of insoluble aggregates. Stably folded globular proteins and intrinsically disordered proteins (IDPs) use alternative solutions to the aggregation problem. While in globular proteins folding minimizes the access to aggregation prone regions, IDPs on average display large exposed contact areas. Here, we introduce the concept of average meta-structure correlation maps to analyze sequence space. Using this novel conceptual view we show that representative ensembles of folded and ID proteins show distinct characteristics and respond differently to sequence randomization. By studying the way evolutionary constraints act on IDPs to disable a negative function (aggregation) we might gain insight into the mechanisms by which function-enabling information is encoded in IDPs.     

Protein folding and the order/disorder paradox

Most proteins encoded by the nuclear genome are synthesized in the cytoplasm and fold into precise 3D structures. During synthesis, the nascent polypeptide begins to fold as it traverses the large subunit of the ribosome and is assisted by molecular chaperones in attaining its precise folded/highly ordered state efficiently and in a biologically relevant timescale. Proteins that are misfolded are culled, re-routed, and marked by mechanisms such as ubiquitinylation for degradation ensuring strict quality control (QC). In addition to the highly ordered "globular" proteins, emerging evidence indicates that a large fraction of the proteome also comprises the so-called "Intrinsically Disordered Proteins" (IDPs). IDPs are proteins that lack rigid 3D structures and instead, exist as dynamic ensembles. The dynamic structures in the IDPs have many similarities with "normal" globular proteins such as the native (ordered), and non-native (molten globule, pre-molten globule, and coil-like) states seen during folding of "normal" globular proteins. However, unlike the case of the nascent globular proteins, IDPs evade being detected as "misfolded" and degraded by the cell's QC system. We refer to this paradox as the order/disorder paradox and postulate that the IDPs capitalize on their intrinsic promiscuity and ability to undergo disorder-to-order transitions upon binding to biological targets (coupled folding and binding) to escape the cell's surveillance machinery. Understanding the mechanism by which the IDPs evade the quality check has wide implications from protein folding to disease biology since the aggregation of misfolded proteins underlies several debilitating illnesses such as many neurodegenerative diseases and cancer.     

Are amyloid fibrils molecular spandrels?

Amyloid-β, the protein implicated in Alzheimer’s disease, along with a number of other proteins, has been shown to form amyloid fibrils. Fibril forming proteins share no common primary structure and have little known function. Furthermore, all proteins have the ability to form amyloid fibrils under certain conditions as the fibrillar structure lies at the global free energy minimum of proteins. This raises the question of the mechanism of the evolution of the amyloid fibril structure. Experimental evidence supports the hypothesis that the fibril structure is a by-product of the forces of protein folding and lies outside the bounds of evolutionary pressures.     

Repeat-protein folding: new insights into origins of cooperativity, stability, and topology

Although our understanding of globular protein folding continues to advance, the irregular tertiary structures and high cooperativity of globular proteins complicates energetic dissection. Recently, proteins with regular, repetitive tertiary structures have been identified that sidestep limitations imposed by globular protein architecture. Here we review recent studies of repeat-protein folding. These studies uniquely advance our understanding of both the energetics and kinetics of protein folding. Equilibrium studies provide detailed maps of local stabilities, access to energy landscapes, insights into cooperativity, determination of nearest-neighbor interaction parameters using statistical thermodynamics, relationships between consensus sequences and repeat-protein stability. Kinetic studies provide insight into the influence of short-range topology on folding rates, the degree to which folding proceeds by parallel (versus localized) pathways, and the factors that select among multiple potential pathways. The recent application of force spectroscopy to repeat-protein unfolding is providing a unique route to test and extend many of these findings.     

A Folding Insulator Defines Cryptic Domains in Tropomyosin

Multidomain proteins are the product of evolutionary selection for diversity of function through concatenation and repurposing of existing modular units of structures. In structures of proteins with multiple domains, components are often globular units stitched together with flexible linkers. Multidomain proteins often fold as multiple distinct order–disorder transitions. However, the relationship between structure and folding is not always straightforward. Tropomyosin binds to actin in muscle and cytoskeletal filaments. The structure is that of a continuous ɑ-helix lacking domain boundaries, but unfolding shows distinct transitions suggesting at least three possible domains do exist. To explore how domains might occur in a continuous structure, we used Lifson-Roig helix-coil models with sequence domains of varying helical nucleation propensities. Of these models, ones with a central folding insulator, separating folding of N- and C-terminal domains, are most consistent with experimental folding studies. The positions of domain boundaries are identified by hydrogen–deuterium exchange mass spectrometry. The presence of structurally cryptic folding domains in tropomyosin could relate to its evolution and explain the uneven distribution of deleterious mutations that lead to various cardiomyopathies.     

A structural dissection of amino acid substitutions in helical transmembrane proteins

The evolution of protein folds is under strong constraints from their surrounding environment. Although folding in water‐soluble proteins is driven primarily by hydrophobic forces, the nature of the forces that determine the folding and stability of transmembrane proteins are still not fully understood. Furthermore, the chemically heterogeneous lipid bilayer has a non‐uniform effect on protein structure. In this article, we attempt to get an insight into the nature of this effect by examining the impact of various types of local structure environment on amino acid substitution, based on alignments of high‐resolution structures of polytopic helical transmembrane proteins combined with sequences of close homologs. Compared to globular proteins, burying amino acid sidechains, especially hydrophilic ones, led to a lower increase in conservation in both the lipid‐water interface region and the hydrocarbon core region. This observation is due to surface residues in HTM proteins especially in the HC region being relatively highly conserved, suggesting higher evolutionary constraints from their specific interactions with the surrounding lipid molecules. Polar and small residues, particularly Pro and Gly, show a noticeable increase in conservation as they are positioned more towards the centre of the membrane, which is consistent with their recognized key roles in structural stability. In addition, the examination of hydrogen bonds in the membrane environment identified some exposed hydrophilic residues being better conserved when not hydrogen‐bonded to other residues, supporting the importance of lipid‐protein sidechain interactions. The conclusions presented in this study highlight the distinct features of substitution matrices that take into account the membrane environment, and their potential role in improving sequence‐structure alignments of transmembrane proteins. Proteins 2010; © 2010 Wiley‐Liss, Inc.     

Single amino acid substitutions producing instability of globular proteins. Calculation of their frequencies in the entire mutational spectra of the α- and β-subunits of human hemoglobin

Summary The frequencies of substitutions resulting in protein instability were calculated by a method estimating changes in stability produced by amino acid substitutions. The method takes into account the accessibility of an amino acid position to a solvent and changes in the specificity of amino acid interactions. When tested on human mutant hemoglobins, the method yielded predictions with a preciseness of 80%. The consideration of the evolutionary homologous proteins in the analysis allowed us to estimate the evolutionary constraints imposed on stability of their spatial structure. With these limitations, approximately 50% of amino acid substitutions in the entire mutational spectra of the - and -subunits of human hemoglobin were found to damage the spatial structure of the globular proteins.     

The evolutionary landscape of functional model proteins.

To study the distinct influences of structure and function on evolution, we propose a minimalist model for proteins with binding pockets, called functional model proteins, based on a shifted-HP model on a two-dimensional square lattice. These model proteins are not maximally compact and contain an empty lattice site surrounded by at least three nearest neighbors, thus providing a binding pocket. Functional model proteins possess a unique native state, cooperative folding and tolerance to mutation. Due to the explicit functionality in these models (by design), we have been able to explore their fitness or evolutionary landscapes, as characterized by the size and distribution of homologous families and by the complexity of the inter-relatedness of the functional model proteins. Mindful that these minimalist models are highly idealized and two-dimensional, functional model proteins should nevertheless provide a useful means for exploring the constraints of maintaining structure and function on the evolution of proteins.     

Spatial constraints and group behaviour in globular proteins

A comparative analysis has been attempted on the spatial placement of amino acid residues derived from radial, ellipsoidal and exposure arrangements. The group behaviour of residues and their restraining influence in protein folding have been brought out. A study is also made on the geometry of proteins, the exposure arrangement of residues and the spatial distribution of the physical properties of the residues in globular proteins. It has been shown that the group constraints along with the information on the shape of the globular proteins would be highly useful in assigning the spatial and exposure arrangements of residues in globular proteins.     

Expression Differentiation Is Constrained to Low-Expression Proteins over Ecological Timescales

Protein expression level is one of the strongest predictors of protein sequence evolutionary rate, with high-expression protein sequences evolving at slower rates than low-expression protein sequences largely because of constraints on protein folding and function. Expression evolutionary rates also have been shown to be negatively correlated with expression level across human and mouse orthologs over relatively long divergence times (i.e., ∼100 million years). Long-term evolutionary patterns, however, often cannot be extrapolated to microevolutionary processes (and vice versa), and whether this relationship holds for traits evolving under directional selection within a single species over ecological timescales (i.e., <5000 years) is unknown and not necessarily expected. Expression is a metabolically costly process, and the expression level of a particular protein is predicted to be a tradeoff between the benefit of its function and the costs of its expression. Selection should drive the expression level of all proteins close to values that maximize fitness, particularly for high-expression proteins because of the increased energetic cost of production. Therefore, stabilizing selection may reduce the amount of standing expression variation for high-expression proteins, and in combination with physiological constraints that may place an upper bound on the range of beneficial expression variation, these constraints could severely limit the availability of beneficial expression variants. To determine whether rapid-expression evolution was restricted to low-expression proteins owing to these constraints on highly expressed proteins over ecological timescales, we compared venom protein expression levels across mainland and island populations for three species of pit vipers. We detected significant differentiation in protein expression levels in two of the three species and found that rapid-expression differentiation was restricted to low-expression proteins. Our results suggest that various constraints on high-expression proteins reduce the availability of beneficial expression variants relative to low-expression proteins, enabling low-expression proteins to evolve and potentially lead to more rapid adaptation.     

Inter-residue interactions in protein folding and stability

During the process of protein folding, the amino acid residues along the polypeptide chain interact with each other in a cooperative manner to form the stable native structure. The knowledge about inter-residue interactions in protein structures is very helpful to understand the mechanism of protein folding and stability. In this review, we introduce the classification of inter-residue interactions into short, medium and long range based on a simple geometric approach. The features of these interactions in different structural classes of globular and membrane proteins, and in various folds have been delineated. The development of contact potentials and the application of inter-residue contacts for predicting the structural class and secondary structures of globular proteins, solvent accessibility, fold recognition and ab initio tertiary structure prediction have been evaluated. Further, the relationship between inter-residue contacts and protein-folding rates has been highlighted. Moreover, the importance of inter-residue interactions in protein-folding kinetics and for understanding the stability of proteins has been discussed. In essence, the information gained from the studies on inter-residue interactions provides valuable insights for understanding protein folding and de novo protein design.     

Dynamic polyhedral models of globular proteins

We have devised several mechanical models of globular proteins by approximating them to various polyhedra (dodecahedron, truncated octahedron, icosahedron, truncated icosahedron). The models comprise hollow blocks linked together in a flexible chain. Between blocks there is a set of several reversible, weak magnetic interactions such that when the chain is agitated, it will fold into a stable polyhedral structure about the size of a hand. Folding may be followed in real time with a video camera. Key to the success of the folding process is the lightness of the chain. Several side chains may also be added to the blocks such that they come together to create a polyhedral core when the chain folds. The models have a number of similarities to globular proteins: each chain folds into a unique, but dynamic, three-dimensional structure; the instructions that determine this structure are built into the configuration of blocks; and it is difficult to predict this structure given the unfolded block configuration. Furthermore, the chains fold quickly, generally in less than a minute, several pathways are involved, and these pathways progress through elements of "native" structure. In particular, the models emphasize the importance of restricted conformational mobility in assisting the chain to fold, and also in eliminating undesirable interactions. Because of these similarities to globular proteins, we believe that the polyhedral models will, with continued development, be helpful in understanding the protein folding process, while at the same time acting as valuable educational visual aids. They might also inspire the construction of new types of microscopic, self-assembling devices.     

To Investigate Protein Evolution by Detecting Suppressed Epitope Structures

Material remains of ancestor nucleotides and proteins are largely unavailable, thus sequence comparison among homologous genes in present-day organisms forms the core of current knowledge of molecular evolution. Variation in protein three-dimensional structure is a basis for functional diversity. To study the evolution of three-dimensional structures in related proteins would significantly improve our understanding of protein evolution and function. A protein may contain ancestor conformations that have been allosterically suppressed by evolutionarily additive structures. Using monoclonal antibody probes to detect such conformation in proteins after removing the suppressor structure, our study demonstrates three-dimensional structure evidence for the evolutionary relationship between troponin I and troponin T, two subunits of the troponin complex in the Ca²⁺-regulatory system of striated muscle, and among their muscle type-specific isoforms. The experimental data show the feasibility of detecting evolutionarily suppressed history-telling structural states in proteins by removing conformational modulator segments added during evolution. In addition to identifying structural modifications that were critical to the emergence of diverged proteins, investigating this novel mode of evolution will help us to understand the origin and functional potential of protein structures.     

Evolution of Protein Targeting into “Complex” Plastids: The “Secretory Transport Hypothesis”

Abstract: In algae different types of plastids are known, which vary in pigment content and ultrastructure, providing an opportunity to study their evolutionary origin. One interesting feature is the number of envelope membranes surrounding the plastids. Red algae, green algae and glaucophytes have plastids with two membranes. They are thought to originate from a primary endocytobiosis event, a process in which a prokaryotic cyanobacterium was engulfed by a eukaryotic host cell and transformed into a plastid. Several other algal groups, like euglenophytes and heterokont algae (diatoms, brown algae, etc.), have plastids with three or four surrounding membranes, respectively, probably reflecting the evolution of these organisms by so‐called secondary endocytobiosis, which is the uptake of a eukaryotic alga by a eukaryotic host cell. A prerequisite for the successful establishment of primary or secondary endocytobiosis must be the development of suitable protein targeting machineries to allow the transport of nucleus‐encoded plastid proteins across the various plastid envelope membranes. Here, we discuss the possible evolution of such protein transport systems. We propose that the secretory system of the respective host cell might have been the essential tool to establish protein transport into primary as well as into secondary plastids.     

Amyloid Formation by Human Carboxypeptidase D Transthyretin-like Domain under Physiological Conditions

Protein aggregation is linked to a growing list of diseases, but it is also an intrinsic property of polypeptides, because the formation of functional globular proteins comes at the expense of an inherent aggregation propensity. Certain proteins can access aggregation-prone states from native-like conformations without the need to cross the energy barrier for unfolding. This is the case of transthyretin (TTR), a homotetrameric protein whose dissociation into its monomers initiates the aggregation cascade. Domains with structural homology to TTR exist in a number of proteins, including the M14B subfamily carboxypeptidases. We show here that the monomeric transthyretin-like domain of human carboxypeptidase D aggregates under close to physiological conditions into amyloid structures, with the population of folded but aggregation-prone states being controlled by the conformational stability of the domain. We thus confirm that the TTR fold keeps a generic residual aggregation propensity upon folding, resulting from the presence of preformed amyloidogenic β-strands in the native state. These structural elements should serve for functional/structural purposes, because they have not been purged out by evolution, but at the same time they put proteins like carboxypeptidase D at risk of aggregation in biological environments and thus can potentially lead to deposition diseases.     

Protein Folding and Evolution are Driven by the Maxwell Demon Activity of Proteins

In this paper we propose a theoretical model of protein folding and protein evolution in which a polypeptide (sequence/structure) is assumed to behave as a Maxwell Demon or Information Gathering and Using System (IGUS) that performs measurements aiming at the construction of the native structure. Our model proposes that a physical meaning to Shannon information (H) and Chaitin's algorithmic information (K) parameters can be both defined and referred from the IGUS standpoint. Our hypothesis accounts for the interdependence of protein folding and protein evolution through mutual influencing relationships mediated by the IGUS. In brief, IGUS activity in protein folding determines long term tendencies that emerge at the evolutionary time-scale.Thus, protein evolution is a consequence of measurements executed by proteins at the cellular level, where the IGUS imposes a tendency to attain a highly unique stable native form that promotes the updating of the information content. The folding kinetics observed is, thus, the outcome of an evolutionary process where the polypeptide-IGUS drives the evolution of its linear sequence. Finally, we describe protein evolution as an entropic process that tends to increase the content of mutual algorithmic information between the sequence and the structure. This model enables one: 1. To comprehend that full determination of the three-dimensional structure by the linear sequence is a tendency where satisfaction is only possible at thermodynamic equilibrium.2. To account for the observed randomness of the amino acid sequences. 3. To predict an alternation of periods of selection and neutral diffusion during protein evolutionary time.