How random are intrinsically disordered proteins? A small angle scattering perspective

While the crucial role of intrinsically disordered proteins (IDPs) in the cell cycle is now recognized, deciphering their molecular mode of action at the structural level still remains highly challenging and requires a combination of many biophysical approaches. Among them, small angle X-ray scattering (SAXS) has been extremely successful in the last decade and has become an indispensable technique for addressing many of the fundamental questions regarding the activities of IDPs. After introducing some experimental issues specific to IDPs and in relation to the latest technical developments, this article presents the interest of the theory of polymer physics to evaluate the flexibility of fully disordered proteins. The different strategies to obtain 3-dimensional models of IDPs, free in solution and associated in a complex, are then reviewed. Indeed, recent computational advances have made it possible to readily extract maximum information from the scattering curve with a special emphasis on highly flexible systems, such as multidomain proteins and IDPs. Furthermore, integrated computational approaches now enable the generation of ensembles of conformers to translate the unique flexible characteristics of IDPs by taking into consideration the constraints of more and more various complementary experiment. In particular, a combination of SAXS with high-resolution techniques, such as x-ray crystallography and NMR, allows us to provide reliable models and to gain unique structural insights about the protein over multiple structural scales. The latest neutron scattering experiments also promise new advances in the study of the conformational changes of macromolecules involving more complex systems.     

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.     

Quantifying polypeptide conformational space: sensitivity to conformation and ensemble definition

Quantifying the density of conformations over phase space (the conformational distribution) is needed to model important macromolecular processes such as protein folding. In this work, we quantify the conformational distribution for a simple polypeptide (N-mer polyalanine) using the cumulative distribution function (CDF), which gives the probability that two randomly selected conformations are separated by less than a "conformational" distance and whose inverse gives conformation counts as a function of conformational radius. An important finding is that the conformation counts obtained by the CDF inverse depend critically on the assignment of a conformation's distance span and the ensemble (e.g., unfolded state model): varying ensemble and conformation definition (1 --> 2 A) varies the CDF-based conformation counts for Ala(50) from 10(11) to 10(69). In particular, relatively short molecular dynamics (MD) relaxation of Ala(50)'s random-walk ensemble reduces the number of conformers from 10(55) to 10(14) (using a 1 A root-mean-square-deviation radius conformation definition) pointing to potential disconnections in comparing the results from simplified models of unfolded proteins with those from all-atom MD simulations. Explicit waters are found to roughen the landscape considerably. Under some common conformation definitions, the results herein provide (i) an upper limit to the number of accessible conformations that compose unfolded states of proteins, (ii) the optimal clustering radius/conformation radius for counting conformations for a given energy and solvent model, (iii) a means of comparing various studies, and (iv) an assessment of the applicability of random search in protein folding.     

Segmental,Domain‐Selective Perdeuteration and Small‐Angle Neutron Scattering for Structural Analysis of Multi‐Domain Proteins

Multi‐domain proteins play critical roles in fine‐tuning essential processes in cellular signaling and gene regulation. Typically, multiple globular domains that are connected by flexible linkers undergo dynamic rearrangements upon binding to protein, DNA or RNA ligands. RNA binding proteins (RBPs) represent an important class of multi‐domain proteins, which regulate gene expression by recognizing linear or structured RNA sequence motifs. Here, we employ segmental perdeuteration of the three RNA recognition motif (RRM) domains in the RBP TIA‐1 using Sortase A mediated protein ligation. We show that domain‐selective perdeuteration combined with contrast‐matched small‐angle neutron scattering (SANS), SAXS and computational modeling provides valuable information to precisely define relative domain arrangements. The approach is generally applicable to study conformational arrangements of individual domains in multi‐domain proteins and changes induced by ligand binding.     

Segmental,Domain‐Selective Perdeuteration and Small‐Angle Neutron Scattering for Structural Analysis of Multi‐Domain Proteins

Multi‐domain proteins play critical roles in fine‐tuning essential processes in cellular signaling and gene regulation. Typically, multiple globular domains that are connected by flexible linkers undergo dynamic rearrangements upon binding to protein, DNA or RNA ligands. RNA binding proteins (RBPs) represent an important class of multi‐domain proteins, which regulate gene expression by recognizing linear or structured RNA sequence motifs. Here, we employ segmental perdeuteration of the three RNA recognition motif (RRM) domains in the RBP TIA‐1 using Sortase A mediated protein ligation. We show that domain‐selective perdeuteration combined with contrast‐matched small‐angle neutron scattering (SANS), SAXS and computational modeling provides valuable information to precisely define relative domain arrangements. The approach is generally applicable to study conformational arrangements of individual domains in multi‐domain proteins and changes induced by ligand binding.     

Relationship among ligand conformations in solution, in the solid state, and at the Hsp90 binding site: geldanamycin and radicicol

The unknown effects of a receptor's environment on a ligand's conformation presents a difficult challenge in predicting feasible bioactive conformations, particularly if the receptor is ill-defined. The primary hypothesis of this work is that a structure's conformational ensemble in solution presents viable candidates for protein binding. The experimental solution profile can be achieved with the NAMFIS (NMR analysis of molecular flexibility in solution) method, which deconvolutes the average NMR spectrum of small flexible molecules into individual contributing conformations with varying populations. Geldanamycin and radicicol are structurally different macrocycles determined by X-ray crystallography to bind to a common site on the cellular chaperone heat shock protein 90 (Hsp90). Without benefit of a receptor structure, NAMFIS has identified the bioactive conformers of geldanamycin and radicicol in CDCl3 solution with populations of 4% and 21%, respectively. Conversely, docking the set of NAMFIS conformers into the unliganded proteins with GLIDE followed by MM-GBSA scoring reproduces the experimental crystallographic binding poses.     

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.     

Refinement of multidomain protein structures by combination of solution small-angle X-ray scattering and NMR data

Determination of the 3D structures of multidomain proteins by solution NMR methods presents a number of unique challenges related to their larger molecular size and the usual scarcity of constraints at the interdomain interface, often resulting in a decrease in structural accuracy. In this respect, experimental information from small-angle scattering of X-ray radiation in solution (SAXS) presents a suitable complement to the NMR data, as it provides an independent constraint on the overall molecular shape. A computational procedure is described that allows incorporation of such SAXS data into the mainstream high-resolution macromolecular structure refinement. The method is illustrated for a two-domain 177-amino-acid protein, gammaS crystallin, using an experimental SAXS data set fitted at resolutions from approximately 200 A to approximately 30 A. Inclusion of these data during structure refinement decreases the backbone coordinate root-mean-square difference between the derived model and the high-resolution crystal structure of a 54% homologous gammaB crystallin from 1.96 +/- 0.07 A to 1.31 +/- 0.04 A. Combining SAXS data with NMR restraints can be accomplished at a moderate computational expense and is expected to become useful for multidomain proteins, multimeric assemblies, and tight macromolecular complexes.     

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.     

Probing dynamic protein ensembles with atomic proximity measures

The emerging role of internal dynamics in protein fold and function requires new avenues of structure analysis. We analyzed the dynamically restrained conformational ensemble of ubiquitin generated from residual dipolar coupling data, in terms of protruding and buried atoms as well as interatomic distances, using four proximity-based algorithms, CX, DPX, PRIDE and PRIDE-NMR (http://hydra.icgeb.trieste.it/protein/). We found that Ubiquitin, this relatively rigid molecule has a highly diverse dynamic ensemble. The environment of protruding atoms is highly variable across conformers, on the other hand, only a part of buried atoms tends to fluctuate. The variability of the ensemble cautions against the use of single conformers when explaining functional phenomena. We also give a detailed evaluation of PRIDE-NMR on a wide dataset and discuss its usage in the light of the features of available NMR distance restraint sets in public databases.     

Novel inhibitor discovery through virtual screening against multiple protein conformations generated via ligand-directed modeling: a maternal embryonic leucine zipper kinase example

Kinase targets have been demonstrated to undergo major conformational reorganization upon ligand binding. Such protein conformational plasticity remains a significant challenge in structure-based virtual screening methodology and may be approximated by screening against an ensemble of diverse protein conformations. Maternal embryonic leucine zipper kinase (MELK), a member of serine-threonine kinase family, has been recently found to be involved in the tumerogenic state of glioblastoma, breast, ovarian, and colon cancers. We therefore modeled several conformers of MELK utilizing the available chemogenomic and crystallographic data of homologous kinases. We carried out docking pose prediction and virtual screening enrichment studies with these conformers. The performances of the ensembles were evaluated by their ability to reproduce known inhibitor bioactive conformations and to efficiently recover known active compounds early in the virtual screen when seeded with decoy sets. A few of the individual MELK conformers performed satisfactorily in reproducing the native protein-ligand pharmacophoric interactions up to 50% of the cases. By selecting an ensemble of a few representative conformational states, most of the known inhibitor binding poses could be rationalized. For example, a four conformer ensemble is able to recover 95% of the studied actives, especially with imperfect scoring function(s). The virtual screening enrichment varied considerably among different MELK conformers. Enrichment appears to improve by selection of a proper protein conformation. For example, several holo and unliganded active conformations are better to accommodate diverse chemotypes than ATP-bound conformer. These results prove that using an ensemble of diverse conformations could give a better performance. Applying this approach, we were able to screen a commercially available library of half a million compounds against three conformers to discover three novel inhibitors of MELK, one from each template. Among the three compounds validated via experimental enzyme inhibition assays, one is relatively potent (15; K(d) = 0.37 μM), one moderately active (12; K(d) = 3.2 μM), and one weak but very selective (9; K(d) = 18 μM). These novel hits may be utilized to assist in the development of small molecule therapeutic agents useful in diseases caused by deregulated MELK, and perhaps more importantly, the approach demonstrates the advantages of choosing an appropriate ensemble of a few conformers in pursuing compound potency, selectivity, and novel chemotypes over using single target conformation for structure-based drug design in general.     

Sulfonated molecules that bind a partially structured species of beta2-microglobulin also influence refolding and fibrillogenesis

Human beta2-microglobulin (beta2-m) is a small amyloidogenic protein responsible for dialysis-related amyloidosis, which represents a severe complication of long-term hemodialysis. A therapeutic approach for this amyloidosis could be based on the stabilization of beta2-m through the binding to a small molecule, to possibly inhibit protein misfolding and amyloid fibril formation. The search of a strong ligand of this protein is extremely challenging: by using CE in affinity and refolding experiments we study the effect that previously selected sulfonated molecules have on the equilibrium between the native form and an ensemble of conformers populating the slow phase of beta2-m folding. These data are correlated with the effect that the same molecules exert on in vitro fibrillogenesis experiments.     

Effect of Linker Length and Flexibility on the <Emphasis Type="Italic">Clostridium thermocellum</Emphasis> Esterase Displayed on <Emphasis Type="Italic">Bacillus subtilis</Emphasis> Spores

In fusion protein design strategies, the flexibility and length of linkers are important parameters affecting the bioactivity of multifunctional proteins. A series of fusion proteins with different linkers were constructed. The effect of temperature, pH, and organic solvents was investigated on the enzymatic activity. Fusion proteins with P1(PTPTPT) and P2((PTPTPT)₂) linkers remained highly active with wide temperature range. At pH 9.6, the relative activity of fusion proteins with (PTPTPT)₂ and S2(EGKSSGSGSESKST) linkers was 70 and 62 % (1.75 and 1.5 times of that of non-linker ones). Fusion proteins with S3((GGGGS)₄) linker retained 55 % activity after 5 h of incubation at 80 °C (1.2-fold of that of non-linker fusion proteins and 1.9-fold of GGGGS-linker fusion proteins). Finally, the relative activity of fusion proteins having different linkers was increased with 20 % dimethyl sulfoxide (DMSO) and methanol; relative activity of fusion proteins with EGKSSGSGSESKST linkers was enhanced 1.5- and 2.2-fold, respectively. These results suggest that longer flexible linker can enhance the activity and stability of displayed esterase than shorter flexible linker. Optimizing peptide linkers with length, flexibility, and amino acid composition could improve the thermostability and activity of the displayed enzyme.     

Dependence of effective molarity on linker length for an intramolecular protein-ligand system

This paper reports dissociation constants and "effective molarities" (M(eff)) for the intramolecular binding of a ligand covalently attached to the surface of a protein by oligo(ethylene glycol) (EG(n)) linkers of different lengths (n = 0, 2, 5, 10, and 20) and compares these experimental values with theoretical estimates from polymer theory. As expected, the value of M(eff) is lowest when the linker is too short (n = 0) to allow the ligand to bind noncovalently at the active site of the protein without strain, is highest when the linker is the optimal length (n = 2) to allow such binding to occur, and decreases monotonically as the length increases past this optimal value (but only by a factor of approximately 8 from n = 2 to n = 20). These experimental results are not compatible with a model in which the single bonds of the linker are completely restricted when the ligand has bound noncovalently to the active site of the protein, but they are quantitatively compatible with a model that treats the linker as a random-coil polymer. Calorimetry revealed that enthalpic interactions between the linker and the protein are not important in determining the thermodynamics of the system. Taken together, these results suggest that the manifestation of the linker in the thermodynamics of binding is exclusively entropic. The values of M(eff) are, theoretically, intrinsic properties of the EG(n) linkers and can be used to predict the avidities of multivalent ligands with these linkers for multivalent proteins. The weak dependence of M(eff) on linker length suggests that multivalent ligands containing flexible linkers that are longer than the spacing between the binding sites of a multivalent protein will be effective in binding, and that the use of flexible linkers with lengths somewhat greater than the optimal distance between binding sites is a justifiable strategy for the design of multivalent ligands.     

Insights into the interplay between electronic structure and protein dynamics: The case of ubiquitin

A quantum mechanical analysis of an experimental ensemble comprising 128 conformers of the protein ubiquitin has been carried out with the aid of LMO–SCF–COSMO calculations. The permanent dipole moment of the protein fluctuates in the range from 131 to 283 D while the energy-weighted average dipole has a magnitude of 197 D. The HOMO–LUMO energy gap of the conformational ensemble ranges from 7.389 to 8.397 eV and appears to being affected mainly by fluctuations in the HOMO energy. An inspection of the frontier orbitals of the 128 conformers indicates that their localization is not affected by the protein dynamics.     

Water‐in‐Oil Micro‐Emulsion Enhances the Secondary Structure of a Protein by Confinement

A scheme is presented in which an organic solvent environment in combination with surfactants is used to confine a natively unfolded protein inside an inverse microemulsion droplet. This type of confinement allows a study that provides unique insight into the dynamic structure of an unfolded, flexible protein which is still solvated and thus under near‐physiological conditions. In a model system, the protein osteopontin (OPN) is used. It is a highly phosphorylated glycoprotein that is expressed in a wide range of cells and tissues for which limited structural analysis exists due to the high degree of flexibility and large number of post‐translational modifications. OPN is implicated in tissue functions, such as inflammation and mineralisation. It also has a key function in tumour metastasis and progression. Circular dichroism measurements show that confinement enhances the secondary structural features of the protein. Small‐angle X‐ray scattering and dynamic light scattering show that OPN changes from being a flexible protein in aqueous solution to adopting a less flexible and more compact structure inside the microemulsion droplets. This novel approach for confining proteins while they are still hydrated may aid in studying the structure of a wide range of natively unfolded proteins.     

Monte Carlo simulation algorithm for B‐DNA

Understanding the structure–function relationship of biomolecules containing DNA has motivated experiments aimed at determining molecular structure using methods such as small‐angle X‐ray and neutron scattering (SAXS and SANS). SAXS and SANS are useful for determining macromolecular shape in solution, a process which benefits by using atomistic models that reproduce the scattering data. The variety of algorithms available for creating and modifying model DNA structures lack the ability to rapidly modify all‐atom models to generate structure ensembles. This article describes a Monte Carlo algorithm for simulating DNA, not with the goal of predicting an equilibrium structure, but rather to generate an ensemble of plausible structures which can be filtered using experimental results to identify a sub‐ensemble of conformations that reproduce the solution scattering of DNA macromolecules. The algorithm generates an ensemble of atomic structures through an iterative cycle in which B‐DNA is represented using a wormlike bead–rod model, new configurations are generated by sampling bend and twist moves, then atomic detail is recovered by back mapping from the final coarse‐grained configuration. Using this algorithm on commodity computing hardware, one can rapidly generate an ensemble of atomic level models, each model representing a physically realistic configuration that could be further studied using molecular dynamics. © 2016 Wiley Periodicals, Inc.