Similarities and Differences among Protein Dynamics Studied by Variable Temperature Nuclear Magnetic Resonance Relaxation

Understanding and describing the dynamics of proteins is one of the major challenges in biology. Here, we use multifield variable-temperature NMR longitudinal relaxation (R₁) measurements to determine the hierarchical activation energies of motions of four different proteins: two small globular proteins (GB1 and the SH3 domain of α-spectrin), an intrinsically disordered protein (the C-terminus of the nucleoprotein of the Sendai virus, Sendai Ntail), and an outer membrane protein (OmpG). The activation energies map the motions occurring in the side chains, in the backbone, and in the hydration shells of the proteins. We were able to identify similarities and differences in the average motions of the proteins. We find that the NMR relaxation properties of the four proteins do share similar features. The data characterizing average backbone motions are found to be very similar, the same for methyl group rotations, and similar activation energies are measured. The main observed difference occurs for the intrinsically disordered Sendai Ntail, where we observe much lower energy of activation for motions of protons associated with the protein-solvent interface as compared to the others. We also observe variability between the proteins regarding side chain ¹⁵N relaxation of lysine residues, with a higher activation energy observed in OmpG. This hints at strong interactions with negatively charged lipids in the bilayer and provides a possible mechanistic clue for the "positive-inside" rule for helical membrane proteins. Overall, these observations refine the understanding of the similarities and differences between hierarchical dynamics in proteins.

Rotamer Dynamics: Analysis of Rotamers in Molecular Dynamics Simulations of Proteins

Given by χ torsional angles, rotamers describe the side-chain conformations of amino acid residues in a protein based on the rotational isomers (hence the word rotamer). Constructed rotamer libraries, based on either protein crystal structures or dynamics studies, are the tools for classifying rotamers (torsional angles) in a way that reflect their frequency in nature. Rotamer libraries are routinely used in structure modeling and evaluation. In this perspective article, we would like to encourage researchers to apply rotamer analyses beyond their traditional use. Molecular dynamics (MD) of proteins highlight the in silico behavior of molecules in solution and thus can identify favorable side-chain conformations. In this article, we used simple computational tools to study rotamer dynamics (RD) in MD simulations. First, we isolated each frame in the MD trajectories in separate Protein Data Bank files via the cpptraj module in AMBER. Then, we extracted torsional angles via the Bio3D module in R language. The classification of torsional angles was also done in R according to the penultimate rotamer library. RD analysis is useful for various applications such as protein folding, study of rotamer-rotamer relationship in protein-protein interaction, real-time correlation between secondary structures and rotamers, study of flexibility of side chains in binding site for molecular docking preparations, use of RD as guide in functional analysis and study of structural changes caused by mutations, providing parameters for improving coarse-grained MD accuracy and speed, and many others. Major challenges facing RD to emerge as a new scientific field involve the validation of results via easy, inexpensive wet-lab methods. This realm is yet to be explored.

Dynamics, Conformational Entropy, and Frustration in Protein-Protein Interactions Involving an Intrinsically Disordered Protein Domain

Intrinsically disordered proteins (IDPs) are abundant in the eukaryotic proteome. However, little is known about the role of subnanosecond dynamics and the conformational entropy that it represents in protein-protein interactions involving IDPs. Using nuclear magnetic resonance side chain and backbone relaxation, stopped-flow kinetics, isothermal titration calorimetry, and computational studies, we have characterized the interaction between the globular TAZ1 domain of the CREB binding protein and the intrinsically disordered transactivation domain of STAT2 (TAD-STAT2). We show that the TAZ1/TAD-STAT2 complex retains considerable subnanosecond motions, with TAD-STAT2 undergoing only a partial disorder-to-order transition. We report here the first experimental determination of the conformational entropy change for both binding partners in an IDP binding interaction and find that the total change even exceeds in magnitude the binding enthalpy and is comparable to the contribution from the hydrophobic effect, demonstrating its importance in the binding energetics. Furthermore, we show that the conformational entropy change for TAZ1 is also instrumental in maintaining a biologically meaningful binding affinity. Strikingly, a spatial clustering of very high amplitude motions and a cluster of more rigid sites in the complex exist, which through computational studies we found to overlap with regions that experience energetic frustration and are less frustrated, respectively. Thus, the residual dynamics in the bound state could be necessary for faster dissociation, which is important for proteins that interact with multiple binding partners.

Dynamics of the Peripheral Membrane Protein P2 from Human Myelin Measured by Neutron Scattering--A Comparison between Wild-Type Protein and a Hinge Mutant

Myelin protein P2 is a fatty acid-binding structural component of the myelin sheath in the peripheral nervous system, and its function is related to its membrane binding capacity. Here, the link between P2 protein dynamics and structure and function was studied using elastic incoherent neutron scattering (EINS). The P38G mutation, at the hinge between the β barrel and the α-helical lid, increased the lipid stacking capacity of human P2 in vitro, and the mutated protein was also functional in cultured cells. The P38G mutation did not change the overall structure of the protein. For a deeper insight into P2 structure-function relationships, information on protein dynamics in the 10 ps to 1 ns time scale was obtained using EINS. Values of mean square displacements mainly from protein H atoms were extracted for wild-type P2 and the P38G mutant and compared. Our results show that at physiological temperatures, the P38G mutant is more dynamic than the wild-type P2 protein, especially on a slow 1-ns time scale. Molecular dynamics simulations confirmed the enhanced dynamics of the mutant variant, especially within the portal region in the presence of bound fatty acid. The increased softness of the hinge mutant of human myelin P2 protein is likely related to an enhanced flexibility of the portal region of this fatty acid-binding protein, as well as to its interactions with the lipid bilayer surface requiring conformational adaptations.

Is a malleable protein necessarily highly dynamic? The hydrophobic core of the nuclear coactivator binding domain is well ordered

The nuclear coactivator binding domain of CREB binding protein folds into remarkably different structures in complex with different ligands. To understand the mechanism of the structural adaptability in the nuclear coactivator binding domain (NCBD), we have compared the dynamics of the hydrophobic core of NCBD in the ligand-free state and in a well-folded complex with the ligand activator for thyroid hormone and retinoid receptors using multiple NMR methods including methyl chemical shifts, coupling constants, and methyl order parameters. From all NMR measures, the aliphatic side chains in the hydrophobic core are slightly more dynamic in the free protein than in the complex, but have mobility comparable to the hydrophobic cores of average folded proteins. Urea titration monitored by NMR reveals that all parts of the protein, including the side-chain packing in the hydrophobic core, denatures in a single cooperative process. The molten globule characteristics of NCBD are thus restricted to a slowly fluctuating tertiary structure. Consequently, the conformational plasticity of the protein is most likely related to its low overall stability rather than an intrinsically flexible protein structure. The well-defined structure supports a model of molecular recognition dominated by conformational selection, whereas only minor structural adjustments are necessary after the association.