Biomedical subjects

Joan-Emma Shea

Publications and source records attributed to Joan-Emma Shea.

17 recordsLinked to original sources

PubMed2006

The structure of the Alzheimer amyloid beta 10-35 peptide probed through replica-exchange molecular dynamics simulations in explicit solvent.

The conformational states sampled by the Alzheimer amyloid beta (10-35) (Abeta 10-35) peptide were probed using replica-exchange molecular dynamics (REMD) simulations in explicit solvent. The Abeta 10-35 peptide is a fragment of the full-length Abeta 40/42 peptide that possesses many of the amyloidogenic properties of its full-length counterpart. Under physiological temperature and pressure, our simulations reveal that the Abeta 10-35 peptide does not possess a single unique folded state. Rather, this peptide exists as a mixture of collapsed globular states that remain in rapid dynamic equilibrium with each other. This conformational ensemble is dominated by random coil and bend structures with insignificant presence of an alpha-helical or beta-sheet structure. The 3D structure of Abeta 10-35 is seen to be defined by a salt bridge formed between the side-chains of K28 and D23. This salt bridge is also observed in Abeta fibrils and our simulations suggest that monomeric conformations of Abeta 10-35 contain pre-folded structural motifs that promote rapid aggregation of this peptide.

Algorithms↗

PubMed2006

Folding landscapes of the Alzheimer amyloid-beta(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.

Alzheimer Disease↗

PubMed2006

Effects of solvent on the structure of the Alzheimer amyloid-beta(25-35) peptide.

The free energy landscape for folding of the Alzheimer's amyloid-beta(25-35) peptide is explored using replica exchange molecular dynamics in both pure water and in HFIP/water cosolvent. This amphiphilic peptide is a natural by-product of the Alzheimer's amyloid-beta(1-40) peptide and retains the toxicity of its full-length counterpart as well as the ability to aggregate into beta-sheet-rich fibrils. Our simulations reveal that the peptide preferentially populates a helical structure in apolar organic solvent, while in pure water, the peptide adopts collapsed coil conformations and to a lesser extent beta-hairpin conformations. The beta-hairpin is characterized by a type II' beta-turn involving residues G29 and A30 and two short beta-strands involving residues N27, K28, I31, and I32. The hairpin is stabilized by backbone hydrogen-bonding interactions between residues K28 and I31; S26 and G33; and by side-chain-to-side-chain interactions between N27 and I32. Implications regarding the mechanism of aggregation of this peptide into fibrils and the role of the environment in modulating secondary structure are discussed.

Alzheimer Disease↗

PubMed2006

Effects of surface tethering on protein folding mechanisms.

The folding mechanisms of proteins are increasingly being probed through single-molecule experiments in which the protein is immobilized on a surface. Nevertheless, a clear understanding of how the surface might affect folding, and whether or not it changes folding from its bulk behavior, is lacking. In this work, we use molecular dynamics simulations of a model beta-barrel protein tethered to a surface to systematically investigate how the surface impacts folding. In the bulk, this protein folds in a three-state manner through a compact intermediate state, and its transition state (TS) has a well formed hydrophobic core. Upon tethering, we find that folding rates and stability are impacted differently by the surface, with dependencies on both the length and location of the tether. Significant changes in folding times are observed for tether points that do not alter the folding temperature. Tethering also locally enhances the formation of structure for residues proximal to the tether point. We find that neither the folding mechanism nor the TS of this protein are altered if the tether is in a fully structured or completely unstructured region of the TS. By contrast, tethering in a partially structured region of the TS leads to dramatic changes. For one such tether point, the intermediate present in bulk folding is eliminated, leading to a two-state folding process with a heterogeneous, highly unstructured TS ensemble. These results have implications for both the design of single-molecule experiments and biotechnological applications of tethered proteins.

Models, Molecular↗

PubMed2006

Aggregation of polyalanine in a hydrophobic environment.

The dimerization of polyalanine peptides in a hydrophobic environment was explored using replica exchange molecular dynamics simulations. A nonpolar solvent (cyclohexane) was used to mimic, among other hydrophobic environments, the hydrophobic interior of a membrane in which the peptides are fully embedded. Our simulations reveal that while the polyalanine monomer preferentially adopts a beta-hairpin conformation, dimeric phases exist in an equilibrium between random coil, alpha-helical, beta-sheet, and beta-hairpin states. A thermodynamic characterization of the dimeric phases reveals that electric dipole-dipole interactions and optimal side-chain packing stabilize alpha-helical conformations, while hydrogen bond interactions favor beta-sheet conformations. Possible pathways leading to the formation of alpha-helical and beta-sheet dimers are discussed.

Alanine↗

PubMed2006

Elucidating amyloid beta-protein folding and assembly: A multidisciplinary approach.

Oligomeric, neurotoxic amyloid protein assemblies are thought to be causative agents in Alzheimer's and other neurodegenerative diseases. Development of oligomer-specific therapeutic agents requires a mechanistic understanding of the oligomerization process. This is a daunting task because amyloidogenic protein oligomers often are metastable and comprise structurally heterogeneous populations in equilibrium with monomers and fibrils. A single methodological approach cannot elucidate the entire protein assembly process. An integrated multidisciplinary program is required. We discuss here the synergistic application of in hydro, in vacuo, and in silico methods to the study of the amyloid beta-protein, the key pathogenetic agent in Alzheimer's disease.

Amino Acid Sequence↗

PubMed2006

Amyloid beta-protein monomer structure: a computational and experimental study.

The structural properties of the Abeta42 peptide, a main constituent of the amyloid plaques formed in Alzheimer's disease, were investigated through a combination of ion-mobility mass spectrometry and theoretical modeling. Replica exchange molecular dynamics simulations using a fully atomic description of the peptide and implicit water solvent were performed on the -3 charge state of the peptide, its preferred state under experimental conditions. Equilibrated structures at 300 K were clustered into three distinct families with similar structural features within a family and with significant root mean square deviations between families. An analysis of secondary structure indicates the Abeta42 peptide conformations are dominated by loops and turns but show some helical structure in the C-terminal hydrophobic tail. A second calculation on Abeta42 in a solvent-free environment yields compact structures turned "inside out" from the solution structures (hydrophobic parts on the outside, polar parts on the inside). Ion mobility experiments on the Abeta42 -3 charge state electrosprayed from solution yield a bimodal arrival time distribution. This distribution can be quantitatively fit using cross-sections from dehydrated forms of the three families of calculated solution structures and the calculated solvent-free family of structures. Implications of the calculations on the early stages of aggregation of Abeta42 are discussed.

Amyloid beta-Peptides↗

PubMed2006

Structure of the 21-30 fragment of amyloid beta-protein.

Folding and self-assembly of the 42-residue amyloid beta-protein (Abeta) are linked to Alzheimer's disease (AD). The 21-30 region of Abeta, Abeta(21-30), is resistant to proteolysis and is believed to nucleate the folding of full-length Abeta. The conformational space accessible to the Abeta(21-30) peptide is investigated by using replica exchange molecular dynamics simulations in explicit solvent. Conformations belonging to the global free energy minimum (the "native" state) from simulation are in good agreement with reported NMR structures. These conformations possess a bend motif spanning the central residues V24-K28. This bend is stabilized by a network of hydrogen bonds involving the side chain of residue D23 and the amide hydrogens of adjacent residues G25, S26, N27, and K28, as well as by a salt bridge formed between side chains of K28 and E22. The non-native states of this peptide are compact and retain a native-like bend topology. The persistence of structure in the denatured state may account for the resistance of this peptide to protease degradation and aggregation, even at elevated temperatures.

Amyloid beta-Peptides↗

PubMed2005

The influence of different treatments of electrostatic interactions on the thermodynamics of folding of peptides.

Replica exchange molecular dynamics simulations were performed to investigate the effects of different electrostatic treatments on the structure and thermodynamics of a small beta-hairpin forming peptide. Three different electrostatic schemes were considered: regular cutoffs, generalized reaction field (GRF), and particle mesh Ewald (PME), with the peptide modeled using OPLS/AA all-atom force field with explicit TIP3P water. Both the GRF and PME methods yielded results consistent with experiment, with free energy surfaces displaying a single minimum corresponding to the native beta-hairpin structure. In contrast, use of straight cutoffs led to the population of an additional local minimum corresponding to nonhairpin conformations that compete with the formation of the native beta-hairpin at low temperatures. This extra minimum would not be apparent in conventional constant-temperature molecular dynamics simulations run for a few nanoseconds. This result points to the critical need of careful sampling of conformational space to assess the quality of different numerical treatments of long-range forces. While differences emerged in the nature of the unfolded states populated using PME and GRF approaches, simulations on the beta-hairpin forming peptide and on two additional control peptides indicate that the GRF treatment of electrostatics offers a satisfactory compromise between accuracy and computational speed for the identification of low-energy conformations. A GRF-based approach emerges as a viable means for treating larger biological systems that would be prohibitively costly to simulate using PME methods.

Computer Simulation↗

PubMed2005

Effects of frustration, confinement, and surface interactions on the dimerization of an off-lattice beta-barrel protein.

We study the effects of confinement, sequence frustration, and surface interactions on the thermodynamics of dimerization of an off-lattice minimalist beta-barrel protein using replica exchange molecular dynamics. We vary the degree of frustration of the protein by tuning the specificity of the hydrophobic interactions and investigate dimerization in confining spheres of different radii. We also investigate surface effects by tethering the first residue of one of the proteins to a uniformly repulsive surface. We find that increasing the confinement and frustration stabilize the dimer, while adding a repulsive surface decreases its stability. Different ensembles of structures, including properly dimerized and various partially dimerized states, are observed at the association transition temperature T(a), depending on the amount of frustration and whether a surface is present. The presence of a surface is predicted to alter the morphology of larger aggregates formed from partially unfolded dimeric conformations.

Chemistry, Physical↗

PubMed2005

Reconstruction of the src-SH3 protein domain transition state ensemble using multiscale molecular dynamics simulations.

We use an integrated computational approach to reconstruct accurately the transition state ensemble (TSE) for folding of the src-SH3 protein domain. We first identify putative TSE conformations from free energy surfaces generated by importance sampling molecular dynamics for a fully atomic, solvated model of the src-SH3 protein domain. These putative TSE conformations are then subjected to a folding analysis using a coarse-grained representation of the protein and rapid discrete molecular dynamics simulations. Those conformations that fold to the native conformation with a probability (P(fold)) of approximately 0.5, constitute the true transition state. Approximately 20% of the putative TSE structures were found to have a P(fold) near 0.5, indicating that, although correct TSE conformations are populated at the free energy barrier, there is a critical need to refine this ensemble. Our simulations indicate that the true TSE conformations are compact, with a well-defined central beta sheet, in good agreement with previous experimental and theoretical studies. A structured central beta sheet was found to be present in a number of pre-TSE conformations, however, indicating that this element, although required in the transition state, does not define it uniquely. An additional tight cluster of contacts between highly conserved residues belonging to the diverging turn and second beta-sheet of the protein emerged as being critical elements of the folding nucleus. A number of commonly used order parameters to identify the transition state for folding were investigated, with the number of native Cbeta contacts displaying the most satisfactory correlation with P(fold) values.

Models, Molecular↗

PubMed2005

Amyloid beta-protein: monomer structure and early aggregation states of Abeta42 and its Pro19 alloform.

The amyloid beta-protein (Abeta) is a seminal neuropathic agent in Alzheimer's disease (AD). Recent evidence points to soluble Abeta oligomers as the probable neurotoxic species. Among the naturally occurring Abeta peptides, the 42-residue form Abeta42 is linked particularly strongly with AD, even though it is produced at approximately 10% of the levels of the more abundant 40-residue form Abeta40. Here, we apply mass spectrometry and ion mobility to the study of Abeta42 and its Pro19 alloform. The Phe19 --> Pro19 substitution blocks fibril formation by [Pro19]Abeta42. Evidence indicates that solution-like structures of Abeta monomers are electrosprayed and characterized. Unfiltered solutions of Abeta42 produce only monomers and large oligomers, whereas [Pro19]Abeta42 solutions produce abundant monomers, dimers, trimers, and tetramers but no large oligomers. When passed through a 10,000 amu filter and immediately sampled, Abeta42 solutions produce monomers, dimers, tetramers, hexamers, and an aggregate of two hexamers that may be the first step in protofibril formation. These results are consistent with recently published photochemical cross-linking data and lend support to recent aggregation mechanisms proposed by Bitan, Teplow, and co-workers [J. Biol. Chem. 2003, 278, 34882-34889].

Amino Acid Sequence↗

PubMed2005

The physics of the interactions governing folding and association of proteins.

The review discusses the molecular origins of the forces and free energies that determine several things about proteins, and how experiment and theory reveal this information. The first subject is the stability of the folded, native structures. The second is the range of molecular mechanisms by which proteins find their way to those folded structures in laboratory environments. The third is the much more complex problem of how folding occurs in the cellular environment. This topic includes a discussion of crowding and of the roles of chaperone molecules. The review concludes with a discussion of protein aggregation and fibril formation and of misfolding and therapies associated with it.

Amino Acids↗

PubMed2004

Temperature dependence of the free energy landscape of the src-SH3 protein domain.

The temperature dependence of the free energy landscape of the src-SH3 protein domain is investigated through fully atomic simulations in explicit solvent. Simulations are performed above and below the folding transition temperature, enabling an analysis of both protein folding and unfolding. The transition state for folding and unfolding, identified from the free energy surfaces, is found to be very similar, with structure in the central hydrophobic sheet and little structure throughout the rest of the protein. This is a result of a polarized folding (unfolding) mechanism involving early formation (late loss) of the central hydrophobic sheet at the transition state. Unfolding simulations map qualitatively well onto low-temperature free energy surfaces but appear, however, to miss important features observed in folding simulations. In particular, details of the folding mechanism involving the opening and closing of the hydrophobic core are not captured by unfolding simulations performed under strongly denaturing conditions. In addition, free energy surfaces at high temperatures do not display a desolvation barrier found at lower temperatures, involving the expulsion of water molecules from the hydrophobic core.

Computer Simulation↗

PubMed2004

Self-assembly of peptides into a beta-barrel motif.

We report the results of a study of the self-assembly of four minimalist peptide strands with a native beta-barrel structure. Using a soft-well potential to mimic cellular crowding, molecular dynamics simulations were performed in confining spheres of varying radii. By utilizing a previously introduced scaling factor lambda for the non-native hydrophobic interactions (0<lambda<1), we were able to study models with varying degrees of frustration. Both the thermodynamics and kinetics of a Go-like model (lambda=0) and a highly frustrated model (lambda=0.9) were studied. Additionally, we used an extrapolation technique to investigate the thermodynamics of assembly at intermediate values of lambda. As in our earlier work [J. Chem. Phys. 118, 8106 (2003)] on a connected Go-like model beta-barrel protein, we find that the stability of the assembled protein increases with decreasing sphere size, and that larger confining spheres result in increased assembly times. Additionally, the lambda=0 model seems to undergo distinct phase transitions during the assembly process. In contrast, the more frustrated model (lambda=0.9) appears to undergo a glasslike transition at temperatures comparable to the assembly temperature of the Go model, and that this transition is relatively nonspecific. Our results suggest the assembly process is dependent on both sequence and environment, with implications for the formation of misassembled aggregates.

Computer Simulation↗

PubMed2003

Posttransition state desolvation of the hydrophobic core of the src-SH3 protein domain.

The folding thermodynamics of the src-SH3 protein domain were characterized under refolding conditions through biased fully atomic molecular dynamics simulations with explicit solvent. The calculated free energy surfaces along several reaction coordinates revealed two barriers. The first, larger barrier was identified as the transition state barrier for folding, associated with the formation of the first hydrophobic sheet of the protein. phi values calculated from structures residing at the transition state barrier agree well with experimental phi values. The microscopic information obtained from our simulations allowed us to unambiguously assign intermediate phi values as the result of multiple folding pathways. The second, smaller barrier occurs later in the folding process and is associated with the cooperative expulsion of water molecules between the hydrophobic sheets of the protein. This posttransition state desolvation barrier cannot be observed through traditional folding experiments, but is found to be critical to the correct packing of the hydrophobic core in the final stages of folding. Hydrogen exchange and NMR experiments are suggested to probe this barrier.

Binding Sites↗

PubMed2002

Probing the folding free energy landscape of the Src-SH3 protein domain.

The mechanism and thermodynamics of folding of the Src homology 3 (SH3) protein domain are characterized at an atomic level through molecular dynamics with importance sampling. This methodology enables the construction of the folding free energy landscape of the protein as a function of representative reaction coordinates. We observe that folding proceeds in a downhill manner under native conditions, with early compaction and structure formation in the hydrophobic sheet consisting of the three central beta strands of the protein. This state bears considerable resemblance to the experimentally determined transition state for folding. Folding proceeds further with the formation of the second hydrophobic sheet consisting of the terminal strands and the RT loop. The final stages of folding appear to involve the formation of the hydrophobic core through the expulsion of water molecules bridging the two hydrophobic sheets. This work sheds new light on the complementary roles of sequence and topology in governing the folding mechanism of small proteins and provides further support for the role of water in facilitating the late stages in folding.

Animals↗