Characterization of a possible amyloidogenic precursor in glutamine-repeat neurodegenerative diseases

The role of α-sheet structure in amyloidogenesis: characterization and implications

Amyloid diseases are linked to protein misfolding whereby the amyloidogenic protein undergoes a conformational change, aggregates and eventually forms amyloid fibrils. While the amyloid fibrils and plaques are hallmarks of these diseases, they typically form late in the disease process and do not correlate with disease. Instead, there is growing evidence that smaller, soluble toxic oligomers form prior and appear to be early triggers of the molecular pathology underlying these diseases. Nearly 20 years ago, we proposed the α-sheet hypothesis after discovering that the early conformational changes observed during atomistic molecular dynamics simulations involve the formation of a non-standard protein structure, α-sheet. Furthermore, we proposed that toxic oligomers contain α-sheet structure and that preferentially targeting this structure could neutralize the toxicity, prevent further aggregation and serve as the basis for early detection of disease. Here, we present the origin of the α-sheet hypothesis and describe α-sheet structure and the corresponding mechanisms of conversion. We discuss experimental studies demonstrating that both mammalian and bacterial amyloid systems form α-sheet oligomers before converting to conventional β-sheet fibrils. Furthermore, we show that the process can be inhibited with de novo designed α-sheet peptides complementary to the structure in the toxic oligomers.

A Robust Assay to Monitor Ataxin-3 Amyloid Fibril Assembly

Spinocerebellar ataxia type 3 (SCA3) is caused by the expansion of a glutamine repeat in the protein ataxin-3, which is deposited as intracellular aggregates in affected brain regions. Despite the controversial role of ataxin-3 amyloid structures in SCA3 pathology, the identification of molecules with the capacity to prevent aberrant self-assembly and stabilize functional conformation(s) of ataxin-3 is a key to the development of therapeutic solutions. Amyloid-specific kinetic assays are routinely used to measure rates of protein self-assembly in vitro and are employed during screening for fibrillation inhibitors. The high tendency of ataxin-3 to assemble into oligomeric structures implies that minor changes in experimental conditions can modify ataxin-3 amyloid assembly kinetics. Here, we determine the self-association rates of ataxin-3 and present a detailed study of the aggregation of normal and pathogenic ataxin-3, highlighting the experimental conditions that should be considered when implementing and validating ataxin-3 amyloid progress curves in different settings and in the presence of ataxin-3 interactors. This assay provides a unique and robust platform to screen for modulators of the first steps of ataxin-3 aggregation—a starting point for further studies with cell and animal models of SCA3.

Tumorigenic p53 mutants undergo common structural disruptions including conversion to α-sheet structure

The p53 protein is a commonly studied cancer target because of its role in tumor suppression. Unfortunately, it is susceptible to mutation-associated loss of function; approximately 50% of cancers are associated with mutations to p53, the majority of which are located in the central DNA-binding domain. Here, we report molecular dynamics simulations of wild-type (WT) p53 and 20 different mutants, including a stabilized pseudo-WT mutant. Our findings indicate that p53 mutants tend to exacerbate latent structural-disruption tendencies, or vulnerabilities, already present in the WT protein, suggesting that it may be possible to develop cancer therapies by targeting a relatively small set of structural-disruption motifs rather than a multitude of effects specific to each mutant. In addition, α-sheet secondary structure formed in almost all of the proteins. α-Sheet has been hypothesized and recently demonstrated to play a role in amyloidogenesis, and its presence in the reported p53 simulations coincides with the recent re-consideration of cancer as an amyloid disease.

Aggregation Mechanism of Alzheimer's Amyloid β-Peptide Mediated by α-Strand/α-Sheet Structure

Alzheimer's disease (AD) is one of the most common neurodegenerative diseases and a widespread form of dementia. Aggregated forms of the amyloid β-peptide (Aβ) are identified as a toxic species responsible for neuronal damage in AD. Extensive research has been conducted to reveal the aggregation mechanism of Aβ. However, the structure of pathological aggregates and the mechanism of aggregation are not well understood. Recently, experimental studies have confirmed that the α-sheet structure in Aβ drives aggregation and toxicity in AD. However, how the α-sheet structure is formed in Aβ and how it contributes to Aβ aggregation remains elusive. In the present study, molecular dynamics simulations suggest that Aβ adopts the α-strand conformation by peptide-plane flipping. Multiple α-strands interact through hydrogen bonding to form α-sheets. This structure acts as a nucleus that initiates and promotes aggregation and fibrillation of Aβ. Our findings are supported by previous experimental as well as theoretical studies. This study provides valuable structural insights for the design of anti-AD drugs exploiting the α-strand/α-sheet structure.

Molecular Dynamics Simulations Applied to Structural and Dynamical Transitions of the Huntingtin Protein: A Review

Over the recent years, Huntington's disease (HD) has become widely discussed in the scientific literature especially because at the mutant level there are several contradictions regarding the aggregation mechanism. The specific role of the physiological huntingtin protein remains unknown, due to the lack of characterization of its entire crystallographic structure, making the experimental and theoretical research even harder when taking into consideration its involvement in multiple biological functions and its high affinity for different interacting partners. Different types of models, containing fewer (not more than 35 Qs) polyglutamine residues for the WT structure and above 35 Qs for the mutants, were subjected to classical or advanced MD simulations to establish the proteins' structural stability by evaluating their conformational changes. Outside the polyQ tract, there are two other regions of interest (the N17 domain and the polyP rich domain) considered to be essential for the aggregation kinetics at the mutant level. The polymerization process is considered to be dependent on the polyQ length. As the polyQ tract's dimension increases, the structures present more β-sheet conformations. Contrarily, it is also considered that the aggregation stability is not necessarily dependent on the number of Qs, while the initial stage of the aggregation seed might play the decisive role. A general assumption regarding the polyP domain is that it might preserve the polyQ structures soluble by acting as an antagonist for β-sheet formation.

Local structural motifs in proteins: Detection and characterization of fragments inserted in helices

The global/local fold of protein structures is stabilized by a variety of specific interactions. A primary role in this context is played by hydrogen bonds. In order to identify novel motifs in proteins, we searched Protein Data Bank structures looking for backbone H-bonds formed by NH groups of two (or more) consecutive residues with consecutive CO groups of distant residues in the sequence. The present analysis unravels the occurrence of recurrent structural motifs that, to the best of our knowledge, had not been characterized in literature. Indeed, these H-bonding patterns are found (i) in a specific parallel β-sheet capping, (ii) in linking of β-hairpins to α-helices, and (iii) in α-helix insertions. Interestingly, structural analyses of these motifs indicate that Gly residues frequently occupy prominent positions. The formation of these motifs is likely favored by the limited propensity of Gly to be embodied in helices/sheets. Of particular interest is the motif corresponding to insertions in helices that was detected in 1% of analyzed structures. Inserted fragments may assume different structures and aminoacid compositions and usually display diversified evolutionary conservation. Since inserted regions are physically separated from the rest of the protein structure, they represent hot spots for ad-hoc protein functionalization.

Identification of a novel site of interaction between ataxin-3 and the amyloid aggregation inhibitor polyglutamine binding peptide 1

Amyloid diseases represent a growing social and economic burden in the developed world. Understanding the assembly pathway and the inhibition of amyloid formation is key to developing therapies to treat these diseases. The neurodegenerative condition Machado-Joseph disease is characterised by the self-aggregation of the protein ataxin-3. Ataxin-3 consists of a globular N-terminal Josephin domain, which can aggregate into curvilinear protofibrils, and an unstructured, dynamically disordered C-terminal domain containing three ubiquitin interacting motifs separated by a polyglutamine stretch. Upon expansion of the polyglutamine region above 50 residues, ataxin-3 undergoes a second stage of aggregation in which long, straight amyloid fibrils form. A peptide inhibitor of polyglutamine aggregation, known as polyQ binding peptide 1, has been shown previously to prevent the maturation of ataxin-3 fibrils. However, the mechanism of this inhibition remains unclear. Using nanoelectrospray ionisation-mass spectrometry, we demonstrate that polyQ binding peptide 1 binds to monomeric ataxin-3. By investigating the ability of polyQ binding peptide 1 to bind to truncated ataxin-3 constructs lacking one or more domains, we localise the site of this interaction to a 39-residue sequence immediately C-terminal to the Josephin domain. The results suggest a new mechanism for the inhibition of polyglutamine aggregation by polyQ binding peptide 1 in which binding to a region outside of the polyglutamine tract can prevent fibril formation, highlighting the importance of polyglutamine flanking regions in controlling aggregation and disease.

Possible Existence of α-Sheets in the Amyloid Fibrils Formed by a TTR105-115 Mutant

Herein, we combine several methods to characterize the fibrils formed by a TTR105-115 mutant in which Leu111 is replaced by the unnatural amino acid aspartic acid 4-methyl ester. We find that this mutant peptide exhibits significantly different aggregation behavior than the wild-type peptide: (1) it forms fibrils with a much faster rate, (2) its fibrils lack the long-range helical twists observed in TTR105-115 fibrils, (3) its fibrils exhibit a giant far-UV circular dichroism signal, and (4) its fibrils give rise to an unusual amide I' band consisting of four distinct and sharp peaks. On the basis of these results and also several previous computational studies, we hypothesize that the fibrils formed by this TTR mutant peptide contain both β- and α-sheets.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

Insights into Unfolded Proteins from the Intrinsic ϕ/ψ Propensities of the AAXAA Host-Guest Series

Various host-guest peptide series are used by experimentalists as reference conformational states. One such use is as a baseline for random-coil NMR chemical shifts. Comparison to this random-coil baseline, through secondary chemical shifts, is used to infer protein secondary structure. The use of these random-coil data sets rests on the perception that the reference chemical shifts arise from states where there is little or no conformational bias. However, there is growing evidence that the conformational composition of natively and nonnatively unfolded proteins fail to approach anything that can be construed as random coil. Here, we use molecular dynamics simulations of an alanine-based host-guest peptide series (AAXAA) as a model of unfolded and denatured states to examine the intrinsic propensities of the amino acids. We produced ensembles that are in good agreement with the experimental NMR chemical shifts and confirm that the sampling of the 20 natural amino acids in this peptide series is be far from random. Preferences toward certain regions of conformational space were both present and dependent upon the environment when compared under conditions typically used to denature proteins, i.e., thermal and chemical denaturation. Moreover, the simulations allowed us to examine the conformational makeup of the underlying ensembles giving rise to the ensemble-averaged chemical shifts. We present these data as an intrinsic backbone propensity library that forms part of our Structural Library of Intrinsic Residue Propensities to inform model building, to aid in interpretation of experiment, and for structure prediction of natively and nonnatively unfolded states.

New Dynamic Rotamer Libraries: Data-Driven Analysis of Side-Chain Conformational Propensities

Most rotamer libraries are generated from subsets of the PDB and do not fully represent the conformational scope of protein side chains. Previous attempts to rectify this sparse coverage of conformational space have involved application of weighting and smoothing functions. We resolve these limitations by using physics-based molecular dynamics simulations to determine more accurate frequencies of rotameric states. This work forms part of our Dynameomics initiative and uses a set of 807 proteins selected to represent 97% of known autonomous protein folds, thereby eliminating the bias toward common topologies found within the PDB. Our Dynameomics derived rotamer libraries encompass 4.8 × 10(9) rotamers, sampled from at least 51,000 occurrences of each of 93,642 residues. Here, we provide a backbone-dependent rotamer library, based on secondary structure ϕ/ψ regions, and an update to our 2011 backbone-independent library that addresses the doubling of our dataset since its original publication.

Peptides Composed of Alternating L- and D-Amino Acids Inhibit Amyloidogenesis in Three Distinct Amyloid Systems Independent of Sequence

There is now substantial evidence that soluble oligomers are primary toxic agents in amyloid diseases. The development of an antibody recognizing the toxic soluble oligomeric forms of different and unrelated amyloid species suggests a common conformational intermediate during amyloidogenesis. We previously observed a common occurrence of a novel secondary structure element, which we call α-sheet, in molecular dynamics (MD) simulations of various amyloidogenic proteins, and we hypothesized that the toxic conformer is composed of α-sheet structure. As such, α-sheet may represent a conformational signature of the misfolded intermediates of amyloidogenesis and a potential unique binding target for peptide inhibitors. Recently, we reported the design and characterization of a novel hairpin peptide (α1 or AP90) that adopts stable α-sheet structure and inhibits the aggregation of the β-Amyloid Peptide Aβ42 and transthyretin. AP90 is a 23-residue hairpin peptide featuring alternating D- and L-amino acids with favorable conformational propensities for α-sheet formation, and a designed turn. For this study, we reverse engineered AP90 to identify which of its design features is most responsible for conferring α-sheet stability and inhibitory activity. We present experimental characterization (CD and FTIR) of seven peptides designed to accomplish this. In addition, we measured their ability to inhibit aggregation in three unrelated amyloid species: Aβ42, transthyretin, and human islet amylin polypeptide. We found that a hairpin peptide featuring alternating L- and D-amino acids, independent of sequence, is sufficient for conferring α-sheet structure and inhibition of aggregation. Additionally, we show a correlation between α-sheet structural stability and inhibitory activity.

Huntingtin exon 1 fibrils feature an interdigitated β-hairpin-based polyglutamine core

Polyglutamine expansion within the exon1 of huntingtin leads to protein misfolding, aggregation, and cytotoxicity in Huntington's disease. This incurable neurodegenerative disease is the most prevalent member of a family of CAG repeat expansion disorders. Although mature exon1 fibrils are viable candidates for the toxic species, their molecular structure and how they form have remained poorly understood. Using advanced magic angle spinning solid-state NMR, we directly probe the structure of the rigid core that is at the heart of huntingtin exon1 fibrils and other polyglutamine aggregates, via measurements of long-range intramolecular and intermolecular contacts, backbone and side-chain torsion angles, relaxation measurements, and calculations of chemical shifts. These experiments reveal the presence of β-hairpin-containing β-sheets that are connected through interdigitating extended side chains. Despite dramatic differences in aggregation behavior, huntingtin exon1 fibrils and other polyglutamine-based aggregates contain identical β-strand-based cores. Prior structural models, derived from X-ray fiber diffraction and computational analyses, are shown to be inconsistent with the solid-state NMR results. Internally, the polyglutamine amyloid fibrils are coassembled from differently structured monomers, which we describe as a type of "intrinsic" polymorphism. A stochastic polyglutamine-specific aggregation mechanism is introduced to explain this phenomenon. We show that the aggregation of mutant huntingtin exon1 proceeds via an intramolecular collapse of the expanded polyglutamine domain and discuss the implications of this observation for our understanding of its misfolding and aggregation mechanisms.

An Analysis of Biomolecular Force Fields for Simulations of Polyglutamine in Solution

Polyglutamine (polyQ) peptides are a useful model system for biophysical studies of protein folding and aggregation, both for their intriguing aggregation properties and their own relevance to human disease. The genetic expansion of a polyQ tract triggers the formation of amyloid aggregates associated with nine neurodegenerative diseases. Several clearly identifiable and separable factors, notably the length of the polyQ tract, influence the mechanism of aggregation, its associated kinetics, and the ensemble of structures formed. Atomistic simulations are well positioned to answer open questions regarding the thermodynamics and kinetics of polyQ folding and aggregation. The additional, explicit representation of water permits deeper investigation of the role of solvent dynamics, and it permits a direct comparison of simulation results with infrared spectroscopy experiments. The generation of meaningful simulation results hinges on satisfying two essential criteria: achieving sufficient conformational sampling to draw statistically valid conclusions, and accurately reproducing the intermolecular forces that govern system structure and dynamics. In this work, we examine the ability of 12 biomolecular force fields to reproduce the properties of a simple, 30-residue polyQ peptide (Q30) in explicit water. In addition to secondary and tertiary structure, we consider generic structural properties of polymers that provide additional dimensions for analysis of the highly degenerate disordered states of the molecule. We find that the 12 force fields produce a wide range of predictions. We identify AMBER ff99SB, AMBER ff99SB*, and OPLS-AA/L to be most suitable for studies of polyQ folding and aggregation.

Structural determinants of polyglutamine protofibrils and crystallites

Nine inherited neurodegenerative diseases are associated with the expansion of the CAG codon. Once the translated polyglutamine expansion becomes longer than ~36 residues, it triggers the formation of intraneural protein aggregates that often display the signature of cross-β amyloid fibrils. Here, we use fully atomistic molecular dynamics simulations to probe the structural stability and conformational dynamics of both previously proposed and new polyglutamine aggregate models. We test the relative stability of parallel and antiparallel β sheets, and characterize possible steric interfaces between neighboring sheets and the effects of different alignments of the side-chain carboxamide dipoles. Results indicate that (i) different initial oligomer structures converge to crystals consistent with available diffraction data, after undergoing cooperative side-chain rotational transitions and quarter-stagger displacements on a microsecond time scale, (ii) structures previously deemed stable on a hundred nanosecond time scale are unstable over the microsecond time scale, and (iii) conversely, structures previously deemed unstable did not account for the correct side-chain packing and once the correct symmetry is considered the structures become stable for over a microsecond, due to tightly interdigitated side chains, which lock into highly regular polar zippers with inter-side-chain and backbone-side-chain hydrogen bonds. With these insights, we built Q40 monomeric models with different combinations of arc and hairpin turns and tested them for stability. The stable monomers were further probed as a function of repeat length. Our results are consistent with the aggregation threshold. These results explain and reconcile previously reported experimental and model discrepancies about polyglutamine aggregate structures.

Structural and dynamic properties of the human prion protein

Prion diseases involve the conformational conversion of the cellular prion protein (PrP(C)) to its misfolded pathogenic form (PrP(Sc)). To better understand the structural mechanism of this conversion, we performed extensive all-atom, explicit-solvent molecular-dynamics simulations for three structures of the wild-type human PrP (huPrP) at different pH values and temperatures. Residue 129 is polymorphic, being either Met or Val. Two of the three structures have Met in position 129 and the other has Val. Lowering the pH or raising the temperature induced large conformational changes of the C-terminal globular domain and increased exposure of its hydrophobic core. In some simulations, HA and its preceding S1-HA loop underwent large displacements. The C-terminus of HB was unstable and sometimes partially unfolded. Two hydrophobic residues, Phe-198 and Met-134, frequently became exposed to solvent. These conformational changes became more dramatic at lower pH or higher temperature. Furthermore, Tyr-169 and the S2-HB loop, or the X-loop, were different in the starting structures but converged to common conformations in the simulations for the Met-129, but not the Val-129, protein. α-Strands and β-strands formed in the initially unstructured N-terminus. α-Strand propensity in the N-terminus was different between the Met-129 and Val129 proteins, but β-strand propensity was similar. This study reveals detailed structural and dynamic properties of huPrP, providing insight into the mechanism of the conversion of PrP(C) to PrP(Sc).

Designed α-sheet peptides inhibit amyloid formation by targeting toxic oligomers

Previous studies suggest that the toxic soluble-oligomeric form of different amyloid proteins share a common backbone conformation, but the amorphous nature of this oligomer prevents its structural characterization by experiment. Based on molecular dynamics simulations we proposed that toxic intermediates of different amyloid proteins adopt a common, nonstandard secondary structure, called α-sheet. Here we report the experimental characterization of peptides designed to be complementary to the α-sheet conformation observed in the simulations. We demonstrate inhibition of aggregation in two different amyloid systems, β-amyloid peptide (Aβ) and transthyretin, by these designed α-sheet peptides. When immobilized the α-sheet designs preferentially bind species from solutions enriched in the toxic conformer compared with non-aggregated, nontoxic species or mature fibrils. The designs display characteristic spectroscopic signatures distinguishing them from conventional secondary structures, supporting α-sheet as a structure involved in the toxic oligomer stage of amyloid formation and paving the way for novel therapeutics and diagnostics.

DOI:http://dx.doi.org/10.7554/eLife.01681.001

The build up of very thin fibres called amyloid fibrils is known to lead to more than 40 different human diseases, including Parkinson’s disease and rheumatoid arthritis. These diseases involve soluble proteins or peptides joining other proteins or peptides to form the fibrils, which are not soluble. However, the damage is done by the time the fibrils form because soluble intermediate structures formed by the proteins and peptides are toxic. The development of methods that can detect these toxic intermediate structures could lead to earlier interventions before significant damage.

Amyloid fibrils are known to have a beta-sheet structure that is found in many protein systems. In 2004, based on computer simulations, researchers predicted that proteins and peptides that go on to form amyloid fibrils would pass through a related but less stable structure called an alpha-sheet, and that this structure would be toxic. Now Hopping et al., including some of the researchers involved in the 2004 work, have confirmed that the alpha-sheet structure is indeed involved in the formation of amyloid fibrils.

To do this Hopping et al. designed peptides with alpha-sheet structures that could bind to the alpha-sheet structures predicted by their simulations. When these complementary designed peptides were added to a solution of peptide that causes Alzheimer’s Disease, or a protein that causes systemic amyloid disease, the designed peptides bound the toxic peptides or proteins and prevented the formation of fibrils.

The results of Hopping et al. suggest that designed alpha-sheet compounds might be able to capture peptides and proteins that are implicated in a wide variety of amyloid diseases, independent of their composition and native structure, by targeting the intermediate alpha-sheet structure. Future challenges include showing that most proteins and peptides pass through this intermediate structure as they form fibrils, and improving the sensitivity of the binding in the hope of developing diagnostics for amyloid diseases.

DOI:http://dx.doi.org/10.7554/eLife.01681.002

NMR spectroscopy reveals a preferred conformation with a defined hydrophobic cluster for polyglutamine binding peptide 1

Several important human inherited neurodegenerative diseases are caused by "polyQ expansions", which are aberrant long repeats of glutamine residues in proteins. PolyQ binding peptide 1 (QBP1), whose minimal active core sequence is Trp-Lys-Trp-Trp-Pro-Gly-Ile-Phe, binds to expanded polyQs and blocks their β-structure transition, aggregation and in vivo neurodegeneration. Whereas QBP1 is a widely used, commercially available product, its structure is unknown. Here, we have characterized the conformations of QBP1 and a scrambled peptide (Trp-Pro-Ile-Trp-Lys-Gly-Trp-Phe) in aqueous solution by CD, fluorescence and NMR spectroscopies. A CD maximum at 227 nm suggests the presence of rigid Trp side chains in QBP1. Based on 41 NOE-derived distance constraints, the 3D structure of QBP1 was determined. The side chains of Trp 4 and Ile 7, and to a lesser extent, those of Lys 2, Trp 3 and Phe 8, form a small hydrophobic cluster. Pro 5 and Gly 6 adopt a type II tight turn and Lys 2's ζ-NH3(+) is positioned to form a favorable cation-π interaction with Trp 4's indole ring. In contrast, the scrambled QBP1 peptide, which lacks inhibitory activity, does not adopt a preferred structure. These results provide a basis for future structure-based design approaches to further optimize QBP1 for therapy.

Computational approaches to understanding protein aggregation in neurodegeneration

The generation of toxic non-native protein conformers has emerged as a unifying thread among disorders such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. Atomic-level detail regarding dynamical changes that facilitate protein aggregation, as well as the structural features of large-scale ordered aggregates and soluble non-native oligomers, would contribute significantly to current understanding of these complex phenomena and offer potential strategies for inhibiting formation of cytotoxic species. However, experimental limitations often preclude the acquisition of high-resolution structural and mechanistic information for aggregating systems. Computational methods, particularly those combine both all-atom and coarse-grained simulations to cover a wide range of time and length scales, have thus emerged as crucial tools for investigating protein aggregation. Here we review the current state of computational methodology for the study of protein self-assembly, with a focus on the application of these methods toward understanding of protein aggregates in human neurodegenerative disorders.

Molecular dynamics simulations capture the misfolding of the bovine prion protein at acidic pH

Bovine spongiform encephalopathy (BSE), or mad cow disease, is a fatal neurodegenerative disease that is transmissible to humans and that is currently incurable. BSE is caused by the prion protein (PrP), which adopts two conformers; PrPC is the native innocuous form, which is α-helix rich; and PrPSc is the β-sheet rich misfolded form, which is infectious and forms neurotoxic species. Acidic pH induces the conversion of PrPC to PrPSc. We have performed molecular dynamics simulations of bovine PrP at various pH regimes. An acidic pH environment induced conformational changes that were not observed in neutral pH simulations. Putative misfolded structures, with nonnative β-strands formed in the flexible N-terminal domain, were found in acidic pH simulations. Two distinct pathways were observed for the formation of nonnative β-strands: at low pH, hydrophobic contacts with M129 nucleated the nonnative β-strand; at mid-pH, polar contacts involving Q168 and D178 facilitated the formation of a hairpin at the flexible N-terminus. These mid- and low pH simulations capture the process of nonnative β-strand formation, thereby improving our understanding of how PrPC misfolds into the β-sheet rich PrPSc and how pH factors into the process.

Structural consequences of mutations to the α-tocopherol transfer protein associated with the neurodegenerative disease ataxia with vitamin E deficiency

The α-tocopherol transfer protein (α-TTP) is a liver protein that transfers α-tocopherol (vitamin E) to very-low-density lipoproteins (VLDLs). These VLDLs are then circulated throughout the body to maintain blood α-tocopherol levels. Mutations to the α-TTP gene are associated with ataxia with vitamin E deficiency, a disease characterized by peripheral nerve degeneration. In this study, molecular dynamics simulations of the E141K and R59W disease-associated mutants were performed. The mutants displayed disruptions in and around the ligand-binding pocket. Structural analysis and ligand docking to the mutant structures predicted a decreased affinity for α-tocopherol. To determine the detailed mechanism of the mutation-related changes, we developed a new tool called ContactWalker that analyzes contact differences between mutant and wild-type proteins and highlights pathways of altered contacts within the mutant proteins. Taken together, our findings are in agreement with experiment and suggest structural explanations for the weakened ability of the mutants to bind and carry α-tocopherol.

Assessing polyglutamine conformation in the nucleating event by molecular dynamics simulations

Polyglutamine (polyQ) diseases comprise a group of dominantly inherited pathology caused by an expansion of an unstable polyQ stretch which is presumed to form β-sheets. Similar to other amyloid pathologies, polyQ amyloidogenesis occurs via a nucleated polymerization mechanism, and proceeds through energetically unfavorable nucleus whose existence and structure are difficult to detect. Here, we use atomistic molecular dynamics simulations in explicit solvent to assess the conformation of the polyQ stretch in the nucleus that initiates polyQ fibrillization. Comparison of the kinetic stability of various structures of polyQ peptide with a Q-length in the pathological range (Q40) revealed that steric zipper or nanotube-like structures (β-nanotube or β-pseudohelix) are not kinetically stable enough to serve as a template to initiate polyQ fibrillization as opposed to β-hairpin-based (β-sheet and β-sheetstack) or α-helical conformations. The selection of different structures of the polyQ stretch in the aggregation-initiating event may provide an alternative explanation for polyQ aggregate polymorphism.

Common features at the start of the neurodegeneration cascade

A single-molecule study reveals that neurotoxic proteins share common structural features that may trigger neurodegeneration, thus identifying new targets for therapy and diagnosis.

Amyloidogenic neurodegenerative diseases are incurable conditions with high social impact that are typically caused by specific, largely disordered proteins. However, the underlying molecular mechanism remains elusive to established techniques. A favored hypothesis postulates that a critical conformational change in the monomer (an ideal therapeutic target) in these “neurotoxic proteins” triggers the pathogenic cascade. We use force spectroscopy and a novel methodology for unequivocal single-molecule identification to demonstrate a rich conformational polymorphism in the monomer of four representative neurotoxic proteins. This polymorphism strongly correlates with amyloidogenesis and neurotoxicity: it is absent in a fibrillization-incompetent mutant, favored by familial-disease mutations and diminished by a surprisingly promiscuous inhibitor of the critical monomeric β-conformational change, neurotoxicity, and neurodegeneration. Hence, we postulate that specific mechanostable conformers are the cause of these diseases, representing important new early-diagnostic and therapeutic targets. The demonstrated ability to inhibit the conformational heterogeneity of these proteins by a single pharmacological agent reveals common features in the monomer and suggests a common pathway to diagnose, prevent, halt, or reverse multiple neurodegenerative diseases.

Neurodegenerative diseases like Alzheimer's or Parkinson's are currently incurable. They are caused by different proteins that, under certain circumstances, aggregate and become toxic as we grow older, but the molecular events underlying this process remain unclear. The lack of a well-defined structure, and the tendency of these “neurotoxic proteins” to aggregate make them difficult to study using conventional techniques. Here, we use an established single-molecule manipulation technique combined with a new protein-engineering strategy to show that all these proteins can adopt a rich collection of structures (conformers) that includes a high proportion of mechanostable conformers, which are associated with toxicity and disease. We also find that a known drug can block the formation of these mechanostable structures in different neurotoxic proteins. We suggest that the most mechanostable conformers, or their precursors, may trigger the pathogenic cascade that results in toxicity. We thus propose that these mechanostable structures are ideal targets for early diagnosis, prevention, and treatment of these fatal diseases.

Are long-range structural correlations behind the aggregration phenomena of polyglutamine diseases?

We have characterized the conformational ensembles of polyglutamine peptides of various lengths (ranging from to ), both with and without the presence of a C-terminal polyproline hexapeptide. For this, we used state-of-the-art molecular dynamics simulations combined with a novel statistical analysis to characterize the various properties of the backbone dihedral angles and secondary structural motifs of the glutamine residues. For (i.e., just above the pathological length for Huntington's disease), the equilibrium conformations of the monomer consist primarily of disordered, compact structures with non-negligible -helical and turn content. We also observed a relatively small population of extended structures suitable for forming aggregates including - and -strands, and - and -hairpins. Most importantly, for we find that there exists a long-range correlation (ranging for at least residues) among the backbone dihedral angles of the Q residues. For polyglutamine peptides below the pathological length, the population of the extended strands and hairpins is considerably smaller, and the correlations are short-range (at most residues apart). Adding a C-terminal hexaproline to suppresses both the population of these rare motifs and the long-range correlation of the dihedral angles. We argue that the long-range correlation of the polyglutamine homopeptide, along with the presence of these rare motifs, could be responsible for its aggregation phenomena.

Nine neurodegenerative diseases are caused by polyglutamine (polyQ) expansions greater than a given threshold in proteins with little or no homology except for the polyQ regions. The diseases all share a common feature: the formation of polyQ aggregates and eventual neuronal death. Using molecular dynamics simulations, we have explored the conformations of polyQ peptides. Results indicate that for peptides (i.e., just above the pathological length for Hungtington's disease), the equilibrium conformations were found to consist primarily of disordered, compact structures with a non-negligible -helical and turn content. We also observed a small population of extended structures suitable for forming aggregates. For peptides below the pathological length, the population of these structures was found to be considerably lower. For longer peptides, we found evidence for long-range correlations among the dihedral angles. This correlation turns out to be short-range for the smaller polyQ peptides, and is suppressed (along with the extended structural motifs) when a C-terminal polyproline tail is added to the peptides. We believe that the existence of these long-range correlations in above-threshold polyQ peptides, along with the presence of rare motifs, could be responsible for the experimentally observed aggregation phenomena associated with polyQ diseases.

Secondary structures of native and pathogenic huntingtin N-terminal fragments

Huntington's disease is a neurodegenerative disorder caused by a polyglutamine (polyQ) expansion in the N-terminal fragment of the Huntingtin (Htt) protein. Structural properties of Htt N-terminal regions and the molecular mechanism leading to protein aggregation have not been fully explained yet. We performed all-atom replica exchange molecular dynamics to investigate the structures of Htt N-terminal parts with polyQ tracts of nonpathogenic and pathogenic lengths. The monomers were composed of the headpiece (17 N-terminal residues), a polyQ tract (polyQ(17) for native and polyQ(55) for pathogenic sequence), and a polyP(11) region, followed by 17 amino acids of mixed sequence. We found that corresponding regions in both fragments fold to similar secondary structures; the headpiece and polyQ stretch adopt mainly α-helical conformations, and polyP(11) forms the PP II-type helix. The native N-terminal fragment is more compact and stabilized by hydrophobic interactions between the surface of polyP(11) and the amphipathic helix of the headpiece. In the pathogenic fragment the headpiece is solvent exposed and does not interact with polyP(11). The predicted structure of the native N-terminal fragment agrees with the X-ray structure of the Htt first exon containing polyQ(17). The structure of the pathogenic fragment adheres to an aggregation model that is mediated by the Htt headpiece.

The Aggregation Inhibitor Peptide QBP1 as a Therapeutic Molecule for the Polyglutamine Neurodegenerative Diseases

Misfolding and abnormal aggregation of proteins in the brain are implicated in the pathogenesis of various neurodegenerative diseases including Alzheimer's, Parkinson's, and the polyglutamine (polyQ) diseases. In the polyQ diseases, an abnormally expanded polyQ stretch triggers misfolding and aggregation of the disease-causing proteins, eventually resulting in neurodegeneration. In this paper, we introduce our therapeutic strategy against the polyQ diseases using polyQ binding peptide 1 (QBP1), a peptide that we identified by phage display screening. We showed that QBP1 specifically binds to the expanded polyQ stretch and inhibits its misfolding and aggregation, resulting in suppression of neurodegeneration in cell culture and animal models of the polyQ diseases. We further demonstrated the potential of protein transduction domains (PTDs) for in vivo delivery of QBP1. We hope that in the near future, chemical analogues of aggregation inhibitor peptides including QBP1 will be developed against protein misfolding-associated neurodegenerative diseases.

Malleability of folding intermediates in the homeodomain superfamily

Members of the homeodomain superfamily are three-helix bundle proteins whose second and third helices form a helix-turn-helix motif (HTH). Their folding mechanism slides from the ultrafast, three-state framework mechanism for the engrailed homeodomain (EnHD), in which the HTH motif is independently stable, to an apparent two-state nucleation-condensation model for family members with an unstable HTH motif. The folding intermediate of EnHD has nearly native HTH structure, but it is not docked with helix1. The determinant of whether two- or three-state folding was hypothesized to be the stability of the HTH substructure. Here, we describe a detailed Φ-value analysis of the folding of the Pit1 homeodomain, which has similar ultrafast kinetics to that of EnHD. Formation of helix1 was strongly coupled with formation of HTH, which was initially surprising because they are uncoupled in the EnHD folding intermediate. However, we found a key difference between Pit1 and EnHD: The isolated peptide corresponding to the HTH motif in Pit1 was not folded in the absence of H1. Independent molecular dynamics simulations of Pit1 unfolding found an intermediate with H1 misfolded onto the HTH motif. The Pit1 folding pathway is the connection between that of EnHD and the slower folding homeodomains and provides a link in the transition of mechanisms from two- to three-state folding in this superfamily. The malleability of folding intermediates can lead to unstable substructures being stabilized by a variety of nonnative interactions, adding to the continuum of folding mechanisms.

Manifestations of native topology in the denatured state ensemble of Rhodopseudomonas palustris cytochrome c'

To provide insight into the role of local sequence in the nonrandom coil behavior of the denatured state, we have extended our measurements of histidine-heme loop formation equilibria for cytochrome c' to 6 M guanidine hydrochloride. We observe that there is some reduction in the scatter about the best fit line of loop stability versus loop size data in 6 M versus 3 M guanidine hydrochloride, but the scatter is not eliminated. The scaling exponent, ν(3), of 2.5 ± 0.2 is also similar to that found previously in 3 M guanidine hydrochloride (2.6 ± 0.3). Rates of histidine-heme loop breakage in the denatured state of cytochrome c' show that some histidine-heme loops are significantly more persistent than others at both 3 and 6 M guanidine hydrochloride. Rates of histidine-heme loop formation more closely approximate random coil behavior. This observation indicates that heterogeneity in the denatured state ensemble results mainly from contact persistence. When mapped onto the structure of cytochrome c', the histidine-heme loops with slow breakage rates coincide with chain reversals between helices 1 and 2 and between helices 2 and 3. Molecular dynamics simulations of the unfolding of cytochrome c' at 498 K show that these reverse turns persist in the unfolded state. Thus, these portions of the primary structure of cytochrome c' set up the topology of cytochrome c' in the denatured state, predisposing the protein to fold efficiently to its native structure.

Diverse effects on the native β-sheet of the human prion protein due to disease-associated mutations

Prion diseases are fatal neurodegenerative disorders that involve the conversion of the normal cellular form of the prion protein (PrP(C)) to a misfolded pathogenic form (PrP(Sc)). There are many genetic mutations of PrP associated with human prion diseases. Three of these point mutations are located at the first strand of the native β-sheet in human PrP: G131V, S132I, and A133V. To understand the underlying structural and dynamic effects of these disease-causing mutations on the human PrP, we performed molecular dynamics of wild-type and mutated human PrP. The results indicate that the mutations induced different effects but they were all related to misfolding of the native β-sheet: G131V caused the elongation of the native β-sheet, A133V disrupted the native β-sheet, and S132I converted the native β-sheet to an α-sheet. The observed changes were due to the reorientation of side chain-side chain interactions upon introducing the mutations. In addition, all mutations impaired a structurally conserved water site at the native β-sheet. Our work suggests various misfolding pathways for human PrP in response to mutation.

Quantum mechanical studies on model alpha-pleated sheets

Pauling and Corey proposed a pleated-sheet configuration, now called alpha-sheet, as one of the protein secondary structures in addition to alpha-helix and beta-sheet. Recently, it has been suggested that alpha-sheet is a common feature of amyloidogenic intermediates. We have investigated the stability of antiparallel beta-sheet and two conformations of alpha-sheet in solution phase using the density functional theoretical method. The peptides are modeled as two-strand acetyl-(Ala)(2)-N-methylamine. Using stages of geometry optimization and single point energy calculation at B3LYP/cc-pVTZ//B3LYP/6-31G* level and including zero-point energies, thermal, and entropic contribution, we have found that beta-sheet is the most stable conformation, while the alpha-sheet proposed by Pauling and Corey has 13.6 kcal/mol higher free energy than the beta-sheet. The alpha-sheet that resembles the structure observed in molecular dynamics simulations of amyloidogenic proteins at low pH becomes distorted after stages of geometry optimization in solution. Whether the alpha-sheets with longer chains would be increasingly favorable in water relative to the increase in internal energy of the chain needs further investigation. Different from the quantum mechanics results, AMBER parm94 force field gives small difference in solution phase energy between alpha-sheet and beta-sheet. The predicted amide I IR spectra of alpha-sheet shows the main band at higher frequency than beta-sheet.

Temperature dependence of the flexibility of thermophilic and mesophilic flavoenzymes of the nitroreductase fold

A widely held hypothesis regarding the thermostability of thermophilic proteins states asserts that, at any given temperature, thermophilic proteins are more rigid than their mesophilic counterparts. Many experimental and computational studies have addressed this question with conflicting results. Here, we compare two homologous enzymes, one mesophilic (Escherichia coli FMN-dependent nitroreductase; NTR) and one thermophilic (Thermus thermophilus NADH oxidase; NOX), by multiple molecular dynamics simulations at temperatures from 5 to 100 degrees C. We find that the global rigidity/flexibility of the two proteins, assessed by a variety of metrics, is similar on the time scale of our simulations. However, the thermophilic enzyme retains its native conformation to a much greater degree at high temperature than does the mesophilic enzyme, both globally and within the active site. The simulations identify the helix F-helix G 'arm' as the region with the greatest difference in loss of native contacts between the two proteins with increasing temperature. In particular, a network of electrostatic interactions holds helix F to the body of the protein in the thermophilic protein, and this network is absent in the mesophilic counterpart.

Secondary structure of Huntingtin amino-terminal region

Huntington's disease is a genetic neurodegenerative disorder resulting from polyglutamine (polyQ) expansion (>36Q) within the first exon of Huntingtin (Htt) protein. We applied X-ray crystallography to determine the secondary structure of the first exon (EX1) of Htt17Q. The structure of Htt17Q-EX1 consists of an amino-terminal alpha helix, poly17Q region, and polyproline helix formed by the proline-rich region. The poly17Q region adopts multiple conformations in the structure, including alpha helix, random coil, and extended loop. The conformation of the poly17Q region is influenced by the conformation of neighboring protein regions, demonstrating the importance of the native protein context. We propose that the conformational flexibility of the polyQ region observed in our structure is a common characteristic of many amyloidogenic proteins. We further propose that the pathogenic polyQ expansion in the Htt protein increases the length of the random coil, which promotes aggregation and facilitates abnormal interactions with other proteins in cells.

Structural changes to monomeric CuZn superoxide dismutase caused by the familial amyotrophic lateral sclerosis-associated mutation A4V

Amyotrophic lateral sclerosis (ALS) is a progressive motor neuron degenerative disease, and the inherited form, familial ALS (fALS), has been linked to over 100 different point mutations scattered throughout the Cu-Zn superoxide dismutase protein (SOD1). The disease is likely due to a toxic gain of function caused by the misfolding, oligomerization, and eventual aggregation of mutant SOD1, but it is not yet understood how the structurally diverse mutations result in a common disease phenotype. The behavior of the apo-monomer fALS-associated mutant protein A4V was explored using molecular-dynamics simulations to elucidate characteristic structural changes to the protein that may allow the mutant form to improperly associate with other monomer subunits. Simulations showed that the mutant protein is less stable than the WT protein overall, with shifts in residue-residue contacts that lead to destabilization of the dimer and metal-binding sites, and stabilization of nonnative contacts that leads to a misfolded state. These findings provide a unifying explanation for disparate experimental observations, allow a better understanding of alterations of residue contacts that accompany loss of SOD1 structural integrity, and suggest sites where compensatory changes may stabilize the mutant structure.

A hotspot of inactivation: The A22S and V108M polymorphisms individually destabilize the active site structure of catechol O-methyltransferase

Human catechol O-methyltransferase (COMT) contains three common polymorphisms (A22S, A52T, and V108M), two of which (A22S and V108M) render the protein susceptible to deactivation by temperature or oxidation. We have performed multiple molecular dynamics simulations of the wild-type, A22S, A52T, and V108M COMT proteins to explore the structural consequences of these mutations. In total, we have amassed more than 1.4 micros of simulation time, representing the largest set of simulations detailing the effects of polymorphisms on a protein system to date. The A52T mutation had no significant effect on COMT structure in accord with experiment, thereby serving as a good negative control for the simulation set. Residues 22 (alpha2) and 108 (alpha5) interact with each other throughout the simulations and are located in a polymorphic hotspot approximately 20 A from the active site. Introduction of either the larger Ser (22) or Met (108) tightens this interaction, pulling alpha2 and alpha5 toward each other and away from the protein core. The V108M polymorphism rearranges active-site residues in alpha5, beta3, and alpha6, increasing the S-adenosylmethionine site solvent exposure. The A22S mutation reorients alpha2, moving critical catechol-binding residues away from the substrate-binding pocket. The A22S and V108M polymorphisms evolved independently in Northern European and Asian populations. While the decreased activities of both A22S and V108M COMT are associated with an increased risk for schizophrenia, the V108M-induced destabilization is also linked with improved cognitive function. These results suggest that polymorphisms within this hotspot may have evolved to regulate COMT activity and that heterozygosity for either mutation may be advantageous.

Potential aggregation prone regions in biotherapeutics: A survey of commercial monoclonal antibodies

Aggregation of a biotherapeutic is of significant concern and judicious process and formulation development is required to minimize aggregate levels in the final product. Aggregation of a protein in solution is driven by intrinsic and extrinsic factors. In this work we have focused on aggregation as an intrinsic property of the molecule. We have studied the sequences and Fab structures of commercial and non-commercial antibody sequences for their vulnerability towards aggregation by using sequence based computational tools to identify potential aggregation-prone motifs or regions. The mAbs in our dataset contain 2 to 8 aggregation-prone motifs per heavy and light chain pair. Some of these motifs are located in variable domains, primarily in CDRs. Most aggregation-prone motifs are rich in beta branched aliphatic and aromatic residues. Hydroxyl-containing Ser/Thr residues are also found in several aggregation-prone motifs while charged residues are rare. The motifs found in light chain CDR3 are glutamine (Q)/asparagine (N) rich. These motifs are similar to the reported aggregation promoting regions found in prion and amyloidogenic proteins that are also rich in Q/N, aliphatic and aromatic residues. The implication is that one possible mechanism for aggregation of mAbs may be through formation of cross-beta structures and fibrils. Mapping on the available Fab-receptor/antigen complex structures reveals that these motifs in CDRs might also contribute significantly towards receptor/antigen binding. Our analysis identifies the opportunity and tools for simultaneous optimization of the therapeutic protein sequence for potency and specificity while reducing vulnerability towards aggregation.

Biophysical characterization of intrinsically disordered proteins

The challenges associated with the structural characterization of disordered proteins have resulted in the application of a host of biophysical methods to such systems. NMR spectroscopy is perhaps the most readily suited technique for providing high-resolution structural information on disordered protein states in solution. Optical methods, solid state NMR, ESR and X-ray scattering can also provide valuable information regarding the ensemble of conformations sampled by disordered states. Finally, computational studies have begun to assume an increasingly important role in interpreting and extending the impact of experimental data obtained for such systems. This article discusses recent advances in the applications of these methods to intrinsically disordered proteins.

Folding of polyglutamine chains

Long polyglutamine chains have been associated with a number of neurodegenerative diseases. These include Huntington's disease, where expanded polyglutamine (PolyQ) sequences longer than 36 residues are correlated with the onset of symptoms. In this paper we study the folding pathway of a 54-residue PolyQ chain into a beta-helical structure. Transition path sampling Monte Carlo simulations are used to generate unbiased reactive pathways between unfolded configurations and the folded beta-helical structure of the polyglutamine chain. The folding process is examined in both explicit water and an implicit solvent. Both models reveal that the formation of a few critical contacts is necessary and sufficient for the molecule to fold. Once the primary contacts are formed, the fate of the protein is sealed and it is largely committed to fold. We find that, consistent with emerging hypotheses about PolyQ aggregation, a stable beta-helical structure could serve as the nucleus for subsequent polymerization of amyloid fibrils. Our results indicate that PolyQ sequences shorter than 36 residues cannot form that nucleus, and it is also shown that specific mutations inferred from an analysis of the simulated folding pathway exacerbate its stability.

Amide I infrared spectral features characteristic of some untypical conformations appearing in the structures suggested for amyloids

Amide I infrared (IR) spectral features are studied, by using the density functional theoretical method, for two untypical (but possibly rather prevalent) structures inspired from those recently suggested for amyloids: a structure consisting of loop regions in the (alpha L, alpha R) conformation stacked to form an alpha-sheet, and a structure involving some main-chain peptide groups (of any residues) and some side-chain amide groups of glutamine and asparagine residues closely located with each other. The amide I vibrational (off-diagonal) coupling constants are examined by extracting them from the calculated Cartesian-based force constants with the average partial vector method and by comparing them with those estimated on the basis of the transition dipole coupling mechanism. It is suggested that the amide I IR band characteristic of the alpha-sheet conformation in dry environment (without hydrogen bonding to solvent water molecules) is located in a high-frequency region (approximately >1670 cm(-1), somewhat higher than that of alpha-helix), because of the dependence of the diagonal (uncoupled) frequency and the off-diagonal coupling constant on the Phi and Psi dihedral angles. It is also shown that the amide I vibrations of the closely located peptide and amide groups are strongly coupled through-space with each other, and in the presence of this type of strong vibrational coupling, a noticeable change in the IR intensity upon (13)C=O substitution may occur even for a mode that arises mainly from an unsubstituted group and is not much shifted in frequency. The meaning of these results in the interpretation of observed amide I spectral profiles, especially the possible usefulness of IR spectroscopic measurements for detecting those untypical structures in the process of amyloid formation, is also discussed.

Insights into structure, stability, and toxicity of monomeric and aggregated polyglutamine models from molecular dynamics simulations

Nine genetically inherited neurodegenerative diseases are linked to abnormal expansions of a polyglutamine (polyQ) encoding region. Over the years, several structural models for polyQ regions have been proposed and confuted. The cross-beta-spine steric zipper motif, identified recently for the GNNQQNY peptide, represents an attractive model for amyloid fibers formed by polyQ fragments. Here we report a detailed molecular dynamics investigation of polyQ models assembled by cross-beta-spine steric zipper motifs. Our simulations indicate clearly that these assemblies are very stable. Glutamine side chains contribute strongly to the overall stability of the models by fitting perfectly within the zipper. In contrast to GNNQQNY zipper motifs, hydrogen bonding interactions provide a significant contribution to the overall stability of polyQ models. Molecular dynamics simulations carried out on monomeric polyQ forms (composed by 40-60 residues) show clearly that they can also assume structures stabilized by steric zipper motifs. Based on these findings, we build monomeric polyQ models that can explain recent data on the toxicity exerted by these species. In a more general context, our data suggests that polyQ models with interdigitated side chains can provide a structural rationale to several literature experiments on polyQ formation, stability, and toxicity.

The role of the turn in beta-hairpin formation during WW domain folding

The folding of WW domains is rate limited by formation of a beta-hairpin comprising residues from strands 1 and 2. Residues in the turn of this hairpin have reported Phi-values for folding close to 1 and have been proposed to nucleate folding. High Phi-values do not necessarily imply that the energetics of formation are a driving force for initiating folding. We demonstrate by NMR studies and molecular dynamics simulations that the first turn of the hYAP, FBP28, and PIN1 WW domains is structurally dynamic and solvent exposed in the native and folding transition states. It is, therefore, unlikely that the formation of the beta-turn per se provides the energetic driving force for hairpin folding. It is more likely that the turn acts as an easily formed hinge that facilitates the formation of the hairpin; it is a nucleus as defined by the nucleation-condensation mechanism whereby a diffuse nucleus is stabilized by associated interactions.