Structural basis of p38α regulation by hematopoietic tyrosine phosphatase

Linear motif specificity in signaling through p38α and ERK2 mitogen-activated protein kinases

Mitogen-activated protein kinase (MAPK) cascades are essential for eukaryotic cells to integrate and respond to diverse stimuli. Maintaining specificity in signaling through MAPK networks is key to coupling distinct inputs to appropriate cellular responses. Docking sites-short linear motifs found in MAPK substrates, regulators, and scaffolds-can promote signaling specificity through selective interactions, but how they do so remains unresolved. Here, we screened a proteomic library for sequences interacting with the MAPKs extracellular signal-regulated kinase 2 (ERK2) and p38α, identifying selective and promiscuous docking motifs. Sequences specific for p38α had high net charge and lysine content, and selective binding depended on a pair of acidic residues unique to the p38α docking interface. Finally, we validated a set of full-length proteins harboring docking sites selected in our screens to be authentic MAPK interactors and substrates. This study identifies features that help define MAPK signaling networks and explains how specific docking motifs promote signaling integrity.

Investigating Mammalian Formins with SMIFH2 Fifteen Years in: Novel Targets and Unexpected Biology

The mammalian formin family comprises fifteen multi-domain proteins that regulate actin dynamics and microtubules in vitro and in cells. Evolutionarily conserved formin homology (FH) 1 and 2 domains allow formins to locally modulate the cell cytoskeleton. Formins are involved in several developmental and homeostatic processes, as well as human diseases. However, functional redundancy has long hampered studies of individual formins with genetic loss-of-function approaches and prevents the rapid inhibition of formin activities in cells. The discovery of small molecule inhibitor of formin homology 2 domains (SMIFH2) in 2009 was a disruptive change that provided a powerful chemical tool to explore formins' functions across biological scales. Here, I critically discuss the characterization of SMIFH2 as a pan-formin inhibitor, as well as growing evidence of unexpected off-target effects. By collating the literature and information hidden in public repositories, outstanding controversies and fundamental open questions about the substrates and mechanism of action of SMIFH2 emerge. Whenever possible, I propose explanations for these discrepancies and roadmaps to address the paramount open questions. Furthermore, I suggest that SMIFH2 be reclassified as a multi-target inhibitor for its appealing activities on proteins involved in pathological formin-dependent processes. Notwithstanding all drawbacks and limitations, SMIFH2 will continue to prove useful in studying formins in health and disease in the years to come.

Proteome-wide screening for mitogen-activated protein kinase docking motifs and interactors

Essential functions of mitogen-activated protein kinases (MAPKs) depend on their capacity to selectively phosphorylate a limited repertoire of substrates. MAPKs harbor a conserved groove located outside of the catalytic cleft that binds to short linear sequence motifs found in substrates and regulators. However, the weak and transient nature of these "docking" interactions poses a challenge to defining MAPK interactomes and associated sequence motifs. Here, we describe a yeast-based genetic screening pipeline to evaluate large collections of MAPK docking sequences in parallel. Using this platform, we analyzed a combinatorial library based on the docking sequences from the MAPK kinases MKK6 and MKK7, defining features critical for binding to the stress-activated MAPKs JNK1 and p38α. Our screen of a library consisting of ~12,000 sequences from the human proteome revealed multiple MAPK-selective interactors, including many that did not conform to previously defined docking motifs. Analysis of p38α/JNK1 exchange mutants identified specific docking groove residues that mediate selective binding. Last, we verified that docking sequences identified in the screen functioned in substrate recruitment in vitro and in cultured cells. Together, these studies establish an approach to characterize MAPK docking sequences and provide a resource for future investigation of signaling downstream of p38 and JNK.

Comprehensive analysis of PTPN family expression and prognosis in acute myeloid leukemia

Background: Tyrosyl phosphorylation is carried out by a group of enzymes known as non-receptor protein tyrosine phosphatases (PTPNs). In the current investigation, it is hoped to shed light on the relationships between the expression patterns of PTPN family members and the prognosis of acute myeloid leukemia (AML). Methods: PTPN expression was examined using GEPIA and GEO databases. To investigate the connection between PTPN expression and survival in AML patients, we downloaded data from the Broad TCGA Firehose and Clinical Proteomic Tumor Analysis (CPTAC) of the Cancer Genome Atlas (TCGA). We used quantitative real-time PCR (qRT-PCR) to confirm that essential genes were performed in clinical samples and cell lines. We then used western blot to verify that the genes expressed in the above databases were positive in normal tissues, AML patient samples, and AML cell lines. Next, we investigated associations between genome-wide expression profiles and PTPN6 expression using the GEO datasets. We investigated the interactive exploration of multidimensional cancer genomics using the cBioPortal datasets. Using the DAVID database, a study of gene ontology enrichment was performed. The protein-protein interaction (PPI) network was created using the STRING portal, and the gene-gene interaction network was performed using GeneMANIA. Results: Data from GEO and GEPIA revealed that most PTPN family members were linked to AML. Patients with leukemia have elevated levels of several PTPN members. All of the AML patients' poor overall survival (OS, p < .05) was significantly linked with higher expression of PTPN1, PTPN6, and PTPN7. Additionally, clinical samples showed that the expression of PTPN 6, PTPN 7, PTPN 13, and PTPN 14 was higher than normal in AML patients (p = .0116, p = .0034, p = .0092, and p = .0057, respectively) and AML cell lines (p = .0004, p = .0035, p = .0357, and p = .0177, respectively). Western blotting results showed that the expression of PTPN6 in AML samples and AML cell lines was significantly higher than that in normal control samples. Conclusion: Differentially expressed PTPN family members were found in AML. The prognosis of patients and PTPN gene expression were shown to be correlated. PTPN6 is one of these members and may be used as an AML diagnostic and prognostic marker.

Mapping the deformability of natural and designed cellulosomes in solution

BACKGROUND

Natural cellulosome multi-enzyme complexes, their components, and engineered 'designer cellulosomes' (DCs) promise an efficient means of breaking down cellulosic substrates into valuable biofuel products. Their broad uptake in biotechnology relies on boosting proximity-based synergy among the resident enzymes, but the modular architecture challenges structure determination and rational design.

RESULTS

We used small angle X-ray scattering combined with molecular modeling to study the solution structure of cellulosomal components. These include three dockerin-bearing cellulases with distinct substrate specificities, original scaffoldins from the human gut bacterium Ruminococcus champanellensis (ScaA, ScaH and ScaK) and a trivalent cohesin-bearing designer scaffoldin (Scaf20L), followed by cellulosomal complexes comprising these components, and the nonavalent fully loaded Clostridium thermocellum CipA in complex with Cel8A from the same bacterium. The size analysis of Rg and Dmax values deduced from the scattering curves and corresponding molecular models highlight their variable aspects, depending on composition, size and spatial organization of the objects in solution.

CONCLUSIONS

Our data quantifies variability of form and compactness of cellulosomal components in solution and confirms that this native plasticity may well be related to speciation with respect to the substrate that is targeted. By showing that scaffoldins or components display enhanced compactness compared to the free objects, we provide new routes to rationally enhance their stability and performance in their environment of action.

The interaction of p38 with its upstream kinase MKK6

Mitogen-activated protein kinase (MAPK; p38, ERK, and JNK) cascades are evolutionarily conserved signaling pathways that regulate the cellular response to a variety of extracellular stimuli, such as growth factors and interleukins. The MAPK p38 is activated by its specific upstream MAPK kinases, MKK6 and MKK3. However, a comprehensive molecular understanding of how these cognate upstream kinases bind and activate p38 is still missing. Here, we combine NMR spectroscopy and isothermal titration calorimetry to define the binding interface between full-length MKK6 and p38. It was shown that p38 engages MKK6 not only via its hydrophobic docking groove, but also influences helix αF, a secondary structural element that plays a key role in organizing the kinase core. It was also shown that, unlike MAPK phosphatases, the p38 conserved docking (CD) site is much less affected by MKK6 binding. Finally, it was demonstrated that these interactions with p38 are conserved independent of the MKK6 activation state. Together, the results revealed differences between specificity markers of p38 regulation by upstream kinases, which do not effectively engage the CD site, and downstream phosphatases, which require the CD site for productive binding.

The methyl 13C-edited/13C-filtered transferred NOE for studying protein interactions with short linear motifs

Many proteins interact with their ligand proteins by recognition of short linear motifs that are often intrinsically disordered. These interactions are usually weak and are characterized by fast exchange. NMR spectroscopy is a powerful tool to study weak interactions. The methods that have been commonly used are analysis of chemicals shift perturbations (CSP) upon ligand binding and saturation transfer difference spectroscopy. These two methods identify residues at the binding interface between the protein and its ligand. In the present study, we used a combination of transferred-NOE, specific methyl-labeling and an optimized isotope-edited/isotope-filtered NOESY experiment to study specific interactions between the 42 kDa p38α mitogen-activated protein kinase and the kinase interaction motif (KIM) on the STEP phosphatase. These measurements distinguished between residues that both exhibit CSPs upon ligand binding and interact with the KIM peptide from residues that exhibit CSPs but do not interact with the peptide. In addition, these results provide information about pairwise interactions that is important for a more reliable docking of the KIM peptide into its interacting surface on p38α. This combination of techniques should be applicable for many protein-peptide complexes up to 80 kDa for which methyl resonance assignment can be achieved.

Capturing the Conformational Ensemble of the Mixed Folded Polyglutamine Protein Ataxin-3

Ataxin-3 is a deubiquitinase involved in protein quality control and other essential cellular functions. It preferentially interacts with polyubiquitin chains of four or more units attached to proteins delivered to the ubiquitin-proteasome system. Ataxin-3 is composed of an N-terminal Josephin domain and a flexible C terminus that contains two or three ubiquitin-interacting motifs (UIMs) and a polyglutamine tract, which, when expanded beyond a threshold, leads to protein aggregation and misfolding and causes spinocerebellar ataxia type 3. The high-resolution structure of the Josephin domain is available, but the structural and dynamical heterogeneity of ataxin-3 has so far hindered the structural description of the full-length protein. Here, we characterize non-expanded and expanded variants of ataxin-3 in terms of conformational ensembles adopted by the proteins in solution by jointly using experimental data from nuclear magnetic resonance and small-angle X-ray scattering with coarse-grained simulations. Our results pave the way to a molecular understanding of polyubiquitin recognition.

The expression patterns and the diagnostic/prognostic roles of PTPN family members in digestive tract cancers

Background

Non-receptor protein tyrosine phosphatases (PTPNs) are a set of enzymes involved in the tyrosyl phosphorylation. The present study intended to clarify the associations between the expression patterns of PTPN family members, and diagnosis as well as the prognosis of digestive tract cancers.

Methods

Oncomine and Ualcan were used to analyze PTPN expressions. Data from The Cancer Genome Atlas (TCGA) were downloaded through UCSC Xena for validation and to explore the relationship of the PTPN expression with diagnosis, clinicopathological parameters and survival of digestive tract cancers. Gene ontology enrichment analysis was conducted using the DAVID database. The gene–gene interaction network was performed by GeneMANIA and the protein–protein interaction (PPI) network was built using STRING portal coupled with Cytoscape. The expression of differentially expressed PTPNs in cancer cell lines were explored using CCLE. Moreover, by histological verification, the expression of four PTPNs in digestive tract cancers were further analyzed.

Results

Most PTPN family members were associated with digestive tract cancers according to Oncomine, Ualcan and TCGA data. Several PTPN members were differentially expressed in digestive tract cancers. For esophageal carcinoma (ESCA), PTPN1 and PTPN12 levels were correlated with incidence; PTPN20 was associated with poor prognosis. For stomach adenocarcinoma (STAD), PTPN2 and PTPN12 levels were correlated with incidence; PTPN3, PTPN5, PTPN7, PTPN11, PTPN13, PTPN14, PTPN18 and PTPN23 were correlated with pathological grade; PTPN20 expression was related with both TNM stage and N stage; PTPN22 was associated with T stage and pathological grade; decreased expression of PTPN5 and PTPN13 implied worse overall survival of STAD, while elevated PTPN6 expression indicated better prognosis. For colorectal cancer (CRC), PTPN2, PTPN21 and PTPN22 levels were correlated with incidence; expression of PTPN5, PTPN12, and PTPN14 was correlated with TNM stage and N stage; high PTPN5 or PTPN7 expression was associated with increased hazards of death. CCLE analyses showed that in esophagus cancer cell lines, PTPN1, PTPN4 and PTPN12 were highly expressed; in gastric cancer cell lines, PTPN2 and PTPN12 were highly expressed; in colorectal cancer cell lines, PTPN12 was highly expressed while PTPN22 was downregulated. Results of histological verification experiment showed differential expressions of PTPN22 in CRC, and PTPN12 in GC and CRC.

Conclusions

Members of PTPN family were differentially expressed in digestive tract cancers. Correlations were found between PTPN genes and clinicopathological parameters of patients. Expression of PTPN12 was upregulated in both STAD and CRC, and thus could be used as a diagnostic biomarker. Differential expression of PTPN12 in GC and CRC, and PTPN22 in CRC were presented in our histological verification experiment.

Combining molecular dynamics simulations with small-angle X-ray and neutron scattering data to study multi-domain proteins in solution

Many proteins contain multiple folded domains separated by flexible linkers, and the ability to describe the structure and conformational heterogeneity of such flexible systems pushes the limits of structural biology. Using the three-domain protein TIA-1 as an example, we here combine coarse-grained molecular dynamics simulations with previously measured small-angle scattering data to study the conformation of TIA-1 in solution. We show that while the coarse-grained potential (Martini) in itself leads to too compact conformations, increasing the strength of protein-water interactions results in ensembles that are in very good agreement with experiments. We show how these ensembles can be refined further using a Bayesian/Maximum Entropy approach, and examine the robustness to errors in the energy function. In particular we find that as long as the initial simulation is relatively good, reweighting against experiments is very robust. We also study the relative information in X-ray and neutron scattering experiments and find that refining against the SAXS experiments leads to improvement in the SANS data. Our results suggest a general strategy for studying the conformation of multi-domain proteins in solution that combines coarse-grained simulations with small-angle X-ray scattering data that are generally most easy to obtain. These results may in turn be used to design further small-angle neutron scattering experiments that exploit contrast variation through 1H/2H isotope substitutions.

How to learn from inconsistencies: Integrating molecular simulations with experimental data

Molecular simulations and biophysical experiments can be used to provide independent and complementary insights into the molecular origin of biological processes. A particularly useful strategy is to use molecular simulations as a modeling tool to interpret experimental measurements, and to use experimental data to refine our biophysical models. Thus, explicit integration and synergy between molecular simulations and experiments is fundamental for furthering our understanding of biological processes. This is especially true in the case where discrepancies between measured and simulated observables emerge. In this chapter, we provide an overview of some of the core ideas behind methods that were developed to improve the consistency between experimental information and numerical predictions. We distinguish between situations where experiments are used to refine our understanding and models of specific systems, and situations where experiments are used more generally to refine transferable models. We discuss different philosophies and attempt to unify them in a single framework. Until now, such integration between experiments and simulations have mostly been applied to equilibrium data, and we discuss more recent developments aimed to analyze time-dependent or time-resolved data.

A Dynamic Switch in Inactive p38γ Leads to an Excited State on the Pathway to an Active Kinase

The inactive state of mitogen-activated protein kinases (MAPKs) adopts an open conformation while the active state exists in a compact form stabilized by phosphorylation. In the active state, eukaryotic kinases undergo breathing motions related to substrate binding and product release that have not previously been detected in the inactive state. However, docking interactions of partner proteins with inactive MAPK kinases exhibit allostery in binding of activating kinases. Interactions at a site distant from the activation loop are coupled to the configuration of the activation loop, suggesting that the inactive state may also undergo concerted dynamics. X-ray crystallographic studies of nonphosphorylated, inactive p38γ reveal differences in domain orientations and active site structure in the two molecules in the asymmetric unit. One molecule resembles an inactive kinase with an open active site. The second molecule has a rotation of the N-lobe that leads to partial compaction of the active site, resulting in a conformation that is intermediate between the inactive open state and the fully closed state of the activated kinase. Although the compact state of apo p38γ displays several of the features of the activated enzyme, it remains catalytically inert. In solution, the kinase fluctuates on a millisecond time scale between the open ground state and a weakly populated excited state that is similar in structure to the compact state observed in the crystal. The nuclear magnetic resonance and crystal structure data imply that interconversion between the open and compact states involves a molecular switch associated with the DFG loop.

Structural insights into Acyl-coenzyme A binding domain containing 3 (ACBD3) protein hijacking by picornaviruses

Many picornaviruses hijack the Golgi resident Acyl-coenzyme A binding domain containing 3 (ACBD3) protein in order to recruit the phosphatidylinositol 4-kinase B (PI4KB) to viral replication organelles (ROs). PI4KB, once recruited and activated by ACBD3 protein, produces the lipid phosphatidylinositol 4-phosphate (PI4P), which is a key step in the biogenesis of viral ROs. To do so, picornaviruses use their small nonstructural protein 3A that binds the Golgi dynamics domain of the ACBD3 protein. Here, we present the analysis of the highly flexible ACBD3 proteins and the viral 3A protein in solution using small-angle X-ray scattering and computer simulations. Our analysis revealed that both the ACBD3 protein and the 3A:ACBD3 protein complex have an extended and flexible conformation in solution.

Preparation of Phosphorylated Proteins for NMR Spectroscopy

Phosphorylation is a ubiquitous posttranslational modification that is essential for the regulation of many cellular processes. The human genome consists of more than 200,000 phosphorylation sites, whose phosphorylation is tightly controlled by ≥500 kinases and ~200 phosphatases. Given the large number of phosphorylation sites and the key role phosphorylation plays in regulating cellular processes, it is essential to characterize the impact of phosphorylation on substrate structure, dynamics, and function. However, a major challenge is the large-scale production of phosphorylated proteins in vitro for these structural, functional, and dynamic studies. Here, we describe an efficient protocol used routinely in our laboratory for the production of phosphorylated proteins. We also describe the methods used for identifying, characterizing, and separating the resulting phosphorylated proteins for subsequent studies.

Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications

Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.

Efficient Ensemble Refinement by Reweighting

Ensemble refinement produces structural ensembles of flexible and dynamic biomolecules by integrating experimental data and molecular simulations. Here we present two efficient numerical methods to solve the computationally challenging maximum-entropy problem arising from a Bayesian formulation of ensemble refinement. Recasting the resulting constrained weight optimization problem into an unconstrained form enables the use of gradient-based algorithms. In two complementary formulations that differ in their dimensionality, we optimize either the log-weights directly or the generalized forces appearing in the explicit analytical form of the solution. We first demonstrate the robustness, accuracy, and efficiency of the two methods using synthetic data. We then use NMR J-couplings to reweight an all-atom molecular dynamics simulation ensemble of the disordered peptide Ala-5 simulated with the AMBER99SB*-ildn-q force field. After reweighting, we find a consistent increase in the population of the polyproline-II conformations and a decrease of α-helical-like conformations. Ensemble refinement makes it possible to infer detailed structural models for biomolecules exhibiting significant dynamics, such as intrinsically disordered proteins, by combining input from experiment and simulation in a balanced manner.

Inferring Structural Ensembles of Flexible and Dynamic Macromolecules Using Bayesian, Maximum Entropy, and Minimal-Ensemble Refinement Methods

The flexible and dynamic nature of biomolecules and biomolecular complexes is essential for many cellular functions in living organisms but poses a challenge for experimental methods to determine high-resolution structural models. To meet this challenge, experiments are combined with molecular simulations. The latter propose models for structural ensembles, and the experimental data can be used to steer these simulations and to select ensembles that most likely underlie the experimental data. Here, we explain in detail how the "Bayesian Inference Of ENsembles" (BioEn) method can be used to refine such ensembles using a wide range of experimental data. The "Ensemble Refinement of SAXS" (EROS) method is a special case of BioEn, inspired by the Gull-Daniell formulation of maximum entropy image processing and focused originally on X-ray solution scattering experiments (SAXS) and then extended to integrative structural modeling. We also briefly sketch the "minimum ensemble method," a maximum-parsimony refinement method that seeks to represent an ensemble with a minimal number of representative structures.

Dynamic activation and regulation of the mitogen-activated protein kinase p38

Mitogen-activated protein kinases, which include p38, are essential for cell differentiation and autophagy. The current model for p38 activation involves activation-loop phosphorylation with subsequent substrate binding leading to substrate phosphorylation. Despite extensive efforts, the molecular mechanism of activation remains unclear. Here, using NMR spectroscopy, we show how the modulation of protein dynamics across timescales activates p38. We find that activation-loop phosphorylation does not change the average conformation of p38; rather it quenches the loop ps-ns dynamics. We then show that substrate binding to nonphosphorylated and phosphorylated p38 results in uniform µs-ms backbone dynamics at catalytically essential regions and across the entire molecule, respectively. Together, these results show that phosphorylation and substrate binding flatten the energy landscape of the protein, making essential elements of allostery and activation dynamically accessible. The high degree of structural conservation among ser/thr kinases suggests that elements of this mechanism may be conserved across the kinase family.

1H, 15N and 13C sequence specific backbone assignment of the vanadate inhibited hematopoietic tyrosine phosphatase

The sequence-specific backbone assignment of hematopoietic protein tyrosine phosphatase (HePTP; PTPN7) in presence of vanadate has been determined, based on triple-resonance experiments using uniformly [¹³C,¹⁵N]-labeled protein. These assignments facilitate further studies of HePTP in the presence of inhibitors to target leukemia and provide further insights into the function of protein tyrosine phosphatases.

Homing in: Mechanisms of Substrate Targeting by Protein Kinases

Protein phosphorylation is the most common reversible post-translational modification in eukaryotes. Humans have over 500 protein kinases, of which more than a dozen are established targets for anticancer drugs. All kinases share a structurally similar catalytic domain, yet each one is uniquely positioned within signaling networks controlling essentially all aspects of cell behavior. Kinases are distinguished from one another based on their modes of regulation and their substrate repertoires. Coupling specific inputs to the proper signaling outputs requires that kinases phosphorylate a limited number of sites to the exclusion of hundreds of thousands of off-target phosphorylation sites. Here, we review recent progress in understanding mechanisms of kinase substrate specificity and how they function to shape cellular signaling networks.

The length but not the sequence of peptide linker modules exerts the primary influence on the conformations of protein domains in cellulosome multi-enzyme complexes

Cellulosomes are large multi-protein catalysts produced by various anaerobic microorganisms to efficiently degrade plant cell-wall polysaccharides down into simple sugars. X-ray and physicochemical structural characterisations show that cellulosomes are composed of numerous protein domains that are connected by unstructured polypeptide segments, yet the properties and possible roles of these 'linker' peptides are largely unknown. We have performed coarse-grained and all-atom molecular dynamics computer simulations of a number of cellulosomal linkers of different lengths and compositions. Our data demonstrates that the effective stiffness of the linker peptides, as quantified by the equilibrium fluctuations in the end-to-end distances, depends primarily on the length of the linker and less so on the specific amino acid sequence. The presence of excluded volume - provided by the domains that are connected - dampens the motion of the linker residues and reduces the effective stiffness of the linkers. Simultaneously, the presence of the linkers alters the conformations of the protein domains that are connected. We demonstrate that short, stiff linkers induce significant rearrangements in the folded domains of the mini-cellulosome composed of endoglucanase Cel8A in complex with scaffoldin ScafT (Cel8A-ScafT) of Clostridium thermocellum as well as in a two-cohesin system derived from the scaffoldin ScaB of Acetivibrio cellulolyticus. We give experimentally testable predictions on structural changes in protein domains that depend on the length of linkers.

Investigating energy-based pool structure selection in the structure ensemble modeling with experimental distance constraints: The example from a multidomain protein Pub1

The structural variations of multidomain proteins with flexible parts mediate many biological processes, and a structure ensemble can be determined by selecting a weighted combination of representative structures from a simulated structure pool, producing the best fit to experimental constraints such as interatomic distance. In this study, a hybrid structure-based and physics-based atomistic force field with an efficient sampling strategy is adopted to simulate a model di-domain protein against experimental paramagnetic relaxation enhancement (PRE) data that correspond to distance constraints. The molecular dynamics simulations produce a wide range of conformations depicted on a protein energy landscape. Subsequently, a conformational ensemble recovered with low-energy structures and the minimum-size restraint is identified in good agreement with experimental PRE rates, and the result is also supported by chemical shift perturbations and small-angle X-ray scattering data. It is illustrated that the regularizations of energy and ensemble-size prevent an arbitrary interpretation of protein conformations. Moreover, energy is found to serve as a critical control to refine the structure pool and prevent data overfitting, because the absence of energy regularization exposes ensemble construction to the noise from high-energy structures and causes a more ambiguous representation of protein conformations. Finally, we perform structure-ensemble optimizations with a topology-based structure pool, to enhance the understanding on the ensemble results from different sources of pool candidates.

Deciphering the Dynamic Interaction Profile of an Intrinsically Disordered Protein by NMR Exchange Spectroscopy

Intrinsically disordered proteins (IDPs) display a large number of interaction modes including folding-upon-binding, binding without major structural transitions, or binding through highly dynamic, so-called fuzzy, complexes. The vast majority of experimental information about IDP binding modes have been inferred from crystal structures of proteins in complex with short peptides of IDPs. However, crystal structures provide a mainly static view of the complexes and do not give information about the conformational dynamics experienced by the IDP in the bound state. Knowledge of the dynamics of IDP complexes is of fundamental importance to understand how IDPs engage in highly specific interactions without concomitantly high binding affinity. Here, we combine rotating-frame R_1ρ, Carr-Purcell-Meiboom Gill relaxation dispersion as well as chemical exchange saturation transfer to decipher the dynamic interaction profile of an IDP in complex with its partner. We apply the approach to the dynamic signaling complex formed between the mitogen-activated protein kinase (MAPK) p38α and the intrinsically disordered regulatory domain of the MAPK kinase MKK4. Our study demonstrates that MKK4 employs a subtle combination of interaction modes in order to bind to p38α, leading to a complex displaying significantly different dynamics across the bound regions.

NMR Spectroscopy to Study MAP Kinase Binding to MAP Kinase Phosphatases

NMR spectroscopy and other solution methods are increasingly being used to obtain novel insights into the mechanisms by which MAPK regulatory proteins bind and direct the activity of MAPKs. Here, we describe how interactions between the MAPK p38α and its regulatory proteins are studied using NMR spectroscopy, isothermal titration calorimetry, and small angle X-ray scattering (SAXS).

Structures of Dynamic Protein Complexes: Hybrid Techniques to Study MAP Kinase Complexes and the ESCRT System

The integration of complementary molecular methods (including X-ray crystallography, NMR spectroscopy, small angle X-ray/neutron scattering, and computational techniques) is frequently required to obtain a comprehensive understanding of dynamic macromolecular complexes. In particular, these techniques are critical for studying intrinsically disordered protein regions (IDRs) or intrinsically disordered proteins (IDPs) that are part of large protein:protein complexes. Here, we explain how to prepare IDP samples suitable for study using NMR spectroscopy, and describe a novel SAXS modeling method (ensemble refinement of SAXS; EROS) that integrates the results from complementary methods, including crystal structures and NMR chemical shift perturbations, among others, to accurately model SAXS data and describe ensemble structures of dynamic macromolecular complexes.

Design, Synthesis, and Biological Evaluation of Novel Type I1/2 p38α MAP Kinase Inhibitors with Excellent Selectivity, High Potency, and Prolonged Target Residence Time by Interfering with the R-Spine

We recently reported 1a (skepinone-L) as a type I p38α MAP kinase inhibitor with high potency and excellent selectivity in vitro and in vivo. However, as a type I inhibitor, it is entirely ATP-competitive and shows just a moderate residence time. Thus, the scope was to develop a new class of advanced compounds maintaining the structural binding features of skepinone-L scaffold like inducing a glycine flip at the hinge region and occupying both hydrophobic regions I and II. Extending this scaffold with suitable residues resulted in an interference with the kinase's R-Spine. By synthesizing 69 compounds, we could significantly prolong the target residence time with one example to 3663 s, along with an excellent selectivity score of 0.006 and an outstanding potency of 1.0 nM. This new binding mode was validated by cocrystallization, showing all binding interactions typifying type I¹/₂ binding. Moreover, microsomal studies showed convenient metabolic stability of the most potent, herein reported representatives.

Redox-dependent dimerization of p38α mitogen-activated protein kinase with mitogen-activated protein kinase kinase 3

The kinase p38α MAPK (p38α) plays a pivotal role in many biological processes. p38α is activated by canonical upstream kinases that phosphorylate the activation region. The purpose of our study was to determine whether such activation may depend on redox-sensing cysteines within p38α. p38α was activated and formed a disulfide-bound heterodimer with MAP2K3 (MKK3) in rat cardiomyocytes and isolated hearts exposed to H₂O₂. This disulfide heterodimer was sensitive to reduction by mercaptoethanol and was enhanced by the thioredoxin-reductase inhibitor auranofin. We predicted that Cys-119 or Cys-162 of p38α, close to the known MKK3 docking domain, were relevant for these redox characteristics. The C119S mutation decreased whereas the C162S mutation increased the dimer formation, suggesting that these two Cys residues act as vicinal thiols, consistent with C119S/C162S being incapable of sensing H₂O₂. Similarly, disulfide heterodimer formation was abolished in H9C2 cells expressing both MKK3 and p38α C119S/C162S and subjected to simulated ischemia and reperfusion. However, the p38α C119S/C162S mutants did not exhibit appreciable alteration in activating dual phosphorylation. In contrast, the anti-inflammatory agent 10-nitro-oleic acid (NO₂-OA), a component of the Mediterranean diet, reduced p38α activation and covalently modified Cys-119/Cys-162, probably obstructing MKK3 access. Moreover, NO₂-OA reduced the dephosphorylation of p38α by hematopoietic tyrosine phosphatase (HePTP). Furthermore, steric obstruction of Cys-119/Cys-162 by NO₂-OA pretreatment in Langendorff-perfused murine hearts prevented the p38-MKK3 disulfide dimer formation and attenuated H₂O₂-induced contractile dysfunction. Our findings suggest that cysteine residues within p38α act as redox sensors that can dynamically regulate the association between p38 and MKK3.

Stiffness of the C-terminal disordered linker affects the geometry of the active site in endoglucanase Cel8A

Cellulosomes are complex multi-enzyme machineries which efficiently degrade plant cell-wall polysaccharides. The multiple domains of the cellulosome proteins are often tethered together by intrinsically disordered regions. The properties and functions of these disordered linkers are not well understood. In this work, we study endoglucanase Cel8A, which is a relevant enzymatic component of the cellulosomes of Clostridium thermocellum. We use both all-atom and coarse-grained simulations to investigate how the conformations of the catalytic domain of Cel8A are affected by the disordered linker at its C terminus. We find that when the endoglucanase is bound to its substrate, the effective stiffness of the linker can influence the distances between groups of amino-acid residues throughout the entire enzymatic domain. In particular, variations in the linker stiffness can lead to small changes in the geometry of the active-site cleft. We suggest that such geometrical changes may have an effect on the catalytic activity of the enzyme.

The KIM-family protein-tyrosine phosphatases use distinct reversible oxidation intermediates: Intramolecular or intermolecular disulfide bond formation

The kinase interaction motif (KIM) family of protein-tyrosine phosphatases (PTPs) includes hematopoietic protein-tyrosine phosphatase (HePTP), striatal-enriched protein-tyrosine phosphatase (STEP), and protein-tyrosine phosphatase receptor type R (PTPRR). KIM-PTPs bind and dephosphorylate mitogen-activated protein kinases (MAPKs) and thereby critically modulate cell proliferation and differentiation. PTP activity can readily be diminished by reactive oxygen species (ROS), e.g. H₂O₂, which oxidize the catalytically indispensable active-site cysteine. This initial oxidation generates an unstable sulfenic acid intermediate that is quickly converted into either a sulfinic/sulfonic acid (catalytically dead and irreversible inactivation) or a stable sulfenamide or disulfide bond intermediate (reversible inactivation). Critically, our understanding of ROS-mediated PTP oxidation is not yet sufficient to predict the molecular responses of PTPs to oxidative stress. However, identifying distinct responses will enable novel routes for PTP-selective drug design, important for managing diseases such as cancer and Alzheimer's disease. Therefore, we performed a detailed biochemical and molecular study of all KIM-PTP family members to determine their H₂O₂ oxidation profiles and identify their reversible inactivation mechanism(s). We show that despite having nearly identical 3D structures and sequences, each KIM-PTP family member has a unique oxidation profile. Furthermore, we also show that whereas STEP and PTPRR stabilize their reversibly oxidized state by forming an intramolecular disulfide bond, HePTP uses an unexpected mechanism, namely, formation of a reversible intermolecular disulfide bond. In summary, despite being closely related, KIM-PTPs significantly differ in oxidation profiles. These findings highlight that oxidation protection is critical when analyzing PTPs, for example, in drug screening.

The D Domain of LRRC4 anchors ERK1/2 in the cytoplasm and competitively inhibits MEK/ERK activation in glioma cells

Background

As a well-characterized key player in various signal transduction networks, extracellular-signal-regulated kinase (ERK1/2) has been widely implicated in the development of many malignancies. We previously found that Leucine-rich repeat containing 4 (LRRC4) was a tumor suppressor and a negative regulator of the ERK/MAPK pathway in glioma tumorigenesis. However, the precise molecular role of LRRC4 in ERK signal transmission is unclear.

Methods

The interaction between LRRC4 and ERK1/2 was assessed by co-immunoprecipitation and GST pull-down assays in vivo and in vitro. We also investigated the interaction of LRRC4 and ERK1/2 and the role of the D domain in ERK activation in glioma cells.

Results

Here, we showed that LRRC4 and ERK1/2 interact via the D domain and CD domain, respectively. Following EGF stimuli, the D domain of LRRC4 anchors ERK1/2 in the cytoplasm and abrogates ERK1/2 activation and nuclear translocation. In glioblastoma cells, ectopic LRRC4 expression competitively inhibited the interaction of endogenous mitogen-activated protein kinase (MEK) and ERK1/2. Mutation of the D domain decreased the LRRC4-mediated inhibition of MAPK signaling and its anti-proliferation and anti-invasion roles.

Conclusions

Our results demonstrated that the D domain of LRRC4 anchors ERK1/2 in the cytoplasm and competitively inhibits MEK/ERK activation in glioma cells. These findings identify a new mechanism underlying glioblastoma progression and suggest a novel therapeutic strategy by restoring the activity of LRRC4 to decrease MAPK cascade activation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13045-016-0355-1) contains supplementary material, which is available to authorized users.

Structural Basis for the Subversion of MAP Kinase Signaling by an Intrinsically Disordered Parasite Secreted Agonist

The causative agent of toxoplasmosis, the intracellular parasite Toxoplasma gondii, delivers a protein, GRA24, into the cells it infects that interacts with the mitogen-activated protein (MAP) kinase p38α (MAPK14), leading to activation and nuclear translocation of the host kinase and a subsequent inflammatory response that controls the progress of the parasite. The purification of a recombinant complex of GRA24 and human p38α has allowed the molecular basis of this activation to be determined. GRA24 is shown to be intrinsically disordered, binding two kinases that act independently, and is the only factor required to bypass the canonical mitogen-activated protein kinase activation pathway. An adapted kinase interaction motif (KIM) forms a highly stable complex that competes with cytoplasmic regulatory partners. In addition, the recombinant complex forms a powerful in vitro tool to evaluate the specificity and effectiveness of p38α inhibitors that have advanced to clinical trials, as it provides a hitherto unavailable stable and highly active form of p38α.

•

Toxoplasmosis controls its host immune response via a protein effector, GRA24

•

A recombinant complex of GRA24 and MAPK p38α demonstrates how the protein works

•

An adapted KIM domain ensures activation and a sustained inflammatory response

•

The recombinant complex is useful in the evaluation of p38 inhibitors

Bayesian ensemble refinement by replica simulations and reweighting

We describe different Bayesian ensemble refinement methods, examine their interrelation, and discuss their practical application. With ensemble refinement, the properties of dynamic and partially disordered (bio)molecular structures can be characterized by integrating a wide range of experimental data, including measurements of ensemble-averaged observables. We start from a Bayesian formulation in which the posterior is a functional that ranks different configuration space distributions. By maximizing this posterior, we derive an optimal Bayesian ensemble distribution. For discrete configurations, this optimal distribution is identical to that obtained by the maximum entropy "ensemble refinement of SAXS" (EROS) formulation. Bayesian replica ensemble refinement enhances the sampling of relevant configurations by imposing restraints on averages of observables in coupled replica molecular dynamics simulations. We show that the strength of the restraints should scale linearly with the number of replicas to ensure convergence to the optimal Bayesian result in the limit of infinitely many replicas. In the "Bayesian inference of ensembles" method, we combine the replica and EROS approaches to accelerate the convergence. An adaptive algorithm can be used to sample directly from the optimal ensemble, without replicas. We discuss the incorporation of single-molecule measurements and dynamic observables such as relaxation parameters. The theoretical analysis of different Bayesian ensemble refinement approaches provides a basis for practical applications and a starting point for further investigations.

Investigating the Role of Large-Scale Domain Dynamics in Protein-Protein Interactions

Intrinsically disordered linkers provide multi-domain proteins with degrees of conformational freedom that are often essential for function. These highly dynamic assemblies represent a significant fraction of all proteomes, and deciphering the physical basis of their interactions represents a considerable challenge. Here we describe the difficulties associated with mapping the large-scale domain dynamics and describe two recent examples where solution state methods, in particular NMR spectroscopy, are used to investigate conformational exchange on very different timescales.

NMR Characterization of Information Flow and Allosteric Communities in the MAP Kinase p38γ

The intramolecular network structure of a protein provides valuable insights into allosteric sites and communication pathways. However, a straightforward method to comprehensively map and characterize these pathways is not currently available. Here we present an approach to characterize intramolecular network structure using NMR chemical shift perturbations. We apply the method to the mitogen activated protein kinase (MAPK) p38γ. p38γ contains allosteric sites that are conserved among eukaryotic kinases as well as unique to the MAPK family. How these regulatory sites communicate with catalytic residues is not well understood. Using our method, we observe and characterize for the first time information flux between regulatory sites through a conserved kinase infrastructure. This network is accessed, reinforced, and broken in various states of p38γ, reflecting the functional state of the protein. We demonstrate that the approach detects critical junctions in the network corresponding to biologically significant allosteric sites and pathways.

Molecular Basis of the Interaction of the Human Protein Tyrosine Phosphatase Non-receptor Type 4 (PTPN4) with the Mitogen-activated Protein Kinase p38γ

The human protein tyrosine phosphatase non-receptor type 4 (PTPN4) prevents cell death induction in neuroblastoma and glioblastoma cell lines in a PDZ·PDZ binding motifs-dependent manner, but the cellular partners of PTPN4 involved in cell protection are unknown. Here, we described the mitogen-activated protein kinase p38γ as a cellular partner of PTPN4. The main contribution to the p38γ·PTPN4 complex formation is the tight interaction between the C terminus of p38γ and the PDZ domain of PTPN4. We solved the crystal structure of the PDZ domain of PTPN4 bound to the p38γ C terminus. We identified the molecular basis of recognition of the C-terminal sequence of p38γ that displays the highest affinity among all endogenous partners of PTPN4. We showed that the p38γ C terminus is also an efficient inducer of cell death after its intracellular delivery. In addition to recruiting the kinase, the binding of the C-terminal sequence of p38γ to PTPN4 abolishes the catalytic autoinhibition of PTPN4 and thus activates the phosphatase, which can efficiently dephosphorylate the activation loop of p38γ. We presume that the p38γ·PTPN4 interaction promotes cellular signaling, preventing cell death induction.

Applications of isothermal titration calorimetry - the research and technical developments from 2011 to 2015

Isothermal titration calorimetry is a widely used biophysical technique for studying the formation or dissociation of molecular complexes. Over the last 5 years, much work has been published on the interpretation of isothermal titration calorimetry (ITC) data for single binding and multiple binding sites. As over 80% of ITC papers are on macromolecules of biological origin, this interpretation is challenging. Some researchers have attempted to link the thermodynamics constants to events at the molecular level. This review highlights work carried out using binding sites characterized using x-ray crystallography techniques that allow speculation about individual bond formation and the displacement of individual water molecules during ligand binding and link these events to the thermodynamic constants for binding. The review also considers research conducted with synthetic binding partners where specific binding events like anion-π and π-π interactions were studied. The revival of assays that enable both thermodynamic and kinetic information to be collected from ITC data is highlighted. Lastly, published criticism of ITC research from a physical chemistry perspective is appraised and practical advice provided for researchers unfamiliar with thermodynamics and its interpretation. Copyright © 2016 John Wiley & Sons, Ltd.

Low potency toxins reveal dense interaction networks in metabolism

BACKGROUND

The chemicals of metabolism are constructed of a small set of atoms and bonds. This may be because chemical structures outside the chemical space in which life operates are incompatible with biochemistry, or because mechanisms to make or utilize such excluded structures has not evolved. In this paper I address the extent to which biochemistry is restricted to a small fraction of the chemical space of possible chemicals, a restricted subset that I call Biochemical Space. I explore evidence that this restriction is at least in part due to selection again specific structures, and suggest a mechanism by which this occurs.

RESULTS

Chemicals that contain structures that our outside Biochemical Space (UnBiological groups) are more likely to be toxic to a wide range of organisms, even though they have no specifically toxic groups and no obvious mechanism of toxicity. This correlation of UnBiological with toxicity is stronger for low potency (millimolar) toxins. I relate this to the observation that most chemicals interact with many biological structures at low millimolar toxicity. I hypothesise that life has to select its components not only to have a specific set of functions but also to avoid interactions with all the other components of life that might degrade their function.

CONCLUSIONS

The chemistry of life has to form a dense, self-consistent network of chemical structures, and cannot easily be arbitrarily extended. The toxicity of arbitrary chemicals is a reflection of the disruption to that network occasioned by trying to insert a chemical into it without also selecting all the other components to tolerate that chemical. This suggests new ways to test for the toxicity of chemicals, and that engineering organisms to make high concentrations of materials such as chemical precursors or fuels may require more substantial engineering than just of the synthetic pathways involved.

Two hydrophobic residues can determine the specificity of mitogen-activated protein kinase docking interactions

MAPKs bind to many of their upstream regulators and downstream substrates via a short docking motif (the D-site) on their binding partner. MAPKs that are in different families (e.g. ERK, JNK, and p38) can bind selectively to D-sites in their authentic substrates and regulators while discriminating against D-sites in other pathways. Here we demonstrate that the short hydrophobic region at the distal end of the D-site plays a critical role in determining the high selectivity of JNK MAPKs for docking sites in their cognate MAPK kinases. Changing just 1 or 2 key hydrophobic residues in this submotif is sufficient to turn a weak JNK-binding D-site into a strong one, or vice versa. These specificity-determining differences are also found in the D-sites of the ETS family transcription factors Elk-1 and Net. Moreover, swapping two hydrophobic residues between these D-sites switches the relative efficiency of Elk-1 and Net as substrates for ERK versus JNK, as predicted. These results provide new insights into docking specificity and suggest that this specificity can evolve rapidly by changes to just 1 or 2 amino acids.

Structure-Based Assignment of Ile, Leu, and Val Methyl Groups in the Active and Inactive Forms of the Mitogen-Activated Protein Kinase Extracellular Signal-Regulated Kinase 2

Resonance assignments are the first step in most NMR studies of protein structure, function, and dynamics. Standard protein assignment methods employ through-bond backbone experiments on uniformly (13)C/(15)N-labeled proteins. For larger proteins, this through-bond assignment procedure often breaks down due to rapid relaxation and spectral overlap. The challenges involved in studies of larger proteins led to efficient methods for (13)C labeling of side chain methyl groups, which have favorable relaxation properties and high signal-to-noise. These methyls are often still assigned by linking them to the previously assigned backbone, thus limiting the applications for larger proteins. Here, a structure-based procedure is described for assignment of (13)C(1)H3-labeled methyls by comparing distance information obtained from three-dimensional methyl-methyl nuclear Overhauser effect (NOE) spectroscopy with the X-ray structure. The Ile, Leu, or Val (ILV) methyl type is determined by through-bond experiments, and the methyl-methyl NOE data are analyzed in combination with the known structure. A hierarchical approach was employed that maps the largest observed "NOE-methyl cluster" onto the structure. The combination of identification of ILV methyl type with mapping of the NOE-methyl clusters greatly simplifies the assignment process. This method was applied to the inactive and active forms of the 42-kDa ILV (13)C(1)H3-methyl labeled extracellular signal-regulated kinase 2 (ERK2), leading to assignment of 60% of the methyls, including 90% of Ile residues. A series of ILV to Ala mutants were analyzed, which helped confirm the assignments. These assignments were used to probe the local and long-range effects of ligand binding to inactive and active ERK2.

Structure and dynamics of the MKK7-JNK signaling complex

Signaling specificity in the mitogen-activated protein kinase (MAPK) pathways is controlled by disordered domains of the MAPK kinases (MKKs) that specifically bind to their cognate MAPKs via linear docking motifs. MKK7 activates the c-Jun N-terminal kinase (JNK) pathway and is the only MKK containing three motifs within its regulatory domain. Here, we characterize the conformational behavior and interaction mechanism of the MKK7 regulatory domain. Using NMR spectroscopy, we develop an atomic resolution ensemble description of MKK7, revealing highly diverse intrinsic conformational propensities of the three docking sites, suggesting that prerecognition sampling of the bound-state conformation is not prerequisite for binding. Although the different sites exhibit similar affinities for JNK1, interaction kinetics differ considerably. Importantly, we determine the crystal structure of JNK1 in complex with the second docking site of MKK7, revealing two different binding modes of the docking motif correlating with observations from NMR exchange spectroscopy. Our results provide unique insight into how signaling specificity is regulated by linear motifs and, in general, into the role of conformational disorder in MAPK signaling.

Ensemble Calculation for Intrinsically Disordered Proteins Using NMR Parameters

Intrinsically disordered proteins (IDPs) perform their function despite their lack of well-defined tertiary structure. Residual structure has been observed in IDPs, commonly described as transient/dynamic or expressed in terms of fractional populations. In order to understand how the protein primary sequence dictates the dynamic and structural properties of IDPs and in general to understand how IDPs function, atomic-level descriptions are needed. Nuclear magnetic resonance spectroscopy provides information about local and long-range structure in IDPs at amino acid specific resolution and can be used in combination with ensemble descriptions to represent the dynamic nature of IDPs. In this chapter we describe sample-and-select approaches for ensemble modelling of local structural propensities in IDPs with specific emphasis on validation of these ensembles.

Application of SAXS for the Structural Characterization of IDPs

Small-angle X-ray scattering (SAXS) is a powerful structural method allowing one to study the structure, folding state and flexibility of native particles and complexes in solution and to rapidly analyze structural changes in response to variations in external conditions. New high brilliance sources and novel data analysis methods significantly enhanced resolution and reliability of structural models provided by the technique. Automation of the SAXS experiment, data processing and interpretation make solution SAXS a streamline tool for large scale structural studies in molecular biology. The method provides low resolution macromolecular shapes ab initio and is readily combined with other structural and biochemical techniques in integrative studies. Very importantly, SAXS is sensitive to macromolecular flexibility being one of the few structural techniques applicable to flexible systems and intrinsically disordered proteins (IDPs). A major recent development is the use of SAXS to study particle dynamics in solution by ensemble approaches, which allow one to quantitatively characterize flexible systems. Of special interest is the joint use of SAXS with solution NMR, given that both methods yield highly complementary structural information, in particular, for IDPs. In this chapter, we present the basics of SAXS and also consider protocols of the experiment and data analysis for different scenarios depending on the type of the studied object. These include ab initio shape reconstruction, validation of available high resolution structures and rigid body modelling for folded macromolecules and also characterisation of flexible proteins with the ensemble methods. The methods are illustrated by examples of recent applications and further perspectives of the integrative use of SAXS with NMR in the studies of IDPs are discussed.

Methods for SAXS-based structure determination of biomolecular complexes

Measurements from small-angle X-ray scattering (SAXS) are highly informative to determine the structures of bimolecular complexes in solution. Here, current and recent SAXS-driven developments are described, with an emphasis on computational modeling. In particular, accurate methods to computing one theoretical scattering profile from a given structure model are discussed, with a key focus on structure factor coarse-graining and hydration contribution. Methods for reconstructing topological structures from an experimental SAXS profile are currently under active development. We report on several modeling tools designed for conformation generation that make use of either atomic-level or coarse-grained representations. Furthermore, since large, flexible biomolecules can adopt multiple well-defined conformations, a traditional single-conformation SAXS analysis is inappropriate, so we also discuss recent methods that utilize the concept of ensemble optimization, weighing in on the SAXS contributions of a heterogeneous mixture of conformations. These tools will ultimately posit the usefulness of SAXS data beyond a simple space-filling approach by providing a reliable structure characterization of biomolecular complexes under physiological conditions.

Large, dynamic, multi-protein complexes: a challenge for structural biology

Structural biology elucidates atomic structures of macromolecules such as proteins, DNA, RNA, and their complexes to understand the basic mechanisms of their functions. Among proteins that pose the most difficult problems to current efforts are those which have several large domains connected by long, flexible polypeptide segments. Although abundant and critically important in biological cells, such proteins have proven intractable by conventional techniques. This gap has recently led to the advancement of hybrid methods that use state-of-the-art computational tools to combine complementary data from various high- and low-resolution experiments. In this review, we briefly discuss the individual experimental techniques to illustrate their strengths and limitations, and then focus on the use of hybrid methods in structural biology. We describe how representative structures of dynamic multi-protein complexes are obtained utilizing the EROS hybrid method that we have co-developed.

Recovering a representative conformational ensemble from underdetermined macromolecular structural data

Structural analysis of proteins and nucleic acids is complicated by their inherent flexibility, conferred, for example, by linkers between their contiguous domains. Therefore, the macromolecule needs to be represented by an ensemble of conformations instead of a single conformation. Determining this ensemble is challenging because the experimental data are a convoluted average of contributions from multiple conformations. As the number of the ensemble degrees of freedom generally greatly exceeds the number of independent observables, directly deconvolving experimental data into a representative ensemble is an ill-posed problem. Recent developments in sparse approximations and compressive sensing have demonstrated that useful information can be recovered from underdetermined (ill-posed) systems of linear equations by using sparsity regularization. Inspired by these advances, we designed the Sparse Ensemble Selection (SES) method for recovering multiple conformations from a limited number of observations. SES is more general and accurate than previously published minimum-ensemble methods, and we use it to obtain representative conformational ensembles of Lys48-linked diubiquitin, characterized by the residual dipolar coupling data measured at several pH conditions. These representative ensembles are validated against NMR chemical shift perturbation data and compared to maximum-entropy results. The SES method reproduced and quantified the previously observed pH dependence of the major conformation of Lys48-linked diubiquitin, and revealed lesser-populated conformations that are preorganized for binding known diubiquitin receptors, thus providing insights into possible mechanisms of receptor recognition by polyubiquitin. SES is applicable to any experimental observables that can be expressed as a weighted linear combination of data for individual states.