Molecular Insight into Lipid Nanoparticle Assembly from NMR Spectroscopy and Molecular Dynamics Simulation

Lipid nanoparticles (LNPs) have emerged as the premier drug delivery system for oligonucleotide vaccines and therapeutics in recent years. Despite their prosperous advancement in research and clinical applications, there is a significant lack of mechanistic understanding of the assembly of lipid particles at the molecular level. In our study, we utilized a combination of solution and solid-state NMR, together with molecular dynamics simulations, to elucidate local structures and interactions of chemical components across multiple motional regimes. Our results comprehensively evaluated the impact of formulation components and engineering process factors on the particle formation and identified the interplay of phospholipids (DSPC), poly(ethylene glycol) (PEG) lipid conjugates, and cholesterol in governing the particle size and lipid dynamics from a structural perspective, using static 31P NMR techniques. These studies provide novel insights into the impact of particle engineering on the molecular properties of the LNP envelope membrane. Additionally, molecular interactions and compositional distribution play a critical role in particle engineering and the consequent stability and potency. In this study, we have identified intermolecular contacts among the lipid components using one-dimensional 1H-13C cross-polarization magic angle spinning experiments, 1H relaxation measurements, and two-dimensional 1H-1H correlation methods, providing a structural basis for the lipid assembly. Interestingly, the cationic and ionizable lipids, conventionally regarded as stabilizing agents primarily located within the core of LNPs, were found to interact with PEG lipids and coexist in the outer layer of the particles. We suggest that LNPs examined here are comprised of an outer layer rich in lipid components surrounding a core region. Our high-resolution findings offer insightful structural and dynamic details pertaining to the individual chemical components in the lipid particles and their interactions influence lipid complex structure and stability in particle engineering.

Pitfalls in measurements of R1 relaxation rates of protein backbone 15N nuclei

The dynamics of the backbone and side-chains of protein are routinely studied by interpreting experimentally determined 15N spin relaxation rates. R1(15N), the longitudinal relaxation rate, reports on fast motions and encodes, together with the transverse relaxation R2, structural information about the shape of the molecule and the orientation of the amide bond vectors in the internal diffusion frame. Determining error-free 15N longitudinal relaxation rates remains a challenge for small, disordered, and medium-sized proteins. Here, we show that mono-exponential fitting is sufficient, with no statistical preference for bi-exponential fitting up to 800 MHz. A detailed comparison of the TROSY and HSQC techniques at medium and high fields showed no statistically significant differences. The least error-prone DD/CSA interference removal technique is the selective inversion of amide signals while avoiding water resonance. The exchange of amide with solvent deuterons appears to affect the rate R1 of solvent-exposed amides in all fields tested and in each DD/CSA interference removal technique in a statistically significant manner. In summary, the most accurate R1(15N) rates in proteins are achieved by selective amide inversion, without the addition of D2O. Importantly, at high magnetic fields stronger than 800 MHz, when non-mono-exponential decay is involved, it is advisable to consider elimination of the shortest delays (typically up to 0.32 s) or bi-exponential fitting.

Nanosecond Molecular Motion in pHP1α Liquid-Liquid Phase Separation Captured by Solid-State NMR

The relationship among protein structure, function, and dynamics is fundamental to biological activity, particularly in more complex biomolecular systems. Solid-state and solution-state NMR techniques offer powerful means to probe these dynamics across various time scales. However, standard assumptions about molecular motion are often challenged in phase-separated systems like phosphorylated heterochromatin protein 1 alpha (pHP1α), which exhibit both solid- and solution-like characteristics. This study investigates the nanosecond molecular motions in pHP1α liquid-liquid phase separation (LLPS) using relaxation in hetNOE-filtered HSQC signals. By systematically analyzing motions captured by hetNOE-filtered HSQC and conventional HSQC, we characterize the global dynamics site-specifically in pHP1α LLPS. Our findings reveal ∼15 ns motion in the pHP1α LLPS system, suggesting the coexistence of different dynamic phases, and support previous observations on its role in chromatin organization. This work contributes to the expanding literature on phase-separated biomolecular behavior, with implications for understanding the molecular basis of chromatin compaction and genomic stability.

Molecular dynamics as an efficient process to predict 15N chemical shift anisotropy at very high NMR magnetic fields

The emergence of very high NMR magnetic fields will certainly encourage the study of larger biological systems with their dynamics and interactions. NMR spin relaxation allows probing the dynamical properties of proteins where the 15N longitudinal (R1) and transverse (R2) relaxation rates in addition to the 1H-15N heteronuclear NOE describe the ps-ns time scale. Their analytical representation involves the chemical shift anisotropy (CSA) effect that represents the major contribution at a very high magnetic field above 18.8 T. An accurate analysis of the latter parameters in terms of model free (MF) requires considering its effect. Until now, a uniform value of -160 ppm for the CSA has been widely used to derive the backbone order parameters (S2), giving rise to a large fluctuation of its value at very high magnetic fields. Conversely, the use of a site-specific CSA improves the accurate analysis of protein dynamics but requires a cost-effective experimental multi-field approach. In the present paper, we show how the CSA mainly contributes to the relaxation parameters at 28.2 T compared to lower magnetic fields and may bias the determination of S2. We propose to replace the time-consuming measurement of spin relaxation at multiple fields by a combination of molecular dynamics (MD) and the measurement of spin relaxation at one very high magnetic field only. We applied this strategy to three well-folded proteins (ubiquitin, GB3 and ribonuclease H) to show that the determined order parameters are in good agreement with the ones obtained by means of experimental data only.

Structural dynamics of calcium and integrin-binding protein 2 (CIB2) reveal uncommon flexibility and heterogeneous calcium and magnesium loading

Calcium- and Integrin-Binding protein 2 (CIB2) is a widely expressed protein with an uncertain biological role. Two of its four EF-hand motifs bind Mg(II) and/or Ca(II), thus triggering conformational changes. Although previous studies suggested that CIB2 preferentially binds Mg(II) over Ca(II) under physiological conditions, an atomic level characterization of CIB2 in the presence of both cations was lacking. Based on a combination of solution NMR, exhaustive molecular dynamics simulations and isothermal titration and differential scanning calorimetry, we characterized the interaction of CIB2 with both Ca(II) and Mg(II) ions and elucidated the protein regions involved in the interaction with the α7B integrin target. Analysis of experimental amide nitrogen relaxation rates shows that the EF4 motif exhibits high mobility regardless of the specific bound metal ion and demonstrates that the Mg(II)- and Ca(II)-bound state of CIB2 is relatively floppy, with pico-nanosecond motions induced in a region involved in target recognition. Overall, our data indicate a preferential, thermodynamically stable but structurally flexible state for CIB2, in which a Mg(II) ion is bound to EF3 and a Ca(II) ion to EF4. These results unveil the role of metal binding events in CIB2 and offer new insights into the dynamic regulation of target recognition.

Structural basis for regulation of a CBASS-CRISPR-Cas defense island by a transmembrane anti-σ factor and its ECF σ partner

How CRISPR-Cas and cyclic oligonucleotide-based antiphage signaling systems (CBASS) are coordinately deployed against invaders remains unclear. We show that a locus containing two CBASS and one type III-B CRISPR-Cas system, regulated by the transmembrane anti-σ DdvA and its cognate extracytoplasmic function (ECF) σ DdvS, can defend Myxococcus xanthus against a phage. Cryo-electron microscopy reveals DdvA-DdvS pairs assemble as arrow-shaped transmembrane dimers. Each DdvA periplasmic domain adopts a separase/craspase-type tetratricopeptide repeat (TPR)-caspase HetF-associated with TPR (TPR-CHAT) architecture with an incomplete His-Cys active site, lacking three α-helices conserved among CHAT domains. Each active site faces the dimer interface, raising the possibility that signal-induced caspase-like DdvA autoproteolysis in trans precedes RseP-mediated intramembrane proteolysis and DdvS release. Nuclear magnetic resonance reveals a DdvA cytoplasmic CHCC-type zinc-bound three-helix bundle that binds to DdvS σ2 and σ4 domains, undergoing σ4-induced helix extension to trap DdvS. Altogether, we provide structural-mechanistic insights into membrane anti-σ-ECF σ regulation of an antiviral CBASS-CRISPR-Cas defense island.

AsIII Selectively Induces a Disorder-to-Order Transition in the Metalloid Binding Region of the AfArsR Protein

Arsenic is highly toxic and a significant threat to human health, but certain bacteria have developed defense mechanisms initiated by AsIII binding to AsIII-sensing proteins of the ArsR family. The transcriptional regulator AfArsR responds to AsIII and SbIII by coordinating the metalloids with three cysteines, located in a short sequence of the same monomer chain. Here, we characterize the binding of AsIII and HgII to a model peptide encompassing this fragment of the protein via solution equilibrium and spectroscopic/spectrometric techniques (pH potentiometry, UV, CD, NMR, PAC, EXAFS, and ESI-MS) combined with DFT calculations and MD simulations. Coordination of AsIII changes the peptide structure from a random-coil to a well-defined structure of the complex. A trigonal pyramidal AsS3 binding site is formed with almost exactly the same structure as observed in the crystal structure of the native protein, implying that the peptide possesses all of the features required to mimic the AsIII recognition and response selectivity of AfArsR. Contrary to this, binding of HgII to the peptide does not lead to a well-defined structure of the peptide, and the atoms near the metal binding site are displaced and reoriented in the HgII model. Our model study suggests that structural organization of the metal site by the inducer ion is a key element in the mechanism of the metalloid-selective recognition of this protein.

Unraveling motion in proteins by combining NMR relaxometry and molecular dynamics simulations: A case study on ubiquitin

Nuclear magnetic resonance (NMR) relaxation experiments shine light onto the dynamics of molecular systems in the picosecond to millisecond timescales. As these methods cannot provide an atomically resolved view of the motion of atoms, functional groups, or domains giving rise to such signals, relaxation techniques have been combined with molecular dynamics (MD) simulations to obtain mechanistic descriptions and gain insights into the functional role of side chain or domain motion. In this work, we present a comparison of five computational methods that permit the joint analysis of MD simulations and NMR relaxation experiments. We discuss their relative strengths and areas of applicability and demonstrate how they may be utilized to interpret the dynamics in MD simulations with the small protein ubiquitin as a test system. We focus on the aliphatic side chains given the rigidity of the backbone of this protein. We find encouraging agreement between experiment, Markov state models built in the χ1/χ2 rotamer space of isoleucine residues, explicit rotamer jump models, and a decomposition of the motion using ROMANCE. These methods allow us to ascribe the dynamics to specific rotamer jumps. Simulations with eight different combinations of force field and water model highlight how the different metrics may be employed to pinpoint force field deficiencies. Furthermore, the presented comparison offers a perspective on the utility of NMR relaxation to serve as validation data for the prediction of kinetics by state-of-the-art biomolecular force fields.

From Milliseconds to Minutes: Melittin Self-Assembly from Concerted Non-Equilibrium Pressure-Jump and Equilibrium Relaxation Nuclear Magnetic Resonance

Non-equilibrium kinetics techniques like pressure-jump nuclear magnetic resonance (NMR) are powerful in tracking changes in oligomeric populations and are not limited by relaxation rates for the time scales of exchange that can be probed. However, these techniques are less sensitive to minor, transient populations than are Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments. We integrated non-equilibrium pressure-jump and equilibrium CPMG relaxation dispersion data to fully map the kinetic landscape of melittin tetramerization. While monomeric peptides weakly form dimers (Kd,D/M ≈ 26 mM) whose population never exceeds 1.6% at 288 K, dimers associate tightly to form stable tetrameric species (Kd,T/D ≈ 740 nM). Exchange between the monomer and dimer, along with exchange between the dimer and tetramer, occurs on the millisecond time scale. The NMR approach developed herein can be readily applied to studying the folding and misfolding of a wide range of oligomeric assemblies.

Comprehensive analysis of relaxation decays from high-resolution relaxometry

Relaxometry consists in measuring relaxation rates over orders of magnitude of magnetic fields to probe motions of complex systems. High-resolution relaxometry (HRR) experiments can be performed on conventional high-field NMR magnets equipped with a sample shuttle. During the experiment, the sample shuttle transfers the sample between the high-field magnetic center and a chosen position in the stray field for relaxation during a variable delay, thus using the stray field as a variable field. As the relaxation delay occurs outside of the probe, HRR experiments cannot rely on the control of cross-relaxation pathways, which is standard in high-field relaxation pulse sequences. Thus, decay rates are not pure relaxation rates, which may impair a reliable description of the dynamics. Previously, we took into account cross-relaxation effects in the analysis of high-resolution relaxometry data by applying a correction factor to relaxometry decay rates in order to estimate relaxation rates. These correction factors were obtained from the iterative simulation of the relaxation decay while the sample lies outside of the probe and a preceding analysis of relaxation rates which relies on the approximation of a priori multi-exponential decays by mono-exponential functions. However, an analysis protocol matching directly experimental and simulated relaxometry decays should be more self consistent and more generally applicable as it can accommodate deviations from mono-exponential decays. Here, we introduce Matching INtensities for the Optimization of Timescales and Amplitudes of motions Under Relaxometry (MINOTAUR), a framework for the analysis of high-resolution relaxometry that takes as input the intensity decays at all fields. This approach uses the full relaxation matrix to calculate intensity decays, allowing complex relaxation pathways to be taken into account. Therefore, it eliminates the need for a correction of decay rates and for fitting multi-exponential decays with mono-exponential functions. The MINOTAUR software is designed as a flexible framework where relaxation matrices and spectral density functions corresponding to various models of motions can be defined on a case-by-case basis. The agreement with our previous analyses of protein side-chain dynamics from carbon-13 relaxation is excellent, while providing a more robust analysis tool. We expect MINOTAUR to become the tool of choice for the analysis of high-resolution relaxometry.

Studies of proline conformational dynamics in IDPs by 13C-detected cross-correlated NMR relaxation

Intrinsically disordered proteins (IDPs) are significantly enriched in proline residues, which can populate specific local secondary structural elements called PPII helices, characterized by small packing densities. Proline is often thought to promote disorder, but it can participate in specific π·CH interactions with aromatic side chains resulting in reduced conformational flexibilities of the polypeptide. Differential local motional dynamics are relevant for the stabilization of preformed structural elements and can serve as nucleation sites for the establishment of long-range interactions. NMR experiments to probe the dynamics of proline ring systems would thus be highly desirable. Here we present a pulse scheme based on 13C detection to quantify dipole-dipole cross-correlated relaxation (CCR) rates at methylene CH2 groups in proline residues. Applying 13C-CON detection strategy provides exquisite spectral resolution allowing applications also to high molecular weight IDPs even in conditions approaching the physiological ones. The pulse scheme is illustrated with an application to the 220 amino acids long protein Osteopontin, an extracellular cytokine involved in inflammation and cancer progression, and a construct in which three proline-aromatic sequence patches have been mutated.

Fast Motions Dominate Dynamics of Intrinsically Disordered Tau Protein at High Temperatures

Reorientational dynamics of intrinsically disordered proteins (IDPs) contain multiple motions often clustered around three motional modes: ultrafast librational motions of amide groups, fast local backbone conformational fluctuations and slow chain segmental motions. This dynamic picture is mainly based on ¹⁵ N NMR relaxation studies of IDPs at relatively low temperatures where the amide-water proton exchange rates are sufficiently small. Less is known, however, about the dynamics of IDPs at more physiological temperatures. Here, we investigate protein dynamics in a 441-residue long IDP, tau protein, in the temperature range from 0-25 °C, using ¹⁵ N NMR relaxation rates and spectral density analysis. While at these temperatures relaxation rates are still better described in terms of amide group librational motions, local backbone dynamics and chain segmental motions, the temperature-dependent trend of spectral densities suggests that the timescales of fast backbone conformational fluctuations and slower chain segmental motions might become inseparable at higher temperatures. Our data demonstrate the remarkable dynamic plasticity of this prototypical IDP and highlight the need for dynamic studies of IDPs at multiple temperatures.

A Comparison of Bonded and Nonbonded Zinc(II) Force Fields with NMR Data

Classical molecular dynamics (MD) simulations are widely used to inspect the behavior of zinc(II)-proteins at the atomic level, hence the need to properly model the zinc(II) ion and the interaction with its ligands. Different approaches have been developed to represent zinc(II) sites, with the bonded and nonbonded models being the most used. In the present work, we tested the well-known zinc AMBER force field (ZAFF) and a recently developed nonbonded force field (NBFF) to assess how accurately they reproduce the dynamic behavior of zinc(II)-proteins. For this, we selected as benchmark six zinc-fingers. This superfamily is extremely heterogenous in terms of architecture, binding mode, function, and reactivity. From repeated MD simulations, we computed the order parameter (S2) of all backbone N-H bond vectors in each system. These data were superimposed to heteronuclear Overhauser effect measurements taken by NMR spectroscopy. This provides a quantitative estimate of the accuracy of the FFs in reproducing protein dynamics, leveraging the information about the protein backbone mobility contained in the NMR data. The correlation between the MD-computed S2 and the experimental data indicated that both tested FFs reproduce well the dynamic behavior of zinc(II)-proteins, with comparable accuracy. Thus, along with ZAFF, NBFF represents a useful tool to simulate metalloproteins with the advantage of being extensible to diverse systems such as those bearing dinuclear metal sites.

How does it really move? Recent progress in the investigation of protein nanosecond dynamics by NMR and simulation

Nuclear magnetic resonance (NMR) spin relaxation experiments currently probe molecular motions on timescales from picoseconds to nanoseconds. The detailed interpretation of these motions in atomic detail benefits from complementarity with the results from molecular dynamics (MD) simulations. In this mini-review, we describe the recent developments in experimental techniques to study the backbone dynamics from ¹⁵N relaxation and side-chain dynamics from ¹³C relaxation, discuss the different analysis approaches from model-free to dynamics detectors, and highlight the many ways that NMR relaxation experiments and MD simulations can be used together to improve the interpretation and gain insights into protein dynamics.

Explicit Models of Motion to Understand Protein Side-Chain Dynamics

Nuclear magnetic relaxation is widely used to probe protein dynamics. For decades, most analyses of relaxation in proteins have relied successfully on the model-free approach, forgoing mechanistic descriptions of motion. Model-free types of correlation functions cannot describe a large carbon-13 relaxation dataset in protein side chains. Here, we use molecular dynamics simulations to design explicit models of motion and solve Fokker-Planck diffusion equations. These models of motion provide better agreement with relaxation data, mechanistic insight, and a direct link to configuration entropy.

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.

Leveraging the multivalent p53 peptide-MdmX interaction to guide the improvement of small molecule inhibitors

Overexpressed Mdm2 and its 7homolog MdmX impair p53 activity in many cancers. Small molecules mimicking a p53 peptide can effectively inhibit Mdm2 but not MdmX. Here, we show a strategy for improving lead compounds for Mdm2 and MdmX inhibition based on the multivalency of the p53 peptide. Crystal structures of MdmX complexed with nutlin-3a, a strong Mdm2 inhibitor but a weak one for MdmX, reveal that nutlin-3a fits into the ligand binding pocket of MdmX mimicking the p53 peptide. However, due to distinct flexibility around the MdmX ligand binding pocket, the structures are missing many important intermolecular interactions that exist in the MdmX/p53 peptide and Mdm2/nultin-3a complexes. By targeting these flexible regions, we identify allosteric and additive fragments that enhance the binding affinity of nutlin-3a for MdmX, leading to potent Mdm2/MdmX inhibitors with anticancer activity. Our work provides a practical approach to drug design for signal transduction therapy.

Peptide fragments derived from the interfaces of protein-protein interactions (PPIs) provide useful templates for designing small molecule PPI inhibitors. Here, the authors utilize the multivalency of an MdmX-binding p53 peptide to develop a weak inhibitor of MdmX into potent Mdm2/MdmX inhibitors.

The long-standing relationship between paramagnetic NMR and iron-sulfur proteins: the mitoNEET example. An old method for new stories or the other way around?

Paramagnetic NMR spectroscopy and iron-sulfur (Fe-S) proteins have maintained a synergic relationship for decades. Indeed, the hyperfine shifts with their temperature dependencies and the relaxation rates of nuclei of cluster-bound residues have been extensively used as a fingerprint of the type and of the oxidation state of the Fe-S cluster within the protein frame. The identification of NMR signals from residues surrounding the metal cofactor is crucial for understanding the structure-function relationship in Fe-S proteins, but it is generally impaired in standard NMR experiments by paramagnetic relaxation enhancement due to the presence of the paramagnetic cluster(s). On the other hand, the availability of systems of different sizes and stabilities has, over the years, stimulated NMR spectroscopists to exploit iron-sulfur proteins as paradigmatic cases to develop experiments, models, and protocols. Here, the cluster-binding properties of human mitoNEET have been investigated by 1D and 2D 1 H diamagnetic and paramagnetic NMR, in its oxidized and reduced states. The NMR spectra of both oxidation states of mitoNEET appeared to be significantly different from those reported for previously investigated [ Fe 2 S 2 ] 2 + / + proteins. The protocol we have developed in this work conjugates spectroscopic information arising from "classical" paramagnetic NMR with an extended mapping of the signals of residues around the cluster which can be taken, even before the sequence-specific assignment is accomplished, as a fingerprint of the protein region constituting the functional site of the protein. We show how the combined use of 1D NOE experiments, 13 C direct-detected experiments, and double- and triple-resonance experiments tailored using R1 - and/or R2 -based filters significantly reduces the "blind" sphere of the protein around the paramagnetic cluster. This approach provided a detailed description of the unique electronic properties of mitoNEET, which are responsible for its biological function. Indeed, the NMR properties suggested that the specific electronic structure of the cluster possibly drives the functional properties of different [ Fe 2 S 2 ] proteins.

The armadillo-repeat domain of plakophilin 1 binds the C-terminal sterile alpha motif (SAM) of p73

BACKGROUND

Plakophilin 1 (PKP1) is a component of desmosomes, which are key structural components for cell-cell adhesion, and can also be found in other cell locations. The p53, p63 and p73 proteins belong to the p53 family of transcription factors, playing crucial roles in tumour suppression. The α-splice variant of p73 (p73α) has at its C terminus a sterile alpha motif (SAM); such domain, SAMp73, is involved in the interaction with other macromolecules.

METHODS

We studied the binding of SAMp73 with the armadillo domain of PKP1 (ARM-PKP1) in the absence and the presence of 100 mM NaCl, by using several biophysical techniques, namely fluorescence, far-ultraviolet circular dichroism (CD), nuclear magnetic resonance (NMR), isothermal titration calorimetry (ITC), and molecular docking and simulations.

RESULTS

Association was observed between the two proteins, with a dissociation constant of ~5 μM measured by ITC and fluorescence in the absence of NaCl. The binding region of SAMp73 involved residues of the so-called "middle-loop-end-helix" binding region (i.e., comprising the third helix, together with the C terminus of the second one, and the N-cap of the fourth), as shown by ¹⁵N, ¹H- HSQC-NMR spectra. Molecular modelling provided additional information on the possible structure of the binding complex.

CONCLUSIONS

This newly-observed interaction could have potential therapeutic relevance in the tumour pathways where PKP1 is involved, and under conditions when there is a possible inactivation of p53.

GENERAL SIGNIFICANCE

The discovery of the binding between SAMp73 and ARM-PKP1 suggests a functional role for their interaction, including the possibility that SAMp73 could assist PKP1 in signalling pathways.

Investigation of protein-ligand complexes by ligand-based NMR methods

Molecular recognition is at the base of all biological events and its knowledge at atomic level is pivotal in the development of new drug design approaches. NMR spectroscopy is one of the most widely used technique to detect and characterize transient ligand-receptor interactions in solution. In particular, ligand-based NMR approaches, including NOE-based NMR techniques, diffusion experiments and relaxation methods, are excellent tools to investigate how ligands interact with their receptors. Here we describe the key structural information that can be achieved on binding processes thanks to the combined used of advanced NMR and computational methods. Saturation Transfer Difference NMR (STD-NMR), WaterLOGSY, diffusion- and relaxation-based experiments, together with tr-NOE techniques allow, indeed, to investigate the ligand behavior when bound to a receptor, determining, among others, the epitope map of the ligand and its bioactive conformation. The combination of these NMR techniques with computational methods, including docking, molecular dynamics and CORCEMA-ST analysis, permits to define and validate an accurate 3D model of protein-ligand complexes.

Structural transitions in Orb2 prion-like domain relevant for functional aggregation in memory consolidation

The recent structural elucidation of ex vivo Drosophila Orb2 fibrils revealed a novel amyloid formed by interdigitated Gln and His residue side chains belonging to the prion-like domain. However, atomic-level details on the conformational transitions associated with memory consolidation remain unknown. Here, we have characterized the nascent conformation and dynamics of the prion-like domain (PLD) of Orb2A using a nonconventional liquid-state NMR spectroscopy strategy based on 13C detection to afford an essentially complete set of 13Cα, 13Cβ, 1Hα, and backbone 13CO and 15N assignments. At pH 4, where His residues are protonated, the PLD is disordered and flexible, except for a partially populated α-helix spanning residues 55-60, and binds RNA oligos, but not divalent cations. At pH 7, in contrast, His residues are predominantly neutral, and the Q/H segments adopt minor populations of helical structure, show decreased mobility and start to self-associate. At pH 7, the His residues do not bind RNA or Ca2+, but do bind Zn2+, which promotes further association. These findings represent a remarkable case of structural plasticity, based on which an updated model for Orb2A functional amyloidogenesis is suggested.

Molecular dynamics simulations in photosynthesis

Photosynthesis is regulated by a dynamic interplay between proteins, enzymes, pigments, lipids, and cofactors that takes place on a large spatio-temporal scale. Molecular dynamics (MD) simulations provide a powerful toolkit to investigate dynamical processes in (bio)molecular ensembles from the (sub)picosecond to the (sub)millisecond regime and from the Å to hundreds of nm length scale. Therefore, MD is well suited to address a variety of questions arising in the field of photosynthesis research. In this review, we provide an introduction to the basic concepts of MD simulations, at atomistic and coarse-grained level of resolution. Furthermore, we discuss applications of MD simulations to model photosynthetic systems of different sizes and complexity and their connection to experimental observables. Finally, we provide a brief glance on which methods provide opportunities to capture phenomena beyond the applicability of classical MD.

Theoretical and computational framework for the analysis of the relaxation properties of arbitrary spin systems. Application to high-resolution relaxometry

A wide variety of nuclear magnetic resonance experiments rely on the prediction and analysis of relaxation processes. Recently, innovative approaches have been introduced where the sample travels through a broad range of magnetic fields in the course of the experiment, such as dissolution dynamic nuclear polarization or high-resolution relaxometry. Understanding the relaxation properties of nuclear spin systems over orders of magnitude of magnetic fields is essential to rationalize the results of these experiments. For example, during a high-resolution relaxometry experiment, the absence of control of nuclear spin relaxation pathways during the sample transfers and relaxation delays leads to systematic deviations of polarization decays from an ideal mono-exponential decay with the pure longitudinal relaxation rate. These deviations have to be taken into account to describe quantitatively the dynamics of the system. Here, we present computational tools to (1) calculate analytical expressions of relaxation rates for a broad variety of spin systems and (2) use these analytical expressions to correct the deviations arising in high-resolution relaxometry experiments. These tools lead to a better understanding of nuclear spin relaxation, which is required to improve the sensitivity of many pulse sequences, and to better characterize motions in macromolecules.

Dynamic G Protein-Coupled Receptor Signaling Probed by Solution NMR Spectroscopy

Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool for investigating various dynamic features of G protein-coupled receptor (GPCR) signaling. In this Perspective, we focus on NMR techniques to characterize ligand-dependent conformational dynamics of GPCRs as well as the interaction of GPCRs with their environment and ligands. We also describe circumstances under which each technique should be applied, their advantages and disadvantages, and how they can be combined with other strategies to deepen the understanding of GPCR signaling at the molecular level.

Solution NMR insights into dynamic supramolecular assemblies of disordered amyloidogenic proteins

The extraordinary flexibility and structural heterogeneity of intrinsically disordered proteins (IDP) make them functionally versatile molecules. We have now begun to better understand their fundamental role in biology, however many aspects of their behaviour remain difficult to grasp experimentally. This is especially true for the intermolecular interactions which lead to the formation of transient or highly dynamic supramolecular self-assemblies, such as oligomers, aggregation intermediates and biomolecular condensates. Both the emerging functions and pathogenicity of these structures have stimulated great efforts to develop methodologies capable of providing useful insights. Significant progress in solution NMR spectroscopy has made this technique one of the most powerful to describe structural and dynamic features of IDPs within such assemblies at atomic resolution. Here, we review the most recent works that have illuminated key aspects of IDP assemblies and contributed significant advancements towards our understanding of the complex conformational landscape of prototypical disease-associated proteins. We also include a primer on some of the fundamental and innovative NMR methods being used in the discussed studies.

Three-dimensional structure of Megabalanus rosa Cement Protein 20 revealed by multi-dimensional NMR and molecular dynamics simulations

Barnacles employ a protein-based cement to firmly attach to immersed substrates. The cement proteins (CPs) have previously been identified and sequenced. However, the molecular mechanisms of adhesion are not well understood, in particular, because the three-dimensional molecular structure of CPs remained unknown to date. Here, we conducted multi-dimensional nuclear magnetic resonance (NMR) studies and molecular dynamics (MD) simulations of recombinant Megabalanus rosa Cement Protein 20 (rMrCP20). Our NMR results show that rMrCP20 contains three main folded domain regions intervened by two dynamic loops, resulting in multiple protein conformations that exist in equilibrium. We found that 12 out of 32 Cys in the sequence engage in disulfide bonds that stabilize the β-sheet domains owing to their placement at the extremities of β-strands. Another feature unveiled by NMR is the location of basic residues in turn regions that are exposed to the solvent, playing an important role for intermolecular contact with negatively charged surfaces. MD simulations highlight a highly stable and conserved β-motif (β7-β8), which may function as nuclei for amyloid-like nanofibrils previously observed in the cured adhesive cement. To the best of our knowledge, this is the first report describing the tertiary structure of an extracellular biological adhesive protein at the molecular level. This article is part of the theme issue 'Transdisciplinary approaches to the study of adhesion and adhesives in biological systems'.

Protein Dynamics from Accurate Low-Field Site-Specific Longitudinal and Transverse Nuclear Spin Relaxation

Nuclear magnetic relaxation provides invaluable quantitative site-specific information on the dynamics of complex systems. Determining dynamics on nanosecond time scales requires relaxation measurements at low magnetic fields incompatible with high-resolution NMR. Here, we use a two-field NMR spectrometer to measure carbon-13 transverse and longitudinal relaxation rates at a field as low as 0.33 T (proton Larmor frequency 14 MHz) in specifically labeled side chains of the protein ubiquitin. The use of radiofrequency pulses enhances the accuracy of measurements as compared to high-resolution relaxometry approaches, where the sample is moved in the stray field of the superconducting magnet. Importantly, we demonstrate that accurate measurements at a single low magnetic field provide enough information to characterize complex motions on low nanosecond time scales, which opens a new window for the determination of site-specific nanosecond motions in complex systems such as proteins.

Reducing bias in the analysis of solution-state NMR data with dynamics detectors

Nuclear magnetic resonance (NMR) is sensitive to dynamics on a wide range of correlation times. Recently, we have shown that analysis of relaxation rates via fitting to a correlation function with a small number of exponential terms could yield a biased characterization of molecular motion in solid-state NMR due to limited sensitivity of experimental data to certain ranges of correlation times. We introduced an alternative approach based on "detectors" in solid-state NMR, for which detector responses characterize motion for a range of correlation times and reduce potential bias resulting from the use of simple models for the motional correlation functions. Here, we show that similar bias can occur in the analysis of solution-state NMR relaxation data. We have thus adapted the detector approach to solution-state NMR, specifically separating overall tumbling motion from internal motions and accounting for contributions of chemical exchange to transverse relaxation. We demonstrate that internal protein motions can be described with detectors when the overall motion and the internal motions are statistically independent. We illustrate the detector analysis on ubiquitin with typical relaxation data sets recorded at a single high magnetic field or at multiple high magnetic fields and compare with results of model-free analysis. We also compare our methodology to LeMaster's method of dynamics analysis.

Accurate measurement of dipole/dipole transverse cross-correlated relaxation Γ 2 in methylenes and primary amines of uniformly 13 C / 15 N -labeled proteins

Side chains possess a broader conformational space (compared to the backbone) and are directly affected by intra- and intermolecular interactions, hence their dynamics and the corresponding NMR relaxation data are more sensitive and informative. Nevertheless, transverse relaxation inXH 2 (X = 13 C or 15 N ) spin systems is predominantly non-measurable in uniformly 13 C / 15 N -labeled proteins due to cross-correlation effects. In the present publication, we propose a number of pulse sequences for accurate and precise measurement of the dipole-dipole transverse cross-correlated relaxation rateΓ 2XH,XH ′ , which, similarly toR 2 measurements, provides information about the amplitudes of intramolecular dynamics. The suggested approach has allowed us to circumvent a number of obstacles that were limiting earlier applications ofΓ 2XH,XH ′ : (1) impossibility of transmission of the central component of the triplet ofXH 2 group to 1 H -acquisition via INEPT has been solved by transmission of the averaged signal of "inner" and "outer" components of the triplet; (2) direct recording of the entire triplets resulting in substantial overlap of side chain signals has been replaced by recording of individual singlets with the use of1 J CH -modulated approach and constant-time evolution; (3) low sensitivity has been enhanced via proton acquisition which required special attention to a zero-quantum coherence evolution. The proposed method expands the set of "dynamics sensors" covering protein side chains and substantially improves the quality and the level of detail of experimental data describing dynamic processes in proteins and protein complexes.

Chemical Exchange

The phenomenon of chemical or conformational exchange in NMR spectroscopy has enabled detailed characterization of time-dependent aspects of biomolecular function, including folding, molecular recognition, allostery, and catalysis, on timescales from microsecond to second. Importantly, NMR methods based on a variety of spin relaxation parameters have been developed that provide quantitative information on interconversion kinetics, thermodynamic properties, and structural features of molecular states populated to a fraction of a percent at equilibrium and otherwise unobservable by other NMR approaches. The ongoing development of more sophisticated experimental techniques and the necessity to apply these methods to larger and more complex molecular systems engenders a corresponding need for theoretical advances describing such techniques and facilitating data analysis in applications. This review surveys current aspects of the theory of chemical exchange, as utilized in ZZ-exchange; Hahn and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo; and R_1ρ, chemical exchange saturation transfer (CEST), and dark state saturation transfer (DEST) spin-locking experiments. The review emphasizes theoretical results for kinetic topologies with more than two interconverting states, both to obtain compact analytical forms suitable for data analysis and to establish conditions for distinguishability between alternative kinetic schemes.

15N transverse relaxation measurements for the characterization of µs-ms dynamics are deteriorated by the deuterium isotope effect on 15N resulting from solvent exchange

¹⁵N R₂ relaxation measurements are key for the elucidation of the dynamics of both folded and intrinsically disordered proteins (IDPs). Here we show, on the example of the intrinsically disordered protein α-synuclein and the folded domain PDZ2, that at physiological pH and near physiological temperatures amide-water exchange can severely skew Hahn-echo based ¹⁵N R₂ relaxation measurements as well as low frequency data points in CPMG relaxation dispersion experiments. The nature thereof is the solvent exchange with deuterium in the sample buffer, which modulates the ¹⁵N chemical shift tensor via the deuterium isotope effect, adding to the apparent relaxation decay which leads to systematic errors in the relaxation data. This results in an artificial increase of the measured apparent ¹⁵N R₂ rate constants-which should not be mistaken with protein inherent chemical exchange contributions, R_ex, to ¹⁵N R₂. For measurements of ¹⁵N R₂ rate constants of IDPs and folded proteins at physiological temperatures and pH, we recommend therefore the use of a very low D₂O molar fraction in the sample buffer, as low as 1%, or the use of an external D₂O reference along with a modified ¹⁵N R₂ Hahn-echo based experiment. This combination allows for the measurement of R_ex contributions to ¹⁵N R₂ originating from conformational exchange in a time window from µs to ms.

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.

An Evaluation of the Potential of NMR Spectroscopy and Computational Modelling Methods to Inform Biopharmaceutical Formulations

Protein-based therapeutics are considered to be one of the most important classes of pharmaceuticals on the market. The growing need to prolong stability of high protein concentrations in liquid form has proven to be challenging. Therefore, significant effort is being made to design formulations which can enable the storage of these highly concentrated protein therapies for up to 2 years. Currently, the excipient selection approach involves empirical high-throughput screening, but does not reveal details on aggregation mechanisms or the molecular-level effects of the formulations under storage conditions. Computational modelling approaches have the potential to elucidate such mechanisms, and rapidly screen in silico prior to experimental testing. Nuclear Magnetic Resonance (NMR) spectroscopy can also provide complementary insights into excipient–protein interactions. This review will highlight the underpinning principles of molecular modelling and NMR spectroscopy. It will also discuss the advancements in the applications of computational and NMR approaches in investigating excipient–protein interactions.

Simple and Robust Study of Backbone Dynamics of Crystalline Proteins Employing 1H-15N Dipolar Coupling Dispersion

We report a new solid-state multidimensional NMR approach based on the cross-polarization with variable-contact pulse sequence [ Paluch , P. ; Pawlak , T. ; Amoureux , J.-P. ; Potrzebowski , M. J. J. Magn. Reson. 233 , 2013 , 56 ], with ¹H inverse detection and very fast magic angle spinning (ν_R = 60 kHz), dedicated to the measurement of local molecular motions of ¹H-¹⁵N vectors. The introduced three-dimensional experiments, ¹H-¹⁵N-¹H and hCA(N)H, are particularly useful for the study of molecular dynamics of proteins and other complex structures. The applicability and power of this methodology have been revealed by employing as a model sample the GB-1 small protein doped with Na₂CuEDTA. The results clearly prove that the dispersion of ¹H-¹⁵N dipolar coupling constants well correlates with higher order structure of the protein. Our approach complements the conventional studies and offers a fast and reasonably simple method.

Determination of Protein ps-ns Motions by High-Resolution Relaxometry

Many of the functions of biomacromolecules can be rationalized by the characterization of their conformational energy landscapes: the structures of the dominant states, transitions between states and motions within states. Nuclear magnetic resonance (NMR) spectroscopy is the technique of choice to study internal motions in proteins. The determination of motions on picosecond to nanosecond timescales requires the measurement of nuclear spin relaxation rates at multiple magnetic fields. High sensitivity and resolution are obtained only at high magnetic fields, so that, until recently, site-specific relaxation rates in biomolecules were only measured over a narrow range of high magnetic fields. This limitation was particularly striking for the quantification of motions on nanosecond timescales, close to the correlation time for overall rotational diffusion. High-resolution relaxometry is an emerging technique to investigate picosecond-nanosecond motions of proteins. This approach uses a high-field NMR spectrometer equipped with a sample shuttle device, which allows for the measurement of the relaxation rate constants at low magnetic fields, while preserving the sensitivity and resolution of a high-field NMR spectrometer. The combined analysis of high-resolution relaxometry and standard high-field relaxation data provides a more accurate description of the dynamics of proteins, in particular in the nanosecond range. The purpose of this chapter is to describe how to perform high-resolution relaxometry experiments and how to analyze the rates measured with this technique.

Field-cycling NMR experiments in an ultra-wide magnetic field range: relaxation and coherent polarization transfer

An experimental method is described allowing fast field-cycling Nuclear Magnetic Resonance (NMR) experiments over a wide range of magnetic fields from 5 nT to 10 T. The method makes use of a hybrid technique: the high field range is covered by positioning the sample in the inhomogeneous stray field of the NMR spectrometer magnet. For fields below 2 mT a magnetic shield is mounted on top of the spectrometer; inside the shield the magnetic field is controlled by a specially designed coil system. This combination allows us to measure T1-relaxation times and nuclear Overhauser effect parameters over the full range in a routine way. For coupled proton-carbon spin systems relaxation with a common T1 is found at low fields, where the spins are "strongly coupled". In some cases, experiments at ultralow fields provide access to heteronuclear long-lived spin states. Efficient coherent polarization transfer is seen for proton-carbon spin systems at ultralow fields as follows from the observation of quantum oscillations in the polarization evolution. Applications to analysis and the manipulation of heteronuclear spin systems are discussed.

Fast evaluation of protein dynamics from deficient 15N relaxation data

Simple and convenient method of protein dynamics evaluation from the insufficient experimental ¹⁵N relaxation data is presented basing on the ratios, products, and differences of longitudinal and transverse ¹⁵N relaxation rates obtained at a single magnetic field. Firstly, the proposed approach allows evaluating overall tumbling correlation time (nanosecond time scale). Next, local parameters of the model-free approach characterizing local mobility of backbone amide N-H vectors on two different time scales, S² and R _ex , can be elucidated. The generalized order parameter, S², describes motions on the time scale faster than the overall tumbling correlation time (pico- to nanoseconds), while the chemical exchange term, R _ex , identifies processes slower than the overall tumbling correlation time (micro- to milliseconds). Advantages and disadvantages of different methods of data handling are thoroughly discussed.

Structure and Dynamics of an Intrinsically Disordered Protein Region That Partially Folds upon Binding by Chemical-Exchange NMR

Many intrinsically disordered proteins (IDPs) and protein regions (IDRs) engage in transient, yet specific, interactions with a variety of protein partners. Often, if not always, interactions with a protein partner lead to partial folding of the IDR. Characterizing the conformational space of such complexes is challenging: in solution-state NMR, signals of the IDR in the interacting region become broad, weak, and often invisible, while X-ray crystallography only provides information on fully ordered regions. There is thus a need for a simple method to characterize both fully and partially ordered regions in the bound state of IDPs. Here, we introduce an approach based on monitoring chemical exchange by NMR to investigate the state of an IDR that folds upon binding through the observation of the free state of the protein. Structural constraints for the bound state are obtained from chemical shifts, and site-specific dynamics of the bound state are characterized by relaxation rates. The conformation of the interacting part of the IDR was determined and subsequently docked onto the structure of the folded partner. We apply the method to investigate the interaction between the disordered C-terminal region of Artemis and the DNA binding domain of Ligase IV. We show that we can accurately reproduce the structure of the core of the complex determined by X-ray crystallography and identify a broader interface. The method is widely applicable to the biophysical investigation of complexes of disordered proteins and folded proteins.

A Site-Specific Study of the Magnetic Field-Dependent Proton Spin Relaxation of an Iridium N-Heterocyclic Carbene Complex

We report a study of proton spin relaxation of an Iridium N-heterocyclic carbene complex [Ir(COD)(IMes)Cl] complex (where COD=1,5-cyclooctadiene, Imes=1,3-bis(2,4,6-trimethylphenyl)imidazol-2-ylidene). This compound is a pre-catalyst of the most efficient complex allowing the signal amplification by reversible exchange (SABRE) effect, relevant for enhancing weak signals in nuclear magnetic resonance (NMR). An important feature of the study is a combination of relaxation measurements over a wide field range with high-resolution NMR detection. As a result, we are able to measure nuclear magnetic relaxation dispersion (NMRD) curves in the field range 0.1 mT–16.4 T (corresponding to the frequency range 4 kHz–700 MHz) for individual protons in the complex under study. This attractive possibility enables determination of the motional correlation times, τc , for the individual protons by analyzing the features in the NMRD curves (increase of the relaxation times) appearing at the magnetic fields where ωτc ≈1 (here ω is the proton Larmor precession frequency at a given field strength). The following correlation times were determined: (1.3±0.1) ns for the protons of imidazol-2-ylidene, (0.96±0.1) ns for the ortho-protons of two phenyl moieties and (0.95±0.2) ns for the protons of methyl groups. Additionally, we report low-field features coming from “strong coupling” of the protons. One should note that such features must not be misinterpreted by associating them with motional features. From the low-field features we obtain consistent estimates for the proton spin-spin interactions. The analysis of motional correlation times is also of importance for interpretation of spin order transfer from parahydrogen to various substrates in transient organometallic complexes (termed the SABRE effect) at high magnetic field.

High-resolution two-field nuclear magnetic resonance spectroscopy

Two-field NMR provides correlations of nuclear spins at the most favourable magnetic fields in a single experiment.