Thermodynamics of beta-amyloid fibril formation

Integrating Protein Interaction Surface Prediction with a Fragment-Based Drug Design: Automatic Design of New Leads with Fragments on Energy Surfaces

Protein-protein interactions (PPIs) have emerged in the past years as significant pharmacological targets in the development of new therapeutics due to their key roles in determining pathological pathways. Herein, we present fragments on energy surfaces, a simple and general design strategy that integrates the analysis of the dynamic and energetic signatures of proteins to unveil the substructures involved in PPIs, with docking, selection, and combination of drug-like fragments to generate new PPI inhibitor candidates. Specifically, structural representatives of the target protein are used as inputs for the blind physics-based prediction of potential protein interaction surfaces using the matrix of low coupling energy decomposition method. The predicted interaction surfaces are subdivided into overlapping windows that are used as templates to direct the docking and combination of fragments representative of moieties typically found in active drugs. This protocol is then applied and validated using structurally diverse, important PPI targets as test systems. We demonstrate that our approach facilitates the exploration of the molecular diversity space of potential ligands, with no requirement of prior information on the location and properties of interaction surfaces or on the structures of potential lead compounds. Importantly, the hit molecules that emerge from our ab initio design share high chemical similarity with experimentally tested active PPI inhibitors. We propose that the protocol we describe here represents a valuable means of generating initial leads against difficult targets for further development and refinement.

The Subtle Trade-Off between Evolutionary and Energetic Constraints in Protein-Protein Interactions

Life machinery, although overwhelmingly complex, is rooted on a rather limited number of molecular processes. One of the most important is protein-protein interaction. Metabolic regulation, protein folding control, and cellular motility are examples of processes based on the fine-tuned interaction of several protein partners. The region on the protein surface devoted to the recognition of a specific partner is essential for the function of the protein and is, therefore, likely to be conserved during evolution. On the other hand, the physical chemistry of amino acids underlies the mechanism of interactions. Both evolutionary and energetic constraints can then be used to build scoring functions capable of recognizing interaction sites. Our working hypothesis is that residues within the interaction interface tend at the same time to be evolutionarily conserved (to preserve their function) and to provide little contribution to the internal stabilization of the structure of their cognate protein, to facilitate conformational adaptation to the partner. Here, we show that for some classes of protein partners (for example, those involved in signal transduction and in enzymes) evolutionary constraints play the key role in defining the interaction surface. In contrast, energetic constraints emerge as more important in protein partners involved in immune response, in inhibitor proteins, and in structural proteins. Our results indicate that a general-purpose scoring function for protein-protein interaction should not be agnostic of the biological function of the partners.

Connecting macroscopic observables and microscopic assembly events in amyloid formation using coarse grained simulations

The pre-fibrillar stages of amyloid formation have been implicated in cellular toxicity, but have proved to be challenging to study directly in experiments and simulations. Rational strategies to suppress the formation of toxic amyloid oligomers require a better understanding of the mechanisms by which they are generated. We report Dynamical Monte Carlo simulations that allow us to study the early stages of amyloid formation. We use a generic, coarse-grained model of an amyloidogenic peptide that has two internal states: the first one representing the soluble random coil structure and the second one the -sheet conformation. We find that this system exhibits a propensity towards fibrillar self-assembly following the formation of a critical nucleus. Our calculations establish connections between the early nucleation events and the kinetic information available in the later stages of the aggregation process that are commonly probed in experiments. We analyze the kinetic behaviour in our simulations within the framework of the theory of classical nucleated polymerisation, and are able to connect the structural events at the early stages in amyloid growth with the resulting macroscopic observables such as the effective nucleus size. Furthermore, the free-energy landscapes that emerge from these simulations allow us to identify pertinent properties of the monomeric state that could be targeted to suppress oligomer formation.

A number of normally soluble proteins can form amyloid structures in a process associated with neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Mature amyloid structures consist of large fibrils containing thousands of individual proteins aggregated into linear nanostructures; there is increasing evidence, however, that the toxic species responsible for neurodegeneration are not the mature fibrils themselves but rather lower molecular weight precursors commonly known as amyloid oligomers. Unfortunately, these early oligomers are commonly thermodynamically unstable and of nanometer scale dimensions, factors which make them highly challenging to probe in detail in experiments. We have used computer simulations of a model inspired by Alzheimer's Abeta peptide to investigate the early stages of protein aggregation. The results that we obtain were shown to fit Oosawa's polymerization theory, a finding which allows us to provide a connection between the microscopic molecular parameters and macroscopic growth. One crucial parameter is size of the nucleus, i.e. the basic oligomer existing at origin of the formation of each fiber. We have revealed a path for the formation of this nucleus and validate its size by several methods. Our results provide fundamental information for influencing the early stages of amyloid formation in a rational manner.

What drives amyloid molecules to assemble into oligomers and fibrils?

We develop a theory for three states of equilibrium of amyloid peptides: the monomer, oligomer, and fibril. We assume that the oligomeric state is a disordered micellelike collection of a few peptide chains held together loosely by hydrophobic interactions into a spherical hydrophobic core. We assume that fibrillar amyloid chains are aligned and further stabilized by steric zipper interactions-hydrogen bonding, steric packing, and specific hydrophobic side-chain contacts. The model makes a broad set of predictions that are consistent with experimental results: 1), Similar to surfactant micellization, amyloid oligomerization should increase with peptide concentration in solution. 2), The onset of fibrillization limits the concentration of oligomers in the solution. 3), The extent of Aβ fibrillization increases with peptide concentration. 4), The predicted average fibril length versus monomer concentration agrees with data on α-synuclein. 5), Full fibril length distributions agree with data on α-synuclein. 6), Denaturants should melt out fibrils. And finally, 7), added salt should stabilize fibrils by reducing repulsions between amyloid peptide chains. It is of interest that small changes in solvent conditions can tip the equilibrium balance between oligomer and fibril and cause large changes in rates through effects on the transition-state barrier. This model may provide useful insights into the physical processes underlying amyloid diseases.

Self-assembly of protein amyloids: a competition between amorphous and ordered aggregation

Protein aggregation in the form of amyloid fibrils has important biological and technological implications. Although the self-assembly process is highly efficient, aggregates not in the fibrillar form would also occur and it is important to include these disordered species when discussing the thermodynamic equilibrium behavior of the system. Here, we initiate such a task by considering a mixture of monomeric proteins and the corresponding aggregates in the disordered form (micelles) and in the fibrillar form (amyloid fibrils). Starting with a model on the respective binding free energies for these species, we calculate their concentrations at thermal equilibrium. We then discuss how the incorporation of the disordered structure furthers our understanding on the various amyloid promoting factors observed empirically, and on the kinetics of fibrilization.

A universal pathway for amyloid nucleus and precursor formation for insulin

To help identify the etiological agents for amyloid-related diseases, attention is focused here on the fibrillar precursors, also called oligomers and protofibrils, and on modeling the reaction kinetics of the formation of the amyloid nucleus. Insulin is a favored model for amyloid formation, not only because amyloidosis can be a problem in diabetes, but also because aggregation and fibrillation causes problems during production, storage, and delivery. Small angle neutron scattering (SANS) is used to measure the temporal formation of insulin oligomers in H(2)O- and D(2)O-based solvents and obtain consistent evidence of the composition of the insulin nucleus that comprised three dimers or six monomers similar to that recently proposed in the literature. A simple molecular structural model that describes the growth of oligomers under a wide range of environmental conditions is proposed. The model first involves lengthening or end-on-end association of dimers to form three-dimer nuclei, and then exhibits broadening or side-on-side association of nuclei. Using different additives to demonstrate their influence on the kinetics of oligomer formation, we showed that, although the time required to form the nucleus was dependent on a specific system, they all followed a universal pathway for nucleus and precursor formation. The methods and analyses presented here provide the first experimental molecular size description of the details of amyloid nucleus formation and subsequent propagation to fibril precursors independent of kinetics.

Finite size effects in simulations of protein aggregation

It is becoming increasingly clear that the soluble protofibrillar species that proceed amyloid fibril formation are associated with a range of neurodegenerative disorders such as Alzheimer's and Parkinson diseases. Computer simulations of the processes that lead to the formation of these oligomeric species are starting to make significant contributions to our understanding of the determinants of protein aggregation. We simulate different systems at constant concentration but with a different number of peptides and we study the how the finite number of proteins affects the underlying free energy of the system and therefore the relative stability of the species involved in the process. If not taken into account, this finite size effect can undermine the validity of theoretical predictions regarding the relative stability of the species involved and the rates of conversion from one to the other. We discuss the reasons that give rise to this finite size effect form both a probabilistic and energy fluctuations point of view and also how this problem can be dealt by a finite size scaling analysis.

Solutes, but not cells, drain from the brain parenchyma along basement membranes of capillaries and arteries: significance for cerebral amyloid angiopathy and neuroimmunology

Elimination of interstitial fluid and solutes plays a role in homeostasis in the brain, but the pathways are unclear. Previous work suggests that interstitial fluid drains along the walls of arteries.

AIMS

to define the pathways within the walls of capillaries and arteries for drainage of fluid and solutes out of the brain.

METHODS

Fluorescent soluble tracers, dextran (3 kDa) and ovalbumin (40 kDa), and particulate fluospheres (0.02 microm and 1.0 microm in diameter) were injected into the corpus striatum of mice. Brains were examined from 5 min to 7 days by immunocytochemistry and confocal microscopy.

RESULTS

soluble tracers initially spread diffusely through brain parenchyma and then drain out of the brain along basement membranes of capillaries and arteries. Some tracer is takenf up by vascular smooth muscle cells and by perivascular macrophages. No perivascular drainage was observed when dextran was injected into mouse brains following cardiac arrest. Fluospheres expand perivascular spaces between vessel walls and surrounding brain, are ingested by perivascular macrophages but do not appear to leave the brain even following an inflammatory challenge with lipopolysaccharide or kainate.

CONCLUSIONS

capillary and artery basement membranes act as 'lymphatics of the brain' for drainage of fluid and solutes; such drainage appears to require continued cardiac output as it ceases following cardiac arrest. This drainage pathway does not permit migration of cells from brain parenchyma to the periphery. Amyloid-beta is deposited in basement membrane drainage pathways in cerebral amyloid angiopathy, and may impede elimination of amyloid-beta and interstitial fluid from the brain in Alzheimer's disease. Soluble antigens, but not cells, drain from the brain by perivascular pathways. This atypical pattern of drainage may contribute to partial immune privilege of the brain and play a role in neuroimmunological diseases such as multiple sclerosis.

Modeling amyloid fibril formation: a free-energy approach

Amyloid fibrils are structures consisting of many proteins with a well-defined conformation. The formation of these fibrils has been the subject of intense research, largely due to their connection to several diseases. We focus here on the computational studies and discuss these from a free-energy point of view. The fibrillogenic properties of many proteins can be predicted and understood by taking the relevant free energies into account in an appropriate way. This is because both the equilibrium and the kinetic properties of the protein system depend on its free-energy landscape. Advanced simulation techniques can be used to understand the relationship between the free-energy landscape of a protein and its three-dimensional structure and propensity to form amyloid fibrils. We give an overview of existing simulation techniques that operate at a molecular level of detail and that are capable of generating relevant free-energy values. The free energies obtained with these methods can be inserted into a statistical-mechanical or kinetic framework to predict mean fibril properties on length scales and time scales that are inaccessible by molecular-scale simulation methods.

Understanding the Molecular Basis for the Inhibition of the Alzheimer's Aβ-Peptide Oligomerization by Human Serum Albumin Using Saturation Transfer Difference and Off-Resonance Relaxation NMR Spectroscopy

Human serum albumin (HSA) inhibits the formation of amyloid beta-peptide (Abeta) fibrils in human plasma. However, currently it is not known how HSA affects the formation of the highly toxic soluble diffusible oligomers that occur in the initial stages of Abeta fibrillization. We have therefore investigated by solution NMR the interaction of HSA with the Abeta(12-28) peptide, which has been previously shown to provide a reliable and stable model for the early prefibrillar oligomers as well as to contain key determinants for the recognition by albumin. For this purpose we propose a novel NMR approach based on the comparative analysis of Abeta in its inhibited and filtrated states monitored through both saturation transfer difference and recently developed nonselective off-resonance relaxation experiments. This combined NMR strategy reveals a mechanism for the oligomerization inhibitory function of HSA, according to which HSA targets preferentially the soluble oligomers of Abeta(12-28) rather than its monomeric state. Specifically, HSA caps the exposed hydrophobic patches located at the growing and/or transiently exposed sites of the Abeta oligomers, thereby blocking the addition of further monomers and the growth of the prefibrillar assemblies. The proposed model has implications not only for the pharmacological treatment of Alzheimer's disease specifically but also for the inhibition of oligomerization in amyloid-related diseases in general. In addition, the proposed NMR approach is expected to be useful for the investigation of the mechanism of action of other oligomerization inhibitors as well as of other amyloidogenic systems.

Hydrophobic cooperativity as a mechanism for amyloid nucleation

The kinetics of amyloid fibril formation are in most cases explained by classical nucleation theory, yet the mechanisms behind nucleation are not well understood. We show using molecular dynamics simulations that the hydrophobic cooperativity in the self-association of the model amyloidogenic peptide STVIYE is sufficient to allow for nucleation-dependent polymerization with a pentamer critical nucleus. The role of electrostatics was also investigated. Novel considerations of the electrostatic solvation energy using the Born-Onsager equation are put forth to rationalize the aggregation of charged peptides and provide new insight into the energetic differences between parallel and antiparallel beta-sheets. Together these results help explain the influence of molecular charge in the class of fibril-forming hexapeptides recently designed by Serrano and collaborators.

Molecular dynamics simulation suggests possible interaction patterns at early steps of beta2-microglobulin aggregation

Early events in aggregation of proteins are not easily accessible by experiments. In this work, we perform a 5-ns molecular dynamics simulation of an ensemble of 27 copies of beta(2)-microglobulin in explicit solvent. During the simulation, the formation of intermolecular contacts is observed. The simulation highlights the importance of apical residues and, in particular, of those at the N-terminus end of the molecule. The most frequently found pattern of interaction involves a head-to-head contact arrangement of molecules. Hydrophobic contacts appear to be important for the establishment of long-lived (on the simulation timescale) contacts. Although early events on the pathway to aggregation and fibril formation are not directly related to the end-state of the process, which is reached on a much longer timescale, simulation results are consistent with experimental data and in general with a parallel arrangement of intermolecular beta-strand pairs.

Folding Landscapes of the Alzheimer Amyloid-β(12-28) Peptide

The energy landscape for folding of the 12-28 fragment of the Alzheimer amyloid beta (Abeta) peptide is characterized using replica-exchange molecular dynamics simulations with an all-atom peptide model and explicit solvent. At physiological temperatures, the peptide exists mostly as a collapsed random coil, populating a small fraction (less than 10%) of hairpins with a beta-turn at position V18F19, with another 10% of hairpin-like conformations possessing a bend rather than a turn in the central VFFA positions. A small fraction of the populated states, approximately 14%, adopt polyproline II (PPII) conformations. Folding of the structured hairpin states proceeds through the assembly of two locally stable segments, VFFAE and EDVGS. The interactions stabilizing these locally folded structural motifs are in conflict with those stabilizing the global fold of A12-28, a signature of underlying residual frustration in this peptide. At increased temperature, the population of both beta-strand and PPII conformations diminishes in favor of beta-turn and random-coil states. On the basis of the conformational preferences of Abeta 12-28 monomers, two models for the molecular structure of amyloid fibrils formed by this peptide are proposed.

The Effect of Aβ Conformation on the Metal Affinity and Aggregation Mechanism Studied by Circular Dichroism Spectroscopy

The conformational change and associated aggregation of beta amyloid (Abeta) with or without metals is the main cause of Alzheimer's disease (AD). In order to further understand the effects of Abeta and its associated metals on the aggregation mechanism, the influence of Abeta conformation on the metal affinity and aggregation was investigated using circular dichroism (CD) spectroscopy. The Abeta conformation is dependent on pH and trifluoroethanol (TFE). The binding of metals to Abeta was found to be dependent on the Abeta conformation. The aggregation induced by Abeta itself or its associated metals is completely diminished for Abeta in 40% TFE. Only in 5% and 25% TFE can Abeta undergo an alpha-helix to beta-sheet aggregation, which involve a three-state mechanism for the metal-free state, and a two-state transition for the metal-bound state, respectively. The aggregation-inducing activity of metals is in the order, Cu2+ > Fe3+ > or = Al3+ > Zn2+.

Amyloid β-Peptide Oligomerization in Silico: Dimer and Trimer

Soluble oligomers of Alzheimer's amyloid beta protein (Abeta) may act as effectors of neurotoxicity in early stages of Alzheimer's disease. Detailed information about the structure of Abeta in atomistic level and the dynamics of assembly of monomeric Abeta into oligomeric structures is rather elusive. We have performed replica exchange molecular dynamics (REMD) simulations on the formation of the dimer and trimer of Abeta10-35 peptide. We have observed spontaneous formation of several basic structural units that may act as a template or an intermediate for further aggregation of Alzheimer's Abeta protein. Various conformers, including interlocking structures of experimentally known bend double beta strands, are identified.

Sequence Dependence of Amyloid Fibril Formation: Insights from Molecular Dynamics Simulations

The clarification of the physico-chemical determinants underlying amyloid deposition is critical for our understanding of misfolding diseases. With this purpose we have performed a systematic all-atom molecular dynamics (MD) study of a series of single point mutants of the de novo designed amyloidogenic peptide STVIIE. Sixteen different 50ns long simulations using explicit solvent have been carried out starting from four different conformations of a polymeric six-stranded beta-sheet. The simulations have provided evidence for the influence of a small number of site-specific hydrophobic interactions on the packing and stabilization of nascent aggregates, as well as the interplay between side-chain interactions and the net charge of the molecule on the strand arrangement of polymeric beta-sheets. This MD analysis has also shed light into the origin of the position dependence on mutation of beta-sheet polymerization that was found experimentally for this model system. Our results suggest that MD can be applied to detect critical positions for beta-sheet aggregation within a given amyloidogenic stretch. Studies similar to the one presented here can guide site-directed mutations or the design of drugs that specifically disrupt the key stabilizing interactions of beta-sheet aggregates.

Modeling the α-helix to β-hairpin transition mechanism and the formation of oligomeric aggregates of the fibrillogenic peptide A(12–28): insights from all-atom molecular dynamics simulations

In this paper, all-atom molecular dynamics simulations in explicit solvent are used to investigate the structural and dynamical determinants of the alpha-helical to beta-hairpin conformational transition of the 12-28 fragment from the full length Abeta Alzheimer's peptide. The transition from alpha-helical to beta-structure requires the peptide to populate intermediate beta-bend geometries in which several mainly hydrophobic interactions are partially formed. This is followed by the sudden collapse to ordered beta-hairpin structures and the simultaneous disruption of the hydrophobic side-chain interactions with a consequent increase in solvent exposure. The solvent exposure of hydrophobic side-chains belonging to a sequence of five consecutive residues in the beta-hairpin defines a possible starting point for the onset of the aggregation mechanisms. Several different conformations of model oligomeric (dimeric and tetrameric) aggregates are then investigated. These simulations show that while hydrophobic contacts are important to bring together different monomers with a beta-hairpin like conformation, more specific interactions such as hydrogen-bonding and coulombic interactions, should be considered necessary to provide further stabilization and ordering to the nascent fibrillar aggregates.

Replica exchange molecular dynamics simulations of amyloid peptide aggregation

The replica exchange molecular dynamics (REMD) approach is applied to four oligomeric peptide systems. At physiologically relevant temperature values REMD samples conformation space and aggregation transitions more efficiently than constant temperature molecular dynamics (CTMD). During the aggregation process the energetic and structural properties are essentially the same in REMD and CTMD. A condensation stage toward disordered aggregates precedes the beta-sheet formation. Two order parameters, borrowed from anisotropic fluid analysis, are used to monitor the aggregation process. The order parameters do not depend on the peptide sequence and length and therefore allow to compare the amyloidogenic propensity of different peptides

In silico assembly of Alzheimer's Abeta16-22 peptide into beta-sheets

Recent studies suggest that soluble oligomers of amyloid-forming peptides have toxic effects in cell cultures. In this study, the folding of three Alzheimer's A beta(16-22) peptides have been simulated with the activation-relaxation technique and a generic energy model. Starting from randomly chosen states, the predicted lowest energy structure superposes within 1 A rms deviation from its conformation within the fibrils. This antiparallel structure is found to be in equilibrium with several out-of-register antiparallel beta-sheets and mixed parallel-antiparallel beta-sheets, indicating that full structural order in the fibrils requires larger aggregates. Folding involves the formation of dimers followed by the addition of a monomer and proceeds through a generalized mechanism between disordered and native alignments of beta-strands.