Computational approaches in supramolecular chemistry with a special focus on virtual screening

Within this thesis novel computational tools for the rational design of synthetic host-guest complexes (SHGC) were developed and applied that employ the concepts of efficient virtual screening (VS) approaches. The first part describes the development of a fast structure prediction tool for flexible SHGC. The tool was validated on a test dataset comprising crystallographically determined SHGC. In nine of ten cases near-native solutions were generated. The tool can be applied for VS. In the second part of the thesis computational techniques were applied for designing SHGC based on ß-cyclodextrins (ß-CD). We performed a structure-based inverse virtual screening for identifying modified ß-CDs as receptors for the anticancer drug camptothecin (CPT). Six of the proposed receptors exhibited binding affinities which were significantly higher than for any other CPT-receptor. Furthermore, we applied a combination of a similarity-based virtual screening technique with a regression model (RM) for identifying novel high affinity guest molecules of ß-CD. Ten of the proposed guest molecules exhibited a binding free energy of lower than -20 kJ/mol. The last chapter describes a comparison of regression methods regarding their ability to generate predictive RM for thermodynamical parameters (dG, dH and dS) of ß-CD-guest complexes. dG could be predicted in good agreement with experimental values, none of the methods led to comparably good predictive models for dH. dS appears almost unpredictable.Im Rahmen dieser Arbeit wurden rechnergestützte Verfahren (RGV) zum gezielten Entwurf von synthetischen Wirt-Gast Komplexen (SWGK) entwickelt und eingesetzt. Dabei wurde ein Fokus auf schnelle virtuelle Screening (VS) Verfahren gelegt. Der erste Teil beschreibt die Entwicklung eines Programms zur schnellen Strukturvorhersage von flexiblen SWGK. Das Programm wurde auf einem Testdatensatz an kristallographisch vermessenen SWGK validiert. Für neun von zehn SWGK wurden nativ-ähnliche Lösungen gefunden. Das Programm kann für VS eingesetzt werden. Der zweite Teil der Arbeit behandelt RGV zum gezielten Entwurf von ß-Cyclodextrin (ß-CD) Komplexen. Mit Hilfe eines strukturbasierten inversen VS wurden sechs modifizierte ß-CD-Rezeptoren für den Krebsarzneistoff Camptothecin (CPT) gefunden, die deutlich höhere Bindungsaffinitäten zu CPT aufwiesen als alle bislang bekannten CPT-Rezeptoren. Zur Identifizierung neuer hochaffiner Gäste von ß-CD wurde ein ähnlichkeitsbasiertes VS Verfahren in Kombination mit einem Regressionsmodell (RM) eingesetzt. Zehn der mit Hilfe dieses Verfahrens vorgeschlagenen Moleküle wiesen eine Bindungsenergie von unter -20 kJ/mol auf. Das letzte Kapitel beschreibt einen Vergleich von drei Regressionsverfahren. Es wurde die Fähigkeit untersucht, vorhersagekräftige RM für thermodynamische Parameter (dG, dH und dS) von ß-CD-Gast-Komplexen zu generieren. dG konnte mit allen Methoden sehr gut vorhergesagt werden, während dH nur begrenzt und dS unzureichend vorhersagbar war

Bridging Between Protein Dynamics and Evolution Through Simulations and Machine Learning Approaches

Antibiotics resistance posed a serious threat to the public health and caused huge economic cost. β-Lactamases, which are enzymes produced by bacteria to hydrolyze β-lactam based antibiotics, are one of the driving forces behind antibiotic resistance. To explore the antibiotic resistance effect, understanding the mechanistic and dynamical features of β-lactamases through their interactions with antibiotics is critical. In my doctoral research, I applied both molecular dynamic (MD) simulations and machine learning approaches to explore these crucial interactions. Vancomycin is a typical glycopeptide antibiotic, which inhibits the bacterial cell wall through binding with peptidoglycan (PG). The key interactions of vancomycin and cell wall structure are identified by the conformational distributions of vancomycin and its three derivatives with PG complexes. TEM-1 is a serine-based β-lactamase and can hydrolyze the benzyl penicillin antibiotic. The key residues on TEM-1 are identified by random forest classification models. Moreover, the dynamical motions of four antibiotic resistance related proteins TEM-1, TOHO-1, PBP-A and DD-transpeptidase with a benzyl penicillin are analyzed and compared to explore their evolutionary correlation. I also investigated the petroleum thermal cracking mechanism through quantum chemistry calculations, and provided a quantitative and insightful understanding of thermal cracking processes

Development of Computational Methods to Predict Protein Pocket Druggability and Profile Ligands using Structural Data

This thesis presents the development of computational methods and tools using as input three-dimensional structures data of protein-ligand complexes. The tools are useful to mine, profile and predict data from protein-ligand complexes to improve the modeling and the understanding of the protein-ligand recognition. This thesis is divided into five sub-projects. In addition, unpublished results about positioning water molecules in binding pockets are also presented. I developed a statistical model, PockDrug, which combines three properties (hydrophobicity, geometry and aromaticity) to predict the druggability of protein pockets, with results that are not dependent on the pocket estimation methods. The performance of pockets estimated on apo or holo proteins is better than that previously reported in the literature (Publication I). PockDrug is made available through a web server, PockDrug-Server (http://pockdrug.rpbs.univ-paris-diderot.fr), which additionally includes many tools for protein pocket analysis and characterization (Publication II). I developed a customizable computational workflow based on the superimposition of homologous proteins to mine the structural replacements of functional groups in the Protein Data Bank (PDB). Applied to phosphate groups, we identified a surprisingly high number of phosphate non-polar replacements as well as some mechanisms allowing positively charged replacements. In addition, we observed that ligands adopted a U-shape conformation at nucleotide binding pockets across phylogenetically unrelated proteins (Publication III). I investigated the prevalence of salt bridges at protein-ligand complexes in the PDB for five basic functional groups. The prevalence ranges from around 70% for guanidinium to 16% for tertiary ammonium cations, in this latter case appearing to be connected to a smaller volume available for interacting groups. In the absence of strong carboxylate-mediated salt bridges, the environment around the basic functional groups studied appeared enriched in functional groups with acidic properties such as hydroxyl, phenol groups or water molecules (Publication IV). I developed a tool that allows the analysis of binding poses obtained by docking. The tool compares a set of docked ligands to a reference bound ligand (may be different molecule) and provides a graphic output that plots the shape overlap and a Jaccard score based on comparison of molecular interaction fingerprints. The tool was applied to analyse the docking poses of active ligands at the orexin-1 and orexin-2 receptors found as a result of a combined virtual and experimental screen (Publication V). The review of literature focusses on protein-ligand recognition, presenting different concepts and current challenges in drug discovery.Tässä väitöskirjassa esitetään tietokoneavusteisia menetelmiä ja työkaluja, jotka perustuvat proteiini-ligandikompleksien kolmiulotteisiin rakenteisiin. Ne soveltuvat proteiini-ligandikompleksien rakennetiedon louhimiseen, optimointiin ja ennustamiseen. Tavoitteena on parantaa sekä mallinnusta että käsitystä proteiini-liganditunnistuksesta. Väitöskirjassa työkalut kuvataan viitenä eri alahankkeena. Lisäksi esitetään toistaiseksi julkaisemattomia tuloksia vesimolekyylien asemoinnista proteiinien sitoutumistaskuihin. Kehitin PockDrugiksi kutsumani tilastollisen mallin, joka yhdistää kolme ominaisuutta – hydrofobisuuden, geometrian ja aromaattisuuden – proteiinitaskujen lääkekehityskohteeksi soveltuvuuden ennustamista varten siten, että tulokset ovat riippumattomia sitoutumistaskun sijoitusmenetelmästä. Apo- ja holoproteiinien taskujen ennustaminen toimii paremmin kuin alan kirjallisuudessa on aiemmin kuvattu (Julkaisu I). PockDrug on vapaasti käyttäjien saatavilla PockDrug-verkkopalvelimelta (http://pockdrug.rpbs.univ-paris-diderot.fr), jossa on lisäksi useita työkaluja proteiinin sitoutumiskohdan analyysiin ja karakterisointiin (Julkaisu II). Kehitin myös muokattavissa olevan tietokoneavusteisen prosessin, joka perustuu samankaltaisten proteiinien päällekkäin asetteluun, louhiakseni Protein Data Bankista (PDB) toiminnallisten ryhmien rakenteellisia korvikkeita. Tätä fosfaattiryhmiin soveltaessani tunnistin yllättävän paljon poolittomia fosfaattiryhmän korvikkeita ja joitakin positiivisesti varautuneita korvikkeita mahdollistavia mekanismeja. Lisäksi havaitsin, että ligandit omaksuivat U muotoisen konformaation fylogeneettisesti riippumattomien proteiinien nukleotidien sitoutumistaskuissa (Julkaisu III). Tutkin PDB:n proteiini-ligandikompleksien suolasiltojen yleisyyttä viidelle emäksiselle toiminnalliselle ryhmälle. Suolasiltojen yleisyys vaihteli guanidinium-ionin 70 prosentista tertiääristen ammoniumkationien 16 prosenttiin. Jälkimmäisessä tapauksessa suolasiltojen vähäisyys vaikuttaa riippuvan siitä, että vuorovaikuttaville ryhmille on vähemmän tilaa. Mikäli tarkastellut emäksiset ryhmät eivät osallistuneet vahvoihin karboksylaattivälitteisiin suolasiltoihin, niiden ympäristössä vaikutti olevan runsaasti happamia toiminnallisia ryhmiä, kuten hydroksi- ja fenoliryhmiä sekä vesimolekyylejä (Julkaisu IV). Lopuksi kehitin työkalun, joka mahdollistaa telakoinnista saatujen sitoutumisasentojen analyysin. Työkalu vertaa telakoitua ligandisarjaa sitoutuneeseen vertailuligandiin, joka voi olla eri molekyyli. Graafisena tulosteena saadaan diagrammi ligandien muotojen samankaltaisuudesta ja molekyylivuorovaikutusten sormenjälkiin perustuvasta Jaccard-pistemäärästä. Työkalua sovellettiin oreksiini-1- ja oreksiini-2-reseptoreille yhdistetyllä virtuaalisella ja kokeellisella seulonnalla löydettyjen aktiivisten ligandien sitoutumisasentojen analyysiin (Julkaisu V).Cette thèse présente le développement de méthodes et d’outils informatiques basés sur la structure tridimensionnelle des complexes protéine-ligand. Ces différentes méthodes sont utilisées pour extraire, optimiser et prédire des données à partir de la structure des complexes afin d’améliorer la modélisation et la compréhension de la reconnaissance entre une protéine et un ligand. Ce travail de thèse est divisé en cinq projets. En complément, une étude sur le positionnement des molécules d’eau dans les sites de liaisons a aussi été développée et est présentée. Dans une première partie un modèle statistique, PockDrug, a été mis en place. Il combine trois propriétés de poches protéiques (l’hydrophobicité, la géométrie et l’aromaticité) pour prédire la druggabilité des poches protéiques, si une poche protéique peut lier une molécule drug-like. Le modèle est optimisé pour s’affranchir des différentes méthodes d’estimation de poches protéiques. La qualité des prédictions, est meilleure à la fois sur des poches estimées à partir de protéines apo et holo et est supérieure aux autres modèles de la littérature (Publication I). Le modèle PockDrug est disponible sur un serveur web, PockDrug-Server (http://pockdrug.rpbs.univ-paris-diderot.fr) qui inclus d’autres outils pour l’analyse et la caractérisation des poches protéiques. Dans un second temps un protocole, basé sur la superposition de protéines homologues a été développé pour extraire des replacements structuraux de groupements chimiques fonctionnels à partir de la Protein Data Bank (PDB). Appliqué aux phosphates, un grand nombre de remplacements non-polaires ont été identifié pouvant notamment être chargés positivement. Quelques mécanismes de remplacements ont ainsi pu être analysé. Nous avons, par exemple, observé que le ligand adopte une configuration en forme U dans les sites de liaison des nucléotides indépendamment de la phylogénétique des protéines (Publication III). Dans une quatrième partie, la prévalence des ponts salins de cinq groupements chimiques basiques a été étudié dans les complexes protéine-ligand. Ainsi le pourcentage de pont salin fluctue de 70% pour le guanidinium à 16% pour l’amine tertiaire qui a le plus faible volume disponible autour de lui pour accueillir un group pouvant interagir. L’absence d’acide fort comme l’acide carboxylique pour former un pont salin est remplacé par un milieu enrichis en groupement chimiques fonctionnels avec des propriétés acides comme l’hydroxyle, le phénol ou encore les molécules d’eau (Publication IV). Dans un dernier temps un outil permettant l’analyse des poses de ligand obtenues par une méthode d’ancrage moléculaire a été développé. Cet outil compare ces poses à un ligand de référence, qui peut être une molécule différente en combinant l’information du chevauchement de forme de la pose et du ligand de référence et un score de Jaccard basé sur une comparaison des empreintes d’interaction moléculaires du ligand de référence et de la pose. Cette méthode a été utilisé dans l’analyse des résultats d’ancrage moléculaires pour des ligands actifs pour les récepteurs aux orexine 1 et 2. Ces ligands actifs ont été trouvés à partir de résultats combinant un criblage virtuel et expérimental. La revue de la littérature associée est focalisée sur la reconnaissance moléculaire d’un ligand pour une protéine et présente diffèrent concepts et challenges pour la recherche de nouveaux médicaments

Predicting biomolecular function from 3D dynamics : sequence-sensitive coarse-grained elastic network model coupled to machine learning

La dynamique structurelle des biomolécules est intimement liée à leur fonction, mais très coûteuse à étudier expériementalement. Pour cette raison, de nombreuses méthodologies computationnelles ont été développées afin de simuler la dynamique structurelle biomoléculaire. Toutefois, lorsque l'on s'intéresse à la modélisation des effects de milliers de mutations, les méthodes de simulations classiques comme la dynamique moléculaire, que ce soit à l'échelle atomique ou gros-grain, sont trop coûteuses pour la majorité des applications. D'autre part, les méthodes d'analyse de modes normaux de modèles de réseaux élastiques gros-grain (ENM pour "elastic network model") sont très rapides et procurent des solutions analytiques comprenant toutes les échelles de temps. Par contre, la majorité des ENMs considèrent seulement la géométrie du squelette biomoléculaire, ce qui en fait de mauvais choix pour étudier les effets de mutations qui ne changeraient pas cette géométrie. Le "Elastic Network Contact Model" (ENCoM) est le premier ENM sensible à la séquence de la biomolécule à l'étude, ce qui rend possible son utilisation pour l'exploration efficace d'espaces conformationnels complets de variants de séquence. La présente thèse introduit le pipeline computationel ENCoM-DynaSig-ML, qui réduit les espaces conformationnels prédits par ENCoM à des Signatures Dynamiques qui sont ensuite utilisées pour entraîner des modèles d'apprentissage machine simples. ENCoM-DynaSig-ML est capable de prédire la fonction de variants de séquence avec une précision significative, est complémentaire à toutes les méthodes existantes, et peut générer de nouvelles hypothèses à propos des éléments importants de dynamique structurelle pour une fonction moléculaire donnée. Nous présentons trois exemples d'étude de relations séquence-dynamique-fonction: la maturation des microARN, le potentiel d'activation de ligands du récepteur mu-opioïde et l'efficacité enzymatique de l'enzyme VIM-2 lactamase. Cette application novatrice de l'analyse des modes normaux est rapide, demandant seulement quelques secondes de temps de calcul par variant de séquence, et est généralisable à toute biomolécule pour laquelle des données expérimentale de mutagénèse sont disponibles.The dynamics of biomolecules are intimately tied to their functions but experimentally elusive, making their computational study attractive. When modelling the effects of thousands of mutations, time-stepping methods such as classical or enhanced sampling molecular dynamics are too costly for most applications. On the other hand, normal mode analysis of coarse-grained elastic network models (ENMs) provides fast analytical dynamics spanning all timescales. However, the vast majority of ENMs consider backbone geometry alone, making them a poor choice to study point mutations which do not affect the equilibrium structure. The Elastic Network Contact Model (ENCoM) is the first sequence-sensitive ENM, enabling its use for the efficient exploration of full conformational spaces from sequence variants. The present work introduces the ENCoM-DynaSig-ML computational pipeline, in which the ENCoM conformational spaces are reduced to Dynamical Signatures and coupled to simple machine learning algorithms. ENCoM-DynaSig-ML predicts the function of sequence variants with significant accuracy, is complementary to all existing methods, and can generate new hypotheses about which dynamical features are important for the studied biomolecule's function. Examples given are the maturation efficiency of microRNA variants, the activation potential of mu-opioid receptor ligands and the effect of point mutations on VIM-2 lactamase's enzymatic efficiency. This novel application of normal mode analysis is very fast, taking a few seconds CPU time per variant, and is generalizable to any biomolecule on which experimental mutagenesis data exist

Efficient comprehensive scoring of docked proteincomplexes - a machine learning approach

Biological systems and processes rely on a complex network of molecular interactions. The association of biological macromolecules is a fundamental biochemical phenomenon and an unsolved theoretical problem crucial for the understanding of complex living systems. The term protein-protein docking describes the computational prediction of the assembly of protein complexes from the individual subunits. Docking algorithms generally produce a large number of putative protein complexes. In most cases, some of these conformations resemble the native complex structure within an acceptable degree of structural similarity. A major challenge in the field of docking is to extract the near-native structure(s) out of this considerably large pool of solutions, the so called scoring or ranking problem. It has been the aim of this work to develop methods for the efficient and accurate detection of near-native conformations in the scoring or ranking process of docked protein-protein complexes. A series of structural, chemical, biological and physical properties are used in this work to score docked protein-protein complexes. These properties include specialised energy functions, evolutionary relationship, class specific residue interface propensities, gap volume, buried surface area, empiric pair potentials on residue and atom level as well as measures for the tightness of fit. Efficient comprehensive scoring functions have been developed using probabilistic Support Vector Machines in combination with this array of properties on the largest currently available protein-protein docking benchmark. The established scoring functions are shown to be specific for certain types of protein-protein complexes and are able to detect near-native complex conformations from large sets of decoys with high sensitivity. The specific complex classes are Enzyme-Inhibitor/Substrate complexes, Antibody-Antigen complexes and a third class denoted as "Other" complexes which holds all test cases not belonging to either of the two previous classes. The three complex class specific scoring functions were tested on the docking results of 99 complexes in their unbound form for the above mentioned categories. Defining success as scoring a 'true' result with a p-value of better than 0.1, the scoring schemes were found to be successful in 93%, 78% and 63% of the examined cases, respectively. The ranking of near-native structures can be drastically improved, leading to a significant enrichment of near-native complex conformations in the top ranks. It could be shown that the developed scoring schemes outperform five other previously published scoring functions

Ligand recognition by the major urinary protein

Molecular Dynamics (MD) and Quartz Crystal Microbalance (QCM) techniques can provide unique insights into what drives protein-ligand association. The major urinary protein (MUP) binds small ligands in a deeply buried hydrophobic pocket. Detailed calorimetric studies have shown that ligand binding is driven by enthalpic effects, not entropic effects [1]. Previous studies have shown that this is due to 'dewetting' of the binding site cavity even in the absence of ligands, and have also characterised the complex changes in molecular flexibility that accompany ligand binding-features that may be correlated with NMR data [2]. Recent MD revealed the hydration effects of apo-MUP and also shown where certain regions of MUP become more flexible upon ligand binding. They have also shown a water molecule remains close to the tyrosine in the binding pocket [2]. In our current MD studies and OCM experiments we have used wild type and 2 different mutants of MUP to study the binding effects of the ligand IBM. The first mutant has an OH group removed from the binding site of MUP (i.e. tyrosine to phenylalanine (Y120F)). The second mutant has an extra OH group in the binding site (i.e. alanine to serine (A103S)). For all three systems the hydration and flexibility upon ligand binding has been analysed. The hydration analysis from MD reveal (from radial distribution curves and hydration density maps) there is a small density of water that remains even without the presence of the ligand for the WT MUP whereas a larger density of water remains in the binding cavity of the A103S hydrophilic MUP simulation. The results are based on the average structure generated from the 1 mus simulations. The Y120F MUP simulations reveal that there is no water molecules present in the binding cavity. However, as protein molecules are very dynamic in nature, water molecules are observed to hop in and out of the binding pockets for both mutant MUP (but not WT MUP) simulations over the 1 mus simulations. On the other hand the experimental QCM results reveal that on ligand binding no water loss is observed for Y120F mutant MUP whereas A103S and WT MUP have about 2 water molecules which are lost in the binding cavity. The flexibility results from the MD simulations reveal that WT MUP have some residues which increase in flexibility whilst other residues which decrease in flexibility on ligand binding. However, the Y120F hydrophobic MUP show an overall decrease in flexibility whereas the A103S MUP shows an overall increase in flexibility on ligand binding. In contrast the experimental OCM and AFM results reveal that there is an increase in flexibility on ligand binding to all 3 different types of MUP molecules. The experimental and the simulation data have shown a variation in results but it is to be noted that the results cannot be directly compared as the analytical experiments are a surface based techniques whereas the MD simulations do not involve a surface. However, the contrast observed between computer simulation and experiments has revealed important information on the ligand binding effects on MUP. [1] Bingham, R.J., J.B.C. Findlay, S.Y. Hsieh, A.P. Kalverda, A. Kjeliberg, C. Perazzolo, S.E.V. Phillips, K. Seshadri, C.H. Trinh, W. B. TurnbulI, G. Bodenhausen, and S.W. Homans. 2004. Thermodynamics of binding of 2-methoxy-3-lsopropylpyrazlne and 2- methoxy-3-lsobutylpyrazine to the major urinary protein. J. Am. Chem. Soc. 126:1675-1681. [2] Barratt, E., R.J. Bingham. D.J. Warner, C.A. Laughton, S.E.V. Phillips, and S.W. Homans. 2005. Van der Waals interactions dominate ligand-protein association in a protein binding site occluded from solvent water. J. Am. Chem. Soc. 127:11827-11834

Studying protein-ligand interactions using a Monte Carlo procedure

[eng] Biomolecular simulations have been widely used in the study of protein-ligand interactions; comprehending the mechanisms involved in the prediction of binding affinities would have a significant repercussion in the pharmaceutical industry. Notwithstanding the intrinsic difficulty of sampling the phase space, hardware and methodological developments make computer simulations a promising candidate in the resolution of biophysically relevant problems. In this context, the objective of the thesis is the development of a protocol that permits studying protein-ligand interactions, in view to be applied in drug discovery pipelines. The author contributed to the rewriting PELE, our Monte Carlo sampling procedure, using good practices of software development. These involved testing, improving the readability, modularity, encapsulation, maintenance and version control, just to name a few. Importantly, the recoding resulted in a competitive cutting-edge software that is able to integrate new algorithms and platforms, such as new force fields or a graphical user interface, while being reliable and efficient. The rest of the thesis is built upon this development. At this point, we established a protocol of unbiased all-atom simulations using PELE, often combined with Markov (state) Models (MSM) to characterize the energy landscape exploration. In the thesis, we have shown that PELE is a suitable tool to map complex mechanisms in an accurate and efficient manner. For example, we successfully conducted studies of ligand migration in prolyl oligopeptidases and nuclear hormone receptors (NHRs). Using PELE, we could map the ligand migration and binding pathway in such complex systems in less than 48 hours. On the other hand, with this technique we often run batches of 100s of simulations to reduce the wall-clock time. MSM is a useful technique to join these independent simulations in a unique statistical model, as individual trajectories only need to characterize the energy landscape locally, and the global characterization can be extracted from the model. We successfully applied the combination of these two methodologies to quantify binding mechanisms and estimate the binding free energy in systems involving NHRs and tyorsinases. However, this technique represents a significant computational effort. To reduce the computational load, we developed a new methodology to overcome the sampling limitations caused by the ruggedness of the energy landscape. In particular, we used a procedure of iterative simulations with adaptive spawning points based on reinforcement learning ideas. This permits sampling binding mechanisms at a fraction of the cost, and represents a speedup of an order of magnitude in complex systems. Importantly, we show in a proof-of-concept that it can be used to estimate absolute binding free energies. Overall, we hope that the methodologies presented herein help streamline the drug design process.[spa] Las simulaciones biomoleculares se han usado ampliamente en el estudio de interacciones proteína-ligando. Comprender los mecanismos involucrados en la predicción de afinidades de unión tiene una gran repercusión en la industria farmacéutica. A pesar de las dificultades intrínsecas en el muestreo del espacio de fases, mejoras de hardware y metodológicas hacen de las simulaciones por ordenador un candidato prometedor en la resolución de problemas biofísicos con alta relevancia. En este contexto, el objetivo de la tesis es el desarrollo de un protocolo que introduce un estudio más eficiente de las interacciones proteína-ligando, con vistas a diseminar PELE, un procedimiento de muestreo de Monte Carlo, en el diseño de fármacos. Nuestro principal foco ha sido sobrepasar las limitaciones de muestreo causadas por la rugosidad del paisaje de energías, aplicando nuestro protocolo para hacer analsis detallados a nivel atomístico en receptores nucleares de hormonas, receptores acoplados a proteínas G, tirosinasas y prolil oligopeptidasas, en colaboración con una compañía farmacéutica y de varios laboratorios experimentales. Con todo ello, esperamos que las metodologías presentadas en esta tesis ayuden a mejorar el diseño de fármacos