Structure-based rational quest for potential novel inhibitors of human HMG-CoA reductase by combining CoMFA 3D QSAR modeling and virtual screening

3-Hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) catalyzes the formation of mevalonate. In many classes of organisms, this is the committed step leading to the synthesis of essential compounds, such as cholesterol. However, a high level of cholesterol is an important risk factor for coronary heart disease, for which an effective clinical treatment is to block HMGR using inhibitors like statins. Recently the structures of catalytic portion of human HMGR complexed with six different statins have been determined by a delicate crystallography study (Istvan and Deisenhofer Science 2001, 292, 1160-1164), which established a solid basis of structure and mechanism for the rational design, optimization, and development of even better HMGR inhibitors. In this study, three-dimensional quantitative structure-activity relationship (3D QSAR) with comparative molecular field analysis (CoMFA) was performed on a training set of up to 35 statins and statin-like compounds. Predictive models were established by using two different ways: (1) Models-fit, obtained by SYBYL conventional fit-atom molecular alignment rule, has cross-validated coefficients (q2) up to 0.652 and regression coefficients (r2) up to 0.977. (2) Models-dock, obtained by FlexE by docking compounds into the HMGR active site, has cross-validated coefficients (q2) up to 0.731 and regression coefficients (r2) up to 0.947. These models were further validated by an external testing set of 12 statins and statin-like compounds. Integrated with CoMFA 3D QSAR predictive models, molecular surface property (electrostatic and steric) mapping and structure-based (both ligand and receptor) virtual screening have been employed to explore potential novel hits for the HMGR inhibitors. A representative set of eight new compounds of non-statin-like structures but with high pIC(50) values were sorted out in the present study.     

Use of robust classification techniques for the prediction of human cytochrome P450 2D6 inhibition

A new in silico model is developed to predict cytochrome P450 2D6 inhibition from 2D chemical structure. Using a diverse training set of 100 compounds with published inhibition constants, an ensemble approach to recursive partitioning is applied to create a large number of classification trees, each of which yields a yes/no prediction about inhibition for a given compound. These binary classifications are combined to provide an overall prediction, which answers the yes/no question about inhibition and provides a measure of confidence about that prediction. Compared to single-tree models, the ensemble approach is less sensitive to noise in the experimental data as well as to changes in the training set. Internal validation tests indicated an overall classification accuracy of 75%, whereas predictions applied to an external set of 51 compounds yielded 80% accuracy, with all inhibitors correctly identified. The speed and 2D nature of this model make it appropriate for high-throughput processing of large chemical libraries, and the confidence level provides a continuous scale on which to prioritize compounds.     

Prediction of cytochrome P450 3A4, 2D6, and 2C9 inhibitors and substrates by using support vector machines

Statistical learning methods have been used in developing filters for predicting inhibitors of two P450 isoenzymes, CYP3A4 and CYP2D6. This work explores the use of different statistical learning methods for predicting inhibitors of these enzymes and an additional P450 enzyme, CYP2C9, and the substrates of the three P450 isoenzymes. Two consensus support vector machine (CSVM) methods, "positive majority" (PM-CSVM) and "positive probability" (PP-CSVM), were used in this work. These methods were first tested for the prediction of inhibitors of CYP3A4 and CYP2D6 by using a significantly higher number of inhibitors and noninhibitors than that used in earlier studies. They were then applied to the prediction of inhibitors of CYP2C9 and substrates of the three enzymes. Both methods predict inhibitors of CYP3A4 and CYP2D6 at a similar level of accuracy as those of earlier studies. For classification of inhibitors of CYP2C9, the best CSVM method gives an accuracy of 88.9% for inhibitors and 96.3% for noninhibitors. The accuracies for classification of substrates and nonsubstrates of CYP3A4, CYP2D6, and CYP2C9 are 98.2 and 90.9%, 96.6 and 94.4%, and 85.7 and 98.8%, respectively. Both CSVM methods are potentially useful as filters for predicting inhibitors and substrates of P450 isoenzymes. These methods generally give better accuracies than single SVM classification systems, and the performance of the PP-CSVM method is slightly better than that of the PM-CSVM method.     

Discovery of potential M2 channel inhibitors based on the amantadine scaffold via virtual screening and pharmacophore modeling

The M2 channel protein on the influenza A virus membrane has become the main target of the anti-flu drugs amantadine and rimantadine. The structure of the M2 channel proteins of the H3N2 (PDB code 2RLF) and 2009-H1N1 (Genbank accession number GQ385383) viruses may help researchers to solve the drug-resistant problem of these two adamantane-based drugs and develop more powerful new drugs against influenza A virus. In the present study, we searched for new M2 channel inhibitors through a combination of different computational methodologies, including virtual screening with docking and pharmacophore modeling. Virtual screening was performed to calculate the free energies of binding between receptor M2 channel proteins and 200 new designed ligands. After that, pharmacophore analysis was used to identify the important M2 protein-inhibitor interactions and common features of top binding compounds with M2 channel proteins. Finally, the two most potential compounds were determined as novel leads to inhibit M2 channel proteins in both H3N2 and 2009-H1N1 influenza A virus.     

Drug disposition of chiral and achiral drug substrates metabolized by cytochrome P450 2D6 isozyme: case studies, analytical perspectives and developmental implications

The concepts of drug development have evolved over the last few decades. Although number of novel chemical entitities belonging to varied classes have made it to the market, the process of drug development is challenging, intertwined as it is with complexities and uncertainities. The intention of this article is to provide a comprehensive review of novel chemical entities (NCEs) that are substrates to cytochrome P450 (CYP) 2D6 isozyme. Topics covered in this review aim: (1) to provide a framework of the importance of CYP2D6 isozyme in the biotransformation of NCEs as stand-alones and/or in conjunction with other CYP isozymes; (2) to provide several case studies of drug disposition of important drug substrates, (3) to cover key analytical perspectives and key assay considerations to assess the role and involvement of CYP2D6, and (4) to elaborate some important considerations from the development point of view. Additionally, wherever applicable, special emphasis is provided on chiral drug substrates in the various subsections of the review.     

Identification of potential inhibitors of human methionine aminopeptidase (type II) for cancer therapy: Structure-based virtual screening,ADMET prediction and molecular dynamics studies

Methionine Aminopeptidases MetAPs are divalent-cofactor dependent enzymes that are responsible for the cleavage of the initiator Methionine from the nascent polypeptides. MetAPs are classified into two isoforms: namely, MetAP1 and MetAP2. Several studies have revealed that MetAP2 is upregulated in various cancers, and its inhibition has shown to suppress abnormal or excessive blood vessel formation and tumor growth in model organisms. Clinical studies show that the natural product fumagillin, and its analogs are potential inhibitors of MetAP2. However, due to their poor pharmacokinetic properties and neurotoxicities in clinical studies, their further developments have received a great setback. Here, we apply structure-based virtual screening and molecular dynamics methods to identify a new class of potential inhibitors for MetAP2. We screened Otava’s Chemical Library, which consists of about 3 200 000 tangible-chemical compounds, and meticulously selected the top 10 of these compounds based on their inhibitory potentials against MetAP2. The top hit compounds subjected to ADMET predictor using 3 independent ADMET prediction programs, were found to be drug-like. To examine the stability of ligand binding mode, and efficacy, the unbound form of MetAP2, its complexes with fumagillin, spiroepoxytriazole, and the best promising compounds compound-3369841 and compound-3368818 were submitted to 100 ns molecular dynamics simulation. Like fumagillin, spiroepoxytriazole, and both compound-3369841 and compound-3368818 showed stable binding mode over time during the simulations. Taken together, these uninherited-fumagillin compounds may serve as new class of inhibitors or provide scaffolds for further optimization towards the design of more potent MetAP2 inhibitors -development of such inhibitors would be essential strategy against various cancer types.     

Molecular modeling studies of human immunodeficiency virus type 1 protease inhibitors using three‐dimensional quantitative structure‐activity relationship,virtual screening,and docking simulations

In this paper, two 3‐dimensional quantitative structure‐activity relationship models for 60 human immunodeficiency virus (HIV)‐1 protease inhibitors were established using random sampling analysis on molecular surface and translocation comparative molecular field vector analysis (Topomer CoMFA). The non–cross‐validation (r²), cross‐validation (q²), correlation coefficient of external validation (Q²_ext), and F of 2 models were 0.94, 0.80, 0.79, and 198.84 and 0.94, 0.72, 0.75, and 208.53, respectively. The results indicated that 2 models were reasonable and had good prediction ability. Topomer Search was used to search R groups in the ZINC database, 20 new compounds were designed, and the Topomer CoMFA model was used to predicate the biological activity. The results showed that 18 new compounds were more active than the template molecule. So the Topomer Search is effective in screening and can guide the design of new HIV/AIDS drugs. The mechanism of action was studied by molecular docking, and it showed that the protease inhibitors and Ile50, Asp25, and Arg8 sites of HIV‐1 protease have interactions. These results have provided an insight for the design of new potent inhibitors of HIV‐1 protease.     

High-throughput cytochrome P450 (CYP) inhibition screening via a cassette probe-dosing strategy. VI. Simultaneous evaluation of inhibition potential of drugs on human hepatic isozymes CYP2A6, 3A4, 2C9, 2D6 and 2E1

The inhibition potential of drugs towards five major human hepatic cytochrome P450 (CYP) isozymes (CYP2A6, 3A4, 2C9, 2D6, and 2E1) was investigated via cassette dosing of the five probe substrates (coumarin, midazolam, tolbutamide, dextromethorphan, and chlorzoxazone) in human liver microsomes using a 96-well plate format. After microsomal incubations had been terminated with formic acid, the five marker metabolites (7-hydroxycoumarin, 1'-hydroxymidazolam, 4-hydroxytolbutamide, dextrorphan, and 6-hydroxychlorzoxazone) were simultaneously quantified using direct injection/online guard cartridge extraction/tandem mass spectrometry (DI-GCE/MS/MS). Several advantages resulted from the use of a short C(18) guard cartridge (4 mm in length) for DI-GCE/MS/MS, including minimal sample preparation, fast online extraction, short analysis time (2.5 min), and minimal source contamination. In addition, this method demonstrated an inter-day accuracy range from -8.7 - 7.4% with a precision less than 8.3% for the quantification of all the marker metabolites. The inhibition assay for the five CYP isozymes was evaluated using their known selective inhibitors via individual and cassette dosing of the probe substrates. The IC(50) values measured via cassette dosing were consistent with those observed via individual dosing, which were all in agreement with the reported values. In addition, the validated assay was used to evaluate the inhibitory potential of 23 generic drugs (randomly selected) towards the five CYP isozymes. The results suggest the integration of the cassette dosing strategy and the DI-GCE/MS/MS method can provide a reliable in vitro approach to screening the inhibitory potential of new chemical entities, with maximal throughput and cost-effectiveness, in support of drug discovery and development.     

Pharmacophore identification, in silico screening, and virtual library design for inhibitors of the human factor Xa

Factor Xa inhibitors are innovative anticoagulant agents that provide a better safety/efficacy profile compared to other anticoagulative drugs. A chemical feature-based modeling approach was applied to identify crucial pharmacophore patterns from 3D crystal structures of inhibitors bound to human factor Xa (Pdb entries 1fjs, 1kns, 1eqz) using the software LIGANDSCOUT and CATALYST. The complex structures were selected regarding the criteria of high inhibitory potency (i.e. all ligands show K(i) values against factor Xa in the subnanomolar range) and good resolution (i.e. at least 2.2 A) in order to generate selective and high quality pharmacophore models. The resulting chemical-feature based hypotheses were used for virtual screening of commercial molecular databases such as the WDI database. Furthermore, a ligand-based molecular modeling approach was performed to obtain common-feature hypotheses that represent the relevant chemical interactions between 10 bioactive factor Xa inhibitors and the protein, respectively. In a next step a virtual combinatorial library was designed in order to generate new compounds with similar chemical and spatial properties as known inhibitors. The software tool ILIB DIVERSE was used for this procedure in order to provide new scaffolds of this group of anticoagulants. Finally we present the combination of these two techniques, hence virtual screening was performed with selective pharmacophore models in a focused virtual combinatorial database. De novo derived molecular scaffolds that were able to adequately satisfy the pharmacophore criteria are revealed and are promising templates for candidates for further development.     

Novel virtual screening protocol based on the combined use of molecular modeling and electron-ion interaction potential techniques to design HIV-1 integrase inhibitors

HIV-1 integrase (IN) is an essential enzyme for viral replication and represents an intriguing target for the development of new drugs. Although a large number of compounds have been reported to inhibit IN in biochemical assays, no drug active against this enzyme has been approved by the FDA so far. In this study, we report, for the first time, the use of the electron-ion interaction potential (EIIP) technique in combination with molecular modeling approaches for the identification of new IN inhibitors. An innovative virtual screening approach, based on the determination of both short- and long-range interactions between interacting molecules, was employed with the aim of identifying molecules able to inhibit the binding of IN to viral DNA. Moreover, results from a database screening on the commercial Asinex Gold Collection led to the selection of several compounds. One of them showed a significant inhibitory potency toward IN in the overall integration assay. Biological investigations also showed, in agreement with modeling studies, that these compounds prevent recognition of DNA by IN in a fluorescence fluctuation assay, probably by interacting with the DNA binding domain of IN.     

Molecular dynamics, density functional, ADMET predictions, virtual screening, and molecular interaction field studies for identification and evaluation of novel potential CDK2 inhibitors in cancer therapy

In this work, we have used molecular dynamics, density functional theory, virtual screening, ADMET predictions, and molecular interaction field studies to design and propose eight novel potential inhibitors of CDK2. The eight molecules proposed showed interesting structural characteristics that are required for inhibiting the CDK2 activity and show potential as drug candidates for the treatment of cancer. The parameters related to the Rule of Five were calculated, and only one of the molecules violated more than one parameter. One of the proposals and one of the drug-like compounds selected by virtual screening indicated to be promising candidates for CDK2-based cancer therapy.     

Discovery of potential inhibitors targeting the kinase domain of polynucleotide kinase/phosphatase (PNKP): Homology modeling,virtual screening based on multiple conformations,and molecular dynamics simulation

In recent years, the level of interest has been increased in developing the DNA-repair inhibitors, to enhance the cytotoxic effects in the treatment of cancers. Polynucleotide kinase/phosphatase (PNKP) is a critical human DNA repair enzyme that repairs DNA strand breaks by catalyzing the restoration of 5’-phosphate and 3’-hydroxyl termini that are required for subsequent processing by DNA ligases and polymerases. PNKP is the only protein that repairs the 3′-hydroxyl group and 5′-phosphate group, which depicts PNKP as a potential therapeutic target. Besides, PNKP is the only DNA-repair enzyme that contains the 5′-kinase activity, therefore, targeting this kinase domain would motivate the development of novel PNKP-specific inhibitors. However, there are neither crystal structures of human PNKP nor the kinase inhibitors reported so far. Thus, in this present study, a sequential molecular docking-based virtual screening with multiple PNKP conformations integrating homology modeling, molecular dynamics simulation, and binding free energy calculation was developed to discover novel PNKP kinase inhibitors, and the top-scored molecule was finally submitted to molecular dynamics simulation to reveal the binding mechanism between the inhibitor and PNKP. Taken together, the current study could provide some guidance for the molecular docking based-virtual screening of novel PNKP kinase inhibitors.     

Validation of higher-throughput high-performance liquid chromatography/atmospheric pressure chemical ionization tandem mass spectrometry assays to conduct cytochrome P450s CYP2D6 and CYP3A4 enzyme inhibition studies in human liver microsomes

In the early stage of drug discovery, thousands of new chemical entities (NCEs) may be screened before a single drug candidate can be identified for development. In order to accelerate the drug discovery process, we have developed higher-throughput enzyme assays to evaluate the inhibition of cytochrome P450 isoforms 2D6 (CYP2D6) and 3A4 (CYP3A4) in human liver microsomes. The assays are based on high-performance liquid chromatography/tandem mass spectrometry (LC/MS/MS) techniques. The analysis time for each sample was reduced from approximately 20 minutes for the conventional HPLC assay to 30 seconds for the LC/MS/MS assay. For both LC/MS/MS assays, the linearity (r(2) > 0.99), precision (%CV < 15%) and accuracy (% bias <15%) for both inter- and intraday validations were satisfactory. Since the implementation of the LC/MS/MS assays, our sample throughput has increased by over 40-fold.     

Development and application of a method to investigate drug-metabolizing enzyme inhibitors using sparteine for probe of cytochrome P450 2D6 and tris(2,2'-bipyridine)ruthenium(II)-electrogenerated chemiluminescence detection

We studied the detection of drug-metabolizing enzyme inhibitiors using column-switching high performance liquid chromatography with tris(2,2'-bipyridine)ruthenium(II) (Ru(bpy)(3)(2+))-electrogenerated chemiluminescence detection. This can be applied to evaluate the genetic diversity concerning the ability of cytochrome P450 (CYP) 2D6 to metabolize drug in vitro. We demonstrated the ability of CYP2D6 to enable us to examine drugs metabolizing enzyme inhibition with high performance and sensitivity. This method can be applied to investigate metabolite inhibitors of CYP2D6 in vitro and in vivo. Thus, Metixene was found to be a potential CYP2D6 inhibitor.     

QM/MM study of the active species of the human cytochrome P450 3A4, and the influence thereof of the multiple substrate binding

Cytochrome P450 3A4 is involved in the metabolism of 50% of all swallowed drugs. The enzyme functions by means of a high-valent iron-oxo species, called compound I (Cpd I), which is formed after entrance of the substrate to the active site. We explored the features of Cpd I using hybrid quantum mechanical/molecular mechanical calculations on various models that are either substrate-free or containing one and two molecules of diazepam as a substrate. M?ssbauer parameters of Cpd I were computed. Our major finding shows that without the substrate, Cpd I tends to elongate its Fe-S bond, localize the radical on the sulfur, and form hydrogen bonds with A305 and T309, which may hypothetically lead to Cpd I consumption by H-abstraction. However, the positioning of diazepam close to Cpd I, as enforced by the effector molecule, was found to strengthen the NH...S interactions of the conserved I443 and G444 residues with the proximal cysteinate ligand. These interactions are known to stabilize the Fe-S bond, and as such, the presence of the substrate leads to a shorter Fe-S bond and it prevents the localization of the radical on the sulfur. This diazepam-Cpd I stabilization was manifested in the 1W0E conformer. The effector substrate did not influence Cpd I directly but rather by positioning the active substrate close to Cpd I, thus displacing the hydrogen bonds with A305 and T309, and thereby giving preference to substrate oxidation. It is hypothesized that these effects on Cpd I, promoted by the restrained substrate, may be behind the special metabolic behavior observed in cases of multiple substrate binding (also called cooperative binding). This restraint constitutes a mechanism whereby substrates stabilize Cpd I sufficiently long to affect monooxygenation by P450s at the expense of Cpd I destruction by the protein residues.     

High-throughput cytochrome P450 inhibition screening via cassette probe-dosing strategy. IV. Validation of a direct injection on-line guard cartridge extraction/tandem mass spectrometry method for simultaneous CYP3A4, 2D6 and 2E1 inhibition assessment

A highly efficient direct injection on-line guard cartridge extraction/tandem mass spectrometry (DI-GCE/MS/MS) method has been validated for high-throughput evaluation of cytochrome P450 (CYP) 3A4, 2D6 and 2E1 inhibition potential via cassette dosing of midazolam, dextromethorphan and chlorzoxazone using human hepatic microsomes and 96-well microtiter plates. Microsomal incubations were terminated with formic acid, centrifuged, and the resulting supernatants were injected for analysis by DI-GCE/MS/MS. Due to the novel use of an extremely short C(18) guard cartridge (4 mm in length), this method exhibits several advantages such as no sample preparation, excellent on-line extraction, short run time (2.5 min), and minimized source contamination and performance deterioration. The DI-GCE/MS/MS method demonstrates acceptable accuracy and precision for the simultaneous quantification of 1'-hydroxymidazolam, dextrorphan and 6-hydroxychlorzoxazone in microsomal incubations. The inhibition potential of CYP3A4, 2D6 and 2E1 has been evaluated using their known selective inhibitors. The IC(50) values measured by the cassette dosing approach (high-throughput) are consistent with those observed by an individual dosing regimen (conventional) and are all in good agreement with the literature values. The results suggest that the cassette probe-dosing strategy may provide an in vitro approach to minimize cost while maximizing throughput of CYP inhibition evaluation of new chemical entities in support of drug discovery and development.