Concerted motion of a protein-peptide complex: backbone dynamics studies of an (15)N-labeled peptide derived from P(21)-activated kinase bound to Cdc42Hs.GMPPCP.

Cdc42Hs is a member of the Ras superfamily of GTPases which, when active, initiates a cascade beginning with the activation of several kinases, including P(21)-activated kinase (PAK). We previously determined the structure of a complex between a 46 amino acid fragment peptide derived from the PAK binding domain (PBD46) and Cdc42Hs.GMPPCP (Gizachew, D., Guo, W., Chohan, K. K., Sutcliffe, M. J., and Oswald, R. E. (2000) Biochemistry 39, 3963-3971). Previous studies (Loh, A. P., Guo, W., Nicholson, L. K., and Oswald, R. E. (1999) Biochemistry 38, 12547-12557) suggest that the regions of Cdc42Hs that bind effectors and regulators have distinct dynamic properties from the remainder of the protein. Here, we describe the backbone dynamics of PBD46 bound to Cdc42Hs.GMPPCP. T(1), T(2), T(1)(rho), and steady-state nuclear Overhauser effects were measured at 500 and 600 MHz. An extension of the Lipari-Szabo model-free analysis was used to determine the order parameters (S(2)) and local correlation times (tau(e)) of the N-H bond vectors within PBD46. Both Cdc42Hs and PBD46 exhibit increased mobility in the free versus the bound state, suggesting that protein flexibility may be required for high-affinity PBD46 binding and, presumably, the activation of PAK. Different backbone dynamics were observed in different regions of the peptide. The beta-strand region of bound PBD46, which makes contacts with beta2 of Cdc42Hs, exhibits low mobility on the pico- to nanosecond timescale. However, the part of PBD46 that interacts with Switch I of Cdc42Hs exhibits greater mobility. Thus, PBD46 and Cdc42Hs form a tight complex that exhibits concerted dynamics.     

Dynamic and Thermodynamic Response of the Ras Protein Cdc42Hs upon Association with the Effector Domain of PAK3

Human cell division cycle protein 42 (Cdc42Hs) is a small, Rho-type guanosine triphosphatase involved in multiple cellular processes through its interactions with downstream effectors. The binding domain of one such effector, the actin cytoskeleton-regulating p21-activated kinase 3, is known as PBD46. Nitrogen-15 backbone and carbon-13 methyl NMR relaxation was measured to investigate the dynamical changes in activated GMPPCP·Cdc42Hs upon PBD46 binding. Changes in internal motion of the Cdc42Hs, as revealed by methyl axis order parameters, were observed not only near the Cdc42Hs–PBD46 interface but also in remote sites on the Cdc42Hs molecule. The binding-induced changes in side-chain dynamics propagate along the long axis of Cdc42Hs away from the site of PBD46 binding with sharp distance dependence. Overall, the binding of the PBD46 effector domain on the dynamics of methyl-bearing side chains of Cdc42Hs results in a modest rigidification, which is estimated to correspond to an unfavorable change in conformational entropy of approximately − 10 kcal mol^− 1 at 298 K. A cluster of methyl probes closest to the nucleotide-binding pocket of Cdc42Hs becomes more rigid upon binding of PBD46 and is proposed to slow the catalytic hydrolysis of the γ phosphate moiety. An additional cluster of methyl probes surrounding the guanine ring becomes more flexible on binding of PBD46, presumably facilitating nucleotide exchange mediated by a guanosine exchange factor. In addition, the Rho insert helix, which is located at a site remote from the PBD46 binding interface, shows a significant dynamic response to PBD46 binding.     

Main chain and side chain dynamics of the ubiquitin conjugating enzyme variant human Mms2 in the free and ubiquitin-bound States

Protein ubiquitination involves a cascade of enzymatic steps where ubiquitin (Ub) is sequentially transferred as a thiolester intermediate from an E1 enzyme to an E2 enzyme and finally to the protein target with the help of a Ub-protein ligase. Protein ubiquitination brought about by the Ubc13-Mms2 (E2-E2) complex has a unique role in the cell, unrelated to protein degradation. The Mms2-Ubc13 heterodimer links Ub molecules to one another through an isopeptide bond between its own C-terminus and Lys-63 on another Ub. The role of Mms2 is to orient a target-bound Ub molecule such that its Lys-63 is proximal to the C-terminus of the Ub molecule that is covalently linked to the active site of Ubc13. To gain insight into the influence of protein dynamics on the affinity of Ub for Mms2, we have determined pico- to nanosecond time scale fluctuations of the main chain and methyl side chains of human Mms2 in the free and Ub-bound states using solution state (15)N and (2)H nuclear magnetic resonance relaxation measurements. Analysis of the relaxation data allows for a semiquantitative estimation of the conformational entropy change (TDeltaS(NMR)) for the main chain and side chain methyl groups of Mms2 upon binding Ub. The value of TDeltaS(NMR) for the main chain and side chain methyl groups of Mms2 is -8 +/- 2 and -2 +/- 2 kcal mol(-)(1), respectively. The experimental DeltaG(binding) for the Mms2.Ub complex is -6 kcal mol(-)(1). Estimation of DeltaG(binding) using an empirical structure-based approach that does not account for changes in main chain entropy yields a value of -17 +/- 2 kcal mol(-)(1). However, inclusion of TDeltaS(NMR) for the main chain of Mms2 increases the estimated DeltaG(binding) to -9 +/- 3 kcal mol(-)(1). Assuming that changes in Ub main chain dynamics contribute to TDeltaS(NMR) to the same extent as Mms2, the estimated DeltaG(binding) is further reduced to -1 +/- 4 kcal mol(-)(1), a value close to the experimental DeltaG(binding).     

Main chain and side chain dynamics of oxidized flavodoxin from Cyanobacterium anabaena.

Oxidized flavodoxin from Cyanobacterium anabaena PCC 7119 is used as a model system to investigate the fast internal dynamics of a flavin-bearing protein. Virtually complete backbone and side chain resonance NMR assignments of an oxidized flavodoxin point mutant (C55A) have been determined. Backbone and side chain dynamics in flavodoxin (C55A) were investigated using (15)N amide and deuterium methyl NMR relaxation methods. The squared generalized order parameters (S(NH)(2)) for backbone amide N-H bonds are found to be uniformly high ( approximately 0.923 over 109 residues in regular secondary structure), indicating considerable restriction of motion in the backbone of the protein. In contrast, methyl-bearing side chains are considerably heterogeneous in their amplitude of motion, as indicated by obtained symmetry axis squared generalized order parameters (S(axis)(2)). However, in comparison to nonprosthetic group-bearing proteins studied with these NMR relaxation methods, the side chains of oxidized flavodoxin are unusually rigid.     

Backbone dynamics of inactive, active, and effector-bound Cdc42Hs from measurements of (15)N relaxation parameters at multiple field strengths.

Cdc42Hs, a member of the Ras superfamily of GTP-binding proteins, initiates a cascade that begins with the activation of several kinases, including p21-activated kinase (PAK). We have previously determined the structure of Cdc42Hs and found that the regions involved in effector (Switch I) and regulator (Switch II) actions are partially disordered [Feltham, J. L., et al. (1997) Biochemistry 36, 8755-8766]. Recently, we used a 46-amino acid fragment of PAK (PBD46) to define the binding surface on Cdc42Hs, which includes the beta2 strand and a portion of Switch I [Guo, W., et al. (1998) Biochemistry 37, 14030-14037]. Here we describe the backbone dynamics of three constructs of [(15)N]Cdc42Hs (GDP-, GMPPCP-, and GMPPCP- and PBD46-bound) using (15)N-(1)H NMR measurements of T(1), T(1)(rho), and the steady-state NOE at three magnetic field strengths. Residue-specific values of the generalized order parameters (S(s)(2) and S(f)(2)), local correlation time (tau(e)), and exchange rate (R(ex)) were obtained using the Lipari-Szabo model-free formalism. Residues in Switch I were found to exhibit high-amplitude (low-order) motions on a nanosecond time scale, whereas those in Switch II experience low-amplitude motion on the nanosecond time scale and chemical (conformational) exchange on a millisecond time scale. The Insert region of Cdc42Hs-GDP exhibits high-order, nanosecond motions; the time scale of motion in the Insert is reduced in Cdc42Hs-GMPPCP and Cdc42Hs-PBD46. Overall, significant flexibility was observed mainly in the regions of Cdc42Hs that are involved in protein-protein interactions (Switch I, Switch II, and Insert), and flexibility was reduced upon interaction with a protein ligand. These results suggest that protein flexibility is important for high-affinity binding interactions.     

Structure of the complex of Cdc42Hs with a peptide derived from P-21 activated kinase

Cdc42Hs is a member of the Ras superfamily of GTPases and initiates a cascade that begins with the activation of several kinases, including p21-activated kinase (PAK). We have previously used a 46 amino acid fragment of PAK (PBD46) to define the binding surface on Cdc42Hs [Guo et al. (1998) Biochemistry 37, 14030-14037]. Here we describe the three-dimensional solution structure of the Cdc42Hs. GMPPCP-PBD46 complex. Heteronuclear NMR methods were used to assign resonances in the complex, and approximately 2400 distance and dihedral restraints were used to calculate a set of 20 structures using a combination of distance geometry, simulated annealing, and chemical shift and Ramachandran refinement. The overall structure of Cdc42Hs in the complex differs from the uncomplexed structure in two major aspects: (1) the first alpha helix is reoriented to accommodate the binding of the peptide and (2) the regions corresponding to switch I and switch II are less disordered. As suggested by our previous work (Guo et al., 1998) and similar to the complex between Cdc42Hs and fACK [Mott et al. (1999) Nature 399, 384-388], PBD46 forms an intermolecular beta-sheet with beta2 of Cdc42Hs and contacts both switch I and switch II. The extensive binding surface between PBD46 and Cdc42Hs can account for both the high affinity of the complex and the inhibition by PBD46 of GTP hydrolysis.     

Characterization of the backbone and side chain dynamics of the CaM-CaMKIp complex reveals microscopic contributions to protein conformational entropy

Calmodulin is a central mediator of calcium-dependent signal transduction pathways and regulates the activity of a large number of diverse targets. Calcium-dependent interactions of calmodulin with regulated proteins are of generally high affinity but of quite variable thermodynamic origins. Here we investigate the influence of the binding of the calmodulin-binding domain of calmodulin kinase I on the fast internal dynamics of calcium-saturated calmodulin. NMR relaxation was used to probe motion on the backbone (viewed through the backbone amide NH group) and the side chains (viewed through methyl groups). The distribution of the amplitudes of side chain dynamics is trimodal. The microscopic details of side chain motion are compared with those of a thermodynamically and structurally similar complex of calmodulin with the calmodulin-binding domain of the smooth muscle myosin light chain kinase. While there are no significant differences in backbone dynamics and no net change in methyl-bearing side chain dynamics, a large redistribution of the amplitude of methyl dynamics is observed between the two complexes. The variation in dynamics was largely localized to the heterogeneously dynamic target-binding interface, suggesting that differential dynamics of the binding surface plays a functional role in the high-affinity binding interactions of calmodulin. These results begin to reveal a fundamental role for residual protein entropy in molecular recognition by calmodulin.     

Molecular dynamics of anthraquinone DNA intercalators with polyethylene glycol side chains

The interactions of a homologous series of four anthraquinone (AQ) intercalators with increasing lengths of polyethylene glycol (PEG) side chains with DNA have been studied via molecular dynamics (MD) simulations. The geometry, conformation, interactions, and hydration of the complexes were examined. The geometries of the four ligands were similar with parallel stacking to the long axis of the adjacent DNA base pairs. Hydrogen bonding between the AQ amide and DNA led to a preference for the trans-syn conformer. As the side chain lengthened, binding to DNA reduced the conformational space, resulting in an increase in unfavorable entropy. Increased localization of the PEG side chain in the DNA groove, indicating some interaction of the side chain with DNA, also contributed unfavorably to the entropy. The changes in free energy of binding due to entropic considerations (-3.9 to -6.3 kcal/mol) of AQ I-IV were significant. The hydration of the PEG side chain decreased upon binding to DNA. Understanding of side chain conformations, interactions, and hydration changes that accompany the formation of a ligand-DNA complex may be important in the development of new applications of pegylated small molecules that target biological macromolecules.     

Analysis of side chain mobility among protein G B1 domain mutants with widely varying stabilities

"Host-guest" studies of the B1 domain from Streptococcal protein G have been used previously to establish a thermodynamic scale for the beta-sheet-forming propensities of the 20 common amino acids. To investigate the contribution of side chain conformational entropy to the relative stabilities of B1 domain mutants, we have determined the dynamics of side chain methyl groups in 10 of the 20 mutants used in a previous study. Deuterium relaxation rates were measured using two-dimensional NMR techniques for 13CH2D groups. Analysis of the relaxation data using the Lipari-Szabo model-free formalism showed that mutations introduced at the guest position caused small but statistically significant changes in the methyl group dynamics. In addition, there was a low level of covariation of the Lipari-Szabo order parameters among the 10 mutants. The variations in conformational free energy estimated from the order parameters were comparable in magnitude to the variations in global stability of the 10 mutants but did not correlate with the global stability of the domain or with the structural properties of the guest amino acids. The data support the view that conformational entropy in the folded state is one of many factors that can influence the folding thermodynamics of proteins.     

Redistribution and loss of side chain entropy upon formation of a calmodulin-peptide complex

The response of the internal dynamics of calcium-saturated calmodulin to the formation of a complex with a peptide model of the calmodulin-binding domain of the smooth muscle myosin light chain kinase has been studied using NMR relaxation methods. The backbone of calmodulin is found to be unaffected by the binding of the domain, whereas the dynamics of side chains are significantly perturbed. The changes in dynamics are interpreted in terms of a heterogeneous partitioning between structure (enthalpy) and dynamics (entropy). These data provide a microscopic view of the residual entropy of a protein in two functional states and suggest extensive enthalpy/entropy exchange during the formation of a protein-protein interface.     

Model selection for the interpretation of protein side chain methyl dynamics

A number of different dynamics models are considered for fitting ¹³C and ²H side chain methyl relaxation rates. It is shown that in cases where nanosecond time scale dynamics are present the extended Lipari–Szabo model which is explicitly parameterized to include the effects of slow motions can produce wide distributions of fitting parameters even in cases where the errors are relatively small and large numbers of relaxation rates are considered. In contrast, fits of ¹⁵N backbone dynamics using this model are far more robust. The origin of this difference is analyzed and can be explained by the different functional forms of the spectral density in these two cases. The utility of a number of models for the analysis of methyl side chain dynamics is presented.     

Probing protein side chain dynamics via 13C NMR relaxation

Protein side chain dynamics is associated with protein stability, folding, and intermolecular interactions. Detailed dynamics information is crucial for the understanding of protein function and biochemical and biophysical properties, which can be obtained using NMR relaxation techniques. In this review, (13)C relaxation of methine, methylene and methyl groups with and without (1)H decoupling are described briefly for a better understanding of how spin relaxation is associated with motional (dynamics) parameters. Developments in the measurement and interpretation of (13)C auto-relaxation and cross-correlated relaxation data are presented too. Finally, recent progress in the use of (13)C relaxation to probe the dynamics of protein side chains is detailed mainly for the dynamics of non-deuterated proteins on picoseconds-nanosecond timescales.     

Effect of cofactor binding and loop conformation on side chain methyl dynamics in dihydrofolate reductase

Dihydrofolate reductase (DHFR) has several flexible active site loops that facilitate ligand binding and catalysis. Previous studies of backbone dynamics in several complexes of DHFR indicate that the time scale and amplitude of motion depend on the conformation of the active site loops. In this study, information on dynamics is extended to methyl-containing side chains. To understand the role of side chain dynamics in ligand binding and loop conformation, methyl deuterium relaxation rates of Escherichia coli DHFR in binary folate and ternary folate:NADP+ complexes have been measured, together with chi(1) rotamer populations for threonine, isoleucine, and valine residues, determined from measurements of 3J(CgammaCO) and 3J(CgammaN) coupling constants. The results indicate that, in addition to backbone motional restriction in the adenosine-binding site, side chain flexibility in the active site and the surrounding active site loops is diminished upon binding NADP+. Resonances for several methyls in the active site and the surrounding active site loops were severely broadened in the folate:NADP+ ternary complex, suggesting the presence of motion on the chemical shift time scale. The side chains of Ile14 and Ile94, which pack against the nicotinamide and pterin rings of the cofactor and substrate, respectively, exhibit rotamer disorder in the ternary folate:NADP+ complex. Conformational fluctuations of these side chains may play a role in transition state stabilization; the observed line broadening for Ile14 suggests motions on a microsecond/millisecond time scale.     

Frequent side chain methyl carbon‐oxygen hydrogen bonding in proteins revealed by computational and stereochemical analysis of neutron structures

The propensity of backbone Cα atoms to engage in carbon‐oxygen (CH···O) hydrogen bonding is well‐appreciated in protein structure, but side chain CH···O hydrogen bonding remains largely uncharacterized. The extent to which side chain methyl groups in proteins participate in CH···O hydrogen bonding is examined through a survey of neutron crystal structures, quantum chemistry calculations, and molecular dynamics simulations. Using these approaches, methyl groups were observed to form stabilizing CH···O hydrogen bonds within protein structure that are maintained through protein dynamics and participate in correlated motion. Collectively, these findings illustrate that side chain methyl CH···O hydrogen bonding contributes to the energetics of protein structure and folding. Proteins 2015; 83:403–410. © 2014 Wiley Periodicals, Inc.     

On the ability of molecular dynamics force fields to recapitulate NMR derived protein side chain order parameters

Molecular dynamics (MD) simulations have become a central tool for investigating various biophysical questions with atomistic detail. While many different proxies are used to qualify MD force fields, most are based on largely structural parameters such as the root mean square deviation from experimental coordinates or nuclear magnetic resonance (NMR) chemical shifts and residual dipolar couplings. NMR derived Lipari–Szabo squared generalized order parameter (O²) values of amide N? H bond vectors of the polypeptide chain were also often employed for refinement and validation. However, with a few exceptions, side chain methyl symmetry axis order parameters have not been incorporated into experimental reference sets. Using a test set of five diverse proteins, the performance of several force fields implemented in the NAMDD simulation package was examined. It was found that simulations employing explicit water implemented using the TIP3 model generally performed significantly better than those using implicit water in reproducing experimental methyl symmetry axis O² values. Overall the CHARMM27 force field performs nominally better than two implementations of the Amber force field. It appeared that recent quantum mechanics modifications to side chain torsional angles of leucine and isoleucine in the Amber force field have significantly hindered proper motional modeling for these residues. There remained significant room for improvement as even the best correlations of experimental and simulated methyl group Lipari–Szabo generalized order parameters fall below an R² of 0.8.     

Main chain and side chain dynamics of a heme protein: 15N and 2H NMR relaxation studies of R. capsulatus ferrocytochrome c2

A detailed characterization of the main chain and side chain dynamics in R. capsulatus ferrocytochrome c(2) derived from (2)H NMR relaxation of methyl group resonances is presented. (15)N relaxation measurements confirm earlier results indicating that R. capsulatus ferrocytochrome c(2) exhibits minor rotational anisotropy in solution. The current study is focused on the use of deuterium relaxation in side chain methyl groups, which has been shown to provide a detailed and accurate measure of internal dynamics. Results obtained indicate that the side chains of ferrocytochrome c(2) exhibit a wide range of motional amplitudes, but are more rigid than generally found in the interior of nonprosthetic group bearing globular proteins. This unusual rigidity is ascribed to the interactions of the protein with the large heme prosthetic group. This observation has significant implications for the potential of the heme-protein interface to modulate the redox properties of the protein and also points to the need for great precision in the design and engineering of heme proteins.     

Increased rigidity of eglin c at acidic pH: evidence from NMR spin relaxation and MD simulations

To gain physical insights into how proteins respond to changes in pH, the picosecond to nanosecond time scale dynamics of the small serine protease inhibitor eglin c have been studied by NMR spin relaxation experiments and MD simulations under two pH solution conditions, pH 7 and 3. Like many proteins, eglin c is destabilized by a lowering of the pH, although it retains enough stability to maintain its native conformation at pH 3. Backbone (15)N relaxation results show comparable global tumbling times (tau(m)) and model-free order parameters (S(2)) under the two pH conditions, indicating that the molecule maintains its overall molecular shape and structure at low pH, although the backbone rigidity is slightly increased (/ = 0.6%). In contrast, the side-chain methyl dynamics, as measured from (2)H relaxation experiments, show a substantial increase in rigidity at lower pH (/ = 14.8%). Molecular dynamics simulations performed at these pH states produce results consistent with NMR measurements, showing that the two methods are in qualitative agreement. Although a full accounting of the physical basis for the concurrent conformational rigidification and destabilization at low pH requires further investigation, the high level of detail in the MD simulations provides a potential molecular mechanism: the breaking of the hydrogen bond between the side chains of Asp46 and Arg53, and changes in electrostatic interactions, appear to allow the binding loop to move closer to the core part of the protein, resulting in a more compact structure at low pH. This more compact structure may be responsible for the increased level of restriction of molecular motion. As these findings show, the stability of a molecular structure is distinct from its conformational rigidity, and the two can even change in opposite directions, against na?ve expectation.     

Molecular Dynamics of Anthraquinone DNA Intercalators with Polyethylene Glycol Side Chains

Abstract

The interactions of a homologous series of four anthraquinone (AQ) intercalators with increasing lengths of polyethylene glycol (PEG) side chains with DNA have been studied via molecular dynamics (MD) simulations. The geometry, conformation, interactions, and hydration of the complexes were examined. The geometries of the four ligands were similar with parallel stacking to the long axis of the adjacent DNA base pairs. Hydrogen bonding between the AQ amide and DNA led to a preference for the trans-syn conformer. As the side chain lengthened, binding to DNA reduced the conformational space, resulting in an increase in unfavorable entropy. Increased localization of the PEG side chain in the DNA groove, indicating some interaction of the side chain with DNA, also contributed unfavorably to the entropy. The changes in free energy of binding due to entropic considerations (— 3.9 to—6.3 kcal/mol) of AQ I-IV were significant. The hydration of the PEG side chain decreased upon binding to DNA. Understanding of side chain conformations, interactions, and hydration changes that accompany the formation of a ligand-DNA complex may be important in the development of new applications of pegylated small molecules that target biological macromolecules.     

Methyl dynamics of the amyloid-beta peptides Abeta40 and Abeta42

To probe the role of side chain dynamics in Abeta aggregation, we studied the methyl dynamics of native Abeta40 and Abeta42 by measuring cross relaxation rates with interleaved data collection. The methyl groups in the C-terminus are in general more rigid in Abeta42 than in Abeta40, consistent with previous results from backbone (15)N dynamics. This lends support to the hypothesis that a rigid C-terminus in Abeta42 may serve as an internal aggregation seed. Interestingly, two methyl groups of V18 located in the central hydrophobic cluster are more mobile in Abeta42 than in Abeta40, most likely due to the paucity of V18 intra-molecular interactions in Abeta42. V18 may then be more available for inter-molecular interactions to form Abeta42 aggregates. Thus, the side chain mobility of the central hydrophobic cluster may play an important role in Abeta aggregation and may contribute to the difference in aggregation propensity between Abeta40 and Abeta42.     

The dynamical response of hen egg white lysozyme to the binding of a carbohydrate ligand

It has become clear that the binding of small and large ligands to proteins can invoke significant changes in side chain and main chain motion in the fast picosecond to nanosecond timescale. Recently, the use of a "dynamical proxy" has indicated that changes in these motions often reflect significant changes in conformational entropy. These entropic contributions are sometimes of the same order as the total entropy of binding. Thus, it is important to understand the connections amongst motion between the manifold of states accessible to the native state of proteins, the corresponding entropy, and how this impacts the overall energetics of protein function. The interaction of proteins with carbohydrate ligands is central to a range of biological functions. Here, we examine a classic carbohydrate interaction with an enzyme: the binding of wild-type hen egg white lysozyme (HEWL) to the natural, competitive inhibitor chitotriose. Using NMR relaxation experiments, backbone amide and side chain methyl axial order parameters were obtained across apo and chitotriose-bound HEWL. Upon binding, changes in the apparent amplitude of picosecond to nanosecond main chain and side chain motions are seen across the protein. Indeed, binding of chitotriose renders a large contiguous fraction of HEWL effectively completely rigid. Changes in methyl flexibility are most pronounced closest to the binding site, but average to only a small overall change in the dynamics across the protein. The corresponding change in conformational entropy is unfavorable and estimated to be a significant fraction of the total binding entropy.