Optimisation of spatially-encoded diffusion-ordered NMR spectroscopy for the analysis of mixtures

Diffusion-ordered NMR spectroscopy (DOSY NMR) is a widely used method for the analysis of mixtures. It can be used to separate the spectra of a mixture's components and to analyse interactions. The classic implementation of DOSY experiments, based on an incrementation of the diffusion-encoding gradient area, requires several minutes or more to collect a 2D data set. Spatially-encoded (SPEN) DOSY makes it possible to collect a complete data set in less than 1 s, by spatial parallelisation of the effective gradient area. While several short descriptions of SPEN DOSY experiments have been reported, a thorough characterisation of its features and its practical use is missing, and this hinders the use of the method. Here, we present the unusual principles and implementation of the SPEN DOSY experiment, an understanding of which is useful to make optimal use of the method. The encoding and acquisition steps are described, and the parameter relations that govern the setup of SPEN DOSY experiments are discussed. The influence of key parameters, including on sensitivity, is illustrated experimentally on mixtures of small molecules. This study should be useful for the setup of SPEN DOSY experiments, which are particularly useful for systems that evolve in time.

β-Hairpin Peptide Mimics Decrease Human Islet Amyloid Polypeptide (hIAPP) Aggregation

Amyloid diseases are degenerative pathologies, highly prevalent today because they are closely related to aging, that have in common the erroneous folding of intrinsically disordered proteins (IDPs) which aggregate and lead to cell death. Type 2 Diabetes involves a peptide called human islet amyloid polypeptide (hIAPP), which undergoes a conformational change, triggering the aggregation process leading to amyloid aggregates and fibers rich in β-sheets mainly found in the pancreas of all diabetic patients. Inhibiting the aggregation of amyloid proteins has emerged as a relevant therapeutic approach and we have recently developed the design of acyclic flexible hairpins based on peptidic recognition sequences of the amyloid β peptide (Aβ_1–42) as a successful strategy to inhibit its aggregation involved in Alzheimer’s disease. The present work reports the extension of our strategy to hIAPP aggregation inhibitors. The design, synthesis, conformational analyses, and biophysical evaluations of dynamic β-hairpin like structures built on a piperidine-pyrrolidine β-turn inducer are described. By linking to this β-turn inducer three different arms (i) pentapeptide, (ii) tripeptide, and (iii) α/aza/aza/pseudotripeptide, we demonstrate that the careful selection of the peptide-based arms from the sequence of hIAPP allowed to selectively modulate its aggregation, while the peptide character can be decreased. Biophysical assays combining, Thioflavin-T fluorescence, transmission electronic microscopy, capillary electrophoresis, and mass spectrometry showed that the designed compounds inhibit both the oligomerization and the fibrillization of hIAPP. They are also capable to decrease the aggregation process in the presence of membrane models and to strongly delay the membrane-leakage induced by hIAPP. More generally, this work provides the proof of concept that our rational design is a versatile and relevant strategy for developing efficient and selective inhibitors of aggregation of amyloidogenic proteins.

Structural and dynamic characterization of the C-terminal tail of ErbB2: Disordered but not random

ErbB2 (or HER2) is a receptor tyrosine kinase overexpressed in some breast cancers and associated with poor prognosis. Treatments targeting the receptor extracellular and kinase domains have greatly improved disease outcome in the last 20 years. In parallel, the structures of these domains have been described, enabling better mechanistic understanding of the receptor function and targeted inhibition. However, the ErbB2 disordered C-terminal cytoplasmic tail (CtErbB2) remains very poorly characterized in terms of structure, dynamics, and detailed functional mechanism. Yet, it is where signal transduction is triggered via phosphorylation of tyrosine residues and carried out via interaction with adaptor proteins. Here, we report the first description, to our knowledge, of the ErbB2 disordered tail at atomic resolution using NMR, complemented by small-angle x-ray scattering. We show that although no part of CtErbB2 has any fully populated secondary or tertiary structure, it contains several transient α-helices and numerous transient polyproline II helices, populated up to 20 and 40%, respectively, and low but significant compaction. The presence of some structural elements suggests, along the lines of the results obtained for EGFR (ErbB1), that they may have a functional role in ErbB2's autoregulation processes. In addition, the transient formation of polyproline II helices is compliant with previously suggested interactions with SH3 domains. All in all, our in-depth structural study opens perspectives in the mechanistic understanding of ErbB2.

[Formula: see text]H, [Formula: see text]C and [Formula: see text]N assignments of human Grb2 free of ligands

Growth factor receptor-bound 2 (Grb2) is an important link in the receptor tyrosine kinase signaling cascades. It is involved in crucial processes, both physiological (mainly embryogenesis) and pathological (different types of cancer). Several binding partners of all three domains (SH3-SH2-SH3) of this adaptor protein are well described, such as ErbB family members for the SH2 domain and Sos for the SH3 domains. How the different domains interact with each other, both structurally and functionally, is still unclear. These interactions could be essential for regulation processes, and therefore are of great interest. Although a lot of structural data on Grb2 exist, they describe either individual domains, ligand-bound conformations, or frozen pictures of the protein captured by crystallography. Here we report the assignment of backbone and of [Formula: see text] chemical shifts of full-length, apo-Grb2 in solution. In addition to the assigned conformation corresponding to three well-folded domains, a set of peaks compatible with the presence of an unfolded conformation of the N-terminal SH3 domain is observed. This assignment paves the way for future studies of inter-domain interactions and dynamics that have to be taken into account when studying the regulation of Grb2 interactions and signaling pathways.

NMR Characterization of the Influence of Zinc(II) Ions on the Structural and Dynamic Behavior of the New Delhi Metallo-β-Lactamase-1 and on the Binding with Flavonols as Inhibitors

New Delhi metallo-β-lactamase-1 (NDM-1) has recently emerged as a global threat because of its ability to confer resistance to all common β-lactam antibiotics. Understanding the molecular basis of β-lactam hydrolysis by NDM is crucial for designing NDM inhibitors or β-lactams resistant to their hydrolysis. In this study, for the first time, NMR was used to study the influence of Zn(II) ions on the dynamic behavior of NDM-1. Our results highlighted that the binding of Zn(II) in the NDM-1 active site induced several structural and dynamic changes on active site loop 2 (ASL2) and L9 loops and on helix α2. We subsequently studied the interaction of several flavonols: morin, quercetin, and myricetin were identified as natural and specific inhibitors of NDM-1. Quercetin conjugates were also synthesized in an attempt to increase the solubility and bioavailability. Our NMR investigations on NDM-1/flavonol interactions highlighted that both Zn(II) ions and the residues of the NDM-1 ASL1, ASL2, and ASL4 loops are involved in the binding of flavonols. This is the first NMR interaction study of NDM-1/inhibitors, and the models generated using HADDOCK will be useful for the rational design of more active inhibitors, directed against NDM-1.

Structural Characterization of N-WASP Domain V Using MD Simulations with NMR and SAXS Data

Because of their large conformational heterogeneity, structural characterization of intrinsically disordered proteins (IDPs) is very challenging using classical experimental methods alone. In this study, we use NMR and small-angle x-ray scattering (SAXS) data with multiple molecular dynamics (MD) simulations to describe the conformational ensemble of the fully disordered verprolin homology domain of the neural Aldrich syndrome protein involved in the regulation of actin polymerization. First, we studied several back-calculation software of SAXS scattering intensity and optimized the adjustable parameters to accurately calculate the SAXS intensity from an atomic structure. We also identified the most appropriate force fields for MD simulations of this IDP. Then, we analyzed four conformational ensembles of neural Aldrich syndrome protein verprolin homology domain, two generated with the program flexible-meccano with or without NMR-derived information as input and two others generated by MD simulations with two different force fields. These four conformational ensembles were compared to available NMR and SAXS data for validation. We found that MD simulations with the AMBER-03w force field and the TIP4P/2005s water model are able to correctly describe the conformational ensemble of this 67-residue IDP at both local and global level.

1H, 13C and 15N assignments of the C-terminal intrinsically disordered cytosolic fragment of the receptor tyrosine kinase ErbB2

ErbB2 (or HER2) is a receptor tyrosine kinase that is involved in signaling pathways controlling cell division, motility and apoptosis. Though important in development and cell growth homeostasis, this protein, when overexpressed, participates in triggering aggressive HER2+ breast cancers. It is composed of an extracellular part and a transmembrane domain, both important for activation by dimerization, and a cytosolic tyrosine kinase, which activates its intrinsically disordered C-terminal end (CtErbB2). Little is known about this C-terminal part of 268 residues, despite its crucial role in interacting with adaptor proteins involved in signaling. Understanding its structural and dynamic characteristics could eventually lead to the design of new interaction inhibitors, and treatments complementary to those already targeting other parts of ErbB2. Here we report backbone and side-chain assignment of CtErbB2, which, together with structural predictions, confirms its intrinsically disordered nature.

Spatially encoded diffusion-ordered NMR spectroscopy of reaction mixtures in organic solvents

Using a spatial encoding of the diffusion dimension, DOSY NMR data are acquired in less than one second for reaction mixtures in organic solvents.

Spatially encoded 2D and 3D diffusion-ordered NMR spectroscopy

The separation of 2D spectra of components in mixtures is accelerated with spatial encoding.

Efficient simulation of ultrafast magnetic resonance experiments

We present a convenient and powerful simulation formalism for ultrafast NMR spectroscopy. The formalism is based on the Fokker–Planck equation that supports systems with complicated combinations of classical spatial dynamics and quantum mechanical spin dynamics.

Mutations in actin used for structural studies partially disrupt β-thymosin/WH2 domains interaction

Understanding the structural basis of actin cytoskeleton remodeling requires stabilization of actin monomers, oligomers, and filaments in complex with partner proteins, using various biochemical strategies. Here, we report a dramatic destabilization of the dynamic interaction with a model β-thymosin/WH2 domain induced by mutations in actin. This result underlines that mutant actins should be used with prudence to characterize interactions with intrinsically disordered partners as destabilization of dynamic interactions, although identifiable by NMR, may be invisible to other structural techniques. It also highlights how both β-thymosin/WH2 domains and actin tune local structure and dynamics in regulatory processes involving intrinsically disordered domains.

Structure of a GPCR ligand in its receptor-bound state: leukotriene B4 adopts a highly constrained conformation when associated to human BLT2

G protein-coupled receptors (GPCRs) are key players in signal recognition and cell communication and are among the most important targets for drug development. Direct structural information on the conformation of GPCR ligands bound to their receptors is scarce. Using a leukotriene receptor, BLT2, expressed under a perdeuterated form in Escherichia coli , purified in milligram amounts, and folded to its native state using amphipols, we have solved, by (1)H NMR, the structure of receptor-bound leukotriene B4 (LTB4). Upon binding, LTB4 adopts a highly constrained seahorse conformation, at variance with the free state, where it explores a wide range of conformations. This structure provides an experimentally determined template of a pro-inflammatory compound for further pharmacological studies. The novel approach used for its determination could prove powerful to investigate ligand binding to GPCRs and membrane proteins in general.

Insight into the role of dynamics in the conformational switch of the small GTP-binding protein Arf1

Activation of the small GTP-binding protein Arf1, a major regulator of cellular traffic, follows an ordered sequence of structural events, which have been pictured by crystallographic snapshots. Combined with biochemical analysis, these data lead to a model of Arf1 activation, in which opening of its N-terminal helix first translocates Arf1-GDP to membranes, where it is then secured by a register shift of the interswitch β-strands, before GDP is eventually exchanged for GTP. However, how Arf1 rearranges its central β-sheet, an event that involves the loss and re-formation of H-bonds deep within the protein core, is not explained by available structural data. Here, we used Δ17Arf1, in which the N-terminal helix has been deleted, to address this issue by NMR structural and dynamics analysis. We first completed the assignment of Δ17Arf1 bound to GDP, GTP, and GTPγS and established that NMR data are fully consistent with the crystal structures of Arf1-GDP and Δ17Arf1-GTP. Our assignments allowed us to analyze the kinetics of both protein conformational transitions and nucleotide exchange by real-time NMR. Analysis of the dynamics over a very large range of timescale by (15)N relaxation, CPMG relaxation dispersion and H/D exchange reveals that while Δ17Arf1-GTP and full-length Arf1-GDP dynamics is restricted to localized fast motions, Δ17Arf1-GDP features unique intermediate and slow motions in the interswitch region. Altogether, the NMR data bring insight into how that membrane-bound Arf1-GDP, which is mimicked by the truncation of the N-terminal helix, acquires internal motions that enable the toggle of the interswitch.

Robust structure-based resonance assignment for functional protein studies by NMR

High-throughput functional protein NMR studies, like protein interactions or dynamics, require an automated approach for the assignment of the protein backbone. With the availability of a growing number of protein 3D structures, a new class of automated approaches, called structure-based assignment, has been developed quite recently. Structure-based approaches use primarily NMR input data that are not based on J-coupling and for which connections between residues are not limited by through bonds magnetization transfer efficiency. We present here a robust structure-based assignment approach using mainly H(N)-H(N) NOEs networks, as well as (1)H-(15) N residual dipolar couplings and chemical shifts. The NOEnet complete search algorithm is robust against assignment errors, even for sparse input data. Instead of a unique and partly erroneous assignment solution, an optimal assignment ensemble with an accuracy equal or near to 100% is given by NOEnet. We show that even low precision assignment ensembles give enough information for functional studies, like modeling of protein-complexes. Finally, the combination of NOEnet with a low number of ambiguous J-coupling sequential connectivities yields a high precision assignment ensemble. NOEnet will be available under: http://www.icsn.cnrs-gif.fr/download/nmr.

Solution NMR mapping of water-accessible residues in the transmembrane beta-barrel of OmpX

The atomic structure of OmpX, the smallest member of the bacterial outer membrane protein family, has been previously established by X-ray crystallography and NMR spectroscopy. In apparent conflict with electrophysiological studies, the lumen of its transmembrane beta-barrel appears too tightly packed with amino acid side chains to let any solute flow through. In the present study, high-resolution solution NMR spectra were obtained of OmpX kept water-soluble by either amphipol A8-35 or the detergent dihexanoylphosphatidylcholine. Hydrogen/deuterium exchange measurements performed after prolonged equilibration show that, whatever the surfactant used, some of the amide protons of the membrane-spanning region exchange much more readily than others, which likely reflects the dynamics of the barrel.

Inter- and intramolecular contacts in a membrane protein/surfactant complex observed by heteronuclear dipole-to-dipole cross-relaxation

Heteronuclear dipole-to-dipole cross-relaxation has been applied to exploring intermolecular interactions and intramolecular spatial proximities in a large supramolecular structure comprised of a beta-barrel membrane protein, OmpX, in complex with a polymeric surfactant, amphipol A8-35. The experiments, performed in either the laboratory or the rotating frame, reveal the existence of intermolecular contacts between aromatic amino acids and specific groups of the polymer, in addition to intra-protein dipolar interactions, some of them involving carbonyl carbons. This study opens the perspective of collecting by NMR spectroscopy a new kind of through-space structural information involving aromatic and carbonyl (13)C atoms of large proteins.

NOEnet--use of NOE networks for NMR resonance assignment of proteins with known 3D structure

Motivation: A prerequisite for any protein study by NMR is the assignment of the resonances from the ¹⁵N−¹H HSQC spectrum to their corresponding atoms of the protein backbone. Usually, this assignment is obtained by analyzing triple resonance NMR experiments. An alternative assignment strategy exploits the information given by an already available 3D structure of the same or a homologous protein. Up to now, the algorithms that have been developed around the structure-based assignment strategy have the important drawbacks that they cannot guarantee a high assignment accuracy near to 100%.

Results: We propose here a new program, called NOEnet, implementing an efficient complete search algorithm that ensures the correctness of the assignment results. NOEnet exploits the network character of unambiguous NOE constraints to realize an exhaustive search of all matching possibilities of the NOE network onto the structural one. NOEnet has been successfully tested on EIN, a large protein of 28 kDa, using only NOE data. The complete search of NOEnet finds all possible assignments compatible with experimental data that can be defined as an assignment ensemble. We show that multiple assignment possibilities of large NOE networks are restricted to a small spatial assignment range (SAR), so that assignment ensembles, obtained from accessible experimental data, are precise enough to be used for functional proteins studies, like protein–ligand interaction or protein dynamics studies. We believe that NOEnet can become a major tool for the structure-based backbone resonance assignment strategy in NMR.

Availability: The NOEnet program will be available under: http://www.icsn.cnrs-gif.fr/download/nmr

Contact:carine@icsn.cnrs-gif.fr; eric.guittet@icsn.cnrs-gif.fr

Supplementary Information:Supplementary data are available at Bioinformatics online.

(1)H, (13)C and (15)N resonance assignment of the first N-terminal RNA recognition motif (RRM) of the human heterogeneous nuclear ribonucleoprotein H (hnRNP H)

Human heterogeneous nuclear ribonucleoprotein H (hnRNP H) regulates alternative splicing of HIV-1 Tat pre-mRNA. The structure of the first N-terminal domain (residues 1-104) of hnRNP H was solved and its binding to an exonic splicing silencer (pESS2) studied. For this, all backbone and 85% of side-chain resonance frequencies were assigned.

Structure of the histone chaperone ASF1 bound to the histone H3 C-terminal helix and functional insights

Asf1 is a histone chaperone that favors histone H3/H4 assembly and disassembly. We solved the structure of the conserved domain of human ASF1A in complex with the C-terminal helix of histone H3 using nuclear magnetic resonance spectroscopy. This structure is fully compatible with an association of ASF1 with the heterodimeric form of histones H3/H4. In our model, ASF1 substitutes for the second H3/H4 heterodimer that is normally found in heterotetrameric H3/H4 complexes. This result constitutes an essential step in the fundamental understanding of the mechanisms of nucleosome assembly by histone chaperones. Point mutations that perturb the Asf1/histone interface were designed from the structure. The decreased binding affinity of the Asf1-H3/H4 complex correlates with decreased levels of H3-K56 acetylation and phenotypic defects in vivo.

Coupling of Folding and Binding of Thymosin β4 upon Interaction with Monomeric Actin Monitored by Nuclear Magnetic Resonance

Thymosin beta4 is a major actin-sequestering protein, yet the structural basis for its biological function is still unknown. This study provides insight regarding the way this 43-amino acid peptide, mostly unstructured in solution, binds to monomeric actin and prevents its assembly in filaments. We show here that the whole backbone of thymosin beta4 is highly affected upon binding to G-actin. The assignment of all amide protons and nitrogens of thymosin in the bound state, obtained using a combination of NMR experiments and selective labelings, shows that thymosin folds completely upon binding and displays a central extended region flanked by two N- and C-terminal helices. The cleavage of actin by subtilisin in the DNase I binding loop does not modify the structure of thymosin beta4 in the complex, showing that the backbone of the peptide is not in close proximity to segment 42-47 of actin. The combination of our NMR results and previously published mutation and cross-link data allows a better characterization of the binding mode of thymosins on G-actin.

The β-Thymosin/WH2 Domain

The widespread beta-thymosin/WH2 actin binding domain has versatile regulatory properties in actin dynamics and motility. beta-thymosins (isolated WH2 domain) maintain monomeric actin in a "sequestered" nonpolymerizable form. In contrast, when repeated in tandem or inserted in modular proteins, the beta-thymosin/WH2 domain promotes actin assembly at filament barbed ends, like profilin. The structural basis for these opposite functions is addressed using ciboulot, a three beta-thymosin repeat protein. Only the first repeat binds actin and possesses the function of ciboulot. The region that shows the strongest interaction with actin is an amphipathic N-terminal alpha helix, present in all beta-thymosin/WH2 domains, which recognizes the ATP bound actin structure and uses the shear motion of actin linked to ATP hydrolysis to control polymerization. Crystallographic ((1)H, (15)N), NMR, and mutagenetic data reveal that the weaker interaction of the C-terminal region of beta-thymosin/WH2 domain with actin accounts for the switch in function from inhibition to promotion of actin assembly.

Dynamical characterization of residual and non-native structures in a partially folded protein by (15)N NMR relaxation using a model based on a distribution of correlation times

A spectral density model based on a truncated lorentzian distribution of correlation times is used to analyze the nanosecond time-scale dynamics of the partially unfolded domain 2 of annexin I from its (15)N NMR relaxation parameters measured at three magnetic field strengths. The use of a distribution of correlation times enables the characterization of the dynamical features of the NH bonds of the protein in terms of heterogeneity of dynamical states in the nanosecond range. The variation along the sequence of the two dynamical parameters introduced, namely the center and the width of the distribution, points out the different types of residual secondary structures present in the D2 domain. Moreover, it allows a physically sensible interpretation of the dynamical behavior of the different residual helices and of the non-native structures. Also, a striking correspondence is found between the parameters obtained using an extended Lipari and Szabo model and the parameters obtained using the distribution of correlation times. This result led us to propose a specific interpretation of the model-free order parameter for internal motions in the nanosecond range in the case of unfolded states.