Paramagnetic-based NMR restraints lift residual dipolar coupling degeneracy in multidomain detergent-solubilized membrane proteins

Residual dipolar couplings (RDCs) are widely used as orientation-dependent NMR restraints to improve the resolution of the NMR conformational ensemble of biomacromolecules and define the relative orientation of multidomain proteins and protein complexes. However, the interpretation of RDCs is complicated by the intrinsic degeneracy of analytical solutions and protein dynamics that lead to ill-defined orientations of the structural domains (ghost orientations). Here, we illustrate how restraints from paramagnetic relaxation enhancement (PRE) experiments lift the orientational ambiguity of multidomain membrane proteins solubilized in detergent micelles. We tested this approach on monomeric phospholamban (PLN), a 52-residue membrane protein, which is composed of two helical domains connected by a flexible loop. We show that the combination of classical solution NMR restraints (NOEs and dihedral angles) with RDC and PRE constraints resolves topological ambiguities, improving the convergence of the PLN structural ensemble and giving the depth of insertion of the protein within the micelle. The combination of RDCs with PREs will be necessary for improving the accuracy and precision of membrane protein conformational ensembles, where three-dimensional structures are dictated by interactions with the membrane-mimicking environment rather than compact tertiary folds common in globular proteins.     

Conformational Ensemble of Disordered Proteins Probed by Solvent Paramagnetic Relaxation Enhancement (sPRE)

Characterization of the conformational ensemble of disordered proteins is highly important for understanding protein folding and aggregation mechanisms, but remains a computational and experimental challenge owing to the dynamic nature of these proteins. New observables that can provide unique insights into transient residual structures in disordered proteins are needed. Here using denatured ubiquitin as a model system, NMR solvent paramagnetic relaxation enhancement (sPRE) measurements provide an accurate and highly sensitive probe for detecting low populations of residual structure in a disordered protein. Furthermore, a new ensemble calculation approach based on sPRE restraints in conjunction with residual dipolar couplings (RDCs) and small‐angle X‐ray scattering (SAXS) is used to define the conformational ensemble of disordered proteins at atomic resolution. The approach presented should be applicable to a wide range of dynamic macromolecules.     

Calculation of residual dipolar couplings from disordered state ensembles using local alignment

Residual dipolar couplings (RDCs) have been observed in disordered states of several proteins. While their nonuniform values were initially surprising, it has been shown that reasonable approximation of experimental RDCs can be obtained using simple statistical coil models and assuming global alignment of each structure, provided that many thousands of conformers are averaged. Here we show that, by using short local alignment tensors, we can achieve good agreement between experimental and simulated RDCs with far fewer structures than required when using global alignment. This makes the possibility of using RDCs as direct restraints in structural calculations of disordered proteins much more feasible. In addition, it provides insight into the nature of RDCs in disordered states, suggesting that they are primarily reporting on local structure.     

Prediction of protein disorder at the domain level

Intrinsically disordered/unstructured proteins exist in a highly flexible conformational state largely devoid of secondary structural elements and tertiary contacts. Despite their lack of a well defined structure, these proteins often fulfill essential regulatory functions. The intrinsic lack of structure confers functional advantages on these proteins, allowing them to adopt multiple conformations and to bind to different binding partners. The structural flexibility of disordered regions hampers efforts solving structures at high resolution by X-ray crystallography and/or NMR. Removing such proteins/regions from high-throughput structural genomics pipelines would be of significant benefit in terms of cost and success rate. In this paper we outline the theoretical background of structural disorder, and review bioinformatic predictors that can be used to delineate regions most likely to be amenable for structure determination. The primary focus of our review is the interpretation of prediction results in a way that enables segmentation of proteins to separate ordered domains from disordered regions.     

Highly populated turn conformations in natively unfolded tau protein identified from residual dipolar couplings and molecular simulation

Tau, a natively unstructured protein that regulates the organization of neuronal microtubules, is also found in high concentrations in neurofibrillary tangles of Alzheimer's disease and other neurodegenerative disorders. The conformational transition between these vastly different healthy and pathological forms remains poorly understood. We have measured residual dipolar couplings (RDCs), J-couplings, and nuclear Overhauser enhancement (NOE) in construct K18 of tau, containing all four repeat domains R1-R4. NHN RDCs were compared with prediction on the basis of a statistical model describing the intrinsic conformational sampling of unfolded proteins in solution. While local variation and relative amplitude of RDCs agrees with propensity-based prediction for most of the protein, homologous sequences in each repeat domain (DLKN, DLSN, DLSK, and DKFD in repeats R1-R4) show strong disagreement characterized by inversion of the sign of the central couplings. Accelerated molecular dynamic simulations (AMD) in explicit solvent revealed strong tendencies to form turns, identified as type I beta-turns for repeats R1-R3. Incorporation of the backbone dihedral sampling resulting from AMD into the statistical coil model closely reproduces experimental RDC values. These localized sequence-dependent conformational tendencies interrupt the propensity to sample more extended conformations in adjacent strands and are remarkably resistant to local environmental factors, as demonstrated by the persistence of the RDC signature even under harsh denaturing conditions (8 M urea). The role that this specific conformational behavior may play in the transition to the pathological form is discussed.     

Structure and dynamics of an unfolded protein examined by molecular dynamics simulation

The accurate characterization of the structure and dynamics of proteins in disordered states is a difficult problem at the frontier of structural biology whose solution promises to further our understanding of protein folding and intrinsically disordered proteins. Molecular dynamics (MD) simulations have added considerably to our understanding of folded proteins, but the accuracy with which the force fields used in such simulations can describe disordered proteins is unclear. In this work, using a modern force field, we performed a 200 μs unrestrained MD simulation of the acid-unfolded state of an experimentally well-characterized protein, ACBP, to explore the extent to which state-of-the-art simulation can describe the structural and dynamical features of a disordered protein. By comparing the simulation results with the results of NMR experiments, we demonstrate that the simulation successfully captures important aspects of both the local and global structure. Our simulation was ~2 orders of magnitude longer than those in previous studies of unfolded proteins, a length sufficient to observe repeated formation and breaking of helical structure, which we found to occur on a multimicrosecond time scale. We observed one structural feature that formed but did not break during the simulation, highlighting the difficulty in sampling disordered states. Overall, however, our simulation results are in reasonable agreement with the experimental data, demonstrating that MD simulations can already be useful in describing disordered proteins. Finally, our direct calculation of certain NMR observables from the simulation provides new insight into the general relationship between structural features of disordered proteins and experimental NMR relaxation properties.     

Disorder-to-order transition of an intrinsically disordered region of sortase revealed by multiscale enhanced sampling

Molecular functions of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs), such as molecular recognition and cellular signaling, are ascribed to dynamic changes in the conformational space in response to binding of target molecules. Sortase, a transpeptitase in Gram-positive bacteria, has an IDR in a loop which undergoes a disordered-to-ordered transition (called "disordered loop"), accompanying a tilt of another loop ("dynamic loop"), upon binding of a signal peptide and a calcium ion. In this study, all-atom conformational ensembles of sortase were calculated for the four different binding states (with/without the peptide and with/without a calcium ion) by the multiscale enhanced sampling (MSES) simulation to examine how the binding of the peptide and/or calcium influences the conformational ensemble. The MSES is a multiscale and multicopy simulation method that allows an enhanced sampling of the all-atom model of large proteins including explicit solvent. A 100 ns MSES simulation of the ligand-free sortase using 20 replicas (in total 2 μs) demonstrated large flexibility in both the disordered and dynamic loops; however, their distributions were not random but had a clear preference which populates the N-terminal part of the disordered loop near the bound form. The MSES simulations of the three binding states clarified the allosteric mechanism of sortase: the N- and C-terminal parts of the disordered loop undergo a disorder-to-order transition independently of each other upon binding of the peptide and a calcium ion, respectively; however, upon binding of both ligands, the two parts work cooperatively to stabilize the bound peptide.     

The use of residual dipolar coupling for conformational analysis of structurally related natural products

Determining the conformational preferences of molecules in solution remains a considerable challenge. Recently, the use of residual dipolar coupling (RDC) analysis has emerged as a key method to address this. Whilst to date the majority of the applications have focused on biomolecules including proteins and DNA, the use of RDCs for studying small molecules is gaining popularity. Having said that, the method continues to develop, and here, we describe an early case study of the quantification of conformer populations in small molecules using RDC analysis. Having been inspired to study conformational preferences by unexpected differences in the NMR spectra and the reactivity of related natural products, we showed that the use of more established techniques was unsatisfactory in explaining the experimental observations. The use of RDCs provided an improved understanding that, following use of methods to quantify conformer populations using RDCs, culminated in a rationalisation of the contrasting diastereoselectivities observed in a ketone reduction reaction. Copyright © 2015 John Wiley & Sons, Ltd.     

Searching and optimizing structure ensembles for complex flexible sugars

NMR restrictions are suitable to specify the geometry of a molecule when a single well-defined global free energy minimum exists that is significantly lower than other local minima. Carbohydrates are quite flexible, and therefore, NMR observables do not always correlate with a single conformer but instead with an ensemble of low free energy conformers that can be accessed by thermal fluctuations. In this communication, we describe a novel procedure to identify and weight the contribution to the ensemble of local minima conformers based on comparison to residual dipolar couplings (RDCs) or other NMR observables, such as scalar couplings. A genetic algorithm is implemented to globally minimize the R factor comparing calculated RDCs to experiment. This is done by optimizing the weights of different conformers derived from the exhaustive local minima conformational search program, fast sugar structure prediction software (FSPS). We apply this framework to six human milk sugars, LND-1, LNF-1, LNF-2, LNF-3, LNnT, and LNT, and are able to determine corresponding population weights for the ensemble of conformers. Interestingly, our results indicate that in all cases the RDCs can be well represented by only a few most important conformers. This confirms that several, but not all of the glycosidic linkages in histo-blood group "epitopes" are quite rigid.     

Unraveling the nature of the segmentation clock: Intrinsic disorder of clock proteins and their interaction map

Vertebrate segmentation has been proved to be under a strict temporal control governed by a biological clock, known as the segmentation clock. The present experimental evidence suggests that the segmentation clock initiates and maintains its periodic cycle by the periodic activation or inhibition of the Notch signaling pathway as well as the periodic autoregulation of the cyclic genes themselves. In this paper, we investigate the structural and evolutionary properties of the Notch pathway proteins involved in the mice segmentation clock and computationally identify the interaction map within the Notch signaling pathway. The results of our analysis strongly indicate that most of the pathway proteins are intrinsically disordered and that the mechanism of their interaction likely involves helical molecular recognition elements, short loosely structured segments within disordered regions which are directly involved in protein-protein interactions. Predicted interactions are in agreement with gene knock-out studies available in the literature.     

Replica exchange with guided annealing for accelerated sampling of disordered protein conformations

We critically examine a recently proposed convective replica exchange (cRE) method for enhanced sampling of protein conformation based on theoretical and numerical analysis. The results demonstrate that cRE and related replica exchange with guided annealing (RE‐GA) schemes lead to unbalanced exchange attempt probabilities and break detailed balance whenever the system undergoes slow conformational transitions (relative to the temperature diffusion timescale). Nonetheless, numerical simulations suggest that approximate canonical ensembles can be generated for systems with small conformational transition barriers. This suggests that RE‐GA maybe suitable for simulating intrinsically disordered proteins, an important class of newly recognized functional proteins. The efficacy of RE‐GA is demonstrated by calculating the conformational ensembles of intrinsically disordered kinase inducible domain protein. The results show that RE‐GA helps the protein to escape nonspecific compact states more efficiently and provides several fold speedups in generating converged and largely correct ensembles compared to the standard temperature RE. © 2014 Wiley Periodicals, Inc.     

Magnetic Alignment of Polymer Macro‐Nanodiscs Enables Residual‐Dipolar‐Coupling‐Based High‐Resolution Structural Studies by NMR Spectroscopy

Experimentally measured residual dipolar couplings (RDCs) are highly valuable for atomic‐resolution structural and dynamic studies of molecular systems ranging from small molecules to large proteins by solution NMR spectroscopy. Here we demonstrate the first use of magnetic‐alignment behavior of lyotropic liquid‐crystalline polymer macro‐nanodiscs (>20 nm in diameter) as a novel alignment medium for the measurement of RDCs using high‐resolution NMR. The easy preparation of macro‐nanodiscs, their high stability against pH changes and the presence of divalent metal ions, and their high homogeneity make them an efficient tool to investigate a wide range of molecular systems including natural products, proteins, and RNA.     

Intrinsically Disordered Tardigrade Proteins Self-Assemble into Fibrous Gels in Response to Environmental Stress

Tardigrades are remarkable for their ability to survive harsh stress conditions as diverse as extreme temperature and desiccation. The molecular mechanisms that confer this unusual resistance to physical stress remain unknown. Recently, tardigrade-unique intrinsically disordered proteins have been shown to play an essential role in tardigrade anhydrobiosis. Here, we characterize the conformational and physical behaviour of CAHS-8 from Hypsibius exemplaris. NMR spectroscopy reveals that the protein comprises an extended central helical domain flanked by disordered termini. Upon concentration, the protein is shown to successively form oligomers, long fibres, and finally gels constituted of fibres in a strongly temperature-dependent manner. The helical domain forms the core of the fibrillar structure, with the disordered termini remaining highly dynamic within the gel. Soluble proteins can be encapsulated within cavities in the gel, maintaining their functional form. The ability to reversibly form fibrous gels may be associated with the enhanced protective properties of these proteins.     

Elongated Bacterial Pili as a Versatile Alignment Medium for NMR Spectroscopy

In NMR spectroscopy, residual dipolar couplings (RDCs) have emerged as one of the most exquisite probes of biological structure and dynamics. The measurement of RDCs relies on the partial alignment of the molecule of interest, for example by using a liquid crystal as a solvent. Here, we establish bacterial type 1 pili as an alternative liquid-crystalline alignment medium for the measurement of RDCs. To achieve alignment at pilus concentrations that allow for efficient NMR sample preparation, we elongated wild-type pili by recombinant overproduction of the main structural pilus subunit. Building on the extraordinary stability of type 1 pili against spontaneous dissociation and unfolding, we show that the medium is compatible with challenging experimental conditions such as high temperature, the presence of detergents, organic solvents or very acidic pH, setting it apart from most established alignment media. Using human ubiquitin, HIV-1 TAR RNA and camphor as spectroscopic probes, we demonstrate the applicability of the medium for the determination of RDCs of proteins, nucleic acids and small molecules. Our results show that type 1 pili represent a very useful alternative to existing alignment media and may readily assist the characterization of molecular structure and dynamics by NMR.     

Computational strategy for intrinsically disordered protein ligand design leads to the discovery of p53 transactivation domain I binding compounds that activate the p53 pathway

Intrinsically disordered proteins or intrinsically disordered regions (IDPs) have gained much attention in recent years due to their vital roles in biology and prevalence in various human diseases. Although IDPs are perceived as attractive therapeutic targets, rational drug design targeting IDPs remains challenging because of their conformational heterogeneity. Here, we propose a hierarchical computational strategy for IDP drug virtual screening (IDPDVS) and applied it in the discovery of p53 transactivation domain I (TAD1) binding compounds. IDPDVS starts from conformation sampling of the IDP target, then it combines stepwise conformational clustering with druggability evaluation to identify potential ligand binding pockets, followed by multiple docking screening runs and selection of compounds that can bind multi-conformations. p53 is an important tumor suppressor and restoration of its function provides an opportunity to inhibit cancer cell growth. TAD1 locates at the N-terminus of p53 and plays key roles in regulating p53 function. No compounds that directly bind to TAD1 have been reported due to its highly disordered structure. We successfully used IDPDVS to identify two compounds that bind p53 TAD1 and restore wild-type p53 function in cancer cells. Our study demonstrates that IDPDVS is an efficient strategy for IDP drug discovery and p53 TAD1 can be directly targeted by small molecules.

A hierarchical computational strategy for IDP drug virtual screening (IDPDVS) was proposed and successfully applied to identify compounds that bind p53 TAD1 and restore wild-type p53 function in cancer cells.     

Efficacy of independence sampling in replica exchange simulations of ordered and disordered proteins

Recasting temperature replica exchange (T‐RE) as a special case of Gibbs sampling has led to a simple and efficient scheme for enhanced mixing (Chodera and Shirts, J. Chem. Phys., 2011, 135, 194110). To critically examine if T‐RE with independence sampling (T‐REis) improves conformational sampling, we performed T‐RE and T‐REis simulations of ordered and disordered proteins using coarse‐grained and atomistic models. The results demonstrate that T‐REis effectively increase the replica mobility in temperatures space with minimal computational overhead, especially for folded proteins. However, enhanced mixing does not translate well into improved conformational sampling. The convergences of thermodynamic properties interested are similar, with slight improvements for T‐REis of ordered systems. The study re‐affirms the efficiency of T‐RE does not appear to be limited by temperature diffusion, but by the inherent rates of spontaneous large‐scale conformational re‐arrangements. Due to its simplicity and efficacy of enhanced mixing, T‐REis is expected to be more effective when incorporated with various Hamiltonian‐RE protocols. © 2017 Wiley Periodicals, Inc.     

Using Cross-Correlated Spin Relaxation to Characterize Backbone Dihedral Angle Distributions of Flexible Protein Segments

Crucial to the function of proteins is their existence as conformational ensembles sampling numerous and structurally diverse substates. Despite this widely accepted notion there is still a high demand for meaningful and reliable approaches to characterize protein ensembles in solution. As it is usually conducted in solution, NMR spectroscopy offers unique possibilities to address this challenge. Particularly, cross-correlated relaxation (CCR) effects have long been established to encode both protein structure and dynamics in a compelling manner. However, this wealth of information often limits their use in practice as structure and dynamics might prove difficult to disentangle. Using a modern Maximum Entropy (MaxEnt) reweighting approach to interpret CCR rates of Ubiquitin, we demonstrate that these uncertainties do not necessarily impair resolving CCR-encoded structural information. Instead, a suitable balance between complementary CCR experiments and prior information is found to be the most crucial factor in mapping backbone dihedral angle distributions. Experimental and systematic deviations such as oversimplified dynamics appear to be of minor importance. Using Ubiquitin as an example, we demonstrate that CCR rates are capable of characterizing rigid and flexible residues alike, indicating their unharnessed potential in studying disordered proteins.