RESUMO
Predicting the three-dimensional structure of ligand-receptor complexes involving G protein-coupled receptors (GPCRs) is still a challenging task in rational drug design. To evaluate the reliability of the GPCR structural prediction, only a couple of community-wide assessments have been carried out. Our participation in the last edition, DOCK2010, involved the blind prediction of the dopaminergic D(3) receptor in complex with the D(2)/D(3) selective antagonist eticlopride for which the crystal structure has been recently released. Here, we describe a methodology that succeeded to produce a correctly predicted eticlopride-D(3) receptor complex out of three submitted models. Ranking the obtained models in the correct order is the main challenge due to subtle structural differences in the complex that are not sufficiently captured by conventional scoring functions. Importantly, our work reveals that a correct ranking is obtained by including a more sophisticated description of conformational ligand energy on binding. All in all, this case study highlights the current progress in modeling GPCR complexes and underlines that in silico modeling can be a valuable complement in GPCR drug discovery.
Assuntos
Antagonistas de Dopamina/farmacologia , Receptores de Dopamina D3/antagonistas & inibidores , Receptores de Dopamina D3/química , Salicilamidas/farmacologia , Sítios de Ligação , Antagonistas de Dopamina/química , Humanos , Modelos Moleculares , Ligação Proteica , Conformação Proteica , Receptores de Dopamina D3/metabolismo , Salicilamidas/química , Homologia Estrutural de ProteínaRESUMO
The methylerythritol 4-phosphate (MEP) pathway for the biosynthesis of the isoprenoid universal building blocks (isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP)) is present in most of human pathogens and is absent in animals, turning it into a promising therapeutic druggable pathway. Two different strategies, a pharmacophore-directed virtual screening and a protein-protein interaction (PPI)-mimicking cyclic peptide were used to search for compounds that bind to the PPI surface of the 4-(cytidine 5-diphospho)-2C-methyl-D-erythritol kinase (CMK), which catalyzes the fourth step of the MEP pathway. A significant part of the pharmacophore hypothesis used in this study was designed by mimicking water-mediated PPI relevant in the CMK homodimer complex stabilization. After database search and with the aid of docking and molecular dynamics (MD) simulations, a 7H-furo[3,2-g]chromen-7-one derivative and a cyclic peptide were chosen as candidates to be ligands of CMK. Their binding affinities were measured using surface plasmon resonance (SPR) technology.
Assuntos
Desenho de Fármacos , Proteínas de Escherichia coli/metabolismo , Ligantes , Mimetismo Molecular , Fosfotransferases (Aceptor do Grupo Álcool)/metabolismo , Proteínas de Escherichia coli/genética , Conformação Molecular , Simulação de Dinâmica Molecular , Peptídeos Cíclicos/síntese química , Peptídeos Cíclicos/metabolismo , Fosfotransferases (Aceptor do Grupo Álcool)/genética , Ligação Proteica , Solventes/química , Ressonância de Plasmônio de SuperfícieRESUMO
The preclinical assessment of drug-induced ventricular arrhythmia, a major concern for regulators, is typically based on experimental or computational models focused on the potassium channel hERG (human ether-a-go-go-related gene, K(v)11.1). Even if the role of this ion channel in the ventricular repolarization is of critical importance, the complexity of the events involved make the cardiac safety assessment based only on hERG has a high risk of producing either false positive or negative results. We introduce a multiscale simulation system aiming to produce a better cardiotoxicity assessment. At the molecular scale, the proposed system uses a combination of docking simulations on two potassium channels, hERG and KCNQ1, plus three-dimensional quantitative structure-activity relationship modeling for predicting how the tested compound will block the potassium currents IK(r) and IK(s). The obtained results have been introduced in electrophysiological models of the cardiomyocytes and the ventricular tissue, allowing the direct prediction of the drug effects on electrocardiogram simulations. The usefulness of the whole method is illustrated by predicting the cardiotoxic effect of several compounds, including some examples in which classic hERG-based models produce false positive or negative results, yielding correct predictions for all of them. These results can be considered a proof of concept, suggesting that multiscale prediction systems can be suitable for being used for preliminary screening in lead discovery, before the compound is physically available, or in early preclinical development when they can be fed with experimentally obtained data.
Assuntos
Biologia Computacional/métodos , Efeitos Colaterais e Reações Adversas Relacionados a Medicamentos , Coração/efeitos dos fármacos , Canal de Potássio ERG1 , Fenômenos Eletrofisiológicos/efeitos dos fármacos , Canais de Potássio Éter-A-Go-Go/antagonistas & inibidores , Canais de Potássio Éter-A-Go-Go/química , Coração/fisiologia , Humanos , Canal de Potássio KCNQ1/antagonistas & inibidores , Canal de Potássio KCNQ1/química , Modelos Moleculares , Bloqueadores dos Canais de Potássio/efeitos adversos , Bloqueadores dos Canais de Potássio/química , Bloqueadores dos Canais de Potássio/farmacologia , Conformação Proteica , Relação Quantitativa Estrutura-Atividade , Reprodutibilidade dos TestesRESUMO
Glucose-6-phosphate dehydrogenase (G6PDH) is an essential enzyme involved in the first reaction of the oxidative branch of the pentose phosphate pathway (PPP). Recently, G6PDH was suggested as a novel target protein for cancer therapy as one of the final products of the PPP, ribose-5-phosphate, is necessary for nucleic acid synthesis and tumor progression. After analyzing the protein-protein interface of the crystal structure of human G6PDH by means of molecular dynamics simulations, we designed six interface peptides based on the natural sequence of the protein. The three most promising peptides, as predicted by binding free energy calculations, were synthesized and one of them was confirmed as a novel inhibitor of human G6PDH in experimental assays. Together, the active peptide found and its suggested binding mode proposes a new strategy for inhibiting this enzyme and should aid the further design of novel, potent and non-peptidic G6PDH inhibitors.
Assuntos
Inibidores Enzimáticos/química , Inibidores Enzimáticos/farmacologia , Glucosefosfato Desidrogenase/antagonistas & inibidores , Glucosefosfato Desidrogenase/metabolismo , Peptídeos/química , Peptídeos/farmacologia , Humanos , Simulação de Dinâmica Molecular , Estrutura Secundária de Proteína , Estrutura Terciária de ProteínaRESUMO
Survivin, the smallest inhibitor of apoptosis protein (IAP), is a valid target for cancer research. It mediates both the apoptosis pathway and the cell cycle and has been proposed to form a complex with the cyclin-dependent kinase protein CDK4. The resulting complex transports CDK4 from the cytosol to the nucleus, where CDK4 participates in cell division. Survivin has been recognized as a node protein that interacts with several partners; disruption of the formed complexes can lead to new anticancer compounds. We propose a rational model of the survivin/CDK4 complex that fulfills the experimental evidence and that can be used for structure-based design of inhibitors modifying its interface recognition. In particular, the suggested complex involves the alpha helical domain of survivin and resembles the mode of binding of survivin in the survivin/borealin X-ray structure. The proposed model has been obtained by combining protein-protein docking, fractal-based shape complementarity, electrostatics studies and extensive molecular dynamics simulations.
Assuntos
Quinase 4 Dependente de Ciclina/química , Proteínas Inibidoras de Apoptose/química , Proteínas de Neoplasias/química , Sítios de Ligação , Proteínas de Ciclo Celular/química , Cristalografia por Raios X , Quinase 6 Dependente de Ciclina/química , Inibidor p16 de Quinase Dependente de Ciclina/química , Humanos , Simulação de Acoplamento Molecular , Simulação de Dinâmica Molecular , Ligação Proteica , Estrutura Secundária de Proteína , Estrutura Terciária de Proteína , Eletricidade Estática , Homologia Estrutural de Proteína , SurvivinaRESUMO
Transketolase is an enzyme involved in a critical step of the non-oxidative branch of the pentose phosphate pathway whose inhibition could lead to new anticancer drugs. Here, we report new human transketolase inhibitors, based on the phenyl urea scaffold, found by applying structure-based virtual screening. These inhibitors are designed to cover a hot spot in the dimerization interface of the homodimer of the enzyme, providing for the first time compounds with a suggested novel binding mode not based on mimicking the thiamine pyrophosphate cofactor.
Assuntos
Carbanilidas/química , Carbanilidas/farmacologia , Transcetolase/antagonistas & inibidores , Interface Usuário-Computador , Antineoplásicos/química , Antineoplásicos/metabolismo , Antineoplásicos/farmacologia , Carbanilidas/metabolismo , Sobrevivência Celular/efeitos dos fármacos , Cristalografia por Raios X , Avaliação Pré-Clínica de Medicamentos , Inibidores Enzimáticos/química , Inibidores Enzimáticos/metabolismo , Inibidores Enzimáticos/farmacologia , Células HCT116 , Células HT29 , Humanos , Modelos Moleculares , Estrutura Secundária de Proteína , Relação Estrutura-Atividade , Termodinâmica , Transcetolase/química , Transcetolase/metabolismoRESUMO
Protein surface roughness is a structural property associated with ligand-protein and protein-protein binding interfaces. In this work we apply for the first time the concept of surface roughness, expressed as the fractal dimension, to address structure and function of G protein-coupled receptors (GPCRs) which are an important group of drug targets. We calculate the exposure ratio and the fractal dimension for helix-forming residues of the ß(2) adrenergic receptor (ß(2)AR), a model system in GPCR studies, in different conformational states: in complex with agonist, antagonist and partial inverse agonists. We show that both exposure ratio and roughness exhibit periodicity which results from the helical structure of GPCRs. The pattern of roughness and exposure ratio of a protein patch depends on its environment: the residues most exposed to membrane are in general most rough whereas parts of receptors mediating interhelical contacts in a monomer or protein complex are much smoother. We also find that intracellular ends (TM3, TM5, TM6 and TM7) which are relevant for G protein binding and thus receptor signaling, are exposed but smooth. Mapping the values of residual fractal dimension onto receptor 3D structures makes it possible to conclude that the binding sites of orthosteric ligands as well as of cholesterol are characterized with significantly higher roughness than the average for the whole protein. In summary, our study suggests that identification of specific patterns of roughness could be a novel approach to spot possible binding sites which could serve as original drug targets for GPCRs modulation.
Assuntos
Fractais , Receptores Adrenérgicos beta 2/química , Receptores Adrenérgicos beta 2/metabolismo , Sítios de Ligação , Colesterol/metabolismo , Modelos Moleculares , Conformação Proteica , Relação Estrutura-Atividade , Propriedades de SuperfícieRESUMO
Tuberculosis is one of the leading infectious diseases in humans. Discovering new treatments for this disease is urgently required, especially in view of the emergence of multiple drug resistant organisms and to reduce the total duration of current treatments. The synthesis of isoprenoids in Mycobacterium tuberculosis has been reported as an interesting pathway to target, and particular attention has been focused on the methylerythritol phosphate (MEP) pathway comprising the early steps of isoprenoid biosynthesis. In this context we have studied the enzyme 2C-methyl-D-erythritol-4-phosphate cytidylyltransferase (CMS), the third enzyme in the MEP pathway, since the lack of a resolved structure of this protein in M. tuberculosis has seriously limited its use as a drug target. We performed homology modeling of M. tuberculosis CMS in order to provide a reliable model for use in structure-based drug design. After evaluating the quality of the model, we performed a thorough study of the catalytic site and the dimerization interface of the model, which suggested the most important sites (conserved and non-conserved) that could be useful for drug discovery and mutagenesis studies. We found that the metal coordination of CDP-methylerythritol in M. tuberculosis CMS differs substantially with respect to the Escherichia coli variant, consistent with the fact that the former is able to utilize several metal ions for catalysis. Moreover, we propose that electrostatic interactions could explain the higher affinity of the MEP substrate compared with the cytosine 5'-triphosphate substrate in the M. tuberculosis enzyme as reported previously.
Assuntos
Proteínas de Bactérias/química , Modelos Moleculares , Mycobacterium tuberculosis/enzimologia , Nucleotidiltransferases/química , Terpenos/metabolismo , Sequência de Aminoácidos , Proteínas de Bactérias/genética , Proteínas de Bactérias/metabolismo , Sítios de Ligação/genética , Vias Biossintéticas , Cristalografia por Raios X , Cistina Difosfato/química , Cistina Difosfato/metabolismo , Eritritol/análogos & derivados , Eritritol/biossíntese , Eritritol/química , Ligação de Hidrogênio , Metais/química , Metais/metabolismo , Dados de Sequência Molecular , Estrutura Molecular , Mycobacterium tuberculosis/genética , Nucleotidiltransferases/genética , Nucleotidiltransferases/metabolismo , Ligação Proteica , Multimerização Proteica , Estrutura Terciária de Proteína , Homologia de Sequência de Aminoácidos , Fosfatos Açúcares/biossíntese , Fosfatos Açúcares/química , TermodinâmicaRESUMO
Transketolase, the most critical enzyme of the non-oxidative branch of the pentose phosphate pathway, has been reported as a new target protein for cancer research. However, since the crystal structure of human Transketolase is unknown, no structure-based methods can be used to identify new inhibitors. We performed homology modeling of human Transketolase using the crystal structure of yeast as a template, and then refined the model through molecular dynamics simulations. Based on the resulting structure we propose five critical sites containing arginines (Arg 101, Arg 318, Arg 395, Arg 401 and Arg 474) that contribute to dimer stability or catalytic activity. In addition, an interaction analysis of its cofactor (thiamine pyrophosphate) and a binding site description were carried out, suggesting the substrate channel already identified in yeast Transketolase. A binding free energy calculation of its cofactor was performed to establish the main driving forces of binding. In summary, we describe a reliable model of human Transketolase that can be used in structure-based drug design and in the search for new Transketolase inhibitors that disrupt dimer stability and cover the critical sites found.
Assuntos
Desenho de Fármacos , Transcetolase/química , Transcetolase/metabolismo , Sequência de Aminoácidos , Animais , Sítios de Ligação , Domínio Catalítico , Simulação por Computador , Sequência Conservada , Humanos , Modelos Moleculares , Dados de Sequência Molecular , Ligação Proteica , Multimerização Proteica , Estrutura Quaternária de Proteína , Ratos , Saccharomyces cerevisiae/enzimologia , Alinhamento de Sequência , Eletricidade Estática , Homologia Estrutural de ProteínaRESUMO
Apoptosis, also called programmed cell death, is a conserved mechanism inherent to all cells that sentences them to death when they receive the appropriate external stimuli. Inhibitor of apoptosis proteins (IAPs) are a family of regulatory proteins that suppress such cell death. XIAP is the most commonly studied member of the IAP family. It binds to and inhibits Caspases, an important family of apoptotic proteases. In addition, XIAP over-expression has been detected in numerous types of cancer. Smac/DIABLO, a mitochondrial protein that binds to IAPs and promotes Caspase activation, has the opposite action to XIAP and can be considered a key protein in the regulation of IAPs. Survivin, the smallest IAP protein, has received a lot of attention due to its specific expression in many cancer cell lines. It has been shown to interact with Smac/DIABLO, even though the structure of this complex has not yet been reported.We analysed the protein-protein interactions appearing in the Smac/DIABLO-XIAP and Smac/DIABLO-Survivin complexes fully, using molecular dynamics simulations. This information is a first step towards the design of Smac/DIABLO peptidomimetics that could be used as innovative therapeutic agents for the treatment of malignancy. Our results complement the experimental interactions described for the first complex and provide a detailed description for the second. We show that Smac/DIABLO interacts in a similar way with both targets through its amino terminal residues. In addition, we identify a pharmacophore formed by eight stable protein-protein interactions for the XIAP complex and seven stable protein-protein interactions for the Survivin complex, which describe the whole contact surface. This information is used to suggest the binding mode of embelin, the first non-peptidic inhibitor of XIAP, and two of its derivatives. Molecular docking and molecular dynamics simulations were also carried out to describe ligand and receptor flexibility. Finally, an MMGBSA protocol was used to obtain a more quantitative description of the binding in all the complexes studied.
Assuntos
Proteínas Associadas aos Microtúbulos/antagonistas & inibidores , Proteínas de Neoplasias/antagonistas & inibidores , Proteínas Inibidoras de Apoptose Ligadas ao Cromossomo X/antagonistas & inibidores , Sequência de Aminoácidos , Proteínas Reguladoras de Apoptose , Benzoquinonas/química , Simulação por Computador , Humanos , Ligação de Hidrogênio , Proteínas Inibidoras de Apoptose , Peptídeos e Proteínas de Sinalização Intracelular/metabolismo , Ligantes , Proteínas Associadas aos Microtúbulos/química , Proteínas Mitocondriais/metabolismo , Modelos Moleculares , Dados de Sequência Molecular , Proteínas de Neoplasias/química , Ligação Proteica , Eletricidade Estática , Survivina , Termodinâmica , Proteínas Inibidoras de Apoptose Ligadas ao Cromossomo X/química , Zinco/metabolismoRESUMO
Evaluation of binding free energy in receptor-ligand complexes is one of the most important challenges in theoretical drug design. Free energy is directly correlated to the thermodynamic affinity constant, and, as a first step in druglikeness, a lead compound must have this constant in the range of micro- to nanomolar activity. Many efforts have been made to calculate it by rigorous computational approaches, such as free energy perturbation or linear response approximation. However, these methods are still computationally expensive. We focus our work on XIAP, an antiapoptotic protein whose inhibition can lead to new drugs against cancer disease. We report here a comparative evaluation of two completely different methodologies to estimate binding free energy, MMPBSA (a force field based function) and XSCORE (an empirical scoring function), in seven XIAP-peptide complexes using a representative set of structures generated by previous molecular dynamics simulations. Both methods are able to predict the experimental binding free energy with acceptable errors, but if one needs to identify slight differences upon binding, MMPBSA performs better, although XSCORE is not a bad choice taking into account the low computational cost of this method.