Crystal Engineering of an Anti-HIV Drug Based on the Recognition of Assembling Molecular Frameworks

A rational strategy was employed for design of an orthorhombic structure of lamivudine with maleic acid. On the basis of the lamivudine saccharinate structure reported in the literature, maleic acid was chosen to synthesize a salt with the anti-HIV drug because of the structural similarities between the salt formers. Maleic acid has an acid-ionization constant of the anti first proton and an arrangement of their hydrogen bonding functionalities similar to those of saccharin. Likewise, there is a saccharin-like conformational rigidity in maleic acid because of the hydrogen-bonded ring formation and the Z-configuration around the C=C double bond. As was conceivably predicted, lamivudine maleate assembles into a structure whose intermolecular architecture is related to that of saccharinate salt of the drug. Therefore, a molecular framework responsible for crystal assembly into a lamivudine saccharinate-like structure could be recognized in the salt formers. Furthermore, structural correlations and structure-solubility relationships were established for lamivudine maleate and saccharinate. Although there is a same molecular framework in maleic acid and saccharin, these salt formers are Structurally different in some aspects. When compared to saccharin, neither out-of-plane SO(2) oxygens nor a benzene group occur in maleic acid. Both features could be related to higher solubility of lamivudine maleate. Here, we also anticipate that multicomponent molecular crystals of lamivudine with other salt formers possessing the molecular framework responsible for crystal assembly can be engineered successfully.

Quantitative structure-activity relationships for a series of inhibitors of cruzain from Trypanosoma cruzi: Molecular modeling, CoMFA and CoMSIA studies

Human parasitic diseases are the foremost threat to human health and welfare around the world. Trypanosomiasis is a very serious infectious disease against which the currently available drugs are limited and not effective. Therefore, there is an urgent need for new chemotherapeutic agents. One attractive drug target is the major cysteine protease from Trypanosoma cruzi, cruzain. In the present work, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted on a series of thiosemicarbazone and semicarbazone derivatives as inhibitors of cruzain. Molecular modeling studies were performed in order to identify the preferred binding mode of the inhibitors into the enzyme active site, and to generate structural alignments for the three-dimensional quantitative structure-activity relationship (3D QSAR) investigations. Statistically significant models were obtained (CoMFA. r(2) = 0.96 and q(2) = 0.78; CoMSIA, r(2) = 0.91 and q(2) = 0.73), indicating their predictive ability for untested compounds. The models were externally validated employing a test set, and the predicted values were in good agreement with the experimental results. The final QSAR models and the information gathered from the 3D CoMFA and CoMSIA contour maps provided important insights into the chemical and structural basis involved in the molecular recognition process of this family of cruzain inhibitors, and should be useful for the design of new structurally related analogs with improved potency. (C) 2009 Elsevier Inc. All rights reserved.

Fragment-Based QSAR and Molecular Modeling Studies on a Series of Discodermolide Analogs as Microtubule-Stabilizing Anticancer Agents

Inhibition of microtubule function is an attractive rational approach to anticancer therapy. Although taxanes are the most prominent among the microtubule-stabilizers, their clinical toxicity, poor pharmacokinetic properties, and resistance have stimulated the search for new antitumor agents having the same mechanism of action. Discodermolide is an example of nontaxane natural product that has the same mechanism of action, demonstrating superior antitumor efficacy and therapeutic index. The extraordinary chemical and biological properties have qualified discodermolide as a lead structure for the design of novel anticancer agents with optimized therapeutic properties. In the present work, we have employed a specialized fragment-based method to develop robust quantitative structure - activity relationship models for a series of synthetic discodermolide analogs. The generated molecular recognition patterns were combined with three-dimensional molecular modeling studies as a fundamental step on the path to understanding the molecular basis of drug-receptor interactions within this important series of potent antitumoral agents.

Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: Evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17 beta-estradiol (E-2) and the selective ER modulator raloxifene (RAL)] from the human ER alpha ligand-binding domain in monomeric and dimeric forms. E-2 dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II`). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I`). Remarkably, ER dimerization strongly suppresses Paths II and II` for E-2 dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

Structural basis for selective inhibition of trypanosomatid glyceraldehyde-3-phosphate dehydrogenase: Molecular docking and 3D QSAR studies

The glycolytic enzyme glyceraldehyde-3 -phosphate dehydrogenase (GAPDH) is as an attractive target for the development of novel antitrypanosomatid agents. In the present work, comparative molecular field analysis and comparative molecular similarity index analysis were conducted on a large series of selective inhibitors of trypanosomatid GAPDH. Four statistically significant models were obtained (r(2) > 0.90 and q(2) > 0.70), indicating their predictive ability for untested compounds. The models were then used to predict the potency of an external test set, and the predicted values were in good agreement with the experimental results. Molecular modeling studies provided further insight into the structural basis for selective inhibition of trypanosomatid GAPDH.

Comparative molecular field analysis of a series of inhibitors of HIV-1 protease

Several protease inhibitors have reached the world market in the last fifteen years, dramatically improving the quality of life and life expectancy of millions of HIV-infected patients. In spite of the tremendous research efforts in this area, resistant HIV-1 variants are constantly decreasing the ability of the drugs to efficiently inhibit the enzyme. As a consequence, inhibitors with novel frameworks are necessary to circumvent resistance to chemotherapy. In the present work, we have created 3D QSAR models for a series of 82 HIV-1 protease inhibitors employing the comparative molecular field analysis (CoMFA) method. Significant correlation coefficients were obtained (q(2) = 0.82 and r(2) = 0.97), indicating the internal consistency of the best model, which was then used to evaluate an external test set containing 17 compounds. The predicted values were in good agreement with the experimental results, showing the robustness of the model and its substantial predictive power for untested compounds. The final QSAR model and the information gathered from the CoMFA contour maps should be useful for the design of novel anti-HIV agents with improved potency.

Molecular characterization of human polyomavirus JC in Brazilian AIDS patients with and without progressive multifocal leukoencephalopathy

Background: JC virus (JCV), the causative agent of progressive multifocal leukoencephalopathy (PML), is classified in 8 different genotypes. Previous reports have suggested a positive association between specific genotypes and PML. Objective: To compare genotypes and adaptive mutations of JCV strains from Brazilian AIDS patients with and without PML. Study design: The VP1 region of JCV was amplified by polymerase chain reaction from cerebrospinal fluid samples from 51 patients with PML and from urine samples of 47 patients with AIDS without central nervous system disease. Genotyping was done by phylogenetic analysis. Amino acid replacement and selection pressures were also investigated. Results: JCV genotype frequency distributions showed that genotypes 2 (32.7%), 1 (26.5%) and 3 (23.5%) were the most prevalent. Genotype 1 had a positive association (p < 0.0001) and genotype 3 showed an inverse association (p < 0.001) with PML. A previously undescribed point mutation at residue 91 (L/I or L/V) and (L/P), non-genotype-associated, was found in 5/49 (10.2%) and 2/47 (4.3%) JCV sequences from PML and non-PML patients, respectively. This mutation was under positive selection only in PML patients. A previously described substitution of T-A in position 128 showed a significant difference between PML and non-PML cases (70% versus 16%, respectively, p < 0.0005). Conclusion: In Brazilian patients with AIDS, JCV genotype 1 showed a strong association with PML (p < 0.0001) and JCV genotype 3 showed an inverse association with PML. The possible association of aminoacids substitution in residues 91 and 128 with PML in patients with AIDS must be further investigated. (C) 2010 Elsevier B.V. All rights reserved.

Covalent attachment of 3,4,9,10-perylenediimides onto the walls of mesoporous molecular sieves MCM-41 and SBA-15

This work describes the covalent grafting of 3,4,9,10-perylenediimides (PDI), which are fluorescent dyes with very interesting optical properties, onto the walls of mesoporous molecular sieves MCM-41 and SBA-15. The mesoporous materials were first treated with 3-aminopropyltriethoxysilane (APTES) in anhydrous toluene, generating amine-containing surfaces. The amine-containing materials were then reacted with 3,4,9,10-perylenetetracarboxylic dianhydride (PTCA), generating surface-grafted PDI. Infrared spectra of the materials showed that the reaction with amino groups took place at both anhydride ends of the PTCA molecule, resulting in surface attached diimides. No sign of unreacted anhydride groups were found. The new materials, designated as MCMN2PDI and SBAN(2)PDI, presented absorption and emission spectra corresponding to weakly coupled PDI chromophores, in contrast to the strongly coupled rings usually found in solid PDI samples. The materials showed a red fluorescence, which could be observed by the naked eye under UV irradiation or with a fluorescence microscope. The PDI-modified mesoporous materials showed electrical conductivity when pressed into a pellet. The results presented here show that the new materials are potentially useful in the design of nanowires. (C) 2007 Elsevier Inc. All rights reserved.

Design of a Modern Liposome and Bee Venom Formulation for the Traditional VIT-Venom Immunotherapy

Traditional venom immunotherapy uses injections of whole bee venom in buffer or adsorbed in Al (OH)(3) in an expensive, time-consuming way. New strategies to improve the safety and efficacy of this treatment with a reduction of injections would, therefore, be of general interest. It would improve patient compliance and provide socio-economic benefits. Liposomes have a long tradition in drug delivery because they increase the therapeutic index and avoid drug degradation and secondary effects. However, bee venom melittin (Mel) and phospholipase (PLA(2)) destroy the phospholipid membranes. Our central idea was to inhibit the PLA(2) and Mel activities through histidine alkylation and or tryptophan oxidation (with pbb, para-bromo-phenacyl bromide, and/or NBSN-bromosuccinimide, respectively) to make their encapsulations possible within stabilized liposomes. We strongly believe that this formulation will be nontoxic but immunogenic. In this paper, we present the whole bee venom conformation characterization during and after chemical modification and after interaction with liposome by ultraviolet, circular dichroism, and fluorescence spectroscopies. The PLA(2) and Mel activities were, measured indirectly by changes in turbidity at 400(nm), rhodamine leak-out, and hemolysis. The native whole bee venom (BV) presented 78.06% of alpha-helical content. The alkylation (A-BV) and succynilation (S-BV) of BV increased 0.44 and 0.20% of its alpha-helical content. The double-modified venom (S-A-BV) had a 0.74% increase of alpha-helical content. The BV chemical modification induced another change on protein conformations observed by Trp that became buried with respect to the native whole BV. It was demonstrated that the liposomal membranes must contain pbb (SPC:Cho:pbb, 26:7:1) as a component to protect them from aggregation and/or fusion. The membranes containing pbb maintained the same turbidity (100%) after incubation with modified venom, in contrast with pbb-free membranes that showed a 15% size decrease. This size decrease was interpreted as membrane degradation and was corroborated by a 50% rhodamine leak-out. Another fact that confirmed our interpretation was the observed 100% inhibition of the hemolytic activity after venom modification with pbb and NBS (S-A-BV). When S-A-BV interacted with liposomes, other protein conformational changes were observed and characterized by the increase of 1.93% on S-A-BV alpha-helical content and the presence of tryptophan residues in a more hydrophobic environment. In other words, the S-A-BV interacted with liposomal membranes, but this interaction was not effective to cause aggregation, leak-out, or fusion. A stable formulation composed by S-A-BV encapsulated within liposomes composed by SPC:Cho:pbb, at a ratio of 26:7:1, was devised. Large unilamellar vesicles of 202.5 nm with a negative surface charge (-24.29 mV) encapsulated 95% of S-A-BV. This formulation can, now, be assayed on VIT.

Approaches for multicopper oxidases in the design of electrochemical sensors for analytical applications

We report an effective approach for the construction of a biomimetic sensor of multicopper oxidases by immobilizing a cyclic-tetrameric copper(II) species, containing the ligand (4-imidazolyl)ethylene-2-amino-1-ethylpyridine (apyhist), in the Nafion (R) membrane on a vitreous carbon electrode surface. This complex provides a tetranuclear arrangement of copper ions that allows an effective reduction of oxygen to water, in a catalytic cycle involving four electrons. The electrochemical reduction of oxygen was studied at pH 9.0 buffer solution by using cyclic voltammetry, chronoamperometry, rotating disk electrode voltammetry and scanning electrochemical microscopy techniques. The mediator shows good electrocatalytic ability for the reduction of O(2) at pH 9.0, with reduction of overpotential (350 mV) and increased current response in comparison with results obtained with a bare glassy carbon electrode. The heterogeneous rate constant (k(ME)`) for the reduction of O(2) at the modified electrode was determined by using a Koutecky-Levich plot. In addition, the charge transport rate through the coating and the apparent diffusion coefficient of O(2) into the modifier film were also evaluated. The overall process was found to be governed by the charge transport through the coating, occurring at the interface or at a finite layer at the electrode/coating interface. The proposed study opens up the way for the development of bioelectronic devices based on molecular recognition and self-organization. (C) 2010 Elsevier Ltd. All rights reserved.

Design of novel antituberculosis compounds using graph-theoretical and substructural approaches

The increasing resistance of Mycobacterium tuberculosis to the existing drugs has alarmed the worldwide scientific community. In an attempt to overcome this problem, two models for the design and prediction of new antituberculosis agents were obtained. The first used a mixed approach, containing descriptors based on fragments and the topological substructural molecular design approach (TOPS-MODE) descriptors. The other model used a combination of two-dimensional (2D) and three-dimensional (3D) descriptors. A data set of 167 compounds with great structural variability, 72 of them antituberculosis agents and 95 compounds belonging to other pharmaceutical categories, was analyzed. The first model showed sensitivity, specificity, and accuracy values above 80% and the second one showed values higher than 75% for these statistical indices. Subsequently, 12 structures of imidazoles not included in this study were designed, taking into account the two models. In both cases accuracy was 100%, showing that the methodology in silico developed by us is promising for the rational design of antituberculosis drugs.

Insights into the Molecular Requirements for the Anti-obesity Activity of a Series of CB1 Ligands

Two-dimensional and 3D quantitative structure-activity relationships studies were performed on a series of diarylpyridines that acts as cannabinoid receptor ligands by means of hologram quantitative structure-activity relationships and comparative molecular field analysis methods. The quantitative structure-activity relationships models were built using a data set of 52 CB1 ligands that can be used as anti-obesity agents. Significant correlation coefficients (hologram quantitative structure-activity relationships: r 2 = 0.91, q 2 = 0.78; comparative molecular field analysis: r 2 = 0.98, q 2 = 0.77) were obtained, indicating the potential of these 2D and 3D models for untested compounds. The models were then used to predict the potency of an external test set, and the predicted (calculated) values are in good agreement with the experimental results. The final quantitative structure-activity relationships models, along with the information obtained from 2D contribution maps and 3D contour maps, obtained in this study are useful tools for the design of novel CB1 ligands with improved anti-obesity potency.

Docking and molecular dynamics simulation of quinone compounds with trypanocidal activity

In this work, two different docking programs were used, AutoDock and FlexX, which use different types of scoring functions and searching methods. The docking poses of all quinone compounds studied stayed in the same region in the trypanothione reductase. This region is a hydrophobic pocket near to Phe396, Pro398 and Leu399 amino acid residues. The compounds studied displays a higher affinity in trypanothione reductase (TR) than glutathione reductase (GR), since only two out of 28 quinone compounds presented more favorable docking energy in the site of human enzyme. The interaction of quinone compounds with the TR enzyme is in agreement with other studies, which showed different binding sites from the ones formed by cysteines 52 and 58. To verify the results obtained by docking, we carried out a molecular dynamics simulation with the compounds that presented the highest and lowest docking energies. The results showed that the root mean square deviation (RMSD) between the initial and final pose were very small. In addition, the hydrogen bond pattern was conserved along the simulation. In the parasite enzyme, the amino acid residues Leu399, Met400 and Lys402 are replaced in the human enzyme by Met406, Tyr407 and Ala409, respectively. In view of the fact that Leu399 is an amino acid of the Z site, this difference could be explored to design selective inhibitors of TR.

Desenvolvimento de redes de polímeros interpenetrantes (IPN) para a aplicação como eletrólito polimérico

Uma nova rede de polímeros interpenetrantes (IPN) baseada em poliuretana de óleo de mamona e poli(etileno glicol) e poli(metacrilato de metila) foi preparada para ser utilizada como eletrólito polimérico. Os seguintes parâmetros de polimerização foram avaliados: massa molecular do poli(etileno glicol) (PEG), concentração de PEG e concentração de metacrilato de metila. As membranas de IPN foram caracterizadas por calorimetria diferencial de varredura (DSC) e espectroscopia de infravermelho por transformada de Fourier (FT-IR). Os eletrólitos de redes de polímeros interpenetrantes (IPNE) foram preparados a partir da dopagem com sal de lítio através do inchamento numa solução de 10% em massa de LiClO4 na mistura de carbonato de etileno e carbonato de propileno na razão mássica de 50:50. As IPNEs foram caracterizadas por espectroscopia de impedância eletroquímica e Raman. As IPNEs foram testadas como eletrólito polimérico em supercapacitores. As células capacitivas foram preparadas utilizando eletrodos de polipirrol (PPy). Os valores de capacitância e eficiência foram calculados por impedância eletroquímica, voltametria cíclica e ciclos galvonostáticos de carga e descarga. Os valores de capacitância obtidos foram em torno de 90 F.g-1 e eficiência variou no intervalo de 88 a 99%. Os valores de densidade de potência foram superiores a 250 W.kg-1 enquanto que a densidade de energia variou de 10 a 33 W.h.kg-1, dependendo da composição da IPNE. As características eletroquímicas do eletrólito formado pela IPN-LiClO4 (IPNE) foram comparadas aos eletrólitos poliméricos convencionais, tais como poli(difluoreto de vinilideno)-(hexafluorpropileno) ((PVDF-HFP/LiClO4) e poliuretana comercial (Bayer desmopan 385) (PU385/LiClO4). As condutividades na temperatura ambiente foram da ordem de 10-3 S.cm-1. A capacitância da célula utilizando eletrodos de PPy com eletrólito de PVDFHFP foi de 115 F.g-1 (30 mF.cm-2) e 110 F.g-1 (25 mF.cm-2) para a célula com PU385 comparadas a 90 F.g-1 (20 mF.cm-2) para a IPNE. Os capacitores preparados com eletrólito de IPNE apresentaram valores de capacitância inferior aos demais, entretanto provaram ser mais estáveis e mais resistentes aos ciclos de carga/descarga. A interpenetração de duas redes poliméricas, PU e PMMA produziu um eletrólito com boa estabilidade mecânica e elétrica. Um protótipo de supercapacitor de estado sólido foi produzindo utilizando eletrodos impressos de carbono ativado (PCE) e o eletrólito polimérico de IPNE. A técnica de impressão de carbono possui várias vantagens em relação aos outros métodos de manufatura de eletrodos de carbono, pois a área do eletrodo, espessura e composição são variáveis que podem ser controladas experimentalmente. As células apresentaram uma larga janela eletroquímica (4V) e valores da capacitância da ordem de 113 mF.cm-2 (16 F.g-1). Métodos alternativos de preparação do PCE investigados incluem o uso de IPNE como polímero de ligação ao carbono ativado, estes eletrodos apresentaram valores de capacitância similares aos produzidos com PVDF. A influência do número de camadas de carbono usadas na produção do PCE também foi alvo de estudo. Em relação ao eletrólito polimérico, o plastificante e o sal de lítio foram adicionados durante a síntese, formando a IPNGel. As células apresentaram alta capacitância e boa estabilidade após 4000 ciclos de carga e descarga. As membranas de IPN foram testadas também como reservatório de medicamento em sistemas de transporte transdérmico por iontoforese. Os filmes, mecanicamente estáveis, formaram géis quando inchado em soluções saturadas de lidocaina.HCl, anestésico local, em propileno glicol (PG), poli(etileno glicol) (PEG400) e suas misturas. O grau de inchamento em PG foi de 15% e 35% em PEG400. Agentes químicos de penetração foram utilizados para diminuir a resistência da barreira causada pela pele, dentre eles o próprio PG, a 2-pirrolidinona (E1) e a 1-dodecil-2-pirrolidinona (E2). Os géis foram caracterizados por espectroscopia de impedância eletroquímica e transporte passivo e por iontoforese através de uma membrana artificial (celofane). O sistema IPN/ lidocaina.HCl apresentou uma correlação linear entre medicamento liberado e a corrente aplicada. Os melhores resultados de transporte de medicamento foram obtidos utilizando o PG como solvente.