The serine protease α-thrombin is a dual-action protein that mediates the blood-clotting cascade. Thrombin alone is a procoagulant, cleaving fibrinogen to make the fibrin clot, but the thrombinthrombomodulin (TM) complex initiates the anticoagulant pathway by cleaving protein C. A TM fragment consisting of only the fifth and sixth EGF-like domains (TM56) is sufficient to bind thrombin, but the presence of the fourth EGF-like domain (TM456) is critical to induce the anticoagulant activity of thrombin. Crystallography of the thrombin-TM456 complex revealed no significant structural changes in thrombin, suggesting that TM4 may only provide a scaffold for optimal alignment of protein C for its cleavage by thrombin. However, a variety of experimental data have suggested that the presence of TM4 may affect the dynamic properties of the active site loops. In the present work, we have used both conventional and accelerated molecular dynamics simulation to study the structural dynamic properties of thrombin, thrombin: TM56, and thrombin:TM456 across a broad range of time scales.Two distinct yet interrelated allosteric pathways are identified that mediate both the pro-and anticoagulant activities of thrombin. One allosteric pathway, which is present in both thrombin: TM56 and thrombin:TM456, directly links the TM5 domain to the thrombin active site. The other allosteric pathway, which is only present on slow time scales in the presence of the TM4 domain, involves an extended network of correlated motions linking the TM4 and TM5 domains and the active site loops of thrombin.allostery | thrombosis

The backbone dynamics of human α-thrombin inhibited at the active site serine were analyzed using R(1), R(2), and heteronuclear NOE experiments, variable temperature TROSY 2D [(1)H-(15)N] correlation spectra, and R(ex) measurements. The N-terminus of the heavy chain, which is formed upon zymogen activation and inserts into the protein core, is highly ordered, as is much of the double beta-barrel core. Some of the surface loops, by contrast, remain very dynamic with order parameters as low as 0.5 indicating significant motions on the ps-ns timescale. Regions of the protein that were thought to be dynamic in the zymogen and to become rigid upon activation, in particular the γ-loop, the 180s loop, and the Na(+) binding site have order parameters below 0.8. Significant R(ex) was observed in most of the γ-loop, in regions proximal to the light chain, and in the β-sheet core. Accelerated molecular dynamics simulations yielded a molecular ensemble consistent with measured residual dipolar couplings that revealed dynamic motions up to milliseconds. Several regions, including the light chain and two proximal loops, did not appear highly dynamic on the ps-ns timescale, but had significant motions on slower timescales.

It is well known that high hydrostatic pressures can induce the unfolding of proteins. The physical underpinnings of this phenomenon have been investigated extensively but remain controversial. Changes in solvation energetics have been commonly proposed as a driving force for pressure-induced unfolding. Recently, the elimination of void volumes in the native folded state has been argued to be the principal determinant. Here we use the cavity-containing L99A mutant of T 4 lysozyme to examine the pressure-induced destabilization of this multidomain protein by using solution NMR spectroscopy. The cavity-containing C-terminal domain completely unfolds at moderate pressures, whereas the N-terminal domain remains largely structured to pressures as high as 2.5 kbar. The sensitivity to pressure is suppressed by the binding of benzene to the hydrophobic cavity. These results contrast to the pseudo-WT protein, which has a residual cavity volume very similar to that of the L99A-benzene complex but shows extensive subglobal reorganizations with pressure. Encapsulation of the L99A mutant in the aqueous nanoscale core of a reverse micelle is used to examine the hydration of the hydrophobic cavity. The confined space effect of encapsulation suppresses the pressure-induced unfolding transition and allows observation of the filling of the cavity with water at elevated pressures. This indicates that hydration of the hydrophobic cavity is more energetically unfavorable than global unfolding. Overall, these observations point to a range of cooperativity and energetics within the T 4 lysozyme molecule and illuminate the fact that small changes in physical parameters can significantly alter the pressure sensitivity of proteins.protein stability | protein folding and cooperativity | protein hydration | high-pressure NMR | reverse micelle encapsulation

Cetyltrimethylammonium bromide (CTAB)/hexanol reverse micelles have found a variety of applications that demand control over physical parameters. Water content or loading is among the most basic tunable components and is the major driver of the physical properties of these systems. This study uses small-angle scattering with contrast variation to characterize these systems as a function of water loading. The scattering data were analyzed with a variety of approaches, resulting in converging physical specifications. Equations that describe basic physical parameters were determined that allow for characterization and manipulation of the CTAB/hexanol reverse micelle surfactant system. The shape of the reverse micelles was revealed to be slightly ellipsoidal and varies slightly through the water loading range. The surfactant shell is shown to contain a higher fraction of hexanol upon addition of water. Analysis reveals that the size, shape, and surfactant/cosurfactant composition are directly tunable by variation of the water content and that these properties are consequences of the balance of forces present in the reverse micelles.

Thrombin is the central protease in the cascade of blood coagulation proteases. The structure of thrombin consists of a double β-barrel core surrounded by connecting loops and helices. Compared to chymotrypsin, thrombin has more extended loops that are thought to have arisen from insertions in the serine protease that evolved to impart greater specificity. Previous experiments showed thermodynamic coupling between ligand binding at the active site and distal exosites. We present a combined approach of molecular dynamics (MD), accelerated molecular dynamics (AMD), and analysis of the residual local frustration of apo-thrombin and active-site-bound (PPACK-thrombin). Community analysis of the MD ensembles identified changes upon active site occupation in groups of residues linked through correlated motions and physical contacts. AMD simulations, calibrated on measured residual dipolar couplings, reveal that upon active site ligation, correlated loop motions are quenched, but new ones connecting the active site with distal sites where allosteric regulators bind emerge. Residual local frustration analysis reveals a striking correlation between frustrated contacts and regions undergoing slow time scale dynamics. The results elucidate a motional network that probably evolved through retention of frustrated contacts to provide facile conversion between ensembles of states.

Encapsulation of small molecules, proteins, and other macromolecules within the protective water core of reverse micelles is emerging as a powerful strategy for a variety of applications. The cationic surfactant cetyltrimethylammonium bromide (CTAB) in combination with hexanol as a cosurfactant is particularly useful in the context of solution NMR spectroscopy of encapsulated proteins. Small-angle X-ray and neutron scattering is employed to investigate the internal structure of the CTAB/hexanol reverse micelle particle under conditions appropriate for high-resolution NMR spectroscopy. The scattering profiles are used to benchmark extensive molecular dynamics simulations of this reverse micelle system and indicate that the parameters used in these simulations recapitulate experimental results. Scattering profiles and simulations indicate formation of homogeneous solutions of small approximately spherical reverse micelle particles at a water loading of 20 composed of ∼150 CTAB and 240 hexanol molecules. The 3000 waters comprising the reverse micelle core show a gradient of translational diffusion that reaches that of bulk water at the center. Rotational diffusion is slowed relative to bulk throughout the water core, with the greatest slowing near the CTAB headgroups. The 5 Å thick interfacial region of the micelle consists of overlapping layers of Br(-) enriched water, CTAB headgroups, and hexanol hydroxyl groups, containing about one-third of the total water. This study employs well-parametrized MD simulations, X-ray and neutron scattering, and electrostatic theory to illuminate fundamental properties of CTAB/hexanol reverse micelle size, shape, partitioning, and water behavior.

Reverse micelle (RM) encapsulation of proteins for NMR spectroscopy has many advantages over standard NMR methods such as enhanced tumbling and improved sensitivity. It has opened many otherwise difficult lines of investigation including the study of membrane associated proteins, large soluble proteins, unstable protein states, and the study of protein surface hydration dynamics. Recent technological developments have extended the ability of RM encapsulation with high structural fidelity for nearly all proteins and thereby allow high-quality state-of-the-art NMR spectroscopy. Optimal conditions are achieved using a streamlined screening protocol, which is described here. Commonly studied proteins spanning a range of molecular weights are used as examples. Very low-viscosity alkane solvents, such as propane or ethane, are useful for studying very large proteins but require the use of specialized equipment to permit preparation and maintenance of well-behaved solutions under elevated pressure. The procedures for the preparation and use of solutions of reverse micelles in liquefied ethane and propane are described. The focus of this chapter is to provide procedures to optimally encapsulate proteins in reverse micelles for modern NMR applications.

No abstract