IDPConformerGenerator: A Flexible Software Suite for Sampling the Conformational Space of Disordered Protein States.

The power of structural information for informing biological mechanisms is clear for stable folded macromolecules, but similar structure-function insight is more difficult to obtain for highly dynamic systems such as intrinsically disordered proteins (IDPs) which must be described as structural ensembles. Here, we present IDPConformerGenerator, a flexible, modular open-source software platform for generating large and diverse ensembles of disordered protein states that builds conformers that obey geometric, steric, and other physical restraints on the input sequence. IDPConformerGenerator samples backbone phi (φ), psi (ψ), and omega (ω) torsion angles of relevant sequence fragments from loops and secondary structure elements extracted from folded protein structures in the RCSB Protein Data Bank and builds side chains from robust Monte Carlo algorithms using expanded rotamer libraries. IDPConformerGenerator has many user-defined options enabling variable fractional sampling of secondary structures, supports Bayesian models for assessing the agreement of IDP ensembles for consistency with experimental data, and introduces a machine learning approach to transform between internal and Cartesian coordinates with reduced error. IDPConformerGenerator will facilitate the characterization of disordered proteins to ultimately provide structural insights into these states that have key biological functions.

Identification of a molecular locus for normalizing dysregulated GABA release from interneurons in the Fragile X brain.

Principal neurons encode information by varying their firing rate and patterns precisely fine-tuned through GABAergic interneurons. Dysregulation of inhibition can lead to neuropsychiatric disorders, yet little is known about the molecular basis underlying inhibitory control. Here, we find that excessive GABA release from basket cells (BCs) attenuates the firing frequency of Purkinje neurons (PNs) in the cerebellum of Fragile X Mental Retardation 1 (Fmr1) knockout (KO) mice, a model of Fragile X Syndrome (FXS) with abrogated expression of the Fragile X Mental Retardation Protein (FMRP). This over-inhibition originates from increased excitability and Ca2+ transients in the presynaptic terminals, where Kv1.2 potassium channels are downregulated. By paired patch-clamp recordings, we further demonstrate that acutely introducing an N-terminal fragment of FMRP into BCs normalizes GABA release in the Fmr1-KO synapses. Conversely, direct injection of an inhibitory FMRP antibody into BCs, or membrane depolarization of BCs, enhances GABA release in the wild type synapses, leading to abnormal inhibitory transmission comparable to the Fmr1-KO neurons. We discover that the N-terminus of FMRP directly binds to a phosphorylated serine motif on the C-terminus of Kv1.2; and that loss of this interaction in BCs exaggerates GABA release, compromising the firing activity of PNs and thus the output from the cerebellar circuitry. An allosteric Kv1.2 agonist, docosahexaenoic acid, rectifies the dysregulated inhibition in vitro as well as acoustic startle reflex and social interaction in vivo of the Fmr1-KO mice. Our results unravel a novel molecular locus for targeted intervention of FXS and perhaps autism.

Folding of an intrinsically disordered protein by phosphorylation as a regulatory switch.

Intrinsically disordered proteins play important roles in cell signalling, transcription, translation and cell cycle regulation. Although they lack stable tertiary structure, many intrinsically disordered proteins undergo disorder-to-order transitions upon binding to partners. Similarly, several folded proteins use regulated order-to-disorder transitions to mediate biological function. In principle, the function of intrinsically disordered proteins may be controlled by post-translational modifications that lead to structural changes such as folding, although this has not been observed. Here we show that multisite phosphorylation induces folding of the intrinsically disordered 4E-BP2, the major neural isoform of the family of three mammalian proteins that bind eIF4E and suppress cap-dependent translation initiation. In its non-phosphorylated state, 4E-BP2 interacts tightly with eIF4E using both a canonical YXXXXLΦ motif (starting at Y54) that undergoes a disorder-to-helix transition upon binding and a dynamic secondary binding site. We demonstrate that phosphorylation at T37 and T46 induces folding of residues P18-R62 of 4E-BP2 into a four-stranded ß-domain that sequesters the helical YXXXXLΦ motif into a partly buried ß-strand, blocking its accessibility to eIF4E. The folded state of pT37pT46 4E-BP2 is weakly stable, decreasing affinity by 100-fold and leading to an order-to-disorder transition upon binding to eIF4E, whereas fully phosphorylated 4E-BP2 is more stable, decreasing affinity by a factor of approximately 4,000. These results highlight stabilization of a phosphorylation-induced fold as the essential mechanism for phospho-regulation of the 4E-BP:eIF4E interaction and exemplify a new mode of biological regulation mediated by intrinsically disordered proteins.

Conformational Ensembles of an Intrinsically Disordered Protein Consistent with NMR, SAXS, and Single-Molecule FRET.

Intrinsically disordered proteins (IDPs) have fluctuating heterogeneous conformations, which makes their structural characterization challenging. Although challenging, characterization of the conformational ensembles of IDPs is of great interest, since their conformational ensembles are the link between their sequences and functions. An accurate description of IDP conformational ensembles depends crucially on the amount and quality of the experimental data, how it is integrated, and if it supports a consistent structural picture. We used integrative modeling and validation to apply conformational restraints and assess agreement with the most common structural techniques for IDPs: Nuclear Magnetic Resonance (NMR) spectroscopy, Small-angle X-ray Scattering (SAXS), and single-molecule Förster Resonance Energy Transfer (smFRET). Agreement with such a diverse set of experimental data suggests that details of the generated ensembles can now be examined with a high degree of confidence. Using the disordered N-terminal region of the Sic1 protein as a test case, we examined relationships between average global polymeric descriptions and higher-moments of their distributions. To resolve apparent discrepancies between smFRET and SAXS inferences, we integrated SAXS data with NMR data and reserved the smFRET data for independent validation. Consistency with smFRET, which was not guaranteed a priori, indicates that, globally, the perturbative effects of NMR or smFRET labels on the Sic1 ensemble are minimal. Analysis of the ensembles revealed distinguishing features of Sic1, such as overall compactness and large end-to-end distance fluctuations, which are consistent with biophysical models of Sic1's ultrasensitive binding to its partner Cdc4. Our results underscore the importance of integrative modeling and validation in generating and drawing conclusions from IDP conformational ensembles.

Role of CBS and Bateman Domains in Phosphorylation-Dependent Regulation of a CLC Anion Channel.

Eukaryotic CLC anion channels and transporters are homodimeric proteins composed of multiple α-helical membrane domains and large cytoplasmic C-termini containing two cystathionine-ß-synthase domains (CBS1 and CBS2) that dimerize to form a Bateman domain. The Bateman domains of adjacent CLC subunits interact to form a Bateman domain dimer. The functions of CLC CBS and Bateman domains are poorly understood. We utilized the Caenorhabditis elegans CLC-1/2/Ka/Kb anion channel homolog CLH-3b to characterize the regulatory roles of CLC cytoplasmic domains. CLH-3b activity is reduced by phosphorylation or deletion of a 14-amino-acid activation domain (AD) located on the linker connecting CBS1 and CBS2. We demonstrate here that phosphorylation-dependent reductions in channel activity require an intact Bateman domain dimer and concomitant phosphorylation or deletion of both ADs. Regulation of a CLH-3b AD deletion mutant is reconstituted by intracellular perfusion with recombinant 14-amino-acid AD peptides. The sulfhydryl reactive reagent 2-(trimethylammonium)ethyl methanethiosulfonate bromide (MTSET) alters in a phosphorylation-dependent manner the activity of channels containing single cysteine residues that are engineered into the short intracellular loop connecting membrane α-helices H and I (H-I loop), the AD, CBS1, and CBS2. In contrast, MTSET has no effect on channels in which cysteine residues are engineered into intracellular regions that are dispensable for regulation. These studies together with our previous work suggest that binding and unbinding of the AD to the Bateman domain dimer induces conformational changes that are transduced to channel membrane domains via the H-I loop. Our findings provide new, to our knowledge, insights into the roles of CLC Bateman domains and the structure-function relationships that govern the regulation of CLC protein activity by diverse ligands and signaling pathways.

Regulatory R region of the CFTR chloride channel is a dynamic integrator of phospho-dependent intra- and intermolecular interactions.

Intrinsically disordered proteins play crucial roles in regulatory processes and often function as protein interaction hubs. Here, we present a detailed characterization of a full-length disordered hub protein region involved in multiple dynamic complexes. We performed NMR, CD, and fluorescence binding studies on the nonphosphorylated and highly PKA-phosphorylated human cystic fibrosis transmembrane conductance regulator (CFTR) regulatory region, a ∼200-residue disordered segment involved in phosphorylation-dependent regulation of channel trafficking and gating. Our data provide evidence for dynamic, phosphorylation-dependent, multisite interactions of various segments of the regulatory region for its intra- and intermolecular partners, including the CFTR nucleotide binding domains 1 and 2, a 42-residue peptide from the C terminus of CFTR, the SLC26A3 sulphate transporter and antisigma factor antagonist (STAS) domain, and 14-3-3ß. Because of its large number of binding partners, multivalent binding of individually weak sites facilitates rapid exchange between free and bound states to allow the regulatory region to engage with different partners and generate a graded or rheostat-like response to phosphorylation. Our results enrich the understanding of how disordered binding segments interact with multiple targets. We present structural models consistent with our data that illustrate this dynamic aspect of phospho-regulation of CFTR by the disordered regulatory region.

Characterization of disordered proteins with ENSEMBLE.

UNLABELLED: ENSEMBLE is a computational approach for determining a set of conformations that represents the structural ensemble of a disordered protein based on input experimental data. The disordered protein can be an unfolded or intrinsically disordered state. Here, we introduce the latest version of the program, which has been enhanced to facilitate its general release and includes an intuitive user interface, as well as new approaches to treat data and analyse results. AVAILABILITY AND IMPLEMENTATION: ENSEMBLE is a program implemented in C and embedded in a Perl wrapper. It is supported on main Linux distributions. Source codes and installation files, including a detailed example, can be freely downloaded at http://abragam.med.utoronto.ca/∼JFKlab.

Human galectin-3 (Mac-2 antigen): defining molecular switches of affinity to natural glycoproteins, structural and dynamic aspects of glycan binding by flexible ligand docking and putative regulatory sequences in the proximal promoter region.

BACKGROUND: Human galectin-3 (Mac-2 antigen) is a cell-type-specific multifunctional effector owing to selective binding of distinct cell-surface glycoconjugates harboring ß-galactosides. The structural basis underlying the apparent preferences for distinct glycoproteins and for expression is so far unknown. METHODS: We strategically combined solid-phase assays on 43 natural glycoproteins with a new statistical approach to fully flexible computational docking and also processed the proximal promoter region in silico. RESULTS: The degree of branching in N-glycans and clustering of core 1 O-glycans are positive modulators for avidity. Sialylation of N-glycans in α2-6 linkage and of core 1 O-glycans in α2-3 linkage along with core 2 branching was an unfavorable factor, despite the presence of suited glycans in the vicinity. The lectin-ligand contact profile was scrutinized for six natural di- and tetrasaccharides enabling a statistical grading by analyzing flexible docking trajectories. The computational analysis of the proximal promoter region delineated putative sites for Lmo2/c-Ets-1 binding and new sites with potential for RUNX binding. GENERAL SIGNIFICANCE: These results identify new features of glycan selectivity and ligand contact by combining solid-phase assays with in silico work as well as of reactivity potential of the promoter.

Toward bioinspired galectin mimetics: identification of ligand-contacting peptides by proteolytic-excision mass spectrometry.

Clinically relevant bioactivities of human galectins (adhesion/growth-regulatory galactoside-specific lectins) inspired the design of peptides as new tools to elicit favorable effects (e.g., in growth control) or block harmful binding (e.g., in tissue invasion). To obtain the bioinspired lead compounds, we combined a proteolytic fragmentation approach without/with ligand contact (excision) with mass spectrometric identification of affinity-bound protein fragments, using galectin-1 and -3 as models. Two peptides from the carbohydrate recognition domains were obtained in each case in experimental series rigorously controlled for specificity, and the [157-162] peptide of galectin-3 proved to be active in blocking lectin binding to a neoglycoprotein and to tumor cell surfaces. This approach affords peptide sequences for structural optimization and intrafamily/phylogenetic galectin comparison at the binding-site level with a minimal requirement of protein quantity, and it is even amenable to mixtures.

Recessive NOS1AP variants impair actin remodeling and cause glomerulopathy in humans and mice.

Nephrotic syndrome (NS) is a leading cause of chronic kidney disease. We found recessive NOS1AP variants in two families with early-onset NS by exome sequencing. Overexpression of wild-type (WT) NOS1AP, but not cDNA constructs bearing patient variants, increased active CDC42 and promoted filopodia and podosome formation. Pharmacologic inhibition of CDC42 or its effectors, formin proteins, reduced NOS1AP-induced filopodia formation. NOS1AP knockdown reduced podocyte migration rate (PMR), which was rescued by overexpression of WT Nos1ap but not by constructs bearing patient variants. PMR in NOS1AP knockdown podocytes was also rescued by constitutively active CDC42Q61L or the formin DIAPH3 Modeling a NOS1AP patient variant in knock-in human kidney organoids revealed malformed glomeruli with increased apoptosis. Nos1apEx3-/Ex3- mice recapitulated the human phenotype, exhibiting proteinuria, foot process effacement, and glomerulosclerosis. These findings demonstrate that recessive NOS1AP variants impair CDC42/DIAPH-dependent actin remodeling, cause aberrant organoid glomerulogenesis, and lead to a glomerulopathy in humans and mice.

SAMPLEX: automatic mapping of perturbed and unperturbed regions of proteins and complexes.

BACKGROUND: The activity of proteins within the cell is characterized by their motions, flexibility, interactions or even the particularly intriguing case of partially unfolded states. In the last two cases, a part of the protein is affected either by binding or unfolding and the detection of the respective perturbed and unperturbed region(s) is a fundamental part of the structural characterization of these states. This can be achieved by comparing experimental data of the same protein in two different states (bound/unbound, folded/unfolded). For instance, measurements of chemical shift perturbations (CSPs) from NMR 1H-15N HSQC experiments gives an excellent opportunity to discriminate both moieties. RESULTS: We describe an innovative, automatic and unbiased method to distinguish perturbed and unperturbed regions in a protein existing in two distinct states (folded/partially unfolded, bound/unbound). The SAMPLEX program takes as input a set of data and the corresponding three-dimensional structure and returns the confidence for each residue to be in a perturbed or unperturbed state. Its performance is demonstrated for different applications including the prediction of disordered regions in partially unfolded proteins and of interacting regions in protein complexes. CONCLUSIONS: The proposed approach is suitable for partially unfolded states of proteins, local perturbations due to small ligands and protein-protein interfaces. The method is not restricted to NMR data, but is generic and can be applied to a wide variety of information.

Extended Experimental Inferential Structure Determination Method in Determining the Structural Ensembles of Disordered Protein States.

Proteins with intrinsic or unfolded state disorder comprise a new frontier in structural biology, requiring the characterization of diverse and dynamic structural ensembles. We introduce a comprehensive Bayesian framework, the Extended Experimental Inferential Structure Determination (X-EISD) method, that calculates the maximum log-likelihood of a disordered protein ensemble. X-EISD accounts for the uncertainties of a range of experimental data and back-calculation models from structures, including NMR chemical shifts, J-couplings, Nuclear Overhauser Effects (NOEs), paramagnetic relaxation enhancements (PREs), residual dipolar couplings (RDCs), hydrodynamic radii (R h ), single molecule fluorescence Förster resonance energy transfer (smFRET) and small angle X-ray scattering (SAXS). We apply X-EISD to the joint optimization against experimental data for the unfolded drkN SH3 domain and find that combining a local data type, such as chemical shifts or J-couplings, paired with long-ranged restraints such as NOEs, PREs or smFRET, yields structural ensembles in good agreement with all other data types if combined with representative IDP conformers.

MINOES: a new approach to select a representative ensemble of structures in NMR studies of (partially) unfolded states. Application to Delta25-PYP.

In nature, some proteins partially unfold under specific environmental conditions. These unfolded states typically consist of a large ensemble of conformations; their proper description is therefore a challenging problem. NMR spectroscopy is particularly well suited for this task: information on conformational preferences can be derived, for example, from chemical shifts or residual dipolar couplings. This information, which is measured as a time- and ensemble-average, can be used to model these states by generating large ensembles of conformations. The challenge is then to select a minimum representative set of conformations out of a large ensemble to represent the unfolded state. We have developed for this purpose an algorithm called MINOES (MINimum Optimal Ensemble Selection), which is based on an iterative procedure based on a driven expansion/contraction selection process. MINOES aims at selecting an optimal and minimal ensemble of conformations that, on average, maximizes the agreement between back-calculated and experimental (NMR) data, without any a-priori assumption about the required ensemble size. This approach is demonstrated by modeling the partially unfolded state of a deletion mutant of the Photoactive Yellow Protein, Delta25-PYP, which has been previously characterized by NMR (Bernard et al., Structure 2005;13:953-962).

HADDOCK versus HADDOCK: new features and performance of HADDOCK2.0 on the CAPRI targets.

Here we present version 2.0 of HADDOCK, which incorporates considerable improvements and new features. HADDOCK is now able to model not only protein-protein complexes but also other kinds of biomolecular complexes and multi-component (N > 2) systems. In the absence of any experimental and/or predicted information to drive the docking, HADDOCK now offers two additional ab initio docking modes based on either random patch definition or center-of-mass restraints. The docking protocol has been considerably improved, supporting among other solvated docking, automatic definition of semi-flexible regions, and inclusion of a desolvation energy term in the scoring scheme. The performance of HADDOCK2.0 is evaluated on the targets of rounds 4-11, run in a semi-automated mode using the original information we used in our CAPRI submissions. This enables a direct assessment of the progress made since the previous versions. Although HADDOCK performed very well in CAPRI (65% and 71% success rates, overall and for unbound targets only, respectively), a substantial improvement was achieved with HADDOCK2.0.

An evolutionary switch in ND2 enables Src kinase regulation of NMDA receptors.

The non-receptor tyrosine kinase Src is a key signalling hub for upregulating the function of N-methyl D-aspartate receptors (NMDARs). Src is anchored within the NMDAR complex via NADH dehydrogenase subunit 2 (ND2), a mitochondrially encoded adaptor protein. The interacting regions between Src and ND2 have been broadly identified, but the interaction between ND2 and the NMDAR has remained elusive. Here we generate a homology model of ND2 and dock it onto the NMDAR via the transmembrane domain of GluN1. This interaction is enabled by the evolutionary loss of three helices in bilaterian ND2 proteins compared to their ancestral homologues. We experimentally validate our model and demonstrate that blocking this interaction with an ND2 fragment identified in our experimental studies prevents Src-mediated upregulation of NMDAR currents in neurons. Our findings establish the mode of interaction between an NMDAR accessory protein with one of the core subunits of the receptor.

Structural changes of CFTR R region upon phosphorylation: a plastic platform for intramolecular and intermolecular interactions.

Chloride channel gating and trafficking of the cystic fibrosis transmembrane conductance regulator (CFTR) are regulated by phosphorylation. Intrinsically disordered segments of the protein are responsible for phospho-regulation, particularly the regulatory (R) region that is a target for several kinases and phosphatases. The R region remains disordered following phosphorylation, with different phosphorylation states sampling various conformations. Recent studies have demonstrated the crucial role that intramolecular and intermolecular interactions of the R region play in CFTR regulation. Different partners compete for the same binding segment, with the R region containing multiple overlapping binding elements. The non-phosphorylated R region interacts with the nucleotide binding domains and inhibits channel activity by blocking heterodimerization. Phosphorylation shifts the equilibrium such that the R region is excluded from the dimer interface, facilitating gating and processing by stimulating R region interactions with other domains and proteins. The dynamic conformational sampling and transient binding of the R region to multiple partners enables complex control of CFTR channel activity and trafficking.

Bioinformatics and functional analysis of an Entamoeba histolytica mannosyltransferase necessary for parasite complement resistance and hepatical infection.

The glycosylphosphatidylinositol (GPI) moiety is one of the ways by which many cell surface proteins, such as Gal/GalNAc lectin and proteophosphoglycans (PPGs) attach to the surface of Entamoeba histolytica, the agent of human amoebiasis. It is believed that these GPI-anchored molecules are involved in parasite adhesion to cells, mucus and the extracellular matrix. We identified an E. histolytica homolog of PIG-M, which is a mannosyltransferase required for synthesis of GPI. The sequence and structural analysis led to the conclusion that EhPIG-M1 is composed of one signal peptide and 11 transmembrane domains with two large intra luminal loops, one of which contains the DXD motif, involved in the enzymatic catalysis and conserved in most glycosyltransferases. Expressing a fragment of the EhPIG-M1 encoding gene in antisense orientation generated parasite lines diminished in EhPIG-M1 levels; these lines displayed reduced GPI production, were highly sensitive to complement and were dramatically inhibited for amoebic abscess formation. The data suggest a role for GPI surface anchored molecules in the survival of E. histolytica during pathogenesis.

Activity-structure correlations in divergent lectin evolution: fine specificity of chicken galectin CG-14 and computational analysis of flexible ligand docking for CG-14 and the closely related CG-16.

Gene duplication and sequence divergence are driving forces toward establishing protein families. To examine how sequence changes affect carbohydrate specificity, the two closely related proto-type chicken galectins CG-14 and CG-16 were selected as models. Binding properties were analyzed using a highly sensitive solid-phase assay. We tested 56 free saccharides and 34 well-defined glycoproteins. The two galectins share preference for the II (Galbeta1-4GlcNAc) versus I (Galbeta1-3GlcNAc) version of beta-galactosides. A pronounced difference is found owing to the reactivity of CG-14 with histo-blood group ABH active oligosaccharides and A/B active glycoproteins. These experimental results prompted to determine activity-structure correlations by modeling. Computational analysis included consideration of the flexibility of binding partners and the presence of water molecules. It provided a comparative description of complete carbohydrate recognition domains, which had so far not been characterized in animal galectins. The structural models assigned II, I selectivity to a region downstream of the central Trp moiety. Docking revealed that the tetrasaccharides can be accommodated in their free-state low-energy conformations. CG-14's preference for A versus B epitopes could be attributed to a contact between His124 and the N-acetyl group of GalNAc. Regarding intergalectin comparison, the Ala53/Cys51 exchange affects the interaction potential of His54/His52. Close inspection of simulated dynamic interplay revealed reorientation of His124 at the site of the His124/Glu123 substitution, with potential impact on ligand dissociation. In summary, this study identifies activity differences and provides information on their relation to structural divergence, epitomizing the value of this combined approach beyond galectins.

Entamoeba histolytica: inhibition of cellular functions by overexpression of EhGEF1, a novel Rho/Rac guanine nucleotide exchange factor.

The molecular, biochemical, and cellular characterization of EhGEF1 protein is described. Complete cDNA sequence of 1890 bp revealed an open reading frame that encodes a protein of 69 kDa. EhGEF1 is constituted of Dbl homology domain, pleckstrin homology domain, and several putative regulation sites. Studies of guanine nucleotide exchange activity of EhGEF1 on several GTPases from Entamoeba histolytica and Homo sapiens showed preferential activation on EhRacG, suggesting that EhGEF1 protein could be involved in mechanisms related to actin cytoskeleton activation, cytokinesis, capping, and uroid formation in trophozoite. Confocal microscopy studies of pExEhNeo/HSV-tagged-EhGEF1-transfected cells showed that trophozoites stimulated with ConA, EhGEF1, and EhRacG were localized at plasma membrane. Cellular studies showed that F-actin content of pExEhNeo/HSV-tagged-EhGEF1-transfected trophozoites as well as cellular migration and cell damage capacity were significantly altered. The observations suggest that EhRacG was the principal target of EhGEF1 and that EhGEF1 may provide a link between F-actin dynamics and EhRacG signaling.