[Discovering new antiretroviral compounds in "Big Data" chemical space of the SAVI library]. / Poisk novykh antiretrovirusnykh soedinenii v khimicheskom prostranstve "bol'shikh dannykh" biblioteki SAVI.

Despite significant advances in the application of highly active antiretroviral therapy, the development of new drugs for the treatment of HIV infection remains an important task because the existing drugs do not provide a complete cure, cause serious side effects and lead to the emergence of resistance. In 2015, a consortium of American and European scientists and specialists launched a project to create the SAVI (Synthetically Accessible Virtual Inventory) library. Its 2016 version of over 283 million structures of new easily synthesizable organic molecules, each annotated with a proposed synthetic route, were generated in silico for the purpose of searching for safer and more potent pharmacological substances. We have developed an algorithm for comparing large chemical databases (DB) based on the representation of structural formulas in SMILES codes, and evaluated the possibility of detecting new antiretroviral compounds in the SAVI database. After analyzing the intersection of SAVI with 97 million structures of the PubChem database, we found that only a small part of the SAVI (~0.015%) is represented in PubChem, which indicates a significant novelty of this virtual library. However, among those structures, 632 compounds tested for anti-HIV activity were detected, 41 of which had the desired activity. Thus, our studies for the first time demonstrated that SAVI is a promising source for the search for new anti-HIV compounds.

Improving (Q)SAR predictions by examining bias in the selection of compounds for experimental testing.

Existing data on structures and biological activities are limited and distributed unevenly across distinct molecular targets and chemical compounds. The question arises if these data represent an unbiased sample of the general population of chemical-biological interactions. To answer this question, we analyzed ChEMBL data for 87,583 molecules tested against 919 protein targets using supervised and unsupervised approaches. Hierarchical clustering of the Murcko frameworks generated using Chemistry Development Toolkit showed that the available data form a big diffuse cloud without apparent structure. In contrast hereto, PASS-based classifiers allowed prediction whether the compound had been tested against the particular molecular target, despite whether it was active or not. Thus, one may conclude that the selection of chemical compounds for testing against specific targets is biased, probably due to the influence of prior knowledge. We assessed the possibility to improve (Q)SAR predictions using this fact: PASS prediction of the interaction with the particular target for compounds predicted as tested against the target has significantly higher accuracy than for those predicted as untested (average ROC AUC are about 0.87 and 0.75, respectively). Thus, considering the existing bias in the data of the training set may increase the performance of virtual screening.

Internet resources integrating many small-molecule databases.

()New data, tools and services recently made available on the web server (http://cactus.nci.nih.gov) of the Computer-Aided Drug Design (CADD) Group, NCI, NIH, developed in the context of chemoinformatics and drug development work, are presented. These tools are designed for searching for structures in very large databases of small molecules. One of them is a web service-the Chemical Structure Look-up Service (CSLS)-for very rapid structure look-up in an aggregated collection of more than 80 databases comprising more than 27 million unique structures at the time of this writing. CSLS contains pointers to the entries in toxicology-related databases, catalogues of commercially available samples, drugs, assay results data sets, and databases in several other categories. CSLS allows the user to find out very rapidly in which one(s) of all these databases a given structure occurs independent of the representation of the input structure, by making use of InChIs as well as new CACTVS hashcode-based identifiers. These latter, calculable, identifiers are designed to take into account tautomerism, different resonance structures drawn for charged species, and presence of additional fragments. They make possible fine-tunable yet rapid compound identification and database overlap analyses in very large compound collections.

Protein kinase C. Modeling of the binding site and prediction of binding constants.

A detailed examination of the mode of binding of phorbol esters to protein kinase C led to a model of the phorbol binding site in the enzyme. The efficacy with which various synthetic diacylglycerol analogs and ribonolactones are able to bind to this site was determined by means of semiempirical quantum mechanical calculations using PM3, and an estimate of the binding energy was made in each case. Sixteen synthetic analogs of 1,2-diacylglycerol and two natural products were studied, and their calculated energies of binding to this model were correlated with the measured Ki values. The binding energies calculated for this receptor model, together with solubility and entropy considerations, allow prediction through regressive fit of free energies of binding which correlate very well with the measured binding constants.

Conformationally locked nucleoside analogues. Synthesis of dideoxycarbocyclic nucleoside analogues structurally related to neplanocin C.

The glycon moiety of nucleosides in solution is known to exist in a rapid dynamic equilibrium between extreme northern and southern conformations as defined by the pseudorotation cycle. The concept of preparing rigid nucleoside analogues with the glycon conformation locked in one of these two extremes was tested with the synthesis of some cyclopropane-fused dideoxycarbocyclic nucleosides, similar to the well-known class of anti-HIV active dideoxynucleosides. The new compounds described here are dideoxynucleoside analogues of the fermentation product neplanocin C (6) which exhibits a typical northern geometry for its 6-oxabicyclo[3.1.0]hexane pseudosugar moiety. However, in view of the lability of the epoxide ring in this system, the equivalent cyclopropane-fused bicyclo[3.1.0]hexane system was used instead to prepare the corresponding dideoxynucleoside analogues bearing all the common bases [(+/-)-9-13]. Due to the well-documented preference of unrestricted bicyclo[3.1.0]hexane systems to exist exclusively in a boat conformation, the resulting nucleosides are structurally locked in a typical northern conformation similar to that of neplanocin C. The locked northern conformation in these nucleosides remained unchanged in solution in the 20-80 degrees C temperature range according to variable temperature 1H NMR studies. For the synthesis of these compounds, racemic trans-1-[(benzyloxy)methyl]-4-hydroxybicyclo[3.1.0]hexane [(+/-)-18] was prepared by a samarium-promoted cyclopropanation reaction with the antecedent cyclopentenol. All of the bases were incorporated under Mitsunobu conditions and converted to the desired final products following a standard methodology. Anti-HIV evaluation revealed that only the adenosine analogue (+/-)-9 possessed enough activity to warrant resolution into its optical antipodes. This was realized by chiral HPLC chromatography to give the individual enantiomers (-)-32 and (+)-33. Adenosine deaminase was used to identify isomer (+)-33 as the enantiomer with the "natural" configuration which was solely responsible for the observed biological activity and toxicity of (+/-)-9. It is possible that the exclusive northern conformation adopted by these nucleosides reduces their substrate affinity for the various activating kinases, except in the case of the adenosine analogue.

HIV-1 integrase pharmacophore: discovery of inhibitors through three-dimensional database searching.

Starting from a known inhibitor of human immunodeficiency virus type 1 (HIV-1) integrase (IN); caffeic acid phenethyl ester (CAPE), a putative three-point pharmacophore for binding of inhibitors to IN was derived. This pharmacophore was used to search the National Cancer Institute three-dimensional (3D) structural database. Out of the open, nonproprietary part of this database, comprising approximately 200000 compounds, 267 structures were found to match the pharmacophore in at least one conformation, and 60 of those were tested in an in vitro assay against HIV-1 IN. Out of these, 19 were found to inhibit both the 3'-processing and strand transfer of IN at micromolar concentrations. In order to test the validity of this pharmacophore, a small 3D database of 152 published IN inhibitors was built. A search in this database yielded a statistically significant correlation of the presence of this pharmacophore and the potency of the compounds. An automated pharmacophore identification procedure performed on this set of compounds provided additional support for the importance of this pharmacophore for binding of inhibitors to IN and hinted at a possible second pharmacophore. The role of aromatic moieties in the binding of ligands to HIV-1 IN through interactions with divalent metal cations, which are known to be necessary for activity of the enzyme, was explored in ab initio calculations.

Discovery of HIV-1 integrase inhibitors by pharmacophore searching.

Based upon a class of known HIV-1 integrase inhibitors, several pharmacophore models were proposed from molecular modeling studies and validated using a 3D database of 152, compounds for which integrase assay data are known. Using the most probable pharmacophore model as the query, the NCI 3D database of 206,876 compounds was searched, and 340 compounds that contain the pharmacophore query were identified. Twenty-nine of these compounds were selected and tested in the HIV-1 integrase assay. This led to the discovery of 10 novel, structurally diverse HIV-1 integrase inhibitors, four of which have an IC50 value less than 30 microM and are promising lead compounds for further HIV-1 integrase inhibitor development.

Depsides and depsidones as inhibitors of HIV-1 integrase: discovery of novel inhibitors through 3D database searching.

Seventeen lichen acids comprising despides, depsidones, and their synthetic derivatives have been examined for their inhibitory activity against HIV-1 integrase, and two pharmacophores associated with inhibition of this enzyme have been identified. A search of the NCI 3D database of approximately 200,000 structures yielded some 800 compounds which contain one or the other pharmacophore. Forty-two of these compounds were assayed for HIV-1 integrase inhibition, and of these, 27 had inhibitory IC50 values of less than 100 microM; 15 were below 50 microM. Several of these compounds were also examined for their activity against HIV-2 integrase and mammalian topoisomerase I.

An optimized protein kinase C activating diacylglycerol combining high binding affinity (Ki) with reduced lipophilicity (log P).

A small, focused combinatorial library encompassing all possible permutations of acyl branched alkyl chains-small and large, saturated and unsaturated-was generated from the active diacylglycerol enantiomer (S-DAG) to help identify the analogue with the highest binding affinity (lowest Ki) for protein kinase C (PK-C) combined with the minimum lipophilicity (log P). The selected ligand (3B) activated PK-C more effectively than sn-1,2-dioctanoylglycerol (diC8) despite being 1.4 log units more hydrophilic. Compound 3B indeed represents the most potent, hydrophilic DAG ligand to date. With the help of a green fluorescent protein (GFP)-tagged PK-Calpha, 3B was able to translocate the full length protein to the membrane with an optimal dose of 100 microM in CHO-K1 cells, while diC8 failed to achieve translocation even at doses 3-fold higher. Molecular modeling of 3B into an empty C1b domain of PK-Cdelta clearly showed the existence of a preferred binding orientation. In addition, molecular dynamic simulations suggest that binding discrimination could result from a favorable van der Waals (VDW) interaction between the large, branched sn-1 acyl group of 3B and the aromatic rings of Trp252 (PK-Cdelta) or Tyr252 (PK-Calpha). The DAG analogue of 3B in which the acyl groups are reversed (2C) showed a decrease in binding affinity reflecting the capacity of PK-C to effectively discriminate between alternative orientations of the acyl chains.

Antiretroviral agents as inhibitors of both human immunodeficiency virus type 1 integrase and protease.

The human immunodeficiency virus type one integrase (HIV-1 integrase) is required for integration of a double-stranded DNA copy of the viral RNA genome into a host chromosome and for HIV replication. We have previously reported that phenolic moieties in compounds such as flavones, caffeic acid phenethyl ester (CAPE), tyrphostins, and curcumin confer inhibitory activity against HIV-1 integrase. We have investigated the actions of several recently described protease inhibitors, possessing novel structural features, on HIV-1 integrase. NSC 158393, which contains four 4-hydroxycoumarin residues, was found to exhibit antiviral, antiprotease, and antiintegrase activity. Both the DNA binding and catalytic activities (3'-processing and strand transfer) of integrase were inhibited at micromolar concentrations. Disintegration catalyzed by an integrase mutant containing only the central catalytic domain was also inhibited, indicating that the binding site for these compounds resides in the central 50-212 amino acids of HIV-1 integrase. Binding at or near the integrase catalytic site was also suggested by a global inhibition of the choice of attacking nucleophile in the 3'-processing reaction. NSC 158393 inhibited HIV-2, feline, and simian immunodeficiency virus integrases while eukaryotic topoisomerase I was inhibited at higher concentrations, suggesting selective inhibition of retroviral integrases. Molecular modeling studies revealed that the two hydroxyls and two carbonyl moieties in NSC 158393 may represent essential elements of the pharmacophore. Antiviral efficacy was observed with NSC 158393 derivatives that inhibited both HIV protease and integrase, and the most potent integrase inhibitors also inhibited HIV protease. Hydroxycoumarins may provide lead compounds for development of novel antiviral agents based upon the concurrent inhibition of two viral targets, HIV-1 integrase and protease.

Discovery of novel, non-peptide HIV-1 protease inhibitors by pharmacophore searching.

Fifteen novel non-peptide HIV-1 protease inhibitors were identified by flexible 3D database pharmacophore searching of the NCI DIS 3D database. The pharmacophore query used in the search was derived directly from the X-ray determined structures of protease/inhibitor complexes. These 15 inhibitors, belonging to nine different chemical classes, are promising leads for further development. The two best inhibitors found, NSC 32180, a "dimer" of 4-hydroxycoumarin, and NSC 117027, a "tetramer" of 2-hydroxy quinone, had ID50 values of 0.32 and 0.75 microM for HIV-1 protease inhibition, respectively, and two other inhibitors had ID50 values close to 1 microM. Among the potent inhibitors, NSC 158393 not only demonstrated activity against HIV-1 protease (ID50 1.7 microM) but also exhibited promising antiviral activity in HIV-1-infected CEM-SS cells (EC50 = 11.5 microM). Validation of the pharmacophore used in the search was accomplished by conformational analysis. The binding modes of the most potent inhibitor found in our studies, NSC 32180, were predicted employing docking and molecular dynamics techniques.

Conformationally constrained analogues of diacylglycerol. 18. The incorporation of a hydroxamate moiety into diacylglycerol-lactones reduces lipophilicity and helps discriminate between sn-1 and sn-2 binding modes to protein kinase C (PK-C). Implications for isozyme specificity.

An approach to reduce the log P in a series of diacylglycerol (DAG)-lactones known for their high binding affinity for protein kinase C (PK-C) is presented. Branched alkyl groups with reduced lipophilicity were selected and combined with the replacement of the ester or lactone oxygens by NH or NOH groups. Compound 6a with an isosteric N-hydroxyl amide arm represents the most potent and least lipophilic DAG analogue known to date.

Inhibition of (cytosine C5)-methyltransferase by oligonucleotides containing flexible (cyclopentane) and conformationally constrained (bicyclo[3.1.0]hexane) abasic sites.

Pseudorotationally locked sugar analogues based on bicyclo[3.1.0]-hexane templates were placed in DNA duplexes as abasic target sites in the M. HhaI recognition sequence. The binding affinity of the enzyme increases when the abasic site is constrained to the South conformation and decreases when it is constrained to the North conformation. A structural understanding of these differences is provided.

Pharmacophores in drug design and discovery.

Compounds containing a specific pharmacophore--the minimum structural features necessary for enzyme binding--can be retrieved from a database such as the National Cancer Institute repository by means of three-dimensional (3D) searching, which allows the retrieval of all compounds possessing a specified set of atoms with a given 3D geometry. The ways in which pharmacophores can be found and characterized and the details of the 3D searching methods are described. Elaboration of compounds found in such searches and their subsequent development as lead drugs is also discussed.

Computer-aided design and discovery of protein-protein interaction inhibitors as agents for anti-HIV therapy.

Protein-protein interactions (PPI) are involved in most of the essential processes that occur in organisms. In recent years, PPI have become the object of increasing attention in drug discovery, particularly for anti-HIV drugs. Although the use of combinations of existing drugs, termed highly active antiretroviral therapy (HAART), has revolutionized the treatment of HIV/AIDS, problems with these agents, such as the rapid emergence of drug-resistant HIV-1 mutants and serious adverse effects, have highlighted the need for further discovery of new drugs and new targets. Numerous investigations have shown that PPI play a key role in the virus's life cycle and that blocking or modulating them has a significant therapeutic potential. Here we summarize the recent progress in computer-aided design of PPI inhibitors, mainly focusing on the selection of the drug targets (HIV enzymes and virus entry machinery) and the utilization of peptides and small molecules to prevent a variety of protein-protein interactions (viral-viral or viral-host) that play a vital role in the progression of HIV infection.

QSAR of conformationally flexible molecules: comparative molecular field analysis of protein-tyrosine kinase inhibitors.

Comparative Molecular Field Analysis (CoMFA) has been applied to a study of quantitative structure-activity relationships (QSAR) of conformationally flexible molecules. The relationship between three-dimensional structure and activity of 20 styrene derivatives which inhibit protein-tyrosine kinase was determined. A technique was developed that allows accurate prediction of the inhibitory activity of these molecules and identification in each case of the active conformation. The problem of multiple energetically acceptable conformations was approached in an iterative procedure. Use was made of the varying degrees of symmetry among the molecules. First, CoMFA QSAR models were developed using only those compounds that possess a symmetrical substituent pattern on the phenyl ring. These CoMFA models were then used to select the active conformers of the less symmetrical compounds in the set. Allowing multiple conformers for each compound in the dataset yielded higher crossvalidated r2 values and better predictivity of the QSAR models. Different probe atoms (C+, O-, neutral C) were explored, the O- probe atom exhibiting the highest selectivity in the conformer selection process.

Chem-X and CAMBRIDGE. Comparison of computer generated chemical structures with X-ray crystallographic data.

The structures of a number of molecules as determined by X-ray crystallography have been compared with the structures for the same molecules as calculated by the 3D structure generation program, Chem-X. In the group of molecules examined, ChemModel produced structures that were essentially identical to those based upon X-ray data in 57% of the cases. The corresponding figure for the widely used alternative model builder, CONCORD, was 38%. The superior performance of ChemModel was due entirely to that program's ability to generate multiple structures covering the entire conformational space.

Molecular modeling in the discovery of drug leads.

The National Cancer Institute of the U.S.A. maintains a repository of about 500,000 chemicals which it has tested at some time for anticancer activity. As new chemotherapeutic targets present themselves, methods have been developed by which this large database can be re-examined without resorting to expensive high-volume biological screening. Electronic screening, the method described in this paper, begins with the identification of a target enzyme. The pharmacophore used by the enzyme in binding to substrates is identified, and then all compounds in the database that contain the pharmacophore are found. These compounds are then further filtered, for example, by physical properties such as solubility, and the relatively small number of compounds that survive are submitted for biological testing. This use of a primary electronic screen in the search for ligands of protein kinase C is described. The screen is very fast, and the method is quite generally applicable to different enzymes.

Conformational changes of small molecules binding to proteins.

Flexible molecules change their conformation upon binding to a protein. This was shown by the analysis of small molecules whose structures have been determined by X-ray crystallography of both the pure compound and the compound bound to a protein. Thirty-three compounds present both in the Cambridge Structural Database and the Brookhaven Protein Data Bank were analyzed, and both were compared with the global energy minimum conformation calculated by the molecular mechanics program CHARMm. It was found that the conformation bound to the protein differs from that in the crystal structure and also from that of the global energy minimum, and the degree of deformation depends upon the number of freely rotatable bonds in the molecule. Analysis of the conformational energies of the flexible molecules showed that, for most of those compounds, both the crystal and the protein-bound conformations are energetically well above the global minimum, and, in many cases, not even in any local energy minimum. Semi-empirical calculations performed for a select number of structures, using both the AM1 and PM3 hamiltonians, confirmed these results. These findings are discussed as to their impact upon contemporary methods of drug design.

Effects of tyrphostins, protein kinase inhibitors, on human immunodeficiency virus type 1 integrase.

Efficient replication of HIV-1 requires establishment of the proviral state, i.e., the integration of a DNA copy of the viral genome, synthesized by reverse transcriptase, into a chromosome of the host cell. Integration is catalyzed by the viral integrase protein. We have previously reported that phenolic moieties in compounds such as napthoquinones, flavones, caffeic acid phenethyl ester (CAPE), and curcumin confer inhibitory activity against HIV-1 integrase. We have extended these findings by examining the effects of tryphostins, tyrosine kinase inhibitors. The catalytic activities of HIV-1 integrase and the formation of enzyme-DNA complexes using photocross-linking were examined. Both steps of the integration reaction, 3'-processing and strand transfer, were inhibited by tyrphostins at micromolar concentrations. The DNA binding activity of integrase was inhibited at higher concentrations of tryphostins. Disintegration, an apparent reversal of the strand transfer reaction, catalyzed by an integrase mutant lacking the N-terminal zinc finger and C-terminal DNA binding domains is also inhibited by tyrphostins, indicating that the binding site for these compounds resides in the central catalytic core of HIV-1 integrase. Binding of tyrphostins at or near the integrase catalytic site was also suggested by experiments showing a global inhibition of the choice of attacking nucleophile in the 3'-processing reaction. None of the tryphostins tested inhibited eukaryotic topoisomerase I, even at 100 microM, suggesting selectivity for integrase inhibition. Molecular-modeling studies have revealed that, after energy minimization, several tyrphostins may adopt folded conformations. The similarity of the tyrphostin family to other families of inhibitors is discussed. Tyrphostins may provide lead compounds for development of novel antiviral agents for the treatment of acquired immunodeficiency syndrome based upon inhibition of HIV-1 integrase.