Application of Mixed-Solvent Molecular Dynamics Simulations for Prediction of Allosteric Sites on G Protein-Coupled Receptors.

The development of small molecule allosteric modulators acting at G protein-coupled receptors (GPCRs) is becoming increasingly attractive. Such compounds have advantages over traditional drugs acting at orthosteric sites on these receptors, in particular target specificity. However, the number and locations of druggable allosteric sites within most clinically relevant GPCRs are unknown. In the present study, we describe the development and application of a mixed-solvent molecular dynamics (MixMD)-based method for the identification of allosteric sites on GPCRs. The method employs small organic probes with druglike qualities to identify druggable hotspots in multiple replicate short-timescale simulations. As proof of principle, we first applied the method retrospectively to a test set of five GPCRs (cannabinoid receptor type 1, C-C chemokine receptor type 2, M2 muscarinic receptor, P2Y purinoceptor 1, and protease-activated receptor 2) with known allosteric sites in diverse locations. This resulted in the identification of the known allosteric sites on these receptors. We then applied the method to the µ-opioid receptor. Several allosteric modulators for this receptor are known, although the binding sites for these modulators are not known. The MixMD-based method revealed several potential allosteric sites on the mu-opioid receptor. Implementation of the MixMD-based method should aid future efforts in the structure-based drug design of drugs targeting allosteric sites on GPCRs. SIGNIFICANCE STATEMENT: Allosteric modulation of G protein-coupled receptors (GPCRs) has the potential to provide more selective drugs. However, there are limited structures of GPCRs bound to allosteric modulators, and obtaining such structures is problematic. Current computational methods utilize static structures and therefore may not identify hidden or cryptic sites. Here we describe the use of small organic probes and molecular dynamics to identify druggable allosteric hotspots on GPCRs. The results reinforce the importance of protein dynamics in allosteric site identification.

Molnupiravir; molecular and functional descriptors of mitochondrial safety.

Molnupiravir is an orally active nucleoside analog antiviral drug that recently was approved by the U.S. FDA for emergency treatment of adult patients infected with the SARS-CoV-2 (COVID-19) virus and at risk for severe progression. The active form of the drug, N-hydroxycytidine (NHC) triphosphate competes for incorporation by RNA-dependent RNA-polymerase (RdRp) into the replicating viral genome resulting in mutations and arrest of the replicating virus. Historically, some nucleoside analog antiviral drugs have been found to lack specificity for the virus and also inhibit replication and/or expression of the mitochondrial genome. The objective of the present study was to test whether molnupiravir and/or NHC also target mitochondrial DNA polymerase gamma (PolG) or RNA polymerase (POLRMT) activity to inhibit the replication and/or expression of the mitochondrial genome leading to impaired mitochondrial function. Human-derived HepG2 cells were exposed for 48 h in culture to increasing concentrations of either molnupiravir or NHC after which cytotoxicity, mtDNA copy number and mitochondrial gene expression were determined. The phenotypic endpoint, mitochondrial respiration, was measured with the Seahorse® XF96 Extracellular Flux Analyzer. Both molnupiravir and NHC were cytotoxic at concentrations of ≥10 µM. However, at non-cytotoxic concentrations, neither significantly altered mitochondrial gene dose or transcription, or mitochondrial respiration. From this we conclude that mitochondrial toxicity is not a primary off target in the mechanism of cytotoxicity for either molnupiravir or its active metabolite NHC in the HepG2 cell line.

Computational Identification of Possible Allosteric Sites and Modulators of the SARS-CoV-2 Main Protease.

In this study, we target the main protease (Mpro) of the SARS-CoV-2 virus as it is a crucial enzyme for viral replication. Herein, we report three plausible allosteric sites on Mpro that can expand structure-based drug discovery efforts for new Mpro inhibitors. To find these sites, we used mixed-solvent molecular dynamics (MixMD) simulations, an efficient computational protocol that finds binding hotspots through mapping the surface of unbound proteins with 5% cosolvents in water. We have used normal mode analysis to support our claim of allosteric control for these sites. Further, we have performed virtual screening against the sites with 361 hits from Mpro screenings available through the National Center for Advancing Translational Sciences (NCATS). We have identified the NCATS inhibitors that bind to the remote sites better than the active site of Mpro, and we propose these molecules may be allosteric regulators of the system. After identifying our sites, new X-ray crystal structures were released that show fragment molecules in the sites we found, supporting the notion that these sites are accurate and druggable.

Mixed-solvent molecular dynamics simulation-based discovery of a putative allosteric site on regulator of G protein signaling 4.

Regulator of G protein signaling 4 (RGS4) is an intracellular protein that binds to the Gα subunit ofheterotrimeric G proteins and aids in terminating G protein coupled receptor signaling. RGS4 has been implicated in pain, schizophrenia, and the control of cardiac contractility. Inhibitors of RGS4 have been developed but bind covalently to cysteine residues on the protein. Therefore, we sought to identify alternative druggable sites on RGS4 using mixed-solvent molecular dynamics simulations, which employ low concentrations of organic probes to identify druggable hotspots on the protein. Pseudo-ligands were placed in consensus hotspots, and perturbation with normal mode analysis led to the identification and characterization of a putative allosteric site, which would be invaluable for structure-based drug design of non-covalent, small molecule inhibitors. Future studies on the mechanism of this allostery will aid in the development of novel therapeutics targeting RGS4.

Remdesivir; molecular and functional measures of mitochondrial safety.

Remdesivir is one of a few antiviral drugs approved for treating severe cases of coronavirus 2 (SARS-CoV-2) infection in hospitalized patients. The prodrug is a nucleoside analog that interferes with viral replication by inhibiting viral RNA-dependent RNA polymerase. The drug has also been shown to be a weak inhibitor of human mitochondrial RNA polymerase, leaving open the possibility of mitochondrial off-targets and toxicity. The investigation was designed to explore whether remdesivir causes mitochondrial toxicity, using both genomic and functional parameters in the assessment. Human-derived HepG2 liver cells were exposed for up to 48 h in culture to increasing concentrations of remdesivir. At sub-cytotoxic concentrations (<1 µM), the drug failed to alter either the number of copies or the expression of the mitochondrial genome. mtDNA copy number was unaffected as was the relative rates of expression of mtDNA-encoded and nuclear encoded subunits of complexes I and IV of the mitochondrial respiratory chain. Consistent with this is the observation that remdesivir was without effect on mitochondrial respiration, including basal respiration, proton leak, maximum uncoupled respiration, spare respiratory capacity or coupling efficiency. We conclude that although remdesivir has weak inhibitory activity towards mitochondrial RNA polymerase, mitochondria are not primary off-targets for the mechanism of cytotoxicity of the drug.

Precise modelling and interpretation of bioactivities of ligands targeting G protein-coupled receptors.

MOTIVATION: Accurate prediction and interpretation of ligand bioactivities are essential for virtual screening and drug discovery. Unfortunately, many important drug targets lack experimental data about the ligand bioactivities; this is particularly true for G protein-coupled receptors (GPCRs), which account for the targets of about a third of drugs currently on the market. Computational approaches with the potential of precise assessment of ligand bioactivities and determination of key substructural features which determine ligand bioactivities are needed to address this issue. RESULTS: A new method, SED, was proposed to predict ligand bioactivities and to recognize key substructures associated with GPCRs through the coupling of screening for Lasso of long extended-connectivity fingerprints (ECFPs) with deep neural network training. The SED pipeline contains three successive steps: (i) representation of long ECFPs for ligand molecules, (ii) feature selection by screening for Lasso of ECFPs and (iii) bioactivity prediction through a deep neural network regression model. The method was examined on a set of 16 representative GPCRs that cover most subfamilies of human GPCRs, where each has 300-5000 ligand associations. The results show that SED achieves excellent performance in modelling ligand bioactivities, especially for those in the GPCR datasets without sufficient ligand associations, where SED improved the baseline predictors by 12% in correlation coefficient (r2) and 19% in root mean square error. Detail data analyses suggest that the major advantage of SED lies on its ability to detect substructures from long ECFPs which significantly improves the predictive performance. AVAILABILITY AND IMPLEMENTATION: The source code and datasets of SED are freely available at https://zhanglab.ccmb.med.umich.edu/SED/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Homologous G Protein-Coupled Receptors Boost the Modeling and Interpretation of Bioactivities of Ligand Molecules.

G protein-coupled receptors (GPCRs) are one of the most important drug targets, accounting for ∼34% of drugs on the market. For drug discovery, accurate modeling and explanation of bioactivities of ligands is critical for the screening and optimization of hit compounds. Homologous GPCRs are more likely to interact with chemically similar ligands, and they tend to share common binding modes with ligand molecules. The inclusion of homologous GPCRs in learning bioactivities of ligands potentially enhances the accuracy and interpretability of models due to utilizing increased training sample size and the existence of common ligand substructures that control bioactivities. Accurate modeling and interpretation of bioactivities of ligands by combining homologous GPCRs can be formulated as multitask learning with joint feature learning problem and naturally matched with the group lasso learning algorithm. Thus, we proposed a multitask regression learning with group lasso (MTR-GL) implemented by l2,1-norm regularization to model bioactivities of ligand molecules and then tested the algorithm on a series of thirty-five representative GPCRs datasets that cover nine subfamilies of human GPCRs. The results show that MTR-GL is overall superior to single-task learning methods and classic multitask learning with joint feature learning methods. Moreover, MTR-GL achieves better performance than state-of-the-art deep multitask learning based methods of predicting ligand bioactivities on most datasets (31/35), where MTR-GL obtained an average improvement of 38% on correlation coefficient (r2) and 29% on root-mean-square error over the DeepNeuralNet-QSAR predictors.

WDL-RF: predicting bioactivities of ligand molecules acting with G protein-coupled receptors by combining weighted deep learning and random forest.

Motivation: Precise assessment of ligand bioactivities (including IC50, EC50, Ki, Kd, etc.) is essential for virtual screening and lead compound identification. However, not all ligands have experimentally determined activities. In particular, many G protein-coupled receptors (GPCRs), which are the largest integral membrane protein family and represent targets of nearly 40% drugs on the market, lack published experimental data about ligand interactions. Computational methods with the ability to accurately predict the bioactivity of ligands can help efficiently address this problem. Results: We proposed a new method, WDL-RF, using weighted deep learning and random forest, to model the bioactivity of GPCR-associated ligand molecules. The pipeline of our algorithm consists of two consecutive stages: (i) molecular fingerprint generation through a new weighted deep learning method, and (ii) bioactivity calculations with a random forest model; where one uniqueness of the approach is that the model allows end-to-end learning of prediction pipelines with input ligands being of arbitrary size. The method was tested on a set of twenty-six non-redundant GPCRs that have a high number of active ligands, each with 200-4000 ligand associations. The results from our benchmark show that WDL-RF can generate bioactivity predictions with an average root-mean square error 1.33 and correlation coefficient (r2) 0.80 compared to the experimental measurements, which are significantly more accurate than the control predictors with different molecular fingerprints and descriptors. In particular, data-driven molecular fingerprint features, as extracted from the weighted deep learning models, can help solve deficiencies stemming from the use of traditional hand-crafted features and significantly increase the efficiency of short molecular fingerprints in virtual screening. Availability and implementation: The WDL-RF web server, as well as source codes and datasets of WDL-RF, is freely available at https://zhanglab.ccmb.med.umich.edu/WDL-RF/ for academic purposes. Supplementary information: Supplementary data are available at Bioinformatics online.

GLASS: a comprehensive database for experimentally validated GPCR-ligand associations.

MOTIVATION: G protein-coupled receptors (GPCRs) are probably the most attractive drug target membrane proteins, which constitute nearly half of drug targets in the contemporary drug discovery industry. While the majority of drug discovery studies employ existing GPCR and ligand interactions to identify new compounds, there remains a shortage of specific databases with precisely annotated GPCR-ligand associations. RESULTS: We have developed a new database, GLASS, which aims to provide a comprehensive, manually curated resource for experimentally validated GPCR-ligand associations. A new text-mining algorithm was proposed to collect GPCR-ligand interactions from the biomedical literature, which is then crosschecked with five primary pharmacological datasets, to enhance the coverage and accuracy of GPCR-ligand association data identifications. A special architecture has been designed to allow users for making homologous ligand search with flexible bioactivity parameters. The current database contains ∼500 000 unique entries, of which the vast majority stems from ligand associations with rhodopsin- and secretin-like receptors. The GLASS database should find its most useful application in various in silico GPCR screening and functional annotation studies. AVAILABILITY AND IMPLEMENTATION: The website of GLASS database is freely available at http://zhanglab.ccmb.med.umich.edu/GLASS/. CONTACT: zhng@umich.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The assembly of proline-rich membrane anchor (PRiMA)-linked acetylcholinesterase enzyme: glycosylation is required for enzymatic activity but not for oligomerization.

Acetylcholinesterase (AChE) anchors onto cell membranes by a transmembrane protein PRiMA (proline-rich membrane anchor) as a tetrameric form in vertebrate brain. The assembly of AChE tetramer with PRiMA requires the C-terminal "t-peptide" in AChE catalytic subunit (AChE(T)). Although mature AChE is well known N-glycosylated, the role of glycosylation in forming the physiologically active PRiMA-linked AChE tetramer has not been studied. Here, several lines of evidence indicate that the N-linked glycosylation of AChE(T) plays a major role for acquisition of AChE full enzymatic activity but does not affect its oligomerization. The expression of the AChE(T) mutant, in which all N-glycosylation sites were deleted, together with PRiMA in HEK293T cells produced a glycan-depleted PRiMA-linked AChE tetramer but with a much higher K(m) value as compared with the wild type. This glycan-depleted enzyme was assembled in endoplasmic reticulum but was not transported to Golgi apparatus or plasma membrane.

In silico analysis of SARS-CoV-2 proteins as targets for clinically available drugs.

The ongoing pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) requires treatments with rapid clinical translatability. Here we develop a multi-target and multi-ligand virtual screening method to identify FDA-approved drugs with potential activity against SARS-CoV-2 at traditional and understudied viral targets. 1,268 FDA-approved small molecule drugs were docked to 47 putative binding sites across 23 SARS-CoV-2 proteins. We compared drugs between binding sites and filtered out compounds that had no reported activity in an in vitro screen against SARS-CoV-2 infection of human liver (Huh-7) cells. This identified 17 "high-confidence", and 97 "medium-confidence" drug-site pairs. The "high-confidence" group was subjected to molecular dynamics simulations to yield six compounds with stable binding poses at their optimal target proteins. Three drugs-amprenavir, levomefolic acid, and calcipotriol-were predicted to bind to 3 different sites on the spike protein, domperidone to the Mac1 domain of the non-structural protein (Nsp) 3, avanafil to Nsp15, and nintedanib to the nucleocapsid protein involved in packaging the viral RNA. Our "two-way" virtual docking screen also provides a framework to prioritize drugs for testing in future emergencies requiring rapidly available clinical drugs and/or treating diseases where a moderate number of targets are known.

The genetic heterogeneity and drug resistance mechanisms of relapsed refractory multiple myeloma.

Multiple myeloma is the second most common hematological malignancy. Despite significant advances in treatment, relapse is common and carries a poor prognosis. Thus, it is critical to elucidate the genetic factors contributing to disease progression and drug resistance. Here, we carry out integrative clinical sequencing of 511 relapsed, refractory multiple myeloma (RRMM) patients to define the disease's molecular alterations landscape. The NF-κB and RAS/MAPK pathways are more commonly altered than previously reported, with a prevalence of 45-65% each. In the RAS/MAPK pathway, there is a long tail of variants associated with the RASopathies. By comparing our RRMM cases with untreated patients, we identify a diverse set of alterations conferring resistance to three main classes of targeted therapy in 22% of our cohort. Activating mutations in IL6ST are also enriched in RRMM. Taken together, our study serves as a resource for future investigations of RRMM biology and potentially informs clinical management.

The PRiMA-linked cholinesterase tetramers are assembled from homodimers: hybrid molecules composed of acetylcholinesterase and butyrylcholinesterase dimers are up-regulated during development of chicken brain.

Acetylcholinesterase (AChE) is anchored onto cell membranes by the transmembrane protein PRiMA (proline-rich membrane anchor) as a tetrameric globular form that is prominently expressed in vertebrate brain. In parallel, the PRiMA-linked tetrameric butyrylcholinesterase (BChE) is also found in the brain. A single type of AChE-BChE hybrid tetramer was formed in cell cultures by co-transfection of cDNAs encoding AChE(T) and BChE(T) with proline-rich attachment domain-containing proteins, PRiMA I, PRiMA II, or a fragment of ColQ having a C-terminal GPI addition signal (Q(N-GPI)). Using AChE and BChE mutants, we showed that AChE-BChE hybrids linked with PRiMA or Q(N-GPI) always consist of AChE(T) and BChE(T) homodimers. The dimer formation of AChE(T) and BChE(T) depends on the catalytic domains, and the assembly of tetramers with a proline-rich attachment domain-containing protein requires the presence of C-terminal "t-peptides" in cholinesterase subunits. Our results indicate that PRiMA- or ColQ-linked cholinesterase tetramers are assembled from AChE(T) or BChE(T) homodimers. Moreover, the PRiMA-linked AChE-BChE hybrids occur naturally in chicken brain, and their expression increases during development, suggesting that they might play a role in cholinergic neurotransmission.

Targeting acetylcholinesterase to membrane rafts: a function mediated by the proline-rich membrane anchor (PRiMA) in neurons.

In the mammalian brain, acetylcholinesterase (AChE) is anchored in cell membranes by a transmembrane protein PRiMA (proline-rich membrane anchor). We present evidence that at least part of the PRiMA-linked AChE is integrated in membrane microdomains called rafts. A significant proportion of PRiMA-linked AChE tetramers from rat brain was recovered in raft fractions; this proportion was markedly higher at low rather than at high concentrations of cold Triton X-100. The detergent-resistant fraction increased during brain development. In NG108-15 neuroblastoma cells transfected with cDNAs encoding AChE(T) and PRiMA, PRiMA-linked G(4) AChE was found in membrane rafts and showed the same sensitivity to cold Triton X-100 extraction as in the brain. The association of PRiMA-linked AChE with rafts was weaker than that of glycosylphosphatidylinositol-anchored G(2) AChE or G(4) Q(N)-H(C)-linked AChE. It was found to depend on the presence of a cholesterol-binding motif, called CRAC (cholesterol recognition/interaction amino acid consensus), located at the junction of transmembrane and cytoplasmic domains of both PRiMA I and II isoforms. The cytoplasmic domain of PRiMA, which differs between PRiMA I and PRiMA II, appeared to play some role in stabilizing the raft localization of G(4) AChE, because the Triton X-100-resistant fraction was smaller with the shorter PRiMA II isoform than that with the longer PRiMA I isoform.

Ontological modeling and analysis of experimentally or clinically verified drugs against coronavirus infection.

Our systematic literature collection and annotation identified 106 chemical drugs and 31 antibodies effective against the infection of at least one human coronavirus (including SARS-CoV, SAR-CoV-2, and MERS-CoV) in vitro or in vivo in an experimental or clinical setting. A total of 163 drug protein targets were identified, and 125 biological processes involving the drug targets were significantly enriched based on a Gene Ontology (GO) enrichment analysis. The Coronavirus Infectious Disease Ontology (CIDO) was used as an ontological platform to represent the anti-coronaviral drugs, chemical compounds, drug targets, biological processes, viruses, and the relations among these entities. In addition to new term generation, CIDO also adopted various terms from existing ontologies and developed new relations and axioms to semantically represent our annotated knowledge. The CIDO knowledgebase was systematically analyzed for scientific insights. To support rational drug design, a "Host-coronavirus interaction (HCI) checkpoint cocktail" strategy was proposed to interrupt the important checkpoints in the dynamic HCI network, and ontologies would greatly support the design process with interoperable knowledge representation and reasoning.

Virtual Screening of Human Class-A GPCRs Using Ligand Profiles Built on Multiple Ligand-Receptor Interactions.

G protein-coupled receptors (GPCRs) are a large family of integral membrane proteins responsible for cellular signal transductions. Identification of therapeutic compounds to regulate physiological processes is an important first step of drug discovery. We proposed MAGELLAN, a novel hierarchical virtual-screening (VS) pipeline, which starts with low-resolution protein structure prediction and structure-based binding-site identification, followed by homologous GPCR detections through structure and orthosteric binding-site comparisons. Ligand profiles constructed from the homologous ligand-GPCR complexes are then used to thread through compound databases for VS. The pipeline was first tested in a large-scale retrospective screening experiment against 224 human Class A GPCRs, where MAGELLAN achieved a median enrichment factor (EF) of 14.38, significantly higher than that using individual ligand profiles. Next, MAGELLAN was examined on 5 and 20 GPCRs from two public VS databases (DUD-E and GPCR-Bench) and resulted in an average EF of 9.75 and 13.70, respectively, which compare favorably with other state-of-the-art docking- and ligand-based methods, including AutoDock Vina (with EF = 1.48/3.16 in DUD-E and GPCR-Bench), DOCK 6 (2.12/3.47 in DUD-E and GPCR-Bench), PoLi (2.2 in DUD-E), and FINDSITECcomb2.0 (2.90 in DUD-E). Detailed data analyses show that the major advantage of MAGELLAN is attributed to the power of ligand profiling, which integrates complementary methods for ligand-GPCR interaction recognition and thus significantly improves the coverage and sensitivity of VS models. Finally, cases studies on opioid and motilin receptors show that new connections between functionally related GPCRs can be visualized in the minimum spanning tree built on the similarities of predicted ligand-binding ensembles, suggesting a novel use of MAGELLAN for GPCR deorphanization.

A conceptual framework for predicting the toxicity of reactive chemicals: modeling soft electrophilicity.

Although the literature is replete with QSAR models developed for many toxic effects caused by reversible chemical interactions, the development of QSARs for the toxic effects of reactive chemicals lacks a consistent approach. While limitations exit, an appropriate starting-point for modeling reactive toxicity is the applicability of the general rules of organic chemical reactions and the association of these reactions to cellular targets of importance in toxicology. The identification of plausible "molecular initiating events" based on covalent reactions with nucleophiles in proteins and DNA provides the unifying concept for a framework for reactive toxicity. This paper outlines the proposed framework for reactive toxicity. Empirical measures of the chemical reactivity of xenobiotics with a model nucleophile (thiol) are used to simulate the relative rates at which a reactive chemical is likely to bind irreversibly to cellular targets. These measures of intrinsic reactivity serve as correlates to a variety of toxic effects; what's more they appear to be more appropriate endpoints for QSAR modeling than the toxicity endpoints themselves.

Cumulative and irreversible cardiac mitochondrial dysfunction induced by doxorubicin.

Interference with mitochondrial calcium regulation is proposed to be a primary causative event in the mechanism of doxorubicin-induced cardiotoxicity. We previously reported disruption of mitochondrial calcium homeostasis after chronic doxorubicin administration (Solen et al. Toxicol. Appl. Pharmacol, 129: 214-222, 1994). The present study was designed to characterize the dose-dependent and cumulative interference with mitochondrial calcium regulation and to assess the reversibility of this functional lesion. Sprague Dawley rats were treated with 2 mg/kg/week doxorubicin s.c. for 4-8 weeks. With succinate as substrate, cardiac mitochondria isolated from rats after 4 weeks of treatment with doxorubicin expressed a lower calcium loading capacity compared with control. This suppression of calcium loading capacity increased with successive doses to 8 weeks of treatment (P < 0.05) and persisted for 5 weeks after the last doxorubicin injection, and was corroborated by dose-dependent and irreversible histopathological changes. Preincubation of mitochondria with tamoxifen, DTT, or monobromobimane did not reverse the diminished calcium loading capacity caused by doxorubicin. In contrast, incubation with cyclosporin A abolished any discernible difference in mitochondrial calcium loading capacity between doxorubicin-treated and saline-treated rats. The decrease in cardiac mitochondrial calcium loading capacity was not attributable to bioenergetic changes in the electron transport chain, because the mitochondrial coupling efficiency was not altered by doxorubicin treatment. However, the ADP/ATP translocase content was significantly lower in mitochondria from rats that received 8 weeks of doxorubicin treatment. These data indicate that doxorubicin treatment in vivo causes a dose-dependent and irreversible decrease in mitochondrial calcium loading capacity. Suppression of adenine nucleotide translocase content may be a key factor altering the calcium-dependent regulation of the mitochondrial permeability transition pore, which may account for the cumulative and irreversible loss of myocardial function in patients receiving doxorubicin chemotherapy.

Preferential oxidation of cardiac mitochondrial DNA following acute intoxication with doxorubicin.

The purpose of this investigation was to determine whether acute doxorubicin intoxication causes a preferential accumulation of 8-hydroxydeoxyguanosine (8OHdG) adducts to mitochondrial DNA (mtDNA) as opposed to nuclear DNA (nDNA), particularly in cardiac tissue. Adult male rats received a single i.p. bolus of doxorubicin (15 mg/kg) and were killed 1-14 days later. Acute intoxication with doxorubicin caused a 2-fold greater increase in 8OHdG adducts to mtDNA compared to nDNA, the concentration of adducts to both nDNA and mtDNA being 20%-40% greater for heart as opposed to liver. For both tissues, the relative abundance of adducts was highest at the earliest time-point examined (24 h) and decreased to control values by 2 weeks. The temporal dilution of 8OHdG adducts was not the result of cell hyperplasia and was only partially due to amplification of the mitochondrial genome, most probably via an increase in DNA copy number rather than a stimulation of mitochondrial biogenesis.

Cardioselective and cumulative oxidation of mitochondrial DNA following subchronic doxorubicin administration.

We recently reported the preferential accumulation of 8-hydroxydeoxyguanosine (8OHdG) adducts in cardiac mitochondrial DNA (mtDNA) following acute intoxication of rats with doxorubicin (C.M. Palmeira et al., Biochim. Biophys. Acta, 1321 (1997) 101-106). The concentration of 8OHdG adducts decreased to control values within 2 weeks. Since conventional antineoplastic therapy entails repeated administration of small doses of doxorubicin, it was of interest to characterize the kinetics for the accumulation and repair of 8OHdG adducts in the various DNA fractions. Weekly injections of doxorubicin (2 mg/kg, i.p.) to adult male Sprague-Dawley rats caused a cumulative dose-dependent increase in the concentration of 8OHdG adducts in both mtDNA and nuclear DNA (nDNA) from heart and liver. Following six weekly injections, the concentration of 8OHdG in cardiac mtDNA was 50% higher than liver mtDNA and twice that of cardiac nDNA. In contrast to the rapid repair of 8OHdG observed during the first days following an acute intoxicating dose of doxorubicin, the concentration of 8OHdG adducts remained constant between 1 and 5 weeks following the last injection. This was true for all DNA fractions examined. The cardioselective accumulation and persistence of 8OHdG adducts to mtDNA is consistent with the implication of mitochondrial dysfunction in the cumulative and irreversible cardiotoxicity observed clinically in patients receiving doxorubicin cancer chemotherapy.