Trp-cage is a synthetic 20-residue miniprotein which folds rapidly and spontaneously to a well-defined globular structure more typical of larger proteins. Due to its small size and fast folding, it is an ideal model system for experimental and theoretical investigations of protein folding mechanisms. However, Trp-cage's exact folding mechanism is still a matter of debate. Here we investigate Trp-cage's relaxation dynamics in the amide I' spectral region (1530-1700 cm(-1)) using time-resolved infrared spectroscopy. Residue-specific information was obtained by incorporating an isotopic label ((13)C═(18)O) into the amide carbonyl group of residue Gly11, thereby spectrally isolating an individual 310-helical residue. The folding-unfolding equilibrium is perturbed using a nanosecond temperature-jump (T-jump), and the subsequent re-equilibration is probed by observing the time-dependent vibrational response in the amide I' region. We observe bimodal relaxation kinetics with time constants of 100 ± 10 and 770 ± 40 ns at 322 K, suggesting that the folding involves an intermediate state, the character of which can be determined from the time- and frequency-resolved data. We find that the relaxation dynamics close to the melting temperature involve fast fluctuations in the polyproline II region, whereas the slower process can be attributed to conformational rearrangements due to the global (un)folding transition of the protein. Combined analysis of our T-jump data and molecular dynamics simulations indicates that the formation of a well-defined α-helix precedes the rapid formation of the hydrophobic cage structure, implying a native-like folding intermediate, that mainly differs from the folded conformation in the orientation of the C-terminal polyproline II helix relative to the N-terminal part of the backbone. We find that the main free-energy barrier is positioned between the folding intermediate and the unfolded state ensemble, and that it involves the formation of the α-helix, the 310-helix, and the Asp9-Arg16 salt bridge. Our results suggest that at low temperature (T ≪ Tm) a folding path via formation of α-helical contacts followed by hydrophobic clustering becomes more important.

We have applied the recently developed multiple state transition interface sampling approach to alanine dipeptide in explicit water. We extract the rate constant matrix for configurational changes between each pair of metastable states. The results are comparable with values from previous literature and show that the method is applicable to biomolecular systems.

The lipid composition of cell membranes has increasingly been recognized as playing an important role in the function of various membrane proteins, including G Protein-Coupled Receptors (GPCRs). For instance, experimental and computational evidence has pointed to lipids influencing receptor oligomerization directly, by physically interacting with the receptor, and/or indirectly, by altering the bulk properties of the membrane. While the exact role of oligomerization in the function of class A GPCRs such as the μ-opioid receptor (MOR) is still unclear, insight as to how these receptors oligomerize and the relevance of the lipid environment to this phenomenon is crucial to our understanding of receptor function. To examine the effect of lipids and different MOR conformations on receptor oligomerization we carried out extensive coarse-grained molecular dynamics simulations of crystal structures of inactive and/or activated MOR embedded in an idealized mammalian plasma membrane composed of 63 lipid types asymmetrically distributed across the two leaflets. The results of these simulations point, for the first time, to specific direct and indirect effects of the lipids, as well as the receptor conformation, on the spatio-temporal organization of MOR in the plasma membrane. While sphingomyelin-rich, high-order lipid regions near certain transmembrane (TM) helices of MOR induce an effective long-range attractive force on individual protomers, both long-range lipid order and interface formation are found to be conformation dependent, with a larger number of different interfaces formed by inactive MOR compared to active MOR.

Although protein folding is typically studied in dilute solution, folding in a cell will be affected by interactions with other biomolecules and excluded volume effects. Here, we examine the effect of hydrophobic confinement on folding of the Trp-cage miniprotein. We used replica exchange molecular dynamics simulations to probe the differences between folding in the bulk, on a hydrophobic surface, and confined between two hydrophobic walls. In addition to promotion of helix formation due to reduced conformational entropy of the unfolded state upon confinement, adsorption of Trp-cage to a hydrophobic surface stabilizes intermediate structures not present in the bulk. These new intermediate structures may alter the folding mechanism and kinetics and show the importance of including environmental effects when studying protein folding.

Protein kinases play a key role in regulating cellular processes. Kinase dysfunction can lead to disease, making them an attractive target for drug design. The B-Raf kinase is a key target for the treatment of melanoma since a single mutation (V600E) is found in more than 50% of all malignant melanomas. Despite the importance of B-Raf in melanoma treatment, the molecular mechanism by which the mutation increases kinase activity remains elusive. Since kinases are tightly regulated by a conformational transition between an active and inactive state, which is difficult to capture experimentally, large-scale enhanced-sampling simulations are performed to examine the mechanism by which the V600E mutation enhances the activity of the B-Raf monomer. The results reveal that the mutation has a twofold effect. First, the mutation increases the barrier of the active to inactive transition trapping B-Raf in the active state. The mutation also increases the flexibility of the activation loop which might speed-up the rate-limiting step of phosphorylation. Both effects can be explained by the formation of salt-bridges with the Glu600 residue.

Equation ͑8͒should read:There is an extra T in the first exponential in the published paper.The value of ប was erroneously used instead of Planck's constant, h ͑h =2ប͒, when calculating vibrational entropies using Eq. ͑9͒. The consequence of this error is that the values given for the entropic factors, e ͑−⌬S/k B ͒ , and diffusion constant pre-exponential factors, D 0 , in Tables I, IIb, IIIb, IVb, and V are too large by a factor of 2. This error does not affect the overall trends or conclusions presented in the paper, but the agreement of the pre-exponential factors with experimental data is poorer than previously represented.

No abstract

This article is part of a themed section on Emerging Areas of Opioid Pharmacology. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v175.14/issuetoc.