PubMed Health⌕ Search

Biomedical subjects

V Daggett

Publications and source records attributed to V Daggett.

At least 19 recordsLinked to original sources
PubMed1999

Engineering out motion: introduction of a de novo disulfide bond and a salt bridge designed to close a dynamic cleft on the surface of cytochrome b5.

A previous molecular dynamics (MD) simulation of cytochrome b5 (cyt b5) at 25 degrees C displayed localized dynamics on the surface of the protein giving rise to the periodic formation of a cleft that provides access to the heme through a protected hydrophobic channel [Storch and Daggett (1995) Biochemistry 34, 9682]. Here we describe the production and testing of mutants designed to prevent the cleft from opening using a combination of experimental and theoretical techniques. Two mutants have been designed to close the surface cleft: S18D to introduce a salt bridge and S18C:R47C to incorporate a disulfide bond. The putative cleft forms between two separate cores of the protein: one is structural in nature and can be monitored through the fluorescence of Trp 22, and the other binds the heme prosthetic group and can be tracked via heme absorbance. An increase in motion localized to the cleft region was observed for each protein, except for the disulfide-containing variant, in MD simulations at 50 degrees C compared to simulations at 25 degrees C. For the disulfide-containing variant, the cleft remained closed. Both urea and temperature denaturation curves were nearly identical for wild-type and mutant proteins when heme absorbance was monitored. In contrast, fluorescence studies revealed oxidized S18C:R47C to be considerably more stable based on the midpoints of the denaturation transitions, Tm and U1/2. Moreover, the fluorescence changes for each protein were complete at approximately 50 degrees C and a urea concentration of approximately 3.9 M, significantly below the temperature and urea concentration (62 degrees C, 5 M urea) required to observe heme release. In addition, solvent accessibility based on acrylamide quenching of Trp 22 was lower in the S18C:R47C mutant, particularly at 50 degrees C, before heme release [presented in the accompanying paper (58)]. The results suggest that a constraining disulfide bond can be designed to inhibit dynamic cleft formation on the surface of cyt b5. Located near the heme, the native dynamics of the cleft may be functionally important for protein-protein recognition and/or complex stabilization.

Animals↗
PubMed1999

Engineering out motion: a surface disulfide bond alters the mobility of tryptophan 22 in cytochrome b5 as probed by time-resolved fluorescence and 1H NMR experiments.

In the accompanying paper [Storch et al. (1999) Biochemistry 38, 5054-5064] equilibrium denaturation studies and molecular dynamics (MD) simulations were used to investigate localized dynamics on the surface of cytochrome b5 (cyt b5) that result in the formation of a cleft. In those studies, an S18C:R47C disulfide mutant was engineered to inhibit cleft mobility. Temperature- and urea-induced denaturation studies revealed significant differences in Trp 22 fluorescence between the wild-type and mutant proteins. On the basis of the results, it was proposed that wild type populates a conformational ensemble that is unavailable to the disulfide mutant and is mediated by cleft mobility. As a result, the solvent accessibility of Trp 22 is decreased in S18C:R47C, suggesting that the local environment of this residue is less mobile due to the constraining effects of the disulfide on cleft dynamics. To further probe the structural effects on the local environment of Trp 22 caused by inhibition of cleft formation, we report here the results of steady-state and time-resolved fluorescence quenching, differential phase/modulation fluorescence anisotropy, and 1H NMR studies. In Trp fluorescence experiments, the Stern-Volmer quenching constant increases in wild type versus the oxidized disulfide mutant with increasing temperature. At 50 degrees C, KSV is nearly 1.5-fold greater in wild type compared to the oxidized disulfide mutant. In the reduced disulfide mutant, KSV was the same as wild type. The bimolecular collisional quenching constant, kq, for acrylamide quenching of Trp 22 increases 2.7-fold for wild type and only 1.8-fold for S18C:R47C, upon increasing the temperature from 25 to 50 degrees C. The time-resolved anisotropy decay at 25 degrees C was fit to a double-exponential decay for both the wild type and S18C:R47C. Both proteins exhibited a minor contribution from a low-amplitude fast decay, consistent with local motion of Trp 22. This component was more prevalent in the wild type, and the fractional contribution increased significantly upon raising the temperature. The fast rotational component of the S18C:R47C mutant was less sensitive to increasing temperature. A comparison of the 1H NMR monitored temperature titration of the delta-methyl protons of Ile 76 for wild type and oxidized disulfide mutant, S18C:R47C, showed a significantly smaller downfield shift for the mutant protein, suggesting that Trp 22 in the mutant protein experiences comparatively decreased cleft dynamics in core 2 at higher temperatures. Furthermore, comparison of the delta'-methyl protons of Leu 25 in the two proteins revealed a difference in the ratio of the equilibrium heme conformers of 1.2:1 for S18C:R47C versus 1.5:1 for wild type at 40 degrees C. The difference in equilibrium heme orientations between wild type and S18C:R47C suggests that the disulfide bond affects heme binding within core 1, possibly through damped cleft fluctuations. Taken together, the NMR and fluorescence studies support the proposal that an engineered disulfide bond inhibits the formation of a dynamic cleft on the surface of cyt b5.

Amino Acid Substitution↗
PubMed1998

Non-native interactions in protein folding intermediates: molecular dynamics simulations of hen lysozyme.

Molecular dynamics simulations of protein denaturation can complement and extend experimental studies of protein folding by providing atomic-level structural information about conformational transitions and any conformational states along the unfolding pathway. Previous unfolding simulations of hen egg-white lysozyme have resulted in intermediate structures with an unfolded alpha-domain and a structured beta-domain, which is inconsistent with experiment. In contrast, the beta-domain unfolded first in the two simulations presented here leaving a structured alpha-domain. Following this, intermediate states were identified that differ with respect to the packing of the helices and the elements of non-native structure adopted. The non-native structure is critical for explaining many of the experimental observations. Overall, the pooled ensemble of these intermediates is in agreement with the experimental data for the major kinetic intermediate, suggesting that the kinetic intermediate may be made up of distinct, but rapidly interconverting, partially folded conformations distinguished primarily by differences in helix packing.

Amino Acid Sequence↗
PubMed1998

Barstar has a highly dynamic hydrophobic core: evidence from molecular dynamics simulations and nuclear magnetic resonance relaxation data.

The dynamic behavior of the ribonuclease inhibitor barstar has been investigated by molecular dynamics (MD) simulations in explicit water. Two 2.5 ns MD simulations were performed, and an ensemble of 25 000 structures was generated. This ensemble reproduces the solution structures and is consistent with the experimental structural restraints from NMR spectroscopy. Reorientation of the backbone NH bond vectors and side chain methyl groups was monitored by calculation of autocorrelation functions and the generalized S2 order parameters. Order parameters derived for motion in the approximately 100 ps time scale were compared with those obtained from NMR relaxation measurements. Consistent with experiment, the backbone NH bond vectors were relatively rigid. In contrast, the side chain methyl groups exhibited a wide dynamic range, from restricted motion comparable to that of the backbone to rapid unrestricted motion. The order parameters for the methyl groups correlate well with their spatial separation from the backbone and are residue-type dependent. Smaller S2axis values were observed for leucine methyl groups, in part due to side chain hopping between two predominant rotamers (g+t and tg-). Motions such as the flipping of aromatic rings and the hopping of leucine side chains were prevalent within the hydrophobic core, suggesting that the core is fluid-like with low energy barriers between native conformational substates. Thus, our studies suggest that the entropy of the native state can be significant and should not be discounted in thermodynamic considerations of protein folding. On the basis of our results, the side chain motion represents the primary source of the residual entropy of the native state and entropic considerations based solely on backbone dynamics would be incomplete.

Bacillus↗
PubMed1998

Synergy between simulation and experiment in describing the energy landscape of protein folding.

Experimental data from protein engineering studies and NMR spectroscopy have been used by theoreticians to develop algorithms for helix propensity and to benchmark computer simulations of folding pathways and energy landscapes. Molecular dynamic simulations of the unfolding of chymotrypsin inhibitor 2 (CI2) have provided detailed structural models of the transition state ensemble for unfolding/folding of the protein. We now have used the simulated transition state structures to design faster folding mutants of CI2. The models pinpoint a number of unfavorable local interactions at the carboxyl terminus of the single alpha-helix and in the protease-binding loop region of CI2. By removing these interactions or replacing them with stabilizing ones, we have increased the rate of folding of the protein up to 40-fold (tau = 0.4 ms). This correspondence, and other examples of agreement between experiment and theory in general, Phi-values and molecular dynamics simulations, in particular, suggest that significant progress has been made toward describing complete folding pathways at atomic resolution by combining experiment and simulation.

Kinetics↗
PubMed1998

Simulations of the structural and dynamical properties of denatured proteins: the "molten coil" state of bovine pancreatic trypsin inhibitor.

The dynamic nature of denatured, unfolded proteins makes it difficult to characterize their structures experimentally. To complement experiment and to obtain more detailed information about the structure and dynamic behavior of the denatured state, we have performed eleven 2.5 ns molecular dynamics simulations of reduced bovine pancreatic trypsin inhibitor (BPTI) at high temperature in water and a control simulation at 298 K, for a total of 30 ns of simulation time. In a neutral pH environment (acidic residues ionized), the unfolded protein structures were compact with an average radius of gyration 9% greater than the native state. The compact conformations resulted from the transient formation of non-native hydrophobic clusters, turns and salt bridges. However, when the acidic residues were protonated, the protein periodically expanded to a radius of gyration of 18 to 20 A. The early steps in unfolding were similar in the different simulations until passing through the major transition state of unfolding. Afterwards, unfolding proceeded through one of two general pathways with respect to secondary structure: loss of the C-terminal helix followed by loss of beta-structure or the opposite. To determine whether the protein preferentially sampled particular conformational substates in the denatured state, pairwise Calpha root-mean-square deviations were measured between all structures, but similar structures were found between only two trajectories. Yet, similar composite properties (secondary structure content, side-chain and water contacts, solvent accessible surface area, etc.) were observed for the structures that unfolded through different pathways. Somewhat surprisingly, the unfolded structures are in agreement with both past experiments suggesting that reduced BPTI is a random coil and more recent experiments providing evidence for non-random structure, demonstrating how ensembles of fluctuating structures can give rise to experimental observables that are seemingly at odds.

Aprotinin↗
PubMed1998

Molecular dynamics simulation of the unfolding of barnase: characterization of the major intermediate.

The folding/unfolding pathway of barnase has been studied extensively using the protein engineering method, which has provided indirect structural information for the transition state and the major folding intermediate. To further characterize the structural properties of the intermediate, we have simulated the thermal denaturation of barnase beginning from the average NMR structure. Our results indicate that there are at least two intermediates on the unfolding pathway. The three hydrophobic cores are partially formed in the major intermediate (I1), with core1 and core3 being slightly stronger than core2. Helix alpha 1 is substantially formed, with the center being stronger than the termini. The first turn of alpha 2 is lost and alpha 3 is unfolded. The center of the beta-sheet is substantially formed, but the edges are disrupted. These structural characteristics are in good qualitative agreement with the experimental data. For semi-quantitative comparison with experimental data, the extent of native structure of individual residues is characterized by a structure index, S, that reflects both secondary and tertiary structure. There is good agreement between S and the experimentally measured phi values, which are based on energetics, except for three residues. These residues are polar and non-conservative mutations were made to obtain phi values, which can complicate structural interpretations. These residues make strong side-chain interactions in I1, but the backbone structure is disrupted, leading to low S values. Thus, this discordance highlights possible limitations in both the phi value and S value analyses: strong polar interactions in the intermediate may give rise to high phi values that are not reflective of structure per se; however, due to sampling limitations, any one simulation is not expected to capture all of the features of the true conformational ensemble. In any case, these simulations provide an experimentally testable, atomic-level structural model for the major folding intermediate of barnase, as well as the detailed pathway from the native to the intermediate state.

Bacterial Proteins↗
PubMed1998

Molecular dynamics simulations of hydrophobic collapse of ubiquitin.

Nine nonnative conformations of ubiquitin, generated during two different thermal denaturation trajectories, were simulated under nearly native conditions (62 degrees C). The simulations included all protein and solvent atoms explicitly, and simulation times ranged from 1-2.4 ns. The starting structures had alpha-carbon root-mean-square deviations (RMSDs) from the crystal structure of 4-12 A and radii of gyration as high as 1.3 times that of the native state. In all but one case, the protein collapsed when the temperature was lowered and sampled conformations as compact as those reached in a control simulation beginning from the crystal structure. In contrast, the protein did not collapse when simulated in a 60% methanol:water mixture. The behavior of the protein depended on the starting structure: during simulation of the most native-like starting structures (<5 A RMSD to the crystal structure) the RMSD decreased, the number of native hydrogen bonds increased, and the secondary and tertiary structure increased. Intermediate starting structures (5-10 A RMSD) collapsed to the radius of gyration of the control simulation, hydrophobic residues were preferentially buried, and the protein acquired some native contacts. However, the protein did not refold. The least native starting structures (10-12 A RMSD) did not collapse as completely as the more native-like structures; instead, they experienced large fluctuations in radius of gyration and went through cycles of expansion and collapse, with improved burial of hydrophobic residues in successive collapsed states.

Computer Simulation↗
PubMed1998

Structure-function aspects of prion proteins.

Prions diseases are fatal neurodegenerative disorders resulting from conformational changes in the prion protein from the normal cellular form, PrPC, to the infectious scrapie isoform, PrPSc. High resolution structures for PrPC are now available, and biochemical investigations are shedding light on the nature and determinants of the conformational transition. Together, these studies are beginning to provide a framework to describe structure-function relationships of the prion protein.

Animals↗
PubMed1997

Characterization of residual structure in the thermally denatured state of barnase by simulation and experiment: description of the folding pathway.

Residual structure in the denatured state of a protein may contain clues about the early events in folding. We have simulated by molecular dynamics the denatured state of barnase, which has been studied by NMR spectroscopy. An ensemble of 10(4) structures was generated after 2 ns of unfolding and following for a further 2 ns. The ensemble was heterogeneous, but there was nonrandom, residual structure with persistent interactions. Helical structure in the C-terminal portion of helix alpha1 (residues 13-17) and in helix alpha2 as well as a turn and nonnative hydrophobic clustering between beta3 and beta4 were observed, consistent with NMR data. In addition, there were tertiary contacts between residues in alpha1 and the C-terminal portion of the beta-sheet. The simulated structures allow the rudimentary NMR data to be fleshed out. The consistency between simulation and experiment inspires confidence in the methods. A description of the folding pathway of barnase from the denatured to the native state can be constructed by combining the simulation with experimental data from phi value analysis and NMR.

Amino Acid Sequence↗
PubMed1997

Selective neuronal targeting in prion disease.

The pattern of scrapie prion protein (PrP(Sc)) accumulation in the brain is different for each prion strain. We tested whether the PrP(Sc) deposition pattern is influenced by the Asn-linked oligosaccharides of PrP(C) in transgenic mice. Deletion of the first oligosaccharide altered PrP(C) trafficking and prevented infection with two prion strains. Deletion of the second did not alter PrP(C) trafficking, permitted infection with one prion strain, and had a profound effect on the PrP(Sc) deposition pattern. Our data raise the possibility that glycosylation can modify the conformation of PrP(C). Glycosylation could affect the affinity of PrP(C) for a particular conformer of PrP(Sc), thereby determining the rate of nascent PrP(Sc) formation and the specific patterns of PrP(Sc) deposition.

Animals↗
PubMed1997

The prion folding problem.

Prion diseases are neurodegenerative disorders in which dramatic conformational change in the structure of the prion protein is the fundamental event. This structural transition involves the loss of substantial alpha-helical content and the acquisition of beta-sheet structure. A convergence of recent biological and structural studies argues that the mechanism underlying the prion diseases is truly unprecedented.

Biopolymers↗
PubMed1996

Structural consequences of heme removal: molecular dynamics simulations of rat and bovine apocytochrome b5.

Molecular dynamics simulations of rat and bovine apocytochrome b5 were performed to investigate the structural and dynamical consequences of heme removal. A crystal structure is available for the bovine holoprotein, while experimental studies of apocytochrome b5 have focused on the rat protein. The rat and bovine proteins are 93% homologous by sequence, and the sequence differences (six residues) appear to have no effect on the structure of the native holoprotein, as seen by the correlation of a bovine simulation with rat holocytochrome b5 experimental data (Storch & Daggett, 1995). There was a marked effect, however, on the structure and dynamics of the apo form. The bovine apocytochrome b5 simulation displayed subtle inconsistencies when compared to the experimental results on the rat apoprotein. Therefore, the rat protein was constructed from the bovine crystal structure coordinates. The MD simulation of the rat apoprotein displayed greater deviations from the crystal structure, yet it was in much closer agreement to the experimental data for the apoprotein. Additionally, the six variant residues fall in the regions where the bovine protein deviated from experiment. The two hydrophobic cores of the rat protein behaved very differently. Core 2 was well maintained, retained native-like structure, and is in good agreement with NMR data (Moore & Lecomte, 1990). Conversely, core 1, which is normally constrained by the prosthetic heme group, exhibited conformational heterogeneity, increased mobility, and some loss of secondary structure. Thus, the model of rat apocytochrome b5 complements past studies by providing structural information about core 1 that has proved difficult to obtain by experiment. The bovine simulation serves as a prediction, since little to no experimental data exist for this form of the apoprotein.

Animals↗
PubMed1996

Identification and characterization of the unfolding transition state of chymotrypsin inhibitor 2 by molecular dynamics simulations.

Temperature-induced unfolding of chymotrypsin inhibitor 2 (CI2) in water has been investigated using molecular dynamics simulations. One simulation (2.2 ns) has been analyzed in detail and three additional simulations (each > or = 1 ns) were performed to check the generality of the results. Concurrent loss of secondary and tertiary structure during unfolding was observed in all the simulations. For each simulation, the major transition state of unfolding was identified based on conformational analysis of protein structures along the unfolding trajectory. The transition state has a considerably weakened hydrophobic core and disrupted secondary structure. Nevertheless, the overall structure of the transition state is closer to the native state than to the unfolded state. The disruption of the hydrophobic core appears to be rate limiting. However, other energy barriers have to be overcome before reaching the major transition state. A method is described to quantitatively compare the structure of the simulated transition state with that characterized by protein engineering experiments. Good agreement with the experimental data is obtained for all four transition state models (the correlation coefficient R = 0.80 to 0.93) and the average over all four models gives the best correlation (R = 0.94). These simulations provide the first comprehensive atomic-level view of what the unfolding transition state of C12 may look like.

Chymotrypsin↗
PubMed1996

Structure of the transition state for folding of a protein derived from experiment and simulation.

Independent experimental and theoretical studies of the unfolding of barley chymotrypsin inhibitor 2 (CI2) are compared in an attempt to derive plausible three-dimensional structural models of the transition states. A very simple structure index is calculated along the sequence for the molecular dynamics-generated transition state models to facilitate comparison with the phi F values. The two are in good agreement overall (correlation coefficient = 0.87), which suggests that the theoretical models should provide a structural framework for interpretation of the phi F values. Both experiment and simulation indicate that the transition state is a distorted form of the native state in which the alpha-helix is weakened but partially intact and the beta-sheet is quite disrupted. As inferred from the phi f values and observed directly in the simulations, the unfolding of CI2 is cooperative and there is a "folding core" comprising a patch on the alpha-helix and a portion of the beta-sheet, nucleated by interactions between Ala16, Ile49 and other neighbouring residues. The protein becomes less structured radiating away from this core. Overall the data indicate that CI2 folds by a nucleation-collapse mechanism. In the absence of experimental information, we have little confidence that the molecular dynamics simulations are correct, especially when only one or a few simulations are performed. On the other hand, even though the experimentally derived phi values may reflect the extent of overall structure formation, they do not provide an actual atomic-resolution three dimensional structure of the transition state. By combining the two approaches, however, we have a framework for interpreting phi F values and can hopefully arrive at a more trustworthy model of the transition state. The process is in some ways similar to the combination of molecular dynamics and NMR data to solve the tertiary structure of proteins.

Amino Acid Sequence↗
PubMed1996

Molecular dynamics simulations of apocytochrome b562--the highly ordered limit of molten globules.

BACKGROUND: Cytochrome B562 is a heme-containing, four-helix bundle. It has been proposed that the apo form of the protein is a molten globule. We present a molecular dynamics study of apocytochrome b562 to investigate its structural and dynamic properties. RESULTS: Our simulations suggest that all four helices are essentially intact and confirm that the experimental difficulties of assigning helical NOEs in the C-terminal helix are not due to structural disorder. The increased 'moltenness' of the apoprotein is due to an increased mobility of the sidechains. The small observed increase in compressibility for the apoprotein is proposed to be the result of an increase in the intrinsic protein compressibility, which is opposed by the increase in the size of the protein hydration shell. CONCLUSIONS: Apocytochrome b562 is postulated to be near the highly ordered limit of the molten globule state, a structure whose molten character is due primarily to increased sidechain mobility with concurrent loss in tertiary contacts between the helices, rather than changes in the folding topology or substantially increased disorder of the secondary structure.

Animals↗
PubMed1995

Molecular dynamics simulation of cytochrome b5: implications for protein-protein recognition.

Cytochrome b5 participates in electron-transfer reactions with a variety of different proteins. To explore how this protein might discern between structurally varied proteins, we have performed a molecular dynamics simulation focusing on its structural stability and dynamic behavior in solution. The protein was simulated in water at 298 K and pH 6.9 for 2.5 ns. The protein deviated significantly from the crystal structure midway through the simulation, but ultimately the crystalline conformation was regained. The simulation was at all times well behaved as judged by comparison to structural NMR data obtained in solution. One region of the protein backbone that deviated from the crystal conformation contains acidic residues implicated in electrostatic-based protein-protein recognition. The mobility in this region caused the protein to display different patterns of residues at the surface with time, as well as the formation of a large cleft partially exposing the hydrophobic core lining the heme pocket. Furthermore, the position and cyclical formation of this cleft suggest that hydrophobic interactions may be important in protein-protein recognition events and possibly even electron transfer, as the cleft allows for easy access to the heme group. These results indicate that thermal motion could provide a low-energy mechanism for controlling recognition events. Thus, the dynamical behavior observed through the varying solution conformations sampled may be important in influencing the diverse range of protein-protein interactions in which cytochrome b5 participates.

Computer Simulation↗