Cytochrome P450 3A4 inhibition by ketoconazole: tackling the problem of ligand cooperativity using molecular dynamics simulations and free-energy calculations

Unraveling the Ligand-Binding Sites of CYP3A4 by Molecular Dynamics Simulations with Solvent Probes

Cytochrome P450 3A4 (CYP3A4) is one of the most important drug-metabolizing enzymes in the human body and is well known for its complicated, atypical kinetic characteristics. The existence of multiple ligand-binding sites in CYP3A4 has been widely recognized as being capable of interfering with the active pocket through allosteric effects. The identification of ligand-binding sites other than the canonical active site above the heme is especially important for understanding the atypical kinetic characteristics of CYP3A4 and the intriguing association between the ligand and the receptor. In this study, we first employed mixed-solvent molecular dynamics (MixMD) simulations coupled with the online computational predictive tools to explore potential ligand-binding sites in CYP3A4. The MixMD approach demonstrates better performance in dealing with the receptor flexibility compared with other computational tools. From the sites identified by MixMD, we then picked out multiple sites for further exploration using ensemble docking and conventional molecular dynamics (cMD) simulations. Our results indicate that three extra sites are suitable for ligand binding in CYP3A4, including one experimentally confirmed site and two novel sites.

Spatial engineering of single-atom Fe adjacent to Cu-assisted nanozymes for biomimetic O2 activation

The precise design of single-atom nanozymes (SAzymes) and understanding of their biocatalytic mechanisms hold great promise for developing ideal bio-enzyme substitutes. While considerable efforts have been directed towards mimicking partial bio-inspired structures, the integration of heterogeneous SAzymes configurations and homogeneous enzyme-like mechanism remains an enormous challenge. Here, we show a spatial engineering strategy to fabricate dual-sites SAzymes with atomic Fe active center and adjacent Cu sites. Compared to planar Fe-Cu dual-atomic sites, vertically stacked Fe-Cu geometry in FePc@2D-Cu-N-C possesses highly optimized scaffolds, favorable substrate affinity, and fast electron transfer. These characteristics of FePc@2D-Cu-N-C SAzyme induces biomimetic O2 activation through homogenous enzymatic pathway, resembling functional and mechanistic similarity to natural cytochrome c oxidase. Furthermore, it presents an appealing alternative of cytochrome P450 3A4 for drug metabolism and drug-drug interaction. These findings are expected to deepen the fundamental understanding of atomic-level design in next-generation bio-inspired nanozymes.

Polyfluoroalkyl Chain-Based Assemblies for Biomimetic Catalysis

Amphiphobic fluoroalkyl chains are exploited for creating robust and diverse self-assembled biomimetic catalysts. Long terminal perfluoroalkyl chains (Cn F2n+1 with n=6, 8, and 10) linked with a short perhydroalkyl chains (Cm H2m with m=2 and 3) were used to synthesize several 1,4,7-triazacyclononane (TACN) derivatives, Cn F2n+1 -Cm H2m -TACN. In the presence of an equimolar amount of Zn2+ ions that coordinate the TACN moiety and drive the self-assembly into micelle-like aggregates, the critical aggregation concentration of polyfluorinated Cn F2n+1 -Cm H2m -TACN⋅Zn2+ was lowered by ∼1 order of magnitude compared to the traditional perhyroalkyl counterpart with identical carbon number of alkyl chain. When 2'-hydroxypropyl-4-nitrophenyl phosphate was used as the model phosphate substrate, polyfluorinated Cn F2n+1 -Cm H2m -TACN⋅Zn2+ assemblies showed higher affinity and catalytic activity, compared to its perhyroalkyl chain-based counterpart. Coarse-grained molecular dynamic simulations have been introduced to explore the supramolecular assembly of polyfluoroalkyl chains in the presence of Zn2+ ions and to better understand their enhanced catalytic activity.

Evaluation of novel pyrazol-4-yl pyridine derivatives possessing arylsulfonamide tethers as c-Jun N-terminal kinase (JNK) inhibitors in leukemia cells

A series of 36 pyrazol-4-yl pyridine derivatives (8a-i, 9a-i, 10a-i, and 11a-i) was designed, synthesized, and evaluated for its antiproliferative activity over NCI-60 cancer cell line panel and inhibitory effect against JNK isoforms (JNK1, JNK2, and JNK3). All the synthesized compounds were tested against the NCI-60 cancer cell line panel. Compounds 11b, 11c, 11g, and 11i were selected to determine their GI50s and exerted a superior potency over the reference standard SP600125 against the tested cell lines. 11c showed a GI50 of 1.28 μM against K562 leukemic cells. Vero cells were used to assess 11c cytotoxicity compared to the tested cancer cells. The target compounds were tested against hJNK isoforms in which compound 11e exhibited the highest potency against JNK isoforms with IC50 values of 1.81, 12.7, and 10.5 nM against JNK1, JNK2, and JNK3, respectively. Kinase profiling of 11e showed higher JNK selectivity in 50 kinase panels. Compounds 11c and 11e showed cell population arrest at the G2/M phase, induced early apoptosis, and slightly inhibited beclin-1 production at higher concentrations in K562 leukemia cells relative to SP600125. NanoBRET assay of 11e showed intracellular JNK1 inhibition with an IC50 of 2.81 μM. Also, it inhibited CYP2D6 and 3A4 with different extent and its hERG activity showed little cardiac toxicity with an IC50 of 4.82 μM. hJNK3 was used as a template to generate the hJNK1 crystal structure to explore the binding mode of 11e (PDB ID: 8ENJ) with a resolution of 2.8 °A and showed a typical type I kinase inhibition against hJNK1. Binding energy scores showed that selectivity of 11e towards JNK1 could be attributed to additional hydrophobic interactions relative to JNK3.

Revisiting Supersaturation of a Biopharmaceutical Classification System IIB Drug: Evaluation via a Multi-Cup Dissolution Approach and Molecular Dynamic Simulation

As a subclass of the biopharmaceutical classification system (BCS) class II, basic drugs (BCS IIB) exhibit pH-dependent solubility and tend to generate supersaturation in the gastrointestinal tract, leading to less qualified in vitro-in vivo correlation (IVIVC). This study aims to develop a physiologically based multi-cup dissolution approach to improve the evaluation of the supersaturation for a higher quality of IVIVC and preliminarily explores the molecular mechanism of supersaturation and precipitation of ketoconazole affected by Polyvinylpyrrolidone-vinyl acetate copolymer (PVPVA) and hydroxypropyl methyl-cellulose (HPMC). The concentration of ketoconazole in each cup of the dynamic gastrointestinal model (DGIM) was measured using fiber optical probes. Molecular interactions between ketoconazole and PVPVA or HPMC were simulated by Materials Studio. The results demonstrated that PVPVA and HPMC improved and maintained the supersaturation of ketoconazole. PVPVA exhibited superior precipitation inhibitory effect on ketoconazole molecule aggregation due to slightly stronger van der Waals forces as well as unique electrostatic forces, thereby further enhancing in vitro drug absorption, which correlated well with in vivo drug absorption. Compared with a conventional dissolution apparatus paddle method, the DGIM improved the mean prediction error through the IVIVC from 19.30% to 9.96%, reaching the qualification criteria. In conclusion, the physiologically based multi-cup dissolution approach enables improved evaluation of supersaturation in gastrointestinal transportation of BCS IIB drug ketoconazole, enabling screening screen precipitation inhibitors and achieving qualified IVIVC for drug formulation studies.

Distinct effects of form selective cytochrome P450 inhibitors on cytochrome P450-mediated monooxygenase and hydrogen peroxide generating NADPH oxidase

A characteristic of cytochrome P450 (CYP) enzymes is their ability to generate H₂O₂, either directly or indirectly via superoxide anion, a reaction referred to as "NADPH oxidase" activity. H₂O₂ production by CYPs can lead to the accumulation of cytotoxic reactive oxygen species which can compromise cellular functioning and contribute to tissue injury. Herein we determined if form selective CYP inhibitors could distinguish between the activities of the monooxygenase and NADPH oxidase activities of rat recombinant CYP1A2, CYP2E1, CYP3A1 and CYP3A2 and CYP1A1/2-enriched β-naphthoflavone-induced rat liver microsomes, CYP2E1-enriched isoniazide-induced rat liver microsomes and CYP3A subfamily-enriched dexamethasone-induced rat liver microsomes. In the presence of 7,8-benzoflavone (2.0 μM) for CYP1A2 and 4-methylpyrazole (32 μM) or DMSO (16 mM) for CYP2E1, monooxygenase activity was blocked without affecting NADPH oxidase activity for both the recombinant enzymes and microsomal preparations. Ketoconazole (1.0 μM), a form selective inhibitor for CYP3A subfamily enzymes, completely inhibited monooxygenase activity of rat recombinant CYP3A1/3A2 and CYP3A subfamily in rat liver microsomes; it also partially inhibited NADPH oxidase activity. 7,8-benzoflavone is a type I ligand, which competes with substrate binding, while 4-methylpyrazole and DMSO are type II heme binding ligands. Interactions of heme with these type II ligands was not sufficient to interfere with oxygen activation, which is required for NADPH oxidase activity. Ketoconazole, a type II ligand known to bind multiple sites on CYP3A subfamily enzymes in close proximity to heme, also interfered, at least in part, with oxygen activation. These data indicate that form specific inhibitors can be used to distinguish between monooxygenase reactions and H₂O₂ generating NADPH oxidase of CYP1A2 and CYP2E1. Mechanisms by which ketoconazole inhibits CYP3A NADPH oxidase remain to be determined.

Computer-Aided (In Silico) Modeling of Cytochrome P450-Mediated Food–Drug Interactions (FDI)

Modifications of the activity of Cytochrome 450 (CYP) enzymes by compounds in food might impair medical treatments. These CYP-mediated food–drug interactions (FDI) play a major role in drug clearance in the intestine and liver. Inter-individual variation in both CYP expression and structure is an important determinant of FDI. Traditional targeted approaches have highlighted a limited number of dietary inhibitors and single-nucleotide variations (SNVs), each determining personal CYP activity and inhibition. These approaches are costly in time, money and labor. Here, we review computational tools and databases that are already available and are relevant to predicting CYP-mediated FDIs. Computer-aided approaches such as protein–ligand interaction modeling and the virtual screening of big data narrow down hundreds of thousands of items in databanks to a few putative targets, to which the research resources could be further directed. Structure-based methods are used to explore the structural nature of the interaction between compounds and CYP enzymes. However, while collections of chemical, biochemical and genetic data are available today and call for the implementation of big-data approaches, ligand-based machine-learning approaches for virtual screening are still scarcely used for FDI studies. This review of CYP-mediated FDIs promises to attract scientists and the general public.

Structural characterization of the homotropic cooperative binding of azamulin to human cytochrome P450 3A5

Cytochrome P450 3A4 and 3A5 catalyze the metabolic clearance of a large portion of therapeutic drugs. Azamulin is used as a selective inhibitor for 3A4 and 3A5 to define their roles in metabolism of new chemical entities during drug development. In contrast to 3A4, 3A5 exhibits homotropic cooperativity for the sequential binding of two azamulin molecules at concentrations used for inhibition. To define the underlying sites and mechanisms for cooperativity, an X-ray crystal structure of 3A5 was determined with two azamulin molecules in the active site that are stacked in an antiparallel orientation. One azamulin resides proximal to the heme in a pose similar to the 3A4-azamulin complex. Comparison to the 3A5 apo structure indicates that the distal azamulin in 3A5 ternary complex causes a significant induced fit that excludes water from the hydrophobic surfaces of binding cavity and the distal azamulin, which is augmented by the stacking interaction with the proximal azamulin. Homotropic cooperativity was not observed for the binding of related pleuromutilin antibiotics, tiamulin, retapamulin, and lefamulin, to 3A5, which are larger and unlikely to bind in the distal site in a stacked orientation. Formation of the 3A5 complex with two azamulin molecules may prevent time-dependent inhibition that is seen for 3A4 by restricting alternate product formation and/or access of reactive intermediates to vulnerable protein sites. These results also contribute to a better understanding of sites for cooperative binding and the differential structural plasticity of 3A5 and 3A4 that contribute to differential substrate and inhibitor binding.

Computational Insight into Biotransformation Profiles of Organophosphorus Flame Retardants to Their Diester Metabolites by Cytochrome P450

Biotransformation of organophosphorus flame retardants (OPFRs) mediated by cytochrome P450 enzymes (CYPs) has a potential correlation with their toxicological effects on humans. In this work, we employed five typical OPFRs including tris(1,3-dichloro-2-propyl) phosphate (TDCIPP), tris(1-chloro-2-propyl) phosphate (TCIPP), tri(2-chloroethyl) phosphate (TCEP), triethyl phosphate (TEP), and 2-ethylhexyl diphenyl phosphate (EHDPHP), and performed density functional theory (DFT) calculations to clarify the CYP-catalyzed biotransformation of five OPFRs to their diester metabolites. The DFT results show that the reaction mechanism consists of Cα-hydroxylation and O-dealkylation steps, and the biotransformation activities of five OPFRs may follow the order of TCEP ≈ TEP ≈ EHDPHP > TCIPP > TDCIPP. We further performed molecular dynamics (MD) simulations to unravel the binding interactions of five OPFRs in the CYP3A4 isoform. Binding mode analyses demonstrate that CYP3A4-mediated metabolism of TDCIPP, TCIPP, TCEP, and TEP can produce the diester metabolites, while EHDPHP metabolism may generate para-hydroxyEHDPHP as the primary metabolite. Moreover, the EHDPHP and TDCIPP have higher binding potential to CYP3A4 than TCIPP, TCEP, and TEP. This work reports the biotransformation profiles and binding features of five OPFRs in CYP, which can provide meaningful clues for the further studies of the metabolic fates of OPFRs and toxicological effects associated with the relevant metabolites.

Discovery of first-in-class imidazothiazole-based potent and selective ErbB4 (HER4) kinase inhibitors

This article reports on novel imidazothiazole derivatives as first-in-class potent and selective ErbB4 (HER4) inhibitors. There are no other reported selective inhibitors of this kinase in the literature, that's why they are considered as first-in-class. In addition, none of the reported non-selective ErbB4 inhibitors possesses imidazothiazole nucleus in its structure. Therefore, there is novelty in this work in both kinase selectivity and chemical structure. Compounds Ik and IIa are the most potent ErbB4 kinase inhibitor (IC₅₀ = 15.24 and 17.70 nM, respectively). Compound Ik showed promising antiproliferative activity. It is selective towards cancer cell lines than normal cells. Its ability to penetrate T-47D cell membrane and inhibit ErbB4 kinase inside the cells has been confirmed. Moreover, both compound Ik and IIa have additional merits such as weak potency against hERG ion channels and against CYP 3A4 and 2D6. Molecular docking and dynamic simulation studies were carried out to explain binding interactions.

Crystal Structures of Drug-Metabolizing CYPs

The complex enzyme kinetics displayed by drug-metabolizing cytochrome P450 enzymes (CYPs) (see Chapter 9 ) can, in part, be explained by an examination of their crystallographic protein structures. Fortunately, despite low sequence similarity between different families of drug-metabolizing CYPs, there exists a high degree of structural homology within the superfamily. This similarity in the protein fold allows for a direct comparison of the structural features of CYPs that contribute toward differences in substrate binding, heterotropic and homotropic cooperativity, and genetic variability in drug metabolism. In this chapter, we first provide an overview of the nomenclature and the role of structural features that are common in all CYPs. We then apply these definitions to understand the different substrate specificities and functions in the CYP3A, CYP2C, and CYP2D families of enzymes.

In Silico Exploration of Conformational Dynamics and Novel Inhibitors for Targeting MEF2-Associated Transcriptional Activity

Myocyte enhancer factor 2 (MEF2; MEF2A-MEF2D) transcription factors regulate gene expression in a variety of developmental processes by binding to AT-rich DNA motifs via highly conserved N-terminal extensions known as MADS-box and MEF2 domains. Despite the fact that MEF2 proteins exhibit high similarity at their N-terminal regions and share a common consensus DNA binding motif, their functional preferences may vary significantly in the adjacent regions to the DNA binding core segment. The current study delineates the conformational paradigm, clustered recognition, and comparative DNA binding preferences for MEF2A and MEF2B-specific MADS-box/MEF2 domains at the YTA(A/T)4TAR consensus motif. In both MEF2A and MEF2B proteins, α1-helix plays a crucial role through acquiring more flexibility by attaining loop conformation. In comparison to apo-MEF2, an outward disposition of the distal portion of α1-helix and movement of its proximal part to β1 allows synergistic repositioning of the α1-α2 linker, C-terminal region, and MEF2 domain, resulting in the formation of a hydrophobic groove for DNA binding. In both instances, conformational switching of the helical content is the main contributing factor while preserving the overall β-topology to maintain the inside-out conformation of subdivided α1-helix flip. Multivariate statistical analysis reveals that MEF2B obscures less accessible conformational space for DNA binding as compared to the MEF2A-DNA complex. The presence of similar structural requirements and conserved residues including Arg10, Phe21, and Arg24 in accentuating the MEF2-specific DNA recognition mechanism led us to perform structure-based virtual screening for isolating novel inhibitors that are able to target MEF2-DNA binding regions. The top hits (acetamide, benzamide, carboxamide, and enamide) obtained through preliminary assay were scrutinized to binding potential analysis at the MEF2-DNA binding groove, energy values, absorption, distribution, toxicity, and Lipinski's rule of five assessments. Based on these findings, we propose valuable active drug-like molecules for selective applications against MEF2A and MEF2B. The current study may help in uncovering the atomistic-level mechanistic DNA binding patterns of MEF2 proteins, and data may be valuable in devising effective therapeutic strategies for MEF2-associated disorders.

Molecular Basis for Metabolic Regioselectivity and Mechanism of Cytochrome P450s toward Carcinogenic 4-(Methylnitrosamino)-(3-pyridyl)-1-butanone

As an abundantly present tobacco component, carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) has also been detected in atmospheric particulate matter, suggesting the ineluctable exposure risk of this contaminant. NNK metabolic activation by cytochrome P450 enzymes (CYPs) is a prerequisite to exerting its genotoxicity, but the metabolic regioselectivity and mechanism are still unknown. Here the binding feature and regioselectivity of CYPs 1A1, 1A2, 2A6, 2A13, 2B6, and 3A4 toward NNK are unraveled through molecular docking and molecular dynamics (MD) simulations. Binding mode analyses reveal that 1A2 and 2B6 have definite preferences for NNK α-methyl hydroxylation, while the other four CYPs preferentially catalyze α-methylene hydroxylation. The binding affinities between NNK and CYPs evaluated by the binding free energies follow the order 2A13 > 2B6 > 1A2 > 2A6 > 1A1 > 3A4. Density functional theory (DFT) calculations are further performed to characterize the mechanism of NNK biotransformation. Results show that the α-hydroxyNNK generated from α-hydroxylation may undergo nonenzymatic decomposition to form genotoxic diazohydroxide and aldehyde, and further oxidation by P450 to yield nitrosamide, which mainly contributes to NNK toxification capacity. Meanwhile the pyridine N-oxidation and denitrosation of C_α-radical intermediate play an important role in detoxifying NNK. Overall, the present study provides the molecular basis for CYP-catalyzed regioselectivity and mechanism of NNK biotransformation, which can enable the identification of metabolites for assessing the health risk of individual NNK exposure.

Binding Modes and Metabolism of Caffeine

A correct estimate of ligand binding modes and a ratio of their occupancies is crucial for calculations of binding free energies. The newly developed method BLUES combines molecular dynamics with nonequilibrium candidate Monte Carlo. Nonequilibrium candidate Monte Carlo generates a plethora of possible binding modes and molecular dynamics enables the system to relax. We used BLUES to investigate binding modes of caffeine in the active site of its metabolizing enzyme Cytochrome P450 1A2 with the aim of elucidating metabolite-formation profiles at different concentrations. Because the activation energies of all sites of metabolism do not show a clear preference for one metabolite over the others, the orientations in the active site must play a key role. In simulations with caffeine located in a spacious pocket above the I-helix, it points N3 and N1 to the heme iron, whereas in simulations where caffeine is in close proximity to the heme N7 and C8 are preferably oriented toward the heme iron. We propose a mechanism where at low caffeine concentrations caffeine binds to the upper part of the active site, leading to formation of the main metabolite paraxanthine. On the other hand, at high concentrations two molecules are located in the active site, forcing one molecule into close proximity to the heme and yielding metabolites theophylline and trimethyluretic acid. Our results offer an explanation of previously published experimental results.

Prediction of Drug-Drug Interactions Related to Inhibition or Induction of Drug-Metabolizing Enzymes

Drug-drug interaction (DDI) is the phenomenon of alteration of the pharmacological activity of a drug(s) when another drug(s) is co-administered in cases of so-called polypharmacy. There are three types of DDIs: pharmacokinetic (PK), pharmacodynamic, and pharmaceutical. PK is the most frequent type of DDI, which often appears as a result of the inhibition or induction of drug-metabolising enzymes (DME). In this review, we summarise in silico methods that may be applied for the prediction of the inhibition or induction of DMEs and describe appropriate computational methods for DDI prediction, showing the current situation and perspectives of these approaches in medicinal and pharmaceutical chemistry. We review sources of information on DDI, which can be used in pharmaceutical investigations and medicinal practice and/or for the creation of computational models. The problem of the inaccuracy and redundancy of these data are discussed. We provide information on the state-of-the-art physiologically- based pharmacokinetic modelling (PBPK) approaches and DME-based in silico methods. In the section on ligand-based methods, we describe pharmacophore models, molecular field analysis, quantitative structure-activity relationships (QSAR), and similarity analysis applied to the prediction of DDI related to the inhibition or induction of DME. In conclusion, we discuss the problems of DDI severity assessment, mention factors that influence severity, and highlight the issues, perspectives and practical using of in silico methods.

How Can Interleukin-1 Receptor Antagonist Modulate Distinct Cell Death Pathways?

Multiple mechanisms of cell death exist (apoptosis, necroptosis, pyroptosis) and the subtle balance of several distinct proteins and inhibitors tightly regulates the cell fate toward one or the other pathway. Here, by combining coimmunoprecipitation, enzyme assays, and molecular simulations, we ascribe a new role, within this entangled regulatory network, to the interleukin-1 receptor antagonist (IL-1Ra). Our study enlightens that IL-1Ra, which usually inhibits the inflammatory effects of IL-1α/β by binding to IL-1 receptor, under advanced pathological states prevents apoptosis and/or necroptosis by noncompetitively inhibiting the activity of caspase-8 and -9. Consensus docking, followed by cumulative 10 μs of molecular dynamics simulations unprecedentedly reveal that IL-1Ra binds both caspases at their dimeric interface, preventing, in this manner, the formation of their catalytically/signaling active form. The resulting IL-1Ra/caspase-8(9) adducts are stabilized by hydrophobic and by few key hydrogen bonding interactions, formed by residues fully conserved across distinct caspases (-3, -6, -7, -8, and -9), and closely resemble the binding mode of the caspases inhibitors XIAP (X-linked inhibitor of apoptosis) and c-FLIP (cellular FLICE-like inhibitory protein). Tight regulation of the different forms of cell death has a major impact on distinct human illnesses (i.e., cancer, neurodegeneration, ischemic injury, atherosclerosis, viral/bacterial infections, and immune reaction). Hence, our study, pinpointing IL-1Ra as new actor of the intricate cell death regulatory network and gaining an atomic-scale understanding of its mechanism may open new avenues toward innovative therapeutic strategies to tackle major human diseases.

Molecular dynamics investigation of stereoselective inhibition mechanism of HIF-2α/ARNT heterodimer

Hypoxia-inducible factors (HIFs) are heterodimeric transcription factors related with the onset and progression of solid tumors. Studies demonstrated a class of tetrazole containing chiral inhibitors could stereoselectively disrupt the HIF-2 dimerization and reduce the target gene expression. However, the dynamical features and structural motifs of the HIF-2 heterodimer caused by the binding of enantiomers have not been rationalized at the atomistic level. In this work, molecular dynamics (MD) simulations combined with adaptive steered MD (ASMD) simulations were used to investigate stereoselective interrupting mechanism of HIF-2. Our results decipher that the binding of ligand A (S, R)-24 begets the significant conformation changes of β-sheets and interrupts the HIF-2α/ARNT heterodimerization, which may be attributed to the disruption of the hydrogen bond and salt bridge interactions formed by the 4 foremost residues (Asp240, Arg247, Glu362, and Arg366) and the destruction of hydrophobic interactions on the binding interface. By contrast, the binding of ligand B (R, S)-24 does not disrupt protein dimerization and causes the motion of Fα helix in HIF-2α PAS-B domain to further change the major tunnel for ligand ingress and engress. The present work provides important molecular-level insight into the effect of the binding enantiomers on HIF-2 heterodimerization and bridges the gap between theory and the experimental results, which may conduce to develop highly potent antagonists for intervening the HIF-2-driven tumors.

Warfarin and Flavonoids Do Not Share the Same Binding Region in Binding to the IIA Subdomain of Human Serum Albumin

Human serum albumin (HSA) binds a variety of xenobiotics, including flavonoids and warfarin. The binding of another ligand to the IIA binding site on HSA can cause warfarin displacement and potentially the elevation of its free concentration in blood. Studies dealing with flavonoid-induced warfarin displacement from HSA provided controversial results: estimated risk of displacement ranged from none to serious. To resolve these controversies, in vitro study of simultaneous binding of warfarin and eight different flavonoid aglycons and glycosides to HSA was carried out by fluorescence spectroscopy as well as molecular docking. Results show that warfarin and flavonoids do not share the same binding region in binding to HSA. Interactions were only observed at high warfarin concentrations not attainable under recommended dosing regimes. Docking experiments show that flavonoid aglycons and glycosides do not bind at warfarin high affinity sites, but rather to different regions within the IIA HSA subdomain. Thus, the risk of clinically significant warfarin-flavonoid interaction in binding to HSA should be regarded as negligible.

Exploring Molecular Mechanisms of Paradoxical Activation in the BRAF Kinase Dimers: Atomistic Simulations of Conformational Dynamics and Modeling of Allosteric Communication Networks and Signaling Pathways

The recent studies have revealed that most BRAF inhibitors can paradoxically induce kinase activation by promoting dimerization and enzyme transactivation. Despite rapidly growing number of structural and functional studies about the BRAF dimer complexes, the molecular basis of paradoxical activation phenomenon is poorly understood and remains largely hypothetical. In this work, we have explored the relationships between inhibitor binding, protein dynamics and allosteric signaling in the BRAF dimers using a network-centric approach. Using this theoretical framework, we have combined molecular dynamics simulations with coevolutionary analysis and modeling of the residue interaction networks to determine molecular determinants of paradoxical activation. We have investigated functional effects produced by paradox inducer inhibitors PLX4720, Dabrafenib, Vemurafenib and a paradox breaker inhibitor PLX7904. Functional dynamics and binding free energy analyses of the BRAF dimer complexes have suggested that negative cooperativity effect and dimer-promoting potential of the inhibitors could be important drivers of paradoxical activation. We have introduced a protein structure network model in which coevolutionary residue dependencies and dynamic maps of residue correlations are integrated in the construction and analysis of the residue interaction networks. The results have shown that coevolutionary residues in the BRAF structures could assemble into independent structural modules and form a global interaction network that may promote dimerization. We have also found that BRAF inhibitors could modulate centrality and communication propensities of global mediating centers in the residue interaction networks. By simulating allosteric communication pathways in the BRAF structures, we have determined that paradox inducer and breaker inhibitors may activate specific signaling routes that correlate with the extent of paradoxical activation. While paradox inducer inhibitors may facilitate a rapid and efficient communication via an optimal single pathway, the paradox breaker may induce a broader ensemble of suboptimal and less efficient communication routes. The central finding of our study is that paradox breaker PLX7904 could mimic structural, dynamic and network features of the inactive BRAF-WT monomer that may be required for evading paradoxical activation. The results of this study rationalize the existing structure-functional experiments by offering a network-centric rationale of the paradoxical activation phenomenon. We argue that BRAF inhibitors that amplify dynamic features of the inactive BRAF-WT monomer and intervene with the allosteric interaction networks may serve as effective paradox breakers in cellular environment.

Elucidating a chemical defense mechanism of Antarctic sponges: A computational study

In 2000, a novel secondary metabolite (erebusinone, Ereb) was isolated from the Antarctic sea sponge, Isodictya erinacea. The bioactivity of Ereb was investigated, and it was found to inhibit molting when fed to the arthropod species Orchomene plebs. Xanthurenic acid (XA) is a known endogenous molt regulator present in arthropods. Experimental studies have confirmed that XA inhibits molting by binding to either (or both) of two P450 enzymes (CYP315a1 or CYP314a1) that are responsible for the final two hydroxylations in the production of the molt-inducing hormone, 20-hydroxyecdysone (20E). The lack of crystal structures and biochemical assays for CYP315a1 or CYP314a1, has prevented further experimental exploration of XA and Ereb's molt inhibition mechanisms. Herein, a wide array of computational techniques - homology modeling, molecular dynamics simulations, binding site bioinformatics, flexible receptor-flexible ligand docking, and molecular mechanics-generalized Born surface area calculations - have been employed to elucidate the structure-function relationships between the aforementioned P450s and the two described small molecule inhibitors (Ereb and XA). Results indicate that Ereb likely targets CYP315a1 by interacting with a network of aromatic residues in the binding site, while XA may inhibit both CYP315a1 and CYP314a1 because of its aromatic, as well as charged nature.

Polyphenols: Extraction Methods, Antioxidative Action, Bioavailability and Anticarcinogenic Effects

Being secondary plant metabolites, polyphenols represent a large and diverse group of substances abundantly present in a majority of fruits, herbs and vegetables. The current contribution is focused on their bioavailability, antioxidative and anticarcinogenic properties. An overview of extraction methods is also given, with supercritical fluid extraction highlighted as a promising eco-friendly alternative providing exceptional separation and protection from degradation of unstable polyphenols. The protective role of polyphenols against reactive oxygen and nitrogen species, UV light, plant pathogens, parasites and predators results in several beneficial biological activities giving rise to prophylaxis or possibly even to a cure for several prevailing human diseases, especially various cancer types. Omnipresence, specificity of the response and the absence of or low toxicity are crucial advantages of polyphenols as anticancer agents. The main problem represents their low bioavailability and rapid metabolism. One of the promising solutions lies in nanoformulation of polyphenols that prevents their degradation and thus enables significantly higher concentrations to reach the target cells. Another, more practiced, solution is the use of mixtures of various polyphenols that bring synergistic effects, resulting in lowering of the required therapeutic dose and in multitargeted action. The combination of polyphenols with existing drugs and therapies also shows promising results and significantly reduces their toxicity.

Drug repurposing to target Ebola virus replication and virulence using structural systems pharmacology

Background

The recent outbreak of Ebola has been cited as the largest in history. Despite this global health crisis, few drugs are available to efficiently treat Ebola infections. Drug repurposing provides a potentially efficient solution to accelerating the development of therapeutic approaches in response to Ebola outbreak. To identify such candidates, we use an integrated structural systems pharmacology pipeline which combines proteome-scale ligand binding site comparison, protein-ligand docking, and Molecular Dynamics (MD) simulation.

Results

One thousand seven hundred and sixty-six FDA-approved drugs and 259 experimental drugs were screened to identify those with the potential to inhibit the replication and virulence of Ebola, and to determine the binding modes with their respective targets. Initial screening has identified a number of promising hits. Notably, Indinavir; an HIV protease inhibitor, may be effective in reducing the virulence of Ebola. Additionally, an antifungal (Sinefungin) and several anti-viral drugs (e.g. Maraviroc, Abacavir, Telbivudine, and Cidofovir) may inhibit Ebola RNA-directed RNA polymerase through targeting the MTase domain.

Conclusions

Identification of safe drug candidates is a crucial first step toward the determination of timely and effective therapeutic approaches to address and mitigate the impact of the Ebola global crisis and future outbreaks of pathogenic diseases. Further in vitro and in vivo testing to evaluate the anti-Ebola activity of these drugs is warranted.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0941-9) contains supplementary material, which is available to authorized users.

Quantitative Correlation of Conformational Binding Enthalpy with Substrate Specificity of Serine Proteases

Members of the same protease family show different substrate specificity, even if they share identical folds, depending on the physiological processes they are part of. Here, we investigate the key factors for subpocket and global specificity of factor Xa, elastase, and granzyme B which despite all being serine proteases and sharing the chymotrypsin-fold show distinct substrate specificity profiles. We determined subpocket interaction potentials with GRID for static X-ray structures and an in silico generated ensemble of conformations. Subpocket interaction potentials determined for static X-ray structures turned out to be insufficient to explain serine protease specificity for all subpockets. Therefore, we generated conformational ensembles using molecular dynamics simulations. We identified representative binding site conformations using distance-based hierarchical agglomerative clustering and determined subpocket interaction potentials for each representative conformation of the binding site. Considering the differences in subpocket interaction potentials for these representative conformations as well as their abundance allowed us to quantitatively explain subpocket specificity for the nonprime side for all three example proteases on a molecular level. The methods to identify key regions determining subpocket specificity introduced in this study are directly applicable to other serine proteases, and the results provide starting points for new strategies in rational drug design.

Co-solvation effect on the binding mode of the α-mangostin/β-cyclodextrin inclusion complex

Cyclodextrins (CDs) have been extensively utilized as host molecules to enhance the solubility, stability and bioavailability of hydrophobic drug molecules through the formation of inclusion complexes. It was previously reported that the use of co-solvents in such studies may result in ternary (host:guest:co-solvent) complex formation. The objective of this work was to investigate the effect of ethanol as a co-solvent on the inclusion complex formation between α-mangostin (α-MGS) and β-CD, using both experimental and theoretical studies. Experimental phase-solubility studies were carried out in order to assess complex formation, with the mechanism of association being probed using a mathematical model. It was found that α-MGS was poorly soluble at low ethanol concentrations (0–10% v/v), but higher concentrations (10–40% v/v) resulted in better α-MGS solubility at all β-CD concentrations studied (0–10 mM). From the equilibrium constant calculation, the inclusion complex is still a binary complex (1:1), even in the presence of ethanol. The results from our theoretical study confirm that the binding mode is binary complex and the presence of ethanol as co-solvent enhances the solubility of α-MGS with some effects on the binding affinity with β-CD, depending on the concentration employed.

Reaction of pyranose dehydrogenase from Agaricus meleagris with its carbohydrate substrates

Monomeric Agaricus meleagris pyranose dehydrogenase (AmPDH) belongs to the glucose-methanol-choline family of oxidoreductases. An FAD cofactor is covalently tethered to His103 of the enzyme. AmPDH can double oxidize various mono- and oligosaccharides at different positions (C1 to C4). To study the structure/function relationship of selected active-site residues of AmPDH pertaining to substrate (carbohydrate) turnover in more detail, several active-site variants were generated, heterologously expressed in Pichia pastoris, and characterized by biochemical, biophysical and computational means. The crystal structure of AmPDH shows two active-site histidines, both of which could take on the role as the catalytic base in the reductive half-reaction. Steady-state kinetics revealed that His512 is the only catalytic base because H512A showed a reduction in (kcat /KM )glucose by a factor of 10(5) , whereas this catalytic efficiency was reduced by two or three orders of magnitude for His556 variants (H556A, H556N). This was further corroborated by transient-state kinetics, where a comparable decrease in the reductive rate constant was observed for H556A, whereas the rate constant for the oxidative half-reaction (using benzoquinone as substrate) was increased for H556A compared to recombinant wild-type AmPDH. Steady-state kinetics furthermore indicated that Gln392, Tyr510, Val511 and His556 are important for the catalytic efficiency of PDH. Molecular dynamics (MD) simulations and free energy calculations were used to predict d-glucose oxidation sites, which were validated by GC-MS measurements. These simulations also suggest that van der Waals interactions are the main driving force for substrate recognition and binding.

Entropy in bimolecular simulations: A comprehensive review of atomic fluctuations-based methods

Entropy of binding constitutes a major, and in many cases a detrimental, component of the binding affinity in biomolecular interactions. While the enthalpic part of the binding free energy is easier to calculate, estimating the entropy of binding is further more complicated. A precise evaluation of entropy requires a comprehensive exploration of the complete phase space of the interacting entities. As this task is extremely hard to accomplish in the context of conventional molecular simulations, calculating entropy has involved many approximations. Most of these golden standard methods focused on developing a reliable estimation of the conformational part of the entropy. Here, we review these methods with a particular emphasis on the different techniques that extract entropy from atomic fluctuations. The theoretical formalisms behind each method is explained highlighting its strengths as well as its limitations, followed by a description of a number of case studies for each method. We hope that this brief, yet comprehensive, review provides a useful tool to understand these methods and realize the practical issues that may arise in such calculations.

Prediction of cytochrome P450 mediated metabolism

Cytochrome P450 enzymes (CYPs) form one of the most important enzyme families involved in the metabolism of xenobiotics. CYPs comprise many isoforms, which catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be formed. However, it is often hard to rationalize what metabolites these enzymes generate. In recent years, many different in silico approaches have been developed to predict binding or regioselective product formation for the different CYP isoforms. These comprise ligand-based methods that are trained on experimental CYP data and structure-based methods that consider how the substrate is oriented in the active site or/and how reactive the part of the substrate that is accessible to the heme group is. We will review key aspects for various approaches that are available to predict binding and site of metabolism (SOM), what outcome can be expected from the predictions, and how they could potentially be improved.

Pi-pi Stacking Mediated Cooperative Mechanism for Human Cytochrome P450 3A4

Human Cytochrome P450 3A4 (CYP3A4) is an important member of the cytochrome P450 superfamily with responsibility for metabolizing ~50% of clinical drugs. Experimental evidence showed that CYP3A4 can adopt multiple substrates in its active site to form a cooperative binding model, accelerating substrate metabolism efficiency. In the current study, we constructed both normal and cooperative binding models of human CYP3A4 with antifungal drug ketoconazoles (KLN). Molecular dynamics simulation and free energy calculation were then carried out to study the cooperative binding mechanism. Our simulation showed that the second KLN in the cooperative binding model had a positive impact on the first one binding in the active site by two significant pi-pi stacking interactions. The first one was formed by Phe215, functioning to position the first KLN in a favorable orientation in the active site for further metabolism reactions. The second one was contributed by Phe304. This pi-pi stacking was enhanced in the cooperative binding model by the parallel conformation between the aromatic rings in Phe304 and the dioxolan moiety of the first KLN. These findings can provide an atomic insight into the cooperative binding in CYP3A4, revealing a novel pi-pi stacking mechanism for drug-drug interactions.

Transformation pathways of MeO-PBDEs catalyzed by active center of P450 enzymes: a DFT investigation employing 6-MeO-BDE-47 as a case

Recent in vivo and in vitro experiments indicated that methoxylated polybrominated diphenyl ethers (MeO-PBDEs) can be biotransformed into hydroxylated PBDEs (HO-PBDEs) that are more toxic than PBDEs and MeO-PBDEs. Nevertheless, the enzymatic transformation mechanism is not clear. We hypothesized that cytochrome P450 enzymes (CYPs) play a key role in the transformation and employed the density functional theory calculations to unveil the mechanism. The transformation of a model compound, 6-MeO-BDE-47, catalyzed by the active center of CYPs (Compound I), was computed. For the first time, our results show that the energy barriers for the addition of Compound I to the C atoms on the phenyl of 6-MeO-BDE-47 are much higher than that for hydroxylation of the methoxyl, indicating that O-demethylation is a dominating metabolic pathway. This is in line with experimental observations performed by others. The pathways for the transformation of 6-MeO-BDE-47 catalyzed by Compound I were clarified. A C-H bond of the methoxyl is activated by Compound I, followed by radical rebound to form carbinol intermediates, then the carbinols decompose to form 6-HO-BDE-47 with the assistance of water molecules. The computational method can be potentially employed to develop models that predict biotransformation of xenobiotics catalyzed by CYPs.