Search

Non-steroidal anti-inflammatory drugs (NSAIDs) are a powerful class of inhibitors targeting two isoforms of the family of cyclooxygenase enzymes (COX-1 and COX-2). While NSAIDs are widely used in the management of pain, in particular as a treatment for osteo- and rheumatoid arthritis, their long-term use has been associated with numerous on- and off-target effects. As the carboxylic acid moiety present in common NSAIDs is responsible for some of their adverse effects, but is not required for their anti-inflammatory activity, we sought to mask this group through direct coupling to glucosamine, which is thought to prevent cartilage degradation. We report herein the conjugation of commonly prescribed NSAIDs to glucosamine hydrochloride and the use of molecular docking to show that addition of the carbohydrate moiety to the parent NSAID can enhance binding in the active site of COX-2. In a preliminary, in vitro screening assay, the diclofenac-glucosamine bioconjugate exhibited 10-fold greater activity toward COX-2, making it an ideal candidate for future in vivo studies. Furthermore, in an intriguing result, we observed that the mefenamic acid-glucosamine bioconjugate displayed enhanced activity toward COX-1 rather than COX-2.

To maintain the lipid asymmetry of the cell envelope in Gram-negative bacteria, the MlaC protein serves as a lipid transfer factor and delivers phospholipids from the outer to the inner membrane. A strategy of antibiotic discovery is to design a proper compound that can tightly bind to the MlaC protein and inhibit the MlaC function. In this study, we performed virtual screening on multiple MlaC structures obtained from molecular dynamics simulations to identify potential MlaC binders. Our results suggested that clorobiocin is a compound that could bind to the MlaC protein. Through the comparison of the bound geometry between clorobiocin and novobiocin, we pointed out that the methyl-pyrrole group of the noviose sugar in clorobiocin forms hydrophobic interactions with amino acids in the phospholipid binding pocket, which allows the compound to bind deep in the active site. This also explains why clorobiocin shows a tighter binding affinity than novobiocin. Our study highlights a practical path of antibiotic development against Gram-negative bacteria.

Malaria, mainly caused by Plasmodium falciparum and Plasmodium vivax, has been a growing cause of morbidity and mortality. P. falciparum is more lethal than is P. vivax, but there is a vital need for effective drugs against both species. Geranylgeranyl diphosphate synthase (GGPPS) is an enzyme involved in the biosynthesis of quinones and in protein prenylation and has been proposed to be a malaria drug target. However, the structure of P. falciparumGGPPS (PfGGPPS) has not been determined, due to difficulties in crystallization. Here, we created a PfGGPPS model using the homologous P.vivaxGGPPS X-ray structure as a template. We simulated the modeled PfGGPPS as well as PvGGPPS using conventional and Gaussian accelerated molecular dynamics in both apo- and GGPP-bound states. The MD simulations revealed a striking similarity in the dynamics of both enzymes with loop 9-10 controlling access to the active site. We also found that GGPP stabilizes PfGGPPS and PvGGPPS into closed conformations and via similar mechanisms. Shape-based analysis of the binding sites throughout the simulations suggests that the two enzymes will be readily targeted by the same inhibitors. Finally, we produced three MD-validated conformations of PfGGPPS to be used in future virtual screenings for potential new antimalarial drugs acting on both PvGGPPS and PfGGPPS.

The A2A adenosine receptor (A2A AR) is a G protein-coupled receptor that is pharmacologically targeted for the treatment of inflammation, sepsis, cancer, neurodegeneration, and Parkinson's disease. Recently, we applied long-timescale molecular dynamics simulations on two ligand-free receptor conformations, starting from the agonist-bound (PDB ID: 3QAK) and antagonist-bound (PDB ID: 3EML) X-ray structures. This analysis revealed four distinct conformers of the A2A AR: the active, intermediate 1, intermediate 2, and inactive. In this study, we apply the fragment-based mapping algorithm, FTMap, on these receptor conformations to uncover five non-orthosteric sites on the A2A AR. Two sites that are identified in the active conformation are located in the intracellular region of the transmembrane helices (TM) 3/TM4 and the G protein-binding site in the intracellular region between TM2/TM3/TM6/TM7. Three sites are identified in the intermediate 1 and intermediate 2 conformations, annexing a site in the lipid interface of TM5/TM6. Five sites are identified in the inactive conformation, comprising a site in the intracellular region of TM1/TM7 and in the extracellular region of TM3/TM4 of the A2A AR. We postulate that these sites on the A2A AR be screened for allosteric modulators for the treatment of inflammatory and neurological diseases.

Staphylococcus aureus is a leading cause of hospital-acquired infections in the USA and is a major health concern as methicillin-resistant S. aureus and other antibiotic-resistant strains are common. Compounds that inhibit the S. aureus sortase (SrtA) cysteine transpeptidase may function as potent anti-infective agents as this enzyme attaches virulence factors to the bacterial cell wall. While a variety of SrtA inhibitors have been discovered, the vast majority of these small molecules have not been optimized using structure-based approaches. Here we have used NMR spectroscopy to determine the molecular basis through which pyridazinone-based small molecules inhibit SrtA. These inhibitors covalently modify the active cysteine thiol and partially mimic the natural substrate of SrtA by inducing the closure of an active site loop. Computational and synthetic chemistry methods led to second-generation analogues that are ~70-fold more potent than the lead molecule. These optimized molecules exhibit broad-spectrum activity against other types of class A sortases, have reduced cytotoxicity, and impair SrtA-mediated protein display on S. aureus cell surface. Our work shows that pyridazinone analogues are attractive candidates for further development into anti-infective agents, and highlights the utility of employing NMR spectroscopy and solubility-optimized small molecules in structure-based drug discovery.

Intrinsically disordered proteins are proteins which lack of specific tertiary structure and unable to fold spontaneously without the partner binding. These intrinsically disordered proteins are found to associate with various diseases, such as diabetes, cancer, and neurodegenerative diseases. However, current widely used force fields, such as ff99SB, ff14SB, OPLS/AA, and Charmm27, are insufficient in sampling the conformational characters of intrinsically disordered proteins. In this study, the CMAP method was used to correct the φ/ψ distributions of disorder-promoting amino acids. The simulation results show that the force filed parameters (ff14IDPs) can improve the φ/ψ distributions of the disorder-promoting amino acids, with RMSD less than 0.10% relative to the benchmark data of intrinsically disordered proteins. Further test suggests that the calculated secondary chemical shifts under ff14IDPs are in quantitative agreement with the data of NMR experiment for five tested systems. In addition, the simulation results show that ff14IDPs can still be used to model structural proteins, such as tested lysozyme and ubiquitin, with better performance in coil regions than the original general Amber force field ff14SB. These findings confirm that the newly developed Amber ff14IDPs is a robust model for improving the conformation sampling of intrinsically disordered proteins.

Continuing our search to find more potent and selective 5-LOX inhibitors, we present now the enzymatic evaluation of seventeen isoflavones (IR) and nine isoflavans (HIR), and their in vitro and in cellulo potency against human leukocyte 5-LOX. Of the 26 compounds tested, 10 isoflavones and 9 isoflavans possessed micromolar potency, but only three were selective against 5-LOX (IR-2, HIR-303, and HIR-309), with IC50 values at least 10 times lower than those of 12-LOX, 15-LOX-1, and 15-LOX-2. Of these three, IR-2 (6,7-dihydroxy-4-methoxy-isoflavone, known as texasin) was the most selective 5-LOX inhibitor, with over 80-fold potency difference compared with other isozymes; Steered Molecular Dynamics (SMD) studies supported these findings. The presence of the catechol group on ring A (6,7-dihydroxy versus 7,8-dihydroxy) correlated with their biological activity, but the reduction of ring C, converting the isoflavones to isoflavans, and the substituent positions on ring B did not affect their potency against 5-LOX. Two of the most potent/selective inhibitors (HIR-303 and HIR-309) were reductive inhibitors and were potent against 5-LOX in human whole blood, indicating that isoflavans can be potent and selective inhibitors against human leukocyte 5-LOX in vitro and in cellulo.

With the rise in antibiotic resistance, there is interest in discovering new drugs active against new targets. Here, we investigate the dynamic structures of three isoprenoid synthases from Mycobacterium tuberculosis using molecular dynamics (MD) methods with a view to discovering new drug leads. Two of the enzymes, cis-farnesyl diphosphate synthase (cis-FPPS) and cis-decaprenyl diphosphate synthase (cis-DPPS), are involved in bacterial cell wall biosynthesis, while the third, tuberculosinyl adenosine synthase (Rv3378c), is involved in virulence factor formation. The MD results for these three enzymes were then compared with previous results on undecaprenyl diphosphate synthase (UPPS) by means of active site volume fluctuation and principal component analyses. In addition, an analysis of the binding of prenyl diphosphates to cis-FPPS, cis-DPPS, and UPPS utilizing the new MD results is reported. We also screened libraries of inhibitors against cis-DPPS, finding ~1 μm inhibitors, and used the receiver operating characteristic-area under the curve (ROC-AUC) method to test the predictive power of X-ray and MD-derived cis-DPPS receptors. We found that one compound with potent M. tuberculosis cell growth inhibition activity was an IC(50) ~0.5- to 20-μm inhibitor (depending on substrate) of cis-DPPS, a ~660-nm inhibitor of Rv3378c as well as a 4.8-μm inhibitor of cis-FPPS, opening up the possibility of multitarget inhibition involving both cell wall biosynthesis and virulence factor formation.

In the fight against heart failure, therapeutics that have the ability to increase the contractile power of the heart are urgently needed. One possible route of action to improve heart contractile power is increasing the calcium sensitivity of the thin filament. From a pharmaceutical standpoint, calcium sensitizers have the distinct advantage of not altering cardiomyocyte calcium levels and thus have lower potential for side-effects. Small chemical molecules have been shown to bind to the interface between cTnC and the cTnI switch peptide and exhibit calcium-sensitizing properties, possibly by stabilizing cTnC in an open conformation. Building on existing structural data of a known calcium sensitizer bound to cardiac troponin, we combined computational structure-based virtual screening drug discovery methods and solution NMR titration assays to identify a novel calcium sensitizer 4-(4-(2,5-dimethylphenyl)-1-piperazinyl)-3-pyridinamine (NSC147866) which binds to cTnC and the cTnC-cTnI147-163 complex. Its presence increases the affinity of switch peptide to cTnC by approximately a factor of two. This action is comparable to that of known levosimendan analogues.

Non-nucleoside reverse transcriptase inhibitors (NNRTIs) are potent anti-HIV chemotherapeutics. Although there are FDA-approved NNRTIs, challenges such as the development of resistance have limited their utility. Here, we describe the identification of novel NNRTIs through a combination of computational and experimental approaches. Based on the known plasticity of the NNRTI binding pocket (NNIBP), we adopted an ensemble-based virtual screening strategy: coupling receptor conformations from 10 X-ray crystal structures with 120 snapshots from a total of 480 ns of molecular dynamics (MD) trajectories. A screening library of 2864 National Cancer Institute (NCI) compounds was built and docked against the ensembles in a hierarchical fashion. Sixteen diverse compounds were tested for their ability to block HIV infection in human tissue cultures using a luciferase-based reporter assay. Three promising compounds were further characterized, using a HIV-1 RT-based polymerase assay, to determine the specific mechanism of inhibition. We found that 2 of the three compounds inhibited the polymerase activity of RT (with potency similar to the positive control, the FDA-approved drug nevirapine). Through a computational approach, we were able to discover two compounds which inhibit HIV replication and block the activity of RT, thus offering the potential for optimization into mature inhibitors.

Proteins with expanded polyglutamine (polyQ) segments cause a number of fatal neurodegenerative disorders, including Huntington's disease (HD). Previous high-throughput screens in cellular and biochemical models of HD have revealed compounds that mitigate polyQ aggregation and proteotoxicity, providing insight into the mechanisms of disease and leads for potential therapeutics. However, the structural diversity of natural products has not yet been fully mobilized toward these goals. Here, we have screened a collection of ~11 000 natural product extracts for the ability to recover the slow growth of ΔProQ103-expressing yeast cells in 384-well plates (Z' ~ 0.7, CV ~ 8%). This screen identified actinomycin D as a strong inhibitor of polyQ aggregation and proteotoxicity at nanomolar concentrations (~50-500 ng/mL). We found that a low dose of actinomycin D increased the levels of the heat-shock proteins Hsp104, Hsp70 and Hsp26 and enhanced binding of Hsp70 to the polyQ in yeast. Actinomycin also suppressed aggregation of polyQ in mammalian cells, suggesting a conserved mechanism. These results establish natural products as a rich source of compounds with interesting mechanisms of action against polyQ disorders.

G-protein-coupled receptors (GPCRs) are key cellular signaling proteins and have been targeted by approximately 30-40% of marketed drugs for treating many human diseases including cancer and heart failure. Recently, we directly observed activation of the M2 muscarinic receptor through long-timescale accelerated molecular dynamics (aMD) simulation, which revealed distinct inactive, intermediate and active conformers of the receptor. Here, FTMAP is applied to search for 'hot spots' in these activation-associated conformers using a library of 16 organic probe molecules that represent fragments of potential drugs. Seven allosteric (non-orthosteric) binding sites are identified in the M2 receptor through the FTMAP analysis. These sites are distributed in the solvent-exposed extracellular and intracellular mouth regions, as well as the lipid-exposed pockets formed by the transmembrane α helices TM3-TM4, TM5-TM6 and TM7-TM1/TM2. They serve as promising target sites for designing novel allosteric modulators as receptor-selective drugs.

Staphylococcus aureus is the leading cause of hospital-acquired infections in the United States. The emergence of multidrug-resistant strains of S. aureus has created an urgent need for new antibiotics. Staphylococcus aureus uses the sortase A enzyme to display surface virulence factors suggesting that compounds that inhibit its activity will function as potent anti-infective agents. Here, we report the identification of several inhibitors of sortase A using virtual screening methods that employ the relaxed complex scheme, an advanced computer-docking methodology that accounts for protein receptor flexibility. Experimental testing validates that several compounds identified in the screen inhibit the activity of sortase A. A lead compound based on the 2-phenyl-2,3-dihydro-1H-perimidine scaffold is particularly promising, and its binding mechanism was further investigated using molecular dynamics simulations and conducting preliminary structure-activity relationship studies.

Following sequencing and assembly of the human genome, the preferred methods for identification of new drug targets have changed dramatically. Modern tactics such as genome-wide association studies (GWAS) and deep sequencing are fundamentally different from the pharmacology-guided approaches used previously, in which knowledge of small molecule ligands acting at their cellular targets was the primary discovery engine. A consequence of the 'target-first, pharmacology-second' strategy is that many predicted drug targets are non-enzymes, such as scaffolding, regulatory or structural proteins, and their activities are often dependent on protein-protein interactions (PPIs). These types of targets create unique challenges to drug discovery efforts because enzymatic turnover cannot be used as a convenient surrogate for compound potency. Moreover, it is often challenging to predict how ligand binding to non-enzymes might affect changes in protein function and/or pathobiology. Thus, in the postgenomic era, targets might be strongly implicated by molecular biology-based methods, yet they often later earn the designation of 'undruggable'. Can the scope of available targets be widened to include these promising, but challenging, non-enzymes? In this review, we discuss advances in high-throughput screening (HTS) technology and chemical library design that are emerging to deal with these challenges.

The protozoan parasite Trypanosoma cruzi, the etiological agent of Chagas' disease, affects millions of individuals and continues to be an important global health concern. The poor efficacy and unfavorable side effects of current treatments necessitate novel therapeutics. Cruzain, the major cysteine protease of T. cruzi, is one potential novel target. Recent advances in a class of vinyl sulfone inhibitors are encouraging; however, as most potential therapeutics fail in clinical trials and both disease progression and resistance call for combination therapy with several drugs, the identification of additional classes of inhibitory molecules is essential. Using an exhaustive virtual-screening and experimental validation approach, we identify several additional small-molecule cruzain inhibitors. Further optimization of these chemical scaffolds could lead to the development of novel drugs useful in the treatment of Chagas' disease.

The relaxed complex scheme, a virtual-screening methodology that accounts for protein receptor flexibility, was used to identify a low-micromolar, non-bisphosphonate inhibitor of farnesyl diphosphate synthase. Serendipitously, we also found that several predicted farnesyl diphosphate synthase inhibitors were low-micromolar inhibitors of undecaprenyl diphosphate synthase. These results are of interest because farnesyl diphosphate synthase inhibitors are being pursued as both anti-infective and anticancer agents, and undecaprenyl diphosphate synthase inhibitors are antibacterial drug leads.

Matrix metalloproteinases are zinc-containing enzymes capable of degrading all components of the extracellular matrix. Owing to their role in human disease, matrix metalloproteinase have been the subject of extensive study. A bioinorganic approach was recently used to identify novel inhibitors based on a maltol zinc-binding group, but accompanying molecular-docking studies failed to explain why one of these inhibitors, AM-6, had approximately 2500-fold selectivity for MMP-3 over MMP-2. A number of studies have suggested that the matrix-metalloproteinase active site is highly flexible, leading some to speculate that differences in active-site flexibility may explain inhibitor selectivity. To extend the bioinorganic approach in a way that accounts for MMP-2 and MMP-3 dynamics, we here investigate the predicted binding modes and energies of AM-6 docked into multiple structures extracted from matrix-metalloproteinase molecular dynamics simulations. Our findings suggest that accounting for protein dynamics is essential for the accurate prediction of binding affinity and selectivity. Additionally, AM-6 and other similar inhibitors likely select for and stabilize only a subpopulation of all matrix-metalloproteinase conformations sampled by the apo protein. Consequently, when attempting to predict ligand affinity and selectivity using an ensemble of protein structures, it may be wise to disregard protein conformations that cannot accommodate the ligand.

In most eubacteria, apicomplexans, and most plants, including the causal agents for diseases such as malaria, leprosy, and tuberculosis, the methylerythritol phosphate pathway is the route for the biosynthesis of the C(5) precursors to the essential isoprenoid class of compounds. Owing to their absence in humans, the enzymes of the methylerythritol phosphate pathway have become attractive targets for drug discovery. This work investigates a new class of inhibitors against the second enzyme of the pathway, 1-deoxy-D-xylulose 5-phosphate reductoisomerase. Inhibition of this enzyme may involve the chelation of a crucial active site Mn ion, and the metal-chelating moieties studied here have previously been shown to be successful in application to the zinc-dependent metalloproteinases. Quantum mechanics and docking calculations presented in this work suggest the transferability of these metal-chelating compounds to Mn-containing 1-deoxy-D-xylulose 5-phosphate reductoisomerase enzyme, as a promising starting point to the development of potent inhibitors.

To address the problem of specificity in G-protein coupled receptor (GPCR) drug discovery, there has been tremendous recent interest in allosteric drugs that bind at sites topographically distinct from the orthosteric site. Unfortunately, structure-based drug design of allosteric GPCR ligands has been frustrated by the paucity of structural data for allosteric binding sites, making a strong case for predictive computational methods. In this work, we map the surfaces of the beta1 (beta1AR) and beta2 (beta2AR) adrenergic receptor structures to detect a series of five potentially druggable allosteric sites. We employ the FTMAP algorithm to identify 'hot spots' with affinity for a variety of organic probe molecules corresponding to drug fragments. Our work is distinguished by an ensemble-based approach, whereby we map diverse receptor conformations taken from molecular dynamics (MD) simulations totaling approximately 0.5 micros. Our results reveal distinct pockets formed at both solvent-exposed and lipid-exposed cavities, which we interpret in light of experimental data and which may constitute novel targets for GPCR drug discovery. This mapping data can now serve to drive a combination of fragment-based and virtual screening approaches for the discovery of small molecules that bind at these sites and which may offer highly selective therapies.