RESUMO
DCL1 is the ribonuclease that carries out miRNA biogenesis in plants. Substrate pri-miRNA recognition by DCL1 requires two double stranded RNA binding domains located at the C-terminus of the protein. We have previously shown that the first of these domains, DCL1-A, is intrinsically disordered and folds upon binding pri-miRNA. Integrating NMR and SAXS data, we study here the conformational landscape of free DCL1-A through an ensemble description. Our results reveal that secondary structure elements, corresponding to the folded form of the protein, are transiently populated in the unbound state. The conformation of one of the dsRNA binding regions in the free protein shows that, at a local level, RNA recognition proceeds through a conformational selection mechanism. We further explored the stability of the preformed structural elements via temperature and urea destabilization. The C-terminal helix is halfway on the folding pathway in free DCL1-A, constituting a potential nucleation site for the final folding of the protein. In contrast, the N-terminal helix adopts stable non-native structures that could hinder the correct folding of the protein in the absence of RNA. This description of the unfolded form allows us to understand details of the mechanism of binding-induced folding of the protein.
Assuntos
Proteínas de Arabidopsis/metabolismo , Proteínas de Ciclo Celular/metabolismo , Proteínas Intrinsicamente Desordenadas/química , MicroRNAs/metabolismo , Ribonuclease III/metabolismo , Arabidopsis , Proteínas de Arabidopsis/química , Proteínas de Ciclo Celular/química , Dicroísmo Circular , Espectroscopia de Ressonância Magnética , MicroRNAs/química , Modelos Químicos , Ligação Proteica , Conformação Proteica , Domínios Proteicos , Dobramento de Proteína , Ribonuclease III/química , Espalhamento a Baixo Ângulo , Temperatura , Difração de Raios XRESUMO
Recognition and complex formation between proteins and carbohydrates is a key issue in many important biological processes. Determination of the three-dimensional structure of such complexes is thus most relevant, but particularly challenging because of their usually low binding affinity. In silico docking methods have a long-standing tradition in predicting protein-ligand complexes, and allow a potentially fast exploration of a number of possible protein-carbohydrate complex structures. However, determining which of these predicted complexes represents the correct structure is not always straightforward. In this work, we present a modification of the scoring function provided by AutoDock4, a widely used docking software, on the basis of analysis of the solvent structure adjacent to the protein surface, as derived from molecular dynamics simulations, that allows the definition and characterization of regions with higher water occupancy than the bulk solvent, called water sites. They mimic the interaction held between the carbohydrate -OH groups and the protein. We used this information for an improved docking method in relation to its capacity to correctly predict the protein-carbohydrate complexes for a number of tested proteins, whose ligands range in size from mono- to tetrasaccharide. Our results show that the presented method significantly improves the docking predictions. The resulting solvent-structure-biased docking protocol, therefore, appears as a powerful tool for the design and optimization of development of glycomimetic drugs, while providing new insights into protein-carbohydrate interactions. Moreover, the achieved improvement also underscores the relevance of the solvent structure to the protein carbohydrate recognition process.
Assuntos
Carboidratos/química , Simulação de Dinâmica Molecular , Proteínas/química , Solventes/química , Sítios de Ligação , Galectinas/química , Ligantes , Simulação de Acoplamento Molecular , Ligação Proteica , Conformação Proteica , Estrutura Terciária de Proteína , Software , Água/química , Água/metabolismoRESUMO
Galectins, a family of evolutionarily conserved animal lectins, have been shown to modulate signaling processes leading to inflammation, apoptosis, immunoregulation, and angiogenesis through their ability to interact with poly-N-acetyllactosamine-enriched glycoconjugates. To date 16 human galectin carbohydrate recognition domains have been established by sequence analysis and found to be expressed in several tissues. Given the divergent functions of these lectins, it is of vital importance to understand common and differential features in order to search for specific inhibitors of individual members of the human galectin family. In this work we performed an integrated computational analysis of all individual members of the human galectin family. In the first place, we have built homology-based models for galectin-4 and -12 N-terminus, placental protein 13 (PP13) and PP13-like protein for which no experimental structural information is available. We have then performed classical molecular dynamics simulations of the whole 15 members family in free and ligand-bound states to analyze protein and protein-ligand interaction dynamics. Our results show that all galectins adopt the same fold, and the carbohydrate recognition domains are very similar with structural differences located in specific loops. These differences are reflected in the dynamics characteristics, where mobility differences translate into entropy values which significantly influence their ligand affinity. Thus, ligand selectivity appears to be modulated by subtle differences in the monosaccharide binding sites. Taken together, our results may contribute to the understanding, at a molecular level, of the structural and dynamical determinants that distinguish individual human galectins.
Assuntos
Galectina 4/análise , Galectinas/análise , Polissacarídeos/metabolismo , Proteínas da Gravidez/análise , Transdução de Sinais/fisiologia , Biologia de Sistemas/métodos , Sequência de Aminoácidos , Sítios de Ligação , Cristalografia por Raios X , Bases de Dados de Proteínas , Entropia , Epitopos , Galectina 4/química , Galectina 4/imunologia , Galectina 4/metabolismo , Galectinas/química , Galectinas/imunologia , Galectinas/metabolismo , Humanos , Ligantes , Espectroscopia de Ressonância Magnética , Modelos Moleculares , Simulação de Dinâmica Molecular , Dados de Sequência Molecular , Filogenia , Polissacarídeos/imunologia , Proteínas da Gravidez/química , Proteínas da Gravidez/imunologia , Proteínas da Gravidez/metabolismo , Ligação Proteica , Estrutura Terciária de Proteína , Homologia de Sequência de AminoácidosRESUMO
Lectins are able to recognize specific carbohydrate structures through their carbohydrate recognition domain (CRD). The lectin from the mushroom Agaricus bisporus (ABL) has the remarkable ability of selectively recognizing the TF-antigen, composed of Galß1-3GalNAc, Ser/Thr linked to proteins, specifically exposed in neoplastic tissues. Strikingly, the recently solved crystal structure of tetrameric ABL in the presence of TF-antigen and other carbohydrates showed that each monomer has two CRDs, each being able to bind specifically to different monosaccharides that differ only in the configuration of a single hydroxyl, like N-acetyl-d-galactosamine (GalNAc) and N-acetyl-d-glucosamine (GlcNAc). Understanding how lectin CRDs bind and discriminate mono and/or (poly)-saccharides is an important issue in glycobiology, with potential impact in the design of better and selective lectin inhibitors with potential therapeutic properties. In this work, and based on the unusual monosaccharide epimeric specificity of the ABL CRDs, we have performed molecular dynamics simulations of the natural (crystallographic) and inverted (changing GalNAc for GlcNAc and vice-versa) ABL-monosaccharide complexes in order to understand the selective ligand recognition properties of each CRD. We also performed a detailed analysis of the CRD local solvent structure, using previously developed methodology, and related it with the recognition mechanism. Our results provide a detailed picture of each ABL CRD specificity, allowing a better understanding of the carbohydrate selective recognition process in this particular lectin.
Assuntos
Lectinas/química , Acetilgalactosamina/química , Acetilglucosamina/química , Agaricus/química , Agaricus/metabolismo , Antígenos Glicosídicos Associados a Tumores/metabolismo , Sítios de Ligação , Configuração de Carboidratos , Ligação de Hidrogênio , Lectinas/metabolismo , Simulação de Dinâmica Molecular , Ligação Proteica , TermodinâmicaRESUMO
Formation of protein ligand complexes is a fundamental phenomenon in biochemistry. During the process, significant solvent reorganization is produced along the contact surface and many water molecules strongly bound to the protein's ligand binding site must be displaced. Both the thermodynamics and kinetics of this process are complex and a clear understanding at the microscopic level has been not achieved so far. Special attention has been paid to the structure of water molecules on carbohydrate recognition sites of various proteins, and many studies support the idea that displacement of these water molecules should have a crucial effect on the binding free energy. Molecular dynamics (MD) simulations in explicit water solvent is a very promising approach for this type of studies. Using MD simulations combined with statistical mechanics analysis, thermodynamic properties of these water molecules can be computed and analyzed in a comparative view. Using this idea, we developed a set of analysis tools to link solvation with ligand binding in a key carbohydrate binding protein, human galectin-1 (hGal-1). Specifically, we defined water sites (WS) in terms of the thermodynamic properties of water molecules strongly bound to protein surfaces. In the present work, we selected a group of proteins whose ligand bound complexes have been already structurally characterized in order to extend the analysis of the role of the surface associated water molecules in the ligand binding and recognition process. The selected proteins are concanavalin-A (Con-A), galectin-3 (Gal-3), cyclophilin-A (Cyp-A), and two modules CBM40 and CBM32 of the multimodular bacterial sialidase. Our results show that the probability of finding water molecules inside the WS, p(v), with respect to the bulk density is directly correlated to the likeliness of finding an hydroxyl group of the ligand in the protein-ligand complex. This information can be used to analyze in detail the solvation structure of the carbohydrate recognition domain (CRD) and its relation to the possible protein ligand complexes and suggests addition of OH-containing functional groups to displace water from high p(v) WS to enhance drugs, specially glycomimetic-drugs, protein affinity, and/or specificity.