Interdomain dynamics explored by paramagnetic NMR

Encoded loop-lanthanide-binding tags for long-range distance measurements in proteins by NMR and EPR spectroscopy

We recently engineered encodable lanthanide binding tags (LBTs) into proteins and demonstrated their applicability in Nuclear Magnetic Resonance (NMR) spectroscopy, X-ray crystallography and luminescence studies. Here, we engineered two-loop-LBTs into the model protein interleukin-1β (IL1β) and measured (1)H, (15)N-pseudocontact shifts (PCSs) by NMR spectroscopy. We determined the Δχ-tensors associated with each Tm(3+)-loaded loop-LBT and show that the experimental PCSs yield structural information at the interface between the two metal ion centers at atomic resolution. Such information is very valuable for the determination of the sites of interfaces in protein-protein-complexes. Combining the experimental PCSs of the two-loop-LBT construct IL1β-S2R2 and the respective single-loop-LBT constructs IL1β-S2, IL1β-R2 we additionally determined the distance between the metal ion centers. Further, we explore the use of two-loop LBTs loaded with Gd(3+) as a novel tool for distance determination by Electron Paramagnetic Resonance spectroscopy and show the NMR-derived distances to be remarkably consistent with distances derived from Pulsed Electron-Electron Dipolar Resonance.

Information content of long-range NMR data for the characterization of conformational heterogeneity

Long-range NMR data, namely residual dipolar couplings (RDCs) from external alignment and paramagnetic data, are becoming increasingly popular for the characterization of conformational heterogeneity of multidomain biomacromolecules and protein complexes. The question addressed here is how much information is contained in these averaged data. We have analyzed and compared the information content of conformationally averaged RDCs caused by steric alignment and of both RDCs and pseudocontact shifts caused by paramagnetic alignment, and found that, despite the substantial differences, they contain a similar amount of information. Furthermore, using several synthetic tests we find that both sets of data are equally good towards recovering the major state(s) in conformational distributions.

Identification of inhibitors of biological interactions involving intrinsically disordered proteins

Protein-protein interactions involving disordered partners have unique features and represent prominent targets in drug discovery processes. Intrinsically Disordered Proteins (IDPs) are involved in cellular regulation, signaling and control: they bind to multiple partners and these high-specificity/low-affinity interactions play crucial roles in many human diseases. Disordered regions, terminal tails and flexible linkers are particularly abundant in DNA-binding proteins and play crucial roles in the affinity and specificity of DNA recognizing processes. Protein complexes involving IDPs are short-lived and typically involve short amino acid stretches bearing few "hot spots", thus the identification of molecules able to modulate them can produce important lead compounds: in this scenario peptides and/or peptidomimetics, deriving from structure-based, combinatorial or protein dissection approaches, can play a key role as hit compounds. Here, we propose a panoramic review of the structural features of IDPs and how they regulate molecular recognition mechanisms focusing attention on recently reported drug-design strategies in the field of IDPs.

Visualizing transient dark states by NMR spectroscopy

Myriad biological processes proceed through states that defy characterization by conventional atomic-resolution structural biological methods. The invisibility of these 'dark' states can arise from their transient nature, low equilibrium population, large molecular weight, and/or heterogeneity. Although they are invisible, these dark states underlie a range of processes, acting as encounter complexes between proteins and as intermediates in protein folding and aggregation. New methods have made these states accessible to high-resolution analysis by nuclear magnetic resonance (NMR) spectroscopy, as long as the dark state is in dynamic equilibrium with an NMR-visible species. These methods - paramagnetic NMR, relaxation dispersion, saturation transfer, lifetime line broadening, and hydrogen exchange - allow the exploration of otherwise invisible states in exchange with a visible species over a range of timescales, each taking advantage of some unique property of the dark state to amplify its effect on a particular NMR observable. In this review, we introduce these methods and explore two specific techniques - paramagnetic relaxation enhancement and dark state exchange saturation transfer - in greater detail.

FANTEN: a new web-based interface for the analysis of magnetic anisotropy-induced NMR data

Pseudocontact shifts (PCSs) and residual dipolar couplings (RDCs) arising from the presence of paramagnetic metal ions in proteins as well as RDCs due to partial orientation induced by external orienting media are nowadays routinely measured as a part of the NMR characterization of biologically relevant systems. PCSs and RDCs are becoming more and more popular as restraints (1) to determine and/or refine protein structures in solution, (2) to monitor the extent of conformational heterogeneity in systems composed of rigid domains which can reorient with respect to one another, and (3) to obtain structural information in protein-protein complexes. The use of both PCSs and RDCs proceeds through the determination of the anisotropy tensors which are at the origin of these NMR observables. A new user-friendly web tool, called FANTEN (Finding ANisotropy TENsors), has been developed for the determination of the anisotropy tensors related to PCSs and RDCs and has been made freely available through the WeNMR ( http://fanten-enmr.cerm.unifi.it:8080 ) gateway. The program has many new features not available in other existing programs, among which the possibility of a joint analysis of several sets of PCS and RDC data and the possibility to perform rigid body minimizations.

Towards understanding the molecular recognition process in prokaryotic zinc-finger domain

Eukaryotic Cys2His2 zinc finger domain is one of the most common and important structural motifs involved in protein-DNA interaction. The recognition motif is characterized by the tetrahedral coordination of a zinc ion by conserved cysteine and histidine residues. We have characterized the prokaryotic Cys2His2 zinc finger motif, included in the DNA binding region (Ros87) of Ros protein from Agrobacterium tumefaciens, demonstrating that, although possessing a similar zinc coordination sphere, this domain presents significant differences from its eukaryotic counterpart. Furthermore, basic residues flanking the zinc binding region on either side have been demonstrated, by Electrophoretic Mobility Shift Assay (EMSA) experiments, to be essential for Ros DNA binding. In spite of this wealth of knowledge, the structural details of the mechanism through which the prokaryotic zinc fingers recognize their target genes are still unclear. Here, to gain insights into the molecular DNA recognition process of prokaryotic zinc finger domains we applied a strategy in which we performed molecular docking studies using a combination of Nuclear Magnetic Resonance (NMR) and Molecular Dynamics (MD) simulations data. The results demonstrate that the MD ensemble provides a reasonable picture of Ros87 backbone dynamics in solution. The Ros87-DNA model indicates that the interaction involves the first two residue of the first α-helix, and several residues located in the basic regions flanking the zinc finger domain. Interestingly, the prokaryotic zinc finger domain, mainly with the C-terminal tail that is wrapped around the DNA, binds a more extended recognition site than the eukaryotic counterpart. Our analysis demonstrates that the introduction of the protein flexibility in docking studies can improve, in terms of accuracy, the quality of the obtained models and could be particularly useful for protein showing high conformational heterogeneity as well as for computational drug design applications.

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

The present work demonstrates that NMR spin relaxation rate constants for molecules interconverting between states with different diffusion tensors can be modeled theoretically by combining orientational correlation functions for exchanging spherical molecules with locally isotropic approximations for the diffusion anisotropic tensors. The resulting expressions are validated by comparison with correlation functions obtained by Monte Carlo simulations and are accurate for moderate degrees of diffusion anisotropy typically encountered in investigations of globular proteins. The results are complementary to an elegant, but more complex, formalism that is accurate for all degrees of diffusion anisotropy [Ryabov, Y.; Clore, G. M.; Schwieters, C. D. J. Chem. Phys. 2012, 136, 034108].

Exploring regions of conformational space occupied by two-domain proteins

The presence of heterogeneity in the interdomain arrangement of several biomolecules is required for their function. Here we present a method to obtain crucial clues to distinguish between different kinds of protein conformational distributions based on experimental NMR data. The method explores subregions of the conformational space and provides both upper and lower bounds of probability for the system to be in each subregion.

Transient electrostatic interactions dominate the conformational equilibrium sampled by multidomain splicing factor U2AF65: a combined NMR and SAXS study

Multidomain proteins containing intrinsically disordered linkers exhibit large-scale dynamic modes that play key roles in a multitude of molecular recognition and signaling processes. Here, we determine the conformational space sampled by the multidomain splicing factor U2AF65 using complementary nuclear magnetic resonance spectroscopy and small-angle scattering data. Available degrees of conformational freedom are initially stochastically sampled and experimental data then used to delineate the potential energy landscape in terms of statistical probability. The spatial distribution of U2AF65 conformations is found to be highly anisotropic, comprising significantly populated interdomain contacts that appear to be electrostatic in origin. This hypothesis is supported by the reduction of signature PREs reporting on expected interfaces with increasing salt concentration. The described spatial distribution reveals the complete spectrum of the unbound forms of U2AF65 that coexist with the small percentage of a preformed RNA-bound domain arrangement required for polypyrimidine-tract recognition by conformational selection. More generally, the proposed approach to describing conformational equilibria of multidomain proteins can be further combined with other experimental data that are sensitive to domain dynamics.

Determination of the individual roles of the linker residues in the interdomain motions of calmodulin using NMR chemical shifts

Many protein molecules are formed by two or more domains whose structures and dynamics are closely related to their biological functions. It is thus important to develop methods to determine the structural properties of these multidomain proteins. Here, we characterize the interdomain motions in the calcium-bound state of calmodulin (Ca(2+)-CaM) using NMR chemical shifts as replica-averaged structural restraints in molecular dynamics simulations. We find that the conformational fluctuations of the interdomain linker, which are largely responsible for the overall interdomain motions of CaM, can be well described by exploiting the information provided by chemical shifts. We thus identify 10 residues in the interdomain linker region that change their conformations upon substrate binding. Five of these residues (Met76, Lys77, Thr79, Asp80 and Ser81) are highly flexible and cover the range of conformations observed in the substrate-bound state, while the remaining five (Arg74, Lys75, Asp78, Glu82 and Glu83) are much more rigid and do not populate conformations typical of the substrate-bound form. The ensemble of conformations representing the Ca(2+)-CaM state obtained in this study is in good agreement with residual dipolar coupling, paramagnetic resonance enhancement, small-angle X-ray scattering and fluorescence resonance energy transfer measurements, which were not used as restraints in the calculations. These results provide initial evidence that chemical shifts can be used to characterize the conformational fluctuations of multidomain proteins.

The IR-¹⁵N-HSQC-AP experiment: a new tool for NMR spectroscopy of paramagnetic molecules

A crucial factor for the understanding of structure-function relationships in metalloproteins is the identification of NMR signals from residues surrounding the metal cofactor. When the latter is paramagnetic, the NMR information in the proximity of the metal center may be scarce, because fast nuclear relaxation quenches signal intensity and coherence transfer efficiency. To identify residues at a short distance from a paramagnetic center, we developed a modified version of the ¹⁵N-HSQC experiment where (1) an inversion recovery filter is added prior to HSQC, (2) the INEPT period has been optimized according to fast relaxation of interested spins, (3) the inverse INEPT has been eliminated and signals acquired as antiphase doublets. The experiment has been successfully tested on a human [Fe₂S₂] protein which is involved in the biogenesis of iron-sulfur proteins. Thirteen HN resonances, unobserved with conventional HSQC experiments, could be identified. The structural arrangement of the protein scaffold in the proximity of the Fe/S cluster is fundamental to comprehend the molecular processes responsible for the transfer of Fe/S groups in the iron-sulfur protein assembly machineries.