Santiago Di Lella

When Galectins Recognize Glycans: From Biochemistry to Physiology and Back Again

In the past decade, increasing efforts have been devoted to the study of galectins, a family of evolutionarily conserved glycan-binding proteins with multifunctional properties. Galectins function, either intracellularly or extracellularly, as key biological mediators capable of monitoring changes occurring on the cell surface during fundamental biological processes such as cellular communication, inflammation, development, and differentiation. Their highly conserved structures, exquisite carbohydrate specificity, and ability to modulate a broad spectrum of biological processes have captivated a wide range of scientists from a wide spectrum of disciplines, including biochemistry, biophysics, cell biology, and physiology. However, in spite of enormous efforts to dissect the functions and properties of these glycan-binding proteins, limited information about how structural and biochemical aspects of these proteins can influence biological functions is available. In this review, we aim to integrate structural, biochemical, and functional aspects of this bewildering and ancient family of glycan-binding proteins and discuss their implications in physiologic and pathologic settings.

Modeling heme proteins using atomistic simulations

Heme proteins are found in all living organisms, and perform a wide variety of tasks ranging from electron transport, to the oxidation of organic compounds, to the sensing and transport of small molecules. In this work we review the application of classical and quantum-mechanical atomistic simulation tools to the investigation of several relevant issues in heme proteins chemistry: (i) conformational analysis, ligand migration, and solvation effects studied using classical molecular dynamics simulations; (ii) electronic structure and spin state energetics of the active sites explored using quantum-mechanics (QM) methods; (iii) the interaction of heme proteins with small ligands studied through hybrid quantum mechanics-molecular mechanics (QM-MM) techniques; (iv) and finally chemical reactivity and catalysis tackled by a combination of quantum and classical tools.

Carbohydrate-binding proteins: Dissecting ligand structures through solvent environment occupancy.

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 proteins 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.

Protein dynamics and ligand migration interplay as studied by computer simulation.

Since proteins are dynamic systems in living organisms, the employment of methodologies contemplating this crucial characteristic results fundamental to allow revealing several aspects of their function. In this work, we present results obtained using classical mechanical atomistic simulation tools applied to understand the connection between protein dynamics and ligand migration. Firstly, we will present a review of the different sampling schemes used in the last years to obtain both ligand migration pathways and the thermodynamic information associated with the process. Secondly, we will focus on representative examples in which the schemes previously presented are employed, concerning the following: i) ligand migration, tunnels, and cavities in myoglobin and neuroglobin; ii) ligand migration in truncated hemoglobin members; iii) NO escape and conformational changes in nitrophorins; iv) ligand selectivity in catalase and hydrogenase; and v) larger ligand migration: the P450 and haloalkane dehalogenase cases. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Solvent structure improves docking prediction in Lectin- Carbohydrate complexes

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.

Critical Role of the Solvent Environment in Galectin-1 Binding to the Disaccharide Lactose †

Galectin-1 (Gal-1), a member of a family of evolutionarily conserved glycan-binding proteins, binds specifically to poly-N-acetyllactosamine-enriched glycoconjugates. Through interactions with these glycoconjugates, this protein modulates inflammatory responses and contributes to tumor progression and immune cell homeostasis. The carbohydrate recognition domain includes the single protein tryptophan (Trp68). UV resonance Raman spectroscopy and molecular dynamic simulation were used to examine the change in the environment of the Trp on ligand binding. The UV Raman spectra and the calculated water radial distribution functions show that, while no large structural changes in the protein follow lactose binding, substantial solvent reorganization occurs. These new insights into the microscopic role of water molecules in Gal-1 binding to its specific carbohydrate ligands provides a better understanding of the physicochemical properties of Gal-1-saccharide interactions, which will be useful for the design of synthetic inhibitors for therapeutic purposes.

An Integrated Computational Analysis of the Structure, Dynamics, and Ligand Binding Interactions of the Human Galectin Network

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.

Structural basis for ligand recognition in a mushroom lectin: solvent structure as specificity predictor

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.

Linking the Structure and Thermal Stability of β-Galactoside-Binding Protein Galectin-1 to Ligand Binding and Dimerization Equilibria

The stability of proteins involves a critical balance of interactions of different orders of magnitude. In this work, we present experimental evidence of an increased thermal stability of galectin-1, a multifunctional beta-galactoside-binding protein, upon binding to the disaccharide lactose. Analysis of structural changes occurring upon binding of lectin to its specific glycans and thermal denaturation of the protein and the complex were analyzed by circular dichroism. On the other hand, we studied dimerization as another factor that may induce structural and thermal stability changes. The results were then complemented with molecular dynamics simulations followed by a detailed computation of thermodynamic properties, including the internal energy, solvation free energy, and conformational entropy. In addition, an energetic profile of the binding and dimerization processes is also presented. Whereas binding and cross-linking of lactose do not alter galectin-1 structure, this interaction leads to substantial changes in the flexibility and internal energy of the protein which confers increased thermal stability to this endogenous lectin. Given that an improved understanding of the physicochemical properties of galectin-glycan lattices may contribute to the dissection of their biological functions and prediction of their therapeutic applications, our study suggests that galectin binding to specific disaccharide ligands may increase the thermal stability of this glycan-binding protein, an effect that could influence its critical biological functions.