Search

(Chymo)trypsin-like serine fold proteases belong to the serine/cysteine proteases found in eukaryotes, prokaryotes, and viruses. Their catalytic activity is carried out using a triad of amino acids, a nucleophile, a base, and an acid. For this superfamily of proteases, we propose the existence of a universal 3D structure comprising 11 amino acids near the catalytic nucleophile and base – Nucleophile-Base Catalytic Zone (NBCZone). The comparison of NBCZones among 169 eukaryotic, prokaryotic, and viral (chymo)trypsin-like proteases suggested the existence of 15 distinct groups determined by the combination of amino acids located at two “key” structure-functional positions 54T and 55T near the catalytic base His57T. Most eukaryotic and prokaryotic proteases fell into two major groups, [ST]A and TN. Usually, proteases of [ST]A group contain a disulfide bond between cysteines Cys42T and Cys58T of the NBCZone. In contrast, viral proteases were distributed among seven groups, and lack this disulfide bond. Furthermore, only the [ST]A group of eukaryotic proteases contains glycine at position 43T, which is instrumental for activation of these enzymes. In contrast, due to the side chains of residues at position 43T prokaryotic and viral proteases do not have the ability to carry out the structural transition of the eukaryotic zymogen-zyme type.

The ongoing pandemic caused by SARS-CoV-2 has gone through various phases. From the initial outbreak the virus has mutated several times, with some lineages showing even stronger infectivity and faster spread than the original virus. Among all the variants, beta, gamma, delta and the latest (omicron) are currently classified as variants of concern (VOC) while the remaining are labelled either as variants of interest (VOI) or variants under monitoring (VUM). In this work, we have focused on the mutations observed in important variants, particularly at the receptor-binding domain (RBD) of the spike protein that is responsible for the interactions with the host ACE2 receptor and binding of antibodies. Studying these mutations is particularly important for understanding the viral infectivity, spread of the disease and for tracking the escape routes of this virus from antibodies. Molecular dynamics (MD) based alchemical free energy calculations have been shown to be very accurate in predicting the free energy change due to a mutation that could have a deleterious or a stabilising effect on the protein itself or its binding affinity to another protein. Here, we investigated the significance of six commonly observed spike RBD mutations on the stability of the spike protein binding to ACE2 by free energy calculations using high throughput MD simulations. For comparison, we also used other (rigorous and non-rigorous) binding free energy prediction methods and compared our results with the experimental data if available. The alchemical free energy-based method consistently predicted the free-energy changes with an accuracy close to ±1.0 kcal/mol when compared with the available experimental values. As per our simulation data the most significant mutations responsible for stabilising the spike RBD interactions with human ACE2 are N501Y and L452R.

Due to the ever-increasing antimicrobial resistance there is an urgent need to continuously design and develop novel antimicrobial agents. Inspired by the broad antibacterial activities of various heterocyclic compounds such as 2-quinolone derivatives, we designed and synthesized new methyl-(2-oxo-1,2-dihydroquinolin-4-yl)-L-alaninate-1,2,3-triazole derivatives via 1,3-dipolar cycloaddition reaction of 1-propargyl-2-quinolone-L-alaninate with appropriate azide groups. The synthesized compounds were obtained in good yield ranging from 75 to 80 %. The chemical structures of these novel hybrid molecules were determined by spectroscopic methods and the antimicrobial activity of the compounds was investigated against both bacterial and fungal strains. The tested compounds showed significant antimicrobial activity and weak to moderate antifungal activity. Despite the evident similarity of the quinolone moiety of our compounds with fluoroquinolones, our compounds do not function by inhibiting DNA gyrase. Computational characterization of the compounds shows that they have attractive physicochemical and pharmacokinetic properties and could serve as templates for developing potential antimicrobial agents for clinical use.

Coagulation, complement, and innate immunity are tightly interwoven and form an alliance that can be traced back to early eukaryotic evolution. Here we employed an ecoimmunological approach using Tissue Factor Pathway Inhibitor (TFPI)-1-derived peptides from the different classes of vertebrates (i.e. fish, reptile, bird, and mammals) and tested whether they can boost killing of various human bacterial pathogens in plasma. We found signs of species-specific conservation and diversification during evolution in these peptides that significantly impact their antibacterial activity. Though all peptides tested executed bactericidal activity in mammalian plasma (with the exception of rodents), no killing was observed in plasma from birds, reptiles, and fish, pointing to a crucial role for the classical pathway of the complement system. We also observed an interference of these peptides with the human intrinsic pathway of coagulation though, unlike complement activation, this mechanism appears not to be evolutionary conserved.

Nanoemulsions are a rapidly growing drug delivery technology capable of increasing a drug's aqueous solubility and stability. A novel oil‐in‐water nanoemulsion using a polymer, poly(decalactone) (PDL), instead of a conventional oil was recently reported. The amount of drug loading in a polymer‐based formulation is mainly governed by the drug's solubility in the polymer. Thus, herein the power of molecular dynamics simulations (MDS) for the calculation of the Hildebrand solubility parameter to predict PDL–drug miscibility is utilized. The MDS results are subsequently verified by formulating a PDL nanoemulsion with a dispersed droplet size of less than 200 nm by using a block copolymer of PDL (mPEG‐b‐PDL) as a surfactant, with seven different hydrophobic drug molecules. The MDS results are consistent with the experimental findings in terms of increment in the aqueous solubilities of the drugs in PDL nanoemulsion, where celecoxib demonstrated the highest while methotrexate exhibited the lowest solubility increment. Consequently, the reported MDS method can be utilized to predict a drug's solubility/miscibility in PDL to estimate the level of drug loading. The MDS facilitated screening of drugs could consequently emerge as an efficient approach in designing PDL nanoemulsions stabilized by mPEG‐b‐PDL or other similar systems.

Drugs that target human thymidylate synthase (hTS) are widely used in anti-cancer therapy. However, treatment with classical substrate-site-directed TS inhibitors induces its over-expression and the development of drug resistance. We thus pursued an alternative strategy that led to the discovery of TS-dimer disrupters that bind at the monomer-monomer interface and shift the dimerization equilibrium of both the recombinant and the intracellular protein toward the inactive monomers. We performed a structural, spectroscopic and kinetic investigation of the effects of these small molecules andthe best one, E7, accelerates the proteasomal degradation of hTS in cancer cells. E7 showed a superior anticancer profile to fluorouracil in a mouse model of human pancreatic and ovarian cancer. Thus, over sixty years after the discovery of the first TS prodrug inhibitor, fluorouracil, E7 breaks the link between TS inhibition and enhanced expression in response, providing a strategy to fight drug-resistant cancers.

Approaches to increase the efficiency in developing drugs and diagnostics tools, including new drug delivery and diagnostic technologies, are needed for improved diagnosis and treatment of major diseases and health problems such as cancer, inflammatory diseases, chronic wounds, and antibiotic resistance. Development within several areas of research ranging from computational sciences, material sciences, bioengineering to biomedical sciences and bioimaging is needed to realize innovative drug development and diagnostic (DDD) approaches. Here, an overview of recent progresses within key areas that can provide customizable solutions to improve processes and the approaches taken within DDD is provided. Due to the broadness of the area, unfortunately all relevant aspects such as pharmacokinetics of bioactive molecules and delivery systems cannot be covered. Tailored approaches within (i) bioinformatics and computer-aided drug design, (ii) nanotechnology, (iii) novel materials and technologies for drug delivery and diagnostic systems, and (iv) disease models to predict safety and efficacy of medicines under development are focused on. Current developments and challenges ahead are discussed. The broad scope reflects the multidisciplinary nature of the field of DDD and aims to highlight the convergence of biological, pharmaceutical, and medical disciplines needed to meet the societal challenges of the 21st century.

Ras GTPases act as "molecular switches", alternating between inactive GDP-bound and active GTP-bound conformation. Ras-oncogenes were discovered over three decades ago, but there are still no effective therapies for Ras-driven cancers. So far, drug discovery strategies have been unsuccessful, because of a lack of suitable screening methodologies and well-defined binding pockets on the Ras proteins. Here, we addressed the former by introducing a homogeneous quenching resonance energy transfer (QRET) technique-based screening strategy for Ras interfacial and competitive inhibitors. We demonstrate that using a unique GTP-specific antibody fragment to monitor GTPase cycling in the presence of a guanine nucleotide exchange factor (GEF) and a GTPase activating protein (GAP) is an efficient method for Ras inhibitor high-throughput screening. When compared to a conventional GEF-stimulated nucleotide exchange assay in a proof-of-concept screen, we identified an overlapping set of potential inhibitor compounds but also compounds found exclusively with the new GTP hydrolysis monitoring-based GTPase cycling assay.

Summary Species tropism constitutes a serious problem for developing relevant animal models of infection. Human pathogens can express virulence factors that show specific selectivity to human proteins, while their affinity for orthologs from other species can vary significantly. Suitable animal species must be used to analyse whether virulence factors are potential targets for drug development. We developed an assay that rapidly predicts applicable animal species for studying virulence factors binding plasma proteins. We used two well‐characterized Staphylococcus aureus proteins, SSL7 and Efb, to develop an ELISA‐based inhibition assay using plasma from different animal species. The interaction between SSL7 and human C5 and the binding of Efb to human fibrinogen and human C3 was studied. Affinity experiments and Western blot analyses were used to validate the assay. Human, monkey and cat plasma interfered with binding of SSL7 to human C5. Binding of Efb to human fibrinogen was blocked in human, monkey, gerbil and pig plasma, while human, monkey, gerbil, rabbit, cat and guinea pig plasma inhibited the binding of Efb to human C3. These results emphasize the importance of choosing correct animal models, and thus, our approach is a rapid and cost‐effective method that can be used to prevent unnecessary animal experiments.

The SARS-CoV-2 was confirmed to cause the global pandemic of coronavirus disease 2019 (COVID-19). The 3-chymotrypsin-like protease (3CLpro), an essential enzyme for viral replication, is a valid target to combat SARS-CoV and MERS-CoV. In this work, we present a structure-based study to identify potential covalent inhibitors containing a variety of chemical warheads. The targeted Asinex Focused Covalent (AFCL) library was screened based on different reaction types and potential covalent inhibitors were identified. In addition, we screened FDA-approved protease inhibitors to find candidates to be repurposed against SARS-CoV-2 3CLpro. A number of compounds with significant covalent docking scores were identified. These compounds were able to establish a covalent bond (C-S) with the reactive thiol group of Cys145 and to form favorable interactions with residues lining the substrate-binding site. Moreover, paritaprevir and simeprevir from FDA-approved protease inhibitors were identified as potential inhibitors of SARS-CoV-2 3CLpro. The mechanism and dynamic stability of binding between the identified compounds and SARS-CoV-2 3CLpro were characterized by molecular dynamics (MD) simulations. The identified compounds are potential inhibitors worthy of further development as COVID-19 drugs. Importantly, the identified FDA-approved anti-hepatitis-C virus (HCV) drugs paritaprevir and simeprevir could be ready for clinical trials to treat infected patients and help curb COVID-19. Communicated by Ramaswamy H. Sarma. ispartof: JOURNAL OF BIOMOLECULAR STRUCTURE & DYNAMICS vol:39 issue:13 pages:4936-4948 ispartof: location:England status: published

A new series of amino acid derivatives of quinolines was synthesized through the hydrolysis of amino acid methyl esters of quinoline carboxamides with alkali hydroxide. The compounds were purified on silica gel by column chromatography and further characterized by TLC, NMR and ESI-TOF mass spectrometry. All compounds were screened for in vitro antimicrobial activity against different bacterial strains using the microdilution method. Most of the synthesized amino acid-quinolines show more potent or equipotent inhibitory action against the tested bacteria than their correspond esters. In addition, many of them exhibit fluorescent properties and could possibly be utilized as fluorophores. Molecular docking and simulation studies of the compounds at putative bacterial target enzymes suggest that the antimicrobial potency of these synthesized analogues could be due to enzyme inhibition via their favorable binding at the fluoroquinolone binding site at the GyrA subunit of DNA gyrase and/or the ParC subunit of topoisomerase-IV.

Chemokine receptors are key regulators of cell migration in terms of immunity and inflammation. Among these, CCR5 and CXCR4 play pivotal roles in cancer metastasis and HIV-1 transmission and infection. They act as essential co-receptors for HIV and furnish a route to the cell entry. In particular, inhibition of either CCR5 or CXCR4 leads very often the virus to shift to a more virulent dual-tropic strain. Therefore, dual receptor inhibition might improve the therapeutic strategies against HIV. In this study, we aimed to discover selective CCR5, CXCR4, and dual CCR5/CXCR4 antagonists using both receptor- and ligand-based computational methods. We employed this approach to fully incorporate the interaction attributes of the binding pocket together with molecular dynamics (MD) simulations and binding free energy calculations. The best hits were evaluated for their anti-HIV-1 activity against CXCR4- and CCR5-specific NL4.3 and BaL strains. Moreover, the Ca2+ mobilization assay was used to evaluate their antagonistic activity. From the 27 tested compounds, three were identified as inhibitors: compounds 27 (CCR5), 6 (CXCR4) and 3 (dual) with IC50 values ranging from 10.64 to 64.56 μM. The binding mode analysis suggests that the active compounds form a salt bridge with the glutamates and π-stacking interactions with the aromatic side chains binding site residues of the respective co-receptor. The presented hierarchical virtual screening approach provides essential aspects in identifying potential antagonists in terms of selectivity against a specific co-receptor. The compounds having multiple heterocyclic nitrogen atoms proved to be relatively more specific towards CXCR4 inhibition as compared to CCR5. The identified compounds serve as a starting point for further development of HIV entry inhibitors through synthesis and quantitative structure-activity relationship studies., Graphical abstract Image, graphical abstract

Human thymidylate synthase (hTS), a target for antiproliferative drugs, is an obligate homodimer. Single-point mutations to alanine at the monomer–monomer interface may enable the identification of specific residues that delineate sites for drugs aimed at perturbing the protein–protein interactions critical for activity. We computationally identified putative hotspot residues at the interface and designed mutants to perturb the intersubunit interaction. Dimer dissociation constants measured by a FRET-based assay range from 60 nM for wild-type hTS up to about 1 mM for single-point mutants and agree with computational predictions of the effects of these mutations. Mutations that are remote from the active site retain full or partial activity, although the substrate KM values were generally higher and the dimer was less stable. The lower dimer stability of the mutants can facilitate access to the dimer interface by small molecules and thereby aid the design of inhibitors that bind at the dimer interface.Graphical abstract

Sterol 14α-demethylase (cytochrome P450 family 51 (CYP51)) is an essential enzyme occurring in all biological kingdoms. In eukaryotes, it is located in the membrane of the endoplasmic reticulum. Selective inhibitors of trypanosomal CYP51s that do not affect the human CYP51 have been discovered in vitro and found to cure acute and chronic mouse Chagas disease without severe side effects in vivo. Crystal structures indicate that CYP51 may be more rigid than most CYPs, and it has been proposed that this property may facilitate antiparasitic drug design. Therefore, to investigate the dynamics of trypanosomal CYP51, we built a model of membrane-bound Trypanosoma brucei CYP51 and then performed molecular dynamics simulations of T. brucei CYP51 in membrane-bound and soluble forms. We compared the dynamics of T. brucei CYP51 with those of human CYP51, CYP2C9, and CYP2E1. In the simulations, the CYP51s display low mobility in the buried active site although overall mobility is similar in all the CYPs studied. The simulations suggest that in CYP51, pathway 2f serves as the major ligand access tunnel, and both pathways 2f (leading to membrane) and S (leading to solvent) can serve as ligand egress tunnels. Compared with the other CYPs, the residues at the entrance of the ligand access tunnels in CYP51 have higher mobility that may be necessary to facilitate the passage of its large sterol ligands. The water (W) tunnel is accessible to solvent during most of the simulations of CYP51, but its width is affected by the conformations of the heme's two propionate groups. These differ from those observed in the other CYPs studied because of differences in their hydrogen-bonding network. Our simulations give insights into the dynamics of CYP51 that complement the available experimental data and have implications for drug design against CYP51 enzymes.

The recent outbreak of Zika virus (ZIKV) infection in Brazil has developed to a global health concern due to its likely association with birth defects (primary microcephaly) and neurological complications. Consequently, there is an urgent need to develop a vaccine to prevent or a medicine to treat the infection. In this study, immunoinformatics approach was employed to predict antigenic epitopes of Zika viral proteins to aid in development of a peptide vaccine against ZIKV. Both linear and conformational B-cell epitopes as well as cytotoxic T-lymphocyte (CTL) epitopes were predicted for ZIKV Envelope (E), NS3 and NS5 proteins. We further investigated the binding interactions of altogether 15 antigenic CTL epitopes with three class I major histocompatibility complex (MHC I) proteins after docking the peptides to the binding groove of the MHC I proteins. The stability of the resulting peptide-MHC I complexes was further studied by molecular dynamics simulations. The simulation results highlight the limits of rigid-body docking methods. Some of the antigenic epitopes predicted and analyzed in this work might present a preliminary set of peptides for future vaccine development against ZIKV.

Human thymidylate synthase (hTS) is an established anticancer target. It catalyses the production of 2′-deoxythymidine-5′-monophosphate, an essential building block for DNA synthesis. Because of the development of cellular drug resistance against current hTS inhibitors, alternative inhibition strategies are needed. hTS exists in two forms, active and inactive, defined by the conformation of the active-site (AS) loop, which carries the catalytic cysteine, C195. To investigate the mechanism of activation and inactivation, targeted molecular dynamics (TMD) simulations of the transitions between active and inactive states of hTS were performed. Analysis of changes in the dihedral angles in the AS loop during different TMD simulations revealed complex conformational transitions. Despite hTS being a homodimeric enzyme and the conformational transition significantly involving the dimer interface, the transition occurs in an asymmetric, sequential manner via an ensemble of pathways. In addition to C195, which reoriented during the simulations, other key residues in the rotation of the AS loop included W182 and R185. The interactions of the cognate bulky W182 residues at the dimer interface hindered the simultaneous twist of the AS loops in the hTS dimer. Interactions of R185, which is unique for hTS, with ligands at different allosteric sites affected the activation transition. In addition to providing insights into the activation/inactivation mechanism of hTS and how conformational transitions can occur in homodimeric proteins, our observations suggest that blocking of AS loop rotation by ligands binding in the large cavity between the loops could be one way to stabilize inactive hTS and inhibit the enzyme. Proteins 2011; © 2011 Wiley-Liss, Inc.

This Special Issue contains thirteen articles that provide a vivid snapshot of the state-of-the-art of molecular modeling in drug design, illustrating recent advances and critically discussing important challenges [...]

Advances in understanding the physiology and pharmacology of the endogenous cannabinoid system have potentiated the interest of cannabinoid receptors as potential therapeutic targets. Cannabinoids have been shown to modulate a variety of immune cell functions and have therapeutic implications on central nervous system (CNS) inflammation, chronic inflammatory conditions such as arthritis, and may be therapeutically useful in treating autoimmune conditions such as multiple sclerosis. Many of these drug effects occur through cannabinoid receptor signalling mechanisms and the modulation of cytokines and other gene products. Further, endocannabinoids have been found to have many physiological and patho-physiological functions, including mood alteration and analgesia, control of energy balance, gut motility, motor and co-ordination activities, as well as alleviation of neurological, psychiatric and eating disorders. Plants offer a wide range of chemical diversity and have been a growing domain in the search for effective cannabinoid ligands. Cannabis sativa L. with the known plant cannabinoid, Delta(9-)tetrahydrocannabinol (THC) and Echinacea species with the cannabinoid (CB) receptor-binding lipophilic alkamides are the best known herbal cannabimimetics. This review focuses on the state of the art in CB ligands from plants, as well their possible therapeutic and immunomodulatory effects.

Sialic acid-binding immunoglobulin-like lectin-9 (Siglec-9) on leukocyte surface is a counter-receptor for endothelial cell surface adhesin, human primary amine oxidase (hAOC3), a target protein for anti-inflammatory agents. This interaction can be used to detect inflammation and cancer in vivo, since the labeled peptides derived from the second C2 domain (C22) of Siglec-9 specifically bind to the inflammation-inducible hAOC3. As limited knowledge on the interaction between Siglec-9 and hAOC3 has hampered both hAOC3-targeted drug design and in vivo imaging applications, we have now produced and purified the extracellular region of Siglec-9 (Siglec-9-EC) consisting of the V, C21 and C22 domains, modeled its 3D structure and characterized the hAOC3-Siglec-9 interactions using biophysical methods and activity/inhibition assays. Our results assign individual, previously unknown roles for the V and C22 domains. The V domain is responsible for the unusually tight Siglec-9-hAOC3 interactions whereas the intact C22 domain of Siglec-9 is required for modulating the enzymatic activity of hAOC3, crucial for the hAOC3-mediated leukocyte trafficking. By characterizing the Siglec-9-EC mutants, we could conclude that R120 in the V domain likely interacts with the terminal sialic acids of hAOC3 attached glycans whereas residues R284 and R290 in C22 are involved in the interactions with the active site channel of hAOC3. Furthermore, the C22 domain binding enhances the enzymatic activity of hAOC3 although the sialic acid-binding capacity of the V domain of Siglec-9 is abolished by the R120S mutation. To conclude, our results prove that the V and C22 domains of Siglec-9-EC interact with hAOC3 in a multifaceted and unique way, forming both glycan-mediated and direct protein-protein interactions, respectively. The reported results on the mechanism of the Siglec-9-hAOC3 interaction are valuable for the development of hAOC3-targeted therapeutics and diagnostic tools.

The recent expiration of several blockbuster pharmaceutical patents offers new opportunities for generic drug production in Africa. Moreover, 2015 marked a critical juncture; a transition from the Millennium Development Goals to Sustainable Development Goals. The implications for African economies in the area of generic drug production and global health outcomes are vast given the potential to increase access to medicines for neglected diseases and other emerging health crises. This issue-based article analyses the extent to which several coeval variables of governance and macro-economic nature can potentially create the market and institutional conditions to spur innovations for improving access to medicine via cross-sector social partnerships. Proposals for solving grand challenges in Africa’s pharmaceutical markets often fail to address the most fundamental impediments to innovation, besides being mostly donor-driven. Through document analysis, we problematize conventional formulae for healthcare governance with a measured critique of prevailing orthodoxies by offering implementable alternatives. We propose inclusive, innovative models for marshalling sustainable access to high-quality affordable medicine. We identify bottom-up and entrepreneurially viable strategic reversal of decades of systematic damages that have contributed to the underdeveloped pharmaceutical market, whilst striking a reasonable balance between what the desirable future is and what can set the stage for a durable change through game-changing market and innovative mechanisms.

Human thymidylate synthase is a homodimeric enzyme that plays a key role in DNA synthesis and is a target for several clinically important anticancer drugs that bind to its active site. We have designed peptides to specifically target its dimer interface. Here we show through X-ray diffraction, spectroscopic, kinetic, and calorimetric evidence that the peptides do indeed bind at the interface of the dimeric protein and stabilize its di-inactive form. The “LR” peptide binds at a previously unknown binding site and shows a previously undescribed mechanism for the allosteric inhibition of a homodimeric enzyme. It inhibits the intracellular enzyme in ovarian cancer cells and reduces cellular growth at low micromolar concentrations in both cisplatin-sensitive and -resistant cells without causing protein overexpression. This peptide demonstrates the potential of allosteric inhibition of hTS for overcoming platinum drug resistance in ovarian cancer.

Many enzymes and proteins are regulated by their quaternary structure and/or by their association in homo- and/or hetero-oligomer complexes. Thus, these protein-protein interactions can be good targets for blocking or modulating protein function therapeutically. The large number of oligomeric structures in the Protein Data Bank (http://www.rcsb.org/) reflects growing interest in proteins that function as multimeric complexes. In this review, we consider the particular case of homodimeric enzymes as drug targets. There is intense interest in drugs that inhibit dimerization of a functionally obligate homodimeric enzyme. Because amino acid conservation within enzyme interfaces is often low compared to conservation in active sites, it may be easier to achieve drugs that target protein interfaces selectively and specifically. Two main types of dimerization inhibitors have been developed: peptides or peptidomimetics based on sequences involved in protein-protein interactions, and small molecules that act at hot spots in protein-protein interfaces. Examples include inhibitors of HIV protease and HIV integrase. Studying the mechanisms of action and locating the binding sites of such inhibitors requires different techniques for different proteins. For some enzymes, ligand binding is only detectable in vivo or after unfolding of the complexes. Here, we review the structural features of dimeric enzymes and give examples of inhibition through interference in dimer stability. Several techniques for studying these complex phenomena will be presented.

Owing to its central role in DNA synthesis, human thymi-dylate synthase (hTS) is a well-established target for che-motherapeutic agents, such as fluoropyrimidines. The useof hTS inhibitors in cancer therapy is limited by theirtoxicity and the development of cellular drug resistance.Here, with the aim of shedding light on the structuralrole of the A-helix in fluoropyrimidine resistance, we havecreated a fluoropyrimidine-resistant mutant by making asingle point mutation, Glu30Trp. We postulated thatresidue 30, which is located in the A-helix, close to butoutside the enzyme active site, could have a long-rangeeffect on inhibitor binding. The mutant shows 100 timeslower specific activity with respect to the wild-type hTSand is resistant to the classical inhibitor, FdUMP, asshown by a 6-fold higher inhibition constant. Circulardichroism experiments show that the mutant is folded.The results of molecular modeling and simulation suggestthat the Glu30Trp mutation gives rise to resistance byaltering the hydrogen-bond network between residue 30and the active site.Keywords: antitumor agents/drug resistance/humanthymidylate synthase/protein design/protein dynamicsIntroductionHuman thymidylate synthase (hTS) (EC 2.1.1.45) catalyzesthe reaction in which dUMP (2

Given the three-dimensional structure of a protein, how can one find the sites where other molecules might bind to it? Do these sites have the properties necessary for high affinity binding? Is this protein a suitable target for drug design? Here, we discuss recent developments in computational methods to address these and related questions. Geometric methods to identify pockets on protein surfaces have been developed over many years but, with new algorithms, their performance is still improving. Simulation methods show promise in accounting for protein conformational variability to identify transient pockets but lack the ease of use of many of the (rigid) shape-based tools. Sequence and structure comparison approaches are benefiting from the constantly increasing size of sequence and structure databases. Energetic methods can aid identification and characterization of binding pockets, and have undergone recent improvements in the treatment of solvation and hydrophobicity. The "druggability" of a binding site is still difficult to predict with an automated procedure. The methodologies available for this purpose range from simple shape and hydrophobicity scores to computationally demanding free energy simulations.

The endocannabinoid system consists of two cannabinoid receptors (CB1 and CB2), endogenous ligands (endocannabinoids), and the enzymes involved in the metabolism of the endocannabinoids, including fatty acid amide hydrolase (FAAH) and monoglyceride lipase (MGL). In the present study, virtual screening of MGL inhibitors was performed by utilizing a comparative model of the human MGL enzyme. All hit molecules were tested for their potential MGL inhibitory activity, but no compounds were found capable of inhibiting MGL-like enzymatic activity in rat cerebellar membranes. However, these compounds were also tested for their potential FAAH inhibitory activity and five compounds (2-6) inhibiting FAAH were found with IC50 values between 4 and 44 microM. In addition, the hit molecules from the virtual screening of CB2 receptor ligands (reported previously in Salo et al. J. Med. Chem. 2005, 48, 7166) were also tested in our FAAH assay, and four active compounds (7-10) were found with IC50 values between 0.52 and 22 microM. Additionally, compound 7 inhibited MGL-like enzymatic activity with an IC50 value of 31 microM.