[en] This thesis has conducted a comprehensive study on the sequence-structure-function relations for glycosyltransferase proteins with GTA fold, through a bioinformatic computational biology approach and. The thesis is divided into four main chapters. In the first two, is approached this study by an essential enzyme in Mycoplasma genitalium involved in synthesis membrane glycolipids (MG517). In the third chapter, conformational changes of a catalytic loop are studied in the mechanism of a protein of Mycobacterium tuberculosis (GpgS), which starts the 6-O-methyl glucose biosynthetic pathway of lipopolysaccharide (MGLPs) in this organism. Finally, in chapter 4, the relationship between a specific region of the glycosyltransferase with GTA folding and substrate specificity is addressed. Starting with the study of specific proteins, such as Mycoplasma genitalium MG517 protein or Mycobacterium tuberculosis GpgS, through the structures of all proteins crystallized with GTA folding, and the use of all existing sequences for this protein superfamily (over 100000), it has found common characteristics to them all linking sequence, structure and function of each protein. All this, has been achieved using bioinformatics and computational techniques and methods, among which are: homology modeling, simulations of Molecular Dynamics and Metadinámica, the Docking of ligands, overlapping three-dimensional structures, multiple alignments and phylogenetic tree building. Thanks to this study, it has been possible to identify a consensus topology common to all GTAs, with which it has built a three-dimensional model of the N-terminus of the protein MG517, region which hitherto was known structure. The three-dimensional model has been validated experimentally by molecular and dynamic, which has identified key catalytic MG517 positions as general base E193, D40, Y126, Y169, Y218 I170 and substrate binding. Furthermore, for it has been proposed a MG517 monotopic membrane interaction model by an amphipathic helix in the C-terminal region. The same topology consensus has permitted the refinement of a multiple alignment of GTAs, with which we generate a profile Hidden Makov Model (HMM) for the N-terminal region of this group of proteins. This profile facilitates the alignment of GTAs new proteins and identifying their secondary structures. It has also enabled us to identify, within the consensus topology, a region of highly variable sequence and structure, even for proteins of the same family, where each protein specific acceptor is positioned that we have called "Variable Region". Have been described conformational changes in the catalytic loop GpgS protein, by long duration Molecular Dynamics simulations and different calculations of metadinámica using many collective variables. It has been shown that different conformations, are an intrinsic property of the protein, but are displaced towards its inactive form in the absence of ligands. The presence of the donor substrate plus metal in the active site, promotes the movement of the side chains of two residues of the loop, Arg256 and His258, which shifts the equilibrium to the active form and stabilizes by the interaction of His258 with metal, proposing an induced fit mechanism for this protein. Completing these results with docking of ligands calculations, it has been possible to propose the order of entry of the ligands, being the acceptor the first to reach the active site, followed by the donor, when that happens the conformation loop changes. Following these simulations, it has been observed that the interaction of the metal with a histidine residue in the catalytic loop, is a feature common to the vast majority of GTAs families, with the already known motif DXD or tetrad aspartates proposal in literature D-DXD-D, proposed switch to D-DXD-DH. We have studied the evolution of GTAs superfamily and their relationship with the substrates, finding that the overall folding of GTA domain, defines the acceptor specificity and the homology between sequences is more influenced by the acceptor molecule of sugar than the molecule donor. The variable region seems to suffer less evolutionary pressure than the rest of the sequence, which explains its greater variability, however, contain highly conserved residues that interact with the specific acceptor for each protein. This region has been used for profiling specific HMM for each acceptor and protein family. These profiles have been successfully used for screening GTA proteins of unknown function and predicting their acceptor.

[es] En esta tesis se ha realizado un estudio global, acerca de las relaciones secuencia-estructura-función de proteínas glicosiltransferasa con plegamiento GTA, mediante un enfoque bioinformático y de biología computacional. La tesis está estructurada en cuatro capítulos principales. En los dos primeros se aborda este estudio para una enzima esencial de Mycoplasma genitalium involucrada en la síntesis de glicolípidos de membrana (MG517). En el tercer capítulo se estudian los cambios conformacionales de un bucle catalítico, en el mecanismo de una proteína de Micobacterium tuberculosis (GpgS), que inicia la ruta biosintética de los lipopolisacáridos 6-O-metilglucosa (MGLPs) en este organismo. Por último, en el capítulo 4 se aborda la relación entre una región específica de las glicosiltransferasas con plegamiento GTA, y la especificidad por el sustrato. Comenzando por el estudio de proteínas específicas, como la proteína MG517 de Mycoplasma genitalium o GpgS de Micobacterium tuberculosis, pasando por el de las estructuras de todas las proteínas cristalizadas con plegamiento GTA, hasta la utilización de todas las secuencias existentes para esta superfamilia de proteínas (más de 100000), se han encontrado características comunes a todas ellas que relacionan la secuencia, con la estructura y la función de cada proteína. Todo ello se ha conseguido utilizando técnicas y métodos bioinformáticos y computacionales, entre los que destacan: el modelado por homología, las simulaciones de Dinámica Molecular y Metadinámica, el Docking de ligandos, la superposición de estructuras tridimensionales, alineamientos múltiples y la construcción de árboles filogenéticos. Gracias a este estudio, se ha podido identificar una topología consenso común a todas las GTAs, con la que se ha construido un modelo tridimensional de la región N-terminal de la proteína MG517, de la que hasta ahora no existía estructura conocida. El modelo tridimensional ha sido validado por dinámica molecular y experimentalmente, lo cual ha permitido identificar posiciones catalíticas clave en MG517 como E193 base general, D40, Y126, Y169, I170 y Y218 de unión a sustratos. Además, para MG517 se ha propuesto un modelo de interacción monotópica con la membrana, mediante una hélice anfipática en su región C-terminal. La misma topología consenso, ha permitido el refinamiento de un alineamiento múltiple de GTAs, con el que se ha generado un perfil Hidden Makov Model (HMM) para la región N-terminal de este grupo de proteínas. Este perfil facilita el alineamiento de nuevas proteínas GTAs y la identificación de sus estructuras secundarias. También ha permitido identificar, dentro de la topología consenso, una región de secuencia y estructura muy variable, incluso para proteínas de la misma familia, donde se posiciona el aceptor específico de cada proteína y que hemos denominado "Región Variable". Se han descrito los cambios conformacionales del bucle catalítico en la proteína GpgS mediante simulaciones de dinámica molecular de larga duración y distintos cálculos de metadinámica utilizando multitud de variables colectivas. Se ha demostrado que las distintas conformaciones son una propiedad intrínseca de la proteína, pero están desplazadas hacia su forma inactiva en ausencia de ligandos. La presencia del sustrato dador más el metal en el centro activo, promueve el movimiento de las cadenas laterales de dos residuos del bucle, Arg256 e His258, que desplaza el equilibrio hacia la forma activa y la estabiliza mediante la interacción de His258 con el metal, proponiéndose un mecanismo de ajuste inducido para esta proteína. Completando estos resultados con cálculos de docking de ligandos, se ha podido proponer el orden de entrada de los ligandos, siendo el aceptor el primero en llegar al centro activo, seguido por el dador, momento en que sucede el cambio conformacional en el bucle. A raíz de estas simulaciones, se ha observado que la interacción del metal con un residuo de histidina en el bucle catalítico, es una característica común a la gran mayoría de familias GTAs, junto a las ya conocidas del motivo DXD o la tétrada de aspartatos propuesta en literatura D-DXD-D, que se propone cambiar a D-DXD-D-H. Se ha estudiado la evolución de la superfamilia GTAs y su relación con los sustratos, encontrando que el plegamiento global del dominio GTA, define la especificidad del aceptor y que la homología entre secuencias está más influida por esta molécula aceptora del azúcar que por la molécula dadora. La región variable parece sufrir una presión evolutiva menor que el resto de la secuencia, lo que explica su mayor variabilidad, sin embargo, contiene residuos altamente conservados que interaccionan con el aceptor específico de cada proteína. Se ha utilizado esta región para la generación de perfiles HMM específicos para cada aceptor y familia de proteínas. Estos perfiles se han utilizado con éxito para el screening de proteínas GTA de función desconocida y la predicción de su aceptor.