Dynamic and thermodynamic response of the Ras protein Cdc42Hs upon association with the effector domain of PAK3

Pressure, motion, and conformational entropy in molecular recognition by proteins

The thermodynamics of molecular recognition by proteins is a central determinant of complex biochemistry. For over a half-century, detailed cryogenic structures have provided deep insight into the energetic contributions to ligand binding by proteins. More recently, a dynamical proxy based on NMR-relaxation methods has revealed an unexpected richness in the contributions of conformational entropy to the thermodynamics of ligand binding. Here, we report the pressure dependence of fast internal motion within the ribonuclease barnase and its complex with the protein barstar. In what we believe is a first example, we find that protein dynamics are conserved along the pressure-binding thermodynamic cycle. The femtomolar affinity of the barnase-barstar complex exists despite a penalty by -TΔS_conf of +11.7 kJ/mol at ambient pressure. At high pressure, however, the overall change in side-chain dynamics is zero, and binding occurs with no conformational entropy penalty, suggesting an important role of conformational dynamics in the adaptation of protein function to extreme environments. Distinctive clustering of the pressure sensitivity is observed in response to both pressure and binding, indicating the presence of conformational heterogeneity involving less efficiently packed alternative conformation(s). The structural segregation of dynamics observed in barnase is striking and shows how changes in both the magnitude and the sign of regional contributions of conformational entropy to the thermodynamics of protein function are possible.

Early-stage structure-based drug discovery for small GTPases by NMR spectroscopy

Small GTPase or Ras superfamily, including Ras, Rho, Rab, Ran and Arf, are fundamental in regulating a wide range of cellular processes such as growth, differentiation, migration and apoptosis. They share structural and functional similarities for binding guanine nucleotides and hydrolyzing GTP. Dysregulations of Ras proteins are involved in the pathophysiology of multiple human diseases, however there is still a stringent need for effective treatments targeting these proteins. For decades, small GTPases were recognized as 'undruggable' targets due to their complex regulatory mechanisms and lack of deep pockets for ligand binding. NMR has been critical in deciphering the structural and dynamic properties of the switch regions that are underpinning molecular switch functions of small GTPases, which pave the way for developing new effective inhibitors. The recent progress of drug or lead molecule development made for small GTPases profoundly delineated how modern NMR techniques reshape the field of drug discovery. In this review, we will summarize the progress of structural and dynamic studies of small GTPases, the NMR techniques developed for structure-based drug screening and their applications in early-stage drug discovery for small GTPases.

Inter-Active Site Communication Mediated by the Dimer Interface β-Sheet in the Half-the-Sites Enzyme, Thymidylate Synthase

Thymidylate synthase (TS) is a dimeric enzyme conserved in all life forms that exhibits the allosteric feature of half-the-sites activity. Neither the reason for nor the mechanism of this phenomenon is understood. We used a combined nuclear magnetic resonance (NMR) and molecular dynamics approach to study a stable intermediate preceding hydride transfer, which is the rate-limiting and half-the-sites step. In NMR titrations with ligands leading to this intermediate, we measured chemical shifts of the apoenzyme (lig₀), the saturated holoenzyme (lig₂), and the typically elusive singly bound (lig₁) states. Approximately 40 amides showed quartet patterns providing direct NMR evidence of coupling between the active site and probes >30 Å away in the distal subunit. Quartet peak patterns have symmetrical character, indicating reciprocity in communicating the first and second binding events to the distal protomer. Quartets include key catalytic residues and map to the dimer interface β-sheet, which also represents the shortest path between the two active sites. Simulations corroborate the coupling observed in solution in that there is excellent overlap between quartet residues and main-chain atoms having intersubunit cross-correlated motions. Simulations identify five hot spot residues, three of which lie at the kink in the unique β-bulge abutting the active sites on either end of the sheet. Interstrand cross-correlated motions become more organized and pronounced as the enzyme progresses from lig₀ to lig₁ and ultimately lig₂. Coupling in the apparently symmetrical complex has implications for half-the-sites reactivity and potentially resolves the paradox of inequivalent TS active sites despite the vast majority of X-ray structures appearing to be symmetrical.

Allostery and dynamics in small G proteins

The Ras family of small guanine nucleotide-binding proteins behave as molecular switches: they are switched off and inactive when bound to GDP but can be activated by GTP binding in response to signal transduction pathways. Early structural analysis showed that two regions of the protein, which change conformation depending on the nucleotide present, mediate this switch. A large number of X-ray, NMR and simulation studies have shown that this is an over-simplification. The switch regions themselves are highly dynamic and can exist in distinct sub-states in the GTP-bound form that have different affinities for other proteins. Furthermore, regions outside the switches have been found to be sensitive to the nucleotide state of the protein, indicating that allosteric change is more widespread than previously thought. Taken together, the accrued knowledge about small G protein structures, allostery and dynamics will be essential for the design and testing of the next generation of inhibitors, both orthosteric and allosteric, as well as for understanding their mode of action.

Site-directed spin label electron paramagnetic resonance spectroscopy as a probe of conformational dynamics in the Fe(III) "locked-off" state of the CO-sensing transcription factor CooA

The transcriptional activator CooA belongs to the CRP/FNR (cAMP receptor protein/fumarate and nitrate reductase) superfamily of transcriptional regulators and uses heme to sense carbon monoxide (CO). Effector-driven allosteric activation is well understood in CRP, a CooA homologue. A structural allosteric activation model for CooA exists which parallels that of CRP; however, the role of protein dynamics, which is crucial in CRP, is not well understood in CooA. We employed site-directed spin labeling electron paramagnetic resonance spectroscopy to probe CooA motions on the μs-ms timescale. We created a series of Cys substitution variants, each with a cysteine residue introduced into a key functional region of the protein: K26C, E60C, F132C, D134C, and S175C. The heme environment and DNA binding affinity of each variant were comparable to those of wild-type CooA, with the exception of F132C, which displayed reduced DNA binding affinity. This observation confirms a previously hypothesized role for Phe132 in transmitting the allosteric CO binding signal. Osmolyte perturbation studies of Fe(III) "locked-off" CooA variants labeled with either MTSL or MAL-6 nitroxide spin labels revealed that multicomponent EPR spectra report on conformational flexibility on the μs-ms timescale. Multiple dynamic populations exist at every site examined in the structurally uncharacterized Fe(III) "locked-off" CooA. This observation suggests that, in direct contrast to effector-free CRP, Fe(III) "locked-off" CooA undergoes conformational exchange on the μs-ms timescale. Importantly, we establish MAL-6 as a spin label with a redox-stable linkage that may be utilized to compare conformational dynamics between functional states of CooA.

Entropy in molecular recognition by proteins

Molecular recognition by proteins is fundamental to molecular biology. Dissection of the thermodynamic energy terms governing protein-ligand interactions has proven difficult, with determination of entropic contributions being particularly elusive. NMR relaxation measurements have suggested that changes in protein conformational entropy can be quantitatively obtained through a dynamical proxy, but the generality of this relationship has not been shown. Twenty-eight protein-ligand complexes are used to show a quantitative relationship between measures of fast side-chain motion and the underlying conformational entropy. We find that the contribution of conformational entropy can range from favorable to unfavorable, which demonstrates the potential of this thermodynamic variable to modulate protein-ligand interactions. For about one-quarter of these complexes, the absence of conformational entropy would render the resulting affinity biologically meaningless. The dynamical proxy for conformational entropy or "entropy meter" also allows for refinement of the contributions of solvent entropy and the loss in rotational-translational entropy accompanying formation of high-affinity complexes. Furthermore, structure-based application of the approach can also provide insight into long-lived specific water-protein interactions that escape the generic treatments of solvent entropy based simply on changes in accessible surface area. These results provide a comprehensive and unified view of the general role of entropy in high-affinity molecular recognition by proteins.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

General order parameter based correlation analysis of protein backbone motions between experimental NMR relaxation measurements and molecular dynamics simulations

Internal backbone dynamic motions are essential for different protein functions and occur on a wide range of time scales, from femtoseconds to seconds. Molecular dynamic (MD) simulations and nuclear magnetic resonance (NMR) spin relaxation measurements are valuable tools to gain access to fast (nanosecond) internal motions. However, there exist few reports on correlation analysis between MD and NMR relaxation data. Here, backbone relaxation measurements of (15)N-labeled SH3 (Src homology 3) domain proteins in aqueous buffer were used to generate general order parameters (S(2)) using a model-free approach. Simultaneously, 80 ns MD simulations of SH3 domain proteins in a defined hydrated box at neutral pH were conducted and the general order parameters (S(2)) were derived from the MD trajectory. Correlation analysis using the Gromos force field indicated that S(2) values from NMR relaxation measurements and MD simulations were significantly different. MD simulations were performed on models with different charge states for three histidine residues, and with different water models, which were SPC (simple point charge) water model and SPC/E (extended simple point charge) water model. S(2) parameters from MD simulations with charges for all three histidines and with the SPC/E water model correlated well with S(2) calculated from the experimental NMR relaxation measurements, in a site-specific manner.