Explore the vast range of titles available.

It is currently the consensus belief that protective osmolytes such as trimethylamine N-oxide (TMAO) favor protein folding by being excluded from the vicinity of a protein, whereas denaturing osmolytes such as urea lead to protein unfolding by strongly binding to the surface. Despite there being consensus on how TMAO and urea affect proteins as a whole, very little is known as to their effects on the individual mechanisms responsible for protein structure formation, especially hydrophobic association. In the present study, we use single-molecule atomic force microscopy and molecular dynamics simulations to investigate the effects of TMAO and urea on the unfolding of the hydrophobic homopolymer polystyrene. Incorporated with interfacial energy measurements, our results show that TMAO and urea act on polystyrene as a protectant and a denaturant, respectively, while complying with Tanford–Wyman preferential binding theory. We provide a molecular explanation suggesting that TMAO molecules have a greater thermodynamic binding affinity with the collapsed conformation of polystyrene than with the extended conformation, while the reverse is true for urea molecules. Results presented here from both experiment and simulation are in line with earlier predictions on a model Lennard–Jones polymer while also demonstrating the distinction in the mechanism of osmolyte action between protein and hydrophobic polymer. This marks, to our knowledge, the first experimental observation of TMAO-induced hydrophobic collapse in a ternary aqueous system.

Single-molecule force spectroscopies are remarkable tools for studying protein folding and unfolding, but force unfolding explores protein configurations that are potentially very different from the ones traditionally explored in chemical or thermal denaturation. Understanding these differences is crucial because such configurations serve as starting points of folding studies, and thus can affect both the folding mechanism and the kinetics. Here we provide a detailed comparison of both chemically induced and force-induced unfolded state ensembles of ubiquitin based on extensive, all-atom simulations of the protein either extended by force or denatured by urea. As expected, the respective unfolded states are very different on a macromolecular scale, being fully extended under force with no contacts and partially extended in urea with many nonnative contacts. The amount of residual secondary structure also differs: A significant population of α-helices is found in chemically denatured configurations but such helices are absent under force, except at the lowest applied force of 30 pN where short helices form transiently. We see that typical-size helices are unstable above this force, and β-sheets cannot form. More surprisingly, we observe striking differences in the backbone dihedral angle distributions for the protein unfolded under force and the one unfolded by denaturant. A simple model based on the dialanine peptide is shown to not only provide an explanation for these striking differences but also illustrates how the force dependence of the protein dihedral angle distributions give rise to the worm-like chain behavior of the chain upon force.

The thermodynamic properties and phase behavior of water in confined regions can vary significantly from that observed in the bulk. This is particularly true for systems in which the confinement is on the molecular-length scale. In this study, we use molecular dynamics simulations and a powerful solvent analysis technique based on inhomogenous solvation theory to investigate the properties of water molecules that solvate the confined regions of protein active sites. Our simulations and analysis indicate that the solvation of protein active sites that are characterized by hydrophobic enclosure and correlated hydrogen bonds induce atypical entropic and enthalpic penalties of hydration. These penalties apparently stabilize the protein-ligand complex with respect to the independently solvated ligand and protein, which leads to enhanced binding affinities. Our analysis elucidates several challenging cases, including the super affinity of the streptavidin-biotin system.

We performed molecular dynamics simulations of the collapse of a two-domain protein, the BphC enzyme, into a globular structure to examine how water molecules mediate hydrophobic collapse of proteins. In the interdomain region, liquid water persists with a density 10 to 15% lower than in the bulk, even at small domain separations. Water depletion and hydrophobic collapse occur on a nanosecond time scale, which is two orders of magnitude slower than that found in the collapse of idealized paraffin-like plates. When the electrostatic protein-water forces are turned off, a dewetting transition occurs in the interdomain region and the collapse speeds up by more than an order of magnitude. When attractive van der Waals forces are turned off as well, the dewetting in the interdomain region is more profound, and the collapse is even faster.

The folding free energy landscape of the C-terminal β hairpin of protein G has been explored in this study with explicit solvent under periodic boundary condition and Oplsaa force field. A highly parallel replica exchange method that combines molecular dynamics trajectories with a temperature exchange Monte Carlo process is used for sampling with the help of a new efficient algorithm P3ME/RESPA. The simulation results show that the hydrophobic core and the β strand hydrogen bond form at roughly the same time. The free energy landscape with respect to various reaction coordinates is found to be rugged at low temperatures and becomes a smooth funnel-like landscape at about 360 K. In contrast to some very recent studies, no significant helical content has been found in our simulation at all temperatures studied. The β hairpin population and hydrogen-bond probability are in reasonable agreement with the experiment at biological temperature, but both decay more slowly than the experiment with temperature.

The elastic restoring force of tissues must be able to operate over the very wide range of loading rates experienced by living organisms. It is surprising that even the fastest events involving animal muscle tissues do not surpass a few hundred hertz. We propose that this limit is set in part by the elastic dynamics of tethered proteins extending and relaxing under a changing load. Here we study the elastic dynamics of tethered proteins using a fast force spectrometer with sub-millisecond time resolution, combined with Brownian and Molecular Dynamics simulations. We show that the act of tethering a polypeptide to an object, an inseparable part of protein elasticity in vivo and in experimental setups, greatly reduces the attempt frequency with which the protein samples its free energy. Indeed, our data shows that a tethered polypeptide can traverse its free-energy landscape with a surprisingly low effective diffusion coefficient D ₑff ∼ 1,200 nm ²/s. By contrast, our Molecular Dynamics simulations show that diffusion of an isolated protein under force occurs at D ₑff ∼ 10 ⁸ nm ²/s. This discrepancy is attributed to the drag force caused by the tethering object. From the physiological time scales of tissue elasticity, we calculate that tethered elastic proteins equilibrate in vivo with D ₑff ∼ 10 ⁴–10 ⁶ nm ²/s which is two to four orders magnitude smaller than the values measured for untethered proteins in bulk.

The widely used Arrhenius equation describes the kinetics of simple two-state reactions, with the implicit assumption of a single transition state with a well-defined activation energy barrier ΔE, as the rate-limiting step. However, it has become increasingly clear that the saddle point of the free-energy surface in most reactions is populated by ensembles of conformations, leading to nonexponential kinetics. Here we present a theory that generalizes the Arrhenius equation to include static disorder of conformational degrees of freedom as a function of an external perturbation to fully account for a diverse set of transition states. The effect of a perturbation on static disorder is best examined at the single-molecule level. Here we use force-clamp spectroscopy to study the nonexponential kinetics of single ubiquitin proteins unfolding under force. We find that the measured variance in ΔE shows both force-dependent and independent components, where the force-dependent component scales with F², in excellent agreement with our theory. Our study illustrates a novel adaptation of the classical Arrhenius equation that accounts for the microscopic origins of nonexponential kinetics, which are essential in understanding the rapidly growing body of single-molecule data.

We propose a mechanism, based on a ≥10-μs molecular dynamics simulation, for the surprising misfolding of hen egg-white lysozyme caused by a single mutation (W62G). Our simulations of the wild-type and mutant lysozymes in 8 M urea solution at biological temperature (with both pH 2 and 7) reveal that the mutant structure is much less stable than that of the wild type, with the mutant showing larger fluctuations and less native-like contacts. Analysis of local contacts reveals that the Trp-62 residue is the key to a cooperative long-range interaction within the wild type, where it acts like a bridge between two neighboring basic residues. Thus, a native-like cluster or nucleation site can form near these residues in the wild type but not in the mutant. The time evolution of the secondary structure also exhibits a quicker loss of the β-sheets in the mutant than in the wild type, whereas some of the α-helices persist during the entire simulation in both the wild type and the mutant in 8 M urea (even though the tertiary structures are basically all gone). These findings, while supporting the general conclusions of a recent experimental study by Dobson and coworkers [Klein-Seetharam J, Oikama M, Grimshaw SB, Wirmer J, Duchardt E, Ueda T, Imoto T, Smith LJ, Dobson CM, Schwalbe H (2002) Science 295:1719-1722], provide a detailed but different molecular picture of the misfolding mechanism.

We unfold and extend single proteins at a high force and then linearly relax the force to probe their collapse mechanisms. We observe a large variability in the extent of their recoil. Although chain entropy makes a small contribution, we show that the observed variability results from hydrophobic interactions with randomly varying magnitude from protein to protein. This collapse mechanism is common to highly extended proteins, including nonfolding elastomeric proteins like PEVK from titin. Our observations explain the puzzling differences between the folding behavior of highly extended proteins, from those folding after chemical or thermal denaturation. Probing the collapse of highly extended proteins with force spectroscopy allows separation of the different driving forces in protein folding.

Thioredoxins (Trxs) reduce disulfide bonds via a Michaelis-Menten mechanism. Upon substrate stretching at high forces, an S N 2 reaction can be used by bacterial Trxs. A third mechanism, single-electron transfer, is now revealed in Trxs of either bacterial or eukaryotic origin, and is correlated with the depth of the Trx substrate-binding groove. Thioredoxins (Trxs) are oxidoreductase enzymes, present in all organisms, that catalyze the reduction of disulfide bonds in proteins. By applying a calibrated force to a substrate disulfide, the chemical mechanisms of Trx catalysis can be examined in detail at the single-molecule level. Here we use single-molecule force-clamp spectroscopy to explore the chemical evolution of Trx catalysis by probing the chemistry of eight different Trx enzymes. All Trxs show a characteristic Michaelis-Menten mechanism that is detected when the disulfide bond is stretched at low forces, but at high forces, two different chemical behaviors distinguish bacterial-origin from eukaryotic-origin Trxs. Eukaryotic-origin Trxs reduce disulfide bonds through a single-electron transfer reaction (SET), whereas bacterial-origin Trxs show both nucleophilic substitution (S N 2) and SET reactions. A computational analysis of Trx structures identifies the evolution of the binding groove as an important factor controlling the chemistry of Trx catalysis.