Molecular models of NS3 protease variants of the Hepatitis C virus

Background: Hepatitis C virus (HCV) currently infects approximately three percent of the world population. In view of the lack of vaccines against HCV, there is an urgent need for an efficient treatment of the disease by an effective antiviral drug. Rational drug design has not been the primary way for discovering major therapeutics. Nevertheless, there are reports of success in the development of inhibitor using a structure-based approach. One of the possible targets for drug development against HCV is the NS3 protease variants. Based on the three-dimensional structure of these variants we expect to identify new NS3 protease inhibitors. In order to speed up the modeling process all NS3 protease variant models were generated in a Beowulf cluster. The potential of the structural bioinformatics for development of new antiviral drugs is discussed.Results: the atomic coordinates of crystallographic structure 1CU1 and 1DY9 were used as starting model for modeling of the NS3 protease variant structures. The NS3 protease variant structures are composed of six subdomains, which occur in sequence along the polypeptide chain. The protease domain exhibits the dual beta-barrel fold that is common among members of the chymotrypsin serine protease family. The helicase domain contains two structurally related beta-alpha-beta subdomains and a third subdomain of seven helices and three short beta strands. The latter domain is usually referred to as the helicase alpha-helical subdomain. The rmsd value of bond lengths and bond angles, the average G-factor and Verify 3D values are presented for NS3 protease variant structures.Conclusions: This project increases the certainty that homology modeling is an useful tool in structural biology and that it can be very valuable in annotating genome sequence information and contributing to structural and functional genomics from virus. The structural models will be used to guide future efforts in the structure-based drug design of a new generation of NS3 protease variants inhibitors. All models in the database are publicly accessible via our interactive website, providing us with large amount of structural models for use in protein-ligand docking analysis.

Shikimate kinase: A potential target for development of novel antitubercular agents

Tuberculosis (TB) remains the leading cause of mortality due to a bacterial pathogen, Mycobacterium tuberculosis. However, no new classes of drugs for TB have been developed in the past 30 years. Therefore there is an urgent need to develop faster acting and effective new antitubercular agents, preferably belonging to new structural classes, to better combat TB, including MDR-TB, to shorten the duration of current treatment to improve patient compliance, and to provide effective treatment of latent tuberculosis infection. The enzymes in the shikimate pathway are potential targets for development of a new generation of antitubercular drugs. The shikimate pathway has been shown by disruption of aroK gene to be essential for the Mycobacterium tuberculosis. The shikimate kinase (SK) catalyses the phosphorylation of the 3-hydroxyl group of shikimic acid (shikimate) using ATP as a co-substrate. SK belongs to family of nucleoside monophosphate (NMP) kinases. The enzyme is an alpha/beta protein consisting of a central sheet of five parallel beta-strands flanked by alpha-helices. The shikimate kinases are composed of three domains: Core domain, Lid domain and Shikimate-binding domain. The Lid and Shikimate-binding domains are responsible for large conformational changes during catalysis. More recently, the precise interactions between SK and substrate have been elucidated, showing the binding of shikimate with three charged residues conserved among the SK sequences. The elucidation of interactions between MtSK and their substrates is crucial for the development of a new generation of drugs against tuberculosis through rational drug design.

OS produtos naturais e a química medicinal moderna

Natural products have been utilized by humans since ancient times and the relief and cure of their diseases was the first purpose for using natural products in medicine. The history of the oriental and occidental civilizations is very rich in examples of the utilization of natural products in medicine and health care. Chinese traditional medicine is one of the most important examples of how natural products can be efficient in the treatment of diseases, and it points to the importance of scientific research on natural products, concerning the discovery of new active chemical entities. The complexity, chemical diversity and biological properties of natural products always fascinated people, and during the last 200 years, this led to the discovery of important new drugs. In the last 30 years, the development of new bioassay techniques, biotechnology methods, bio-guided phytochemical studies, automated high throughput screening and high performance analytical methods, have introduced new concepts and possibilities of rational drug design and drug discovery. In this context, natural products have played an important and decisive role in the development of modern medicinal chemistry.

Thiophenecarboxamide derivatives activated by EthA kill mycobacterium tuberculosis by inhibiting the CTP synthetase PyrG

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Structure-activity relationship and mechanism of action studies of manzamine analogues for the control of neuroinflammation and cerebral infections

Structure-activity relationship studies were carried out by chemical modification of manzamine A (1), 8-hydroxymanzamine A (2), manzamine F (14), and ircinal isolated from the sponge Acanthostrongylophora. The derived analogues were evaluated for antimalarial, antimicrobial, and antineuroinflammatory activities. Several modified products exhibited potent and improved in vitro antineuroinflammatory, antimicrobial, and antimalarial activity. 1 showed improved activity against malaria compared to chloroquine in both multi- and single-dose in vivo experiments. The significant antimalarial potential was revealed by a 100% cure rate of malaria in mice with one administration of 100 mg/kg of 1. The potent antineuroinflammatory activity of the manzamines will provide great benefit for the prevention and treatment of cerebral infections (e.g., Cryptococcus and Plasmodium). In addition, 1 was shown to permeate across the blood-brain barrier (BBB) in an in vitro model using a MDR-MDCK monolayer. Docking studies support that 2 binds to the ATP-noncompetitive pocket of glycogen synthesis kinase-3beta (GSK-3beta), which is a putative target of manzamines. On the basis of the results presented here, it will be possible to initiate rational drug design efforts around this natural product scaffold for the treatment of several different diseases.

Artemisinin, an endoperoxide antimalarial, disrupts the hemoglobin catabolism and heme detoxification systems in malarial parasite

Endoperoxide antimalarials based on the ancient Chinese drug Qinghaosu (artemisinin) are currently our major hope in the fight against drug-resistant malaria. Rational drug design based on artemisinin and its analogues is slow as the mechanism of action of these antimalarials is not clear. Here we report that these drugs, at least in part, exert their effect by interfering with the plasmodial hemoglobin catabolic pathway and inhibition of heme polymerization. In an in vitro experiment we observed inhibition of digestive vacuole proteolytic activity of malarial parasite by artemisinin. These observations were further confirmed by ex vivo experiments showing accumulation of hemoglobin in the parasites treated with artemisinin, suggesting inhibition of hemoglobin degradation. We found artemisinin to be a potent inhibitor of heme polymerization activity mediated by Plasmodium yoelii lysates as well as Plasmodium falciparum histidine-rich protein II. Interaction of artemisinin with the purified malarial hemozoin in vitro resulted in the concentration-dependent breakdown of the malaria pigment. Our results presented here may explain the selective and rapid toxicity of these drugs on mature, hemozoin-containing, stages of malarial parasite. Since artemisinin and its analogues appear to have similar molecular targets as chloroquine despite having different structures, they can potentially bypass the quinoline resistance machinery of the malarial parasite, which causes sublethal accumulation of these drugs in resistant strains.

Bifurcation and singularity analysis of a molecular network for the induction of long-term memory.

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis.

Targeting TRAF6 for Cancer Therapeutical Development

Tumor necrosis factor (TNF)-Receptor Associated Factors (TRAFs) are a family of signal transducer proteins. TRAF6 is a unique member of this family in that it is involved in not only the TNF superfamily, but the toll-like receptor (TLR)/IL-1R (TIR) superfamily. The formation of the complex consisting of Receptor Activator of Nuclear Factor κ B (RANK), with its ligand (RANKL) results in the recruitment of TRAF6, which activates NF-κB, JNK and MAP kinase pathways. TRAF6 is critical in signaling with leading to release of various growth factors in bone, and promotes osteoclastogenesis. TRAF6 has also been implicated as an oncogene in lung cancer and as a target in multiple myeloma. In the hopes of developing small molecule inhibitors of the TRAF6-RANK interaction, multiple steps were carried out. Computational prediction of hot spot residues on the protein-protein interaction of TRAF6 and RANK were examined. Three methods were used: Robetta, KFC2, and HotPoint, each of which uses a different methodology to determine if a residue is a hot spot. These hot spot predictions were considered the basis for resolving the binding site for in silico high-throughput screening using GOLD and the MyriaScreen database of drug/lead-like compounds. Computationally intensive molecular dynamics simulations highlighted the binding mechanism and TRAF6 structural changes upon hit binding. Compounds identified as hits were verified using a GST-pull down assay, comparing inhibition to a RANK decoy peptide. Since many drugs fail due to lack of efficacy and toxicity, predictive models for the evaluation of the LD50 and bioavailability of our TRAF6 hits, and these models can be used towards other drugs and small molecule therapeutics as well. Datasets of compounds and their corresponding bioavailability and LD50 values were curated based, and QSAR models were built using molecular descriptors of these compounds using the k-nearest neighbor (k-NN) method, and quality of these models were cross-validated.

Structure-activity relationships for a new family of sulfonylurea herbicides

Acetohydroxyacid synthase (AHAS; EC 2.2.1.6) catalyzes the first common step in branched-chain amino acid biosynthesis. The enzyme is inhibited by several chemical classes of compounds and this inhibition is the basis of action of the sulfonylurea and imidazolinone herbicides. The commercial sulfonylureas contain a pyrimidine or a triazine ring that is substituted at both meta positions, thus obeying the initial rules proposed by Levitt. Here we assess the activity of 69 monosubstituted sulfonylurea analogs and related compounds as inhibitors of pure recombinant Arabidopsis thaliana AHAS and show that disubstitution is not absolutely essential as exemplified by our novel herbicide, monosulfuron (2-nitro-N-(4'-methyl-pyrimidin-2'-yl) phenyl-sulfonylurea), which has a pyrimidine ring with a single meta substituent. A subset of these compounds was tested for herbicidal activity and it was shown that their effect in vivo correlates well with their potency in vitro as AHAS inhibitors. Three-dimensional quantitative structure-activity relationships were developed using comparative molecular field analysis and comparative molecular similarity indices analysis. For the latter, the best result was obtained when steric, electrostatic, hydrophobic and H-bond acceptor factors were taken into consideration. The resulting fields were mapped on to the published crystal structure of the yeast enzyme and it was shown that the steric and hydrophobic fields are in good agreement with sulfonylurea-AHAS interaction geometry.

Attenuated total reflection-FT-IR spectroscopic imaging of protein crystallization

Protein crystallization is of strategic and commercial relevance in the post-genomic era because of its pivotal role in structural proteomics projects. Although protein structures are crucial for understanding the function of proteins and to the success of rational drug design and other biotechnology applications, obtaining high quality crystals is a major bottleneck to progress. The major means of obtaining crystals is by massive-scale screening of a target protein solution with numerous crystallizing agents. However, when crystals appear in these screens, one cannot easily know if they are crystals of protein, salt, or any other molecule that happens to be present in the trials. We present here a method based on Attenuated Total Reflection (ATR)-FT-IR imaging that reliably identifies protein crystals through a combination of chemical specificity and the visualizing capability of this approach, thus solving a major hurdle in protein crystallization. ATR-FT-IR imaging was successfully applied to study the crystallization of thaumatin and lysozyme in a high-throughput manner, simultaneously from six different solutions. This approach is fast as it studies protein crystallization in situ and provides an opportunity to examine many different samples under a range of conditions.

Extracellular loops 1 and 3 and their associated transmembrane regions of the calcitonin receptor-like receptor are needed for CGRP receptor function

The first and third extracellular loops (ECL) of G protein-coupled receptors (GPCRs) have been implicated in ligand binding and receptor function. This study describes the results of an alanine/leucine scan of ECLs 1 and 3 and loop-associated transmembrane (TM) domains of the secretin-like GPCR calcitonin receptor-like receptor which associates with receptor activity modifying protein 1 to form the CGRP receptor. Leu195Ala, Val198Ala and Ala199Leu at the top of TM2 all reduced aCGRP-mediated cAMP production and internalization; Leu195Ala and Ala199Leu also reduced aCGRP binding. These residues form a hydrophobic cluster within an area defined as the "minor groove" of rhodopsin-like GPCRs. Within ECL1, Ala203Leu and Ala206Leu influenced the ability of aCGRP to stimulate adenylate cyclase. In TM3, His219Ala, Leu220Ala and Leu222Ala have influences on aCGRP binding and cAMP production; they are likely to indirectly influence the binding site for aCGRP as well as having an involvement in signal transduction. On the exofacial surfaces of TMs 6 and 7, a number of residues were identified that reduced cell surface receptor expression, most noticeably Leu351Ala and Glu357Ala in TM6. The residues may contribute to the RAMP1 binding interface. Ile360Ala impaired aCGRP-mediated cAMP production. Ile360 is predicted to be located close to ECL2 and may facilitate receptor activation. Identification of several crucial functional loci gives further insight into the activation mechanism of this complex receptor system and may aid rational drug design.

Rational design of serine protease inhibitors

Proteases regulate a spectrum of diverse physiological processes, and dysregulation of proteolytic activity drives a plethora of pathological conditions. Understanding protease function is essential to appreciating many aspects of normal physiology and progression of disease. Consequently, development of potent and specific inhibitors of proteolytic enzymes is vital to provide tools for the dissection of protease function in biological systems and for the treatment of diseases linked to aberrant proteolytic activity. The studies in this thesis describe the rational design of potent inhibitors of three proteases that are implicated in disease development. Additionally, key features of the interaction of proteases and their cognate inhibitors or substrates are analysed and a series of rational inhibitor design principles are expounded and tested. Rational design of protease inhibitors relies on a comprehensive understanding of protease structure and biochemistry. Analysis of known protease cleavage sites in proteins and peptides is a commonly used source of such information. However, model peptide substrate and protein sequences have widely differing levels of backbone constraint and hence can adopt highly divergent structures when binding to a protease’s active site. This may result in identical sequences in peptides and proteins having different conformations and diverse spatial distribution of amino acid functionalities. Regardless of this, protein and peptide cleavage sites are often regarded as being equivalent. One of the key findings in the following studies is a definitive demonstration of the lack of equivalence between these two classes of substrate and invalidation of the common practice of using the sequences of model peptide substrates to predict cleavage of proteins in vivo. Another important feature for protease substrate recognition is subsite cooperativity. This type of cooperativity is commonly referred to as protease or substrate binding subsite cooperativity and is distinct from allosteric cooperativity, where binding of a molecule distant from the protease active site affects the binding affinity of a substrate. Subsite cooperativity may be intramolecular where neighbouring residues in substrates are interacting, affecting the scissile bond’s susceptibility to protease cleavage. Subsite cooperativity can also be intermolecular where a particular residue’s contribution to binding affinity changes depending on the identity of neighbouring amino acids. Although numerous studies have identified subsite cooperativity effects, these findings are frequently ignored in investigations probing subsite selectivity by screening against diverse combinatorial libraries of peptides (positional scanning synthetic combinatorial library; PS-SCL). This strategy for determining cleavage specificity relies on the averaged rates of hydrolysis for an uncharacterised ensemble of peptide sequences, as opposed to the defined rate of hydrolysis of a known specific substrate. Further, since PS-SCL screens probe the preference of the various protease subsites independently, this method is inherently unable to detect subsite cooperativity. However, mean hydrolysis rates from PS-SCL screens are often interpreted as being comparable to those produced by single peptide cleavages. Before this study no large systematic evaluation had been made to determine the level of correlation between protease selectivity as predicted by screening against a library of combinatorial peptides and cleavage of individual peptides. This subject is specifically explored in the studies described here. In order to establish whether PS-SCL screens could accurately determine the substrate preferences of proteases, a systematic comparison of data from PS-SCLs with libraries containing individually synthesised peptides (sparse matrix library; SML) was carried out. These SML libraries were designed to include all possible sequence combinations of the residues that were suggested to be preferred by a protease using the PS-SCL method. SML screening against the three serine proteases kallikrein 4 (KLK4), kallikrein 14 (KLK14) and plasmin revealed highly preferred peptide substrates that could not have been deduced by PS-SCL screening alone. Comparing protease subsite preference profiles from screens of the two types of peptide libraries showed that the most preferred substrates were not detected by PS SCL screening as a consequence of intermolecular cooperativity being negated by the very nature of PS SCL screening. Sequences that are highly favoured as result of intermolecular cooperativity achieve optimal protease subsite occupancy, and thereby interact with very specific determinants of the protease. Identifying these substrate sequences is important since they may be used to produce potent and selective inhibitors of protolytic enzymes. This study found that highly favoured substrate sequences that relied on intermolecular cooperativity allowed for the production of potent inhibitors of KLK4, KLK14 and plasmin. Peptide aldehydes based on preferred plasmin sequences produced high affinity transition state analogue inhibitors for this protease. The most potent of these maintained specificity over plasma kallikrein (known to have a very similar substrate preference to plasmin). Furthermore, the efficiency of this inhibitor in blocking fibrinolysis in vitro was comparable to aprotinin, which previously saw clinical use to reduce perioperative bleeding. One substrate sequence particularly favoured by KLK4 was substituted into the 14 amino acid, circular sunflower trypsin inhibitor (SFTI). This resulted in a highly potent and selective inhibitor (SFTI-FCQR) which attenuated protease activated receptor signalling by KLK4 in vitro. Moreover, SFTI-FCQR and paclitaxel synergistically reduced growth of ovarian cancer cells in vitro, making this inhibitor a lead compound for further therapeutic development. Similar incorporation of a preferred KLK14 amino acid sequence into the SFTI scaffold produced a potent inhibitor for this protease. However, the conformationally constrained SFTI backbone enforced a different intramolecular cooperativity, which masked a KLK14 specific determinant. As a consequence, the level of selectivity achievable was lower than that found for the KLK4 inhibitor. Standard mechanism inhibitors such as SFTI rely on a stable acyl-enzyme intermediate for high affinity binding. This is achieved by a conformationally constrained canonical binding loop that allows for reformation of the scissile peptide bond after cleavage. Amino acid substitutions within the inhibitor to target a particular protease may compromise structural determinants that support the rigidity of the binding loop and thereby prevent the engineered inhibitor reaching its full potential. An in silico analysis was carried out to examine the potential for further improvements to the potency and selectivity of the SFTI-based KLK4 and KLK14 inhibitors. Molecular dynamics simulations suggested that the substitutions within SFTI required to target KLK4 and KLK14 had compromised the intramolecular hydrogen bond network of the inhibitor and caused a concomitant loss of binding loop stability. Furthermore in silico amino acid substitution revealed a consistent correlation between a higher frequency of formation and the number of internal hydrogen bonds of SFTI-variants and lower inhibition constants. These predictions allowed for the production of second generation inhibitors with enhanced binding affinity toward both targets and highlight the importance of considering intramolecular cooperativity effects when engineering proteins or circular peptides to target proteases. The findings from this study show that although PS-SCLs are a useful tool for high throughput screening of approximate protease preference, later refinement by SML screening is needed to reveal optimal subsite occupancy due to cooperativity in substrate recognition. This investigation has also demonstrated the importance of maintaining structural determinants of backbone constraint and conformation when engineering standard mechanism inhibitors for new targets. Combined these results show that backbone conformation and amino acid cooperativity have more prominent roles than previously appreciated in determining substrate/inhibitor specificity and binding affinity. The three key inhibitors designed during this investigation are now being developed as lead compounds for cancer chemotherapy, control of fibrinolysis and cosmeceutical applications. These compounds form the basis of a portfolio of intellectual property which will be further developed in the coming years.

The rational design of vaccines

This review provides an insight into the various opportunities for vaccine intervention, analysis of strategies for vaccine development, vaccine ability to modulate immune responses and resultant rational vaccine design. In addition, wider aspects are considered, such as biotechnological advances, advances in immunological understanding and host-pathogen interactions. The key question addressed here is, with all our research and understanding, have we reached a new echelon in vaccine development, that of rational design? ©2005 Elsevier Ltd. All rights reserved.

Network approaches to drug discovery

Introduction: Advances in genomics technologies are providing a very large amount of data on genome-wide gene expression profiles, protein molecules and their interactions with other macromolecules and metabolites. Molecular interaction networks provide a useful way to capture this complex data and comprehend it. Networks are beginning to be used in drug discovery, in many steps of the modern discovery pipeline, with large-scale molecular networks being particularly useful for the understanding of the molecular basis of the disease. Areas covered: The authors discuss network approaches used for drug target discovery and lead identification in the drug discovery pipeline. By reconstructing networks of targets, drugs and drug candidates as well as gene expression profiles under normal and disease conditions, the paper illustrates how it is possible to find relationships between different diseases, find biomarkers, explore drug repurposing and study emergence of drug resistance. Furthermore, the authors also look at networks which address particular important aspects such as off-target effects, combination-targets, mechanism of drug action and drug safety. Expert opinion: The network approach represents another paradigm shift in drug discovery science. A network approach provides a fresh perspective of understanding important proteins in the context of their cellular environments, providing a rational basis for deriving useful strategies in drug design. Besides drug target identification and inferring mechanism of action, networks enable us to address new ideas that could prove to be extremely useful for new drug discovery, such as drug repositioning, drug synergy, polypharmacology and personalized medicine.

Feleucin-BO1: A Novel Antimicrobial Non-Apeptide Amide from the Skin Secretion of the Toad, Bombina orientalis, and Design of a Potent Broad-Spectrum Synthetic Analogue, Feleucin-K3

Feleucins-BV1 and -BV2 are recently-described prototypes of a novel antimicrobial nonapeptide (AMP) family identified in the skin secretion of the bombinid toad, Bombina variegata. They are encoded on different precursors that also encode a novel bombinin. Here we describe the identification of feleucin-BO1 (FLGLLGSLLamide) which is co-encoded with a different novel bombinin, named feleucin precursor-associated bombinin (FPA-bombinin-BO), from the skin secretion of Bombina orientalis. Synthetic feleucin-BO1 displayed activity against a reference Gram-positive bacterium. Staphylococcus aureus (MIC 34 μM) but was inactive (> 250 μM) against the Gram-negative bacterium, Escherichia coli, and the yeast, Candida albicans. This pattern of activity was similar to that of the prototypes. Design and synthesis of a cationicity-enhanced analogue, feleucin-K3 (F-K3), in which the amino acid residues at positions 3 (G), 6 (G) and 7 (S) of feleucin-BO1 were substituted with Lys (K) residues, resulted in a peptide with significantly-enhanced potency and spectrum of activity. The MICs of F-K3 against the reference microorganisms were 7 μM (S. aureus), 14 μM (E. coli) and 7 μM (C. albicans). These data indicate that the skin secretions of amphibians can continue to provide novel peptide templates for the rational design of analogues with possible therapeutic utility.