Local isotropic diffusion approximation for coupled internal and overall molecular motions in NMR spin relaxation

Bromodomain Interactions with Acetylated Histone 4 Peptides in the BRD4 Tandem Domain: Effects on Domain Dynamics and Internal Flexibility

The bromodomain and extra-terminal (BET) protein BRD4 regulates gene expression via recruitment of transcriptional regulatory complexes to acetylated chromatin. Like other BET proteins, BRD4 contains two bromodomains, BD1 and BD2, that can interact cooperatively with target proteins and designed ligands, with important implications for drug discovery. Here, we used nuclear magnetic resonance (NMR) spectroscopy to study the dynamics and interactions of the isolated bromodomains, as well as the tandem construct including both domains and the intervening linker, and investigated the effects of binding a tetra-acetylated peptide corresponding to the tail of histone 4. The peptide affinity is lower for both domains in the tandem construct than for the isolated domains. Using ¹⁵N spin relaxation, we determined the global rotational correlation times and residue-specific order parameters for BD1 and BD2. Isolated BD1 is monomeric in the apo state but apparently dimerizes upon binding the tetra-acetylated peptide. Isolated BD2 partially dimerizes in both the apo and peptide-bound states. The backbone order parameters reveal marked differences between BD1 and BD2, primarily in the acetyl-lysine binding site where the ZA loop is more flexible in BD2. Peptide binding reduces the order parameters of the ZA loop in BD1 and the ZA and BC loops in BD2. The AB loop, located distally from the binding site, shows variable dynamics that reflect the different dimerization propensities of the domains. These results provide a basis for understanding target recognition by BRD4.

Explicit models of motions to analyze NMR relaxation data in proteins

Nuclear Magnetic Resonance (NMR) is a tool of choice to characterize molecular motions. In biological macromolecules, pico- to nano-second motions, in particular, can be probed by nuclear spin relaxation rates which depend on the time fluctuations of the orientations of spin interaction frames. For the past 40 years, relaxation rates have been successfully analyzed using the Model Free (MF) approach which makes no assumption on the nature of motions and reports on the effective amplitude and time-scale of the motions. However, obtaining a mechanistic picture of motions from this type of analysis is difficult at best, unless complemented with molecular dynamics (MD) simulations. In spite of their limited accuracy, such simulations can be used to obtain the information necessary to build explicit models of motions designed to analyze NMR relaxation data. Here, we present how to build such models, suited in particular to describe motions of methyl-bearing protein side-chains and compare them with the MF approach. We show on synthetic data that explicit models of motions are more robust in the presence of rotamer jumps which dominate the relaxation in methyl groups of protein side-chains. We expect this work to motivate the use of explicit models of motion to analyze MD and NMR data.

Analysis of NMR Spin-Relaxation Data Using an Inverse Gaussian Distribution Function

Spin relaxation in solution-state NMR spectroscopy is a powerful approach to explore the conformational dynamics of biological macromolecules. Probability distribution functions for overall or internal correlation times have been used previously to model spectral density functions central to spin-relaxation theory. Applications to biological macromolecules rely on transverse relaxation rate constants, and when studying nanosecond timescale motions, sampling at ultralow frequencies is often necessary. Consequently, appropriate distribution functions necessitate spectral density functions that are accurate and convergent as frequencies approach zero. In this work, the inverse Gaussian probability distribution function is derived from general properties of spectral density functions at low and high frequencies for macromolecules in solution, using the principle of maximal entropy. This normalized distribution function is first used to calculate the correlation function, followed by the spectral density function. The resulting model-free spectral density functions are finite at a frequency of zero and can be used to describe distributions of either overall or internal correlation times using the model-free ansatz. To validate the approach, ¹⁵N spin-relaxation data for the bZip transcription factor domain of the Saccharomyces cerevisiae protein GCN4, in the absence of cognate DNA, were analyzed using the inverse Gaussian probability distribution for intramolecular correlation times. The results extend previous models for the conformational dynamics of the intrinsically disordered, DNA-binding region of the bZip transcription factor domain.

A dynamic look backward and forward

The 2015 Gunther Laukien Prize recognized solution NMR studies of protein dynamics and thermodynamics. This Perspective surveys aspects of the development and application of NMR spin relaxation for investigations of protein flexibility and function over multiple time scales in solution. Methods highlighted include analysis of overall rotational diffusion, theoretical descriptions of R1ρ relaxation, and molecular dynamics simulations to interpret NMR spin relaxation. Applications are illustrated for the zinc-finger domain Xfin-31, the calcium-binding proteins calbindin D9k and calmodulin, and the bZip transcription factor of GCN4.

Dynamics of GCN4 facilitate DNA interaction: a model-free analysis of an intrinsically disordered region

Intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered regions (IDRs) are known to play important roles in regulatory and signaling pathways. A critical aspect of these functions is the ability of IDP/IDRs to form highly specific complexes with target molecules. However, elucidation of the contributions of conformational dynamics to function has been limited by challenges associated with structural heterogeneity of IDP/IDRs. Using NMR spin relaxation parameters ((15)N R1, (15)N R2, and {(1)H}-(15)N heteronuclear NOE) collected at four static magnetic fields ranging from 14.1 to 21.1 T, we have analyzed the backbone dynamics of the basic leucine-zipper (bZip) domain of the Saccharomyces cerevisiae transcription factor GCN4, whose DNA binding domain is intrinsically disordered in the absence of DNA substrate. We demonstrate that the extended model-free analysis can be applied to proteins with IDRs such as apo GCN4 and that these results significantly extend previous NMR studies of GCN4 dynamics performed using a single static magnetic field of 11.74 T [Bracken, et al., J. Mol. Biol., 1999, 285, 2133-2146] and correlate well with molecular dynamics simulations [Robustelli, et al., J. Chem. Theory Comput., 2013, 9, 5190-5200]. In contrast to the earlier work, data at multiple static fields allows the time scales of internal dynamics of GCN4 to be reliably quantified. Large amplitude dynamic fluctuations in the DNA-binding region have correlation times (τs ≈ 1.4-2.5 ns) consistent with a two-step mechanism in which partially ordered bZip conformations of GCN4 form initial encounter complexes with DNA and then rapidly rearrange to the high affinity state with fully formed basic region recognition helices.

Coupling between overall rotational diffusion and domain motions in proteins and its effect on dielectric spectra

In this work, we formulate a closed-form solution of the model of a semirigid molecule for the case of fluctuating and reorienting molecular electric dipole moment. We illustrate with numeric calculations the impact of protein domain motions on dielectric spectra using the example of the 128 kDa protein dimer of Enzyme I. We demonstrate that the most drastic effect occurs for situations when the characteristic time of protein domain dynamics is comparable to the time of overall molecular rotational diffusion. We suggest that protein domain motions could be a possible explanation for the high-frequency contribution that accompanies the major relaxation dispersion peak in the dielectric spectra of protein aqueous solutions. We propose that the presented computational methodology could be used for the simultaneous analysis of dielectric spectroscopy and nuclear magnetic resonance data. Proteins 2015; 83:1571-1581. © 2015 Wiley Periodicals, Inc.