Chicken Villin Headpiece Research Articles

Curvature induced structural changes of the chicken villin headpiece subdomain by single walled carbon nanotubes.

Carbon nanotubes (CNTs) are identified as potential candidates for drug and biomolecular loading and delivery. CNTs of different chiralities have different diameters, which may significantly affect their abilities to interact with different types of biomolecules. Herein, we employ classical molecular dynamics simulation to provide insight into the curvature-dependent interactions between a model protein, chicken villin headpiece subdomain (HP36), with CNTs having chiralities (8,8), (12,12), and (20,20). It is revealed that, with increasing radii, the protein encounters more aromatic carbon atoms on the surface of the CNT, leading to its increasing strength of adsorption. However, the extent of adsorption has a limiting magnitude, after which an increase in the radius of the nanotube has practically no effect on the extent of adsorption. Spontaneous encapsulation of the protein was demonstrated using a (28,28) CNT, where the protein is found to undergo insignificant structural perturbation. Finally, steered molecular dynamics simulations have been performed to mimic the force-induced release of the protein from within the nanotube cavity. It has been identified that a minimum force of ∼300 pN and a minimum velocity of 4 Å ns-1 are required to release the protein from the CNT at 300 K. Any external force below the critical magnitude and inducing velocity less than 4 Å ns-1 allows the translocation of the protein through the inner surface of the CNT; however, before being released, the protein undergoes unfolding, thereby losing the secondary structure and biological activity.

Exploring the role of hydrophilic amino acids in unfolding of protein in aqueous ethanol solution.

Hydrophobic association is the key contributor behind the formation of well packed core of a protein which is often believed to be an important step for folding from an unfolded chain to its compact functional form. While most of the protein folding/unfolding studies have evaluated the changes in the hydrophobic interactions during chemical denaturation, the role of hydrophilic amino acids in such processes are not discussed in detail. Here we report the role of the hydrophilic amino acids behind ethanol induced unfolding of protein. Using free energy simulations, we show that chicken villin head piece (HP-36) protein unfolds gradually in presence of water-ethanol binary mixture with increasing composition of ethanol. However, upon mutation of hydrophilic amino acids by glycine while keeping the hydrophobic amino acids intact, the compact state of the protein is found to be stable at all compositions with gradual flattening of the free energy landscape upon increasing compositions. The local environment around the protein in terms of ethanol/water number significantly differs in wild type protein compared to the mutated protein. The calculated Wyman-Tanford preferential binding coefficient of ethanol for wild type protein reveals that a greater number of cosolutes (here ethanol) bind to the unfolded state compared to its folded state. However, no significant increase in binding coefficient of ethanol at the unfolded state is found for mutated protein. Local-bulk partition coefficient calculation also suggests similar scenarios. Our results reveal that the weakening of hydrophobic interactions in aqueous ethanol solution along with larger preferential binding of ethanol to the unfolded state mediated by hydrophilic amino acids combinedly helps unfolding of protein in aqueous ethanol solution.

Defect-assisted protein HP35 denaturation on graphene.

Structural defects in nanomaterials can alter their physical and chemical properties including magnetization, electronic and thermal conductivities, light absorption, and emission capabilities. Here, we investigated the potential impact of these defects on their biological effects through molecular dynamics simulations. By modeling the interaction between a graphene nanosheet and a widely used model protein, the chicken villin headpiece subdomain (HP35), we observed severe protein denaturation upon contact with defective graphene, while the protein remained intact on ideal graphene. The enhanced toxicity of defective graphene was due to the stronger attraction of the surface residues of HP35 from the defect edges (represented by carboxyl groups in our simulations) than from the ideal graphene. Upon binding to defective graphene, the contacting residues were restrained near the defective sites, which acted as "anchors" for the adsorbed protein. The "anchors" subsequently caused the protein to expose its aromatic and hydrophobic core residues to the graphene surface, via strong π-π stacking and hydrophobic interactions, thus leading to the unfolding of the protein. These findings not only highlight the importance of defects in nanomaterials' impact on biological systems, but also provide insights into fine-tuning the potential biological properties of nanomaterials through defect engineering.

The aqueous environment as an active participant in the protein folding process

Existing computational models applied in the protein structure prediction process do not sufficiently account for the presence of the aqueous solvent. The solvent is usually represented by a predetermined number of H2O molecules in the bounding box which contains the target chain. The fuzzy oil drop (FOD) model, presented in this paper, follows an alternative approach, with the solvent assuming the form of a continuous external hydrophobic force field, with a Gaussian distribution. The effect of this force field is to guide hydrophobic residues towards the center of the protein body, while promoting exposure of hydrophilic residues on its surface. This work focuses on the following sample proteins: Engrailed homeodomain (RCSB: 1enh), Chicken villin subdomain hp-35, n68h (RCSB: 1yrf), Chicken villin subdomain hp-35, k65(nle), n68h, k70(nle) (RCSB: 2f4k), Thermostable subdomain from chicken villin headpiece (RCSB: 1vii), de novo designed single chain three-helix bundle (a3d) (RCSB: 2a3d), albumin-binding domain (RCSB: 1prb) and lambda repressor-operator complex (RCSB: 1lmb).

A network of conformational transitions in an unfolding process of HP-35 revealed by high-temperature MD simulation and a Markov state model**Project supported by the National Natural Science Foundation of China (Grant Nos. 11175068 and 11474117) and the Self-determined Research Funds of CCNU from the Colleges Basic Research and Operation of MOE, China (Grant No. 230-20205170054).

An understanding of protein folding/unfolding processes has important implications for all biological processes, including protein degradation, protein translocation, aging, and diseases. All-atom molecular dynamics (MD) simulations are uniquely suitable for it because of their atomic level resolution and accuracy. However, limited by computational capabilities, nowadays even for small and fast-folding proteins, all-atom MD simulations of protein folding still presents a great challenge. An alternative way is to study unfolding process using MD simulations at high temperature. High temperature provides more energy to overcome energetic barriers to unfolding, and information obtained from studying unfolding can shed light on the mechanism of folding. In the present study, a 1000-ns MD simulation at high temperature (500 K) was performed to investigate the unfolding process of a small protein, chicken villin headpiece (HP-35). To infer the folding mechanism, a Markov state model was also built from our simulation, which maps out six macrostates during the folding/unfolding process as well as critical transitions between them, revealing the folding mechanism unambiguously.

Destabilization of Hydrophobic Core of Chicken Villin Headpiece in Guanidinium Chloride Induced Denaturation: Hint of π-Cation Interaction.

Despite their routine use as protein denaturants, the comprehensive understanding of the molecular mechanisms by which urea and guanidinium chloride (GdmCl) disrupts proteins' structure is still lacking. Here, we use steered molecular dynamics simulations along with the umbrella sampling technique to elucidate the mechanism of unfolding of chicken villin headpiece (HP-36) in these two denaturants. We find that while urea denatures protein predominantly by forming hydrogen bonds with the protein backbone, GdmCl commences unfolding by weakening of the hydrophobic interactions present in the core. The potential of mean force calculation indicates the reduction of hydrophobic interactions between two benzene moieties in 6 M GdmCl as compared to 6 M urea. We observe a near parallel orientation between the guanidinium cation and aromatic side chains of the HP-36 suggesting π-cation type stacking interactions which play a crucial role in weakening of the hydrophobic interaction. We use QM/MM optimization calculations to estimate the energetics of this π-cation interaction. Additionally, the consistency of the unfolding paths between high temperature (400 K) unfolding simulations and steered molecular dynamics simulations strengthens the proposed molecular mechanism of unfolding further.

Protein dynamics in the solid state from 2H NMR line shape analysis: a consistent perspective.

Deuterium line shape analysis of CD3 groups has emerged as a particularly useful tool for studying microsecond-millisecond protein motions in the solid state. The models devised so far consist of several independently conceived simple jump-type motions. They are comprised of physical quantities encoded in their simplest form; improvements are only possible by adding yet another simple motion, thereby changing the model. The various treatments developed are case-specific; hence comparison among the different systems is not possible. Here we develop a new methodology for (2)H NMR line shape analysis free of these limitations. It is based on the microscopic-order-macroscopic-disorder (MOMD) approach. In MOMD motions are described by diffusion tensors, spatial restrictions by potentials/ordering tensors, and geometric features by relative tensor orientations. Jump-type motions are recovered in the limit of large orientational potentials. Model improvement is accomplished by monitoring the magnitude, symmetry, and orientation of the various tensors. The generality of MOMD makes possible comparison among different scenarios. CD3 line shapes from the Chicken Villin Headpiece Subdomain and the Streptomyces Subtilisin Inhibitor are used as experimental examples. All of these spectra are reproduced by using rhombic local potentials constrained for simplicity to be given by the L = 2 spherical harmonics, and by axial diffusion tensors. Potential strength and rhombicity are found to be ca. 2-3 k(B)T. The diffusion tensor is tilted at 120° from the C-CD3 axis. The perpendicular (parallel) correlation times for local motion are 0.1-1.0 ms (3.3-30 μs). Activation energies in the 1.1-8.0 kcal/mol range are estimated. Future prospects include extension to the (2)H relaxation limit, application to the (15)N and (13)C NMR nuclei, and accounting for collective motions and anisotropic media.

Sensitivity of polarization fluctuations to the nature of protein-water interactions: study of biological water in four different protein-water systems.

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (∼80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of ∼2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water ⟨δMW(0)δMW(t)⟩ is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (∼50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work.

Multidimensional free energy surface of unfolding of HP-36: microscopic origin of ruggedness.

The protein folding funnel paradigm suggests that folding and unfolding proceed as directed diffusion in a multidimensional free energy surface where a multitude of pathways can be traversed during the protein's sojourn from initial to final state. However, finding even a single pathway, with the detail chronicling of intermediates, is an arduous task. In this work we explore the free energy surface of unfolding pathway through umbrella sampling, for a small globular α-helical protein chicken-villin headpiece (HP-36) when the melting of secondary structures is induced by adding DMSO in aqueous solution. We find that the unfolding proceeds through the initial separation or melting of aggregated hydrophobic core that comprises of three phenylalanine residues (Phe7, Phe11, and Phe18). This separation is accompanied by simultaneous melting of the second helix. Unfolding is found to be a multistage process involving crossing of three consecutive minima and two barriers at the initial stage. At a molecular level, Phe18 is observed to reorient itself towards other hydrophobic grooves to stabilize the intermediate states. We identify the configuration of the intermediates and correlate the intermediates with those obtained in our previous works. We also give an estimate of the barriers for different transition states and observe the softening of the barriers with increasing DMSO concentration. We show that higher concentration of DMSO tunes the unfolding pathway by destabilizing the third minimum and stabilizing the second one, indicating the development of a solvent modified, less rugged pathway. The prime outcome of this work is the demonstration that mixed solvents can profoundly transform the nature of the energy landscape and induce unfolding via a modified route. A successful application of Kramer's rate equation correlating the free energy simulation results shows faster rate of unfolding with increasing DMSO concentration. This work perhaps presents the first systematic theoretical study of the effect of a chemical denaturant on the microscopic free energy surface and rates of unfolding of HP-36.

Salt effects on hydrophobic-core formation in folding of a helical miniprotein studied by molecular dynamics simulations

We have investigated effects of salt ions on folding events of a helical miniprotein chicken villin headpiece subdomain HP36. Low concentrations of ions alter electrostatic interactions between charged groups of a protein and can change the populations of conformers. Here, we compare two data sets of folding simulations of HP36 in explicit water solvent with or without ions. For efficient sampling of the conformational space of HP36, the multicanonical replica-exchange molecular dynamics method was employed. Our analyses suggest that salt alters salt-bridging nature of the protein at later stages of folding at room temperature. Especially, more nonnative, nonlocal salt bridges are formed at near-native conformations in pure water. Our analyses also show that such salt-bridge formation hinders the fully native hydrophobic-core packing at the final stages of folding.

Solvent Sensitivity of Protein Unfolding: Dynamical Study of Chicken Villin Headpiece Subdomain in Water–Ethanol Binary Mixture

We carry out a series of long atomistic molecular dynamics simulations to study the unfolding of a small protein, chicken villin headpiece (HP-36), in water-ethanol (EtOH) binary mixture. The prime objective of this work is to explore the sensitivity of protein unfolding dynamics toward increasing concentration of the cosolvent and unravel essential features of intermediates formed in search of a dynamical pathway toward unfolding. In water-ethanol binary mixtures, HP-36 is found to unfold partially, under ambient conditions, that otherwise requires temperature as high as ∼600 K to denature in pure aqueous solvent. However, an interesting course of pathway is observed to be followed in the process, guided by the formation of unique intermediates. The first step of unfolding is essentially the separation of the cluster formed by three hydrophobic (phenylalanine) residues, namely, Phe-7, Phe-11, and Phe-18, which constitute the hydrophobic core, thereby initiating melting of helix-2 of the protein. The initial steps are similar to temperature-induced unfolding as well as chemical unfolding using DMSO as cosolvent. Subsequent unfolding steps follow a unique path. As water-ethanol shows composition-dependent anomalies, so do the details of unfolding dynamics. With an increase in cosolvent concentration, different partially unfolded intermediates are found to be formed. This is reflected in a remarkable nonmonotonic composition dependence of several order parameters, including fraction of native contacts and protein-solvent interaction energy. The emergence of such partially unfolded states can be attributed to the preferential solvation of the hydrophobic residues by the ethyl groups of ethanol. We further quantify the local dynamics of unfolding by using a Marcus-type theory.

Soliton driven relaxation dynamics and protein collapse in the villin headpiece

Protein collapse from a random chain to the native state involves a dynamical phase transition. During the process, new scales and collective variables become excited while old ones recede and fade away. The presence of different phases and many scales causes formidable computational bottle-necks in approaches that are based on full atomic scale scrutiny. Here we propose a way to describe the folding and unfolding processes effectively, using only the biologically relevant time and distance scales. We merge a coarse grained Landau theory that models the static collapsed protein in the low-temperature limit with a Glauber protocol that describes finite-temperature relaxation dynamics in a statistical system which is out of thermal equilibrium. As an example we inspect the collapse of a HP35 chicken villin headpiece subdomain, a paradigm specimen in protein folding studies. We simulate the folding and unfolding process by repeated heating and cooling cycles between a given low-temperature, i.e. bad solvent, environment where the protein is collapsed and various different high-temperature, i.e. good solvent, environments. We find that as long as the high temperature value stays below a value in the range that separates the random walk phase from the self-avoiding walk phase, we consistently recover the native state upon cooling. But, when heated to sufficiently high temperatures, the native state practically never recurs. Our result confirms Anfinsen’s thermodynamical hypothesis and estimates a temperature range for its validity, in the case of villin.

Glassy Dynamics of Protein Methyl Groups Revealed by Deuteron NMR

We investigated site-specific dynamics of key methyl groups in the hydrophobic core of chicken villin headpiece subdomain (HP36) over the temperature range between 298 and 140 K using deuteron solid-state NMR longitudinal relaxation measurements. The relaxation of the longitudinal magnetization is weakly nonexponential (glassy) at high temperatures and exhibits a stronger degree of nonexponentiality below about 175 K. In addition, the characteristic relaxation times deviate from the simple Arrhenius law. We interpret this behavior via the existence of distribution of activation energy barriers for the three-site methyl jumps, which originates from somewhat different methyl environments within the local energy landscape. The width of the distribution of the activation barriers for methyl jumps is rather significant, about 1.4 kJ/mol. Our experimental results and modeling allow for the description of the apparent change at about 175 K without invoking a specific transition temperature. For most residues in the core, the relaxation behavior at high temperatures points to the existence of conformational exchange between the substates of the landscape, and our model takes into account the kinetics of this process. The observed dynamics are the same for dry and hydrated protein. We also looked at the effect of F58L mutation inside the hydrophobic core on the dynamics of one of the residues and observed a significant increase in its conformational exchange rate constant at high temperatures.

Chemical Unfolding of Chicken Villin Headpiece in Aqueous Dimethyl Sulfoxide Solution: Cosolvent Concentration Dependence, Pathway, and Microscopic Mechanism

Unfolding of a protein often proceeds through partial unfolded intermediate states (PUIS). PUIS have been detected in several experimental and simulation studies. However, complete analyses of transitions between different PUIS and the unfolding trajectory are sparse. To understand such dynamical processes, we study chemical unfolding of a small protein, chicken villin head piece (HP-36), in aqueous dimethyl sulfoxide (DMSO) solution. We carry out molecular dynamics simulations at various solution compositions under ambient conditions. In each concentration, the initial step of unfolding involves separation of two adjacent native contacts, between phenyl alanine residues (11-18 and 7-18). This first step induces, under appropriate conditions, subsequent separation among other hydrophobic contacts, signifying a high degree of cooperativity in the unfolding process. The observed sequence of structural changes in HP-36 on increasing DMSO concentration and the observed sequence of PUIS, are in approximate agreement with earlier simulation results (in pure water) and experimental observations on unfolding of HP-36. Peculiar to water-DMSO mixture, an intervening structural transformation (around 15% of DMSO) in the binary mixture solvent retards the progression of unfolding as composition is increased. This is reflected in a remarkable nonmonotonic composition dependence of RMSD, radius of gyration and the fraction of native contacts. At 30% mole fraction of DMSO, we find the extended randomly coiled structure of the unfolded protein. The molecular mechanism of DMSO induced unfolding process is attributed to the initial preferential solvation of the hydrophobic side chain atoms through the methyl groups of DMSO, followed by the hydrogen bonding of the oxygen atom of DMSO to the exposed backbone NH groups of HP-36.

Correlation between protein secondary structure, backbone bond angles, and side-chain orientations

We investigate the fine structure of the sp3 hybridized covalent bond geometry that governs the tetrahedral architecture around the central C(α) carbon of a protein backbone, and for this we develop new visualization techniques to analyze high-resolution x-ray structures in the Protein Data Bank. We observe that there is a correlation between the deformations of the ideal tetrahedral symmetry and the local secondary structure of the protein. We propose a universal coarse-grained energy function to describe the ensuing side-chain geometry in terms of the C(β) carbon orientations. The energy function can model the side-chain geometry with a subatomic precision. As an example we construct the C(α)-C(β) structure of HP35 chicken villin headpiece. We obtain a configuration that deviates less than 0.4 Å in root-mean-square distance from the experimental x-ray structure.

Conformational Dynamics and Stability of HP35 Studied with 2D IR Vibrational Echoes

Two-dimensional infrared (2D IR) vibrational echo spectroscopy was used to measure the fast dynamics of two variants of chicken villin headpiece 35 (HP35). The CN of cyanophenylalanine residues inserted in the hydrophobic core were used as a vibrational probe. Experiments were performed on both singly (HP35-P) and doubly CN-labeled peptide (HP35-P(2)) within the wild-type sequence, as well as on HP-35 containing a singly labeled cyanophenylalanine and two norleucine mutations (HP35-P NleNle). There is a remarkable similarity between the dynamics measured in singly and doubly CN-labeled HP35, demonstrating that the presence of an additional CN vibrational probe does not significantly alter the dynamics of the small peptide. The substitution of two lysine residues by norleucines markedly improves the stability of HP35 by replacing charged with nonpolar residues, stabilizing the hydrophobic core. The results of the 2D IR experiments reveal that the dynamics of HP35-P are significantly faster than those of HP35-P NleNle. These observations suggest that the slower structural fluctuations in the hydrophobic core, indicating a more tightly structured core, may be an important contributing factor to HP35-P NleNle's increased stability.

Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these "slow" water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.

Slow motions in the hydrophobic core of chicken villin headpiece subdomain and their contributions to configurational entropy and heat capacity from solid-state deuteron NMR measurements.

We have investigated microsecond to millisecond time scale dynamics in several key hydrophobic core methyl groups of chicken villin headpiece subdomain protein (HP36) using a combination of single-site labeling, deuteron solid-state NMR line shape analysis, and computational modeling. Deuteron line shapes of hydrated powder samples are dominated by rotameric jumps and show a large variability of rate constants, activation energies, and rotameric populations. Site-specific activation energies vary from 6 to 38 kJ/mol. An additional mode of diffusion on a restricted arc is significant for some sites. In dry samples, the dynamics is quenched. Parameters of the motional models allow for calculations of configurational entropy and heat capacity, which, together with the rate constants, allow for observation of interplay between thermodynamic and kinetic picture of the landscape. Mutations at key phenylalanine residues at both distal (F47L&F51L) and proximal (F58L) locations to a relatively rigid side chain of L69 have a pronounced effect on alleviating the rigidity of this side chain at room temperature and demonstrate the sensitivity of the hydrophobic core environment to such perturbations.

Temperature dependence of fast carbonyl backbone dynamics in chicken villin headpiece subdomain

Temperature-dependence of protein dynamics can provide information on details of the free energy landscape by probing the characteristics of the potential responsible for the fluctuations. We have investigated the temperature-dependence of picosecond to nanosecond backbone dynamics at carbonyl carbon sites in chicken villin headpiece subdomain protein using a combination of three NMR relaxation rates: (13)C' longitudinal rate, and two cross-correlated rates involving dipolar and chemical shift anisotropy (CSA) relaxation mechanisms, (13)C'/(13)C'-(13)C(α) CSA/dipolar and (13)C'/(13)C'-(15)N CSA/dipolar. Order parameters have been extracted using the Lipari-Szabo model-free approach assuming a separation of the time scales of internal and molecular motions in the 2-16°C temperature range. There is a gradual deviation from this assumption from lower to higher temperatures, such that above 16°C the separation of the time scales is inconsistent with the experimental data and, thus, the Lipari-Szabo formalism can not be applied. While there are variations among the residues, on the average the order parameters indicate a markedly steeper temperature dependence at backbone carbonyl carbons compared to that probed at amide nitrogens in an earlier study. This strongly advocates for probing sites other than amide nitrogen for accurate characterization of the potential and other thermodynamics characteristics of protein backbone.

Comparative Dynamics of Leucine Methyl Groups in FMOC-Leucine and in a Protein Hydrophobic Core Probed by Solid-State Deuteron Nuclear Magnetic Resonance over 7−324 K Temperature Range

Quantitative dynamics of methyl groups in 9-fluorenylmethyloxycarbonyl-leucine (FMOC-leu) have been analyzed and compared with earlier studies of methyl dynamics in chicken villin headpiece subdomain protein (HP36) labeled at L69, a key hydrophobic core position. A combination of deuteron solid-state nuclear magnetic resonance experiments over the temperature range of 7-324 K and computational modeling indicated that while the two compounds show the same modes of motions, there are marked differences in the best-fit parameters of these motions. One of the main results is that the crossover observed in the dynamics of the methyl groups in the HP36 sample at 170 K is absent in FMOC-leu. A second crossover at around 95-88 K is present in both samples. The differences in the behavior of the two compounds suggest that some of the features of methyl dynamics reflect the complexity of the protein hydrophobic core and are not determined solely by local interactions.