PubMed Health⌕ Search

Biomedical subjects

Robert B Best

Publications and source records attributed to Robert B Best.

At least 19 recordsLinked to original sources
PubMed2006

Characterization of the residual structure in the unfolded state of the Delta131Delta fragment of staphylococcal nuclease.

The determination of the conformational preferences in unfolded states of proteins constitutes an important challenge in structural biology. We use inter-residue distances estimated from site-directed spin-labeling NMR experimental measurements as ensemble-averaged restraints in all-atom molecular dynamics simulations to characterise the residual structure of the Delta131Delta fragment of staphylococcal nuclease under physiological conditions. Our findings indicate that Delta131Delta under these conditions shows a tendency to form transiently hydrophobic clusters similar to those present in the native state of wild-type staphylococcal nuclease. Only rarely, however, all these interactions are simultaneously realized to generate conformations with an overall native topology.

Computer Simulation↗
PubMed2006

Relation between native ensembles and experimental structures of proteins.

Different experimental structures of the same protein or of proteins with high sequence similarity contain many small variations. Here we construct ensembles of "high-sequence similarity Protein Data Bank" (HSP) structures and consider the extent to which such ensembles represent the structural heterogeneity of the native state in solution. We find that different NMR measurements probing structure and dynamics of given proteins in solution, including order parameters, scalar couplings, and residual dipolar couplings, are remarkably well reproduced by their respective high-sequence similarity Protein Data Bank ensembles; moreover, we show that the effects of uncertainties in structure determination are insufficient to explain the results. These results highlight the importance of accounting for native-state protein dynamics in making comparisons with ensemble-averaged experimental data and suggest that even a modest number of structures of a protein determined under different conditions, or with small variations in sequence, capture a representative subset of the true native-state ensemble.

Animals↗
PubMed2006

Diffusive model of protein folding dynamics with Kramers turnover in rate.

We study the folding kinetics of a three-helix bundle protein using a coarse polymer model. The folding dynamics can be accurately represented by one-dimensional diffusion along a reaction coordinate selected to capture the transition state. By varying the solvent friction, we show that position-dependent diffusion coefficients are determined by microscopic transitions on a rough energy landscape. A maximum in the folding rate at intermediate friction is explained by "Kramers turnover" in these microscopic dynamics that modulates the rate via the diffusion coefficient; overall folding remains diffusive even close to zero friction. For water friction, we find that the "attempt frequency" (or "speed limit") in a Kramers model of folding is about 2 micros-1, with an activation barrier of about 2kBT, and a folding transition path duration of approximately equal to 100 ns, 2 orders of magnitude less than the folding time of approximately equal to 10 micros.

Diffusion↗
PubMed2006

Structural comparison of the two alternative transition states for folding of TI I27.

TI I27, a beta-sandwich domain from the human muscle protein titin, has been shown to fold via two alternative pathways, which correspond to a change in the folding mechanism. Under physiological conditions, TI I27 folds by a classical nucleation-condensation mechanism (diffuse transition state), whereas at extreme conditions of temperature and denaturant it switches to having a polarized transition state. We have used experimental Phi-values as restraints in ensemble-averaged molecular dynamics simulations to determine the ensembles of structures representing the two transition states. The comparison of these ensembles indicates that when native interactions are substantially weakened, a protein may still be able to fold if it can access an alternative transition state characterized by a much larger entropic contribution. Analysis of the probability distribution of Phi-values derived from ensemble averaged simulations, enables us to identify residues that form contacts in some members of the ensemble but not in others illustrating that many interactions present in transition states are not strictly required for the successful completion of the folding process.

Computer Simulation↗
PubMed2006

Structural interpretation of hydrogen exchange protection factors in proteins: characterization of the native state fluctuations of CI2.

Protection factors obtained from equilibrium hydrogen exchange experiments are an important source of structural information on both native and nonnative states of proteins. We present a method for determining ensembles of protein structures by using hydrogen exchange data as restraints in molecular dynamics simulations in conjunction with an empirical force-field. The method is applied to determine the ensemble of structures representing the native state of chymotrypsin inhibitor 2 (CI2), including the rare, large fluctuations responsible for hydrogen exchange.

Chymotrypsin↗
PubMed2005

Reaction coordinates and rates from transition paths.

The molecular mechanism of a reaction in solution is reflected in its transition-state ensemble and transition paths. We use a Bayesian formula relating the equilibrium and transition-path ensembles to identify transition states, rank reaction coordinates, and estimate rate coefficients. We also introduce a variational procedure to optimize reaction coordinates. The theory is illustrated with applications to protein folding and the dipole reorientation of an ordered water chain inside a carbon nanotube. To describe the folding of a simple model of a three-helix bundle protein, we variationally optimize the weights of a projection onto the matrix of native and nonnative amino acid contacts. The resulting one-dimensional reaction coordinate captures the folding transition state, with formation and packing of helix 2 and 3 constituting the bottleneck for folding.

Adenosine Triphosphate↗
PubMed2005

What contributions to protein side-chain dynamics are probed by NMR experiments? A molecular dynamics simulation analysis.

Molecular dynamics simulations of the structurally homologous proteins TNfn3 and FNfn10 have been used to investigate the contributions to side-chain dynamics measured by NMR relaxation experiments. The results reproduce the variation in core side-chain dynamics observed by NMR and highlight the relevance of anharmonic motion and transitions between local minima for explaining NMR side-chain order parameters. A method is described for calculating converged order parameters by use of replica exchange molecular dynamics in conjunction with an implicit solvent model. These simulations allow the influence of various factors, such as the flexibility of side-chains and their free volume, on the mobility to be tested by perturbing the system. Deletion mutations are found to have the largest effect on the more densely packed FNfn10. Some counterintuitive effects are seen, such as an increase in order parameters close to deletion mutation sites, but these can be rationalized in terms of direct interactions with the modified side-chains. A statistical analysis of published order parameters supports the conclusions drawn from the simulations.

Amino Acids↗
PubMed2005

Simultaneous determination of protein structure and dynamics.

We present a protocol for the experimental determination of ensembles of protein conformations that represent simultaneously the native structure and its associated dynamics. The procedure combines the strengths of nuclear magnetic resonance spectroscopy--for obtaining experimental information at the atomic level about the structural and dynamical features of proteins--with the ability of molecular dynamics simulations to explore a wide range of protein conformations. We illustrate the method for human ubiquitin in solution and find that there is considerable conformational heterogeneity throughout the protein structure. The interior atoms of the protein are tightly packed in each individual conformation that contributes to the ensemble but their overall behaviour can be described as having a significant degree of liquid-like character. The protocol is completely general and should lead to significant advances in our ability to understand and utilize the structures of native proteins.

Computer Simulation↗
PubMed2005

Interpreting dynamically-averaged scalar couplings in proteins.

The experimental determination of scalar three-bond coupling constants represents a powerful method to probe both the structure and dynamics of proteins. The detailed structural interpretation of such coupling constants is usually based on Karplus relationships, which allow the measured couplings to be related to the torsion angles of the molecules. As the measured couplings are sensitive to thermal fluctuations, the parameters in the Karplus relationships are better derived from ensembles representing the distributions of dihedral angles present in solution, rather than from single conformations. We present a method to derive such parameters that uses ensembles of conformations determined through dynamic-ensemble refinement--a method that provides structural ensembles that simultaneously represent both the structure and the associated dynamics of a protein.

Models, Molecular↗
PubMed2005

Slow protein conformational dynamics from multiple experimental structures: the helix/sheet transition of arc repressor.

Conformational transitions underlie the function of many biomolecular systems. Resolving intermediate structural changes, however, is challenging for both experiments and all-atom simulations because the duration of transitions is short relative to the lifetime of the stable species. Simplified descriptions based on a single experimental structure, such as elastic network models or Gō models, are not immediately applicable. Here, we develop a general method that combines multiple coarse-grained models to capture slow conformational transitions. Individually, each model describes one of the experimental structures; together, they approximate the complete energy surface. We demonstrate the method for the helix-to-sheet transition in Arc repressor N11L. We find that the transition involves the partial unfolding of the switch region, and rapid refolding into the alternate structure. Transient local unfolding is consistent with the low hydrogen exchange protection factors of the switch region. Also in agreement with experiment, the isomerization occurs independently of the global folding/dimerization transition.

Dimerization↗
PubMed2004

A mutant chaperonin with rearranged inter-ring electrostatic contacts and temperature-sensitive dissociation.

The chaperonin GroEL assists protein folding through ATP-dependent, cooperative movements that alternately create folding chambers in its two rings. The substitution E461K at the interface between these two rings causes temperature-sensitive, defective protein folding in Escherichia coli. To understand the molecular defect, we have examined the mutant chaperonin by cryo-EM. The normal out-of-register alignment of contacts between subunits of opposing wild-type rings is changed in E461K to an in-register one. This is associated with loss of cooperativity in ATP binding and hydrolysis. Consistent with the loss of negative cooperativity between rings, the cochaperonin GroES binds simultaneously to both E461K rings. These GroES-bound structures were unstable at higher temperature, dissociating into complexes of single E461K rings associated with GroES. Lacking the allosteric signal from the opposite ring, these complexes cannot release their GroES and become trapped, dead-end states.

Adenosine Triphosphatases↗
PubMed2004

Determination of protein structures consistent with NMR order parameters.

Order parameters obtained from NMR experiments characterize distributions of bond vector orientations. Their interpretation, however, usually requires the assumption of a particular motional model. We propose a multiple-copy simulation method in which the experimental order parameters are used as restraints in conjunction with a standard molecular force field. The latter effectively acts as a sophisticated motional model, allowing ensembles of structures consistent with the experimental order parameters to be determined.

Cell Adhesion Molecules↗
PubMed2004

The origin of protein sidechain order parameter distributions.

Previous work by Wand et al. (Nature 2001, 411, 501-504) showed that the NMR order parameters characterizing the amplitude of motion of protein side chains seemed to form a multimodal distribution. At the time, no detailed explanation of this at the molecular level was offered, yet three "classes" of motion were inferred. We have analyzed a larger published data set and found that, although the distribution is multimodal, the evidence for three classes is weak. More significantly, we have been able to provide a simple physical explanation for the distributions based on the results of molecular dynamics simulations. This result will aid in the interpretation of data from NMR dynamics experiments.

Amino Acids↗
PubMed2004

Hydrophobic core fluidity of homologous protein domains: relation of side-chain dynamics to core composition and packing.

The side-chain dynamics of methyl groups in two structurally related proteins from the fibronectin type III (fnIII) superfamily, the third fnIII domain from human tenascin (TNfn3) and the tenth fnIII domain from human fibronectin (FNfn10), have been studied by NMR spectroscopy. Side-chain order parameters reveal that the hydrophobic cores of the two proteins have substantially different mobilities. The core of TNfn3 is very dynamic, with exceptionally low order parameters for the most deeply buried residues, while that of FNfn10 is more like those of other proteins which have been studied with this technique, having a relatively rigid core with uniformly distributed dynamics. The unusually dynamic core of TNfn3 appears to be related to its amino acid composition, which makes it more fluid-like. A further explanation for the mobility of the TNfn3 core may be found in the negative correlation between the order parameter and excess packing volume, which shows that the core of TNfn3 is less densely packed and consequently has lower methyl order parameters for its buried residues. Rotameric transitions, presumably facilitated by the lower packing density, appear to make an important contribution to lowering the order parameters, and have been probed by measuring three-bond scalar couplings. Overall, although backbone dynamics is generally similar for proteins with the same topology on a fast time scale (picoseconds to nanoseconds), this study shows that a single fold can accommodate a wide variation in the dynamic properties of its core.

Amino Acid Sequence↗
PubMed2003

Mechanical unfolding of a titin Ig domain: structure of transition state revealed by combining atomic force microscopy, protein engineering and molecular dynamics simulations.

Titin I27 shows a high resistance to unfolding when subject to external force. To investigate the molecular basis of this mechanical stability, protein engineering Phi-value analysis has been combined with atomic force microscopy to investigate the structure of the barrier to forced unfolding. The results indicate that the transition state for forced unfolding is significantly structured, since highly destabilising mutations in the core do not affect the force required to unfold the protein. As has been shown before, mechanical strength lies in the region of the A' and G-strands but, contrary to previous suggestions, the results indicate clearly that side-chain interactions play a significant role in maintaining mechanical stability. Since Phi-values calculated from molecular dynamics simulations are the same as those determined experimentally, we can, with confidence, use the molecular dynamics simulations to analyse the structure of the transition state in detail, and are able to show loss of interactions between the A' and G-strands with associated A-B and E-F loops in the transition state. The key event is not a simple case of loss of hydrogen bonding interactions between the A' and G-strands alone. Comparison with Phi-values from traditional folding studies shows differences between the force and "no-force" transition states but, nevertheless, the region important for kinetic stability is the same in both cases. This explains the correspondence between hierarchy of kinetic stability (measured in stopped-flow denaturant studies) and mechanical strength in these titin domains.

Computer Simulation↗
PubMed2003

Hidden complexity in the mechanical properties of titin.

Individual molecules of the giant protein titin span the A-bands and I-bands that make up striated muscle. The I-band region of titin is responsible for passive elasticity in such muscle, and contains tandem arrays of immunoglobulin domains. One such domain (I27) has been investigated extensively, using dynamic force spectroscopy and simulation. However, the relevance of these studies to the behaviour of the protein under physiological conditions was not established. Force studies reveal a lengthening of I27 without complete unfolding, forming a stable intermediate that has been suggested to be an important component of titin elasticity. To develop a more complete picture of the forced unfolding pathway, we use mutant titins--certain mutations allow the role of the partly unfolded intermediate to be investigated in more depth. Here we show that, under physiological forces, the partly unfolded intermediate does not contribute to mechanical strength. We also propose a unified forced unfolding model of all I27 analogues studied, and conclude that I27 can withstand higher forces in muscle than was predicted previously.

Biomechanical Phenomena↗
PubMed2002

Mechanical unfolding of a titin Ig domain: structure of unfolding intermediate revealed by combining AFM, molecular dynamics simulations, NMR and protein engineering.

The mechanical unfolding of an immunoglobulin domain from the human muscle protein titin (TI I27) has been shown to proceed via a metastable intermediate in which the A-strand is detached. The structure and properties of this intermediate are characterised in this study. A conservative destabilising mutation in the A-strand has no effect on the unfolding force, nor the dependence of the unfolding force on the pulling speed, indicating that the unfolding forces measured in an AFM experiment are those required for the unfolding of the intermediate and not the native state. A mutant of TI I27 with the A-strand deleted (TI I27-A) is studied by NMR and standard biophysical techniques, combined with protein engineering. Molecular dynamics simulations show TI I27-A to be a good model for the intermediate. It has a structure very similar to the native state, and is surprisingly stable. Comparison with a Phi-value analysis of the unfolding pathway clearly shows that the protein unfolds by a different pathway under an applied force than on addition of denaturant.

Computer Simulation↗