Molecular dynamics reveal the essential role of linker motions in the function of cullin-RING E3 ligases

Impaired nuclear transport induced by juvenile ALS causing P525L mutation in NLS domain of FUS: A molecular mechanistic study

Amyotrophic lateral sclerosis (ALS) and fronto-temporal lobar degeneration (FTLD) are progressive neurological disorders affecting motor neurons. Cellular aggregates of fused in sarcoma (FUS) protein are found in cytoplasm of ALS and FTLD patients. Nuclear localisation signal (NLS) domain of FUS binds to Karyopherin β2 (Kapβ2), which drives nuclear transport of FUS from cytoplasm. Several pathogenic mutations are reported in FUS NLS, which are associated with its impaired nuclear transport and cytoplasmic mis-localisation. P525L mutation in NLS is most commonly found in cases of juvenile ALS (jALS), which affects individuals below 25 years of age. jALS progresses aggressively causing death within a year of its onset. This study elucidates the molecular mechanism behind jALS-causing P525L mutation hindering nuclear transport of FUS. We perform multiple molecular dynamics simulations in aqueous and hydrophobic solvent to understand the effect of the mutation at molecular level. Dynamics of Kapβ2-FUS complex is better captured in hydrophobic solvent compared to aqueous solvent. P525 and Y526 (PY-motif) of NLS exhibit fine-tuned stereochemical arrangement, which is essential for optimum Kapβ2 binding. P525L causes loss of several native contacts at interface leading to weaker binding, which promotes self-aggregation of FUS in cytoplasm. Native complex samples closed conformation, while mutant complex exhibits open conformation exposing hydrophilic residues of Kapβ2 to hydrophobic solvent. Mutant complex also fails to exhibit spring-like motion essential for its transport through nuclear pore complex. This study provides a mechanistic insight of binding affinity between NLS and Kapβ2 that inhibits self-aggregation of FUS preventing the disease condition.

The diversity of molecular interactions involving intrinsically disordered proteins: A molecular modeling perspective

Intrinsically Disordered Proteins and Regions (IDPs/IDRs) are key components of a multitude of biological processes. Conformational malleability enables IDPs/IDRs to perform very specialized functions that cannot be accomplished by globular proteins. The functional role for most of these proteins is related to the recognition of other biomolecules to regulate biological processes or as a part of signaling pathways. Depending on the extent of disorder, the number of interacting sites and the type of partner, very different architectures for the resulting assemblies are possible. More recently, molecular condensates with liquid-like properties composed of multiple copies of IDPs and nucleic acids have been proven to regulate key processes in eukaryotic cells. The structural and kinetic details of disordered biomolecular complexes are difficult to unveil experimentally due to their inherent conformational heterogeneity. Computational approaches, alone or in combination with experimental data, have emerged as unavoidable tools to understand the functional mechanisms of this elusive type of assemblies. The level of description used, all-atom or coarse-grained, strongly depends on the size of the molecular systems and on the timescale of the investigated mechanism. In this mini-review, we describe the most relevant architectures found for molecular interactions involving IDPs/IDRs and the computational strategies applied for their investigation.

Structural Model for Recruitment of RIT1 to the LZTR1 E3 Ligase: Evidences from an Integrated Computational Approach

Leucine-zipper transcription regulator 1 (LZTR1) is a highly mutated tumor suppressor gene, involved in the pathogenesis of several cancer types and developmental disorders. In proteasomal degradation, it acts as an adaptor protein responsible for the recognition and recruitment of substrates to be ubiquitinated in Cullin3-RING ligase E3 (CRL3) machinery. LZTR1 belongs to the BTB-Kelch family, a multi-domain protein where the Kelch propeller plays as the substrate recognition region and for which no experimental structure has been solved. Recently, large effort mutational analyses pointed to the role of disease-associated LZTR1 mutations in the RAS/MAPK signaling pathway and RIT1, a small Ras-related GTPase protein, has been identified by mass spectroscopy to interact with LZTR1. Hence, a better understanding of native structure, molecular mechanism, and substrate specificity would help clarifying the role of LZTR1 in pathological diseases, thus promoting advancement in the development of novel therapeutic strategies. Here, we address the interaction model between adaptor LZTR1 and substrate RIT1 by applying an integrated computational approach, including molecular modeling and docking techniques. We observe that the interaction model LZTR1-RIT1 is stabilized by an electrostatic bond network established between the two protein surfaces, which is reminiscent of homologous ubiquitin ligases complexes. Then, running MD simulations, we characterize differential conformational dynamics of the multi-domain LZTR1, offering interesting implications on the mechanistic role of specific point mutations. We identify G248R and R283Q as damaging mutations involved in the recognition process of the substrate RIT1 and R412C as a possible allosteric mutation from the Kelch to the C-term BTB-domain. Our findings provide important structural insights on targeting CRL3s for drug discovery.

Unusual RNA binding of FUS RRM studied by molecular dynamics simulation and enhanced sampling method

Amyotrophic lateral sclerosis (ALS) and frontotemporal lobe degeneration (FTLD) are two inter-related intractable diseases of motor neuron degeneration. Fused in sarcoma (FUS) is found in cytoplasmic accumulation of ALS and FTLD patients, which readily link the protein with the diseases. The RNA recognition motif (RRM) of FUS has the canonical α-β folds along with an unusual lysine-rich loop (KK-loop) between α1 and β2. This KK-loop is highly conserved among FET family proteins. Another contrasting feature of FUS RRM is the absence of critical binding residues, which are otherwise highly conserved in canonical RRMs. These residues in FUS RRM are Thr286, Glu336, Thr338, and Ser367, which are substitutions of lysine, phenylalanine, phenylalanine, and lysine, respectively, in other RRMs. Considering the importance of FUS in RNA regulation and metabolism, and its implication in ALS and FTLD, it is important to elucidate the underlying molecular mechanism of RNA recognition. In this study, we have performed molecular dynamics simulation with enhanced sampling to understand the conformational dynamics of noncanonical FUS RRM and its binding with RNA. We studied two sets of mutations: one with alanine mutation of KK-loop and another with KK-loop mutations along with critical binding residues mutated back to their canonical form. We find that concerted movement of KK-loop and loop between β2 and β3 facilitates the folding of the partner RNA, indicating an induced-fit mechanism of RNA binding. Flexibility of the RRM is highly restricted upon mutating the lysine residues of the KK-loop, resulting in weaker binding with the RNA. Our results also suggest that absence of the canonical residues in FUS RRM along with the KK-loop is equally important in regulating its binding dynamics. This study provides a significant structural insight into the binding of FUS RRM with its cognate RNA, which may further help in designing potential drugs targeting noncanonical RNA recognition.

High-speed atomic force microscopy directly visualizes conformational dynamics of the HIV Vif protein in complex with three host proteins

Vif (viral infectivity factor) is a protein that is essential for the replication of the HIV-1 virus. The key function of Vif is to disrupt the antiviral activity of host APOBEC3 (apolipoprotein B mRNA-editing enzyme catalytic subunit 3) proteins, which mutate viral nucleic acids. Inside the cell, Vif binds to the host cell proteins Elongin-C, Elongin-B, and core-binding factor subunit β, forming a four-protein complex called VCBC. The structure of VCBC-Cullin5 has recently been solved by X-ray crystallography, and, using molecular dynamics simulations, the dynamics of VCBC have been characterized. Here, we applied time-lapse high-speed atomic force microscopy to visualize the conformational changes of the VCBC complex. We determined the three most favorable conformations of this complex, which we identified as the triangle, dumbbell, and globular structures. Moreover, we characterized the dynamics of each of these structures. Our data revealed the very dynamic behavior of all of them, with the triangle and dumbbell structures being the most dynamic. These findings provide insight into the structure and dynamics of the VCBC complex and may support efforts to improve HIV treatment, because Vif is essential for virus survival in the cell.

Remodeling of the interdomain allosteric linker upon membrane binding of CCTα pulls its active site close to the membrane surface

The rate-limiting step in the biosynthesis of the major membrane phospholipid, phosphatidylcholine, is catalyzed by CTP:phosphocholine cytidylyltransferase (CCT), which is regulated by reversible membrane binding of a long amphipathic helix (domain M). The M domain communicates with the catalytic domain via a conserved ∼20-residue linker, essential for lipid activation of CCT. Previous analysis of this region (denoted as the αE_C/J) using MD simulations, cross-linking, mutagenesis, and solvent accessibility suggested that membrane binding of domain M promotes remodeling of the αE_C/J into a more compact structure that is required for enzyme activation. Here, using tryptophan fluorescence quenching, we show that the allosteric linker lies superficially on the membrane surface. Analyses with truncated CCTs show that the αE_C/J can interact with lipids independently of the M domain. We observed strong FRET between engineered tryptophans in the αE_C/J and vesicles containing dansyl-phosphatidylethanolamine that depended on the native J sequence. These data are incompatible with the extended conformation of the αE helix observed in the previously determined crystal structure of inactive CCT but support a bent αE helix conformation stabilized by J segment interactions. Our results suggest that the membrane-adsorbed, folded allosteric linker may partially cover the active site cleft and pull it close to the membrane surface, where cytidyl transfer can occur efficiently in a relatively anhydrous environment.

Conformational Dynamics of the HIV-Vif Protein Complex

Human immunodeficiency virus-1 viral infectivity factor (Vif) is an intrinsically disordered protein responsible for the ubiquitination of the APOBEC3 (A3) antiviral proteins. Vif folds when it binds Cullin-RING E3 ligase 5 and the transcription cofactor CBF-β. A five-protein complex containing the substrate receptor (Vif, CBF-β, Elongin-B, Elongin-C (VCBC)) and Cullin5 (CUL5) has a published crystal structure, but dynamics of this VCBC-CUL5 complex have not been characterized. Here, we use molecular dynamics (MD) simulations and NMR to characterize the dynamics of the VCBC complex with and without CUL5 and an A3 protein bound. Our simulations show that the VCBC complex undergoes global dynamics involving twisting and clamshell opening of the complex, whereas VCBC-CUL5 maintains a more static conformation, similar to the crystal structure. This observation from MD is supported by methyl-transverse relaxation-optimized spectroscopy NMR data, which indicates that the VCBC complex without CUL5 is dynamic on the μs-ms timescale. Our NMR data also show that the VCBC complex is more conformationally restricted when bound to the antiviral APOBEC3F (one of the A3 proteins), consistent with our MD simulations. Vif contains a flexible linker region located at the hinge of the VCBC complex, which changes conformation in conjunction with the global dynamics of the complex. Like other substrate receptors, VCBC can exist alone or in complex with CUL5 and other proteins in cells. Accordingly, the VCBC complex could be a good target for therapeutics that would inhibit full assembly of the ubiquitination complex by stabilizing an alternate VCBC conformation.

Role of Polyubiquitin Chain Flexibility in the Catalytic Mechanism of Cullin-RING Ubiquitin Ligases

Cullin-RING ubiquitin ligases are a diverse family of ubiquitin ligases that catalyze the synthesis of K48-linked polyubiquitin (polyUb) chains on a variety of substrates, ultimately leading to their degradation by the proteasome. The cullin-RING enzyme scaffold processively attaches a Ub molecule to the distal end of a growing chain up to lengths of eight Ub monomers. However, the molecular mechanism governing how chains of increasing size are built using a scaffold of largely fixed dimensions is not clear. We developed coarse-grained molecular dynamics simulations to describe the dependence of k_cat for cullin-RING ligases on the length and flexibility of the K48-linked polyUb chain attached to the substrate protein, key factors that determine the rate of subsequent Ub attachment to the chain, and therefore, the ensuing biological outcomes of ubiquitination. The results suggest that a number of regulatory mechanisms may lead to variations in the rate of chain elongation for different cullin-RING ligases. Specifically, modulation of the distance between the target lysine and the phosphodegron sequence of the substrate, the distance between the substrate lysine and the active site cysteine of the Ub conjugation enzyme (E2) bound to the cullin-RING scaffold, and flexibility of the bound E2 can lead to significant differences in the processing of K48-linked chains on substrates, potentially leading to differences in biological outcomes.

Cullin neddylation may allosterically tune polyubiquitin chain length and topology

Conjugation of Nedd8 (neddylation) to Cullins (Cul) in Cul-RING E3 ligases (CRLs) stimulates ubiquitination and polyubiquitination of protein substrates. CRL is made up of two Cul-flanked arms: one consists of the substrate-binding and adaptor proteins and the other consists of E2 and Ring-box protein (Rbx). Polyubiquitin chain length and topology determine the substrate fate. Here, we ask how polyubiquitin chains are accommodated in the limited space available between the two arms and what determines the polyubiquitin linkage topology. We focus on Cul5 and Rbx1 in three states: before Cul5 neddylation (closed state), after neddylation (open state), and after deneddylation, exploiting molecular dynamics simulations and the Gaussian Network Model. We observe that regulation of substrate ubiquitination and polyubiquitination takes place through Rbx1 rotations, which are controlled by Nedd8-Rbx1 allosteric communication. Allosteric propagation proceeds from Nedd8 via Cul5 dynamic hinges and hydrogen bonds between the C-terminal domain of Cul5 (Cul5CTD) and Rbx1 (Cul5CTD residues R538/R569 and Rbx1 residue E67, or Cul5CTD E474/E478/N491 and Rbx1 K105). Importantly, at each ubiquitination step (homogeneous or heterogeneous, linear or branched), the polyubiquitin linkages fit into the distances between the two arms, and these match the inherent CRL conformational tendencies. Hinge sites may constitute drug targets.

Crystal structure of pentapeptide-independent chemotaxis receptor methyltransferase (CheR) reveals idiosyncratic structural determinants for receptor recognition

Chemotactic methyltransferase, CheR catalyse methylation of specific glutamate residues in the cytoplasmic domain of methyl-accepting chemotactic protein receptors (MCPRs). The methylation of MCPRs is essential for the chemical sensing and chemotactic bacterial mobility towards favorable chemicals or away from unfavorable ones. In this study, crystal structure of B. subtilis CheR (BsCheR) in complex with S-adenosyl-l-homocysteine (SAH) has been determined to 1.8Å resolution. This is the first report of crystal structure belonging to the pentapeptide-independent CheR (PICheR) class. Till date, only one crystal structure of CheR from S. typhimurium (StCheR) belonging to pentapeptide-dependent CheR (PDCheR) class is available. Structural analysis of BsCheR reveals a helix-X-helix motif (HXH) with Asp53 as the linker residue in the N-terminal domain. The key structural features of the PDCheR β-subdomain involved in the formation of a tight complex with the pentapeptide binding motif in MCPRs were found to be absent in the structure of BsCheR. Additionally, isothermal titration calorimetry (ITC) experiments were performed to investigate S-adenosyl-(l)-methionine (SAM) binding affinity and K_D was determined to be 0.32mM. The structure of BsCheR reveals that mostly residues of the large C-terminal domain contribute to SAH binding, with contributions of few residues from the linker region and the N-terminal domain. Structural investigations and sequence analysis carried out in this study provide critical insights into the distinct receptor recognition mechanism of the PDCheR and PICheR methyltransferase classes.

Molecular dynamics simulations elucidate the mode of protein recognition by Skp1 and the F-box domain in the SCF complex

Polyubiquitination of the target protein by a ubiquitin transferring machinery is key to various cellular processes. E3 ligase Skp1-Cul1-F-box (SCF) is one such complex which plays crucial role in substrate recognition and transfer of the ubiquitin molecule. Previous computational studies have focused on S-phase kinase-associated protein 2 (Skp2), cullin, and RING-finger proteins of this complex, but the roles of the adapter protein Skp1 and F-box domain of Skp2 have not been determined. Using sub-microsecond molecular dynamics simulations of full-length Skp1, unbound Skp2, Skp2-Cks1 (Cks1: Cyclin-dependent kinases regulatory subunit 1), Skp1-Skp2, and Skp1-Skp2-Cks1 complexes, we have elucidated the function of Skp1 and the F-box domain of Skp2. We found that the L16 loop of Skp1, which was deleted in previous X-ray crystallography studies, can offer additional stability to the ternary complex via its interactions with the C-terminal tail of Skp2. Moreover, Skp1 helices H6, H7, and H8 display vivid conformational flexibility when not bound to Skp2, suggesting that these helices can recognize and lock the F-box proteins. Furthermore, we observed that the F-box domain could rotate (5°-129°), and that the binding partner determined the degree of conformational flexibility. Finally, Skp1 and Skp2 were found to execute a domain motion in Skp1-Skp2 and Skp1-Skp2-Cks1 complexes that could decrease the distance between ubiquitination site of the substrate and the ubiquitin molecule by 3 nm. Thus, we propose that both the F-box domain of Skp2 and Skp1-Skp2 domain motions displaying preferential conformational control can together facilitate polyubiquitination of a wide variety of substrates.

Serendipitous SAD Solution for DMSO-Soaked SOCS2-ElonginC-ElonginB Crystals Using Covalently Incorporated Dimethylarsenic: Insights into Substrate Receptor Conformational Flexibility in Cullin RING Ligases

Suppressor of cytokine signalling 2 (SOCS2) is the substrate-binding component of a Cullin-RING E3 ubiquitin ligase (CRL) complex that targets phosphorylated hormone receptors for degradation by the ubiquitin-proteasome system. As a key regulator of the transcriptional response to growth signals, SOCS2 and its protein complex partners are potential targets for small molecule development. We found that crystals of SOCS2 in complex with its adaptor proteins, Elongin C and Elongin B, underwent a change in crystallographic parameters when treated with dimethyl sulfoxide during soaking experiments. To solve the phase problem for the new crystal form we identified the presence of arsenic atoms in the crystals, a result of covalent modification of cysteines by cacodylate, and successfully extracted anomalous signal from these atoms for experimental phasing. The resulting structure provides a means for solving future structures where the crystals must be treated with DMSO for ligand soaking approaches. Additionally, the conformational changes induced in this structure reveal flexibility within SOCS2 that match those postulated by previous molecular dynamics simulations. This conformational flexibility illustrates how SOCS2 can orient its substrates for successful ubiquitination by other elements of the CRL complex.

Principles of allosteric interactions in cell signaling

Linking cell signaling events to the fundamental physicochemical basis of the conformational behavior of single molecules and ultimately to cellular function is a key challenge facing the life sciences. Here we outline the emerging principles of allosteric interactions in cell signaling, with emphasis on the following points. (1) Allosteric efficacy is not a function of the chemical composition of the allosteric pocket but reflects the extent of the population shift between the inactive and active states. That is, the allosteric effect is determined by the extent of preferred binding, not by the overall binding affinity. (2) Coupling between the allosteric and active sites does not decide the allosteric effect; however, it does define the propagation pathways, the allosteric binding sites, and key on-path residues. (3) Atoms of allosteric effectors can act as "driver" or "anchor" and create attractive "pulling" or repulsive "pushing" interactions. Deciphering, quantifying, and integrating the multiple co-occurring events present daunting challenges to our scientific community.

Multiple conformational selection and induced fit events take place in allosteric propagation

The fact that we observe a single conformational selection event during binding does not necessarily mean that only a single conformational selection event takes place, even though this is the common assumption. Here we suggest that conformational selection takes place not once in a given binding/allosteric event, but at every step along the allosteric pathway. This view generalizes conformational selection and makes it applicable also to other allosteric events, such as post-translational modifications (PTMs) and photon absorption. Similar to binding, at each step along a propagation pathway, conformational selection is coupled with induced fit which optimizes the interactions. Thus, as in binding, the allosteric effects induced by PTMs and light relate not only to population shift; but to conformational selection as well. Conformational selection and population shift take place conjointly.

The isolated N terminus of Ring1B is a well-folded, monomeric fragment with native-like structure

The Polycomb group (PcG) proteins assemble into Polycomb repressive complexes (PRCs), PRC1 and PRC2, which act as general transcriptional repressors. PRC1 comprises a variety of biochemical entities endowed with histone H2A monoubiquitylation activity conferred by really interesting new gene (RING) finger E3 ubiquitin ligases Ring1A and Ring1B. All PRC1 complexes contain Ring1 proteins which are essential for Polycomb epigenetic regulation. We have been able to express the isolated N-terminal region of Ring1B, N-Ring1B, comprising the first 221 residues of the 334-residue-long Ring1B. This fragment contains the 41-residue-long RING finger motif, and flanking sequences that form an interacting platform for PcG and non-PcG proteins. We found that the N-Ring1B is a well-folded, monomeric fragment, with native-like structure which unfolds irreversibly. The protein is capable of binding to an ubiquitin-conjugase protein (with an 85% of sequence similarity to the Ring1B physiological partner) with moderate affinity.

The spatial structure of cell signaling systems

The spatial structure of the cell is highly organized at all levels: from small complexes and assemblies, to local nano- and microclusters, to global, micrometer scales across and between cells. We suggest that this multiscale spatial cell organization also organizes signaling and coordinates cellular behavior. We propose a new view of the spatial structure of cell signaling systems. This new view describes cell signaling in terms of dynamic allosteric interactions within and among distinct, spatially organized transient clusters. The clusters vary over time and space and are on length scales from nanometers to micrometers. When considered across these length scales, primary factors in the spatial organization are cell membrane domains and the actin cytoskeleton, both also highly dynamic. A key challenge is to understand the interplay across these multiple scales, link it to the physicochemical basis of the conformational behavior of single molecules and ultimately relate it to cellular function. Overall, our premise is that at these scales, cell signaling should be thought of not primarily as a sequence of diffusion-controlled molecular collisions, but instead transient, allostery-driven cluster re-forming interactions.

Dispersed sites of HIV Vif-dependent polyubiquitination in the DNA deaminase APOBEC3F

APOBEC3F (A3F) and APOBEC3G (A3G) are DNA cytosine deaminases that potently restrict human immunodeficiency virus type 1 replication when the virus is deprived of its accessory protein Vif (virion infectivity factor). Vif counteracts these restriction factors by recruiting A3F and A3G to an E3 ubiquitin (Ub) ligase complex that mediates their polyubiquitination (polyUb) and proteasomal degradation. While previous efforts have identified single amino acid residues in APOBEC3 proteins required for Vif recognition, less is known about the downstream Ub acceptor sites that are targeted. One prior report identified a cluster of polyubiquitinated residues in A3G and proposed an antiparallel model of A3G interaction with the Vif-E3 Ub ligase complex wherein Vif binding at one terminus of A3G orients the opposite terminus for polyUb [Iwatani et al. (2009). Proc. Natl. Acad. Sci. USA, 106, 19539-19544]. To test the generalizability of this model, we carried out a complete mutagenesis of the lysine residues in A3F and used a complementary, unbiased proteomic approach to identify Ub acceptor sites targeted by Vif. Our data indicate that internal lysines are the dominant Ub acceptor sites in both A3F and A3G. In contrast with the proposed antiparallel model, however, we find that the Vif-dependent polyUb of A3F and A3G can occur at multiple acceptor sites dispersed along predicted lysine-enriched surfaces of both the N- and C-terminal deaminase domains. These data suggest an alternative model for binding of APOBEC3 proteins to the Vif-E3 Ub ligase complex and diminish enthusiasm for the amenability of APOBEC3 Ub acceptor sites to therapeutic intervention.

The role of allostery in the ubiquitin-proteasome system

The ubiquitin-proteasome system (UPS) is involved in many cellular processes including protein degradation. Degradation of a protein via this system involves two successive steps: ubiquitination and degradation. Ubiquitination tags the target protein with ubiquitin-like proteins (UBLs), such as ubiquitin, small ubiquitin-like modifier (SUMO) and NEDD8, via a cascade involving three enzymes: activating enzyme E1, conjugating enzyme E2 and E3 ubiquitin ligases. The proteasomes recognize the UBL-tagged substrate proteins and degrade them. Accumulating evidence indicates that allostery is a central player in the regulation of ubiquitination, as well as deubiquitination and degradation. Here, we provide an overview of the key mechanistic roles played by allostery in all steps of these processes, and highlight allosteric drugs targeting them. Throughout the review, we emphasize the crucial mechanistic role played by linkers in allosterically controlling the UPS action by biasing the sampling of the conformational space, which facilitate the catalytic reactions of the ubiquitination and degradation. Finally, we propose that allostery may similarly play key roles in the regulation of molecular machines in the cell, and as such allosteric drugs can be expected to be increasingly exploited in therapeutic regimes.

The underappreciated role of allostery in the cellular network

Allosteric propagation results in communication between distinct sites in the protein structure; it also encodes specific effects on cellular pathways, and in this way it shapes cellular response. One example of long-range effects is binding of morphogens to cell surface receptors, which initiates a cascade of protein interactions that leads to genome activation and specific cellular action. Allosteric propagation results from combinations of multiple factors, takes place through dynamic shifts of conformational ensembles, and affects the equilibria of macromolecular interactions. Here, we (a) emphasize the well-known yet still underappreciated role of allostery in conveying explicit signals across large multimolecular assemblies and distances to specify cellular action; (b) stress the need for quantitation of the allosteric effects; and finally, (c) propose that each specific combination of allosteric effectors along the pathway spells a distinct function. The challenges are colossal; the inspiring reward will be predicting function, misfunction, and outcomes of drug regimes.

A broad view of scaffolding suggests that scaffolding proteins can actively control regulation and signaling of multienzyme complexes through allostery

Enzymes often work sequentially in pathways; and consecutive reaction steps are typically carried out by molecules associated in the same multienzyme complex. Localization confines the enzymes; anchors them; increases the effective concentration of substrates and products; and shortens pathway timescales; however, it does not explain enzyme coordination or pathway branching. Here, we distinguish between metabolic and signaling multienzyme complexes. We argue for a central role of scaffolding proteins in regulating multienzyme complexes signaling and suggest that metabolic multienzyme complexes are less dependent on scaffolding because they undergo conformational control through direct subunit-subunit contacts. In particular, we propose that scaffolding proteins have an essential function in controlling branching in signaling pathways. This new broadened definition of scaffolding proteins goes beyond cases such as the classic yeast mitogen-activated protein kinase Ste5 and encompasses proteins such as E3 ligases which lack active sites and work via allostery. With this definition, we classify the mechanisms of multienzyme complexes based on whether the substrates are transferred through the involvement of scaffolding proteins, and outline the functional merits to metabolic or signaling pathways. Overall, while co-localization topography helps multistep pathways non-specifically, allosteric regulation requires precise multienzyme organization and interactions and works via population shift, either through direct enzyme subunit-subunit interactions or through active involvement of scaffolding proteins. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

G503 is obligatory for coupling of regulatory domains in NCX proteins

In multidomain proteins, interdomain linkers allow an efficient transfer of regulatory information, although it is unclear how the information encoded in the linker structure coins dynamic coupling. Allosteric regulation of NCX proteins involves Ca(2+)-driven tethering of regulatory CBD1 and CBD2 (through a salt bridge network) accompanied by alignment of CBDs and Ca(2+) occlusion at the interface of the two CBDs. Here we investigated "alanine-walk" substitutions in the CBD1-CBD2 linker (501-HAGIFT-506) and found that among all linker residues, only G503 is obligatory for Ca(2+)-induced reorientations of CBDs and slow dissociation of occluded Ca(2+). Moreover, swapping between positions A502 and G503 in the CBD1-CBD2 linker results in a complete loss of slow dissociation of occluded Ca(2+), meaning that dynamic coupling of CBDs requires an exact pose of glycine at position 503. Therefore, accumulating data revealed that position 503 occupied by glycine is absolutely required for Ca(2+)-driven tethering of CBDs, which in turn limits the linker's flexibility and, thus, restricts CBD movements. Because G503 is extremely well conserved in eukaryotic NCX proteins, the information encoded in G503 is essential for dynamic coupling of the two-domain CBD tandem and, thus, for propagation of the allosteric signal.

Loop 7 of E2 enzymes: an ancestral conserved functional motif involved in the E2-mediated steps of the ubiquitination cascade

The ubiquitin (Ub) system controls almost every aspect of eukaryotic cell biology. Protein ubiquitination depends on the sequential action of three classes of enzymes (E1, E2 and E3). E2 Ub-conjugating enzymes have a central role in the ubiquitination pathway, interacting with both E1 and E3, and influencing the ultimate fate of the substrates. Several E2s are characterized by an extended acidic insertion in loop 7 (L7), which if mutated is known to impair the proper E2-related functions. In the present contribution, we show that acidic loop is a conserved ancestral motif in E2s, relying on the presence of alternate hydrophobic and acidic residues. Moreover, the dynamic properties of a subset of family 3 E2s, as well as their binary and ternary complexes with Ub and the cognate E3, have been investigated. Here we provide a model of L7 role in the different steps of the ubiquitination cascade of family 3 E2s. The L7 hydrophobic residues turned out to be the main determinant for the stabilization of the E2 inactive conformations by a tight network of interactions in the catalytic cleft. Moreover, phosphorylation is known from previous studies to promote E2 competent conformations for Ub charging, inducing electrostatic repulsion and acting on the L7 acidic residues. Here we show that these active conformations are stabilized by a network of hydrophobic interactions between L7 and L4, the latter being a conserved interface for E3-recruitment in several E2s. In the successive steps, L7 conserved acidic residues also provide an interaction interface for both Ub and the Rbx1 RING subdomain of the cognate E3. Our data therefore suggest a crucial role for L7 of family 3 E2s in all the E2-mediated steps of the ubiquitination cascade. Its different functions are exploited thank to its conserved hydrophobic and acidic residues in a finely orchestrate mechanism.

Flexible cullins in cullin-RING E3 ligases allosterically regulate ubiquitination

How do the cullins, with conserved structures, accommodate substrate-binding proteins with distinct shapes and sizes? Cullin-RING E3 ubiquitin ligases facilitate ubiquitin transfer from E2 to the substrate, tagging the substrate for degradation. They contain substrate-binding, adaptor, cullin, and Rbx proteins. Previously, we showed that substrate-binding and Rbx proteins are flexible. This allows shortening of the E2-substrate distance for initiation of ubiquitination or increasing the distance to accommodate the polyubiquitin chain. However, the role of the cullin remained unclear. Is cullin a rigid scaffold, or is it flexible and actively assists in the ubiquitin transfer reaction? Why are there different cullins, and how do these cullins specifically facilitate ubiquitination for different substrates? To answer these questions, we performed structural analysis and molecular dynamics simulations based on Cul1, Cul4A, and Cul5 crystal structures. Our results show that these three cullins are not rigid scaffolds but are flexible with conserved hinges in the N-terminal domain. However, the degrees of flexibilities are distinct among the different cullins. Of interest, Cul1 flexibility can also be changed by deletion of the long loop (which is absent in Cul4A) in the N-terminal domain, suggesting that the loop may have an allosteric functional role. In all three cases, these conformational changes increase the E2-substrate distance to a specific range to facilitate polyubiquitination, suggesting that rather than being inert scaffold proteins, cullins allosterically regulate ubiquitination.

Dynamic allostery: linkers are not merely flexible

Most proteins consist of multiple domains. How do linkers efficiently transfer information between sites that are on different domains to activate the protein? Mere flexibility only implies that the conformations would be sampled. For fast timescales between triggering events and cellular response, which often involves large conformational change, flexibility on its own may not constitute a good solution. We posit that successive conformational states along major allosteric propagation pathways are pre-encoded in linker sequences where each state is encoded by the previous one. The barriers between these states that are hierarchically populated are lower, achieving faster timescales even for large conformational changes. We further propose that evolution has optimized the linker sequences and lengths for efficiency, which explains why mutations in linkers may affect protein function and review the literature in this light.

Expanding the proteome: disordered and alternatively folded proteins

Proteins provide much of the scaffolding for life, as well as undertaking a variety of essential catalytic reactions. These characteristic functions have led us to presuppose that proteins are in general functional only when well structured and correctly folded. As we begin to explore the repertoire of possible protein sequences inherent in the human and other genomes, two stark facts that belie this supposition become clear: firstly, the number of apparent open reading frames in the human genome is significantly smaller than appears to be necessary to code for all of the diverse proteins in higher organisms, and secondly that a significant proportion of the protein sequences that would be coded by the genome would not be expected to form stable three-dimensional (3D) structures. Clearly the genome must include coding for a multitude of alternative forms of proteins, some of which may be partly or fully disordered or incompletely structured in their functional states. At the same time as this likelihood was recognized, experimental studies also began to uncover examples of important protein molecules and domains that were incompletely structured or completely disordered in solution, yet remained perfectly functional. In the ensuing years, we have seen an explosion of experimental and genome-annotation studies that have mapped the extent of the intrinsic disorder phenomenon and explored the possible biological rationales for its widespread occurrence. Answers to the question 'why would a particular domain need to be unstructured?' are as varied as the systems where such domains are found. This review provides a survey of recent new directions in this field, and includes an evaluation of the role not only of intrinsically disordered proteins but also of partially structured and highly dynamic members of the disorder-order continuum.

A comparison of the dynamics of pantothenate synthetase from M. tuberculosis and E. coli: computational studies

Pantothenate synthetase (PS) catalyzes the final step of the pantothenate pathway, in which pantothenate is formed from pantoate and β-alanine in an ATP-dependent reaction. Mycobacterium tuberculosis PS (MTB PS) is functionally a dimer and a potential target for novel antitubercular drugs. Molecular dynamics simulations show that the functional dynamics of the enzyme are dominated by motions of a flexible gate loop in the N-terminal domain and of the C-terminal domain. The gate loop motions dominate in MTB PS while the C-terminal domain motion dominates in Escherichia coli PS. Simulations also show that the correlated motions of the domains are severely compromised in the monomeric forms. Mutations that reduce the mobility of the gate loop in MTB PS and increased it in E. coli PS were designed and validated through simulations.

Rbx1 flexible linker facilitates cullin-RING ligase function before neddylation and after deneddylation

In ubiquitination, cullin-RING E3 ubiquitin ligases (CRLs) assist in ubiquitin transfer from ubiquitin-conjugating enzyme E2 to the substrate. Neddylation, which involves NEDD8 transfer from E2 to E3-cullin, stimulates ubiquitination by inducing conformational change in CRLs. However, deneddylation, which removes NEDD8 from cullin, does not suppress ubiquitination in vivo, raising the question of how neddylation/deneddylation exerts its effects. Using molecular-dynamics simulations, we demonstrate that before neddylation occurs, the linker flexibility of Rbx1, a CRL component, leads to conformational changes in CRLs that allow neddylation and initiation of ubiquitination. These large NEDD8-induced conformational changes are retained after deneddylation, allowing both initiation of the ubiquitination process and ubiquitin chain elongation after deneddylation. Furthermore, mutation of lysine, the cullin residue to which NEDD8 covalently attaches, dramatically reduces CRL conformational changes, suggesting that the acceptor lysine allosterically regulates CRLs. Thus, our results imply that neddylation stimulates ubiquitination by CRL conformational control via lysine modification.