Identifying cytochrome p450 functional networks and their allosteric regulatory elements

Predicting phase-I metabolism of piceatannol: an in silico study

Piceatannol is a natural compound found in plants and can be derived from resveratrol. While resveratrol has been extensively researched for its effects and how the body processes it, there are concerns about its use. These concerns include its limited absorption in the body, the need for specific dosages, potential interactions with other drugs, lack of standardization, and limited clinical evidence to support its benefits. Interestingly, Piceatannol, another compound derived from resveratrol, has received less attention from researchers but appears to offer advantages. It has better bioavailability and seems to have a more favorable therapeutic profile compared to resveratrol. Surprisingly, no previous attempts have been made to explore or predict the metabolites of piceatannol when it interacts with the enzyme cytochrome P450. This study aims to fill that gap by predicting how piceatannol is metabolized by cytochrome P450 and assessing any potential toxicity associated with its metabolites. This research is interesting because it's the first of its kind to investigate the metabolic fate of piceatannol, especially in the context of cytochrome P450. The findings have the potential to significantly contribute to the field of piceatannol research, particularly in the food industry where this compound has applications and implications.

Graphical abstract

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.

Exploring the Interactome of Cytochrome P450 2E1 in Human Liver Microsomes with Chemical Crosslinking Mass Spectrometry

Aiming to elucidate the system-wide effects of the alcohol-induced increase in the content of cytochrome P450 2E1 (CYP2E1) on drug metabolism, we explored the array of its protein-protein interactions (interactome) in human liver microsomes (HLM) with chemical crosslinking mass spectrometry (CXMS). Our strategy employs membrane incorporation of purified CYP2E1 modified with photoreactive crosslinkers benzophenone-4-maleimide and 4-(N-succinimidylcarboxy)benzophenone. Exposure of bait-incorporated HLM samples to light was followed by isolating the His-tagged bait protein and its crosslinked aggregates on Ni-NTA agarose. Analyzing the individual bands of SDS-PAGE slabs of thereby isolated protein with the toolset of untargeted proteomics, we detected the crosslinked dimeric and trimeric complexes of CYP2E1 with other drug-metabolizing enzymes. Among the most extensively crosslinked partners of CYP2E1 are the cytochromes P450 2A6, 2C8, 3A4, 4A11, and 4F2, UDP-glucuronosyltransferases (UGTs) 1A and 2B, fatty aldehyde dehydrogenase (ALDH3A2), epoxide hydrolase 1 (EPHX1), disulfide oxidase 1α (ERO1L), and ribophorin II (RPN2). These results demonstrate the exploratory power of the proposed CXMS strategy and corroborate the concept of tight functional integration in the human drug-metabolizing ensemble through protein-protein interactions of the constituting enzymes.

The hydrogen-bond configuration modulates the energy transfer efficiency in helical protein nanotubes

Energy transport in proteins is critical to a variety of physical, chemical, and biological processes in living organisms. While strenuous efforts have been made to study vibrational energy transport in proteins, thermal transport processes across the most fundamental building blocks of proteins, i.e. helices, are not well understood. This work studies energy transport in a group of "isomer" helices. The π-helix is shown to have the highest thermal conductivity, 110% higher than that of the α-helix and 207% higher than that of the 310-helix. The H-bond connectivity is found to govern thermal transport mechanisms including the phonon spectral energy density, dispersion, mode-specific transport, group velocity, and relaxation time. The energy transport is strongly correlated with the H-bond strength which is also modulated by the H-bond connectivity. These fundamental insights provide a novel perspective for understanding energy transfer in proteins and guiding a rational molecule-level design of novel materials with configurable H-bonds.

Changes and sex- and age-related differences in the expression of drug metabolizing enzymes in a KRAS-mutant mouse model of lung cancer

Background

This study aimed to systematically profile the alterations and sex- and age-related differences in the drug metabolizing enzymes (DMEs) in a KRAS-mutant mouse model of lung cancer (KRAS mice).

Methodology

In this study, the LC-MS/MS approach and a probe substrate method were used to detect the alterations in 21 isoforms of DMEs, as well as the enzymatic activities of five isoforms, respectively. Western blotting was applied to study the protein expression of four related receptors.

Results

The proteins contents of CYP2C29 and CYP3A11, were significantly downregulated in the livers of male KRAS mice at 26 weeks (3.7- and 4.4-fold, respectively, p < 0.05). SULT1A1 and SULT1D1 were upregulated by 1.8- to 7.0- fold at 20 (p = 0.015 and 0.017, respectively) and 26 weeks (p = 0.055 and 0.031, respectively). There were positive correlations between protein expression and enzyme activity for CYP2E1, UGT1A9, SULT1A1 and SULT1D1 (r² ≥ 0.5, p < 0.001). Western blotting analysis revealed the downregulation of AHR, FXR and PPARα protein expression in male KRAS mice at 26 weeks. For sex-related differences, CYP2E1 was male-predominant and UGT1A2 was female-predominant in the kidney. UGT1A1 and UGT1A5 expression was female-predominant, whereas UGT2B1 exhibited male-predominant expression in liver tissue. For the tissue distribution of DMEs, 21 subtypes of DMEs were all expressed in liver tissue. In the intestine, the expression levels of CYP2C29, CYP27A1, UGT1A2, 1A5, 1A6a, 1A9, 2B1, 2B5 and 2B36 were under the limitation of quantification. The subtypes of CYP7A1, 1B1, 2E1 and UGT1A1, 2A3, 2B34 were detected in kidney tissue.

Conclusions

This study, for the first time, unveils the variations and sex- and age-related differences in DMEs in C57 BL/6 (WT) mice and KRAS mice.

Expression and functional analysis of cytochrome P450 genes in the wolf spider Pardosa pseudoannulata under cadmium stress

Cytochrome P450 enzymes (CYPs), encoded by Halloween genes, mediate the biosynthesis of molting hormone, ecdysteroids, in arthropods. In this report, the effect of heavy metal cadmium (Cd) stress on the expression of cytochrome P450 genes in the wolf spider Pardosa pseudoannulata was analyzed. The results showed the expression levels of genes encoding for Cd transporters including ABC transporters, zinc transporters, calcium channel proteins and calcium binding proteins were inhibited or induced by Cd stress. In addition, the increase in metallothionein (MT) content and glutathione peroxidase (GPX) activity and decrease in total acetylcholine esterase (AChE) activity were also detected. Apparently, these detoxification methods did not completely protect the spider from the cytotoxicity of Cd stress. Increased mortality of P. pseudoannulata was observed when they were under Cd tress. In total 569 CYP genes belonging to 62 CYP subfamilies were obtained from P. pseudoannulata RNA-seq databases. BlaxtX analysis showed that 150, 161, 11, and 40 CYP genes were similar to the genes dib, phm, sad and shd, respectively, which are thought to catalyze the biosynthesis of ecdysteroids. Gene expression analysis suggested that 25 dib encoding genes, 27 phm encoding genes, 2 sad encoding genes, and 6 shd encoding genes were differentially expressed in TS2 vs. S2 comparison (Cd-treated 2nd instar spider vs. 2nd instar spider), respectively. There were 70 dib, 70 phm and 19 shd encoding genes either upregulated or downregulated, while 3 sad encoding genes were upregulated in TS5 vs. S5 (Cd-treated 5nd instar spider vs. 5nd instar spider). Genes related to heme binding and essential for activating the CYPs were also differentially expressed. Expression levels of cuticle related genes were significant differentially expressed, implying the changes in activities of chitin synthases and chitinase. Therefore we assume that unsuccessful molting process may occur on P. pseudoannulata due to influenced ecdysteroids levels, thus increasing mortality of spider.

On identifying collective displacements in apo-proteins that reveal eventual binding pathways

Binding of small molecules to proteins often involves large conformational changes in the latter, which open up pathways to the binding site. Observing and pinpointing these rare events in large scale, all-atom, computations of specific protein-ligand complexes, is expensive and to a great extent serendipitous. Further, relevant collective variables which characterise specific binding or un-binding scenarios are still difficult to identify despite the large body of work on the subject. Here, we show that possible primary and secondary binding pathways can be discovered from short simulations of the apo-protein without waiting for an actual binding event to occur. We use a projection formalism, introduced earlier to study deformation in solids, to analyse local atomic displacements into two mutually orthogonal subspaces—those which are “affine” i.e. expressible as a homogeneous deformation of the native structure, and those which are not. The susceptibility to non-affine displacements among the various residues in the apo- protein is then shown to correlate with typical binding pathways and sites crucial for allosteric modifications. We validate our observation with all-atom computations of three proteins, T4-Lysozyme, Src kinase and Cytochrome P450.

Designing drugs which target specific proteins involved in diseases consumes a lot of time and effort in the pharmaceutical industry. In recent times, in silico design of drugs using all-atom molecular modelling has started to provide crucial inputs. Even so, discovery of binding pathways of small molecules both at the primary binding site, as well as sites for allosteric control, is time consuming and often fortuitous. We provide here a framework within which critical conformational changes likely to occur during binding are quantified from statistical analysis of configurations of proteins in their apo, or inactive form, greatly simplifying identification of target residues. We illustrate this idea by analysing ligand binding pathways for three proteins T4- Lysozyme, P450 and Src kinase, which are active respectively in the immune system, metabolism and cancer.

Intraindividual Variation and Correlation of Cytochrome P450 Activities in Human Liver Microsomes

We systematically studied and explored intraindividual variation and correlation in the activities of cytochrome P450 (CYP) using 105 normal liver samples. The intraindividual percentage coefficients of variation of relative K_m, V_max, and CL_int were 45.9% (18.3-140%), 44.0% (16.8-127%), and 52.0% (24.2-134%). Overall values of relative K_m, V_max, and CL_int for 10 CYPs were sorted. The order, from highest to lowest, was CYP2D6, 3A4/5, 2C19, 2C9, 2B6, 1A2, 2A6, 2C8, and 2E1 for K_m; CYP3A4/5, 2C19, 2D6, 2C8, 2E1, 2B6, 1A2, 2A6, and 2C9 for V_max; and CYP2D6, 2B6, 3A4/5, 2C19, 2C9, 1A2, 2E1, 2C8, and 2A6 for CL_int. CYP2A6*4, 2B6 785A>G, and 2D6 100C>T contributed to intraindividual variation of CYP activities. The correlations and similarities among 10 CYP activities were analyzed, and it was found that the highest correlation and similarity for V_max, K_m, and CL_int occurred between CYP2C9 and 2C8, CYP2C9 and 1A2, and CYP2C9 and 2C8, respectively. In conclusion, CYP activities exhibit considerable intraindividual variation, with some notable correlations, which might be helpful to comprehensively understand the internal connection among CYP activities and to guide the rational use of drugs.

Toward a systems approach to the human cytochrome P450 ensemble: interactions between CYP2D6 and CYP2E1 and their functional consequences

Functional cross-talk among human drug-metabolizing cytochrome P450 through their association is a topic of emerging importance. Here, we studied the interactions of human CYP2D6, a major metabolizer of psychoactive drugs, with one of the most prevalent human P450 enzymes, ethanol-inducible CYP2E1. Detection of P450-P450 interactions was accomplished through luminescence resonance energy transfer between labeled proteins incorporated into human liver microsomes and the microsomes of insect cells containing NADPH-cytochrome P450 reductase. The potential of CYP2D6 to form oligomers in the microsomal membrane is among the highest observed with human cytochrome P450 studied up to date. We also observed the formation of heteromeric complexes of CYP2D6 with CYP2E1 and CYP3A4, and found a significant modulation of these interactions by 3,4-methylenedioxymethylamphetamine, a widespread drug of abuse metabolized by CYP2D6. Our results demonstrate an ample alteration of the catalytic properties of CYP2D6 and CYP2E1 caused by their association. In particular, we demonstrated that preincubation of microsomes containing co-incorporated CYP2D6 and CYP2E1 with CYP2D6-specific substrates resulted in considerable time-dependent activation of CYP2D6, which presumably occurs via a slow substrate-induced reorganization of CYP2E1-CYP2D6 hetero-oligomers. Furthermore, we demonstrated that the formation of heteromeric complexes between CYP2E1 and CYP2D6 affects the stoichiometry of futile cycling and substrate oxidation by CYP2D6 by means of decreasing the electron leakage through the peroxide-generating pathways. Our results further emphasize the role of P450-P450 interactions in regulatory cross-talk in human drug-metabolizing ensemble and suggest a role of interactions of CYP2E1 with CYP2D6 in pharmacologically important instances of alcohol-drug interactions.

Immobilized Cytochrome P450 for Monitoring of P450-P450 Interactions and Metabolism

Cytochrome P450 (P450) protein-protein interactions have been shown to alter their catalytic activity. Furthermore, these interactions are isoform specific and can elicit activation, inhibition, or no effect on enzymatic activity. Studies show that these effects are also dependent on the protein partner cytochrome P450 reductase (CPR) and the order of protein addition to purified reconstituted enzyme systems. In this study, we use controlled immobilization of P450s to a gold surface to gain a better understanding of P450-P450 interactions between three key drug-metabolizing isoforms (CYP2C9, CYP3A4, and CYP2D6). Molecular modeling was used to assess the favorability of homomeric/heteromeric P450 complex formation. P450 complex formation in vitro was analyzed in real time utilizing surface plasmon resonance. Finally, the effects of P450 complex formation were investigated utilizing our immobilized platform and reconstituted enzyme systems. Molecular modeling shows favorable binding of CYP2C9-CPR, CYP2C9-CYP2D6, CYP2C9-CYP2C9, and CYP2C9-CYP3A4, in rank order.KDvalues obtained via surface plasmon resonance show strong binding, in the nanomolar range, for the above pairs, with CYP2C9-CYP2D6 yielding the lowestKD, followed by CYP2C9-CYP2C9, CYP2C9-CPR, and CYP2C9-CYP3A4. Metabolic incubations show that immobilized CYP2C9 metabolism was activated by homomeric complex formation. CYP2C9 metabolism was not affected by the presence of CYP3A4 with saturating CPR concentrations. CYP2C9 metabolism was activated by CYP2D6 at saturating CPR concentrations in solution but was inhibited when CYP2C9 was immobilized. The order of addition of proteins (CYP2C9, CYP2D6, CYP3A4, and CPR) influenced the magnitude of inhibition for CYP3A4 and CYP2D6. These results indicate isoform-specific P450 interactions and effects on P450-mediated metabolism.

The combination of a novel 2 bp deletion mutation and p.D63H in CYP11B1 cause congenital adrenal hyperplasia due to steroid 11β-hydroxylase deficiency

Deficiency of steroid 11β-hydroxylase activity occurs in 5-8% of patients with congenital adrenal hyperplasia (CAH). The aim of the current study was to identify mutations in the CYP11B1 gene of a patient with CAH due to deficiency of steroid 11β-hydroxylase activity, and to study the functional and structural consequences of these mutations. A molecular genetic analysis of the CYP11B1 gene in this patient and her parents identified a known missense mutation g.5194G>C (p.D63H) and a novel 2 bp deletion mutation (g.9525_9526delCT, corresponding to p.L380V…R420X) in the patient. In vitro expression studies in COS7 cells revealed a decreased 11β-hydroxylase activity in the p.D63H mutant to 2.0±0.8% and in the p.L380V…R420X mutant to 0.2±2.2% for the conversion of 11-deoxycortisol to cortisol. Three dimensional homology models for the normal and mutant proteins were built by using the recently published x-ray structure of the human CYP11B2 as a template. Presumably, the g.9525_9526delCT mutation in CYP11B1 resulted in a truncated protein with a misfolded C-terminal domain that could not efficiently bind heme iron, substrate, and adrenodoxin and had lost its biochemical function. In summary, CAH due to steroid 11β-hydroxylase deficiency can be attributed to both the novel deletion mutation (g.9525_9526delCT, corresponding to p.L380V…R420X) and known missense mutation (g.5194G>C corresponding to p.D63H) in CYP11B1.

Severely Impaired and Dysregulated Cytochrome P450 Expression and Activities in Hepatocellular Carcinoma: Implications for Personalized Treatment in Patients

This study aims to systematically determine the activities and expressions of cytochrome P450s (CYP) in hepatocellular carcinoma (HCC) patients to support their optimal use in personalized treatment of HCC. Activities of seven major drug-metabolizing CYP enzymes (CYP1A2, 2A6, 2C8, 2C9, 2D6, 2E1, and 3A4) were determined in tumors and pericarcinomatous tissues harvested from 26 patients with hepatitis B virus-positive HCC using probe substrates. Protein and mRNA levels of these CYPs were also measured using isotope label-free LC/MS-MS method and real-time PCR, respectively. Maximal metabolic velocity (Vmax) of CYP probe substrates was decreased by 2.5- to 30-fold in tumor microsomes, accompanied by a corresponding decrease in their protein and mRNA expression levels. However, Km values and turnover numbers of substrates in tumor microsomes were not changed. High correlations between activities and CYP protein levels were also observed, but the correlation between activities and mRNA levels was often poor. There was a major decrease in the degree of correlation in CYP expression in tumor tissues, suggesting that CYP expression levels are greatly disrupted by the tumorigenic process. Our unprecedented systemic study of the effects of HCC on CYPs demonstrated that activities of CYPs were seriously impaired and their expression patterns were severely altered by HCC. We proposed that determination of the CYP protein expression profile by LC/MS-MS in each patient is a promising approach that can be clinically used for individualized treatment of HCC.

Modulation of Metabolic Detoxification Pathways Using Foods and Food-Derived Components: A Scientific Review with Clinical Application

Research into human biotransformation and elimination systems continues to evolve. Various clinical and in vivo studies have been undertaken to evaluate the effects of foods and food-derived components on the activity of detoxification pathways, including phase I cytochrome P450 enzymes, phase II conjugation enzymes, Nrf2 signaling, and metallothionein. This review summarizes the research in this area to date, highlighting the potential for foods and nutrients to support and/or modulate detoxification functions. Clinical applications to alter detoxification pathway activity and improve patient outcomes are considered, drawing on the growing understanding of the relationship between detoxification functions and different disease states, genetic polymorphisms, and drug-nutrient interactions. Some caution is recommended, however, due to the limitations of current research as well as indications that many nutrients exert biphasic, dose-dependent effects and that genetic polymorphisms may alter outcomes. A whole-foods approach may, therefore, be prudent.

Allosteric communication pathways and thermal rectification in PDZ-2 protein: a computational study

Allosteric communication in proteins is a fundamental and yet unresolved problem of structural biochemistry. Previous findings, from computational biology ( Ota, N.; Agard, D. A. J. Mol. Biol. 2005 , 351 , 345 - 354 ), have proposed that heat diffuses in a protein through cognate protein allosteric pathways. This work studied heat diffusion in the well-known PDZ-2 protein, and confirmed that this protein has two cognate allosteric pathways and that heat flows preferentially through these. Also, a new property was also observed for protein structures: heat diffuses asymmetrically through the structures. The underling structure of this asymmetrical heat flow was a normal length hydrogen bond (∼2.85 Å) that acted as a thermal rectifier. In contrast, thermal rectification was compromised in short hydrogen bonds (∼2.60 Å), giving rise to symmetrical thermal diffusion. Asymmetrical heat diffusion was due, on a higher scale, to the local, structural organization of residues that, in turn, was also mediated by hydrogen bonds. This asymmetrical/symmetrical energy flow may be relevant for allosteric signal communication directionality in proteins and for the control of heat flow in materials science.

Significantly decreased and more variable expression of major CYPs and UGTs in liver microsomes prepared from HBV-positive human hepatocellular carcinoma and matched pericarcinomatous tissues determined using an isotope label-free UPLC-MS/MS method

PURPOSE

To determine the liver expression of cytochrome P450 (CYPs) and uridine 5'-diphosphate-glucuronosyltransferases (UGTs), the major phase I and II metabolism enzymes responsible for clearance and detoxification of drugs, xenobiotic and endogenous substances.

METHODS

A validated isotope label-free method was established for absolute and simultaneous quantification of 9 CYPs (1A2, 2A6, 2B6, 2C8, 2C9, 2C19, 2D, 2E1 and 3A4) and 5 UGTs (1A1, 1A4, 1A6, 1A9 and 2B7) in human liver microsomes using LC-MS/MS.

RESULTS

The LC-MS/MS method displayed excellent dynamic range (at least 250-fold) and high sensitivity for each of the signature peptides with acceptable recovery, accuracy and precision. The protein expression profile of CYP and UGT isoforms were then determined in match microsomes samples prepared from patients with HBV-positive human hepatocellular carcinoma (HCC). In the tumor microsomes, the average absolute amounts of 8 major CYP isoforms (except CYP2C19) and 3 UGT isoforms (UGT1A1, UGT1A4 and UGT2B7) were decreased significantly (p < 0.05), whereas UGT1A6 and UGT1A9 levels were unchanged (p > 0.05). In addition, among isoforms with altered expression, 6 of 8 CYP isoforms and all three UGT isoforms were much more variable in tumor microsomes. Lastly, the importance of CYP3A4 was greatly diminished whereas the importance of UGT1A6 was enhanced in tumor microsomes.

CONCLUSION

The use of an isotope label-free absolute quantification method for the simultaneous determination of 9 CYPs and 5 UGTs in human liver microsomes reveals that expression levels of CYPs and UGTs in human liver are severely impact by HCC, which could impact drug metabolism, disposition and pharmacotherapy.

Interactions among cytochromes P450 in microsomal membranes: oligomerization of cytochromes P450 3A4, 3A5, and 2E1 and its functional consequences

The body of evidence of physiologically relevant P450-P450 interactions in microsomal membranes continues to grow. Here we probe oligomerization of human CYP3A4, CYP3A5, and CYP2E1 in microsomal membranes. Using a technique based on luminescence resonance energy transfer, we demonstrate that all three proteins are subject to a concentration-dependent equilibrium between the monomeric and oligomeric states. We also observed the formation of mixed oligomers in CYP3A4/CYP3A5, CYP3A4/CYP2E1, and CYP3A5/CYP2E1 pairs and demonstrated that the association of either CYP3A4 or CYP3A5 with CYP2E1 causes activation of the latter enzyme. Earlier we hypothesized that the intersubunit interface in CYP3A4 oligomers is similar to that observed in the crystallographic dimers of some microsomal drug-metabolizing cytochromes P450 (Davydov, D. R., Davydova, N. Y., Sineva, E. V., Kufareva, I., and Halpert, J. R. (2013) Pivotal role of P450-P450 interactions in CYP3A4 allostery: the case of α-naphthoflavone. Biochem. J. 453, 219-230). Here we report the results of intermolecular cross-linking of CYP3A4 oligomers with thiol-reactive bifunctional reagents as well as the luminescence resonance energy transfer measurements of interprobe distances in the oligomers of labeled CYP3A4 single-cysteine mutants. The results provide compelling support for the physiological relevance of the dimer-specific peripheral ligand-binding site observed in certain CYP3A4 structures. According to our interpretation, these results reveal an important general mechanism that regulates the activity and substrate specificity of the cytochrome P450 ensemble through interactions between multiple P450 species. As a result of P450-P450 cross-talk, the catalytic properties of the cytochrome P450 ensemble cannot be predicted by simple summation of the properties of the individual P450 species.

Identifying key residues for protein allostery through rigid residue scan

Allostery is a ubiquitous process for protein regulatory activity in which a binding event can change a protein's function carried out at a distal site. Despite intensive theoretical and experimental investigation of protein allostery in the past five decades, effective methods have yet to be developed that can systematically identify key residues involved in allosteric mechanisms. In this study, we propose the rigid residue scan as a systematic approach to identify important allosteric residues. The third PDZ domain (PDZ3) in the postsynaptic density 95 protein (PSD-95) is used as a model system, and each amino acid residue is treated as a single rigid body during independent molecular dynamics simulations. Various indices based on cross-correlation matrices are used, which allow for two groups of residues with different functions to be identified. The first group is proposed as "switches" that are needed to "turn on" the binding effect of protein allostery. The second group is proposed as "wire residues" that are needed to propagate energy or information from the binding site to distal locations within the same protein. Among the nine residues suggested as important for PDZ3 intramolecular communication in this study, eight have been reported as critical for allostery in PDZ3. Therefore, the rigid residue scan approach is demonstrated to be an effective method for systemically identifying key residues in protein intramolecular communication and allosteric mechanisms.

Gene expression profile of Bombyx mori hemocyte under the stress of destruxin A

Destruxin A (DA) is a cyclo-peptidic mycotoxin from the entomopathogenic fungus Metarhizium anisopliae. To uncover potential genes associated with its molecular mechanisms, a digital gene expression (DGE) profiling analysis was used to compare differentially expressed genes in the hemocytes of silkworm larvae treated with DA. Ten DGE libraries were constructed, sequenced, and assembled, and the unigenes with least 2.0-fold difference were further analyzed. The numbers of up-regulated genes were 10, 20, 18, 74 and 8, as well as the numbers of down-regulated genes were 0, 1, 8, 13 and 3 at 1, 4, 8, 12 and 24 h post treatment, respectively. Totally, the expression of 132 genes were significantly changed, among them, 1, 3 and 12 genes were continually up-regulated at 4, 3 and 2 different time points, respectively, while 1 gene was either up or down-regulated continually at 2 different time points. Furthermore, 68 genes were assigned to one or multiple gene ontology (GO) terms and 89 genes were assigned to specific Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology. In-depth analysis identified that these genes putatively involved in insecticide resistance, cell apoptosis, and innate immune defense. Finally, twenty differentially expressed genes were randomly chosen and validated by quantitative real-time PCR (qRT-PCR). Our studies provide insights into the toxic effect of this microbial insecticide on silkworm's hemocytes, and are helpful to better understanding of the molecular mechanisms of DA as a biological insecticide.