<p>Human Immunodeficiency Virus (HIV) is the causative agent of the pandemic disease Acquired Immune Deficiency Syndrome (AIDS). HIV acts to disrupt the immune system which makes the body susceptible to opportunistic infections. Untreated, AIDS is generally fatal. Twenty years of research by countless scientists around the world has led to the discovery and exploitation of several targets in the replication cycle of HIV. Many lives have been saved, prolonged and improved as a result of this massive effort. One particularly successful target has been the inhibition of HIV protease. In combination with the inhibition of HIV reverse transcriptase, protease inhibitors have helped to reduce viral loads and partially restore the immune system. Unfortunately, viral mutations leading to drug resistance and harmful side-effects of the current medicines have identified the need for new drugs to combat HIV.</p><p>This study presents computational efforts to understand the interaction of inhibitors to HIV protease. The first part of this study has used molecular modelling and Comparative Molecular Field Analysis (CoMFA) to help explain the structure-active relationship of a novel series of protease inhibitors. The inhibitors are sulfamide derivatives structurally similar to the cyclic urea candidate drug mozenavir (DMP-450). The central ring of the sulfamides twists to adopt a nonsymmetrical binding mode distinct from that of the cyclic ureas. The energetics of this twist has been studied with <i>ab initio</i> calculations to develop improved empirical force field parameters for use in molecular modelling.</p><p>The second part of this study has focused on an analysis of the association and dissociation kinetics of a broad collection of HIV protease inhibitors. Quantitative models have been derived using CoMFA which relate the dissociation rate back to the chemical structures. Efforts have also been made to improve the models by systematically varying the parameters used to generate them.</p>

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

Novel and potent ligands to the serotonin7 (5-HT7) receptor have been synthesized. The synthesized compounds include a set of substituted pyrimidines which show high affinity to the 5-HT7 receptor, synthesized by previously described methods [1,2] in high yield. Comparing the affinities of substituted pyrimidines to previously calculated models [3,4] yielded new hypotheses about the nature of interaction between the pyrimidine ligands and the 5-HT7 binding site. Several new series of compounds were synthesized by various methods to validate these hypotheses, including a conjugate addition to vinylpyrimidines [5]. These compounds include benzofurans, oximes, hydrazones, as well as a group of substituted piperazines. All series of compounds show affinity to the 5-HT7 receptor comparable to previously synthesized 5-HT7 ligands. Several of the synthesized ligands show affinity which exceeds that of currently available ligands. The synthesized compounds were evaluated quantitatively by calculating a three-dimensional quantitative structure-affinity relationship (3D-QSAR) for the 5-HT7 receptor. Evaluation of the calculated model validated qualitative assumptions about the data set as well as described regions of interaction in greater detail than previously available. These observations give further insight on the nature of ligand-binding site interactions with highly potent ligands such as 4-(3-furyl)-2-(N-methylpiperazino)pyrimidine which will lead to more potent 5-HT7 receptor ligands. Additionally, a model was calculated for affinity to the 5-HT2a receptor. Comparing this model to that calculated for affinity to the 5-HT7 receptor identified two regions which may be exploited in future sets of ligands to increase selectivity to the 5HT7 receptor.

Human Immunodeficiency Virus (HIV) is the causative agent of the pandemic disease Acquired Immune Deficiency Syndrome (AIDS). HIV acts to disrupt the immune system which makes the body susceptible to opportunistic infections. Untreated, AIDS is generally fatal. Twenty years of research by countless scientists around the world has led to the discovery and exploitation of several targets in the replication cycle of HIV. Many lives have been saved, prolonged and improved as a result of this massive effort. One particularly successful target has been the inhibition of HIV protease. In combination with the inhibition of HIV reverse transcriptase, protease inhibitors have helped to reduce viral loads and partially restore the immune system. Unfortunately, viral mutations leading to drug resistance and harmful side-effects of the current medicines have identified the need for new drugs to combat HIV. This study presents computational efforts to understand the interaction of inhibitors to HIV protease. The first part of this study has used molecular modelling and Comparative Molecular Field Analysis (CoMFA) to help explain the structure-active relationship of a novel series of protease inhibitors. The inhibitors are sulfamide derivatives structurally similar to the cyclic urea candidate drug mozenavir (DMP-450). The central ring of the sulfamides twists to adopt a nonsymmetrical binding mode distinct from that of the cyclic ureas. The energetics of this twist has been studied with ab initio calculations to develop improved empirical force field parameters for use in molecular modelling. The second part of this study has focused on an analysis of the association and dissociation kinetics of a broad collection of HIV protease inhibitors. Quantitative models have been derived using CoMFA which relate the dissociation rate back to the chemical structures. Efforts have also been made to improve the models by systematically varying the parameters used to generate them.

Os inibidores seletivos da recaptação de serotonina (ISRSs) representam a classe mais importante de antidepressivos em uso clínico atualmente. Entretanto, esses medicamentos costumam levar de duas a seis semanas para apresentar os efeitos de sua ação terapêutica. Estudos clínicos mostram que quando um antagonista do receptor 5-HT1A é administrado juntamente com um ISRS, um aumento da concentração extracelular de serotonina é observado nas áreas terminais dos neurônios. Assim, a combinação de um antagonista do receptor 5-HT1A com um ISRS pode acelerar o início da ação antidepressiva, aumentando a eficácia do tratamento farmacológico da depressão. A classe mais importante de ligantes do receptor 5-HT1A são os compostos arilpiperazínicos. O presente estudo teve como objetivo o entendimento das características importantes para as interações entre uma série de compostos arilpiperazínicos com o receptor 5-HT1A. Para tal, foram realizados estudos de Relação Quantitativa entre Estrutura e Atividade (QSAR) bi- e tridimensionais, empregando as seguintes abordagens: métodos quimiométricos baseados em descritores teóricos, QSAR por hologramas (HQSAR) e o método de Análise Comparativa de Campos Moleculares (CoMFA). Essas análises foram complementadas com a modelagem por homologia do receptor 5-HT1A e com estudos de docking ligante-receptor realizados para alguns compostos arilpiperazínicos. Modelos de QSAR com boa consistência interna, habilidade preditiva e estabilidade foram obtidos em todos os casos. Os modos de interação observados apresentaram consistência com dados experimentais disponíveis sobre os resíduos importantes para as interações com ligantes arilpiperazínicos. Os principais resultados indicaram algumas características dos ligantes que são importantes para a afinidade pelo receptor 5-HT1A, tais como a presença de um anel benzotiofeno como substituinte Ar2, substituintes pouco volumosos na posição Z e receptores de ligações de hidrogênio na posição orto do anel Ar1. Esses resultados foram corroborados pelo estudo das interações com o modelo do receptor 5-HT1A, que indicou uma importante interação hidrofóbica do grupo benzotiofeno com o resíduo Trp6.48 do receptor, assim como uma ligação de hidrogênio entre a hidroxila na posição Z e o resíduo Thr3.37 e, ainda, entre o oxigênio do anel Ar1 e o resíduo Asn7.39. As informações obtidas neste estudo podem fornecer subsídios para o planejamento de novos ligantes com afinidade pelo receptor 5-HT1A. / Selective serotonin reuptake inhibitors (SSRIs) are the most important class of antidepressants in current clinical use. However, they present the serious drawback of a delay of two to six weeks in the onset of therapeutic effect. Clinical studies have shown that when a 5-HT1A receptor antagonist is administrated along with a SSRI, an increase of extracellular serotonin concentration in neuronal terminal areas is observed. Thus, the combination of a 5- HT1A receptor antagonist and a SSRI could accelerate the onset of antidepressant action, improving the pharmacological treatment of depression. The most important class of 5-HT1A receptor ligands are arylpiperazine compounds. In the present study, our aim was to understand the main features of the interaction between a series of arylpiperazines and the 5- HT1A receptor. Bi- and Tridimensional Quantitative Structure-Activity Relationship (QSAR) studies were conducted employing the following approaches: chemometric methods based on theoretical descriptors, Hologram QSAR (HQSAR), and Comparative Molecular Field Analysis (CoMFA). These analyses were complemented by 5-HT1A receptor homology modeling and ligand-receptor docking studies. QSAR models presenting good internal consistency, predictive power and stability were obtained in all cases. The observed binding modes are consistent with available experimental data on residues considered crucial for interactions with arylpiperazine compounds. The main results have indicated some important features for optimal binding to the 5-HT1A receptor, such as the presence of a benzothiophene ring as Ar2 substituent, small groups at position Z and hydrogen bond acceptors at the ortho position of Ar1 ring. These results were corroborated by modeling the interactions with the 5- HT1A receptor, which has indicated an important hydrophobic interaction between the benzothiophene group and residue Trp6.48, a hydrogen bond between the OH group at position Z and residue Thr3.37, as well as between the oxygen in Ar1 and residue Asn7.39. The information gathered in these studies can be useful for the design of new ligands displaying affinity to the 5-HT1A receptor.

Computational drug design methods have great potential in drug discovery particularly in lead identification and lead optimization. 3-Phosphoinositide dependent kinase-1 (PDK1) is a protein kinase and a well validated anti-cancer target. Inhibitors of PDK1 have the potential to be developed as anti-cancer drugs. In this work, we have applied various novel computational drug design strategies to design and identify new PDK1 inhibitors with potential anti-cancer activity. We have pursued novel structure-based drug design strategies and identified a new binding mode for celecoxib and its derivatives binding with PDK1. This new binding mode provides a valuable basis for rational design of potent PDK1 inhibitors. In order to understand the structure-activity relationship of indolinone-based PDK1 inhibitors, we have carried out a combined molecular docking and three-dimensional quantitative structure-activity relationship (3D-QSAR) modeling study. The predictive ability of the developed 3D-QSAR models were validated using an external test set of compounds. An efficient strategy of the hierarchical virtual screening with increasing complexity was pursued to identify new hits against PDK1. Our approach uses a combination of ligand-based and structure-based virtual screening including shape-based filtering, rigid docking, and flexible docking. In addition, a more sophisticated molecular dynamics/molecular mechanics- Poisson-Boltzmann surface area (MD/MM-PBSA) analysis was used as the final filter in the virtual screening. Our screening strategy has led to the identification of a new PDK1 inhibitor. The anticancer activities of this compound have been confirmed by the anticancer activity assays of national cancer institute-developmental therapeutics program (NCI-DTP) using 60 cancer cell lines. The PDK1-inhibitor binding mode determined in this study may be valuable in future de novo drug design. The virtual screening approach tested and used in this study could also be applied to lead identification in other drug discovery efforts.

A reemergência de algumas bactérias Gram-positivas, em particular, do gênero Staphylococcus, como principal foco causador de infecções hospitalares, tem se intensificado nas últimas décadas, e, apesar da existência de potentes fármacos voltados para o tratamento de infecções causadas por este gênero de bactéria, as taxas de morbidade e mortalidade prevalecem com perfil crescente. Além disso, um grande problema associado a cepas de MRSA (Methicillin-Resistant Staphylococcus aureus) é o fenótipo de multi-resistência, característica que confere a este microrganismo resistência não apenas à meticilina como também a uma série de outros fármacos, exceto frente à vancomicina e à teicoplanina. Muito tem se feito, mas ainda são poucos os resultados efetivamente aplicáveis no tratamento de infecções com caráter de multi-resistência, justificando, desta forma, a necessidade de desenvolvimento de sucedâneos que sejam consideravelmente mais efetivos para a solução deste problema. Baseado nestes fatos, a proposta deste estudo envolveu o planejamento, síntese, identificação e estudos de QSAR (Quantitative Structure-Activity Relationships) em duas e três dimensões de derivados 5-nitro-2-tiofilidênicos com atividade antimicrobiana frente a cepas padrão e multi-resistente de Staphylococcus aureus. A escolha dos grupos substituintes foi realizada em duas etapas. Na primeira delas seguiu-se metodologia de substituição em anéis aromáticos proposta por Topliss para a otimização da bioatividade de compostos. Em uma segunda etapa, predominantemente quantitativa, foram selecionados mais alguns derivados baseando-se em faixa de hidrofobicidade ótima pré-determinada experimentalmente e na variação de efeito estérico dos grupos substituintes. Quatorze derivados 5-nitro-2-tiofilidênicos foram sintetizados, estruturalmente identificados e avaliados quanto à atividade antimicrobiana frente às cepas padrão (ATCC 25923) e multi-resistente (3SP/R33) de Staphylococcus aureus por determinação da concentração inibitória mínima empregando-se método de macrodiluição sucessiva em tubos. Salienta-se que a cepa 3SP/R33 se mostra resistente a dezenove antibióticos empregados na prática médica e apresenta suscetibilidade apenas à vancomicina. As concentrações inibitória e bactericida mínimas apresentadas pelos compostos sintetizados mostraram sofrer influência significativa da hidrofobicidade sobre as referidas atividades de acordo com os estudos de QSAR-2D e QSAR-3D, sendo os resultados obtidos para a cepa multi-resistente absolutamente compatíveis com os anteriormente determinados para a cepa padrão. Os estudos de QSAR-2D indicaram que a atividade antimicrobiana das 5nitro-2-tiofilideno benzidrazidas substituídas sofre influência significativa de duas propriedades físico-químicas que são a hidrofobicidade e a distribuição eletrônica. A relevância dos descritores estruturais &#963; e efe na determinação da atividade antimicrobiana, sinalizam que a distribuição eletrônica influencia fortemente o aumento da potência antimicrobiana dos compostos em estudo tanto pela influência dos efeitos indutivo e de ressonância na estrutura química do ligante, como também pelos campos moleculares gerados ao redor de grupos substituintes, sugerindo uma possível interação dos mesmos com uma área específica do sítio receptor. Nos estudos de QSAR-3D, foi evidenciado, em concordância com o estudo clássico anteriormente realizado, que a hidrofobicidade prevalece como propriedade de fundamental importância no estabelecimento da atividade antimicrobiana. Foi observada a importância da presença de regiões hidrofílicas pontuais nos compostos de forma a propiciar processos de solvatação e dessolvatação que são críticos na difusão através de membranas biológicas. Pode-se afirmar que a análise de QSAR, considerando os aspectos tridimensionais ligantes, ressaltou a necessidade de um balanço lipofílico-hidrofílico para um bom desempenho das 5-nitro-2-tiofilideno benzidrazidas &#961;-substituídas como agentes antimicrobianos. A partir dos resultados obtidos evidenciou-se, neste estudo, o forte potencial de derivados 5-nitro-2-tiofilidênicos como possível alternativa para o desenvolvimento racional, em nível molecular, de fármacos voltados para o tratamento de infecções causadas por cepas multi-resistentes de Staphylococcus aureus. / In the last decade, there has been a reemergence of Gram-positive bacteria, in particular Staphylococcus, which isconsidered one of the. most causing of nosocomial infections. Although potent antistaphylococcal drugs are available, this infection continues presenting increasing morbidity and mortality rates. Besides, a serious problem associated with MRSA (Methicillin-Resistant Staphylococcus aureus) is the phenotype of multidrug resistance, which is, resistance not only to methicillin but also to many other drugs, except to vancomycin and teicoplanin. Many efforts have been made in a tentative to reduce this problem, nevertheless there is only a few number of alternatives to combat Staphylococcus aureus multidrug-resistant strains, justifying the necessity of development of more effective compounds to the treatment of these infections. Based in these facts, the purpose of this study was the design, synthesis, structural identification and 2D-QSAR and 3D-QSAR (Quantitative Structure-Activity Relationships) studies of 5-nitro-2-thiophylidene derivatives with antimicrobial activity against multidrug-resistant strains of Staphylococcus aureus. The choice of substituent groups was made in two stages. The first stage comprises on application of Topliss operational scheme for aromatic substitution. In a second quantitative stage, more derivatives were selected according by hydrophobicity range previously determined. Other standard considered at the selection of substituent groups was the variation of steric effect. Fourteen 5-nitro-2-thiophylidene derivatives were synthesized, structural identified and tested against standard (A TCC 25923) and multidrug-resistant (3SP/R33) strains of Staphylococcus aureus. The Minimal Inhibitory Concentration, MIC, was determined using the serial dilution tests in two sequential stages. The 3SP/R33 strain is resistant to nineteen antimicrobial agents in use, except to vancomycin. The minimal inhibitory and bactericidal concentrations of synthesized compounds showed, according by 2D-QSAR and 3D-QSAR studies, a significant influence of hydrophobic properties on antimicrobial activity determination and the results obtained for multidrug-resistant strain were consistent with those determined for A TCC 25923 strain. 2D-QSAR studies showed that antimicrobial activity are mainly influenced by two physico-chemical properties: hydrophobicity and electronic distribution. The relevance of &#963; e ephe parameters on antimicrobial activity determination, denotes the contribution of inductive and resonance effects for the polar performed by the substituent groups, probably suggesting an interaction between them and specific receptor site. 3D-QSAR studies showed that hydrophobicity is a essential property to antimicrobial activity determination, sustained the same conclusions previously obtained by Hansch Analysis. It was observed a great concern of small hydrophilic regions distributed on derivatives in order to promote solvation and desolvation process, that have critical importance on diffusion process through the biological membranes. QSAR studies considering three-dimensional properties of ligands indicated the necessity of accurate hydrophilic-hydrophobic balance on nitrothiophene derivatives for their good performance as antimicrobial agents. The results obtained in this preliminary study have shown the potential of synthesized compounds as alternatives to the treatment of infections caused by multidrug-resistant strains of Staphylococcus aureus.

Con el fin de relacionar la estructura y la actividad de series de compuestos, es importante usar descriptores moleculares relevantes. Los descriptores GRIND y VolSurf pertenecen a una nueva familia de descriptores llamado libre de alineamiento. Es decir, que no necesitan alinear los compuestos con el fin de comparar sus campos de interacciones molecular. En este estudio se ha aplicado esos descriptores para la selección de reactivos químicos a partir de una amplia base de datos. La selección se ha echo mediante un protocolo que permite maximizar la diversidad de la muestra y así obtener unos compuestos muy informativos. También se ha desarrollado nuevos descriptores de forma que están basado en los cambios de curvatura de la superficie molecular. Los resultados obtenidos indican que los nuevos descriptores de forma se integran muy bien en los descriptores GRIND originales y que permiten identificar los efectos de forma tanto favorable como desfavorable. Además, se ha desarrollado nuevos descriptores libre de alineamiento llamado 'anchor-GRIND' que usan un átomo de cada molécula como punto de referencia para la comparación de los campos de interacciones molecular. Los descriptores 'anchor-GRIND' permiten una descripción mas precisa y mas sencilla que los descriptores GRIND lo que los hace mas relevante para el análisis de ciertas familias de compuestos. / In order to correlate the differences of structure with the differences of activity of series of compounds, it is important to use relevant molecular descriptors. The GRIND and VolSurf descriptors belong to the so-called alignment-free descriptors family. In other words, they do not require to align the compounds in order to compare its molecular interaction fields. In this study, we applied these descriptors to the selection of chemical reagent from a database of compounds. The selection has been done following a protocol which allows to maximize the diversity of the sample and so to obtain some compounds highly informative. In addition we developed new shape descriptors which are based on the changes of curvature of the molecular surface. The results obtained show that the new shape descriptors are well integrated in the original GRIND descriptors. Furthermore, we designed new alignment-free descriptors called 'anchor-GRIND' which use one atom of each molecule as a reference point for the comparison of the molecular interaction fields. The 'anchor-GRIND' descriptors allow a more precise and more simple description than the GRIND descriptors, which makes them more relevant for the analysis of some families of compounds.