Structure of infectious prions: stabilization by domain swapping

Human FoxP Transcription Factors as Tractable Models of the Evolution and Functional Outcomes of Three-Dimensional Domain Swapping

The association of two or more proteins to adopt a quaternary complex is one of the most widespread mechanisms by which protein function is modulated. In this scenario, three-dimensional domain swapping (3D-DS) constitutes one plausible pathway for the evolution of protein oligomerization that exploits readily available intramolecular contacts to be established in an intermolecular fashion. However, analysis of the oligomerization kinetics and thermodynamics of most extant 3D-DS proteins shows its dependence on protein unfolding, obscuring the elucidation of the emergence of 3D-DS during evolution, its occurrence under physiological conditions, and its biological relevance. Here, we describe the human FoxP subfamily of transcription factors as a feasible model to study the evolution of 3D-DS, due to their significantly faster dissociation and dimerization kinetics and lower dissociation constants in comparison to most 3D-DS models. Through the biophysical and functional characterization of FoxP proteins, relevant structural aspects highlighting the evolutionary adaptations of these proteins to enable efficient 3D-DS have been ascertained. Most biophysical studies on FoxP suggest that the dynamics of the polypeptide chain are crucial to decrease the energy barrier of 3D-DS, enabling its fast oligomerization under physiological conditions. Moreover, comparison of biophysical parameters between human FoxP proteins in the context of their minute sequence differences suggests differential evolutionary strategies to favor homoassociation and presages the possibility of heteroassociations, with direct impacts in their gene regulation function.

Assembly of platforms for signal transduction in the new era: dimerization, helical filament assembly, and beyond

Supramolecular organizing center (SMOC)-mediated signal transduction is an emerging concept in the field of signal transduction that is ushering in a new era. The formation of location-specific, higher-order SMOCs is particularly important for cell death and innate immune signaling processes. Several protein interaction domains, including the death domain (DD) superfamily and the CIDE domain, are representative mediators of SMOC assembly in cell death and innate immune signaling pathways. DD superfamily- and CIDE domain-containing proteins form SMOCs that activate various caspases and provide signaling scaffold platforms. These assemblies can lead to signal transduction and amplification during signaling events. In this review, we summarize recent findings on the molecular basis of DD superfamily- and CIDE domain-mediated SMOC formation.

Improved understanding of large molecular signaling complexes that form during innate (nonspecific) immune responses could help develop treatments for multiple diseases including cancer. Correct cell signaling requires precise protein interactions and binding, which are mediated by specific sites on the surface of the protein molecules involved. Innate immune responses and cell death mechanisms rely on such protein interactions, and defects can cause signaling abnormalities and trigger disease. Hyun Ho Park and co-workers at Chung-Ang University in Seoul, South Korea, reviewed recent insights into the presence of supramolecular organizing centers (SMOCs), localized complexes of signaling proteins that form during immune responses. The researchers highlight existing understanding of SMOC assembly processes. A better understanding of SMOCs will help to explain enzyme activation, signal amplification and cell signaling control mechanisms.

Structural mechanisms of oligomer and amyloid fibril formation by the prion protein

Misfolding and aggregation of the prion protein is responsible for multiple neurodegenerative diseases. Works from several laboratories on folding of both the WT and multiple pathogenic mutant variants of the prion protein have identified several structurally dissimilar intermediates, which might be potential precursors to misfolding and aggregation. The misfolded aggregates themselves are morphologically distinct, critically dependent on the solution conditions under which they are prepared, but always β-sheet rich. Despite the lack of an atomic resolution structure of the infectious pathogenic agent in prion diseases, several low resolution models have identified the β-sheet rich core of the aggregates formed in vitro, to lie in the α2-α3 subdomain of the prion protein, albeit with local stabilities that vary with the type of aggregate. This feature article describes recent advances in the investigation of in vitro prion protein aggregation using multiple spectroscopic probes, with particular focus on (1) identifying aggregation-prone conformations of the monomeric protein, (2) conditions which trigger misfolding and oligomerization, (3) the mechanism of misfolding and aggregation, and (4) the structure of the misfolded intermediates and final aggregates.

Comparing the energy landscapes for native folding and aggregation of PrP

Protein sequences are evolved to encode generally one folded structure, out of a nearly infinite array of possible folds. Underlying this code is a funneled free energy landscape that guides folding to the native conformation. Protein misfolding and aggregation are also a manifestation of free-energy landscapes. The detailed mechanisms of these processes are poorly understood, but often involve rare, transient species and a variety of different pathways. The inherent complexity of misfolding has hampered efforts to measure aggregation pathways and the underlying energy landscape, especially using traditional methods where ensemble averaging obscures important rare and transient events. We recently studied the misfolding and aggregation of prion protein by examining 2 monomers tethered in close proximity as a dimer, showing how the steps leading to the formation of a stable aggregated state can be resolved in the single-molecule limit and the underlying energy landscape thereby reconstructed. This approach allows a more quantitative comparison of native folding versus misfolding, including fundamental differences in the dynamics for misfolding. By identifying key steps and interactions leading to misfolding, it should help to identify potential drug targets. Here we describe the importance of characterizing free-energy landscapes for aggregation and the challenges involved in doing so, and we discuss how single-molecule studies can help test proposed structural models for PrP aggregates.

Structural Biology of PrP Prions

Prion diseases are characterized by the deposition of amyloids, misfolded conformers of the prion protein. The misfolded conformation is self-replicating, by a mechanism solely enciphered in the conformation of the protein. Because of low solubility and heterogeneous aggregate sizes, the detailed atomic structure of the infectious isoform is still unknown. Progress has, however, been made, and has allowed insights into the structural and disease-related mechanisms of prions. Many structural models have been proposed, and a number of them support a consensus trimeric β-helical model, significantly more complex than simple amyloid models. There is evidence that such complexity may be a necessary property of prion structure. Knowledge of the structure of prions will provide a greater understanding of the protein isoform conversion mechanism, and could eventually lead to rationally designed intervention strategies.

Le TN, Aulić S, Bistaffa E, Campagnani I, Virgilio T, Indaco A, Palamara L, Andréoletti O, Tagliavini F, Legname G. Synthetic prions with novel strain-specified properties

Prions are infectious proteins that possess multiple self-propagating structures. The information for strains and structural specific barriers appears to be contained exclusively in the folding of the pathological isoform, PrP^Sc. Many recent studies determined that de novo prion strains could be generated in vitro from the structural conversion of recombinant (rec) prion protein (PrP) into amyloidal structures. Our aim was to elucidate the conformational diversity of pathological recPrP amyloids and their biological activities, as well as to gain novel insights in characterizing molecular events involved in mammalian prion conversion and propagation. To this end we generated infectious materials that possess different conformational structures. Our methodology for the prion conversion of recPrP required only purified rec full-length mouse (Mo) PrP and common chemicals. Neither infected brain extracts nor amplified PrP^Sc were used. Following two different in vitro protocols recMoPrP converted to amyloid fibrils without any seeding factor. Mouse hypothalamic GT1 and neuroblastoma N2a cell lines were infected with these amyloid preparations as fast screening methodology to characterize the infectious materials. Remarkably, a large number of amyloid preparations were able to induce the conformational change of endogenous PrP^C to harbor several distinctive proteinase-resistant PrP forms. One such preparation was characterized in vivo habouring a synthetic prion with novel strain specified neuropathological and biochemical properties.

Prions are infectious proteins capable of acquiring multiple self-propagating structures. The information for strains and structural specific barriers appears to be contained exclusively in the folding of the pathological isoform, designated as PrP^Sc. During propagation, disease-associated conformer PrP^Sc coerces the physiological form, denoted as PrP^C, to adopt the pathological isoform conformation. We describe here the generation of an array of infectious materials with different structural, morphological, biochemical and cell biological characteristics. After producing purified recombinant prion protein of the wild-type mouse full-length sequence in Escherichia coli, we polymerized the protein into various amyloid fibril conformations based on different amyloid preparations. We also applied a build-in methodology for screening amyloid preparations and generate infectious materials using an amyloid-infected cell culture assay. Some of the amyloid fibrils preparations were able to efficiently amplify in PMCA (Protein Misfolding Cyclic Amplification), and to induce endogenous PrP^C to convert into PrP^Sc in both murine hypothalamic GT1 and neuroblastoma N2a cell lines. One such protocol lead to the generation of a novel synthetic prion strain in mice.

The structure of the infectious prion protein: experimental data and molecular models

The structures of the infectious prion protein, PrP(Sc), and that of its proteolytically truncated variant, PrP 27-30, have evaded experimental determination due to their insolubility and propensity to aggregate. Molecular modeling has been used to fill this void and to predict their structures, but various modeling approaches have produced significantly different models. The disagreement between the different modeling solutions indicates the limitations of this method. Over the years, in absence of a three-dimensional (3D) structure, a variety of experimental techniques have been used to gain insights into the structure of this biologically, medically, and agriculturally important isoform. Here, we present an overview of experimental results that were published in recent years, and which provided new insights into the molecular architecture of PrP(Sc) and PrP 27-30. Furthermore, we evaluate all published models in light of these recent, experimental data, and come to the conclusion that none of the models can accommodate all of the experimental constraints. Moreover, this conclusion constitutes an open invitation for renewed efforts to model the structure of PrP(Sc).

Dimer domain swapping versus monomer folding in apo-myoglobin studied by molecular simulations

Using a coarse-grained symmetrized Go model, we performed a series of folding simulations of two apo-myoglobin molecules restrained at a high density, addressing competition of formation of a domain-swapped dimer with folding to two monomer structures.

The diversity and relationship of prion protein self-replicating states

It has become evident that the prion protein (PrP) can form a diverse range of self-replicating structures in addition to bona fide PrP(Sc) or strain-specific PrP(Sc) variants. Some self-replicating states can be only produced in vitro, whereas others can be formed in vivo and in vitro. While transmissible, not all states that replicate in vivo are truly pathogenic. Some of them can replicate silently without causing symptoms or clinical diseases. In the current article we discuss the data on PK-digestion patterns of different self-replicating PrP states in connection with other structural data available to date and assess possible relationships between different self-replicating states. Even though different self-replicating PrP states appear to have significantly different global folding patterns, it seems that the C-terminal region exhibits a cross-β-sheet structure in all self-replicating states, as this region acquires the proteolytically most stable conformation. We also discuss the possibility of the transformation of self-replicating states and triggering of PrP(Sc) formation within the frame of the deformed templating model. The spread of silent self-replicating states is of a particular concern because they can lead to transmissible prion disease. Moreover, examples on how different replication requirements favor different states are discussed. This knowledge can help in designing conditions for selective amplification of a particular PrP state in vitro.

Functional repertoire, molecular pathways and diseases associated with 3D domain swapping in the human proteome

Background

3D domain swapping is a novel structural phenomenon observed in diverse set of protein structures in oligomeric conformations. A distinct structural feature, where structural segments in a protein dimer or higher oligomer were shared between two or more chains of a protein structure, characterizes 3D domain swapping. 3D domain swapping was observed as a key mediator of numerous functional mechanisms and play pathogenic role in various diseases including conformational diseases like amyloidosis, Alzheimer's disease, Parkinson's disease and prion diseases. We report the first study with a focus on identifying functional classes, pathways and diseases mediated by 3D domain swapping in the human proteome.

Methods

We used a panel of four enrichment tools with two different ontologies and two annotations database to derive biological and clinical relevant information associated with 3D domain swapping. Protein domain enrichment analysis followed by Gene Ontology (GO) term enrichment analysis revealed the functional repertoire of proteins involved in swapping. Pathway analysis using KEGG annotations revealed diverse pathway associations of human proteins involved in 3D domain swapping. Disease Ontology was used to find statistically significant associations with proteins in swapped conformation and various disease categories (P-value < 0.05).

Results

We report meta-analysis results of a literature-curated dataset of human gene products involved in 3D domain swapping and discuss new insights about the functional repertoire, pathway associations and disease implications of proteins involved in 3D domain swapping.

Conclusions

Our integrated bioinformatics pipeline comprising of four different enrichment tools, two ontologies and two annotations revealed new insights into the functional and disease correlations with 3D domain swapping. GO term enrichment were used to infer terms associated with three different GO categories. Protein domain enrichment was used to identify conserved domains enriched in swapped proteins. Pathway enrichment analysis using KEGG annotations revealed that proteins with swapped conformations are present in all six classes of KEGG BRITE hierarchy and significantly enriched KEGG pathways were observed in five classes. Five major classes of disease were found to be associated with 3D domain swapping using functional disease ontology based enrichment analysis. Five classes of human diseases: cancer, diseases of the respiratory or pulmonary system, degenerative diseases of the central nervous system, vascular disease and encephalitis were found to be significant. In conclusion, our study shows that bioinformatics based analytical approaches using curated data can enhance the understanding of functional and disease implications of 3D domain swapping.

Implications of 3D domain swapping for protein folding, misfolding and function

Three-dimensional domain swapping is the process by which two identical protein chains exchange a part of their structure to form an intertwined dimer or higher-order oligomer. The phenomenon has been observed in the crystal structures of a range of different proteins. In this chapter we review the experiments that have been performed in order to understand the sequence and structural determinants of domain-swapping and these show how the general principles obtained can be used to engineer proteins to domain swap. We discuss the role of domain swapping in regulating protein function and as one possible mechanism of protein misfolding that can lead to aggregation and disease. We also review a number of interesting pathways of macromolecular assembly involving β-strand insertion or complementation that are related to the domain-swapping phenomenon.

Amyloid-like fibrils from a domain-swapping protein feature a parallel, in-register conformation without native-like interactions

The formation of amyloid-like fibrils is characteristic of various diseases, but the underlying mechanism and the factors that determine whether, when, and how proteins form amyloid, remain uncertain. Certain mechanisms have been proposed based on the three-dimensional or runaway domain swapping, inspired by the fact that some proteins show an apparent correlation between the ability to form domain-swapped dimers and a tendency to form fibrillar aggregates. Intramolecular β-sheet contacts present in the monomeric state could constitute intermolecular β-sheets in the dimeric and fibrillar states. One example is an amyloid-forming mutant of the immunoglobulin binding domain B1 of streptococcal protein G, which in its native conformation consists of a four-stranded β-sheet and one α-helix. Under native conditions this mutant adopts a domain-swapped dimer, and it also forms amyloid-like fibrils, seemingly in correlation to its domain-swapping ability. We employ magic angle spinning solid-state NMR and other methods to examine key structural features of these fibrils. Our results reveal a highly rigid fibril structure that lacks mobile domains and indicate a parallel in-register β-sheet structure and a general loss of native conformation within the mature fibrils. This observation contrasts with predictions that native structure, and in particular intermolecular β-strand interactions seen in the dimeric state, may be preserved in "domain-swapping" fibrils. We discuss these observations in light of recent work on related amyloid-forming proteins that have been argued to follow similar mechanisms and how this may have implications for the role of domain-swapping propensities for amyloid formation.

Optical trapping with high forces reveals unexpected behaviors of prion fibrils

Amyloid fibrils are important in diverse cellular functions, feature in many human diseases and have potential applications in nanotechnology. Here we describe methods that combine optical trapping and fluorescent imaging to characterize the forces that govern the integrity of amyloid fibrils formed by a yeast prion protein. A crucial advance was to use the self-templating properties of amyloidogenic proteins to tether prion fibrils, enabling their manipulation in the optical trap. At normal pulling forces the fibrils were impervious to disruption. At much higher forces (up to 250 pN), discontinuities occurred in force-extension traces before fibril rupture. Experiments with selective amyloid-disrupting agents and mutations demonstrated that such discontinuities were caused by the unfolding of individual subdomains. Thus, our results reveal unusually strong noncovalent intermolecular contacts that maintain fibril integrity even when individual monomers partially unfold and extend fibril length.

Insights into prion protein function from atomistic simulations

Computer simulations are a powerful tool for studies of biological systems. They have often been used to study prion protein (PrP), a protein responsible for neurodegenerative diseases, which include "mad cow disease" in cattle and Creutzfeldt-Jacob disease in humans. An important aspect of the prion protein is its interaction with copper ion, which is thought to be relevant for PrP's yet undetermined function and also potentially play a role in prion diseases. for studies of copper attachment to the prion protein, computer simulations have often been used to complement experimental data and to obtain binding structures of Cu-PrP complexes. This paper summarizes the results of recent ab initio calculations of copper-prion protein interactions focusing on the recently discovered concentration-dependent binding modes in the octarepeat region of this protein. In addition to determining the binding structures, computer simulations were also used to make predictions about PrP's function and the role of copper in prion diseases. The results demonstrate the predictive power and applicability of ab initio simulations for studies of metal-biomolecular complexes.

A rapid coarse residue-based computational method for x-ray solution scattering characterization of protein folds and multiple conformational states of large protein complexes

We present a coarse residue-based computational method to rapidly compute the solution scattering profile from a protein with dynamical fluctuations. The method is built upon a coarse-grained (CG) representation of the protein. This CG representation takes advantage of the intrinsic low-resolution and CG nature of solution scattering data. It allows rapid scattering determination from a large number of conformations that can be extracted from CG simulations to obtain scattering characterization of protein conformations. The method includes several important elements, effective residue structure factors derived from the Protein Data Bank, explicit treatment of water molecules in the hydration layer at the surface of the protein, and an ensemble average of scattering from a variety of appropriate conformations to account for macromolecular flexibility. This simplified method is calibrated and illustrated to accurately reproduce the experimental scattering curve of Hen egg white lysozyme. We then illustrated the applications of this CG method by computing the solution scattering patterns of several representative protein folds and multiple conformational states. The results suggest that solution scattering data, when combined with the reliable computational method that we developed, show great potential for a better structural description of multidomain complexes in different functional states, and for recognizing structural folds when sequence similarity to a protein of known structure is low.

Elastic energy driven polymerization

We present a molecular system where polymerization is controlled externally by tuning the elastic energy of the monomers. The elastic energy, provided by a DNA molecular spring, destabilizes the monomer state through a process analogous to domain swapping. This energy can be large (of approximately 10 kT) and thus drive polymerization at relatively low monomer concentrations. The monomer-dimer equilibrium provides a measurement of the elastic energy of the monomer, which in this construction appears limited by kink formation in the DNA molecular spring, in accord with previous theoretical and experimental investigations of the elasticity of sharply bent DNA.

Structure of mammalian prions

The mammalian prion PrPSc causes fatal neurodegenerative ailments in humans and farm animals. The molecular mechanisms of its propagation and transmissibility between species will not be understood until sufficient knowledge of its structure is gathered. This task has been hampered by the insoluble nature of PrPSc, which renders classic techniques, such as nuclear magnetic resonance or x-ray crystallography, unusable. However, a number of alternative approaches, such as limited proteolysis, electron microscopy and Fourier transform infrared spectroscopy, have yielded valuable low-resolution structural information. Pieced together, these data present PrPSc as a stackable molecule whose core is probably a β-solenoid formed by two or three rungs of short β-strands interspersed with numerous loops and turns. A reasonable understanding of its architecture might soon be achieved.

Analysis of the sequence and structural features of the left-handed beta-helical fold

The left-handed parallel beta-helix (LbetaH) is a structurally repetitive, highly regular, and symmetrical fold formed by coiling of elongated beta-sheets into helical "rungs." This canonical fold has recently received interest as a possible solution to the fibril structure of amyloid and as a building block of self-assembled nanotubular structures. In light of this interest, we aimed to understand the structural requirements of the LbetaH fold. We first sought to determine the sequence characteristics of the repeats by analyzing known structures to identify positional preferences of specific residues types. We then used molecular dynamics simulations to demonstrate the stabilizing effect of successive rungs and the hydrophobic core of the LbetaH. We show that a two-rung structure is the minimally stable LbetaH structure. In addition, we defined the structure-based sequence preference of the LbetaH and undertook a genome-wide sequence search to determine the prevalence of this unique protein fold. This profile-based LbetaH search algorithm predicted a large fraction of LbetaH proteins from microbial origins. However, the relative number of predicted LbetaH proteins per specie was approximately equal across the genomes from prokaryotes to eukaryotes.

Essential role of hydration in aggregation of misfolded prion proteins: quantification by molecular theory of solvation

A statistical-mechanical, three-dimensional molecular theory of solvation (also know as 3D-RISM) and molecular mechanics were used to study the thermodynamics of aggregation of misfolded prion proteins, based on the theoretical molecular models proposed so far. These include the beta-helical prion trimer (BPT) model of Govaerts et al. (2004), the domain-swapped trimeric prion (DSTP) model of Yang et al. (2005), and the model built after the spiral model of DeMarco and Daggett (2004). It is shown that the solvation contribution to the association free energy can overcome the gain in the internal energy upon association of the proteins. The solvation entropic contribution is as important as the energetic term in the total association free energy. Our calculations show that the spiral-like model is thermodynamically less stable, compared to the DSTP and BPT models. Among the latter two models, the DSTP model is more favorable to association. Quantitative assessment of the solvation effects on the association thermodynamics of prion proteins is provided, and explicitly shows that the solvation contribution is a driving force of the association, in particular, for the existing theoretical models of misfolded prion proteins.

Structural characterization of a neurotoxic threonine-rich peptide corresponding to the human prion protein alpha 2-helical 180-195 segment, and comparison with full-length alpha 2-helix-derived peptides

The 173-195 segment corresponding to the helix 2 of the globular PrP domain is a good candidate to be one of the several 'spots' of intrinsic structural flexibility, which might induce local destabilization and concur to protein transformation, leading to aggregation-prone conformations. Here, we report CD and NMR studies on the alpha2-helix-derived peptide of maximal length (hPrP[180-195]) that is able to exhibit a regular structure different from the prevalently random arrangement of other alpha2-helix-derived peptides. This peptide, which has previously been shown to be affected by buffer composition via the ion charge density dependence typical of Hofmeister effects, corresponds to the C-terminal sequence of the PrP(C) full-length alpha2-helix and includes the highly conserved threonine-rich 188-195 segment. At neutral pH, its conformation is dominated by beta-type contributions, which only very strong environmental modifications are able to modify. On TFE addition, an increase of alpha-helical content can be observed, but a fully helical conformation is only obtained in neat TFE. However, linking of the 173-179 segment, as occurring in wild-type and mutant peptides corresponding to the full-length alpha2-helix, perturbs these intrinsic structural propensities in a manner that depends on whether the environment is water or TFE. Overall, these results confirm that the 180-195 parental region in hPrP(C) makes a strong contribution to the chameleon conformational behavior of the segment corresponding to the full-length alpha2-helix, and could play a role in determining structural rearrangements of the entire globular domain.

Left handed beta helix models for mammalian prion fibrils

We propose models for in vitro grown mammalian prion protein fibrils based upon left handed beta helices formed both from the N-terminal and C-terminal regions of the proteinase resistant infectious prion core. The C-terminal threading onto a beta-helical structure is almost uniquely determined by fixing the cysteine disulfide bond on a helix corner. In comparison to known left handed helical peptides, the resulting model structures have similar stability attributes including relatively low root mean square deviations in all atom molecular dynamics, substantial side-chain-to-side-chain hydrogen bonding, good volume packing fraction, and low hydrophilic/hydrophobic frustration. For the N-terminus, we propose a new threading of slightly more than two turns, which improves upon the above characteristics relative to existing three turn beta-helical models. The N-terminal and C-terminal beta helices can be assembled into eight candidate models for the fibril repeat units, held together by large hinge (order 30 residues) domain swapping, with three amenable to fibril promoting domain swapping via a small (five residue) hinge on the N-terminal side. Small concentrations of the metastable C-terminal beta helix in vivo might play a significant role in templating the infectious conformation and in enhancing conversion kinetics for inherited forms of the disease and explain resistance (for canines) involving hypothesized coupling to the methionine 129 sulfur known to play a role in human disease.

Dynamics of the acetylcholinesterase tetramer

Acetylcholinesterase rapidly hydrolyzes the neurotransmitter acetylcholine in cholinergic synapses, including the neuromuscular junction. The tetramer is the most important functional form of the enzyme. Two low-resolution crystal structures have been solved. One is compact with two of its four peripheral anionic sites (PAS) sterically blocked by complementary subunits. The other is a loose tetramer with all four subunits accessible to solvent. These structures lacked the C-terminal amphipathic t-peptide (WAT domain) that interacts with the proline-rich attachment domain (PRAD). A complete tetramer model (AChEt) was built based on the structure of the PRAD/WAT complex and the compact tetramer. Normal mode analysis suggested that AChEt could exist in several conformations with subunits fluctuating relative to one another. Here, a multiscale simulation involving all-atom molecular dynamics and C alpha-based coarse-grained Brownian dynamics simulations was carried out to investigate the large-scale intersubunit dynamics in AChEt. We sampled the ns-mus timescale motions and found that the tetramer indeed constitutes a dynamic assembly of monomers. The intersubunit fluctuation is correlated with the occlusion of the PAS. Such motions of the subunits "gate" ligand-protein association. The gates are open more than 80% of the time on average, which suggests a small reduction in ligand-protein binding. Despite the limitations in the starting model and approximations inherent in coarse graining, these results are consistent with experiments which suggest that binding of a substrate to the PAS is only somewhat hindered by the association of the subunits.

Molecular dynamics with the United-residue force field: ab initio folding simulations of multichain proteins

The implementation of molecular dynamics with the united-residue (UNRES) force field is extended to treat multichain proteins. Constant temperature was maintained in the simulations with Berendsen or Langevin thermostats. The method was tested on three alpha-helical proteins (1G6U and GCN4-p1, each with two chains, and 1C94, with four chains). Simulations were carried out for both the isolated single chains and the multichain complexes. The proteins were folded by starting from the extended conformation with random initial velocities and with the chains parallel to each other. No symmetry constraints or structure information were included for the single chains or the multichain complexes. In the case of single-chain simulations, a high percentage of the trajectories (100% for 1G6U, 90% for GCN4-p1, and 80% for 1C94) converged to nativelike structures (assumed as the experimental structure of a monomer in the multichain complex), showing that, for the proteins studied in this work with the UNRES force field, the interactions between chains are not critical for stabilization of the individual chains. In the case of multichain simulations, the native structures of the 1G6U and GCN4-p1 complexes, but not that of 1C94, are predicted successfully. The association of the subunits does not follow a unique mechanism; the monomers were observed to fold both before and simultaneously with their association.

Probing structural differences in prion protein isoforms by tyrosine nitration

Two conformational isomers of recombinant hamster prion protein (residues 90-232) have been probed by reaction with two tyrosine nitration reagents, peroxynitrite and tetranitromethane. Two conserved tyrosine residues (tyrosines 149 and 150) are not labeled by either reagent in the normal cellular form of the prion protein. These residues become reactive after the protein has been converted to the beta-oligomeric isoform, which is used as a model of the fibrillar form that causes disease. After conversion, a decrease in reactivity is noted for two other conserved residues, tyrosine 225 and tyrosine 226, whereas little to no effect was observed for other tyrosines. Thus, tyrosine nitration has identified two specific regions of the normal prion protein isoform that undergo a change in chemical environment upon conversion to a structure that is enriched in beta-sheet.

Structural properties of prion protein protofibrils and fibrils: an experimental assessment of atomic models

Decades after the prion protein was implicated in transmissible spongiform encephalopathies, the structure of its toxic isoform and its mechanism of toxicity remain unknown. By gathering available experimental data, albeit low resolution, a few pieces of the prion puzzle can be put in place. Currently, there are two fundamentally different models of a prion protofibril. One has its building blocks derived from a molecular dynamics simulation of the prion protein under amyloidogenic conditions, termed the spiral model. The other model was constructed by threading a portion of the prion sequence through a beta-helical structure from the Protein Data Bank. Here we compare and contrast these models with respect to all of the available experimental information, including electron micrographs, symmetries, secondary structure, oligomerization interfaces, enzymatic digestion, epitope exposure, and disaggregation profiles. Much of this information was not available when the two models were introduced. Overall, we find that the spiral model is consistent with all of the experimental results. In contrast, it is difficult to reconcile several of the experimental observables with the beta-helix model. While the experimental constraints are of low resolution, in bringing together the previously disconnected experiments, we have developed a clearer picture of prion aggregates. Both the improved characterization of prion aggregates and the existing atomic models can be used to devise further experiments to better elucidate the misfolding pathway and the structure of prion protofibrils.

Conformational change and assembly through edge beta strands in transthyretin and other amyloid proteins

Numerous diseases are characterized by the formation of insoluble, amyloid protein fibrils. Intensive investigations are beginning to unravel the detailed molecular and structural principles that underlie the spontaneous formation of these fibrils. The amyloid protein transthyretin serves as an excellent system for dissecting the conformational changes and ensuing subunit-subunit associations that lead to amyloid. One working model for tranthyretin amyloid involves the exposure of an "unprotected" edge beta strand, followed by symmetric assembly of subunits to give head-to-head and tail-to-tail protofibrils. The models and principles emerging from studies on transthyretin lead to connections to other amyloid systems.

A mechanism for copper inhibition of infectious prion conversion

We employ ab initio electronic structure calculations to obtain two structural models for copper bound in the strongest binding site of the noninfectious form of the prion protein. The models are compatible with available experimental constraints from electron spin resonance data. The bending of the peptide backbone attendant with the copper binding is not compatible with the requisite straight beta-strand backbone structure for the same sequence contained in two recently proposed models of the prion protein structure in its infectious form. We hypothesize that copper binding at this site is protective against conversion to the infectious form, discuss experimental data that appear to support and conflict with our hypothesis, and propose tests using recombinant prion protein, genetically modified cultured neurons, and transgenic mice.

Prion disease: exponential growth requires membrane binding

A hallmark feature of prions, whether in mammals or yeast and fungi, is exponential growth associated with fission or autocatalysis of protein aggregates. We have employed a rigorous kinetic analysis to recent data from transgenic mice lacking a glycosylphosphatidylinositol membrane anchor to the normal cellular PrP(C) protein, which show that toxicity requires the membrane binding. We find as well that the membrane is necessary for exponential growth of prion aggregates; without it, the kinetics is simply the quadratic-in-time growth characteristic of linear elongation as observed frequently in in vitro amyloid growth experiments with other proteins. This requires both: i), a substantial intercellular concentration of anchorless PrP(C), and ii), a concentration of small scrapies seeding aggregates from the inoculum, which remains relatively constant with time and exceeds the concentration of large polymeric aggregates. We also can explain via this analysis why mice heterozygous for the anchor-full/anchor-free PrP(C) proteins have more rapid incubation than mice heterozygous for anchor-full/null PrP(C), and contrast the mammalian membrane associated fission or autocatalysis with the membrane free fission of yeast and fungal prions.

Annealing prion protein amyloid fibrils at high temperature results in extension of a proteinase K-resistant core

Amyloids are highly ordered, rigid beta-sheet-rich structures that appear to have minimal dynamic flexibility in individual polypeptide chains. Here, we demonstrate that substantial conformational rearrangements occur within mature amyloid fibrils produced from full-length mammalian prion protein. The rearrangement results in a substantial extension of a proteinase K-resistant core and is accompanied by an increase in the beta-sheet-rich conformation. The conformational rearrangement was induced in the presence of low concentrations of Triton X-100 either by brief exposure to 80 degrees C or, with less efficacy, by prolonged incubation at 37 degrees C at pH 7.5 and is referred to here as "annealing." Upon annealing, amyloid fibrils acquired a proteinase K-resistant core identical to that found in bovine spongiform encephalopathy-specific scrapie-associated prion protein. Annealing was also observed when amyloid fibrils were exposed to high temperatures in the absence of detergent but in the presence of brain homogenate. These findings suggest that the amyloid fibrils exist in two conformationally distinct states that are separated by a high energy barrier and that yet unknown cellular cofactors may facilitate transition of the fibrils into thermodynamically more stable state. Our studies provide new insight into the complex behavior of prion polymerization and highlight the annealing process, a previously unknown step in the evolution of amyloid structures.

Protein oligomerization through domain swapping: role of inter-molecular interactions and protein concentration

Domain swapping has been shown to be an important mechanism controlling multiprotein assembly and has been suggested recently as a possible mechanism underlying protein aggregation. Understanding oligomerization via domain swapping is therefore of theoretical and practical importance. By using a symmetrized structure-based (Gō) model, we demonstrate that in the free-energy landscape of domain swapping, a large free-energy barrier separates monomeric and domain-swapped dimeric configurations. We investigate the effect of finite monomer concentration, by implementing a new semi-analytical method, which involves computing the second virial coefficient, a thermodynamic indicator of inter-molecular interactions. This method, together with the symmetrized structure-based (Gō) model, minimizes the need for expensive many-protein simulations, providing a convenient framework to investigate concentration effect. Finally, we perform direct simulations of domain-swapped trimer formation, showing that this modeling approach can be used for higher-order oligomers.