Computer simulations of the refolding of sperm whale apomyoglobin from high-temperature denaturated state

The refolding mechanism of apomyoglobin (apoMb) subsequent to high-temperature unfolding has been examined using computer simulations with atomic level detail. The folding of this protein has been extensively studied experimentally, providing a large database of folding parameters which can be probed using simulations. In the present study, 4-folding trajectories of apoMb were computed starting from coiled structures. A crystal structure of sperm whale myoglobin taken from the Protein Data Bank was used to construct the final native conformation by removal of the heme group followed by energy optimization. The initial unfolded conformations were obtained from high-temperature molecular dynamics simulations. Room-temperature refolding trajectories at neutral pH were obtained using the stochastic difference equation in length algorithm. The folding trajectories were compared with experimental results and two previous molecular dynamics studies at low pH. In contrast to the previous simulations, an extended intermediate with large helical content was not observed. In the present study, a structural collapse occurs without formation of helices or native contacts. Once the protein structure is more compact (radius of gyration<18 A) secondary and tertiary structures appear. These results suggest that apoMb follows a different folding pathway after high-temperature denaturation.     

Solvation dynamics of a protein in the pre molten globule state

The nature of solvent molecules around proteins in native and different non-native states is crucial for understanding the protein folding problem. We have characterized two compact denatured states of glutaminyl-tRNA synthetase (GlnRS) under equilibrium conditions in the presence of a naturally occurring osmolyte, l-glutamate. The solvation dynamics of the compact denatured states and the fully unfolded state has been studied using a covalently attached probe, acrylodan, near the active site. The solvation dynamics progressively becomes faster as the protein goes from the native to the molten globule to the pre molten globule to the fully unfolded state. Anisotropy decay measurements suggest that the pre-molten-globule intermediate is more flexible than the molten globule although the secondary structure is largely similar. Dynamic light scattering studies reveal that both the compact denatured states are aggregated under the measurement conditions. The implications of solvation dynamics in aggregated compact denatured states have been discussed.     

蛋白折叠中的暂态结构与表面水分子慢尺度动力学

蛋白表面水的慢尺度动力学行为往往被认为与蛋白的结构稳定性、功能以及折叠过程有关, 但在分子水平上, 还不清楚水分子的慢尺度动力学如何参与蛋白折叠过程. 以Trp-cage蛋白作为个案, 本文利用40条100 ns(总长4 μs)的全原子分子动力学轨迹,分析了蛋白折叠过程中蛋白表面水分子的停留行为,并探究影响蛋白表面水分子慢尺度行为的微观因素. 结果发现, 即使在蛋白折叠过程中蛋白拓扑结构变化很大, 残基之间也会形成稳定的局部暂态结构. 这些结构为水分子提供饱和、稳定的氢键, 通过与水分子之间的极性相互作用, 以及凹形的几何结构, 约束水分子长时停留, 我们称之为“停留中心”. 停留中心的形成是引起水分子慢尺度行为的重要因素. 另外, 停留中心的分布与蛋白折叠的进程有密切关系, 特别地, 在折叠轨迹中, 疏水核周围的残基组成了一个主要的停留中心. 研究结果不但有助于解释水分子慢尺度特征行为的来源, 还可以为实验中通过研究水分子在蛋白附近的慢尺度行为, 揭示蛋白折叠过程中的关键步骤提供一些启发.     

Dynamic nuclear polarization enhanced nuclear magnetic resonance and electron spin resonance studies of hydration and local water dynamics in micelle and vesicle assemblies

We present a unique analysis tool for the selective detection of local water inside soft molecular assemblies (hydrophobic cores, vesicular bilayers, and micellar structures) suspended in bulk water. Through the use of dynamic nuclear polarization (DNP), the (1)H NMR signal of water is amplified, as it interacts with stable radicals that possess approximately 658 times higher spin polarization. We utilized stable nitroxide radicals covalently attached along the hydrophobic tail of stearic acid molecules that incorporate themselves into surfactant-based micelle or vesicle structures. Here, we present a study of local water content and fluid viscosity inside oleate micelles and vesicles and Triton X-100 micelles to serve as model systems for soft molecular assemblies. This approach is unique because the amplification of the NMR signal is performed in bulk solution and under ambient conditions with site-specific spin labels that only detect the water that is directly interacting with the localized spin labels. Continuous wave (cw) electron spin resonance (ESR) analysis provides rotational dynamics of the spin-labeled molecular chain segments and local polarity parameters that can be related to hydration properties, whereas we show that DNP-enhanced (1)H NMR analysis of fluid samples directly provides translational water dynamics and permeability of the local environment probed by the spin label. Our technique therefore has the potential to become a powerful analysis tool, complementary to cw ESR, to study hydration characteristics of surfactant assemblies, lipid bilayers, or protein aggregates, where water dynamics is a key parameter of their structure and function. In this study, we find that there is significant penetration of water inside the oleate micelles with a higher average local water viscosity ( approximately 1.8 cP) than in bulk water, and Triton X-100 micelles and oleate vesicle bilayers mostly exclude water while allowing for considerable surfactant chain motion and measurable water permeation through the soft structure.     

Parametrization, molecular dynamics simulation, and calculation of electron spin resonance spectra of a nitroxide spin label on a polyalanine alpha-helix

The nitroxide spin label 1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl-methanethiosulfonate (MTSSL), commonly used in site-directed spin labeling of proteins, is studied with molecular dynamics (MD) simulations. After developing force field parameters for the nitroxide moiety and the spin label linker, we simulate MTSSL attached to a polyalanine alpha-helix in explicit solvent to elucidate the factors affecting its conformational dynamics. Electron spin resonance spectra at 9 and 250 GHz are simulated in the time domain using the MD trajectories and including global rotational diffusion appropriate for the tumbling of T4 Lysozyme in solution. Analysis of the MD simulations reveals the presence of significant hydrophobic interactions of the spin label with the alanine side chains.     

Molecular dynamics for surfactant-assisted protein refolding

Surfactants are widely used to refold recombinant proteins that are produced as inclusion bodies in E. Coli. However, the microscopic details of the surfactant-assisted protein refolding processes are yet to be uncovered. In the present work, the authors aim to provide insights into the effect of hydrophobic interactions of a denatured protein with surfactant molecules on the refolding kinetics and equilibrium by using the Langevin dynamics for coarse-grained models. The authors have investigated the folding behavior of a beta-barrel protein in the presence of surfactants of different hydrophobicities and concentrations. It is shown that the protein folding process follows a "collapse-rearrangement" mechanism, i.e., the denatured protein first falls into a collapsed state before acquiring the native conformation. In comparison with the protein folding without surfactants, the protein-surfactant hydrophobic interactions promote the collapse of a denatured protein and, consequently, the formation of a hydrophobic core. However, the surfactants must be released from the hydrophobic core during the rearrangement step, in which the native conformation is formed. The simulation results can be qualitatively reproduced by experiments.     

Selective High‐Resolution Detection of Membrane Protein–Ligand Interaction in Native Membranes Using Trityl–Nitroxide PELDOR

The orchestrated interaction of transmembrane proteins with other molecules mediates several crucial biological processes. Detergent solubilization may significantly alter or even abolish such hetero‐oligomeric interactions, which makes observing them at high resolution in their native environment technically challenging. Dipolar electron paramagnetic resonance (EPR) techniques such as pulsed electro–electron double resonance (PELDOR) can provide very precise distances within biomolecules. To concurrently determine the inter‐subunit interaction and the intra‐subunit conformational changes in hetero‐oligomeric complexes, a combination of different spin labels is required. Orthogonal spin labeling using a triarylmethyl (TAM) label in combination with a nitroxide label is used to detect protein–ligand interactions in native lipid bilayers. This approach provides a higher sensitivity and total selectivity and will greatly facilitate the investigation of multimeric transmembrane complexes employing different spin labels in the native lipid environment.     

A general approach for prediction of motional EPR spectra from Molecular Dynamics (MD) simulations: application to spin labelled protein

A general approach for the prediction of EPR spectra directly and completely from single dynamical trajectories generated from Molecular Dynamics (MD) simulations is described. The approach is applicable to an arbitrary system of electron and nuclear spins described by a general form of the spin-Hamiltonian for the entire motional range. It is shown that for a reliable simulation of motional EPR spectra only a single truncated dynamical trajectory generated until the point when correlation functions of rotational dynamics are completely relaxed is required. The simulation algorithm is based on a combination of the propagation of the spin density matrix in the Liouville space for this initial time interval and the use of well defined parameters calculated entirely from the dynamical trajectory for prediction of the evolution of the spin density matrix at longer times. A new approach is illustrated with the application to a nitroxide spin label MTSL attached to the protein sperm whale myoglobin. It is shown that simulation of the EPR spectrum, which is in excellent agreement with experiment, can be achieved from a single MD trajectory. Calculations reveal the complex nature of the dynamics of a spin label which is a superposition of the fast librational motions within dihedral states, of slow rotameric dynamics among different conformational states of the nitroxide tether and of the slow rotational diffusion of the protein itself. The significance of the slow rotameric dynamics of the nitroxide tether on the overall shape of the EPR spectrum is analysed and discussed.     

Water-protein dynamic coupling and new opportunities for probing it at low to physiological temperatures in aqueous solutions

Both the structure and dynamics of biomolecules are known to be essential for their biological function. In the dehydrated state, the function of biomolecules, such as proteins, is severely impeded, so hydration is required for bioactivity. The dynamics of the hydrated biomolecules and their hydration water are related - but how closely? The problem involves several layers of complexity. Even for water in the bulk state, the contribution from various dynamic components to the overall dynamics is not fully understood. In biological systems, the effects of confinement on the hydration water further complicate the picture. Even if the various components of the hydration water dynamics are properly understood, which of them are coupled to the protein dynamics, and how? The studies of protein dynamics over the wide temperature range, from physiological to low temperatures, provide some answers to these question. At low temperatures, both the protein and its hydration water behave as solids, with only vibrational degrees of freedom. As the temperature is increased, non-vibrational dynamic components start contributing to the measurable dynamics and eventually become dominant at physiological temperatures. Thus, the temperature dependence of the dynamics of protein and its hydration water may allow probing various dynamic components separately. In order to suppress the water freezing, the low-temperature studies of protein rely on either low-hydrated samples (essentially, hydrated protein powders), or cryo-protective solutions. Both approaches introduce the hydration environments not characteristic of the protein environments in living systems, which are typically aqueous protein solutions of various concentrations. In this paper, we discuss the coupling between the dynamic components of the protein and its hydration water by critical examining of the existing literature, and then propose that proteins can be studied in an aqueous solution that is remarkably similar in its dynamic properties to pure water, yet does not freeze down to about 200 K, even in the bulk form. The first experiment of this kind using quasielastic neutron scattering is discussed, and more experiments are proposed.     

Dynamics of water in the hydration layer of a partially unfolded structure of the protein HP-36

Atomistic molecular dynamics simulations of the folded native structure and a partially unfolded molten globule structure of the protein villin headpiece subdomain or HP-36 have been carried out with explicit solvent to explore the effects of unfolding on the dynamical behavior of water present in the hydration layers of different segments (three alpha-helices) of the protein. The calculations revealed that the unfolding of helix-2 influences the translational and rotational motions of water present in the hydration layers of the three helices in a heterogeneous manner. It is observed that a correlation exists between the unfolding of helix-2 and the microscopic kinetics of protein-water hydrogen bonds formed by its residues. This in turn has an influence on the rigidity of the hydration layers of the helices in the unfolded structure versus that in the folded native structure. These results should provide a microscopic explanation to recent solvation dynamics experiments on folded native and unfolded structures of proteins.     

The relative helix and hydrogen bond stability in the B domain of protein A as revealed by integrated tempering sampling molecular dynamics simulation

Molecular dynamics simulations using the integrated tempering sampling method were performed for the folding of wild-type B domain of protein A (BdpA). Starting from random and stretched structures, these simulations allow us to fold this protein into the native-like structure frequently, achieving very small backbone (1.7 A?) and all heavy-atom root-mean-square deviation (2.6 A?). Therefore, the method used here increases the efficiency of configuration sampling and thermodynamics characterization by molecular dynamics simulation. Although inconsistency exists between the calculation and experiments for the absolute stabilities, as a limitation of the force field parameters, the calculated order of helix stability (H3 > H2 > H1) is consistent with that determined by experiments for individual separate helices. The lowest free energy folding pathway of BdpA was found to start with a barrierless and non-cooperative structural collapse from the entirely extended (E) state, which leads to a physiologically unfolded (P) state consisting of multiple stable structures with few native inter-helical hydrophobic interactions formed. In the P state, only H3 is fully structured. The final formation of H1 (and to a lesser extent, H2) in the folded (F) state requires the packing of the inter-helical hydrophobic contacts. In addition, it was found that stabilities of backbone hydrogen bonds are significantly affected by their positions relative to the inter-helical hydrophobic core. As temperature increases, the stability of the hydrogen bonds exposed to the solvent tends to increase while that of the hydrogen bonds buried within the hydrophobic core decreases. Finally, we discuss implications of this study on the general folding mechanism of proteins.     

High-field EPR and ESEEM investigation of the nitrogen quadrupole interaction of nitroxide spin labels in disordered solids: toward differentiation between polarity and proticity matrix effects on protein function

The combination of high-field electron paramagnetic resonance (EPR) with site-directed spin labeling (SDSL) techniques employing nitroxide radicals has turned out to be particularly powerful in revealing subtle changes of the polarity and proticity profiles in proteins enbedded in membranes. This information can be obtained by orientation-selective high-field EPR resolving principal components of the nitroxide Zeeman (g) and hyperfine ( A) tensors of the spin labels attached to specific molecular sites. In contrast to the g- and A-tensors, the (14)N ( I = 1) quadrupole interaction tensor of the nitroxide spin label has not been exploited in EPR for probing effects of the microenvironment of functional protein sites. In this work it is shown that the W-band (95 GHz) high-field electron spin echo envelope modulation (ESEEM) method is well suited for determining with high accuracy the (14)N quadrupole tensor principal components of a nitroxide spin label in disordered frozen solution. By W-band ESEEM the quadrupole components of a five-ring pyrroline-type nitroxide radical in glassy ortho-terphenyl and glycerol solutions have been determined. This radical is the headgroup of the MTS spin label widely used in SDSL protein studies. By DFT calulations and W-band ESEEM experiments it is demonstrated that the Q(yy) value is especially sensitive to the proticity and polarity of the nitroxide environment in H-bonding and nonbonding situations. The quadrupole tensor is shown to be rather insensitive to structural variations of the nitroxide label itself. When using Q(yy) as a testing probe of the environment, its ruggedness toward temperature changes represents an important advantage over the g xx and A(zz) parameters which are usually employed for probing matrix effects on the spin labeled molecular site. Thus, beyond measurenments of g xx and A(zz) of spin labeled protein sites in disordered solids, W-band high-field ESEEM studies of (14)N quadrupole interactions open a new avenue to reliably probe subtle environmental effects on the electronic structure. This is a significant step forward on the way to differentiate between effects from matrix polarity and hydrogen-bond formation.     

The equilibrium properties and folding kinetics of an all-atom Go model of the Trp-cage

The ultrafast-folding 20-residue Trp-cage protein is quickly becoming a new benchmark for molecular dynamics studies. Already several all-atom simulations have probed its equilibrium and kinetic properties. In this work an all-atom Go model is used to accurately represent the side-chain packing and native atomic contacts of the Trp-cage. The model reproduces the hallmark thermodynamics cooperativity of small proteins. Folding simulations observe that in the fast-folding dominant pathway, partial alpha-helical structure forms before hydrophobic core collapse. In the slow-folding secondary pathway, partial core collapse occurs before helical structure. The slow-folding rate of the secondary pathway is attributed to the loss of side-chain rotational freedom, due to the early core collapse, which impedes the helix formation. A major finding is the observation of a low-temperature kinetic intermediate stabilized by a salt bridge between residues Asp-9 and Arg-16. Similar observations [R. Zhou, Proc. Natl. Acad. Sci. U.S.A. 100, 13280 (2003)] were reported in a recent study using an all-atom model of the Trp-cage in explicit water, in which the salt-bridge stabilized intermediate was hypothesized to be the origin of the ultrafast-folding mechanism. A theoretical mutation that eliminates the Asp-9-Arg-16 salt bridge, but leaves the residues intact, is performed. Folding simulations of the mutant Trp-cage observe a two-state free-energy landscape with no kinetic intermediate and a significant decrease in the folding rate, in support of the hypothesis.     

On the interaction of bovine serum albumin with ionic surfactants: temperature induced EPR changes of a maleimide nitroxide reflect local protein dynamics and probe solvent accessibility

The interaction of bovine serum albumin (BSA) with the ionic surfactants sodium dodecylsulfate (SDS, anionic), cetyltrimethylammonium chloride (CTAC, cationic) and N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS, zwitterionic) was studied by electron paramagnetic resonance (EPR) spectroscopy of spin label covalently bound to the single free thiol group of the protein. EPR spectra simulation allows to monitor the protein dynamics at the labeling site and to estimate the changes in standard Gibbs free energy, enthalpy and entropy for transferring the nitroxide side chain from the more motionally restricted to the less restricted component. Whereas SDS and CTAC showed similar increases in the dynamics of the protein backbone for all measured concentrations, HPS presented a smaller effect at concentrations above 1.5mM. At 10mM of surfactants and 0.15 mM BSA, the standard Gibbs free energy change was consistent with protein backbone conformations more expanded and exposed to the solvent as compared to the native protein, but with a less pronounced effect for HPS. In the presence of the surfactants, the enthalpy change, related to the energy required to dissociate the nitroxide side chain from the protein, was greater, suggesting a lower water activity. The nitroxide side chain also detected a higher viscosity environment in the vicinity of the paramagnetic probe induced by the addition of the surfactants. The results suggest that the surfactant-BSA interaction, at higher surfactant concentration, is affected by the affinities of the surfactant to its own micelles and micelle-like aggregates. Complementary DLS data suggests that the temperature induced changes monitored by the nitroxide probe reflects local changes in the vicinity of the single thiol group of Cys-34 BSA residue.     

Mapping the hydration dynamics of ubiquitin

The nature of water's interaction with biomolecules such as proteins has been difficult to examine in detail at atomic resolution. Solution NMR spectroscopy is potentially a powerful method for characterizing both the structural and temporal aspects of protein hydration but has been plagued by artifacts. Encapsulation of the protein of interest within the aqueous core of a reverse micelle particle results in a general slowing of water dynamics, significant reduction in hydrogen exchange chemistry and elimination of contributions from bulk water thereby enabling the use of nuclear Overhauser effects to quantify interactions between the protein surface and hydration water. Here we extend this approach to allow use of dipolar interactions between hydration water and hydrogens bonded to protein carbon atoms. By manipulating the molecular reorientation time of the reverse micelle particle through use of low viscosity liquid propane, the T(1ρ) relaxation time constants of (1)H bonded to (13)C were sufficiently lengthened to allow high quality rotating frame nuclear Overhauser effects to be obtained. These data supplement previous results obtained from dipolar interactions between the protein and hydrogens bonded to nitrogen and in aggregate cover the majority of the molecular surface of the protein. A wide range of hydration dynamics is observed. Clustering of hydration dynamics on the molecular surface is also seen. Regions of long-lived hydration water correspond with regions of the protein that participate in molecular recognition of binding partners suggesting that the contribution of the solvent entropy to the entropy of binding has been maximized through evolution.     

Protein stability at a carbon nanotube interface

The interactions of proteins with solid surfaces occur in a variety of situations. Motivated by the many nanoengineering applications of protein-carbon nanotube hybrids, we investigate the conformational transitions of hen egg white lysozyme adsorbed on a carbon nanotube. Using a C(α) structure-based model and replica exchange molecular dynamics, we show how the folding/unfolding equilibrium of the adsorbed protein varies with the strength of its coupling to the surface. The stability of the native state depends on the balance between the favorable entropy and unfavorable enthalpy change on adsorption. In the case of a weakly attractive surface when the former dominates, the protein is stabilized. In this regime, the protein can fold and unfold while maintaining the same binding fraction. With increasing surface attraction, the unfavorable enthalpic effect dominates, the native state is destabilized, and the protein has to extensively unbind before changing states from unfolded to folded. At the highest surface coupling, the entropic penalty of folding vanishes, and a folding intermediate is strongly stabilized. In this intermediate state, the α-domain of lysozyme is disrupted, while the β-sheet remains fully structured. We rationalize the relative stability of the two domains on the basis of the residue contact order.     

Cooperative alpha-helix formation of beta-lactoglobulin induced by sodium n-alkyl sulfates

It is generally assumed that folding intermediates contain partially formed native-like secondary structures. However, if we consider the fact that the conformational stability of the intermediate state is simpler than that of the native state, it would be expected that the secondary structures in a folding intermediate would not necessarily be similar to those of the native state. beta-Lactoglobulin is a predominantly beta-sheet protein, although it has a markedly high intrinsic preference for alpha-helical structure. The formation of non-native alpha-helical intermediate of beta-lactoglobulin was induced by n-alkyl sulfates including sodium octyl sulfate, SOS; sodium decyl sulfate, SDeS; sodium dodecyl sulfate, SDS; and sodium tetradecyl sulfate, STS at special condition. The effect of n-alkyl sulfates on the structure of native beta-lactoglobulin at pH 2 was utilized to investigate the contribution of hydrophobic interactions to the stability of non-native alpha-helical intermediate. The addition of various concentrations of n-alkyl sulfates to the native state of beta-lactoglobulin (pH 2) appears to support the stabilized form of non-native alpha-helical intermediate at pH 2. The m values of the intermediate state of beta-lactoglobulin by SOS, SDeS, SDS and STS showed substantial variation. The enhancement of m values as the stability criterion of non-native alpha-helical intermediate state corresponded with increasing chain length of the cited n-alkyl sulfates. The present results suggest that the folding reaction of beta-lactoglobulin follows a non-hierarchical mechanism and hydrophobic interactions play important roles in stabilizing the non-native alpha-helical intermediate state.     

Single‐Molecule Detection Reveals Knot Sliding in TrmD Denaturation

An increasing number of proteins are found to contain a knot in their polypeptide chain. Although some studies have looked into the folding mechanism of knotted proteins, why and how these complex topologies form are still far from being fully answered. Moreover, no experimental information about how the knot moves during the protein‐folding process is available. Herein, by combining single‐molecule fluorescence resonance energy transfer (smFRET) experiments with molecular dynamics (MD) simulations, we performed a detailed study to characterize the knot in the denatured state of TrmD, a knotted tRNA (guanosine‐1) methyltransferase from Escherichia coli, as a model system. We found that the knot still existed in the unfolded state of TrmD, consistent with the results for two other knotted proteins, YibK and YbeA. More interestingly, both smFRET experiments and MD simulations revealed that the knot slid towards the C‐terminal during the unfolding process, which could be explained by the relatively strong interactions between the β‐sheet core at the N terminal of the native knot region. The size of the knot in the unfolded state is not larger than that in the native state. In addition, the knot slid in a “downhill” mode with simultaneous chain collapse in the denatured state.     

A Bioresistant Nitroxide Spin Label for In‐Cell EPR Spectroscopy: In Vitro and In Oocytes Protein Structural Dynamics Studies

Approaching protein structural dynamics and protein–protein interactions in the cellular environment is a fundamental challenge. Owing to its absolute sensitivity and to its selectivity to paramagnetic species, site‐directed spin labeling (SDSL) combined with electron paramagnetic resonance (EPR) has the potential to evolve into an efficient method to follow conformational changes in proteins directly inside cells. Until now, the use of nitroxide‐based spin labels for in‐cell studies has represented a major hurdle because of their short persistence in the cellular context. The design and synthesis of the first maleimido‐proxyl‐based spin label (M‐TETPO) resistant towards reduction and being efficient to probe protein dynamics by continuous wave and pulsed EPR is presented. In particular, the extended lifetime of M‐TETPO enabled the study of structural features of a chaperone in the absence and presence of its binding partner at endogenous concentration directly inside cells.