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

Negatively charged, intrinsically disordered regions can accelerate target search by DNA-binding proteins

In eukaryotes, many DNA/RNA-binding proteins possess intrinsically disordered regions (IDRs) with large negative charge, some of which involve a consecutive sequence of aspartate (D) or glutamate (E) residues. We refer to them as D/E repeats. The functional role of D/E repeats is not well understood, though some of them are known to cause autoinhibition through intramolecular electrostatic interaction with functional domains. In this work, we investigated the impacts of D/E repeats on the target DNA search kinetics for the high-mobility group box 1 (HMGB1) protein and the artificial protein constructs of the Antp homeodomain fused with D/E repeats of varied lengths. Our experimental data showed that D/E repeats of particular lengths can accelerate the target association in the overwhelming presence of non-functional high-affinity ligands ('decoys'). Our coarse-grained molecular dynamics (CGMD) simulations showed that the autoinhibited proteins can bind to DNA and transition into the uninhibited complex with DNA through an electrostatically driven induced-fit process. In conjunction with the CGMD simulations, our kinetic model can explain how D/E repeats can accelerate the target association process in the presence of decoys. This study illuminates an unprecedented role of the negatively charged IDRs in the target search process.

Observation of Arginine Side-Chain Motions Coupled to the Global Conformational Exchange Process in Deubiquitinase A

Coupled motions have been demonstrated to be functionally important in a number of enzymes. Noncovalent side-chain interactions play essential roles in coordinating the motions across different structural elements in a protein. However, most of the dynamic studies of proteins are focused on backbone amides or methyl groups in the side chains and little is known about the polar and charged side chains. We have previously characterized the conformational dynamics of deubiquitinase A (DUBA), an isopeptidase, on the microsecond-to-millisecond (μs-ms) time scales with the amide ¹H Carr-Purcell-Meiboom-Gill (CPMG) experiment. We detected a global conformational exchange process on a time scale of approximately 200 μs, which involves most of the structural elements in DUBA, including the active site and the substrate binding interface. Here, we extend our previous study on backbone amides to the arginine side-chain N^ε-H^ε groups using a modified ¹H CPMG pulse sequence that can efficiently detect both backbone amide and arginine side-chain N^ε-H^ε signals in a single experiment. We found that the side chains of three arginines display motions on the same time scale as the backbone amides. Mutations of two of the three arginines to alanines result in a decrease in enzyme activity. One of these two arginines is located in a loop involved in substrate binding. This loop is not visible in the backbone amide-detected experiments due to excess line broadening induced by motions on the μs-ms time scales. These results clearly demonstrate that the motions of some arginine side chains are coupled to the global conformational exchange process and provide an additional probe for motions in a functionally important loop that did not yield visible backbone amide signals, suggesting the value of side-chain experiments on DUBA. The modified ¹H CPMG pulse sequence allows the simultaneous characterization of backbone and arginine side-chain dynamics without any increase in data acquisition time and can be applied to the dynamic studies of any protein that displays measurable amide ¹H relaxation dispersion.

Dynamics of Cations around DNA and Protein as Revealed by 23Na Diffusion NMR Spectroscopy

Counterions are vital for the structure and function of biomolecules. However, the behavior of counterions remains elusive due to the difficulty in characterizing mobile ions. Here, we demonstrate that the dynamics of cations around biological macromolecules can be revealed by ²³Na diffusion nuclear magnetic resonance (NMR) spectroscopy. NMR probe hardware capable of generating strong magnetic field gradients enables ²³Na NMR-based diffusion measurements for Na⁺ ions in solutions of biological macromolecules and their complexes. The dynamic properties of Na⁺ ions interacting with the macromolecules can be investigated using apparent ²³Na diffusion coefficients measured under various conditions. Our diffusion data clearly show that Na⁺ ions retain high mobility within the ion atmosphere around DNA. The ²³Na diffusion NMR method also permits direct observation of the release of Na⁺ ions from nucleic acids upon protein-nucleic acid association. The entropy change due to the ion release can be estimated from the diffusion data.

Hindered Rotations of Protein Asparagine/Glutamine Side-Chain NH2 Groups: Impact of Hydrogen Bonding with DNA

Hindered rotation about an sp² C-N bond is known to occur in arginine (Arg), asparagine (Asn), and glutamine (Gln) side chains of proteins. However, very little is known about the rotational dynamics of Asn and Gln side-chain NH₂ groups. Here, using a unique NMR method, we quantitatively characterized the hindered rotations of protein Asn/Gln side-chain NH₂ groups. This NMR method yields simple NH₂-selective spectra that allow for an accurate determination of the kinetic rate constants for the hindered rotations. Through the NMR measurements at different temperatures, we investigated the energy barriers that restrict the C-N bond rotations of protein side-chain NH₂ groups. Through a comparison of the kinetic data for the free and DNA-bound states of the Antp homeodomain, we also examined the impact of hydrogen bonding on the hindered rotations of the side-chain NH₂ groups. Our data suggest that the hydrogen bonding increases the energy barriers by 1-6 kJ/mol.

De novo determination of near-surface electrostatic potentials by NMR

Electrostatic potentials computed from three-dimensional structures of biomolecules by solving the Poisson-Boltzmann equation are widely used in molecular biophysics, structural biology, and medicinal chemistry. Despite the approximate nature of the Poisson-Boltzmann theory, validation of the computed electrostatic potentials around biological macromolecules is rare and methodologically limited. Here, we present a unique and powerful NMR method that allows for straightforward and extensive comparison with electrostatic models for biomolecules and their complexes. This method utilizes paramagnetic relaxation enhancement arising from analogous cationic and anionic cosolutes whose spatial distributions around biological macromolecules reflect electrostatic potentials. We demonstrate that this NMR method enables de novo determination of near-surface electrostatic potentials for individual protein residues without using any structural information. We applied the method to ubiquitin and the Antp homeodomain-DNA complex. The experimental data agreed well with predictions from the Poisson-Boltzmann theory. Thus, our experimental results clearly support the validity of the theory for these systems. However, our experimental study also illuminates certain weaknesses of the Poisson-Boltzmann theory. For example, we found that the theory predicts stronger dependence of near-surface electrostatic potentials on ionic strength than observed in the experiments. Our data also suggest that conformational flexibility or structural uncertainties may cause large errors in theoretical predictions of electrostatic potentials, particularly for highly charged systems. This NMR-based method permits extensive assessment of near-surface electrostatic potentials for various regions around biological macromolecules and thereby may facilitate improvement of the computational approaches for electrostatic potentials.

Methodological advancements for characterising protein side chains by NMR spectroscopy

The surface of proteins is covered by side chains of polar amino acids that are imperative for modulating protein functionality through the formation of noncovalent intermolecular interactions. However, despite their tremendous importance, the unique structures of protein side chains require tailored approaches for investigation by nuclear magnetic resonance spectroscopy and so have traditionally been understudied compared with the protein backbone. Here, we review substantial recent methodological advancements within nuclear magnetic resonance spectroscopy to address this issue. Specifically, we consider advancements that provide new insight into methyl-bearing side chains, show the potential of using non-natural amino acids and reveal the actions of charged side chains. Combined, the new methods promise unprecedented characterisations of side chains that will further elucidate protein function.

Discrete-state stochastic kinetic models for target DNA search by proteins: Theory and experimental applications

To perform their functions, transcription factors and DNA-repair/modifying enzymes randomly search DNA in order to locate their specific targets on DNA. Discrete-state stochastic kinetic models have been developed to explain how the efficiency of the search process is influenced by the molecular properties of proteins and DNA as well as by other factors such as molecular crowding. These theoretical models not only offer explanations on the relation of microscopic processes to macroscopic behavior of proteins, but also facilitate the analysis and interpretation of experimental data. In this review article, we provide an overview on discrete-state stochastic kinetic models and explain how these models can be applied to experimental investigations using stopped-flow, single-molecule, nuclear magnetic resonance (NMR), and other biophysical and biochemical methods.

Slow Folding of a Helical Protein: Large Barriers, Strong Internal Friction, or a Shallow, Bumpy Landscape?

The rate at which a protein molecule folds is determined by opposing energetic and entropic contributions to the free energy that shape the folding landscape. Delineating the extent to which they impact the diffusional barrier-crossing events, including the magnitude of internal friction and barrier height, has largely been a challenging task. In this work, we extract the underlying thermodynamic and dynamic contributions to the folding rate of an unusually slow-folding helical DNA-binding domain, PurR, which shares the characteristics of ultrafast downhill-folding proteins but nonetheless appears to exhibit an apparent two-state equilibrium. We combine equilibrium spectroscopy, temperature-viscosity-dependent kinetics, statistical mechanical modeling, and coarse-grained simulations to show that the conformational behavior of PurR is highly heterogeneous characterized by a large spread in melting temperatures, marginal thermodynamic barriers, and populated partially structured states. PurR appears to be at the threshold of disorder arising from frustrated electrostatics and weak packing that in turn slows down folding due to a shallow, bumpy landscape and not due to large thermodynamic barriers or strong internal friction. Our work highlights how a strong temperature dependence on the pre-exponential could signal a shallow landscape and not necessarily a slow-folding diffusion coefficient, thus determining the folding timescales of even millisecond folding proteins and hints at possible structural origins for the shallow landscape.

Dynamics of Ionic Interactions at Protein-Nucleic Acid Interfaces

Molecular association of proteins with nucleic acids is required for many biological processes essential to life. Electrostatic interactions via ion pairs (salt bridges) of nucleic acid phosphates and protein side chains are crucial for proteins to bind to DNA or RNA. Counterions around the macromolecules are also key constituents for the thermodynamics of protein–nucleic acid association. Until recently, there had been only a limited amount of experiment-based information about how ions and ionic moieties behave in biological macromolecular processes. In the past decade, there has been significant progress in quantitative experimental research on ionic interactions with nucleic acids and their complexes with proteins. The highly negatively charged surfaces of DNA and RNA electrostatically attract and condense cations, creating a zone called the ion atmosphere. Recent experimental studies were able to examine and validate theoretical models on ions and their mobility and interactions with macromolecules. The ionic interactions are highly dynamic. The counterions rapidly diffuse within the ion atmosphere. Some of the ions are released from the ion atmosphere when proteins bind to nucleic acids, balancing the charge via intermolecular ion pairs of positively charged side chains and negatively charged backbone phosphates. Previously, the release of counterions had been implicated indirectly by the salt-concentration dependence of the association constant.

Recently, direct detection of counterion release by NMR spectroscopy has become possible and enabled more accurate and quantitative analysis of the counterion release and its entropic impact on the thermodynamics of protein–nucleic acid association. Recent studies also revealed the dynamic nature of ion pairs of protein side chains and nucleic acid phosphates. These ion pairs undergo transitions between two major states. In one of the major states, the cation and the anion are in direct contact and form hydrogen bonds. In the other major state, the cation and the anion are separated by water. Transitions between these states rapidly occur on a picosecond to nanosecond time scale. When proteins interact with nucleic acids, interfacial arginine (Arg) and lysine (Lys) side chains exhibit considerably different behaviors. Arg side chains show a higher propensity to form rigid contacts with nucleotide bases, whereas Lys side chains tend to be more mobile at the molecular interfaces. The dynamic ionic interactions may facilitate adaptive molecular recognition and play both thermodynamic and kinetic roles in protein–nucleic acid interactions.

Racemic phosphorothioate as a tool for NMR investigations of protein-DNA complexes

A major driving force for protein-nucleic acid association is electrostatic interactions via ion pairs of the positively charged basic side chains and negatively charged phosphates. For a better understanding of how proteins scan DNA and recognize particular signatures, it is important to gain atomic-level insight into the behavior of basic side chains at the protein-DNA interfaces. NMR spectroscopy is a powerful tool for investigating the structural, dynamic, and kinetic aspects of protein-DNA interactions. However, resonance assignment of basic side-chain cationic moieties at the molecular interfaces remains to be a major challenge. Here, we propose a fast, robust, and inexpensive approach that greatly facilitates resonance assignment of interfacial moieties and also allows for kinetic measurements of protein translocation between two DNA duplexes. This approach utilizes site-specific incorporation of racemic phosphorothioate at the position of a phosphate that interacts with a protein side chain. This modification retains the electric charge of phosphate and therefore is mild, but causes significant chemical shift perturbations for the proximal protein side chains, which facilitates resonance assignment. Due to the racemic nature of the modification, two different chemical shifts are observed for the species with different diastereomers R_P and S_P of the incorporated phosphorothioate group. Kinetic information on the exchange of the protein molecule between R_P and S_P DNA duplexes can be obtained by ¹⁵N_z exchange spectroscopy. We demonstrate the applications of this approach to the Antennapedia homeodomain-DNA complex and the CREB1 basic leucine-zipper (bZIP)-DNA complex.

NMR Observation of Intermolecular Hydrogen Bonds between Protein Tyrosine Side-Chain OH and DNA Phosphate Groups

Hydrogen bonds between protein side-chain hydroxyl (OH) and phosphate groups are one of the most common types of intermolecular hydrogen bonds in protein-DNA/RNA complexes. Using NMR spectroscopy, we identified and characterized the hydrogen bonds between tyrosine side-chain OH and DNA phosphate groups in a protein-DNA complex. These OH groups exhibited relatively slow hydrogen-exchange rates and sizable scalar couplings between hydroxyl 1H and DNA phosphate 31P nuclei across the hydrogen bonds. Information about intermolecular hydrogen bonds facilitates investigations of the DNA/RNA recognition by the protein.

Tunable order-disorder continuum in protein-DNA interactions

DNA-binding protein domains (DBDs) sample diverse conformations in equilibrium facilitating the search and recognition of specific sites on DNA over millions of energetically degenerate competing sites. We hypothesize that DBDs have co-evolved to sense and exploit the strong electric potential from the array of negatively charged phosphate groups on DNA. We test our hypothesis by employing the intrinsically disordered DBD of cytidine repressor (CytR) as a model system. CytR displays a graded increase in structure, stability and folding rate on increasing the osmolarity of the solution that mimics the non-specific screening by DNA phosphates. Electrostatic calculations and an Ising-like statistical mechanical model predict that CytR exhibits features of an electric potential sensor modulating its dimensions and landscape in a unique distance-dependent manner, while DNA plays the role of a non-specific macromolecular chaperone. Accordingly, CytR binds its natural half-site faster than the diffusion-controlled limit and even random DNA conforming to an electrostatic-steering binding mechanism. Our work unravels for the first time the synergistic features of a natural electrostatic potential sensor, a novel binding mechanism driven by electrostatic frustration and disorder, and the role of DNA in promoting distance-dependent protein structural transitions critical for switching between specific and non-specific DNA-binding modes.

Mobility of Histidine Side Chains Analyzed with 15N NMR Relaxation and Cross-Correlation Data: Insight into Zinc-Finger-DNA Interactions

Due to chemical exchange, the mobility of histidine (His) side chains of proteins is typically difficult to analyze by NMR spectroscopy. Using an NMR approach that is uninfluenced by chemical exchange, we investigated internal motions of the His imidazole NH groups that directly interact with DNA phosphates in the Egr-1 zinc-finger-DNA complex. In this approach, the transverse and longitudinal cross-correlation rates for ¹⁵N chemical shift anisotropy and ¹⁵N-¹H dipole-dipole relaxation interference were analyzed together with ¹⁵N longitudinal relaxation rates and heteronuclear Overhauser effect data at two magnetic field strengths. We found that the zinc-coordinating His side chains directly interacting with DNA phosphates are strongly restricted in mobility. This makes a contrast to the arginine and lysine side chains that retain high mobility despite their interactions with DNA phosphates in the same complex. The entropic effects of side-chain mobility on the molecular association are discussed.

Dynamic ion pair behavior stabilizes single α-helices in proteins

Ion pairs are key stabilizing interactions between oppositely charged amino acid side chains in proteins. They are often depicted as single conformer salt bridges (hydrogen-bonded ion pairs) in crystal structures, but it is unclear how dynamic they are in solution. Ion pairs are thought to be particularly important in stabilizing single α-helix (SAH) domains in solution. These highly stable domains are rich in charged residues (such as Arg, Lys, and Glu) with potential ion pairs across adjacent turns of the helix. They provide a good model system to investigate how ion pairs can contribute to protein stability. Using NMR spectroscopy, small-angle X-ray light scattering (SAXS), and molecular dynamics simulations, we provide here experimental evidence that ion pairs exist in a SAH in murine myosin 7a (residues 858-935), but that they are not fixed or long lasting. In silico modeling revealed that the ion pairs within this α-helix exhibit dynamic behavior, rapidly forming and breaking and alternating between different partner residues. The low-energy helical state was compatible with a great variety of ion pair combinations. Flexible ion pair formation utilizing a subset of those available at any one time avoided the entropic penalty of fixing side chain conformations, which likely contributed to helix stability overall. These results indicate the dynamic nature of ion pairs in SAHs. More broadly, thermodynamic stability in other proteins is likely to benefit from the dynamic behavior of multi-option solvent-exposed ion pairs.

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.

Glycosaminoglycan-Protein Interactions by Nuclear Magnetic Resonance (NMR) Spectroscopy

Nuclear magnetic resonance (NMR) spectroscopy is one of the most utilized and informative analytical techniques for investigating glycosaminoglycan (GAG)-protein complexes. NMR methods that are commonly applied to GAG-protein systems include chemical shift perturbation, saturation transfer difference, and transferred nuclear Overhauser effect. Although these NMR methods have revealed valuable insight into the protein-GAG complexes, elucidating high-resolution structural and dynamic information of these often transient interactions remains challenging. In addition, preparation of structurally homogeneous and isotopically enriched GAG ligands for structural investigations continues to be laborious. As a result, understanding of the structure-activity relationship of GAGs is still primitive. To overcome these deficiencies, several innovative NMR techniques have been developed lately. Here, we review some of the commonly used techniques along with more novel methods such as waterLOGSY and experiments to examine structure and dynamic of lysine and arginine side chains to identify GAG-binding sites. We will also present the latest technology that is used to produce isotopically enriched as well as paramagnetically tagged GAG ligands. Recent results that were obtained from solid-state NMR of amyloid’s interaction with GAG are also presented together with a brief discussion on computer assisted modeling of GAG-protein complexes using sparse experimental data.

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.

NMR-based investigations into target DNA search processes of proteins

To perform their function, transcription factors and DNA-repair/modifying enzymes must first locate their targets in the vast presence of nonspecific, but structurally similar sites on genomic DNA. Before reaching their targets, these proteins stochastically scan DNA and dynamically move from one site to another on DNA. Solution NMR spectroscopy provides unique atomic-level insights into the dynamic DNA-scanning processes, which are difficult to gain by any other experimental means. In this review, we provide an introductory overview on the NMR methods for the structural, dynamic, and kinetic investigations of target DNA search by proteins. We also discuss advantages and disadvantages of these NMR methods over other methods such as single-molecule techniques and biochemical approaches.

Impact of two-bond 15N-15N scalar couplings on 15N transverse relaxation measurements for arginine side chains of proteins

NMR relaxation of arginine (Arg) ¹⁵N_ε nuclei is useful for studying side-chain dynamics of proteins. In this work, we studied the impact of two geminal ¹⁵N-¹⁵N scalar couplings on measurements of transverse relaxation rates (R ₂ ) for Arg side-chain ¹⁵N_ε nuclei. For 12 Arg side chains of the DNA-binding domain of the Antp protein, we measured the geminal ¹⁵N-¹⁵N couplings ( ² J _NN ) of the ¹⁵N_ε nuclei and found that the magnitudes of the ² J _NN coupling constants were virtually uniform with an average of 1.2 Hz. Our simulations, assuming ideal 180° rotations for all ¹⁵N nuclei, suggested that the two ² J _NN couplings of this magnitude could in principle cause significant modulation in signal intensities during the Carr-Purcell-Meiboom-Gill (CPMG) scheme for Arg ¹⁵N_ε R ₂ measurements. However, our experimental data show that the expected modulation via two ² J _NN couplings vanishes during the ¹⁵N CPMG scheme. This quenching of J modulation can be explained by the mechanism described in Dittmer and Bodenhausen (Chemphyschem 7:831-836, 2006). This effect allows for accurate measurements of R ₂ relaxation rates for Arg side-chain ¹⁵N_ε nuclei despite the presence of two ² J _NN couplings. Although the so-called recoupling conditions may cause overestimate of R ₂ rates for very mobile Arg side chains, such conditions can readily be avoided through appropriate experimental settings.