A Unique and Simple Approach to Improve Sensitivity in 15N-NMR Relaxation Measurements for NH₃⁺ Groups: Application to a Protein-DNA Complex

Parsing Dynamics of Protein Backbone NH and Side-Chain Methyl Groups using Molecular Dynamics Simulations

Experimental NMR spectroscopy and theoretical molecular dynamics (MD) simulations provide complementary insights into protein conformational dynamics and hence into biological function. The present work describes an extensive set of backbone NH and side-chain methyl group generalized order parameters for the Escherichia coli ribonuclease HI (RNH) enzyme derived from 2-μs microsecond MD simulations using the OPLS4 and AMBER-FF19SB force fields. The simulated generalized order parameters are compared with values derived from NMR 15N and 13CH2D spin relaxation measurements. The squares of the generalized order parameters, S2 for the N-H bond vector and Saxis2 for the methyl group symmetry axis, characterize the equilibrium distribution of vector orientations in a molecular frame of reference. Optimal agreement between simulated and experimental results was obtained by averaging S2 or Saxis2 calculated by dividing the simulated trajectories into 50 ns blocks (∼five times the rotational diffusion correlation time for RNH). With this procedure, the median absolute deviations (MAD) between experimental and simulated values of S2 and Saxis2 are 0.030 (NH) and 0.061 (CH3) for OPLS4 and 0.041 (NH) and 0.078 (CH3) for AMBER-FF19SB. The MAD between OPLS4 and AMBER-FF19SB are 0.021 (NH) and 0.072 (CH3). The generalized order parameters for the methyl group symmetry axis can be decomposed into contributions from backbone fluctuations, between-rotamer dihedral angle transitions, and within-rotamer dihedral angle fluctuations. Analysis of the simulation trajectories shows that (i) backbone and side chain conformational fluctuations exhibit little correlation and that (ii) fluctuations within rotamers are limited and highly uniform with values that depend on the number of dihedral angles considered. Low values of Saxis2, indicative of enhanced side-chain flexibility, result from between-rotamer transitions that can be enhanced by increased local backbone flexibility.

NMR Methods for Characterizing the Basic Side Chains of Proteins: Electrostatic Interactions, Hydrogen Bonds, and Conformational Dynamics

NMR spectroscopy is a powerful tool for studying protein dynamics. Conventionally, NMR studies on protein dynamics have probed motions of protein backbone NH, side-chain aromatic, and CH₃ groups. Recently, there has been remarkable progress in NMR methodologies that can characterize motions of cationic groups in protein side chains. These NMR methods allow investigations of the dynamics of positively charged lysine (Lys) and arginine (Arg) side chains and their hydrogen bonds as well as their electrostatic interactions important for protein function. Here, describing various practical aspects, we provide an overview of the NMR methods for dynamics studies of Lys and Arg side chains. Some example data on protein-DNA complexes are shown. We will also explain how molecular dynamics (MD) simulations can facilitate the interpretation of the NMR data on these basic side chains. Studies combining NMR and MD have revealed the highly dynamic nature of short-range electrostatic interactions via ion pairs, especially those involving Lys side chains.

Lysine Side-Chain Dynamics in the Binding Site of Homeodomain/DNA Complexes As Observed by NMR Relaxation Experiments and Molecular Dynamics Simulations

An important but poorly characterized contribution to the thermodynamics of protein-DNA interactions is the loss of entropy that occurs from restricting the conformational freedom of amino acid side chains. The effect of restricting the flexibility of several side chains at a protein-DNA interface may be comparable in many cases to the other factors that determine the binding thermodynamics and may, therefore, play a key role in dictating the binding affinity and/or specificity. Because the entropic contributions, including the presence and influence of side-chain dynamics, are especially difficult to estimate based on structural information, it is important to pursue experimental and theoretical studies that can provide direct information regarding these issues. We report on studies of a model system, the homeodomain/DNA complex, focusing on the Lys50 class of homeodomains where a key lysine residue in position 50 was shown previously to be critical for binding site specificity. NMR methodology was employed for determining the dynamics of lysine side-chain amino groups via 15N relaxation measurements in the Lys50-class homeodomains from the Drosophila protein Bicoid and the human protein Pitx2. In the case of Pitx2, complexes with both a consensus and a nonconsensus DNA binding site were examined. NMR-derived order parameters indicated moderate to substantial conformational freedom for the lysine NH3+ group in the complexes studied. To complement the experimental NMR measurements, molecular dynamics simulations were performed for the consensus complexes to gain further, detailed insights regarding the dynamics of the Lys50 side chain and other important residues in the protein-DNA interface.

Internal Motions of Basic Side Chains of the Antennapedia Homeodomain in the Free and DNA-Bound States

Basic side chains play crucial roles in protein-DNA interactions. In this study, using NMR spectroscopy, we investigated the dynamics of Arg and Lys side chains of the fruit fly Antennapedia homeodomain in the free state and in the complex with target DNA. We measured 15N relaxation for Arg and Lys side chains at two magnetic fields, from which generalized order parameters for the cationic groups were determined. Mobility of the R5 side chain, which makes hydrogen bonds with a thymine base in the DNA minor groove, was greatly dampened. Several Lys and Arg side chains that form intermolecular ion pairs with DNA phosphates were found to retain high mobility with the order parameter being <0.6 in the DNA-bound state. Interestingly, some of the interfacial cationic groups in the complex were more mobile than in the free protein. The retained or enhanced mobility of the Arg and Lys side chains in the complex should mitigate the overall loss of conformational entropy in the protein-DNA association and allow dynamic molecular recognition.