Search PubMed⌕ Search

Biomedical subjects

Valerie Daggett

Publications and source records attributed to Valerie Daggett.

At least 19 recordsLinked to original sources
PubMed2006

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

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

Amyloid↗
PubMed2006

Direct observation of microscopic reversibility in single-molecule protein folding.

Both folded and unfolded conformations should be observed for a protein at its melting temperature (T(m)), where DeltaG between these states is zero. In an all-atom molecular dynamics simulation of chymotrypsin inhibitor 2 (CI2) at its experimental T(m), the protein rapidly loses its low-temperature native structure; it then unfolds before refolding to a stable, native-like conformation. The initial unfolding follows the unfolding pathway described previously for higher-temperature simulations: the hydrophobic core is disrupted, the beta-sheet pulls apart and the alpha-helix unravels. The unfolded state reached under these conditions maintains a kernel of structure in the form of a non-native hydrophobic cluster. Refolding simply reverses this path, the side-chain interactions shift, the helix refolds, and the native packing and hydrogen bonds are recovered. The end result of this refolding is not the initial crystal structure; it contains the proper topology and the majority of the native contacts, but the structure is expanded and the contacts are long. We believe this to be the native state at elevated temperature, and the change in volume and contact lengths is consistent with experimental studies of other native proteins at elevated temperature and the chemical denaturant equivalent of T(m).

Computer Simulation↗
PubMed2006

Exploring the energy landscape of protein folding using replica-exchange and conventional molecular dynamics simulations.

Two independent replica-exchange molecular dynamics (REMD) simulations with an explicit water model were performed of the Trp-cage mini-protein. In the first REMD simulation, the replicas started from the native conformation, while in the second they started from a nonnative conformation. Initially, the first simulation yielded results qualitatively similar to those of two previously published REMD simulations: the protein appeared to be over-stabilized, with the predicted melting temperature 50-150K higher than the experimental value of 315K. However, as the first REMD simulation progressed, the protein unfolded at all temperatures. In our second REMD simulation, which starts from a nonnative conformation, there was no evidence of significant folding. Transitions from the unfolded to the folded state did not occur on the timescale of these simulations, despite the expected improvement in sampling of REMD over conventional molecular dynamics (MD) simulations. The combined 1.42 micros of simulation time was insufficient for REMD simulations with different starting structures to converge. Conventional MD simulations at a range of temperatures were also performed. In contrast to REMD, the conventional MD simulations provide an estimate of Tm in good agreement with experiment. Furthermore, the conventional MD is a fraction of the cost of REMD and continuous, realistic pathways of the unfolding process at atomic resolution are obtained.

Protein Folding↗
PubMed2006

Phi-analysis at the experimental limits: mechanism of beta-hairpin formation.

The 37-residue Formin-binding protein, FBP28, is a canonical three-stranded beta-sheet WW domain. Because of its small size, it is so insensitive to chemical denaturation that it is barely possible to determine accurately a denaturation curve, as the transition spans 0-7 M guanidinium hydrochloride (GdmCl). It is also only marginally stable, with a free energy of denaturation of just 2.3 kcal/mol at 10 degrees Celsius so only small changes in energy upon mutation can be tolerated. But these properties and relaxation times for folding of 25 micros-400 micros conspire to allow the rapid acquisition of accurate and reproducible kinetic data for Phi-analysis using classical temperature-jump methods. The transition state for folding is highly polarized with some regions having Phi-values of 0 and others 1, as readily seen in chevron plots, with Phi-values of 0 having the refolding arms overlaying and those of 1 the unfolding arms superimposable. Good agreement is seen with transition state structures identified from independent molecular dynamics (MD) simulations at 60, 75, and 100 degrees Celsius, which allows us to explore further the details of the folding and unfolding pathway of FBP28. The first beta-turn is near native-like in the transition state for folding (experimental) and unfolding (MD and experiment). The simulations show that there are transient contacts between the aromatic side-chains of the beta-strands in the denatured state and that these interactions provide the driving force for folding of the first beta-hairpin of this three-stranded sheet. Only after the backbone hydrogen bonds are formed between beta1 and beta2 does a hydrogen bond form to stabilize the intervening turn, or the first beta-turn.

Amino Acid Sequence↗
PubMed2006

Importance of context in protein folding: secondary structural propensities versus tertiary contact-assisted secondary structure formation.

Molecular dynamics simulations can be used to reveal the detailed conformational behaviors of peptides and proteins. By comparing fragment and full-length protein simulations, we can investigate the role of each peptide segment in the folding process. Here, we take advantage of information regarding the helix formation process from our previous simulations of barnase and protein A as well as new simulations of four helical fragments from these proteins at three different temperatures, starting with both helical and extended structures. Segments with high helical propensity began the folding process by tethering the chain through side chain interactions involving either polar interactions, such as salt bridges, or hydrophobic staples. These tethers were frequently nonnative (i.e., not i --> i + 4 spacing) and provided a scaffold for other residues, thereby limiting the conformational search. The helical structure then propagated on both sides of the tether. Segments with low stability and propensity formed later in the folding process and utilized contacts with other portions of the protein when folding. These helices formed via a tertiary contact-assisted mechanism, primarily via hydrophobic contacts between residues distant in sequence. Thus, segments with different helical propensities appear to play different roles during protein folding. Furthermore, the active role of nonlocal side chains in helix formation highlights why we must move beyond simple hierarchical models of protein folding.

Amino Acid Sequence↗
PubMed2006

The folding pathway of spectrin R17 from experiment and simulation: using experimentally validated MD simulations to characterize States hinted at by experiment.

We present an experimental and computational analysis of the folding pathway of the 17th domain of chicken brain alpha-spectrin, R17. Wild-type R17 folds in a two-state manner and the chevron plot (plot of the logarithm of the observed rate constant against concentration of urea) shows essentially linear folding and unfolding arms. A number of mutant proteins, however, show a change in slope of the unfolding arm at high concentration of denaturant, hinting at complexity in the folding landscape. Through a combination of mutational studies and high temperature molecular dynamics simulations we show that the folding of R17 can be described by a model with two sequential transition states separated by an intermediate species. The rate limiting transition state for folding in water has been characterized both through experimental Phi-value analysis and by simulation. In contrast, a detailed analysis of the transition state predicted to dominate under highly denaturing conditions is only possible by simulation.

Amino Acid Sequence↗
PubMed2006

The 108M polymorph of human catechol O-methyltransferase is prone to deformation at physiological temperatures.

The human gene for catechol O-methyltransferase (COMT) contains a common polymorphism that results in substitution of methionine (M) for valine (V) at residue 108 of the soluble form of the protein. While the two proteins have similar kinetic properties, 108M COMT loses activity more rapidly than 108V COMT at 37 degrees C. The cosubstrate S-adenosylmethionine (SAM) stabilizes the activity of 108M COMT at 40 degrees C. The 108M allele has been associated with increased risk for breast cancer, obsessive-compulsive disorder, and aggressive and highly antisocial manifestations of schizophrenia. In the current work, we have constructed homology models for both human COMT polymorphs and performed molecular dynamics simulations of these models at 25, 37, and 50 degrees C to explore the structural consequences of the 108V/M polymorphism. The simulations indicated that replacing valine with the larger methionine residue led to greater solvent exposure of residue 108 and heightened packing interactions between M108 and helices alpha2, alpha4 (especially with R78), and alpha5. These altered packing interactions propagated subtle changes between the polymorphic site and the active site 16 A away, leading to a loosening of the active site. At physiological temperature, 108M COMT sampled a larger distribution of conformations than 108V. 108M COMT was more prone to active-site distortion and had greater overall, and SAM binding site, solvent accessibility than 108V COMT at 37 degrees C. Similar structural perturbations were observed in the 108V protein only at 50 degrees C. Addition of SAM tightened up the cosubstrate pocket in both proteins and prevented the altered packing at the polymorphic site in 108M COMT.

Amino Acid Sequence↗
PubMed2006

Alpha-sheet: The toxic conformer in amyloid diseases?

A novel secondary structure, the alpha-sheet, was identified through molecular dynamics (MD) simulations of various proteins associated with amyloid diseases under amyloidogenic conditions. The structure was first predicted by Pauling and Corey, and it has been directly observed in crystal structures of "nonnatural peptides". There are occurrences of alpha-strands and alpha-sheets in the Protein Data Bank, but they are rare. We propose that alpha-sheet is formed during the conformational changes associated with amyloidosis and that it may represent the toxic conformer. Here, structural properties of the alpha-sheet, background information, and experimental support for this novel structure are presented. Finally we speculate about the possible role of this conformation in disease.

Amyloid↗
PubMed2005

Characterization of two distinct beta2-microglobulin unfolding intermediates that may lead to amyloid fibrils of different morphology.

The self-assembly of beta(2)-microglobulin into fibrils leads to dialysis-related amyloidosis. pH-mediated partial unfolding is required for the formation of the amyloidogenic intermediate that then self-assembles into amyloid fibrils. Two partially folded intermediates of beta(2)-microglobulin have been identified experimentally and linked to the formation of fibrils of distinct morphology, yet it remains difficult to characterize these partially unfolded states at high resolution using experimental approaches. Consequently, we have performed molecular dynamics simulations at neutral and low pH to determine the structures of these partially unfolded amyloidogenic intermediates. In the low-pH simulations, we observed the formation of alpha-sheet structure, which was first proposed by Pauling and Corey. Multiple simulations were performed, and two distinct intermediate state ensembles were identified that may account for the different fibril morphologies. The predominant early unfolding intermediate was nativelike in structure, in agreement with previous NMR studies. The late unfolding intermediate was significantly disordered, but it maintained an extended elongated structure, with hydrophobic clusters and residual alpha-extended chain strands in specific regions of the sequence that map to amyloidogenic peptides. We propose that the formation of alpha-sheet facilitates self-assembly into partially unfolded prefibrillar amyloidogenic intermediates.

Amyloid↗
PubMed2005

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

Several neurodegenerative diseases are linked to expanded repeats of glutamine residues, which lead to the formation of amyloid fibrils and neuronal death. The length of the repeats correlates with the onset of Huntington's disease, such that healthy individuals have <38 residues and individuals with >38 repeats exhibit symptoms. Because it is difficult to obtain atomic-resolution structural information for poly(l-glutamine) (polyQ) in aqueous solution experimentally, we performed molecular dynamics simulations to investigate the conformational behavior of this homopolymer. In simulations of 20-, 40-, and 80-mer polyQ, we observed the formation of the "alpha-extended chain" conformation, which is characterized by alternating residues in the alpha(L) and alpha(R) conformations to yield a sheet. The structural transition from disordered random-coil conformations to the alpha-extended chain conformation exhibits modest length and temperature dependence, in agreement with the experimental observation that aggregation depends on length and temperature. We propose that fibril formation in polyQ may occur through an alpha-sheet structure, which was proposed by Pauling and Corey. Also, we propose an atomic-resolution model of how the inhibitory peptide QBP1 (polyQ-binding peptide 1) may bind to polyQ in an alpha-extended chain conformation to inhibit fibril formation.

Amyloid↗
PubMed2005

Ensemble versus single-molecule protein unfolding.

Molecular dynamics (MD) simulations are the classic single-molecule "experiments," providing atomic-resolution structural and dynamic information. However, the single-molecule nature of the technique has also been its shortcoming, with frequent criticisms of sampling inadequacies and questions regarding the ensemble behavior of large numbers of molecules. Given the increase in computer power, we now address this issue by performing a large number of simulations and comparing individual and ensemble properties. One hundred independent MD simulations of the protein chymotrypsin inhibitor 2 were carried out for 20 ns each at 498 K in water to more fully describe the potentially diverse routes of protein unfolding and investigate how representative a single trajectory can be. Rapid unfolding was observed in all cases with the trajectories distributed about an average "ensemble" path in which secondary and tertiary structure was lost concomitantly, with tertiary structure loss occurring slightly faster. Individual trajectories did, however, sample conformations far from the average path with very heterogeneous time-dependent properties. Nevertheless, all of the simulations but one followed the average ensemble pathway, such that a small number of simulations (5-10) are sufficient to capture the average properties of these states and the unfolding pathway.

Hydrogen Bonding↗
PubMed2005

Simulation and experiment conspire to reveal cryptic intermediates and a slide from the nucleation-condensation to framework mechanism of folding.

There is a change from three-state to two-state kinetics of folding across the homeodomain superfamily of proteins as the mechanism slides from framework to nucleation-condensation. The tendency for framework folding in this family correlates with inherent helical propensity. The cellular myeloblastis protein (c-Myb) falls in the mechanistic transition region. An earlier, preliminary report of protein engineering experiments and molecular dynamics simulations (MD) showed that the folding mechanism for this protein has aspects of both the nucleation-condensation and framework models. In the more in-depth analysis of the MD trajectories presented here, we find that folding may be attributed to both of these mechanisms in different regions of the protein. The folding of the loop, middle helix, and turn is best described by nucleation-condensation, whereas folding of the N and C-terminal helices may be described by the framework model. Experimentally, c-Myb folds by apparent two-state kinetics, but the MD simulations predict that the kinetics hide a high-energy intermediate. We stabilized this hypothetical folding intermediate by deleting a residue (P174) in the loop between its second and third helices, and the mutant intermediate is long-lived in the simulations. Equilibrium and kinetic experiments demonstrate that folding of the DeltaP174 mutant is indeed three-state. The presence and shape of the intermediate observed in the simulations were confirmed by small angle X-ray scattering experiments.

Computer Simulation↗
PubMed2005

The structure of the major transition state for folding of an FF domain from experiment and simulation.

We have analysed the transition state of folding of the four-helix FF domain from HYPA/FBP11 by high-resolution experiment and simulation as part of a continuing effort to understand the principles of folding and the refinement of predictive methods. The major transition state for folding was subjected to a Phi-value analysis utilising 50 mutants. The transition state contained a nucleus for folding centred around the end of helix 1 (H1) and the beginning of helix 2 (H2). Secondary structure in this region was fully formed (PhiF=0.9-1) and tertiary interactions were well developed. Interactions in the distal part of the native structure were weak (PhiF=0-0.2). The hydrophobic core and other parts of the protein displayed intermediate Phi-values, suggesting that interactions coalesce as the end of H1 and beginning of H2 are in the process of being formed. The distribution of Phi-values resembled that of barnase, which folds via an intermediate, rather than that of CI2 which folds by a concerted nucleation-condensation mechanism. The overall picture of the transition state structure identified in molecular dynamics simulations is in qualitative agreement, with the turn connecting H1 and H2 being formed, a loosened core, and H4 partially unfolded and detached from the core. There are some differences in the details and interpretation of specific Phi-values.

Carrier Proteins↗
PubMed2005

Local environmental effects on the structure of the prion protein.

Prion diseases are neurodegenerative diseases causally linked to the partial unfolding and subsequent misfolding and aggregation of the prion protein (PrP). While most proteins fold into a single low energy state, PrP can fold into two distinct isoforms. In its innocuous state, denoted as PrPC, the protein has predominantly alpha-helical secondary structure, however, PrPC can misfold into an isoform rich in extended structure capable of forming toxic and infectious aggregates. While prion disease is believed to be a protein-only disease, one not requiring any non-protein elements for propagation, the different environments the protein finds itself in vivo likely influence its ability to misfold and aggregate. In this review we will examine various molecules, covalent modifications and environments PrP faces in vivo and the effect they have on PrP's local environment and, potentially, conformation. Included in this discussion are: (1) pH, (2) carbohydrates, (3) lipid membranes, (4) metal ions, and (5) small molecules.

Animals↗
PubMed2005

Simulation and experiment at high temperatures: ultrafast folding of a thermophilic protein by nucleation-condensation.

We used Phi-value analysis to characterise the transition state for folding of a thermophilic protein at the relatively high temperature of 325 K. PhiF values for the folding of the three-helix bundle, peripheral subunit binding domain from Bacillus stearothermophilus (E3BD) were determined by temperature-jump experiments in the absence of chemical denaturants. E3BD folded in microseconds through a highly diffuse transition state. Excellent agreement was observed between experiment and the results from eight (independent) molecular dynamics simulations of unfolding at 373 K. We used a combination of heteronuclear NMR experiments and molecular dynamics simulations to characterise the denatured ensemble, and found that it contained very little persistent, residual structure. However, those regions that adopt helical structure in the native state were found by simulation to be poised for helix formation in the denatured state. These regions also had significant structure in the transition state for folding. The overall folding pathway appears to be nucleation-condensation.

Bacterial Proteins↗
PubMed2005

Insight into ribonuclease A domain swapping by molecular dynamics unfolding simulations.

Bovine pancreatic ribonuclease (RNase A) deserves a special place among the numerous proteins that form oligomers by three-dimensional domain swapping. In fact, under destabilizing conditions and at high protein concentrations, it can swap two different domains, the N-terminal alpha-helix or the C-terminal beta-strand, leading to dimers with different quaternary structures. With the change in the unfolding conditions, the relative abundance of the two dimers varies, and the prevalence of one dimer over the other is inverted. To investigate the dynamic behavior of the termini, four independent 10 ns high-temperature molecular dynamics simulations of RNase A were carried out at two different pH values in an attempt to reproduce the experimental conditions of neutral and very low pH that favor the formation of the N- and C-terminal domain-swapped dimers, respectively. In agreement with experimental data, under mild unfolding conditions, a partial or complete opening of the N-terminal arm is observed, whereas the dislocation of the C-terminus away from the core of the structure occurs only during the low-pH simulations. Furthermore, the picture emerging from this study indicates that the same protein can have different pathways for domain swapping. Indeed, in RNase A the C-terminal swapping requires a substantial unfolding of the monomers, whereas the N-terminal swapping can occur through only partial unfolding.

Animals↗
PubMed2005

Cutoff size need not strongly influence molecular dynamics results for solvated polypeptides.

The correct treatment of van der Waals and electrostatic nonbonded interactions in molecular force fields is essential for performing realistic molecular dynamics (MD) simulations of solvated polypeptides. The most computationally tractable treatment of nonbonded interactions in MD utilizes a spherical distance cutoff (typically, 8-12 A) to reduce the number of pairwise interactions. In this work, we assess three spherical atom-based cutoff approaches for use with all-atom explicit solvent MD: abrupt truncation, a CHARMM-style electrostatic shift truncation, and our own force-shifted truncation. The chosen system for this study is an end-capped 17-residue alanine-based alpha-helical peptide, selected because of its use in previous computational and experimental studies. We compare the time-averaged helical content calculated from these MD trajectories with experiment. We also examine the effect of varying the cutoff treatment and distance on energy conservation. We find that the abrupt truncation approach is pathological in its inability to conserve energy. The CHARMM-style shift truncation performs quite well but suffers from energetic instability. On the other hand, the force-shifted spherical cutoff method conserves energy, correctly predicts the experimental helical content, and shows convergence in simulation statistics as the cutoff is increased. This work demonstrates that by using proper and rigorous techniques, it is possible to correctly model polypeptide dynamics in solution with a spherical cutoff. The inherent computational advantage of spherical cutoffs over Ewald summation (and related) techniques is essential in accessing longer MD time scales.

Circular Dichroism↗