Accurate prediction of binding free energies of protein-ligand system has long been a focus area for theoretical and computational studies; with important implications in fields like pharmaceuticals, enzyme-redesign, etc. The aim of this project was to develop such a predictive model for calculating binding free energies of protein-ligand systems based on the LIE-SASA methods. Many models have been successfully fit to experimental data, but a general predictive model, not reliant on experimental values, would make LIE-SASA a more powerful and widely applicable method. The model was developed such that There is no significant increase in computational time No increase in complexity of system setup No increase in the number of empirical parameters. The method was tested on a small number of protein-ligand systems, selected with certain constraints. This was our training set, from which we obtain the complete expression for binding free energy. Expectedly, there was good agreement with experimental values for the training set On applying our model to a similar sized validation set, with the same selection constraints as for the training set, we achieved even better agreement with experimental results, with lower standard errors. Finally, the model was tested by applying it to a set of systems without such selection constraints, and again found good agreement with experimental values. In terms of accuracy, the model was comparable to a system specific empirical fit that was performed on this set. These encouraging results could be an indicator of generality.

<p>Thermodynamic stability measurements on proteins and protein-ligand complexes can offer insights not only into the fundamental properties of protein folding reactions and protein functions, but also into the development of protein-directed therapeutic agents to combat disease. Conventional calorimetric or spectroscopic approaches for measuring protein stability typically require large amounts of purified protein. This requirement has precluded their use in proteomic applications. Stability of Proteins from Rates of Oxidation (SPROX) is a recently developed mass spectrometry-based approach for proteome-wide thermodynamic stability analysis. Since the proteomic coverage of SPROX is fundamentally limited by the detection of methionine-containing peptides, the use of tryptophan-containing peptides was investigated in this dissertation. A new SPROX-like protocol was developed that measured protein folding free energies using the denaturant dependence of the rate at which globally protected tryptophan and methionine residues are modified with dimethyl (2-hydroxyl-5-nitrobenzyl) sulfonium bromide and hydrogen peroxide, respectively. This so-called Hybrid protocol was applied to proteins in yeast and MCF-7 cell lysates and achieved a ~50% increase in proteomic coverage compared to probing only methionine-containing peptides. Subsequently, the Hybrid protocol was successfully utilized to identify and quantify both known and novel protein-ligand interactions in cell lysates. The ligands under study included the well-known Hsp90 inhibitor geldanamycin and the less well-understood omeprazole sulfide that inhibits liver-stage malaria. In addition to protein-small molecule interactions, protein-protein interactions involving Puf6 were investigated using the SPROX technique in comparative thermodynamic analyses performed on wild-type and Puf6-deletion yeast strains. A total of 39 proteins were detected as Puf6 targets and 36 of these targets were previously unknown to interact with Puf6. Finally, to facilitate the SPROX/Hybrid data analysis process and minimize human errors, a Bayesian algorithm was developed for transition midpoint assignment. In summary, the work in this dissertation expanded the scope of SPROX and evaluated the use of SPROX/Hybrid protocols for characterizing protein-ligand interactions in complex biological mixtures.</p> / Dissertation

A central goal of computational biophysics and biochemistry is to understand the behavior, interactions, and reactions of molecules, and to interpret and facilitate experimental design. The objective of this thesis research is to use the molecular modeling and simulation techniques to advance our understanding of principles in molecular structure properties, recognition and interaction at the atomic level. First, a physical molecular mechanics model is built to study the conformational properties of depsipeptide, which shows potential for engineered protein mimetics with controllable structure and function. We explore the possible kinase-substrate binding modes and the likelihood of an [alpha]-helix docking interaction within a kinase active site. Finally, efficient physical models based on a polarizable potential function are developed to describe the structural properties and calculate protein-ligand binding affinities accurately for both trypsin and matrix metalloproteinase. / text

<p>While a number of different proteomic, genomic, and computational approaches exist for the characterization of drug action, each of the experimental approaches developed to date has both strengths and weaknesses. Currently, there is no one "perfect" assay for drug mode-of-action studies. A protocol that could assay all the proteins in the proteome for both direct and indirect binding interactions of drugs would greatly facilitate studies of drug action. Recently, the SPROX (stability of proteins from rates of oxidation) technique was developed as a chemical modification- and mass spectrometry-based strategy for detecting protein-ligand interactions by monitoring the change in thermodynamic stability of proteins upon ligand binding. This is accomplished by monitoring the denaturant dependent oxidation of globally protected methionine residues. The SPROX technique has been interfaced with bottom-up proteomics methods to allow for the proteome-wide analysis of protein-ligand interactions. However, the strategy has been limited by the need to detect and quantify methionine containing peptides in the bottom-up proteomics experiment. </p><p>The work in this dissertation is focused on evaluating the current SPROX protocol, developing modifications to improve proteome coverage, and applying the SPROX platform to two different drug mode-of-action studies. Three main strategies were employed to improve protein coverage. First, a chemo-selective isolation of un-oxidized methionine containing peptides was employed to enrich for methionine containing peptides, and it was found to produce a ~2-fold improvement in proteomic coverage. Second, a pre-fractionation strategy involving the use of isoelectric focusing was employed to decrease sample complexity prior to LC-MS/MS analysis and it was found to generate a ~2-3 fold improvement in proteomic coverage, however when combined with the methionine enrichment strategy the improvement was ~6-fold as the benefits of both were additive. Third, a tryptophan modification strategy was developed that could ultimately expand the number of useful peptides in proteome-wide SPROX experiments to include those that contain tryptophan. Also, investigated was the use of several different mass spectrometer systems (including a bench-top quadrupole and orbitrap system and two different quadrupole time-of-flight systems) in the SPROX protocol. The results of these studies indicate that there is a significant advantage in proteome coverage when faster mass spectrometers are used. The use of high energy collision dissociation (HCD) in the orbitrap system was also more advantageous than the use of collision induced dissociation (CID) in the Q-ToF systems. Regardless of the mass spectrometer used, the major source of error in the SPROX experiment was found to be the random error associated with the LC-MS/MS analysis of isobaric mass tagged peptides. This random error was found to yield a false discovery rate of between 3 and 10% for "hit" peptides in the SPROX experiment. </p><p>The above improvements in the SPROX protocol were used in two protein-ligand binding experiments. One set of experiments involved studies on two small molecules with a specific anti-cancer phenotype in human colon cancer cells. These studies identified 17 proteins as potential "hits" of these two small molecules. After preliminary validation of these proteins, approximately 50% were eliminated as false positives and one protein, p80/nucleophosim, showed consistent data indicating a destabilizing interaction with both small molecules. The destabilization is indicative of an indirect interaction with the small molecules that would be mediated through a protein-protein interaction network. In another set of experiments the breast cancer drug, tamoxifen, and its main, active metabolite, 4-hydroxy tamoxifen, were assayed for binding to the proteins in a yeast cell lysate to better understand its adverse effects on yeast cells. The results of these studies identified ~80 proteins as potential "hits" of these two drugs. After preliminary validation of these proteins, approximately 30% were eliminated as false positives and one protein, SIS1, type II Hsp40, showed consistent data indicative of a direct binding interaction.</p> / Dissertation

Virtually all processes in living organisms are conducted by proteins. Proteins perform their function by binding to other proteins (protein-protein interactions) or small molecules - so called ligands (protein-ligand interactions). Active sites for protein-ligand interactions are pockets in protein structure where ligand can bind. Predicting of ligand binding sites is the first step to study and predict protein functions and structure based drug-design. In this thesis we reviewed current approaches for binding site prediction and proposed our own improvement. We have developed a novel pocket ranking function based on prediction model that predicts ligandability (ability to bind a ligand) of a given point inside of a pocket. Prediction is done considering only a local physicochemical and geometric properties derived from neighbourhood.

Conformationally constrained ligands and their flexible analogues were prepared as inhibitors of the Grb2 SH2 domain in order to study the structural and energetic effects of ligand preorganization in protein-ligand interactions. The compounds were prepared by using trans-cyclopropane-containing amino acid mimics, macrocyclization, or [alpha,alpha]-disubstituted amino acid residues. All trans-cyclopropane containing peptides were more potent than their corresponding succinate containing analogues due to an enthalpic advantage. Surprisingly, the binding of constrained peptides to the domain was entropically disfavored relative to their flexible controls. Effects of proton transfer and desolvation as being the source of the unprecedented entropic penalty for the constrained ligands relative to their respective controls were precluded, and X-ray crystallographic studies revealed that the binding conformations for the respective cyclopropane and succinate containing ligands were similar. This led us to believe that differential changes in protein dynamics may occur upon binding of the constrained and flexible ligands, which could contribute to the observed binding energetics. Two 23-membered macrocyclic ligands were slightly more potent than their corresponding linear controls. The amino acids used to link the N- and C-termini of the linear peptides to form the macrocycles were found to affect the energetics of binding. In one case, the 23-membered macrocycle was more potent than its control due to an entropic advantage, whereas the other 23-membered macrocycle was more potent than its control because it benefited from an enthalpic advantage. [alpha,alpha]-Disubstituted and [alpha]-monosubstituted residues that varied in hydrophobic character were incorporated into Grb2 SH2 domain binding tripeptides, and binding became more favorable as nonpolar surface area increased only for the set of tripeptides possessing cyclic [alpha,alpha]-disubstituted residues. The increase in affinity was due to an increasing enthalplic term, whereas the entropy of binding became less favorable. A total synthesis of (±)-lycopladine A was achieved in five steps from known compounds. The tricyclic core of the natural product was prepared utilizing a novel two-step sequence comprising a conjugate addition of a metalated picoline derivative followed by an intramolecular enolate arylation. It was demonstrated that the natural product existed in a solvent dependent equilibrium with its isomeric lactol. / text

El trabajo presentado en esta tesis cubre varios campos de investigación relacionados con el desarrollo de moléculas bioactivas. Se compone de cinco partes distintas que se resumen aquí. Predicción de la utilidad farmacológica de dianas terapéuticas. El desarrollo de fármacos está generalmente dirigido a inhibir la función de una proteína específica. Pero para validar esta proteína como diana terapéutica, al principio de un proyecto de descubrimiento de fármacos se tiene que saber si una molécula de tipo fármaco puede unirse con suficiente afinidad a la proteína como para alterar su función. Existen métodos que predicen si una potencial diana terapéutica es tratable o no por vía farmacológica, lo que se ha dado en llamar ‘druggability’. El problema es que estos métodos no están accesibles libremente y su validación es discutible. En la primera parte de la tesis se ha compilado un conjunto extensivo de datos de cavidades en proteínas cuyo novel de ‘druggability’ es conocido, haciéndolo accesible en una plataforma web pública (http://fpocket.sourceforge.net/dcd). Estos datos pueden ser modificados por cualquier persona que quiera contribuir al desarrollo de este conjunto de datos, aumentando su volumen o mejorando su calidad. En estudios previos, los sitios druggable se han asociado a cavidades profundas e hidrofóbicas, ignorando la importancia que tiene los grupos polares en el sitio de unión y su posible relación con la ‘druggability’. Utilizando el set de datos compilado previamente, hemos encontrado que aunque las cavidades ‘druggables’ son mas hidrofóbicas, también tienen grupos polares más expuestos pero con poca superficie de interacción. Esta observación es objeto de posteriores investigaciones en la segunda parte de la tesis. Finalmente, se ha utilizado un algoritmo de búsqueda de sitios de unión, fpocket, que empecé a desarrollar como proyecto de master. Este programa se ha utilizado para extraer todas las características de las cavidades ‘druggables’ y no ‘druggables’ y estos parámetros se han utilizado para entrenar un modelo logístico capaz de predecir si un sitio es druggable o no. Demostramos que el algoritmo y la función de puntuación desarrollado durante esta primera parte predice la ‘druggability’ de manera fiable. Los resultados son de igual calidad a los obtenidos con el único otro programa accesible con funcionalidad parecida (SiteFinder, de Schrödinger), pero nuestro programa tiene las siguientes importantes ventajas: 1) es libre; 2) es mucho más eficiente computacionalmente; y 3) trabaja sobre cavidades detectadas automáticamente por el programa, lo que permite aplicaciones a gran escala. En otras partes de la tesis se verá como su aplicación al conjunto del PDB permite novedosas aplicaciones en el área del diseño de fármacos. Análisis de movilidad de cavidades de las proteínas Existen una gran variedad de algoritmos que permiten identificar posibles sitios de unión en las estructuras tridimensionales de las proteínas. El trabajo presentado en la sección anterior de esta tesis permitía extender uno de estos algoritmos para caracterizar la ‘druggability’ de las cavidades. Un problema de gran calado en el estado actual de la técnica, tanto de detección de sitios de unión como en diseño de fármacos en general,4 es que las proteínas se tratan como un cuerpo rígido a pesar de que en realidad gozan de una gran movilidad estructural. El objetivo de esta sección era otorgar todavía otra funcionalidad a fpocket, el programa de predicción de sitios de unión, para permitir también la detección y el análisis de cavidades de proteínas en movimiento. Habitualmente, los movimientos de las proteína se pueden simular usando la dinámica molecular (MD). Una herramienta capaz de analizar conjuntos de estructuras derivados de MD u otras fuentes puede ser, por tanto, extremadamente útil para observar la aparición de cavidades transitorias y su plasticidad. Como resultado final de este trabajo, se presenta un nuevo programa informático, llamado MDpocket y que se enmarca dentro del paquete fpocket. Para cada conformación de la proteína, se ejecuta un ciclo de detección de cavidades con fpocket. Los resultados de este proceso se plasman sobre una malla tridimensional superpuesta a la estructura de la proteína. La malla puede entonces ser visualizada o analizada en mayor detalle. Lo primero se puede llevar a cabo con programas de visualización molecular tales como PyMOL, VMD o Chimera. Otra funcionalidad dentro de MDpocket es la de seguir la evolución de las propiedades de una cavidad o zona de interés (definida por el usuario) a lo largo del tiempo. Cabe destacar que MDpocket es igualmente capaz de identificar sitios de unión de moléculas tipo fármaco como pequeños canales en la matriz proteica que pueden ser importantes para la migración de pequeños ligandos como gases o moléculas de agua. Las posibles aplicaciones de MDpocket se ejemplifican en tres casos distintos. El primero es la capacidad de identificar la apertura transitoria de cavidades en el sitio de unión de ATP en la proteína HSP90. En el segundo ejemplo, se muestra como MDpocket permite identificar un canal de migración de moléculas biatómicas en mioglobina, un sistema de referencia bien conocido. Aquí se demuestra que MDpocket puede, no tan solo identificar los sitios internos de unión a Xenón, sino también los canales que se abren de forma transitoria para permitir a los ligandos migrar de un sitio a otro. En el último ejemplo, las propiedades del sitio de unión a ATP de la proteína cinasa P38 se analizaron a lo largo de una trayectoria de MD, evaluando la capacidad de MDpocket para identificar aquellas conformaciones que pueden ser particularmente útiles para realizar docking molecular. Uno de los principales problemas en docking de proteína-ligando es que el receptor generalmente se considera rígido, mientras que si se utilizan múltiples conformaciones (cristalográficas o derivadas de MD) para representar al receptor es difícil decidir a priori cuales de ellas pueden dar mejores resultados. Aquí mostramos que MDpocket se puede utilizar para seleccionar conformaciones concretas de una trayectoria de MD para usarlas en procesos de docking. Concretamente, hemos observado que la densidad hidrofóbica promedio (previamente identificada como un descriptor importante para predecir ‘druggability’) correlaciona bien con la probabilidad de que el modo de unión de ligandos pueda ser predicho correctamente. Tal como en el trabajo anterior, MDpocket se incluye dentro del proyecto fpocket y se está accesible como una herramienta libre y de código abierto. Relaciones estructura-cinética de unión El control de los tiempos en interacciones moleculares es una propiedad esencial de los sistemas bioquímicos, pero poco se conoce sobre los factores estructurales que gobiernan la cinética de los procesos de reconocimiento molecular. Partiendo de una observación realizada durante el trabajo de predicción de ‘druggability’, aquí se ha investigado el papel que átomos con poca superficie expuesta a solvente pueden jugar en los sitios de unión de la proteína. En particular, encontramos que los átomos polares en los sitios druggable son minoritarios en comparación con los átomos apolares, pero si bien pueden tener poca superficie accesible, tienden a ser mas protuberantes, lo que los hace más accesibles para establecer interacciones. Hemos establecido que esta propiedad puede estar relacionada con la cinética de unión/disociación de un ligando a su cavidad en el receptor. En diseño de fármacos, la vida media del complejo formado entre el fármaco y su diana terapéutica determina en gran medida sus efectos biológicos, pero en ausencia de relaciones estructura cinética, se hace imposible optimizar esta propiedad de forma racional. Aquí se muestra que átomos polares prácticamente enterrados (ABPAs) – un elemento comúnmente encontrado en los sitios de unión de proteínas – tienden a formar puentes de hidrógeno que están protegidos de las moléculas de agua. La formación y ruptura de este tipo de puentes de hidrógeno implica un estado de transición penalizado energéticamente porque ocurre de modo asincrónico con el proceso de deshidratación/rehidratación. En consecuencia, los puentes de hidrógeno protegidos se intercambian a velocidades lentas. Estas conclusiones se basan en el estudio computacional del proceso de unión de un pequeño ligando a un sitio de unión modelo. El receptor modelo se construyó para permitir modular tanto el grado de exposición del átomo polar como la curvatura del entorno apolar. Mediante el uso de dinámicas moleculares con constricciones y la relación de Jarzinsky, se obtuvieron los perfiles de energía libre de unión para cada cavidad. La presencia de un estado de transición (y por tanto menor velocidad de asociación/disociación) puede anticiparse mediante un simple análisis estructural tal como la medición de la superficie accesible del átomo polar o su grado de protrusión. Esto constituye una nueva y valiosa clave para interpretar y predecir relaciones estructura-actividad, que se ha puesto a prueba investigando sistemas reales. En primer lugar, analizando tanto estructuras cristalográficas depositadas en el PDB como trayectorias de dinámica molecular, se ha demostrado que aquellas moléculas de agua que forman puentes de hidrógeno con ABPAs tienden a tener menor movilidad e intercambios más lentos. Posteriormente, la validez del principio se ha demostrado en dos pares de inhibidores de la proteína Hsp90, una diana terapéutica para cáncer, para los que se han obtenido datos estructurales, termodinámicos y cinéticos mediante distintas técnicas experimentales. El acuerdo entre observables macroscópicas y los resultados de simulaciones moleculares confirma la función de los puentes de hidrógeno protegidos de solvente como trampas cinéticas e ilustra como nuestro hallazgo puede ser usado para facilitar el proceso de diseño de fármacos basado en estructura. Base de datos de cavidades: hacia el ‘pocketoma’ El trabajo presentado en el principio de la tesis perseguía un objetivo muy especifico, que se enmarca en un proyecto mayor del grupo de investigación. La herramienta de predicción de ‘druggability’ se desarrolló con el fin de cribar grandes bases de datos estructurales, tales como el PDB, identificando complejos proteína-proteína no obligados (transitorios) que contengan en su interfase una cavidad potencialmente capaz de unir moléculas de tipo fármaco. Con ello se pretende estabilizar selectivamente dicha interacción y conseguir un efecto biológico que pueda ser terapéutico. Dado que la información sobre cavidades y su druggability asociada va a ser explotadas por otras personas en el grupo en este y otro tipo de proyectos encaminados a facilitar la explotación de nuevos mecanismos de acción, es necesario crear una base de datos que contenga esta información y sea fácilmente navegable. Para empezar, se ejecutó el programa para cada una de las estructuras depositadas en el PDB, identificando todas las cavidades y extrayendo sus descriptores, entre los que se incluye la función de puntuación de druggability. Esta información se guarda de forma organizada en una base de datos relacional (pocketDB), que se relaciona con otras bases de datos tales como Uniprot y Uniref, que contienen información sobre secuencias. Igualmente, se incluyeron información sobre la estructura cuaternaria y otros recursos tales como Kegg, haciendo de pocketDB un potente recurso para filtrar de modo eficiente millones de cavidades, identificando aquellas que tengan mayor interés para cada proyecto. De modo particular, cabe destacar la aplicación de pocketDB a la identificación de cavidades ‘druggable’ situadas en la interfase de complejos transitorios proteína-proteína, resultando en 39 complejos candidatos, entre los que se recobraron 3 casos conocidos previamente, lo que valida la metodología. Además otros tres sistemas identificados, después de una inspección minuciosa, han sido seleccionados en el grupo como candidatos para realizar la prueba de concepto que valide esta nueva estrategia farmacológica. Uno de ellos está actualmente en fase de validación experimental. El ‘pocketoma’ La última sección de mi tesis presenta un proyecto que hace un uso extensivo de la base de datos de cavidades (pocketDB) previamente presentada. Dos cuestiones importantes aún persisten hoy día en el proceso de descubrimiento de fármacos y, de algún modo, contribuyen a la alta tasa de fracasos en las fases tardías del desarrollo, que normalmente se explican por la baja eficacia o el exceso de efectos secundarios (normalmente tóxicos) para el organismo. Ambas situaciones pueden ser explicadas por un fenómeno común: la falta de selectividad. Las fármacos interaccionan en la célula con una diversa variedad de macromoléculas además de con la diana terapéutica, induciendo así efectos secundarios imprevistos. Asimismo, considerar solamente una diana para nuestro fármaco puede resultar menos efectivo de lo esperado, puesto que la pérdida de función de una sola proteína diana puede ser fácilmente reemplazada debido a los mecanismos homeostáticos controlados por robustas redes de interacciones, derivando en una pérdida de eficacia de nuestra molécula. Este trabajo pretende sentar las bases para poder establecer relaciones entre macromoléculas biológicas en base a su potencial para interaccionar con una misma entidad química (fármaco). El trabajo se basa en la asunción de que cavidades similares son capaces de unir ligandos similares. Existen varios métodos que calculan la similitud entre cavidades de unión, sin embargo, hasta la fecha no se ha realizado un análisis relacional profundo de todas las cavidades en el PDB que permita establecer relaciones entre ellos. Con el fin de cumplir con los requerimientos técnicos de unos objetivos tan ambiciosos, se ha desarrollado un novedoso método de comparación de cavidades. En este particular método se considera una representación muy abstracta de la cavidad. Dicha representación reúne información tanto sobre la forma de la cavidad cómo sobre la distribución por pares de los puntos de interacción en la superficie. No obstante, la información sobre la topología exacta de la cavidad es ignorada. Tal abstracción intenta relacionar cavidades que estructuralmente están alejadas, pero que se parecen entre ellas en términos de forma y propiedades fisicoquímicas globales. Usando varios ejemplos de validación en un conjunto de cavidades de referencia, el método resulta capaz de recuperar de un gran set de cavidades aquellas más similares y relacionadas entre sí. Del mismo modo, el método demuestra ser útil para encontrar relaciones entre cavidades previamente no relacionadas y que sin embargo unen los mismos ligandos o moléculas muy similares. Esta nueva implementación se ha utilizado en un experimento de exploración a gran escala para encontrar (i) las mismas cavidades en diferentes estructuras, (ii) cavidades relacionadas y (iii) cavidades con la misma estructura de ligando. Se obtuvieron excelentes resultados para estas tres categorías. Seguidamente, se compararon entre ellas todas las cavidades encontradas en el PDB. Se desarrolló una herramienta computacional novedosa para permitir la navegación en el 'pocketoma' resultante de las comparaciones, que será posible descargar gratuitamente. Los resultados presentados muestran por primera vez un espacio global interrelacionando cavidad-ligando para todas las cavidades del PDB. Finalmente, se señalan a continuación dos aplicaciones que ponen de manifiesto el posible impacto del ‘pocketoma’ y la herramienta de navegación. En el primer ejemplo, una cavidad no caracterizada encontrada en GSK3-β fue comparada contra todas las cavidades del PDB, permitiendo obtener una variedad de cavidades, y sistemas, supuestamente relacionadas. Entre los resultados se encontraban las subunidades de unión a ADP dhaL y dhaM de la PTS dependiente di-hidroacetona cinasa. Se encontró que el sitio de unión a ADP era muy similar a la cavidad investigada en GSK3-β con una sorprendente similitud estructural entre ambos. En el último ejemplo, se navegó por el ‘pocketoma’ utilizando la herramienta visual desarrollada para tal tarea. Durante la exploración se encontró una curiosa e importante relación entre el sitio de unión de la hormona del receptor de estrógenos y una cavidad no caracterizada en la proteína Caspasa-3. Seguidamente, se realizó un estudio de docking con ligandos similares a la hormona y se procedió a realizar una extensiva dinámica molecular con el ligando en la mejor posición para verificar la estabilidad del compuesto en la cavidad de Caspasa-3. El complejo proteína-ligando estudiado resultó ser muy estable. Actualmente, el compuesto identificado se encuentra bajo validación experimental por nuestros colaboradores.

The number of experimentally resolved protein structures in the Protein Data Bank has been growing fast in the last 20 years, which motivates the develop- ment of many computational tools for protein-ligand binding sites prediction. Binding sites prediction from protein 3D structure has many important applica- tions; it is an essential step in the complex process of rational drug design, it helps to infer the side-effects of drugs, it provides insight into proteins biological functions and it is helpful in many other fields, such as protein-ligand docking and molecular dynamics. As far as we know, there has not been a study that would systematically investigate general properties of known ligand binding sites on a large scale. In this thesis, we examine these properties using existing experimen- tal and predicted residue-level annotations of protein sequence and structure. We present an automated pipeline for statistical analysis of these annotations, based on hypothesis testing and effect size estimation. It is implemented in Python and it is easily extensible by user-defined annotations. The usage is demonstrated on 33 existing annotations and 4 different datasets. The practical significance of the results is tested with P2Rank prediction method. We hope that the results as well as the pipeline...

Aminoacyl-tRNA synthetases (aaRSs) are at the center of the question of the origin of life and are essential proteins found in all living organisms. AARSs arose early in evolution to interpret genetic code and are believed to be a group of ancient proteins. They constitute a family of enzymes integrating the two levels of cellular organization: nucleic acids and proteins. These enzymes ensure the fidelity of transfer of genetic information from the DNA to the protein. They are responsible for attaching amino acid residues to their cognate tRNA molecules by virtue of matching the nucleotide triplet, which is the first step in the protein synthesis. The translation of genetic code into protein sequence is mediated by tRNA, which accurately picks up the cognate amino acids. The attachment of the cognate amino acid to tRNA is catalyzed by aaRSs, which have binding sites for the anticodon region of tRNA and for the amino acid to be attached. The two binding sites are separated by ≈ 76 Å and experiments have shown that the communication does not go through tRNA (Gale et al., 1996). The problem addressed here is how the information of binding of tRNA anticodon near the anticodon binding site is communicated to the active site through the protein structure. These enzymes are modular with distinct domains on which extensive kinetic and mutational experiments and supported by structural data are available, highlighting the role of inter-domain communication (Alexander and Schimmel, 2001). Hence these proteins present themselves as excellent systems for in-silico studies. Various methods involved for the construction of protein structure networks are well established and analyzed in a variety of ways to gain insights into different aspects of protein structure, stability and function (Kannan and Vishveshwara, 1999; Brinda and Vishveshwara, 2005). In the present study, we have incorporated network parameters for the analysis of molecular dynamics (MD) simulation data, representing the global dynamic behavior of protein in a more elegant way. MD simulations have been performed on the available (and modeled) structures of aaRSs bound to a variety of ligands, and the protein structure networks (PSN) of non-covalent interactions have been characterized in dynamical equilibrium. The changes in the structure networks are used to understand the mode of communication, and the paths between the two sites of interest identified by the analysis of the shortest path. The allosteric concept has played a key role in understanding the biological functions of aaRSs. The rigidity/plasticity and the conformational population are the two important ideas invoked in explaining the allosteric effect. We have explored the conformational changes in the complexes of aaRSs through novel parameters such as cliques and communities (Palla et al., 2005), which identify the rigid regions in the protein structure networks (PSNs) constructed from the non-covalent interactions of amino acid side chains. The thesis consists of 7 chapters. The first chapter constitutes the survey of the literature and also provides suitable background for this study. The aims of the thesis are presented in this chapter. Chapter 2 describes various techniques employed and the new techniques developed for the analysis of PSNs. It includes a brief description of well -known methods of molecular dynamics simulations, essential dynamics, and cross correlation maps. The method used for the construction of graphs and networks is also described in detail. The incorporation of network parameters for the analysis of MD simulation data are done for the first time and has been applied on a well studied protein lysozyme, as described in chapter 3. Chapter 3 focuses on the dynamical behavior of protein structure networks, examined by considering the example of T4-lysozyme. The equilibrium dynamics and the process of unfolding are followed by simulating the protein with explicit water molecules at 300K and at higher temperatures (400K, 500K) respectively. Three simulations of 10ns duration have been performed at 500K to ensure the validity of the results. The snapshots of the protein structure from the simulations are represented as Protein Structure Networks (PSN) of non-covalent interactions. The strength of the non-covalent interaction is evaluated and used as an important criterion in the construction of edges. The profiles of the network parameters such as the degree distribution and the size of the largest cluster (giant component) have been examined as a function of interaction strength (Ghosh et al., 2007). We observe a critical strength of interaction (Icritical) at which there is a transition in the size of the largest cluster. Although the transition profiles at all temperatures show behavior similar to those found in the crystal structures, the 500K simulations show that the non-native structures have lower Icritical values. Based on the interactions evaluated at Icritical value, the folding/unfolding transition region has been identified from the 500K simulation trajectories. Furthermore, the residues in the largest cluster obtained at interaction strength higher than Icritical have been identified to be important for folding. Thus, the compositions of the top largest clusters in the 500K simulations have been monitored to understand the dynamical processes such as folding/unfolding and domain formation/disruption. The results correlate well with experimental findings. In addition, the highly connected residues in the network have been identified from the 300K and 400K simulations and have been correlated with the protein stability as determined from mutation experiments. Based on these analyses, certain residues, on which experimental data is not available, have been predicted to be important for the folding and the stability of the protein. The method can also be employed as a valuable tool in the analysis of MD simulation data, since it captures the details at a global level, which may elude conventional pair-wise interaction analysis. After standardizing the concept of dynamical network analysis using Lysozyme, it was applied to our system of interest, the aaRSs. The investigations carried out on Methionyl-tRNA synthetases (MetRS) are presented in chapter 4. This chapter is divided into three parts: Chapter 4A deals with the introduction to aminoacyl tRNA synthetases (aaRS). Classification and functional insights of aaRSs obtained through various studies are presented. Chapter 4B is again divided into parts: BI and BII. Chapter 4BI elucidates a new technique developed for finding communication pathways essential for proper functioning of aaRS. The enzymes of the family of tRNA synthetases perform their functions with high precision, by synchronously recognizing the anticodon region and the amino acylation region, which is separated by about 70Å in space. This precision in function is brought about by establishing good communication paths between the two regions. We have modelled the structure of E.coli Methionyl tRNA synthetase, which is complexed with tRNA and activated methionine. Molecular dynamics simulations have been performed on the modeled structure to obtain the equilibrated structure of the complex and the cross correlations between the residues in MetRS. Furthermore, the network analysis on these structures has been carried out to elucidate the paths of communication between the aminoacyl activation site and the anticodon recognition site (Ghosh and Vishveshwara, 2007). This study has provided the detailed paths of communication, which are consistent with experimental results. A similar study on the (MetRS + activated methionine) and (MetRS+tRNA) complexes along with ligand free-native enzyme has also been carried out. A comparison of the paths derived from the four simulations has clearly shown that the communication path is strongly correlated and unique to the enzyme complex, which is bound to both the tRNA and the activated methionine. The method developed here could also be utilized to investigate any protein system where the function takes place through long distance communication. The details of the method of our investigation and the biological implications of the results are presented in this chapter. In chapter 4BII, we have explored the conformational changes in the complexes of E.coli Methionyl tRNA synthetase (MetRS) through novel parameters such as cliques and communities, which identify the rigid regions in the protein structure networks (PSNs). The rigidity/plasticity and the conformational population are the two important ideas invoked in explaining the allosteric effect. MetRS belongs to the aminoacyl tRNA Synthetases (aaRSs) family that play a crucial role in initiating the protein synthesis process. The network parameters evaluated here on the conformational ensembles of MetRS complexes, generated from molecular dynamics simulations, have enabled us to understand the inter-domain communication in detail. Additionally, the characterization of conformational changes in terms of cliques/communities has also become possible, which had eluded conventional analyses. Furthermore, we find that most of the residues participating in clique/communities are strikingly different from those that take part in long-range communication. The cliques/communities evaluated here for the first time on PSNs have beautifully captured the local geometries in their detail within the framework of global topology. Here the allosteric effect is revealed at the residue level by identifying the important residues specific for structural rigidity and functional flexibility in MetRS. Chapter 4C focuses on MD simulations of Methionyl tRNA synthetase (AmetRS) from a thermophilic bacterium, Aquifex aeolicus. As describe in Chapter 4B, we have explored the communication pathways between the anticodon binding region and the aminoacylation site, and the conformational changes in the complexes through cliques and communities. The two MetRSs from E.coli and Aquifex aeolicus are structurally and sequentially very close to each other. But the communication pathways between anticodon binding region and the aminoacylation site from A. aeolicus have differed significantly with the communication paths obtained from E.coli. The residue composition and cliques/communities structure participating in communication are not similar in the MetRSs of both these organisms. Furthermore the formation of cliques/communities and hubs in the communication paths are more in A. aeolicus compared to E.coli. The participation of structurally homologous linker peptide, essential for orienting the two domains for efficient communication is same in both the organisms although, the residues composition near domain interface regions including the linker peptide is different. Thus, the diversity in the functioning of two different MetRS has been brought out, by comparing the E.coli and Aquifex aeolicus systems. Protein Structure network analysis of MD simulated trajectories of various ligand bound complexes of Escherichia coli Cysteinyl-tRNA synthetase (CysRS) have been discussed in Chapter 5. The modeling of the complex is done by docking the ligand CysAMP into the tRNA bound structure of E.coli Cysteinyl tRNA synthetase. Molecular dynamics simulations have been performed on the modeled structure and the paths of communications were evaluated using a similar method as used in finding communication paths for MetRS enzymes. Compared to MetRS the evaluation of communication paths in CysRS is complicated due to presence of both direct and indirect readouts. The direct and indirect readouts (DR/IR) involve interaction of protein residues with base-specific functional group and sugar-phosphate backbone of nucleic acids respectively. Two paths of communication between the anticodon region and the activation site has been identified by combining the cross correlation information with the protein structure network constructed on the basis of non-covalent interaction. The complete paths include DR/IR interactions with tRNA. Cliques/communities of non-covalently interacting residues imparting structural rigidity are present along the paths. The reduction of cooperative fluctuation due to the presence of community is compensated by IR/DR interaction and thus plays a crucial role in communication of CysRS. Chapter 6 focuses on free energy calculations of aminoacyl tRNA synthetases with various ligands. The free energy contributions to the binding of the substrates are calculated using a method called MM-PBSA (Massova and Kollman, 2000). The binding free energies were calculated as the difference between the free energy of the enzyme-ligand complex, and the free ligand and protein. The ligand unbinding energy values obtained from the umbrella sampling MD correlates well with the ligand binding energies obtained from MM-PBSA method. Furthermore the essential dynamics was captured from MD simulations trajectories performed on E.coli MetRS, A. aeolius MetRS and E.coli CysRS in terms of the eigenvalues. The top two modes account for more than 50% of the motion in essential space for systems E.coli MetRS, A. aeolius MetRS and E.coli CysRS. Population distribution of protein conformation states are looked at the essential plane defined by the two principal components with highest eigenvalues. This shows how aaRSs existed as a population of conformational states and the variation with the addition of ligands. The population of conformational states is converted into Free energy contour surface. From free energy surfaces, it is evident that the E.coli tRNAMet bound MetRS conformational fluctuations are more, which attributes to less rigidity in the complex. Whereas E.coli tRNACys bound CysRS conformational fluctuations are less and this is reflected in the increase in rigidity of the complex as confirmed by its entropic contribution. Future directions have been discussed in the final chapter (Chapter 7). Specifically, it deals with the ab-initio QM/MM study of the enzymatic reaction involved in the active site of E.coli Methionyl tRNA synthetase. To achieve this, two softwares are integrated: the Quantum Mechanics (QM) part includes small ligands and the Molecular Mechanics (MM) part as protein MetRS are handled using CPMD and Gromacs respectively. The inputs for two reactions pathways are prepared. First reaction involves cyclization reaction of homocysteine in the active site of MetRS and the second reaction deals with charging of methionine in the presence of ATP and magnesium ion. These simulations require very high power computing systems and also time of computation is also very large. With the available computational power we could simulate up to 10ps and it is insufficient for analysis. The future direction will involve the simulations of these systems for longer time, followed by the analysis for reaction pathways.