Probing site-specific conformational distributions in protein folding with solid-state NMR

Probing Chemical Equilibrium in Frozen Sodium Phosphate Buffer Solution by 31P Solid-State NMR

Phosphate buffers are crucial for cryopreservative stability in pharmaceuticals, food processing, biomedical sciences, and biology. However, their freeze concentrates lack quantitative characterization, especially regarding the physicochemical properties of phosphate salt species in equilibrium at subzero temperatures. This study employs 31P solid-state NMR (ssNMR) to analyze frozen sodium phosphate (NaP) solutions, providing insights into phase composition, ionic strength, and pH. For the first time, we have directly quantified phosphate species in frozen NaP buffer, including crystallized disodium phosphate dodecahydrate (Na2HPO4·12H2O) content and the concentrations of H2PO4- and HPO42- in the freeze concentrate. This enabled the calculation of the pH as well as the ionic strength in the freeze concentrate. Trehalose effectively mitigated pH shifts in buffer solutions by preventing the selective crystallization of salt, a spectroscopic phenomenon not previously observed experimentally.

Quantification of Residual Water in Pharmaceutical Frozen Solutions Via 1H Solid-State NMR

Freezing is essential for the stability of biological drug substances and products, particularly in frozen solution formulations and during the primary drying of lyophilized preparations. However, the unfrozen segment within the frozen matrix can alter solute concentration, ionic strength, and stabilizer crystallization, posing risks of increased biophysical instability and faster chemical degradation. While quantifying the unfrozen water content is important for designing stable biopharmaceuticals, there is a lack of analytical techniques for in situ quantitative measurements. In this study, we introduce a 1H magic angle spinning NMR technique to identify the freezing point (Tice) and quantify mobile water content in frozen biologics, applying this method to analyze the freezing of a commercial high-concentration drug product, Dupixent®. Our results demonstrate that water freezing is influenced by buffer salt properties and formulation composition, including the presence of sugar cryoprotectants and protein concentration. Additionally, the 1H chemical shift can probe pH in the unfrozen phase, potentially predicting the microenvironmental acidity in the frozen state. Our proposed methodology provides fresh insights into the analysis of freeze-concentrated solutions, enhancing our understanding of the stability of frozen and lyophilized biopharmaceuticals.

Solid-State NMR Spectroscopy to Probe State and Phase Transitions in Frozen Solutions

Freezing is commonly encountered during the processing and storage of biomacromolecule products. Therefore, understanding the phase and state transitions in pharmaceutical frozen solutions is crucial for the rational development of biopharmaceuticals. Solid-state nuclear magnetic resonance spectroscopy (ssNMR) was used to analyze solutions containing sodium phosphate buffer, histidine, and trehalose. Upon freezing, crystallization of disodium phosphate hydrogen dodecahydrate (Na2HPO4·12H2O, DPDH) and histidine was identified using 31P and 13C ssNMR, respectively, and confirmed by synchrotron X-ray diffractometry (SXRD). Using histidine as a molecular probe and based on the chemical shifts of atoms of interest, the pH of the freeze concentrate was measured. The unfrozen water content in freeze concentrates was quantified by 1H single pulse experiments. 13C-insensitive nuclei enhancement by polarization transfer (INEPT) and cross-polarization (CP) experiments were used as orthogonal tools to characterize the solutes in a "mobile" and a more "solid-like" state in the freeze-concentrated solutions, respectively. The above analyses were applied to a commercial monoclonal antibody (mAb) formulation of dupilumab. This work further establishes ssNMR spectroscopy as a highly capable biophysical tool to investigate the attributes of biopharmaceuticals and thereby provide insights into process optimization and formulation development.

A high-field cellular DNP-supported solid-state NMR approach to study proteins with sub-cellular specificity

Studying the structural aspects of proteins within sub-cellular compartments is of growing interest. Dynamic nuclear polarization supported solid-state NMR (DNP-ssNMR) is uniquely suited to provide such information, but critically lacks the desired sensitivity and resolution. Here we utilize SNAPol-1, a novel biradical, to conduct DNP-ssNMR at high-magnetic fields (800 MHz/527 GHz) inside HeLa cells and isolated cell nuclei electroporated with [13C,15N] labeled ubiquitin. We report that SNAPol-1 passively diffuses and homogenously distributes within whole cells and cell nuclei providing ubiquitin spectra of high sensitivity and remarkably improved spectral resolution. For cell nuclei, physical enrichment facilitates a further 4-fold decrease in measurement time and provides an exclusive structural view of the nuclear ubiquitin pool. Taken together, these advancements enable atomic interrogation of protein conformational plasticity at atomic resolution and with sub-cellular specificity.

Conformational ensembles explain NMR spectra of frozen intrinsically disordered proteins

Protein regions which are intrinsically disordered, exist as an ensemble of rapidly interconverting structures. Cooling proteins to cryogenic temperatures for dynamic nuclear polarization (DNP) magic angle spinning (MAS) NMR studies suspends most of the motions, resulting in peaks that are broad but not featureless. To demonstrate that detailed conformational restraints can be retrieved from the peak shapes of frozen proteins alone, we developed and used a simulation framework to assign peak features to conformers in the ensemble. We validated our simulations by comparing them to spectra of α-synuclein acquired under different experimental conditions. Our assignments of peaks to discrete dihedral angle populations suggest that structural constraints are attainable under cryogenic conditions. The ability to infer ensemble populations from peak shapes has important implications for DNP MAS NMR studies of proteins with regions of disorder in living cells because chemical shifts are the most accessible measured parameter.

Characterization of conformational heterogeneity via higher-dimensionality, proton-detected solid-state NMR

Site-specific heterogeneity of solid protein samples can be exploited as valuable information to answer biological questions ranging from thermodynamic properties determining fibril formation to protein folding and conformational stability upon stress. In particular, for proteins of increasing molecular weight, however, site-resolved assessment without residue-specific labeling is challenging using established methodology, which tends to rely on carbon-detected 2D correlations. Here we develop purely chemical-shift-based approaches for assessment of relative conformational heterogeneity that allows identification of each residue via four chemical-shift dimensions. High dimensionality diminishes the probability of peak overlap in the presence of multiple, heterogeneously broadened resonances. Utilizing backbone dihedral-angle reconstruction from individual contributions to the peak shape either via suitably adapted prediction routines or direct association with a relational database, the methods may in future studies afford assessment of site-specific heterogeneity of proteins without site-specific labeling.

Solid-state NMR methods for the characterization of bioconjugations and protein-material interactions

Protein solid-state NMR has evolved dramatically over the last two decades, with the development of new hardware and sample preparation methodologies. This technique is now ripe for complex applications, among which one can count bioconjugation, protein chemistry and functional biomaterials. In this review, we provide our account on this aspect of protein solid-state NMR.

Time-resolved solid state NMR of biomolecular processes with millisecond time resolution

We review recent efforts to develop and apply an experimental approach to the structural characterization of transient intermediate states in biomolecular processes that involve large changes in molecular conformation or assembly state. This approach depends on solid state nuclear magnetic resonance (ssNMR) measurements that are performed at very low temperatures, typically 25-30 K, with signal enhancements from dynamic nuclear polarization (DNP). This approach also involves novel technology for initiating the process of interest, either by rapid mixing of two solutions or by a rapid inverse temperature jump, and for rapid freezing to trap intermediate states. Initiation by rapid mixing or an inverse temperature jump can be accomplished in approximately-one millisecond. Freezing can be accomplished in approximately 100 microseconds. Thus, millisecond time resolution can be achieved. Recent applications to the process by which the biologically essential calcium sensor protein calmodulin forms a complex with one of its target proteins and the process by which the bee venom peptide melittin converts from an unstructured monomeric state to a helical, tetrameric state after a rapid change in pH or temperature are described briefly. Future applications of millisecond time-resolved ssNMR are also discussed briefly.

Virus Structures and Dynamics by Magic-Angle Spinning NMR

Techniques for atomic-resolution structural biology have evolved during the past several decades. Breakthroughs in instrumentation, sample preparation, and data analysis that occurred in the past decade have enabled characterization of viruses with an unprecedented level of detail. Here we review the recent advances in magic-angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy for structural analysis of viruses and viral assemblies. MAS NMR is a powerful method that yields information on 3D structures and dynamics in a broad range of experimental conditions. After a brief introduction, we discuss recent structural and functional studies of several viruses investigated with atomic resolution at various levels of structural organization, from individual domains of a membrane protein reconstituted into lipid bilayers to virus-like particles and intact viruses. We present examples of the unique information revealed by MAS NMR about drug binding, conduction mechanisms, interactions with cellular host factors, and DNA packaging in biologically relevant environments that are inaccessible by other methods.

Soluble State of Villin Headpiece Protein as a Tool in the Assessment of MD Force Fields

Protein self-assembly plays an important role in cellular processes. Whereas molecular dynamics (MD) represents a powerful tool in studying assembly mechanisms, its predictions depend on the accuracy of underlying force fields, which are known to overly promote protein assembly. We here examine villin headpiece domain, HP36, which remains soluble at concentrations amenable to MD studies. The experimental characterization of soluble HP36 at concentrations of 0.05 to 1 mM reveals concentration-independent 90% monomeric and 10% dimeric populations. Extensive all-atom MD simulations at two protein concentrations, 0.9 and 8.5 mM, probe the HP36 dimer population, stability, and kinetics of dimer formation within two MD force fields, Amber ff14SB and CHARMM36m. MD results demonstrate that whereas CHARMM36m captures experimental HP36 monomer populations at the lower concentration, both force fields overly promote HP36 association at the higher concentration. Moreover, contacts stabilizing HP36 dimers are force-field-dependent. CHARMM36m produces consistently higher HP36 monomer populations, lower association rates, and weaker dependence of these quantities on the protein concentration than Amber ff14SB. Nonetheless, the highest monomer populations and dissociation constants are observed when the TIP3P water model in Amber ff14SB is replaced by TIP4P/2005, showcasing the critical role of the water model in addressing the protein solubility problem in MD.

Quantitative Analysis of Protein Unfolded State Energetics: Experimental and Computational Studies Demonstrate That Non-Native Side-Chain Interactions Stabilize Local Native Backbone Structure

Proteins fold on relatively smooth free energy landscapes which are biased toward the native state, but even simple topologies which fold rapidly can experience roughness on their free energy landscape. The details of these interactions are difficult to elucidate experimentally. Closely related to the problem of deciphering the details of the free energy landscape is the problem of defining the interactions in the denatured state ensemble (DSE) which is populated under native conditions, that is, under conditions where the native state is stable. The DSE of many proteins deviates from random coil models, but quantifying and defining the energetics of the transiently populated interactions in this ensemble is extremely challenging. Characterization of the DSE of proteins which fold to compact structures is also relevant to studies of intrinsically disordered proteins (IDPs) since interactions in the dynamic ensemble populated by IDPs can modulate their behavior. Here we show how experimental thermodynamic and pK_a measurements can be combined with computational thermodynamic integration to quantify interactions in the DSE. We show that non-native side chain interactions can stabilize native backbone structure in the DSE and demonstrate that that even rapidly folding proteins can form energetically significant non-native interactions in their DSE. As an example, we characterize a non-native salt bridge that stabilizes local native backbone structure in the DSE of a widely studied fast-folding protein, the villin headpiece helical domain. The combined computational experimental approach is applicable to other protein unfolded states and provides insight that is impossible to obtain with either method alone.

Advances in studying protein disorder with solid-state NMR

Solution NMR is a key tool to study intrinsically disordered proteins (IDPs), whose importance for biological function is widely accepted. However, disordered proteins are not limited to solution and are also found in non-soluble systems such as fibrils and membrane proteins. In this Trends article, I will discuss how solid-state NMR can be used to study disorder in non-soluble proteins. Techniques based on dipolar couplings can study static protein disorder which either occurs naturally as e.g. in spider silk or can be induced by freeze trapping IDPs or unfolded proteins. In this case, structural ensembles are directly reflected by a static distribution of dihedral angels that can be determined by the distribution of chemical shifts or other methods. Techniques based on J-couplings can detect dynamic protein disorder under MAS. In this case, only average chemical shifts are measured but disorder can be characterized with a variety of data including secondary chemical shifts, relaxation rates, paramagnetic relaxation enhancements, or residual dipolar couplings. I describe both technical aspects and examples of solid-state NMR on protein disorder and end the article with a discussion of challenges and opportunities of this emerging field.

Integrating cryo-EM and NMR data

Single-particle cryo-electron microscopy (cryo-EM) is increasingly used as a technique to determine the atomic structure of challenging biological systems. Recent advances in microscope engineering, electron detection, and image processing have allowed the structural determination of bigger and more flexible targets than possible with the complementary techniques X-ray crystallography and NMR spectroscopy. However, there exist many biological targets for which atomic resolution cannot be currently achieved with cryo-EM, making unambiguous determination of the protein structure impossible. Although determining the structure of large biological systems using solely NMR is often difficult, highly complementary experimental atomic-level data for each molecule can be derived from the spectra, and used in combination with cryo-EM data. We review here strategies with which both techniques can be synergistically combined, in order to reach detail and understanding unattainable by each technique acting alone; and the types of biological systems for which such an approach would be desirable.

Application of millisecond time-resolved solid state NMR to the kinetics and mechanism of melittin self-assembly

Common experimental approaches for characterizing structural conversion processes such as protein folding and self-assembly do not report on all aspects of the evolution from an initial state to the final state. Here, we demonstrate an approach that is based on rapid mixing, freeze-trapping, and low-temperature solid-state NMR (ssNMR) with signal enhancements from dynamic nuclear polarization (DNP). Experiments on the folding and tetramerization of the 26-residue peptide melittin following a rapid pH jump show that multiple aspects of molecular structure can be followed with millisecond time resolution, including secondary structure at specific isotopically labeled sites, intramolecular and intermolecular contacts between specific pairs of labeled residues, and overall structural order. DNP-enhanced ssNMR data reveal that conversion of conformationally disordered melittin monomers at low pH to α-helical conformations at neutral pH occurs on nearly the same timescale as formation of antiparallel melittin dimers, about 6 to 9 ms for 0.3 mM melittin at 24 °C in aqueous solution containing 20% (vol/vol) glycerol and 75 mM sodium phosphate. Although stopped-flow fluorescence data suggest that melittin tetramers form quickly after dimerization, ssNMR spectra show that full structural order within melittin tetramers develops more slowly, in ∼60 ms. Time-resolved ssNMR is likely to find many applications to biomolecular structural conversion processes, including early stages of amyloid formation, viral capsid formation, and protein-protein recognition.

Dynamic Nuclear Polarization Magic-Angle Spinning Nuclear Magnetic Resonance Combined with Molecular Dynamics Simulations Permits Detection of Order and Disorder in Viral Assemblies

We report dynamic nuclear polarization (DNP)-enhanced magic-angle spinning (MAS) NMR spectroscopy in viral capsids from HIV-1 and bacteriophage AP205. Viruses regulate their life cycles and infectivity through modulation of their structures and dynamics. While static structures of capsids from several viruses are now accessible with near-atomic-level resolution, atomic-level understanding of functionally important motions in assembled capsids is lacking. We observed up to 64-fold signal enhancements by DNP, which permitted in-depth analysis of these assemblies. For the HIV-1 CA assemblies, a remarkably high spectral resolution in the 3D and 2D heteronuclear data sets permitted the assignment of a significant fraction of backbone and side-chain resonances. Using an integrated DNP MAS NMR and molecular dynamics (MD) simulation approach, the conformational space sampled by the assembled capsid at cryogenic temperatures was mapped. Qualitatively, a remarkable agreement was observed for the experimental ¹³C/¹⁵N chemical shift distributions and those calculated from substructures along the MD trajectory. Residues that are mobile at physiological temperatures are frozen out in multiple conformers at cryogenic conditions, resulting in broad experimental and calculated chemical shift distributions. Overall, our results suggest that DNP MAS NMR measurements in combination with MD simulations facilitate a thorough understanding of the dynamic signatures of viral capsids.

DNP-Enhanced MAS NMR: A Tool to Snapshot Conformational Ensembles of α-Synuclein in Different States

Intrinsically disordered proteins dynamically sample a wide conformational space and therefore do not adopt a stable and defined three-dimensional conformation. The structural heterogeneity is related to their proper functioning in physiological processes. Knowledge of the conformational ensemble is crucial for a complete comprehension of this kind of proteins. We here present an approach that utilizes dynamic nuclear polarization-enhanced solid-state NMR spectroscopy of sparsely isotope-labeled proteins in frozen solution to take snapshots of the complete structural ensembles by exploiting the inhomogeneously broadened line-shapes. We investigated the intrinsically disordered protein α-synuclein (α-syn), which plays a key role in the etiology of Parkinson’s disease, in three different physiologically relevant states. For the free monomer in frozen solution we could see that the so-called “random coil conformation” consists of α-helical and β-sheet-like conformations, and that secondary chemical shifts of neighboring amino acids tend to be correlated, indicative of frequent formation of secondary structure elements. Based on these results, we could estimate the number of disordered regions in fibrillar α-syn as well as in α-syn bound to membranes in different protein-to-lipid ratios. Our approach thus provides quantitative information on the propensity to sample transient secondary structures in different functional states. Molecular dynamics simulations rationalize the results.

Hyperpolarized MAS NMR of unfolded and misfolded proteins

In this article we give an overview over the use of DNP-enhanced solid-state NMR spectroscopy for the investigation of unfolded, disordered and misfolded proteins. We first provide an overview over studies in which DNP spectroscopy has successfully been applied for the structural investigation of well-folded amyloid fibrils formed by short peptides as well as full-length proteins. Sample cooling to cryogenic temperatures often leads to severe line broadening of resonance signals and thus a loss in resolution. However, inhomogeneous line broadening at low temperatures provides valuable information about residual dynamics and flexibility in proteins, and, in combination with appropriate selective isotope labeling techniques, inhomogeneous linewidths in disordered proteins or protein regions may be exploited for evaluation of conformational ensembles. In the last paragraph we highlight some recent studies where DNP-enhanced MAS-NMR-spectroscopy was applied to the study of disordered proteins/protein regions and inhomogeneous sample preparations.

Heterogeneity in the Folding of Villin Headpiece Subdomain HP36

Small single domain proteins that fold on the microsecond time scale have been the subject of intense interest as models for probing the complexity of folding energy landscapes. The villin headpiece subdomain (HP36) has been extensively studied because of its simple three helix structure, ultrafast folding lifetime of a few microseconds, and stable native fold. We have previously shown that folding as measured by a single ¹³C═¹⁸O isotopic label on residue A57 in helix 2 occurs at a different rate than that measured by global probes of folding, indicating noncooperative complexity in the folding of HP36. In order to determine whether this complexity reflects intermediates or parallel pathways over a small activation barrier, ¹³C═¹⁸O labels were individually incorporated at six different positions in HP36, including into all 3 helices. The equilibrium thermal unfolding transitions and the folding/unfolding dynamics were monitored using the unique IR signature of the ¹³C═¹⁸O label by temperature dependent FTIR and temperature jump IR spectroscopy, respectively. Equilibrium experiments reveal that the ¹³C═¹⁸O labels at different positions in HP36 show drastic differences in the midpoint of their transitions ( T_m), ranging from 45 to 67 °C. Heterogeneity is also observed in the relaxation kinetics; there are differences in the microsecond phase when different labeled positions are probed. At a final temperature of 45 °C, the relaxation rate for ¹³C═¹⁸O A57 is 2.4e + 05 s^-1 whereas for ¹³C═¹⁸O L69 HP36 the relaxation rate is 5.1e + 05 s^-1, two times faster. The observation of site-dependent midpoints for the equilibrium unfolding transitions and differences in the relaxation rates of the labeled positions enables us to probe the progressive accumulation of the folded structure, providing insight into the microscopic details of the folding mechanism.

Hexagonal ice in pure water and biological NMR samples

Ice, in addition to "liquid" water and protein, is an important component of protein samples for NMR spectroscopy at subfreezing temperatures but it has rarely been observed spectroscopically in this context. We characterize its spectroscopic behavior in the temperature range from 100 to 273 K, and find that it behaves like pure water ice. The interference of magic-angle spinning (MAS) as well as rf multiple-pulse sequences with Bjerrum-defect motion greatly influences the ice spectra.

Nuclear magnetic resonance (NMR) applied to membrane–protein complexes

Increasing evidence suggests that most proteins occur and function in complexes rather than as isolated entities when embedded in cellular membranes. Nuclear magnetic resonance (NMR) provides increasing possibilities to study structure, dynamics and assembly of such systems. In our review, we discuss recent methodological progress to study membrane–protein complexes (MPCs) by NMR, starting with expression, isotope-labeling and reconstitution protocols. We review approaches to deal with spectral complexity and limited spectral spectroscopic sensitivity that are usually encountered in NMR-based studies of MPCs. We highlight NMR applications in various classes of MPCs, including G-protein-coupled receptors, ion channels and retinal proteins and extend our discussion to protein–protein complexes that span entire cellular compartments or orchestrate processes such as protein transport across or within membranes. These examples demonstrate the growing potential of NMR-based studies of MPCs to provide critical insight into the energetics of protein–ligand and protein–protein interactions that underlie essential biological functions in cellular membranes.

Low-temperature dynamic nuclear polarization with helium-cooled samples and nitrogen-driven magic-angle spinning

We describe novel instrumentation for low-temperature solid state nuclear magnetic resonance (NMR) with dynamic nuclear polarization (DNP) and magic-angle spinning (MAS), focusing on aspects of this instrumentation that have not been described in detail in previous publications. We characterize the performance of an extended interaction oscillator (EIO) microwave source, operating near 264 GHz with 1.5 W output power, which we use in conjunction with a quasi-optical microwave polarizing system and a MAS NMR probe that employs liquid helium for sample cooling and nitrogen gas for sample spinning. Enhancement factors for cross-polarized (13)C NMR signals in the 100-200 range are demonstrated with DNP at 25K. The dependences of signal amplitudes on sample temperature, as well as microwave power, polarization, and frequency, are presented. We show that sample temperatures below 30K can be achieved with helium consumption rates below 1.3 l/h. To illustrate potential applications of this instrumentation in structural studies of biochemical systems, we compare results from low-temperature DNP experiments on a calmodulin-binding peptide in its free and bound states.

The unusual internal motion of the villin headpiece subdomain

The thermostable 36-residue subdomain of the villin headpiece (HP36) is the smallest known cooperatively folding protein. Although the folding and internal dynamics of HP36 and close variants have been extensively studied, there has not been a comprehensive investigation of side-chain motion in this protein. Here, the fast motion of methyl-bearing amino acid side chains is explored over a range of temperatures using site-resolved solution nuclear magnetic resonance deuterium relaxation. The squared generalized order parameters of methyl groups extensively spatially segregate according to motional classes. This has not been observed before in any protein studied using this methodology. The class segregation is preserved from 275 to 305 K. Motions detected in Helix 3 suggest a fast timescale of conformational heterogeneity that has not been previously observed but is consistent with a range of folding and dynamics studies. Finally, a comparison between the order parameters in solution with previous results based on solid-state nuclear magnetic resonance deuterium line shape analysis of HP36 in partially hydrated powders shows a clear disagreement for half of the sites. This result has significant implications for the interpretation of data derived from a variety of approaches that rely on partially hydrated protein samples.

Successive Stages of Amyloid-β Self-Assembly Characterized by Solid-State Nuclear Magnetic Resonance with Dynamic Nuclear Polarization

Self-assembly of amyloid-β (Aβ) peptides in human brain tissue leads to neurodegeneration in Alzheimer's disease (AD). Amyloid fibrils, whose structures have been extensively characterized by solid state nuclear magnetic resonance (ssNMR) and other methods, are the thermodynamic end point of Aβ self-assembly. Oligomeric and protofibrillar assemblies, whose structures are less well-understood, are also observed as intermediates in the assembly process in vitro and have been implicated as important neurotoxic species in AD. We report experiments in which the structural evolution of 40-residue Aβ (Aβ40) is monitored by ssNMR measurements on frozen solutions prepared at four successive stages of the self-assembly process. Measurements on transient intermediates are enabled by ssNMR signal enhancements from dynamic nuclear polarization (DNP) at temperatures below 30 K. DNP-enhanced ssNMR data reveal a monotonic increase in conformational order from an initial state comprised primarily of monomers and small oligomers in solution at high pH, to larger oligomers near neutral pH, to metastable protofibrils, and finally to fibrils. Surprisingly, the predominant molecular conformation, indicated by (13)C NMR chemical shifts and by side chain contacts between F19 and L34 residues, is qualitatively similar at all stages. However, the in-register parallel β-sheet supramolecular structure, indicated by intermolecular (13)C spin polarization transfers, does not develop before the fibril stage. This work represents the first application of DNP-enhanced ssNMR to the characterization of peptide or protein self-assembly intermediates.

NMR of sedimented, fibrillized, silica-entrapped and microcrystalline (metallo)proteins

Resolution and sensitivity in solid state NMR (SSNMR) can rival the results achieved by solution NMR, and even outperform them in the case of large systems. However, several factors affect the spectral quality in SSNMR samples, and not all systems turn out to be equally amenable for this methodology. In this review we attempt at analyzing the causes of this variable behavior and at providing hints to increase the chances of experimental success.

NMR-based structural biology enhanced by dynamic nuclear polarization at high magnetic field

Dynamic nuclear polarization (DNP) has become a powerful method to enhance spectroscopic sensitivity in the context of magnetic resonance imaging and nuclear magnetic resonance spectroscopy. We show that, compared to DNP at lower field (400 MHz/263 GHz), high field DNP (800 MHz/527 GHz) can significantly enhance spectral resolution and allows exploitation of the paramagnetic relaxation properties of DNP polarizing agents as direct structural probes under magic angle spinning conditions. Applied to a membrane-embedded K(+) channel, this approach allowed us to refine the membrane-embedded channel structure and revealed conformational substates that are present during two different stages of the channel gating cycle. High-field DNP thus offers atomic insight into the role of molecular plasticity during the course of biomolecular function in a complex cellular environment.

COFFDROP: A Coarse-Grained Nonbonded Force Field for Proteins Derived from All-Atom Explicit-Solvent Molecular Dynamics Simulations of Amino Acids

We describe the derivation of a set of bonded and nonbonded coarse-grained (CG) potential functions for use in implicit-solvent Brownian dynamics (BD) simulations of proteins derived from all-atom explicit-solvent molecular dynamics (MD) simulations of amino acids. Bonded potential functions were derived from 1 μs MD simulations of each of the 20 canonical amino acids, with histidine modeled in both its protonated and neutral forms; nonbonded potential functions were derived from 1 μs MD simulations of every possible pairing of the amino acids (231 different systems). The angle and dihedral probability distributions and radial distribution functions sampled during MD were used to optimize a set of CG potential functions through use of the iterative Boltzmann inversion (IBI) method. The optimized set of potential functions—which we term COFFDROP (COarse-grained Force Field for Dynamic Representation Of Proteins)—quantitatively reproduced all of the “target” MD distributions. In a first test of the force field, it was used to predict the clustering behavior of concentrated amino acid solutions; the predictions were directly compared with the results of corresponding all-atom explicit-solvent MD simulations and found to be in excellent agreement. In a second test, BD simulations of the small protein villin headpiece were carried out at concentrations that have recently been studied in all-atom explicit-solvent MD simulations by Petrov and Zagrovic (PLoS Comput. Biol.2014, 5, e1003638). The anomalously strong intermolecular interactions seen in the MD study were reproduced in the COFFDROP simulations; a simple scaling of COFFDROP’s nonbonded parameters, however, produced results in better accordance with experiment. Overall, our results suggest that potential functions derived from simulations of pairwise amino acid interactions might be of quite broad applicability, with COFFDROP likely to be especially useful for modeling unfolded or intrinsically disordered proteins.

Are current atomistic force fields accurate enough to study proteins in crowded environments?

The high concentration of macromolecules in the crowded cellular interior influences different thermodynamic and kinetic properties of proteins, including their structural stabilities, intermolecular binding affinities and enzymatic rates. Moreover, various structural biology methods, such as NMR or different spectroscopies, typically involve samples with relatively high protein concentration. Due to large sampling requirements, however, the accuracy of classical molecular dynamics (MD) simulations in capturing protein behavior at high concentration still remains largely untested. Here, we use explicit-solvent MD simulations and a total of 6.4 µs of simulated time to study wild-type (folded) and oxidatively damaged (unfolded) forms of villin headpiece at 6 mM and 9.2 mM protein concentration. We first perform an exhaustive set of simulations with multiple protein molecules in the simulation box using GROMOS 45a3 and 54a7 force fields together with different types of electrostatics treatment and solution ionic strengths. Surprisingly, the two villin headpiece variants exhibit similar aggregation behavior, despite the fact that their estimated aggregation propensities markedly differ. Importantly, regardless of the simulation protocol applied, wild-type villin headpiece consistently aggregates even under conditions at which it is experimentally known to be soluble. We demonstrate that aggregation is accompanied by a large decrease in the total potential energy, with not only hydrophobic, but also polar residues and backbone contributing substantially. The same effect is directly observed for two other major atomistic force fields (AMBER99SB-ILDN and CHARMM22-CMAP) as well as indirectly shown for additional two (AMBER94, OPLS-AAL), and is possibly due to a general overestimation of the potential energy of protein-protein interactions at the expense of water-water and water-protein interactions. Overall, our results suggest that current MD force fields may distort the picture of protein behavior in biologically relevant crowded environments.

Protein behavior is strongly affected by highly crowded and interaction-rich environments, i.e., typical conditions in both biologically relevant systems, such as the cellular interior, and solution-based structural experiments, including NMR and different spectroscopies. On the other hand, primarily because of limited computational power, molecular dynamics (MD) simulations, a premier high-resolution method for analyzing structure, dynamics and interactions of proteins, have been predominantly used to study individual proteins at infinite dilution. To fill this gap, we use MD simulations to study the behavior of wild-type (aggregation-resistant) and oxidatively damaged (aggregation-prone) forms of villin headpiece at high concentration, and reveal unexpected limitations and inaccuracies of modern-day MD force fields when it comes to modeling proteins at physiologically or experimentally relevant concentrations.

Electronic polarization stabilizes tertiary structure prediction of HP-36

Molecular dynamic (MD) simulations with both implicit and explicit solvent models have been carried out to study the folding dynamics of HP-36 protein. Starting from the extended conformation, the secondary structure of all three helices in HP-36 was formed in about 50 ns and remained stable in the remaining simulation. However, the formation of the tertiary structure was difficult. Although some intermediates were close to the native structure, the overall conformation was not stable. Further analysis revealed that the large structure fluctuation of loop and hydrophobic core regions was devoted mostly to the instability of the structure during MD simulation. The backbone root-mean-square deviation (RMSD) of the loop and hydrophobic core regions showed strong correlation with the backbone RMSD of the whole protein. The free energy landscape indicated that the distribution of main chain torsions in loop and turn regions was far away from the native state. Starting from an intermediate structure extracted from the initial AMBER simulation, HP-36 was found to generally fold to the native state under the dynamically adjusted polarized protein-specific charge (DPPC) simulation, while the peptide did not fold into the native structure when AMBER force filed was used. The two best folded structures were extracted and taken into further simulations in water employing AMBER03 charge and DPPC for 25 ns. Result showed that introducing polarization effect into interacting potential could stabilize the near-native protein structure.

Elucidating molecular motion through structural and dynamic filters of energy-minimized conformer ensembles

Complex RNA structures are constructed from helical segments connected by flexible loops that move spontaneously and in response to binding of small molecule ligands and proteins. Understanding the conformational variability of RNA requires the characterization of the coupled time evolution of interconnected flexible domains. To elucidate the collective molecular motions and explore the conformational landscape of the HIV-1 TAR RNA, we describe a new methodology that utilizes energy-minimized structures generated by the program “Fragment Assembly of RNA with Full-Atom Refinement (FARFAR)”. We apply structural filters in the form of experimental residual dipolar couplings (RDCs) to select a subset of discrete energy-minimized conformers and carry out principal component analyses (PCA) to corroborate the choice of the filtered subset. We use this subset of structures to calculate solution T₁ and T_1ρ relaxation times for ¹³C spins in multiple residues in different domains of the molecule using two simulation protocols that we previously published. We match the experimental T₁ times to within 2% and the T_1ρ times to within less than 10% for helical residues. These results introduce a protocol to construct viable dynamic trajectories for RNA molecules that accord well with experimental NMR data and support the notion that the motions of the helical portions of this small RNA can be described by a relatively small number of discrete conformations exchanging over time scales longer than 1 μs.

Practical considerations over spectral quality in solid state NMR spectroscopy of soluble proteins

Great theoretical and methodological advances are pushing the limits of resolution and sensitivity in solid state NMR (SSNMR). However, sample preparation remains a critical issue for the success of an experiment. The factors affecting spectral quality in SSNMR samples are discussed, examining cases encountered in the literature and presenting new experimental data. A discussion on resolution and sensitivity in sedimented solutes is framed in this context.

Origin of abrupt rise in deuteron NMR longitudinal relaxation times of protein methyl groups below 90 K

In order to examine the origin of the abrupt change in the temperature dependence of (2)H NMR longitudinal relaxation times observed previously for methyl groups of L69 in the hydrophobic core of villin headpiece protein at around 90 K (Vugmeyster et al. J. Am. Chem. Soc. 2010, 132, 4038-4039), we extended the measurements to several other methyl groups in the hydrophobic core. We show that, for all methyl groups, relaxation times experience a dramatic jump several orders of magnitude around this temperature. Theoretical modeling supports the conclusion that the origin of the apparent transition in the relaxation times is due to the existence of the distribution of conformers distinguished by their activation energy for methyl three-site hops. It is also crucial to take into account the differential contribution of individual conformers into overall signal intensity. When a particular conformer approaches the regime at which its three-site hop rate constant is on the order of the quadrupolar coupling interaction constant, the intensity of the signal due to this conformer experiences a sharp drop, thus changing the balance of the contributions of different conformers into the overall signal. As a result, the observed apparent transition in the relaxation rates can be explained without the assumption of an underlying transition in the rate constants. This work in combination with earlier results also shows that the model based on the distribution of conformers explains the relaxation behavior in the entire temperature range between 300 and 70 K.

Secondary structural analysis of proteins based on (13)C chemical shift assignments in unresolved solid-state NMR spectra enhanced by fragmented structure database

Magic-angle-spinning solid-state (13)C NMR spectroscopy is useful for structural analysis of non-crystalline proteins. However, the signal assignments and structural analysis are often hampered by the signal overlaps primarily due to minor structural heterogeneities, especially for uniformly-(13)C,(15)N labeled samples. To overcome this problem, we present a method for assigning (13)C chemical shifts and secondary structures from unresolved two-dimensional (13)C-(13)C MAS NMR spectra by spectral fitting, named reconstruction of spectra using protein local structures (RESPLS). The spectral fitting was conducted using databases of protein fragmented structures related to (13)C(α), (13)C(β), and (13)C' chemical shifts and cross-peak intensities. The experimental (13)C-(13)C inter- and intra-residue correlation spectra of uniformly isotope-labeled ubiquitin in the lyophilized state had a few broad peaks. The fitting analysis for these spectra provided sequence-specific C(α), C(β), and C' chemical shifts with an accuracy of about 1.5 ppm, which enabled the assignment of the secondary structures with an accuracy of 79 %. The structural heterogeneity of the lyophilized ubiquitin is revealed from the results. Test of RESPLS analysis for simulated spectra of five different types of proteins indicated that the method allowed the secondary structure determination with accuracy of about 80 % for the 50-200 residue proteins. These results demonstrate that the RESPLS approach expands the applicability of the NMR to non-crystalline proteins exhibiting unresolved (13)C NMR spectra, such as lyophilized proteins, amyloids, membrane proteins and proteins in living cells.

Glassy dynamics of protein methyl groups revealed by deuteron NMR

We investigated site-specific dynamics of key methyl groups in the hydrophobic core of chicken villin headpiece subdomain (HP36) over the temperature range between 298 and 140 K using deuteron solid-state NMR longitudinal relaxation measurements. The relaxation of the longitudinal magnetization is weakly nonexponential (glassy) at high temperatures and exhibits a stronger degree of nonexponentiality below about 175 K. In addition, the characteristic relaxation times deviate from the simple Arrhenius law. We interpret this behavior via the existence of distribution of activation energy barriers for the three-site methyl jumps, which originates from somewhat different methyl environments within the local energy landscape. The width of the distribution of the activation barriers for methyl jumps is rather significant, about 1.4 kJ/mol. Our experimental results and modeling allow for the description of the apparent change at about 175 K without invoking a specific transition temperature. For most residues in the core, the relaxation behavior at high temperatures points to the existence of conformational exchange between the substates of the landscape, and our model takes into account the kinetics of this process. The observed dynamics are the same for dry and hydrated protein. We also looked at the effect of F58L mutation inside the hydrophobic core on the dynamics of one of the residues and observed a significant increase in its conformational exchange rate constant at high temperatures.

Solid-state NMR spectroscopy of proteins

Solid-state NMR spectroscopy proved to be a versatile tool for characterization of structure and dynamics of complex biochemical systems. In particular, magic angle spinning (MAS) solid-state NMR came to maturity for application towards structural elucidation of biological macromolecules. Current challenges in applying solid-state NMR as well as progress achieved recently will be discussed in the following chapter focusing on conceptual aspects important for structural elucidation of proteins.

Protein linewidth and solvent dynamics in frozen solution NMR

Solid-state NMR of proteins in frozen aqueous solution is a potentially powerful technique in structural biology, especially if it is combined with dynamic nuclear polarization signal enhancement strategies. One concern regarding NMR studies of frozen solution protein samples at low temperatures is that they may have poor linewidths, thus preventing high-resolution studies. To learn more about how the solvent shell composition and temperature affects the protein linewidth, we recorded ¹H, ²H, and ¹³C spectra of ubiquitin in frozen water and frozen glycerol-water solutions at different temperatures. We found that the ¹³C protein linewidths generally increase with decreasing temperature. This line broadening was found to be inhomogeneous and independent of proton decoupling. In pure water, we observe an abrupt line broadening with the freezing of the bulk solvent, followed by continuous line broadening at lower temperatures. In frozen glycerol-water, we did not observe an abrupt line broadening and the NMR lines were generally narrower than for pure water at the same temperature. ¹H and ²H measurements characterizing the dynamics of water that is in exchange with the protein showed that the ¹³C line broadening is relatively independent from the arrest of isotropic water motions.

Structural intermediates during α-synuclein fibrillogenesis on phospholipid vesicles

α-Synuclein (AS) fibrils are the main protein component of Lewy bodies, the pathological hallmark of Parkinson's disease and other related disorders. AS forms helices that bind phospholipid membranes with high affinity, but no atomic level data for AS aggregation in the presence of lipids is yet available. Here, we present direct evidence of a conversion from α-helical conformation to β-sheet fibrils in the presence of anionic phospholipid vesicles and direct conversion to β-sheet fibrils in their absence. We have trapped intermediate states throughout the fibril formation pathways to examine the structural changes using solid-state NMR spectroscopy and electron microscopy. The comparison between mature AS fibrils formed in aqueous buffer and those derived in the presence of anionic phospholipids demonstrates no major changes in the overall fibril fold. However, a site-specific comparison of these fibrillar states demonstrates major perturbations in the N-terminal domain with a partial disruption of the long β-strand located in the 40s and small perturbations in residues located in the "non-β amyloid component" (NAC) domain. Combining all these results, we propose a model for AS fibrillogenesis in the presence of phospholipid vesicles.

Slow motions in the hydrophobic core of chicken villin headpiece subdomain and their contributions to configurational entropy and heat capacity from solid-state deuteron NMR measurements

We have investigated microsecond to millisecond time scale dynamics in several key hydrophobic core methyl groups of chicken villin headpiece subdomain protein (HP36) using a combination of single-site labeling, deuteron solid-state NMR line shape analysis, and computational modeling. Deuteron line shapes of hydrated powder samples are dominated by rotameric jumps and show a large variability of rate constants, activation energies, and rotameric populations. Site-specific activation energies vary from 6 to 38 kJ/mol. An additional mode of diffusion on a restricted arc is significant for some sites. In dry samples, the dynamics is quenched. Parameters of the motional models allow for calculations of configurational entropy and heat capacity, which, together with the rate constants, allow for observation of interplay between thermodynamic and kinetic picture of the landscape. Mutations at key phenylalanine residues at both distal (F47L&F51L) and proximal (F58L) locations to a relatively rigid side chain of L69 have a pronounced effect on alleviating the rigidity of this side chain at room temperature and demonstrate the sensitivity of the hydrophobic core environment to such perturbations.

A Simple Mechanism Underlying the Effect of Protecting Osmolytes on Protein Folding

Osmolytes are small organic compounds that confer to the cell an enhanced adaptability to external conditions. Many osmolytes not only protect the cell from osmotic stress but also stabilize the native structure of proteins. While simplified models able to predict changes to protein stability are available, a general physicochemical explanation of the underlying microscopic mechanism is still missing. Here, we address this issue by performing very long all-atom MD simulations, free energy calculations, and experiments on a well-characterized mini-protein, the villin headpiece. Comparisons between the folding free energy landscapes in pure water and osmolyte solutions, together with experimental validation by means of circular dichroism, unfolding experiments, and NMR, led us to formulate a simple hypothesis for the protecting mechanism. Taken together, our results support a novel mechanistic explanation according to which the main driving force behind native state protection is a change in the solvent rotational diffusion.

Conformational disorder of membrane peptides investigated from solid-state NMR line widths and line shapes

A challenge in the application of solid-state NMR spectroscopy to membrane peptides and proteins is the relatively broad line widths compared to those for solution NMR spectra. To understand the linewidth contributions to membrane protein NMR spectra, we have measured the inhomogeneous and homogeneous line widths of several well-studied membrane peptides under immobilized conditions. (13)C T(2) relaxation times of uniformly (13)C-labeled residues show that the homogeneous line widths of the peptides are comparable to those of crystalline model compounds under identical (1)H decoupling and magic angle spinning conditions, indicating that the homogeneous line widths are determined by conformation-independent factors, including residual dipolar coupling, J-coupling, and intrinsic T(2) relaxation. However, the membrane peptides exhibit larger apparent line widths than the crystalline compounds, indicating conformational disorder. A cationic cell-penetrating peptide, the human immunodeficiency virus TAT, exhibits the largest apparent line widths, which are about five-fold larger than the homogeneous line widths, while the transmembrane helix of the influenza M2 peptide and the β-hairpin antimicrobial peptide PG-1 show moderately larger apparent line widths than the crystalline compounds. These results are consistent with the random coil nature of the TAT peptide, which contrasts with the intramolecularly hydrogen bonded M2 and PG-1. Cross peak line shapes of 2D double-quantum correlation spectra show that the conformational disorder can occur at the residue level and can result from three origins, lipid-peptide interaction, intrinsic conformational disorder encoded in the amino acid sequence, and side-chain rotameric averaging. A particularly important lipid-peptide interaction for cationic membrane peptides is guanidinium-phosphate ion pair interaction. Thus, NMR line widths and line shapes are useful for understanding the conformational disorder of membrane peptides and proteins.

Bridging the Gap between Folding Simulations and Experiments: The Case of the Villin Headpiece

The increasing accuracy of molecular dynamics force fields parameters and the increasing resolution of experimental results allow one to carefully compare and complement in silico data with experimental observations. Here, we study the human villin headpiece C-terminal helical subdomain (HP35) with the recent highly optimized Amber99SB*-ILDN force field and compare the results with recent high resolution triplet-triplet energy transfer (TTET) experiments. The correct reproduction of the main structural features reveals a good agreement between experimental data and simulations.

Temperature dependence of fast carbonyl backbone dynamics in chicken villin headpiece subdomain

Temperature-dependence of protein dynamics can provide information on details of the free energy landscape by probing the characteristics of the potential responsible for the fluctuations. We have investigated the temperature-dependence of picosecond to nanosecond backbone dynamics at carbonyl carbon sites in chicken villin headpiece subdomain protein using a combination of three NMR relaxation rates: (13)C' longitudinal rate, and two cross-correlated rates involving dipolar and chemical shift anisotropy (CSA) relaxation mechanisms, (13)C'/(13)C'-(13)C(α) CSA/dipolar and (13)C'/(13)C'-(15)N CSA/dipolar. Order parameters have been extracted using the Lipari-Szabo model-free approach assuming a separation of the time scales of internal and molecular motions in the 2-16°C temperature range. There is a gradual deviation from this assumption from lower to higher temperatures, such that above 16°C the separation of the time scales is inconsistent with the experimental data and, thus, the Lipari-Szabo formalism can not be applied. While there are variations among the residues, on the average the order parameters indicate a markedly steeper temperature dependence at backbone carbonyl carbons compared to that probed at amide nitrogens in an earlier study. This strongly advocates for probing sites other than amide nitrogen for accurate characterization of the potential and other thermodynamics characteristics of protein backbone.

Kinetic analysis of protein aggregation monitored by real-time 2D solid-state NMR spectroscopy

It is shown that real-time 2D solid-state NMR can be used to obtain kinetic and structural information about the process of protein aggregation. In addition to the incorporation of kinetic information involving intermediate states, this approach can offer atom-specific resolution for all detectable species. The analysis was carried out using experimental data obtained during aggregation of the 10.4 kDa Crh protein, which has been shown to involve a partially unfolded intermediate state prior to aggregation. Based on a single real-time 2D (13)C-(13)C transition spectrum, kinetic information about the refolding and aggregation step could be extracted. In addition, structural rearrangements associated with refolding are estimated and several different aggregation scenarios were compared to the experimental data.

Hydrophobic core formation and dehydration in protein folding studied by generalized-ensemble simulations

Despite its small size, chicken villin headpiece subdomain HP36 folds into the native structure with a stable hydrophobic core within several microseconds. How such a small protein keeps up its conformational stability and fast folding in solution is an important issue for understanding molecular mechanisms of protein folding. In this study, we performed multicanonical replica-exchange simulations of HP36 in explicit water, starting from a fully extended conformation. We observed at least five events of HP36 folding into nativelike conformations. The smallest backbone root mean-square deviation from the crystal structure was 1.1 A. In the nativelike conformations, the stably formed hydrophobic core was fully dehydrated. Statistical analyses of the simulation trajectories show the following sequential events in folding of HP36: 1), Helix 3 is formed at the earliest stage; 2), the backbone and the side chains near the loop between Helices 2 and 3 take nativelike conformations; and 3), the side-chain packing at the hydrophobic core and the dehydration of the core side chains take place simultaneously at the later stage of folding. This sequence suggests that the initial folding nucleus is not necessarily the same as the hydrophobic core, consistent with a recent experimental phi-value analysis.

Protein-ice interaction of an antifreeze protein observed with solid-state NMR

NMR on frozen solutions is an ideal method to study fundamental questions of macromolecular hydration, because the hydration shell of many biomolecules does not freeze together with bulk solvent. In the present study, we present previously undescribed NMR methods to study the interactions of proteins with their hydration shell and the ice lattice in frozen solution. We applied these methods to compare solvent interaction of an ice-binding type III antifreeze protein (AFP III) and ubiquitin a non-ice-binding protein in frozen solution. We measured (1)H-(1)H cross-saturation and cross-relaxation to provide evidence for a molecular contact surface between ice and AFP III at moderate freezing temperatures of -35 °C. This phenomenon is potentially unique for AFPs because ubiquitin shows no such cross relaxation or cross saturation with ice. On the other hand, we detected liquid hydration water and strong water-AFP III and water-ubiquitin cross peaks in frozen solution using relaxation filtered (2)H and HETCOR spectra with additional (1)H-(1)H mixing. These results are consistent with the idea that ubiquitin is surrounded by a hydration shell, which separates it from the bulk ice. For AFP III, the water cross peaks indicate that only a portion of its hydration shell (i.e., at the ice-binding surface) is in contact with the ice lattice. The rest of AFP III's hydration shell behaves similarly to the hydration shell of non-ice-interacting proteins such as ubiquitin and does not freeze together with the bulk water.

Conformational plasticity of the influenza A M2 transmembrane helix in lipid bilayers under varying pH, drug binding, and membrane thickness

Membrane proteins change their conformations to respond to environmental cues, thus conformational plasticity is important for function. The influenza A M2 protein forms an acid-activated proton channel important for the virus lifecycle. Here we have used solid-state NMR spectroscopy to examine the conformational plasticity of membrane-bound transmembrane domain of M2 (M2TM). (13)C and (15)N chemical shifts indicate coupled conformational changes of several pore-facing residues due to changes in bilayer thickness, drug binding, and pH. The structural changes are attributed to the formation of a well-defined helical kink at G34 in the drug-bound state and in thick lipid bilayers, nonideal backbone conformation of the secondary-gate residue V27 in the presence of drug, and nonideal conformation of the proton-sensing residue H37 at high pH. The chemical shifts constrained the (ϕ, ψ) torsion angles for three "basis" states, the equilibrium among which explains the multiple resonances per site in the NMR spectra under different combinations of bilayer thickness, drug binding, and pH conditions. Thus, conformational plasticity is important for the proton conduction and inhibition of M2TM. The study illustrates the utility of NMR chemical shifts for probing the structural plasticity and folding of membrane proteins.

Modelling proteins: conformational sampling and reconstruction of folding kinetics

In the last decades biomolecular simulation has made tremendous inroads to help elucidate biomolecular processes in-silico. Despite enormous advances in molecular dynamics techniques and the available computational power, many problems involve long time scales and large-scale molecular rearrangements that are still difficult to sample adequately. In this review we therefore summarise recent efforts to fundamentally improve this situation by decoupling the sampling of the energy landscape from the description of the kinetics of the process. Recent years have seen the emergence of many advanced sampling techniques, which permit efficient characterisation of the relevant family of molecular conformations by dispensing with the details of the short-term kinetics of the process. Because these methods generate thermodynamic information at best, they must be complemented by techniques to reconstruct the kinetics of the process using the ensemble of relevant conformations. Here we review recent advances for both types of methods and discuss their perspectives to permit efficient and accurate modelling of large-scale conformational changes in biomolecules. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.