Biomedical subjects

Christopher D Snow

Publications and source records attributed to Christopher D Snow.

13 recordsLinked to original sources

PubMed2006

Electric fields at the active site of an enzyme: direct comparison of experiment with theory.

The electric fields produced in folded proteins influence nearly every aspect of protein function. We present a vibrational spectroscopy technique that measures changes in electric field at a specific site of a protein as shifts in frequency (Stark shifts) of a calibrated nitrile vibration. A nitrile-containing inhibitor is used to deliver a unique probe vibration to the active site of human aldose reductase, and the response of the nitrile stretch frequency is measured for a series of mutations in the enzyme active site. These shifts yield quantitative information on electric fields that can be directly compared with electrostatics calculations. We show that extensive molecular dynamics simulations and ensemble averaging are required to reproduce the observed changes in field.

Aldehyde Reductase↗

PubMed2006

Kinetic definition of protein folding transition state ensembles and reaction coordinates.

Using distributed molecular dynamics simulations we located four distinct folding transitions for a 39-residue betabetaalphabeta protein fold. To characterize the nature of each room temperature transition, we calculated the probability of transmission for 500 points along each free energy barrier. We introduced a method for determining transition states by employing the transmission probability, Ptrans, and determined which conformations were transition state ensemble members (Ptrans approximately 0.5). The transmission probability may be used to characterize the barrier in several ways. For example, we ran simulations at 82 degrees C, determined the change in Ptrans with temperature for all 2,000 conformations, and quantified Hammond behavior directly using Ptrans correlation. Additionally, we propose that diffusion along Ptrans may provide the configurational diffusion rate at the top of the barrier. Specifically, given a transition state conformation x0 with estimated Ptrans=0.5, we selected a large set of subsequent conformations from independent trajectories, each exactly a small time deltat after x0 (250 ps). Calculating Ptrans for the new trial conformations, we generated the P(Ptrans|deltat=250 ps) distribution that reflected diffusion. This approach provides a novel perspective on the diffusive nature of a protein folding transition and provides a framework for a quantitative study of activated relaxation kinetics.

Computer Simulation↗

PubMed2005

Direct calculation of the binding free energies of FKBP ligands.

Direct calculations of the absolute free energies of binding for eight ligands to FKBP protein were performed using the Fujitsu BioServer massively parallel computer. Using the latest version of the general assisted model building with energy refinement (AMBER) force field for ligand model parameters and the Bennett acceptance ratio for computing free-energy differences, we obtained an excellent linear fit between the calculated and experimental binding free energies. The rms error from a linear fit is 0.4 kcal/mol for eight ligand complexes. In comparison with a previous study of the binding energies of these same eight ligand complexes, these results suggest that the use of improved model parameters can lead to more predictive binding estimates, and that these estimates can be obtained with significantly less computer time than previously thought. These findings make such direct methods more attractive for use in rational drug design.

Biophysics↗

PubMed2005

Dimerization of the p53 oligomerization domain: identification of a folding nucleus by molecular dynamics simulations.

Dimerization of the p53 oligomerization domain involves coupled folding and binding of monomers. To examine the dimerization, we have performed molecular dynamics (MD) simulations of dimer folding from the rate-limiting transition state ensemble (TSE). Among 799 putative transition state structures that were selected from a large ensemble of high-temperature unfolding trajectories, 129 were identified as members of the TSE via calculation of a 50% transmission coefficient from at least 20 room-temperature simulations. This study is the first to examine the refolding of a protein dimer using MD simulations in explicit water, revealing a folding nucleus for dimerization. Our atomistic simulations are consistent with experiment and offer insight that was previously unobtainable.

Computer Simulation↗

PubMed2005

How well can simulation predict protein folding kinetics and thermodynamics?

Simulation of protein folding has come a long way in five years. Notably, new quantitative comparisons with experiments for small, rapidly folding proteins have become possible. As the only way to validate simulation methodology, this achievement marks a significant advance. Here, we detail these recent achievements and ask whether simulations have indeed rendered quantitative predictions in several areas, including protein folding kinetics, thermodynamics, and physics-based methods for structure prediction. We conclude by looking to the future of such comparisons between simulations and experiments.

Biophysics↗

PubMed2004

Using path sampling to build better Markovian state models: predicting the folding rate and mechanism of a tryptophan zipper beta hairpin.

We propose an efficient method for the prediction of protein folding rate constants and mechanisms. We use molecular dynamics simulation data to build Markovian state models (MSMs), discrete representations of the pathways sampled. Using these MSMs, we can quickly calculate the folding probability (P(fold)) and mean first passage time of all the sampled points. In addition, we provide techniques for evaluating these values under perturbed conditions without expensive recomputations. To demonstrate this method on a challenging system, we apply these techniques to a two-dimensional model energy landscape and the folding of a tryptophan zipper beta hairpin.

Computational Biology↗

PubMed2004

Trp zipper folding kinetics by molecular dynamics and temperature-jump spectroscopy.

We studied the microsecond folding dynamics of three beta hairpins (Trp zippers 1-3, TZ1-TZ3) by using temperature-jump fluorescence and atomistic molecular dynamics in implicit solvent. In addition, we studied TZ2 by using time-resolved IR spectroscopy. By using distributed computing, we obtained an aggregate simulation time of 22 ms. The simulations included 150, 212, and 48 folding events at room temperature for TZ1, TZ2, and TZ3, respectively. The all-atom optimized potentials for liquid simulations (OPLS(aa)) potential set predicted TZ1 and TZ2 properties well; the estimated folding rates agreed with the experimentally determined folding rates and native conformations were the global potential-energy minimum. The simulations also predicted reasonable unfolding activation enthalpies. This work, directly comparing large simulated folding ensembles with multiple spectroscopic probes, revealed both the surprising predictive ability of current models as well as their shortcomings. Specifically, for TZ1-TZ3, OPLS for united atom models had a nonnative free-energy minimum, and the folding rate for OPLS(aa) TZ3 was sensitive to the initial conformation. Finally, we characterized the transition state; all TZs fold by means of similar, native-like transition-state conformations.

Computer Simulation↗

PubMed2003

Surface salt bridges, double-mutant cycles, and protein stability: an experimental and computational analysis of the interaction of the Asp 23 side chain with the N-terminus of the N-terminal domain of the ribosomal protein l9.

Experimental and theoretical double-mutant cycles have been used to investigate a salt bridge in the N-terminal domain of the protein L9. Aspartic acid 23 is the only acidic residue involved in a well-defined pairwise interaction, namely, a partially solvent-exposed salt bridge with the protonated N-terminus of the protein. Mutations were studied in which Asp 23 was substituted by alanine, asparagine, and nitrile alanine. Interactions with the N-terminus were probed by comparisons between proteins with a protonated and acetylated N-terminus. The mutants were all folded, and the structures were unchanged from wild type as judged by CD and 2-D NMR. The coupling free energy between the N-terminus and the side chain of Asp 23 measured through double-mutant cycle analysis was favorable and ranged from -0.7 to -1.7 kcal mol(-)(1), depending upon the set of mutants used. This relatively large coupling free energy for a surface salt bridge likely arises from geometric factors that reduce the entropy loss associated with salt-bridge formation and from structural relaxation in the mutants. Coupling free energies computed with continuum electrostatic calculations agreed well with the experimental values when full account was taken of all potential interactions, particularly those involving Asp 23 and the acetylated N-terminus as well as interactions with solvent. The measured and calculated coupling free energy decreased only slightly when the salt concentration was increased from 100 to 750 mM NaCl. The calculations suggest that the coupling free energy between D23 and the N-terminus measured through the experimental double-mutant cycle analysis is significantly smaller than the actual interaction free energy between the groups in the wild-type structure because of the inapplicability of assumptions frequently used to interpret double-mutant cycles.

Amino Acid Substitution↗

PubMed2003

Atomistic protein folding simulations on the submillisecond time scale using worldwide distributed computing.

Atomistic simulations of protein folding have the potential to be a great complement to experimental studies, but have been severely limited by the time scales accessible with current computer hardware and algorithms. By employing a worldwide distributed computing network of tens of thousands of PCs and algorithms designed to efficiently utilize this new many-processor, highly heterogeneous, loosely coupled distributed computing paradigm, we have been able to simulate hundreds of microseconds of atomistic molecular dynamics. This has allowed us to directly simulate the folding mechanism and to accurately predict the folding rate of several fast-folding proteins and polymers, including a nonbiological helix, polypeptide alpha-helices, a beta-hairpin, and a three-helix bundle protein from the villin headpiece. Our results demonstrate that one can reach the time scales needed to simulate fast folding using distributed computing, and that potential sets used to describe interatomic interactions are sufficiently accurate to reach the folded state with experimentally validated rates, at least for small proteins.

Algorithms↗

PubMed2002

The Trp cage: folding kinetics and unfolded state topology via molecular dynamics simulations.

Using over 75 mus of molecular dynamics simulation, we have generated several thousand folding simulations of the 20-residue Trp cage at experimental temperature and solvent viscosity. A total of 116 independent folding simulations reach RMSDcalpha values below 3 A RMSDcalpha, some as close as 1.4 A RMSDcalpha. We estimate a folding time of 5.5+/-3.5 mus, a rate that is in reasonable agreement with experimental kinetics. Finally, we characterize both the folded and unfolded ensemble under native conditions and note that the average topology of the unfolded ensemble is very similar to the topology of the native state.

Amino Acid Sequence↗

PubMed2002

Simulation of folding of a small alpha-helical protein in atomistic detail using worldwide-distributed computing.

By employing thousands of PCs and new worldwide-distributed computing techniques, we have simulated in atomistic detail the folding of a fast-folding 36-residue alpha-helical protein from the villin headpiece. The total simulated time exceeds 300 micros, orders of magnitude more than previous simulations of a molecule of this size. Starting from an extended state, we obtained an ensemble of folded structures, which is on average 1.7A and 1.9A away from the native state in C(alpha) distance-based root-mean-square deviation (dRMS) and C(beta) dRMS sense, respectively. The folding mechanism of villin is most consistent with the hydrophobic collapse view of folding: the molecule collapses non-specifically very quickly ( approximately 20ns), which greatly reduces the size of the conformational space that needs to be explored in search of the native state. The conformational search in the collapsed state appears to be rate-limited by the formation of the aromatic core: in a significant fraction of our simulations, the C-terminal phenylalanine residue packs improperly with the rest of the hydrophobic core. We suggest that the breaking of this interaction may be the rate-determining step in the course of folding. On the basis of our simulations we estimate the folding rate of villin to be approximately 5micros. By analyzing the average features of the folded ensemble obtained by simulation, we see that the mean folded structure is more similar to the native fold than any individual folded structure. This finding highlights the need for simulating ensembles of molecules and averaging the results in an experiment-like fashion if meaningful comparison between simulation and experiment is to be attempted. Moreover, our results demonstrate that (1) the computational methodology exists to simulate the multi-microsecond regime using distributed computing and (2) that potential sets used to describe interatomic interactions may be sufficiently accurate to reach the folded state, at least for small proteins. We conclude with a comparison between our results and current protein-folding theory.

Amino Acid Sequence↗

PubMed2002

Absolute comparison of simulated and experimental protein-folding dynamics.

Protein folding is difficult to simulate with classical molecular dynamics. Secondary structure motifs such as alpha-helices and beta-hairpins can form in 0.1-10 micros (ref. 1), whereas small proteins have been shown to fold completely in tens of microseconds. The longest folding simulation to date is a single 1- micro s simulation of the villin headpiece; however, such single runs may miss many features of the folding process as it is a heterogeneous reaction involving an ensemble of transition states. Here, we have used a distributed computing implementation to produce tens of thousands of 5-20-ns trajectories (700 micros) to simulate mutants of the designed mini-protein BBA5. The fast relaxation dynamics these predict were compared with the results of laser temperature-jump experiments. Our computational predictions are in excellent agreement with the experimentally determined mean folding times and equilibrium constants. The rapid folding of BBA5 is due to the swift formation of secondary structure. The convergence of experimentally and computationally accessible timescales will allow the comparison of absolute quantities characterizing in vitro and in silico (computed) protein folding.

Amino Acid Sequence↗

PubMed2002

Native-like mean structure in the unfolded ensemble of small proteins.

The nature of the unfolded state plays a great role in our understanding of proteins. However, accurately studying the unfolded state with computer simulation is difficult, due to its complexity and the great deal of sampling required. Using a supercluster of over 10,000 processors we have performed close to 800 micros of molecular dynamics simulation in atomistic detail of the folded and unfolded states of three polypeptides from a range of structural classes: the all-alpha villin headpiece molecule, the beta hairpin tryptophan zipper, and a designed alpha-beta zinc finger mimic. A comparison between the folded and the unfolded ensembles reveals that, even though virtually none of the individual members of the unfolded ensemble exhibits native-like features, the mean unfolded structure (averaged over the entire unfolded ensemble) has a native-like geometry. This suggests several novel implications for protein folding and structure prediction as well as new interpretations for experiments which find structure in ensemble-averaged measurements.

Amino Acid Sequence↗