Flies are a widely distributed pest insect that poses a significant threat to food security. Flight is essential for the dispersal of the adult flies to find new food sources and ideal breeding spots. The supply of metabolic fuel to power the flight muscles of insects is regulated by adipokinetic hormones (AKHs). The fruit fly, Drosophila melanogaster, the flesh fly, Sarcophaga crassipalpis, and the oriental fruit fly, Bactrocera dorsalis all have the same AKH that is present in the blowfly, Phormia terraenovae; this AKH has the code-name Phote-HrTH. Binding of the AKH to the extracellular binding site of a G protein-coupled receptor causes its activation. In this thesis, the structure of Phote-HrTH in SDS micelle solution was determined using NMR restrained molecular dynamics. The peptide was found to bind to the micelle and be reasonably rigid, with an S 2 order parameter of 0.96. The translated protein sequence of the AKH receptor from the fruit fly, Drosophila melanogaster, the flesh fly, Sarcophaga crassipalpis, and the oriental fruit fly, Bactrocera dorsalis were used to construct two models for each receptor: Drome-AKHR, Sarcr-AKHR, and Bacdo-AKHR. It is proposed that these two models represent the active and inactive state of the receptor. The models based on the crystal structure of the β-2 adrenergic receptor were found to bind Phote-HrTH with a predicted binding free energy of –107 kJ mol–1 for Drome-AKHR, –102 kJ mol–1 for Sarcr-AKHR and –102 kJ mol–1 for Bacdo-AKHR. Under molecular dynamics simulation, in a POPC membrane, the β-2AR receptor-like complexes transformed to rhodopsin-like. The identification and characterisation of the ligand-binding site of each receptor provide novel information on ligand-receptor interactions, which could lead to the development of species-specific control substances to use discriminately against these pest flies.

<p>CYP2C9 is one of our major drug metabolising enzymes and belongs to the cytochrome P450 (CYP) super family. The aim of this thesis was to gain an understanding of the quantitative structure–activity relationships (QSAR) of CYP2C9 substrates and inhibitors. This information will be useful in predicting drug metabolism and the potential for drug–drug interactions. To achieve this, a well characterised data set of structurally diverse, competitive CYP2C9 inhibitors was identified in our laboratory. Several computational methodologies, many based on GRID molecular interaction fields, were applied or developed in order to handle issues such as compound alignment and bioactive conformer selection. First, a traditional 3D QSAR was carried out in GOLPE, generating a predictive model. In this model the selection of a bioactive conformer and alignment was based on docking in a homology model of CYP2C9. Secondly, we introduced the concept of alignment independent descriptors from ALMOND. These descriptors were used to generate quantitatively and qualitatively predictive models. We subsequently derived conformation independent descriptors from molecular interaction fields calculated in FlexGRID. This enabled the derivation of 3D QSAR models without taking into account the selection of an alignment or a bioactive conformer. A subsequent programming effort enabled the conversion of this model back to 3D aligned pharmacophores. Similar alignment independent descriptors were also used in the development of the software MetaSite® that predicts the site of metabolism for CYP2C9 ligands. Finally, as crystal information on this isoform emerged, the performance of molecular dynamics simulations and homology models and the flexibility of the protein were evaluated using statistical analyses.</p><p>These modelling efforts have resulted in detailed knowledge of the structural characteristics in ligand interactions with the cytochrome P450 2C9 isoform.</p>

CYP2C9 is one of our major drug metabolising enzymes and belongs to the cytochrome P450 (CYP) super family. The aim of this thesis was to gain an understanding of the quantitative structure–activity relationships (QSAR) of CYP2C9 substrates and inhibitors. This information will be useful in predicting drug metabolism and the potential for drug–drug interactions. To achieve this, a well characterised data set of structurally diverse, competitive CYP2C9 inhibitors was identified in our laboratory. Several computational methodologies, many based on GRID molecular interaction fields, were applied or developed in order to handle issues such as compound alignment and bioactive conformer selection. First, a traditional 3D QSAR was carried out in GOLPE, generating a predictive model. In this model the selection of a bioactive conformer and alignment was based on docking in a homology model of CYP2C9. Secondly, we introduced the concept of alignment independent descriptors from ALMOND. These descriptors were used to generate quantitatively and qualitatively predictive models. We subsequently derived conformation independent descriptors from molecular interaction fields calculated in FlexGRID. This enabled the derivation of 3D QSAR models without taking into account the selection of an alignment or a bioactive conformer. A subsequent programming effort enabled the conversion of this model back to 3D aligned pharmacophores. Similar alignment independent descriptors were also used in the development of the software MetaSite® that predicts the site of metabolism for CYP2C9 ligands. Finally, as crystal information on this isoform emerged, the performance of molecular dynamics simulations and homology models and the flexibility of the protein were evaluated using statistical analyses. These modelling efforts have resulted in detailed knowledge of the structural characteristics in ligand interactions with the cytochrome P450 2C9 isoform.

Cyclotides are small cyclic plant proteins, and this thesis addresses their cytotoxic structure-activity properties and their mode of action on human cancer cell lines. Cyclotides were isolated from Viola odorata and Viola tricolor; three novel cyclotide sequences and two known sequences, but of new origin, were identified using mass spectrometry, amino acid analysis, and Edman degradation. The cyclotide structure includes three disulphide bonds in a knotted arrangement, which forces hydrophobic amino acid residues to be exposed on the surface of the molecule; 3-D homology models of cyclotides have revealed an amphipathic surface and charged residues located at similar positions in the molecules. The charged amino acid residues were shown to play a key role in the cytotoxicity of the cyclotide cycloviolacinO2 on a human lymphoma cell line. Methylation of Glu caused a dramatic change in cytotoxicity, lowering the potency 48 times, whereas concealing the charge of Arg with 1,2-cyclohexanedione caused virtually no change in potency. Acetylation of the two Lys caused a 3-fold reduction in potency, and masking all positive charges caused a 7-fold reduction. Additionally, disturbing the amphipathic structure by reducing and alkylating the disulphide bonds abolished the cytotoxicity. The time dependency of cytotoxicity and cell gross morphology after cyclotide exposure were investigated on the lymphoma cell line. Cells exposed to 4 µM of cycloviolacinO2 showed necrotic characteristics, such as membrane disintegration, within 5 min; a membrane disruptive effect of cycloviolacinO2 was also observed in a functional assay based on liposomes at a peptide-to-lipid molar ratio of 6.5. The anti-tumour properties of cycloviolacinO2 were evaluated on three human cancer cell lines using the hollow fibre assay in vitro and in vivo. The cyclotide exhibited potent anti-tumour activity in the micro-molar concentration range on all cell lines in vitro, but no effect on tumour growth could be established in vivo.

Η ριβονουκλεάση P είναι το ένζυμο το οποίο αναλαμβάνει την δημιουργία του 5´ ώριμου άκρου όλων των πρόδρομων μορίων tRNA. Πρόκειται για ένα ριβονουκλεο-πρωτεϊνικό σύμπλοκο το οποίο εντοπίζεται στα κύτταρα των οργανισμών και από τις τρεις κύριες φυλογενετικές περιοχές, τα Βακτήρια, τα Αρχαία και τους Ευκαρυώτες. Αποτελείται από μια υπομονάδα RNA απαραίτητη για την κατάλυση, ενώ το μέγεθος και ο αριθμός των πρωτεϊνικών υπομονάδων ποικίλλει από μια μικρή στα βακτήρια έως δέκα πρωτεΐνες στο ολοένζυμο που απομονώνεται από τα ανθρώπινα κύτταρα. Η υπομονάδες RNA των βακτηρίων καθώς επίσης και μερικών αρχαίων μπορούν να καταλύσουν την αντίδραση ωρίμανσης του tRNA απουσία της πρωτεΐνης in vitro, είναι δηλαδή ριβοένζυμα. Η ανακάλυψη αυτή διεύρυνε τις αντιλήψεις μας για τις ιδιότητες των βιομορίων και επανέφερε στο προσκήνιο την θεωρία του κόσμου του RNA. Στο ευκαρυωτικό ριβοένζυμο, ο ρόλος των πρωτεϊνών είναι πιο ουσιαστικός, καθώς η υπομονάδα RNA φαίνεται ότι χάνει μεγάλο μέρος της λειτουργικής της ανεξαρτησίας. Η διαλεύκανση των λειτουργών της κάθε υπομονάδας θα δώσει σημαντικές πληροφορίες για την εξέλιξη της RNase P από ένα αρχέγονο ένζυμο σε ένα πολύπλοκο ριβονουκλεοπρωτεϊνικό σύμπλοκο. Η RNase P από το Dictyostelium discoideum διαθέτει μια απαραίτητη για την δραστικότητα υπομονάδα RNA όπως και όλα τα ένζυμα αυτού του είδους. Παράλληλα διαθέτει έντονο πρωτεϊνικό χαρακτήρα καθώς διαθέτει την χαμηλότερη πυκνότητα επιπολής σε σχέση με ένζυμα RNase P από άλλους οργανισμούς. Οι πληροφορίες αυτές προέρχονται από τον αρχικό χαρακτηρισμό του ενζυμικού συμπλόκου, και δεν παρέχουν στοιχεία για την ακριβή σύστασή του. Στην παρούσα μελέτη, πραγματοποιήθηκε κλωνοποίηση και χαρακτηρισμός ενός από τα γονίδια που εντοπίστηκαν στο γονιδίωμα του Dictyostelium, ομόλογα προς χαρακτηρισμένα γονίδια από τον άνθρωπο και άλλους ευκαρυώτες. Το γονίδιο drpp30 κωδικεύει μια πρωτεΐνη 40.7 kDa, σημαντικά μεγαλύτερη από τις ομόλογες της. Η πρωτεΐνη DRpp30 υπερεκφράστηκε σε βακτηριακά κύτταρα, και μετά τον χρωματογραφικό καθαρισμό της χρησιμοποιήθηκε για την παρασκευή πολυκλωνικών αντισωμάτων. Η συμμετοχή της DRpp30 στο μακρομοριακό σύμπλοκο της RNase P πιστοποιήθηκε με ανοσοβιοχημική προσέγγιση, ενώ η ανασυνδυασμένη πρωτεΐνη προσδένει τo pre-tRNA υπόστρωμα του ενζύμου, καθώς και την υπομονάδα RNA in vitrο. Το μοντέλο ομολογίας της DRpp30 βάσει της κρυσταλλικής δομής της ορθόλογης Ph1877 από τα αρχαία, φανερώνει ότι η πρωτεΐνη αποκτά τη δομή αβ βαρελιού (ΤΙΜ barrel fold). Κατά τη διάρκεια της διατριβής, οι προσπάθειες για τον εντοπισμό του γονιδίου της RNA υπομονάδας ήταν σε εξέλιξη, όταν το εν λόγω γονίδιο αναγνωρίστηκε μέσω φυλογενετικών συγκρίσεων από την ομάδα του Norman Pace. Το μετάγραφο του γονιδίου εντοπίστηκε σε ενεργά κλάσματα RNase P, και παράλληλα εντοπίστηκε και ένα μικρότερο μετάγραφο του ίδιου γονιδίου. Προσδιορίστηκαν τα ακριβή 5´και 3´ άκρα των δύο αυτών μορίων και ακολούθησε κλωνοποίηση τους. Τα in vitro μετάγραφα των δύο κλωνοποιημένων αλληλουχιών μπορούν να υποκαθιστούν την ενδογενή RNA υπομονάδα του ολοενζύμου in vitro, ενώ δεν εντοπίστηκε έως τώρα ενζυμική δραστικότητα που να σχετίζεται με τα δύο αυτά μόρια. / Ribonuclease P is a ubiquitus ribonucleoprotein enzyme, responsible for the production of the 5´ mature ends of all precursor tRNA molecules. RNase P endonucleolytic activity has been isolated from organisms representing the three domains of life, namely Bacteria, Archaea and Eukarya. It has been shown to contain an essential RNA subunit and one (Bacteria) or more (Archaea, Eukaryotes) proteins. The RNase P RNA subunits from bacteria and some archaea are catalytically active in vitro, whereas those from eukaryotes and most archaea have lost most of their functionality and require protein subunits for activity. RNase P has been characterized biochemically and genetically in several systems, and structures for both RNA and protein subunits have emerged. The integration of structural and functional data is slowly forming a scenario for the evolution of RNase P from an ancient enzyme to a highly organized ribonucleoprotein complex. Dictyostelium discoideum RNase P harbors an essential RNA subunit, and has high protein content, as judged by its low boyant density. Nevertheless, our knowledge on the exact composition was limited. In the current study, a gene showing significant similarity to human Rpp30 RNase P protein subunit was identified in Dictyostelium genome. The gene encodes a protein (DRpp30) which is significantly larger than its homologues, due to an unusual C-terminus. The gene was cloned, overexpressed, and was used for the production of polyclonal antibodies. The participation of DRpp30 in the macromolecular complex of RNase P was verified by an immunobiochemical approach. The recombinant protein was shown to bind specifically both the RNase P RNA subunit and the pre-tRNA substrate in vitro, thus giving a first insight of its role in the holoenzyme complex. Homology modeling using as a template the archaeal Ph1887p, and molecular dynamics simulations of the modeled structure suggest that DRpp30 adopts a TIM-barrel fold. While our efforts to isolate the gene encoding the RNA subunit of D. discoideum RNase P were in progress, Norman Pace and his group identified it through phylogenetic comparison. The full transcript of the gene was detected in active RNase P samples along with a smaller transcript of the same gene. The exact 5´and 3´ ends of both transcripts were identified and were cloned. Both these transcripts can substitute the endogenous RNA subunit in vitro, but no enzymatic activity associated with these RNA molecules could be detected so far.

Philosophiae Doctor - PhD / Tsetse flies are the biological vectors of human and animal trypanosomiasis and hence representant medical and veterinary importance. The sense of smell plays a significant role in tsetse and its ecological interaction, such as finding blood meal source, resting, and larvicidal sites and for mating. Tsetse olfactory behaviour can be exploited for their management; however, olfactory studies in tsetse flies are still fragmentary. Here in my PhD thesis, using scanning electron microscopy, electrophysiology, behaviour, bioinformatics and molecular biology techniques, I have investigated tsetse flies (Glossina fuscipes fuscipes) olfaction using behaviourally well studied odorants, tsetse repellent by comparing with attractant odour. Insect olfaction is mediated by olfactory sensory neurons (OSNs), located in olfactory sensilla, which are cuticular structures exposed to the environment through pore and create a platform for chemical communication. In the sensilla shaft the dendrite of OSNs are housed, which are protected by called the sensillum lymph produced by support cells and contains a variety of olfactory proteins, including the odorant binding protein (OBP) and chemosensory proteins (CSP). While on the dendrite of OSNs are expressed olfactory receptors. In my PhD, studies I tried to decipher the sense of smell in tsetse fly. In the second chapter, I demonstrated that G. f. fuscipes is equipped with diverse olfactory sensilla, that various from basiconic, trichoid and coeloconic. I also demonstrated, there is shape, length, number difference between sensilla types and sexual dimorphism. There is a major difference between male and female, while male has the unique basiconic sensilla, club shaped found in the pits, which is absent from female pits. In my third chapter, I investigated the odorant receptors which are expressed on the dendrite of the olfactory sensory neurons (OSNs). G. f. fuscipes has 42 ORs, which were not functionally characterised. I used behaviourally well studied odorants, tsetse repellents, composed of four components blend. I demonstrated that tsetse repellent is also a strong antifeedant for both G. pallidipes and G. f. fuscipes using feeding bioassays as compared to the attractant odour, adding the value of tsetse repellent. However, the attractant odour enhanced the feeding index. Using DREAM (deorphanization of receptors based on expression alterations of mRNA levels). I found that in G. f. fuscipes, following a short in vivo exposure to the individual tsetse repellent component as well as an attractant volatile chemical, OSNs that respond to these compounds altered their mRNA expression in two opposite direction, significant downregulation and upregulation in their number of transcripts corresponding to the OR that they expressed and interacted with odorant. Also, I found that the odorants with opposite valence already segregate distinctly at the cellular and molecular target at the periphery, which is the reception of odorants by OSNs, which is the basis of sophisticated olfactory behaviour. Deorphanization of ORs in none model insect is a challenge, here by combining DREAM with molecular dynamics, as docking score, physiology and homology modelling with Drosophila a well-studied model insects, I was able to predict putative receptors of the tsetse repellent components and an attractant odour. However, many ORs were neutral, showing they were not activated by the odorants, demonstrating the selectivity of the technique as well as the receptors. In my fourth chapter, I investigated the OBPs structures and their interaction with odorants molecules. I demonstrated that OBPs are expressed both in the antenna, as well as in other tissues, such as legs. I also demonstrated that there are variations in the expression of OBPs between tissues as well as sexes. I also demonstrated that odorants induced a fast alteration in OBP mRNA expression, some odorants induced a decrease in the transcription of genes corresponding to the activated OBP and others increased the expression by many fold in OBPs in live insect, others were neutral after 5 hours of exposure. Moreover, with subsequent behavioural data showed that the behavioural response of G. f. fuscipes toward 1-octen-3-ol decreased significantly when 1-octen-3-ol putative OBPs were silenced with feeding of double-stranded RNA (dsRNA). In summary, our finding whereby odorant exposure affects the OBPs mRNA, their physiochemical properties and the silencing of these OBPs affected the behavioural response demonstrate that the OBPs are involved in odour detection that affect the percept of the given odorant. The expression of OBPs in olfactory tissues, antenna and their interaction with odorant and their effect on behavioural response when silenced shows their direct involvement in odour detection and reception. Furthermore, their expression in other tissues such as legs indicates they might also have role in other physiological functions, such as taste.

Les mitochondries sont des organites dynamiques dont la morphologie dépend de l’équilibre fusion/fission de leurs membranes. Ce processus essentiel à la survie cellulaire est nommé dynamique mitochondriale et sa dérégulation est associée à des troubles neurologiques. Cependant les mécanismes précis régissant la dynamique mitochondriale ne sont pas élucidés. Cette thèse porte sur la protéine Fzo1p, une grande GTPase de la superfamille des Dynamin-related-Protein. C’est un élément clé impliqué dans la fusion mitochondriale de la membrane externe de la levure. Sa structure et sa dynamique ont été étudiées par modélisation et simulations de dynamiques moléculaires tout-atome dans une bicouche lipidique solvatée. Le modèle structural obtenu tient compte de données expérimentales, de template structuraux, et de modèles ab initio du domaine transmembranaire de Fzo1p. Ce modèle a été validé expérimentalement par mutagenèse dirigée. Des permutations de charges ont confirmé des ponts salins à longue distance prédits dans le modèle. En outre, des mutations ont montré que les domaines coiled-coil de Fzo1p, contrairement à sa partie N-terminale, sont indispensables à sa fonction. L’ensemble des résultats expérimentaux et in silico met en évidence l’implication des domaines charnières dans le changement conformationnel de Fzo1p, ainsi que des résidus critiques affectant sa stabilité. Les précisions atomiques obtenues sur l’interaction de Fzo1p avec le GDP permet de formuler des hypothèses sur le mécanisme moléculaire de la catalyse du GTP pour la fusion membranaire; voire à la compréhension de la dynamique mitochondriale. / Mitochondria are dynamic organelles whose morphology is determined by fusion and fission of their membranes. This essential process is known as mitochondrial dynamics. Defects in mitochondrial dynamics are associated with neurological disorders making the investigation of physiological relevance. However, the precise sequence of events that lead mitochondrial dynamics are still not well characterised. Fzo1p, a large GTPase of the Dynamin-Related Proteins superfamily, is a key component in mitochondrial outer membrane fusion in yeast. During this PhD project I built a model of the protein Fzo1p. The structure and dynamics of the model was investigated through molecular modelling and all-atom molecular dynamics simulation in a fully hydrated lipid bilayer environment. The Fzo1p structural model integrates information from several template structures, experimental knowledge, as well as ab initio models of the transmembrane segments. The model is validated experimentally through directed mutagenesis, for instance charge-swap mutations confirm predicted long-distance salt bridges. A series of mutants indicate that coiled-coil domains are required for protein function at variance with its N-terminal region. Overall, the experimental and in silico approaches pinpoint the hinge domains involved in the putative conformational change and identifies critical residues affecting protein stability. Finally, key Fzo1p-GDP interactions provide insights about the molecular mechanism of membrane fusion catalysis. The model provides insight on atomic level and proposes a structure that will be instructional to understanding mitochondrial membrane fusion.

Philosophiae Doctor - PhD / Tsetse flies are the biological vectors of human and animal trypanosomiasis and hence representant medical and veterinary importance. The sense of smell plays a significant role in tsetse and its ecological interaction, such as finding blood meal source, resting, and larvicidal sites and for mating. Tsetse olfactory behaviour can be exploited for their management; however, olfactory studies in tsetse flies are still fragmentary. Here in my PhD thesis, using scanning electron microscopy, electrophysiology, behaviour, bioinformatics and molecular biology techniques, I have investigated tsetse flies (Glossina fuscipes fuscipes) olfaction using behaviourally well studied odorants, tsetse repellent by comparing with attractant odour. Insect olfaction is mediated by olfactory sensory neurons (OSNs), located in olfactory sensilla, which are cuticular structures exposed to the environment through pore and create a platform for chemical communication. In the sensilla shaft the dendrite of OSNs are housed, which are protected by called the sensillum lymph produced by support cells and contains a variety of olfactory proteins, including the odorant binding protein (OBP) and chemosensory proteins (CSP). While on the dendrite of OSNs are expressed olfactory receptors. In my PhD, studies I tried to decipher the sense of smell in tsetse fly. In the second chapter, I demonstrated that G. f. fuscipes is equipped with diverse olfactory sensilla, that various from basiconic, trichoid and coeloconic. I also demonstrated, there is shape, length, number difference between sensilla types and sexual dimorphism. There is a major difference between male and female, while male has the unique basiconic sensilla, club shaped found in the pits, which is absent from female pits. In my third chapter, I investigated the odorant receptors which are expressed on the dendrite of the olfactory sensory neurons (OSNs). G. f. fuscipes has 42 ORs, which were not functionally characterised. I used behaviourally well studied odorants, tsetse repellents, composed of four components blend. I demonstrated that tsetse repellent is also a strong antifeedant for both G. pallidipes and G. f. fuscipes using feeding bioassays as compared to the attractant odour, adding the value of tsetse repellent. However, the attractant odour enhanced the feeding index. Using DREAM (deorphanization of receptors based on expression alterations of mRNA levels). I found that in G. f. fuscipes, following a short in vivo exposure to the individual tsetse repellent component as well as an attractant volatile chemical, OSNs that respond to these compounds altered their mRNA expression in two opposite direction, significant downregulation and upregulation in their number of transcripts corresponding to the OR that they expressed and interacted with odorant. Also, I found that the odorants with opposite valence already segregate distinctly at the cellular and molecular target at the periphery, which is the reception of odorants by OSNs, which is the basis of sophisticated olfactory behaviour. Deorphanization of ORs in none model insect is a challenge, here by combining DREAM with molecular dynamics, as docking score, physiology and homology modelling with Drosophila a well-studied model insects, I was able to predict putative receptors of the tsetse repellent components and an attractant odour. However, many ORs were neutral, showing they were not activated by the odorants, demonstrating the selectivity of the technique as well as the receptors. In my fourth chapter, I investigated the OBPs structures and their interaction with odorants molecules. I demonstrated that OBPs are expressed both in the antenna, as well as in other tissues, such as legs. I also demonstrated that there are variations in the expression of OBPs between tissues as well as sexes. I also demonstrated that odorants induced a fast alteration in OBP mRNA expression, some odorants induced a decrease in the transcription of genes corresponding to the activated OBP and others increased the expression by many fold in OBPs in live insect, others were neutral after 5 hours of exposure. Moreover, with subsequent behavioural data showed that the behavioural response of G. f. fuscipes toward 1-octen-3-ol decreased significantly when 1-octen-3-ol putative OBPs were silenced with feeding of double-stranded RNA (dsRNA). In summary, our finding whereby odorant exposure affects the OBPs mRNA, their physiochemical properties and the silencing of these OBPs affected the behavioural response demonstrate that the OBPs are involved in odour detection that affect the percept of the given odorant. The expression of OBPs in olfactory tissues, antenna and their interaction with odorant and their effect on behavioural response when silenced shows their direct involvement in odour detection and reception. Furthermore, their expression in other tissues such as legs indicates they might also have role in other physiological functions, such as taste.

No description available.

Accurate predictions of binding free energies from computer simulations are an invaluable resource for understanding biochemical processes and drug action. The primary aim of the work described in the thesis was to predict and understand ligand binding to several proteins of major pharmaceutical importance using computational methods. We report a computational strategy to quantitatively predict the effects of alanine scanning and ligand modifications based on molecular dynamics free energy simulations. A smooth stepwise scheme for free energy perturbation calculations is derived and applied to a series of thirteen alanine mutations of the human neuropeptide Y1 G-protein coupled receptor and a series of eight analogous antagonists. The robustness and accuracy of the method enables univocal interpretation of existing mutagenesis and binding data. We show how these calculations can be used to validate structural models and demonstrate their ability to discriminate against suboptimal ones. Site-directed mutagenesis, homology modelling and docking were further used to characterize agonist binding to the human neuropeptide Y2 receptor, which is important in feeding behavior and an obesity drug target.  In a separate project, homology modelling was also used for rationalization of mutagenesis data for an integron integrase involved in antibiotic resistance. Blockade of the hERG potassium channel by various drug-like compounds, potentially causing serious cardiac side effects, is a major problem in drug development. We have used a homology model of hERG to conduct molecular docking experiments with a series of channel blockers, followed by molecular dynamics simulations of the complexes and evaluation of binding free energies with the linear interaction energy method. The calculations are in good agreement with experimental binding affinities and allow for a rationalization of three-dimensional structure-activity relationships with implications for design of new compounds. Docking, scoring, molecular dynamics, and the linear interaction energy method were also used to predict binding modes and affinities for a large set of inhibitors to HIV-1 reverse transcriptase. Good agreement with experiment was found and the work provides a validation of the methodology as a powerful tool in structure-based drug design. It is also easily scalable for higher throughput of compounds.