The role of the circadian timing system on drug metabolism and detoxification: an update

INTRODUCTION

The 24-hour variations in drug absorption, distribution, metabolism, and elimination, collectively known as pharmacokinetics, are fundamentally influenced by rhythmic physiological processes regulated by the molecular clock. Recent advances have elucidated the intricacies of the circadian timing system and the molecular interplay between biological clocks, enzymes and transporters in preclinical level.

AREA COVERED

Circadian rhythm of the drug metabolizing enzymes and carrier efflux functions possess a major role for drug metabolism and detoxification. The efflux and metabolism function of intestines and liver seems important. The investigations revealed that the ABC and SLC transporter families, along with cytochrome p-450 systems in the intestine, liver, and kidney, play a dominant role in the circadian detoxification of drugs. Additionally, the circadian control of efflux by the blood-brain barrier is also discussed.

EXPERT OPINION

The influence of the circadian timing system on drug pharmacokinetics significantly impacts the efficacy, adverse effects, and toxicity profiles of various drugs. Moreover, the emergence of sex-related circadian changes in the metabolism and detoxification processes has underscored the importance of considering gender-specific differences in drug tolerability and pharmacology. A better understanding of coupling between central clock and circadian metabolism/transport contributes to the development of more rational drug utilization and the implementation of chronotherapy applications.

A cartridge based sensor array platform for multiple coagulation measurements from plasma

This paper proposes a MEMS-based sensor array enabling multiple clot-time tests for plasma in one disposable microfluidic cartridge. The versatile LoC (Lab-on-Chip) platform technology is demonstrated here for real-time coagulation tests (activated Partial Thromboplastin Time (aPTT) and Prothrombin Time (PT)). The system has a reader unit and a disposable cartridge. The reader has no electrical connections to the cartridge. This enables simple and low-cost cartridge designs and avoids reliability problems associated with electrical connections. The cartridge consists of microfluidic channels and MEMS microcantilevers placed in each channel. The microcantilevers are made of electroplated nickel. They are actuated remotely using an external electro-coil and the read-out is also conducted remotely using a laser. The phase difference between the cantilever oscillation and the coil drive is monitored in real time. During coagulation, the viscosity of the blood plasma increases resulting in a change in the phase read-out. The proposed assay was tested on human and control plasma samples for PT and aPTT measurements. PT and aPTT measurements from control plasma samples are comparable with the manufacturer's datasheet and the commercial reference device. The measurement system has an overall 7.28% and 6.33% CV for PT and aPTT, respectively. For further implementation, the microfluidic channels of the cartridge were functionalized for PT and aPTT tests by drying specific reagents in each channel. Since simultaneous PT and aPTT measurements are needed in order to properly evaluate the coagulation system, one of the most prominent features of the proposed assay is enabling parallel measurement of different coagulation parameters. Additionally, the design of the cartridge and the read-out system as well as the obtained reproducible results with 10 μl of the plasma samples suggest an opportunity for a possible point-of-care application.

Enhanced heterotetrameric assembly of potato ADP-glucose pyrophosphorylase using reverse genetics

ADP-glucose pyrophosphorylase (AGPase) is a key allosteric enzyme in plant starch biosynthesis. Plant AGPase is a heterotetrameric enzyme that consists of large (LS) and small subunits (SS), which are encoded by two different genes. Computational and experimental studies have revealed that the heterotetrameric assembly of AGPase is thermodynamically weak. Modeling studies followed by the mutagenesis of the LS of the potato AGPase identified a heterotetramer-deficient mutant, LS(R88A). To enhance heterotetrameric assembly, LS(R88A) cDNA was subjected to error-prone PCR, and second-site revertants were identified according to their ability to restore glycogen accumulation, as assessed with iodine staining. Selected mutations were introduced into the wild-type (WT) LS and co-expressed with the WT SS in Escherichia coli glgC(-). The biochemical characterization of revertants revealed that LS(I90V)SS(WT), LS(Y378C)SS(WT) and LS(D410G)SS(WT) mutants displayed enhanced heterotetrameric assembly with the WT SS. Among these mutants, LS(Y378C)SS(WT) AGPase displayed increased heat stability compared with the WT enzyme. Kinetic characterization of the mutants indicated that the LS(I90V)SS(WT) and LS(Y378C)SS(WT) AGPases have comparable allosteric and kinetic properties. However, the LS(D410G)SS(WT) mutant exhibited altered allosteric properties of being less responsive and more sensitive to 3-phosphoglyceric acid activation and inorganic phosphate inhibition. This study not only enhances our understanding of the interaction between the SS and the LS of AGPase but also enables protein engineering to obtain enhanced assembled heat-stable variants of AGPase, which can be used for the improvement of plant yields.

MEMS biosensor for detection of Hepatitis A and C viruses in serum

Resonant microcantilever arrays are developed for the purpose of label-free and real-time analyte monitoring and biomolecule detection. MEMS cantilevers made of electroplated nickel are functionalized with Hepatitis antibodies. Hepatitis A and C antigens at different concentrations are introduced in undiluted bovine serum. All preparation and measurement steps are carried out in the liquid within a specifically designed flowcell without ever drying the cantilevers throughout the experiment. Both actuation and sensing are done remotely and therefore the MEMS cantilevers have no electrical connections, allowing for easily disposable sensor chips. Actuation is achieved using an electromagnet and the interferometric optical sensing is achieved using laser illumination and embedded diffraction gratings at the tip of each cantilever. Resonant frequency of the cantilevers in dynamic motion is monitored using a self-sustaining closed-loop control circuit and a frequency counter. Specificity is demonstrated by detecting both Hepatitis A and Hepatitis C antigens and their negative controls. This is the first report of Hepatitis antigen detection by resonant cantilevers exposed to undiluted serum. A dynamic range in excess of 1000 and with a minimum detectable concentration limit of 0.1ng/ml (1.66pM) is achieved for both Hepatitis A and C. This result is comparable to labeled detection methods such as ELISA.

Optimization based tumor classification from microarray gene expression data

BACKGROUND

An important use of data obtained from microarray measurements is the classification of tumor types with respect to genes that are either up or down regulated in specific cancer types. A number of algorithms have been proposed to obtain such classifications. These algorithms usually require parameter optimization to obtain accurate results depending on the type of data. Additionally, it is highly critical to find an optimal set of markers among those up or down regulated genes that can be clinically utilized to build assays for the diagnosis or to follow progression of specific cancer types. In this paper, we employ a mixed integer programming based classification algorithm named hyper-box enclosure method (HBE) for the classification of some cancer types with a minimal set of predictor genes. This optimization based method which is a user friendly and efficient classifier may allow the clinicians to diagnose and follow progression of certain cancer types.

METHODOLOGY/PRINCIPAL FINDINGS

We apply HBE algorithm to some well known data sets such as leukemia, prostate cancer, diffuse large B-cell lymphoma (DLBCL), small round blue cell tumors (SRBCT) to find some predictor genes that can be utilized for diagnosis and prognosis in a robust manner with a high accuracy. Our approach does not require any modification or parameter optimization for each data set. Additionally, information gain attribute evaluator, relief attribute evaluator and correlation-based feature selection methods are employed for the gene selection. The results are compared with those from other studies and biological roles of selected genes in corresponding cancer type are described.

CONCLUSIONS/SIGNIFICANCE

The performance of our algorithm overall was better than the other algorithms reported in the literature and classifiers found in WEKA data-mining package. Since it does not require a parameter optimization and it performs consistently very high prediction rate on different type of data sets, HBE method is an effective and consistent tool for cancer type prediction with a small number of gene markers.

Classification of cytochrome P450 inhibitors with respect to binding free energy and pIC50 using common molecular descriptors

Virtual screening of chemical libraries following experimental assays of drug candidates is a common procedure in structure based drug discovery. However, the relationship between binding free energies and biological activities (pIC50) of drug candidates is still an unsolved issue that limits the efficiency and speed of drug development processes. In this study, the relationship between them is investigated based on a common molecular descriptor set for human cytochrome P450 enzymes (CYPs). CYPs play an important role in drug-drug interactions, drug metabolism, and toxicity. Therefore, in silico prediction of CYP inhibition by drug candidates is one of the major considerations in drug discovery. The combination of partial least-squares regression (PLSR) and a variety of classification algorithms were employed by considering this relationship as a classification problem. Our results indicate that PLSR with classification is a powerful tool to predict more than one output such as binding free energy and pIC50 simultaneously. PLSR with mixed-integer linear programming based hyperboxes predicts the binding free energy and pIC50 with a mean accuracy of 87.18% (min: 81.67% max: 97.05%) and 88.09% (min: 79.83% max: 92.90%), respectively, for the cytochrome p450 superfamily using the common 6 molecular descriptors with a 10-fold cross-validation.

Classification of drug molecules considering their IC50 values using mixed-integer linear programming based hyper-boxes method

BACKGROUND

A priori analysis of the activity of drugs on the target protein by computational approaches can be useful in narrowing down drug candidates for further experimental tests. Currently, there are a large number of computational methods that predict the activity of drugs on proteins. In this study, we approach the activity prediction problem as a classification problem and, we aim to improve the classification accuracy by introducing an algorithm that combines partial least squares regression with mixed-integer programming based hyper-boxes classification method, where drug molecules are classified as low active or high active regarding their binding activity (IC50 values) on target proteins. We also aim to determine the most significant molecular descriptors for the drug molecules.

RESULTS

We first apply our approach by analyzing the activities of widely known inhibitor datasets including Acetylcholinesterase (ACHE), Benzodiazepine Receptor (BZR), Dihydrofolate Reductase (DHFR), Cyclooxygenase-2 (COX-2) with known IC50 values. The results at this stage proved that our approach consistently gives better classification accuracies compared to 63 other reported classification methods such as SVM, Naïve Bayes, where we were able to predict the experimentally determined IC50 values with a worst case accuracy of 96%. To further test applicability of this approach we first created dataset for Cytochrome P450 C17 inhibitors and then predicted their activities with 100% accuracy.

CONCLUSION

Our results indicate that this approach can be utilized to predict the inhibitory effects of inhibitors based on their molecular descriptors. This approach will not only enhance drug discovery process, but also save time and resources committed.

Purification and characterization of a type III photolyase from Caulobacter crescentus

The photolyase/cryptochrome family is a large family of flavoproteins that encompasses DNA repair proteins, photolyases, and cryptochromes that regulate blue-light-dependent growth and development in plants, and light-dependent and light-independent circadian clock setting in animals. Phylogenetic analysis has revealed a new class of the family, named type III photolyase, which cosegregates with plant cryptochromes. Here we describe the isolation and characterization of a type III photolyase from Caulobacter crescentus. Spectroscopic analysis shows that the enzyme contains both the methenyl tetrahydrofolate photoantenna and the FAD catalytic cofactor. Biochemical analysis shows that it is a bona fide photolyase that repairs cyclobutane pyrimidine dimers. Mutation of an active site Trp to Arg disrupts FAD binding with no measurable effect on MTHF binding. Using enzyme preparations that contain either both chromophores or only folate, we were able to determine the efficiency and rate of transfer of energy from MTHF to FAD.

Ultrafast dynamics of resonance energy transfer in cryptochrome

In this communication, we report the ultrafast dynamics of resonance energy transfer in a blue-light photoreceptor, Vibrio cholerae cryptochrome. The transfer was observed to occur in 60 ps. We also studied the local rigidity and solvation around the binding site of the photoantenna molecule. The results for the first time show energy transfer in cryptochrome suggesting some mechanistic similarities between photolyase that repairs damaged DNA and cryptochrome that mediates blue-light signaling.

Analysis of the role of intraprotein electron transfer in photoreactivation by DNA photolyase in vivo

Escherichia coli DNA photolyase contains FADH(-) as the catalytic cofactor. The cofactor becomes oxidized to the FADH(*) blue neutral radical during purification. The E-FADH(*) form of the enzyme is catalytically inert but can be converted to the active E-FADH(-) form by a photoreduction reaction that involves intraprotein electron transfer from Trp306. It is thought that the E-FADH(*) form is also transiently generated during pyrimidine dimer repair by photoinduced electron transfer, and it has been suggested that the FADH(*) that is generated after each round of catalysis must be photoreduced before the enzyme can engage in subsequent rounds of repair. In this study, we introduced the Trp306Phe mutation into the chromosomal gene and tested the non-photoreducible W306F mutant for photorepair in vivo. We find that both wild-type and W306F mutant photolyases carry out at least 25 rounds of photorepair at the same rate. We conclude that photoreduction by intraprotein electron transfer is not part of the photolyase photocycle under physiological conditions.

Purification and characterization of three members of the photolyase/cryptochrome family blue-light photoreceptors from Vibrio cholerae

The sequence of Vibrio cholerae genome revealed three genes belonging to the photolyase/cryptochrome blue-light photoreceptor family. The proteins encoded by the three genes were purified and characterized. All three proteins contain folate and flavin cofactors and have absorption peaks in the range of 350-500 nm. Only one of the three, VcPhr, is a photolyase specific for cyclobutane pyrimidine dimers. The other two are cryptochromes and were designated VcCry1 and VcCry2, respectively. Mutation of phr abolishes photoreactivation of UV-induced killing, whereas mutations in cry1 and cry2 do not affect photorepair activity. VcCry1 exhibits some unique features. Of all cryptochromes characterized to date, it is the only one that contains stoichiometric amounts of both chromophores and retains its flavin cofactor in the two-electron reduced FADH2 form. In addition, VcCry1 exhibits RNA binding activity and co-purifies with an RNA of 60-70 nucleotides in length.

Generation, characterization, and heterologous expression of wild-type and up-regulated forms of Arabidopsis thaliana leaf ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase (AGPase), a key enzyme in starch biosynthesis of higher plants, consists of a pair of regulatory large (LS) and catalytically small (SS) subunits. In plants, these subunits are coded by multiple genes resulting in the formation of tissue-specific enzyme forms, which are differentially regulated during plant growth and development. Some AGPase isoforms differ in catalytic and regulatory properties as well as intracellular location. In an effort to gain a better understanding of the role of the leaf AGPase in carbon partitioning and its effect on plant productivity, the Arabidopsis leaf AGPase containing the mature forms of the SS and LS was expressed in a heterologous expression system and characterized enzymatically. The Arabidopsis recombinant AGPase had kinetic values for 3-phosphoglyceric acid, glucose-1-phosphate and Mg(2+) similar to those of the native enzyme. As the N-terminus of the LS has been suggested to be involved in enzyme function, the length of the N-terminal region was extended or shortened. Of the five modified LSs analyzed, only the T5 form lacking six residues of the mature N-terminus was able to form detectable levels of enzyme activity, indicating that the N-terminal region is critical for enzyme function. Two up-regulatory LS mutations that allosterically activate the potato enzyme, a stem isoform, were introduced into the corresponding Arabidopsis LS sequences and co-expressed with wild-type SS. Both modified enzymes showed up-regulatory properties, indicating that these specific residue changes were also operational in the leaf isoform.

Directed molecular evolution of ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase catalyzes a rate-limiting reaction in prokaryotic glycogen and plant starch biosynthesis. Despite sharing similar molecular size and catalytic and allosteric regulatory properties, the prokaryotic and higher plant enzymes differ in higher-order protein structure. The bacterial enzyme is encoded by a single gene whose product of ca. 50,000 Da assembles into a homotetrameric structure. Although the higher plant enzyme has a similar molecular size, it is made up of a pair of large subunits and a pair of small subunits, encoded by different genes. To identify the basis for the evolution of AGPase function and quaternary structure, a potato small subunit homotetrameric mutant, TG-15, was subjected to iterations of DNA shuffling and screened for enzyme variants with up-regulated catalytic and/or regulatory properties. A glycogen selection/screening regimen of buoyant density gradient centrifugation and iodine vapor colony staining on glucose-containing media was used to increase the stringency of selection. This approach led to the isolation of a population of AGPase small subunit homotetramer enzymes with enhanced affinity toward ATP and increased sensitivity to activator and/or greater resistance to inhibition than TG-15. Several enzymes displayed a shift in effector preference from 3-phosphoglycerate to fructose-6 phosphate or fructose-1,6-bis-phosphate, effectors used by specific bacterial AGPases. Our results suggest that evolution of AGPase, with regard to quaternary structure, allosteric effector selectivity, and effector sensitivity, can occur through the introduction of a few point mutations alone with low-level recombination hastening the process.

Analysis of allosteric effector binding sites of potato ADP-glucose pyrophosphorylase through reverse genetics

ADP-glucose pyrophosphorylase (AGPase) is a key regulatory enzyme of bacterial glycogen and plant starch synthesis as it controls carbon flux via its allosteric regulatory behavior. Unlike the bacterial enzyme that is composed of a single subunit type, the plant AGPase is a heterotetrameric enzyme (alpha2beta2) with distinct roles for each subunit type. The large subunit (LS) is involved mainly in allosteric regulation through its interaction with the catalytic small subunit (SS). The LS modulates the catalytic activity of the SS by increasing the allosteric regulatory response of the hetero-oligomeric enzyme. To identify regions of the LS involved in binding of effector molecules, a reverse genetics approach was employed. A potato (Solanum tuberosum L.) AGPase LS down-regulatory mutant (E38A) was subjected to random mutagenesis using error-prone polymerase chain reaction and screened for the capacity to form an enzyme capable of restoring glycogen production in glgC(-) Escherichia coli. Dominant mutations were identified by their capacity to restore glycogen production when the LS containing only the second site mutations was co-expressed with the wild-type SS. Sequence analysis showed that most of the mutations were decidedly nonrandom and were clustered at conserved N- and C-terminal regions. Kinetic analysis of the dominant mutant enzymes indicated that the K(m) values for cofactor and substrates were comparable with the wild-type AGPase, whereas the affinities for activator and inhibitor were altered appreciably. These AGPase variants displayed increased resistance to P(i) inhibition and/or greater sensitivity toward 3-phosphoglyceric acid activation. Further studies of Lys-197, Pro-261, and Lys-420, residues conserved in AGPase sequences, by site-directed mutagenesis suggested that the effectors 3-phosphoglyceric acid and P(i) interact at two closely located binding sites.

Investigation of subunit function in ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase (AGPase), a key regulatory enzyme in higher plant starch biosynthesis, is composed of a pair of large and small subunits (alpha(2)beta(2)). Current evidence suggests that the large subunit has primarily a regulatory function, while the small subunit has both regulatory and catalytic roles. To define the structure-function relationship of the large subunit (LS), the LS of potato AGPase was subjected to chemical mutagenesis and coexpressed with the wild-type (WT) small subunit (SS) cDNA in an AGPase defective Escherichia coli strain. An LS mutant (M143) was isolated, which accumulated very low levels of glycogen compared to the WT recombinant AGPase, but maintained normal catalytic activity when assayed under saturating conditions. Sequence analysis revealed that M143 has a single amino acid change, V463I, which lies adjacent to the C-terminus. This single mutation had no effect on the Km for ATP and Mg(2+), which were similar to the WT enzyme. The K(m) for glucose 1-P, however, was sixfold higher than the WT enzyme. These results suggest that the LS plays a role in binding glucose 1-P through its interaction with the SS.

Transcriptional expression characteristics and subcellular localization of ADP-glucose pyrophosphorylase in the oil plant Perilla frutescens

Three ADP-glucose pyrophosphorylase clones were isolated from the cotyledon cDNA library of the oil plant, Perilla frutescens, and their intracellular localization investigated. Two of three cDNAs (PfagpS1 and PfagpS2) were homologous to the catalytic small subunit of AGPases found in other plants, while the third clone (PfagpL) was highly similar to the large subunit type. Transcripts for PfagpS1 and PfagpS2 were observed in both photosynthetic and non-photosynthetic tissue, showing the highest expression in the stem, while PfagpL transcripts were abundantly expressed in stem and cotyledon. To evaluate the subcellular localization of PfagpS2 and PfagpL as well as the maize BT2, N-terminus-GFP DNA fusion were constructed and transformed into tobacco plants. Immunoblot analysis showed that the expressed PfagpS2- and PfagpL-GFP fusions were targeted to the plastid in the heterologous tobacco system whereas the BT2-GFP remained intact, suggesting a cytoplasmic location. These intracellular assignments were confirmed by direct confocal microscopic examination. GFP signals were localized to the cytoplasm as well as in the nucleus in BT2-GFP plants, and to the plastids in PfagpS2- and PfagpL-GFP plants. Our results indicate that Perilla cotyledons contain multiple AGPase subunits, of which at least two isoforms and very likely the third, are plastidial in nature.

Isolation and characterization of a higher plant ADP-glucose pyrophosphorylase small subunit homotetramer

ADP-glucose pyrophosphorylase (AGPase) is the allosterically regulated gateway for carbon entry into transient and storage starch in plants as well as glycogen in bacteria. This enzyme plays a key role in the modulation of photosynthetic efficiency in source tissues and directly determines the level of storage starch in sink tissues, thus influencing overall crop yield potential. AGPase is a tetrameric enzyme; in higher plants it consists of two regulatory large subunits (LS) and two catalytic small subunits (SS), while in cyanobacteria and prokaryotes the enzyme is homotetrameric. The potato SS gene in pML10 was mutated by hydroxylamine and mutants were screened for elevated homotetrameric activity by iodine vapor staining. This search strategy led to the isolation of SS mutants (SUP-1, TG-15) that had pyrophosphorylase activity in the absence of the LS. TG-15 has a leucine to phenylalanine change at position 48 (L(48)F) that corresponds to a phenylalanine residue at the analogous position in the Escherichia coli homotetrameric AGPase as well as a valine to isoleucine change at position 59 (V(59)I). TG-15 was partially purified and kinetic analysis revealed substrate and effector affinities equal to wild type heterotetrameric enzyme with the exception of ATP binding.

Engineering starch for increased quantity and quality

The characterization and production of starch variants from mutation studies and transgene technology has been invaluable for our understanding of the synthesis of the starch granule. The knowledge gained has allowed for genetic manipulation of the starch biosynthetic pathway in plants. This in vivo approach can be used to generate novel starches and diminishes the need for post-harvest chemically and enzymatically treated starches. Thus, the modification of the starch biosynthetic pathway is a plausible means by which starches with novel properties and applications can be created.

Generation of up-regulated allosteric variants of potato ADP-glucose pyrophosphorylase by reversion genetics

Mutagenesis of the large subunit (LS) of the potato ADP-glucose pyrophosphorylase generated an enzyme, P52L, that was insensitive to 3-phosphoglycerate (3-PGA). To identify additional residues involved in 3-PGA interaction, we subjected P52L LS DNA to a second round of mutagenesis and identified second-site revertants by their ability to restore glycogen accumulation as assessed by iodine (I2) staining. Enzymes from class I revertants with normal I2-staining had an 11- to 49-fold greater affinity for the activator 3-PGA compared with the P52L mutant and a decreased sensitivity to the inhibitor orthophosphate. Sequence analysis of these class I revertants identified a P66L mutation in R4, an E38K mutation in R20, and a G101N mutation in R10 and R32. These mutations appear to restore 3-PGA binding by counteracting the effect of the P52L mutation because introducing E38K or G101N into the wild-type LS led to enzyme variants with higher affinity for the activator 3-PGA and increased resistance to the inhibitor orthophosphate. The generation of these revertant enzymes provides additional structure-function information on the allosteric regulation of higher plant ADP-glucose pyrophosphorylases and validates a strategy for developing novel variants of the enzyme that may be useful in manipulating starch biosynthesis in higher plants.