What common structural features and variations of mammalian P450s are known to date?

Identification of the N-terminal residues responsible for the differential microdomain localization of CYP1A1 and CYP1A2

The endoplasmic reticulum is organized into ordered regions enriched in cholesterol and sphingomyelin, and disordered microdomains characterized by more fluidity. Rabbit CYP1A1 and CYP1A2 localize into disordered and ordered microdomains, respectively. Previously, a CYP1A2 chimera containing the first 109 amino acids of CYP1A1 showed altered microdomain localization. The goal of this study was to identify specific residues responsible for CYP1A microdomain localization. Thus, CYP1A2 chimeras containing substitutions from homologous regions of CYP1A1 were expressed in HEK 293T/17 cells, and the localization was examined after solubilization with Brij 98. A CYP1A2 mutant with the three amino acids from CYP1A1 (VAG) at positions 27 to 29 of CYP1A2 was generated that showed a distribution pattern similar to those of CYP1A1/1A2 chimeras containing both the first 109 amino acids and the first 31 amino acids of CYP1A1 followed by remaining amino acids of CYP1A2. Similarly, the reciprocal substitution of three amino acids from CYP1A2 (AVR) into CYP1A1 resulted in a partial redistribution of the chimera into ordered microdomains. Molecular dynamic simulations indicate that the positive charges of the CYP1A1 and CYP1A2 linker regions between the N termini and catalytic domains resulted in different depths of immersion of the N termini in the membrane. The overlap of the distribution of positively charged residues in CYP1A2 (AVR) and negatively charged phospholipids was higher in the ordered than in the disordered microdomain. These findings identify three residues in the CYP1AN terminus as a novel microdomain-targeting motif of the P450s and provide a mechanistic explanation for the differential microdomain localization of CYP1A.

The Multienzyme Complex Nature of Dehydroepiandrosterone Sulfate Biosynthesis

Dehydroepiandrosterone (DHEA), a precursor of steroid sex hormones, is synthesized by steroid 17-alpha-hydroxylase/17,20-lyase (CYP17A1) with the participation of microsomal cytochrome b5 (CYB5A) and cytochrome P450 reductase (CPR), followed by sulfation by two cytosolic sulfotransferases, SULT1E1 and SULT2A1, for storage and transport to tissues in which its synthesis is not available. The involvement of CYP17A1 and SULTs in these successive reactions led us to consider the possible interaction of SULTs with DHEA-producing CYP17A1 and its redox partners. Text mining analysis, protein-protein network analysis, and gene co-expression analysis were performed to determine the relationships between SULTs and microsomal CYP isoforms. For the first time, using surface plasmon resonance, we detected interactions between CYP17A1 and SULT2A1 or SULT1E1. SULTs also interacted with CYB5A and CPR. The interaction parameters of SULT2A1/CYP17A1 and SULT2A1/CYB5A complexes seemed to be modulated by 3'-phosphoadenosine-5'-phosphosulfate (PAPS). Affinity purification, combined with mass spectrometry (AP-MS), allowed us to identify a spectrum of SULT1E1 potential protein partners, including CYB5A. We showed that the enzymatic activity of SULTs increased in the presence of only CYP17A1 or CYP17A1 and CYB5A mixture. The structures of CYP17A1/SULT1E1 and CYB5A/SULT1E1 complexes were predicted. Our data provide novel fundamental information about the organization of microsomal CYP-dependent macromolecular complexes.

Structural Characterization and Molecular Dynamics Study of the REPI Fusion Protein from Papaver somniferum L

REPI is a pivotal point enzyme in plant benzylisoquinoline alkaloid metabolism as it promotes the evolution of the biosynthetic branch of morphinan alkaloids. Experimental studies of its activity led to the identification of two modules (DRS and DRR) that catalyze two sequential steps of the epimerization of (S)- to (R)-reticuline. Recently, special attention has been paid to its genetic characterization and evolutionary history, but no structural analyses of the REPI protein have been conducted to date. We present here a computational structural characterization of REPI with heme and NADP cofactors in the apo state and in three complexes with substrate (S)-reticuline in DRS and intermediate 1,2-dehydroreticuline in DRS and in DRR. Since no experimental structure exists for REPI, we used its AlphaFold model as a scaffold to build up these four systems, which were submitted to all-atom molecular dynamics (MD) simulations. A comparison of MD results for the four systems revealed key dynamic changes associated with cofactor and ligand binding and provided a dynamic picture of the evolution of their structures and interactions. We also explored the possible dynamic occurrence of tunnels and electrostatic highways potentially involved in alternative mechanisms for channeling the intermediate from DRS to DRR.

Pharmacokinetics during therapeutic hypothermia in neonates: from pathophysiology to translational knowledge and physiologically-based pharmacokinetic (PBPK) modeling

INTRODUCTION

Perinatal asphyxia (PA) still causes significant morbidity and mortality. Therapeutic hypothermia (TH) is the only effective therapy for neonates with moderate to severe hypoxic-ischemic encephalopathy after PA. These neonates need additional pharmacotherapy, and both PA and TH may impact physiology and, consequently, pharmacokinetics (PK) and pharmacodynamics (PD).

AREAS COVERED

This review provides an overview of the available knowledge in PubMed (until November 2022) on the pathophysiology of neonates with PA/TH. In vivo pig models for this setting enable distinguishing the effect of PA versus TH on PK and translating this effect to human neonates. Available asphyxia pig models and methodological considerations are described. A summary of human neonatal PK of supportive pharmacotherapy to improve neurodevelopmental outcomes is provided.

EXPERT OPINION

To support drug development for this population, knowledge from clinical observations (PK data, real-world data on physiology), preclinical (in vitro and in vivo (minipig)) data, and molecular and cellular biology insights can be integrated into a predictive physiologically-based PK (PBPK) framework, as illustrated by the I-PREDICT project (Innovative physiology-based pharmacokinetic model to predict drug exposure in neonates undergoing cooling therapy). Current knowledge, challenges, and expert opinion on the future directions of this research topic are provided.

Characterizing the combined effects of cytochrome P450 missense variation within star allele definitions

Background: Cytochrome P450 (CYP) genetic variation largely impacts drug response. However, many CYP star alleles (haplotypes) lack functional annotation, impeding our understanding of drug metabolism mechanisms. We aimed to investigate the impact of missense variant combinations on CYP protein structures. Methods: Normal mode analysis was conducted on 261 missense variants within 91 CYP haplotypes. CYP2D6*2 and CYP2D6*17 were prioritized for molecular dynamics simulation. Results: Normal mode analysis and molecular dynamics highlight the effects of known CYP missense variants on protein stability and conformational dynamics. Missense variants within haplotypes may have intermodulating effects on protein structure and function. Conclusion: This study highlights the utility of multiscale modeling in interpreting CYP missense variants and particularly their combinations within various star alleles.

Prediction of Sites of Metabolism of CYP3A4 Substrates Utilizing Docking-Derived Geometric Features

Cytochrome P450 3A4 (CYP3A4) is one of the major drug-metabolizing enzymes in the human body and is responsible for the metabolism of ∼50% of clinically used drugs. Therefore, the identification of the compound's sites of metabolism (SOMs) mediated by CYP3A4 is of utmost importance in the early stage of drug discovery and development. Herein, docking-based approaches incorporating geometric features were used for SOMs prediction of CYP3A4 substrates. The cross-docking poses of a relatively large data set containing 474 substrates were analyzed in depth, and a widely observed geometric pattern called the close proximity of SOMs was derived from the poses. On the basis of the close proximity, several structure-based models have been constructed, which demonstrated better performance than those structure-based models using the criterion of Fe-SOM distance. For further improving the prediction performance, the structure-based models were also combined with the well-known ligand-based model SMARTCyp. One combined model exhibited good performance on the SOMs prediction of an external substrate set containing kinase inhibitors, PROTACs, approved drugs, and some lead compounds.

Heterologous Expression of Recombinant Human Cytochrome P450 (CYP) in Escherichia coli: N-Terminal Modification, Expression, Isolation, Purification, and Reconstitution

Cytochrome P450 (CYP) enzymes play important roles in metabolising endogenous and xenobiotic substances. Characterisations of human CYP proteins have been advanced with the rapid development of molecular technology that allows heterologous expression of human CYPs. Among several hosts, bacteria systems such as Escherichia coli (E. coli) have been widely used thanks to their ease of use, high level of protein yields, and affordable maintenance costs. However, the levels of expression in E. coli reported in the literature sometimes differ significantly. This paper aims to review several contributing factors, including N-terminal modifications, co-expression with a chaperon, selections of vectors and E. coli strains, bacteria culture and protein expression conditions, bacteria membrane preparations, CYP protein solubilizations, CYP protein purifications, and reconstitution of CYP catalytic systems. The common factors that would most likely lead to high expression of CYPs were identified and summarised. Nevertheless, each factor may still require careful evaluation for individual CYP isoforms to achieve a maximal expression level and catalytic activity. Recombinant E. coli systems have been evidenced as a useful tool in obtaining the ideal level of human CYP proteins, which ultimately allows for subsequent characterisations of structures and functions.

Comprehensive in vitro and in silico assessments of metabolic capabilities of 24 genomic variants of CYP2C19 using two different substrates

Introduction: Most hepatically cleared drugs are metabolized by cytochromes P450 (CYPs), and Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines provide curated clinical references for CYPs to apply individual genome data for optimized drug therapy. However, incorporating novel pharmacogenetic variants into guidelines takes considerable time. Methods: We comprehensively assessed the drug metabolizing capabilities of CYP2C19 variants discovered through population sequencing of two substrates, S-mephenytoin and omeprazole. Results: Based on established functional assays, 75% (18/24) of the variants not yet described in Pharmacogene Variation (PharmVar) had significantly altered drug metabolizing capabilities. Of them, seven variants with inappreciable protein expression were evaluated as protein damaging by all three in silico prediction algorithms, Sorting intolerant from tolerant (SIFT), Polymorphism Phenotyping v2 (PolyPhen-2), and Combined annotation dependent depletion (CADD). The five variants with decreased metabolic capability (<50%) of wild type for either substrates were evaluated as protein damaging by all three in silico prediction algorithms, except CADD exact score of NM_000769.4:c.593T>C that was 19.68 (<20.0). In the crystal structure of the five polymorphic proteins, each altered residue of all those proteins was observed to affect the key structures of drug binding specificity. We also identified polymorphic proteins indicating different tendencies of metabolic capability between the two substrates (5/24). Discussion: Therefore, we propose a methodology that combines in silico prediction algorithms and functional assays on polymorphic CYPs with multiple substrates to evaluate the changes in the metabolism of all possible genomic variants in CYP genes. The approach would reinforce existing guidelines and provide information for prescribing appropriate medicines for individual patients.

Catch-up Growth and Discontinuation of Fludrocortisone Treatment in Aldosterone Synthase Deficiency

BACKGROUND

Aldosterone synthase deficiency (ASD) caused by mutations in the CYP11B2 gene is characterized by isolated mineralocorticoid deficiency. Data are scarce regarding clinical and biochemical outcomes of the disease in the follow-up.

OBJECTIVE

Assessment of the growth and steroid profiles of patients with ASD at the time of diagnosis and after discontinuation of treatment.

DESIGN AND METHOD

Children with clinical diagnosis of ASD were included in a multicenter study. Growth and treatment characteristics were recorded. Plasma adrenal steroids were measured using liquid chromatography-mass spectrometry. Genetic diagnosis was confirmed by CYP11B2 gene sequencing and in silico analyses.

RESULTS

Sixteen patients from 12 families were included (8 females; median age at presentation: 3.1 months, range: 0.4 to 8.1). The most common symptom was poor weight gain (56.3%). Median age of onset of fludrocortisone treatment was 3.6 months (range: 0.9 to 8.3). Catch-up growth was achieved at median 2 months (range: 0.5 to 14.5) after treatment. Fludrocortisone could be stopped in 5 patients at a median age of 6.0 years (range: 2.2 to 7.6). Plasma steroid profiles revealed reduced aldosterone synthase activity both at diagnosis and after discontinuation of treatment compared to age-matched controls. We identified 6 novel (p.Y195H, c.1200 + 1G > A, p.F130L, p.E198del, c.1122-18G > A, p.I339_E343del) and 4 previously described CYP11B2 variants. The most common variant (40%) was p.T185I.

CONCLUSIONS

Fludrocortisone treatment is associated with a rapid catch-up growth and control of electrolyte imbalances in ASD. Decreased mineralocorticoid requirement over time can be explained by the development of physiological adaptation mechanisms rather than improved aldosterone synthase activity. As complete biochemical remission cannot be achieved, a long-term surveillance of these patients is required.

Uncovering of cytochrome P450 anatomy by SecStrAnnotator

Protein structural families are groups of homologous proteins defined by the organization of secondary structure elements (SSEs). Nowadays, many families contain vast numbers of structures, and the SSEs can help to orient within them. Communities around specific protein families have even developed specialized SSE annotations, always assigning the same name to the equivalent SSEs in homologous proteins. A detailed analysis of the groups of equivalent SSEs provides an overview of the studied family and enriches the analysis of any particular protein at hand. We developed a workflow for the analysis of the secondary structure anatomy of a protein family. We applied this analysis to the model family of cytochromes P450 (CYPs)-a family of important biotransformation enzymes with a community-wide used SSE annotation. We report the occurrence, typical length and amino acid sequence for the equivalent SSE groups, the conservation/variability of these properties and relationship to the substrate recognition sites. We also suggest a generic residue numbering scheme for the CYP family. Comparing the bacterial and eukaryotic part of the family highlights the significant differences and reveals a well-known anomalous group of bacterial CYPs with some typically eukaryotic features. Our workflow for SSE annotation for CYP and other families can be freely used at address https://sestra.ncbr.muni.cz .

Cytochrome P450 expression and regulation in the brain

The regulation of brain cytochrome P450 enzymes (CYPs) is different compared with respective hepatic enzymes. This may result from anatomical bases and physiological functions of the two organs. The brain is composed of a variety of functional structures built of different interconnected cell types endowed with specific receptors that receive various neuronal signals from other brain regions. Those signals activate transcription factors or alter functioning of enzyme proteins. Moreover, the blood-brain barrier (BBB) does not allow free penetration of all substances from the periphery into the brain. Differences in neurotransmitter signaling, availability to endogenous and exogenous active substances, and levels of transcription factors between neuronal and hepatic cells lead to differentiated expression and susceptibility to the regulation of CYP genes in the brain and liver. Herein, we briefly describe the CYP enzymes of CYP1-3 families, their distribution in the brain, and discuss brain-specific regulation of CYP genes. In parallel, a comparison to liver CYP regulation is presented. CYP enzymes play an essential role in maintaining the levels of bioactive molecules within normal ranges. These enzymes modulate the metabolism of endogenous neurochemicals, such as neurosteroids, dopamine, serotonin, melatonin, anandamide, and exogenous substances, including psychotropics, drugs of abuse, neurotoxins, and carcinogens. The role of these enzymes is not restricted to xenobiotic-induced neurotoxicity, but they are also involved in brain physiology. Therefore, it is crucial to recognize the function and regulation of CYP enzymes in the brain to build a foundation for future medicine and neuroprotection and for personalized treatment of brain diseases.

Domain-Swap Dimerization of Acanthamoeba castellanii CYP51 and a Unique Mechanism of Inactivation by Isavuconazole

Cytochromes P450 (P450, CYP) metabolize a wide variety of endogenous and exogenous lipophilic molecules, including most drugs. Sterol 14α-demethylase (CYP51) is a target for antifungal drugs known as conazoles. Using X-ray crystallography, we have discovered a domain-swap homodimerization mode in CYP51 from a human pathogen, Acanthamoeba castellanii CYP51 (AcCYP51). Recombinant AcCYP51 with a truncated transmembrane helix was purified as a heterogeneous mixture corresponding to the dimer and monomer units. Spectral analyses of these two populations have shown that the CO-bound ferrous form of the dimeric protein absorbed at 448 nm (catalytically competent form), whereas the monomeric form absorbed at 420 nm (catalytically incompetent form). AcCYP51 dimerized head-to-head via N-termini swapping, resulting in formation of a nonplanar protein-protein interface exceeding 2000 Å2 with a total solvation energy gain of -35.4 kcal/mol. In the dimer, the protomers faced each other through the F and G α-helices, thus blocking the substrate access channel. In the presence of the drugs clotrimazole and isavuconazole, the AcCYP51 drug complexes crystallized as monomers. Although clotrimazole-bound AcCYP51 adopted a typical CYP monomer structure, isavuconazole-bound AcCYP51 failed to refold 74 N-terminal residues. The failure of AcCYP51 to fully refold upon inhibitor binding in vivo would cause an irreversible loss of a structurally aberrant enzyme through proteolytic degradation. This assumption explains the superior potency of isavuconazole against A. castellanii The dimerization mode observed in this work is compatible with membrane association and may be relevant to other members of the CYP family of biologic, medical, and pharmacological importance. SIGNIFICANCE STATEMENT: We investigated the mechanism of action of antifungal drugs in the human pathogen Acanthamoeba castellanii. We discovered that the enzyme target [Acanthamoeba castellanii sterol 14α-demethylase (AcCYP51)] formed a dimer via an N-termini swap, whereas drug-bound AcCYP51 was monomeric. In the AcCYP51-isavuconazole complex, the protein target failed to refold 74 N-terminal residues, suggesting a fundamentally different mechanism of AcCYP51 inactivation than only blocking the active site. Proteolytic degradation of a structurally aberrant enzyme would explain the superior potency of isavuconazole against A. castellanii.

CYP2J2 Molecular Recognition: A New Axis for Therapeutic Design

Cytochrome P450 (CYP) epoxygenases are a special subset of heme-containing CYP enzymes capable of performing the epoxidation of polyunsaturated fatty acids (PUFA) and the metabolism of xenobiotics. This dual functionality positions epoxygenases along a metabolic crossroad. Therefore, structure-function studies are critical for understanding their role in bioactive oxy-lipid synthesis, drug-PUFA interactions, and for designing therapeutics that directly target the epoxygenases. To better exploit CYP epoxygenases as therapeutic targets, there is a need for improved understanding of epoxygenase structure-function. Of the characterized epoxygenases, human CYP2J2 stands out as a potential target because of its role in cardiovascular physiology. In this review, the early research on the discovery and activity of epoxygenases is contextualized to more recent advances in CYP epoxygenase enzymology with respect to PUFA and drug metabolism. Additionally, this review employs CYP2J2 epoxygenase as a model system to highlight both the seminal works and recent advances in epoxygenase enzymology. Herein we cover CYP2J2's interactions with PUFAs and xenobiotics, its tissue-specific physiological roles in diseased states, and its structural features that enable epoxygenase function. Additionally, the enumeration of research on CYP2J2 identifies the future needs for the molecular characterization of CYP2J2 to enable a new axis of therapeutic design.

A Review of the Important Role of CYP2D6 in Pharmacogenomics

Cytochrome P450 2D6 (CYP2D6) is a critical pharmacogene involved in the metabolism of ~20% of commonly used drugs across a broad spectrum of medical disciplines including psychiatry, pain management, oncology and cardiology. Nevertheless, CYP2D6 is highly polymorphic with single-nucleotide polymorphisms, small insertions/deletions and larger structural variants including multiplications, deletions, tandem arrangements, and hybridisations with non-functional CYP2D7 pseudogenes. The frequency of these variants differs across populations, and they significantly influence the drug-metabolising enzymatic function of CYP2D6. Importantly, altered CYP2D6 function has been associated with both adverse drug reactions and reduced drug efficacy, and there is growing recognition of the clinical and economic burdens associated with suboptimal drug utilisation. To date, pharmacogenomic clinical guidelines for at least 48 CYP2D6-substrate drugs have been developed by prominent pharmacogenomics societies, which contain therapeutic recommendations based on CYP2D6-predicted categories of metaboliser phenotype. Novel algorithms to interpret CYP2D6 function from sequencing data that consider structural variants, and machine learning approaches to characterise the functional impact of novel variants, are being developed. However, CYP2D6 genotyping is yet to be implemented broadly into clinical practice, and so further effort and initiatives are required to overcome the implementation challenges and deliver the potential benefits to the bedside.

Influence of nanoparticles on liver tissue and hepatic functions: A review

Due to the increasing interest in nanotechnology in very large application fields, including biotechnology, electronics and food industries, humans are increasingly exposed to nanoparticles (NPs). Consequently, the question about the safety of these nanomaterials and their impact on human health is a legitimate concern. The liver is the primary organ of detoxification and is one of the tissues that is most exposed to NPs. When they reach the bloodstream, NPs are mainly internalized by liver cells. This review focuses on recent in vitro and in vivo studies addressing the effects of organic and inorganic NPs on the liver. Specifically, the impact of the NPs on hepatic enzyme activities, the inflammatory response and genotoxicity processes will be described. Depending on the physicochemical parameters of the NPs and the conditions of exposure, NPs could lead to global liver injury.

Investigation on drug-binding in heme pocket of CYP2C19 with UV-visible and resonance Raman spectroscopies

Cytochrome P450 (CYP) is a class of heme-containing enzymes which mainly catalyze a monooxygenation reaction of various chemicals, and hence CYP plays a key role in the drug metabolism. Although CYP2C19 isoform is a minor hepatic CYP, it metabolizes clinically important drugs such as omeprazole and S‑mephenytoin. In this work, the interaction of purified CYP2C19 WT (CYP2C19) with seven drugs (phenytoin, S‑mephenytoin, omeprazole, lansoprazole, cimetidine, propranolol, and warfarin) was investigated using spectroscopic methods. The binding of each drug and the induced structural change in the heme distal environment were evaluated. Ferric form of CYP2C19 was revealed to contain a six-coordinate low-spin heme with a water molecule as a sixth ligand in a distal site, and the addition of each drug caused varied minor fraction of five-coordinate heme. It was suggested that the ligated water molecule was partly moved away from the heme distal environment and that the degree of water removal was dependent on the type of drugs. The effect on the coordination was varied with the studied drugs with wide variation in the dissociation constants from 2.6 μM for lansoprazole to 5400 μM for warfarin. Phenytoin and S‑mephenytoin showed that binding to CYP2C19 occurred in a stepwise manner and that the coordination of a water molecule was facilitated in the second binding step. In the ferrous CO-bound state, ν(FeCO) stretching mode was clearly observed at 471 cm^-1 in the absence of drugs. The Raman line was greatly up-shifted by omeprazole (487 cm^-1) and lansoprazole (477 cm^-1) but was minimally affected by propranolol, phenytoin, and S‑mephenytoin. These results indicate that slight chemical modification of a drug greatly affects the heme distal environments upon binding.

Identification of the full 26 cytochrome P450 (CYP) genes and analysis of their expression in response to benzo[α]pyrene in the marine rotifer Brachionus rotundiformis

Cytochrome P450s (CYPs) are a large gene superfamily that are found in all living organisms. CYPs have a key role in detoxification of xenobiotics and endogenous chemicals. Although aquatic invertebrate CYPs and their detoxification mechanisms have been reported, little is known about interspecific comparison of CYPs and their detoxification mechanism in the rotifer Brachionus spp. The aim of this study was to identify the entire CYPs in the rotifer Brachionus rotundiformis (B. rotundiformis) and compare B. rotundiformis-CYPs to the previously reported CYPs in other model Brachionus spp. (B. koreanus, B. plicatilis, and B. calyciflorus). To validate the model, the rotifer, specifically Brachionus rotundiformis was exposed to various concentrations of B[α]P, which is widely used PAH xenobiotic, and analyzed gene expression in response to B[α]P. Here, in silico analysis results showed the total of 26 CYPs from the rotifer B. rotundiformis. Based on the phylogenetic analysis, the 26 B. rotundiformis-CYPs were separated into five different clans: 2, 3, 4, mitochondrial, and 46 clans in comparison to three rotifers species, B. koreanus, B. plicatilis, and B. calyciflorus. To understand the detoxification mechanisms of 26 B. rotundiformis-CYPs, we investigated transcriptional expression of 26 CYPs and found that five CYPs (CYP3045A2, CYP3045B4, CYP3045C10, CYP3049A5, and CYP3049E8) were significantly increased (P < 0.05) in response to 10 and 100 μg B[α]P. In addition, we identified the aryl hydrocarbon receptor (AhR) and aryl hydrocarbon receptor nuclear translocator (ARNT) and observed slight up-regulation of B. rotundiformis-AhR and -ARNT, indicating that these CYPs are likely associated with detoxification mechanism and could be used as potential molecular biomarkers of B[α]P in B. rotundiformis. Overall, this study will be helpful for expanding our knowledge of invertebrate CYPs on detoxification mechanisms associated with AhR signaling pathway in rotifers.

Non-classical 11β-hydroxylase deficiency caused by compound heterozygous mutations: a case study and literature review

BACKGROUND

11β-hydroxylase deficiency (11OHD) is extremely rare, and reports of non-classical 11OHD are even rarer. Non-classical 11OHD usually presents as premature adrenarche, hyperandrogenism, menstrual disorders, and hypertension. Because the symptoms of non-classical 11OHD are mild, delayed diagnosis or misdiagnosis as polycystic ovary syndrome or primary hypertension is common.

CASE PRESENTATION

This paper introduces a case of a young female patient presenting hypertension and menstrual disorders. Laboratory examination revealed increased androgen levels, mild adrenal hyperplasia, mild left ventricular hypertrophy, and mild sclerosis of the lower limb arteries. 11OHD was confirmed by genetic testing, and the patient was found to carry compound heterozygous mutations in CYP11B1 (c.583 T > C and c.1358G > A). The mutation Y195H is located in exon 3 and has not been reported previously. In silico studies indicated that this mutation may cause reduced enzymatic activity. After treatment with hydrocortisone and spironolactone, blood pressure was brought under good control, and menstruation returned to normal. We also conducted a retrospective review of previously reported cases in the literature (over 170 cases since 1991).

CONCLUSIONS

Early diagnosis of non-classical 11OHD is difficult because its symptoms are mild. The possibility of this disease should be considered in patients with early-onset hypertension, menstrual disorders, and hyperandrogenism to provide early treatment and prevent organ damage due to hypertension and hyperandrogenism. CYP11B1 mutations are known to be race-specific and are concentrated in exons 3 and 8, of which mutations in the former are mostly associated with non-classical 11OHD, whereas mutations in the latter are mostly found in classical 11OHD, characterized by severe loss of enzymatic activity.

Microsecond MD simulations of human CYP2D6 wild-type and five allelic variants reveal mechanistic insights on the function

Characterization of cytochrome P450 2D6 (CYP2D6) and the impact of the major identified allelic variants on the activity of one of the most dominating drug-metabolising enzymes is essential to increase drug safety and avoid adverse reactions. Microsecond molecular dynamics simulations have been performed to capture the dynamic signatures of this complex enzyme and five allelic variants with diverse enzymatic activity. In addition to the apo simulations, three substrates (bufuralol, veliparib and tamoxifen) and two inhibitors (prinomastat and quinidine) were included to explore their influence on the structure and dynamical features of the enzyme. Our results indicate that the altered enzyme activity can be attributed to changes in the hydrogen bonding network within the active site, and local structural differences in flexibility, position and shape of the binding pocket. In particular, the increased (CYP2D6*53) or the decreased (CYP2D6*17) activity seems to be related to a change in dynamics of mainly the BC loop due to a modified hydrogen bonding network around this region. In addition, the smallest active site volume was found for CYP2D6*4 (no activity). CYP2D6*2 (normal activity) showed no major differences in dynamic behaviour compared to the wild-type.

Molecular dynamics investigations of regioselectivity of anionic/aromatic substrates by a family of enzymes: a case study of diclofenac binding in CYP2C isoforms

The CYP2C subfamily is of particular importance in the metabolism of drugs, food toxins, and procarcinogens. Like other P450 subfamilies, 2C enzymes share a high sequence identity, but significantly contribute in different ways to hepatic capacity to metabolize drugs. They often metabolize the same substrate to more than one product with different catalytic sites. Because it is challenging to characterize experimentally, much still remains unknown about the reason for why the substrate regioselectivity of these closely related subfamily members is different. Here, we have investigated the structural features of CYP2C8, CYP2C9, and CYP2C19 bound with their shared substrate diclofenac to elucidate the underlying molecular mechanism for the substrate regioselectivity of CYP2C subfamily enzymes. The obtained results demonstrate how a sequence divergence for the active site residues causes heterogeneous variations in the secondary structures and in major tunnel selections, and further affects the shape and chemical properties of the substrate-binding site. Structural analysis and free energy calculations showed that the most important determinants of regioselectivity among the CYP2C isoforms are the geometrical features of the active sites, as well as the hydrogen bonds and the hydrophobic interactions, mainly presenting as the various locations of Arg108 and substitutions of Phe205 for Ile205 in CYP2C8. The MM-GB/SA calculations combined with PMF results accord well with the experimental KM values, bridging the gap between the theory and the experimentally observed results of binding affinity differences. The present study provides important insights into the structure-function relationships of CYP2C subfamily enzymes, the knowledge of ligand binding characteristics and key residue contributions could guide future experimental and computational work on the synthesis of drugs with better pharmacokinetic properties so that CYP interactions could be avoided.

ChannelsDB: database of biomacromolecular tunnels and pores

ChannelsDB (http://ncbr.muni.cz/ChannelsDB) is a database providing information about the positions, geometry and physicochemical properties of channels (pores and tunnels) found within biomacromolecular structures deposited in the Protein Data Bank. Channels were deposited from two sources; from literature using manual deposition and from a software tool automatically detecting tunnels leading to the enzymatic active sites and selected cofactors, and transmembrane pores. The database stores information about geometrical features (e.g. length and radius profile along a channel) and physicochemical properties involving polarity, hydrophobicity, hydropathy, charge and mutability. The stored data are interlinked with available UniProt annotation data mapping known mutation effects to channel-lining residues. All structures with channels are displayed in a clear interactive manner, further facilitating data manipulation and interpretation. As such, ChannelsDB provides an invaluable resource for research related to deciphering the biological function of biomacromolecular channels.

Genome-wide identification of 31 cytochrome P450 (CYP) genes in the freshwater rotifer Brachionus calyciflorus and analysis of their benzo[α]pyrene-induced expression patterns

While marine invertebrate cytochrome P450 (CYP) genes and their roles in detoxification mechanisms have been studied, little information is available regarding freshwater rotifer CYPs and their functions. Here, we used genomic sequences and RNA-seq databases to identify 31 CYP genes in the freshwater rotifer Brachionus calyciflorus. The 31 Bc-CYP genes with a few tandem duplications were clustered into CYP 2, 3, 4, mitochondrial, and 46 clans with two marine rotifers Brachionus plicatilis and Brachionus koreanus. To understand the molecular responses of these 31 Bc-CYP genes, we also examined their expression patterns in response to benzo[α]pyrene (B[α]P). Three Bc-CYP genes (Bc-CYP3044B3, Bc-CYP3049B4, Bc-CYP3049B6) were significantly upregulated (P<0.05) in response to B[α]P, suggesting that these CYP genes can be involved in detoxification in response to B[α]P exposure. These genes might be useful as biomarkers of B[α]P exposure in B. calyciflorus. Overall, our findings expand the repertoire of known CYPs and shed light on their potential roles in xenobiotic detoxification in rotifers.

Effects of water accommodated fractions (WAFs) of crude oil in two congeneric copepods Tigriopus sp

Oil pollution has deleterious effects on marine ecosystems. However, the toxicity of crude oil towards Antarctic marine organisms has not been well studied. We compared the deleterious effects of water accommodated fractions (WAFs) of crude oil on reproduction, intracellular reactive oxygen species (ROS) levels, and antioxidant enzymatic activity in Antarctic (Tigriopus kingsejongensis) and temperate (Tigriopus japonicus) copepods. Reproductive rates of T. kingsejongensis and T. japonicus were significantly reduced (P < 0.05) in response to WAFs. Furthermore, T. kingsejongensis showed elevated levels of ROS and higher antioxidant enzyme (glutathione peroxidase [GPx]) activity than T. japonicus in response to WAFs. CYP genes from congeneric copepods were identified and annotated to better understand molecular detoxification mechanisms. We observed significant up-regulation (P < 0.05) of Tk-CYP3024A3 and Tj-CYP3024A2 in response to WAFs, suggesting that CYP genes may contribute to the detoxification mechanism in response to WAF exposure. These finding also suggest that WAFs may induce oxidative stress, leading to reproductive impairment in copepods. Furthermore, Tk-CYP3024A3 and Tj-CYP3024A2 genes can be considered as potential biomarkers of WAF toxicity in the congeneric copepods T. kingsejongensis and T. japonicus. This study will be helpful for enhancing our knowledge on the harmful effects of WAFs in Antarctic and temperate copepods and provides insight into the underlying molecular mechanisms.

Molecular cloning and 3D model of first cytochrome P450 from CYP3A subfamily in saltwater crocodile (Crocodylus porosus)

Cytochrome P450s (CYPs) play critical role in oxidative metabolism of numerous xenobiotics and endogenous compounds. The first CYP3A subfamily member in saltwater crocodile has been cloned and modelled for three-dimensional (3D) structure. The full-length cDNA was obtained employing reverse transcription polymerase chain reaction (RT-PCR) strategy and rapid amplification of cDNA ends (RACE). The cDNA sequence of 1659 nucleotides includes 132 nucleotides from 5' untranslated region (UTR), an open reading frame of 1527 nucleotides encoding 509 amino acids designated as CYP3A163. The alignment of CYP3A163 sequence with CYP3A subfamily across the lineages exhibit the loss of 1 residue in birds and 7 residues in mammals in comparison to reptiles suggesting the adaptation processes during evolution. The amino acid identity of CYP3A163 with Alligator mississippiensis CYP3A77 and Homo sapiens CYP3A4 is 91% and 62% respectively. The 3D structure of CYP3A163 modelled using human CYP3A4 structure as a template with Phyre² software, represents high similarity with its functionally important motifs and catalytic domain. Both sequence and structure of CYP3A163 display the common and conserved features of CYP3A subfamily. Overall, this study provides primary molecular and structural data of CYP3A163 required to investigate the xenobiotic metabolism in saltwater crocodiles.

Genome-wide identification of 52 cytochrome P450 (CYP) genes in the copepod Tigriopus japonicus and their B[α]P-induced expression patterns

Cytochrome P450s (CYPs) are enzymes with a heme-binding domain that are found in all living organisms. CYP enzymes have important roles associated with detoxification of xenobiotics and endogenous compounds (e.g. steroids, fatty acids, and hormones). Although CYP enzymes have been reported in several invertebrates, including insects, little is known about copepod CYPs. Here, we identified the entire repertoire of CYP genes (n=52) from whole genome and transcriptome sequences of the benthic copepod Tigriopus japonicus, including a tandem duplication (CYP3026A3, CYP3026A4, CYP3026A5), and examined patterns of gene expression over various developmental stages and in response to benzo[α]pyrene (B[α]P) exposure. Through phylogenetic analysis, the 52 T. japonicus CYP genes were assigned to five distinct clans: CYP2 (22 genes), CYP3 (19 genes), CYP4 (two genes), CYP20 (one gene), and mitochondrial (eight genes). Developmental stage and gender-specific expression patterns of the 52 T. japonicus CYPs were analyzed. CYP3022A1 was constitutively expressed during all developmental stages. CYP genes in clans 2 and 3 were induced in response to B[α]P, suggesting that these differentially modulated CYP transcripts are likely involved in defense against exposure to B[α]P and other pollutants. This study enhances our understanding of the repertoire of CYP genes in copepods and of their potential role in development and detoxification in copepods.

Novel and prevalent CYP11B1 gene mutations in Turkish patients with 11-β hydroxylase deficiency

11β-Hydroxylase deficiency is the second most frequent type of congenital adrenal hyperplasia and is more common in those of Turkish descent than in other populations. The purpose of this study is to examine the spectrum of CYP11B1 gene mutations in Turkish patients with 11β-hydroxylase deficiency. Twenty-eight patients from 24 families, ages ranging from 0.1 to 7 years, were included in the study. Clinical diagnosis was based on virilization and high levels of 11-deoxycortisol. Twenty-six cases exhibited the classical and 2 cases the non-classical form. Mutation screening of 9 CYP11B1 exons was performed by direct DNA sequence analysis, specifically amplifying CYP11B1 gene fragments while avoiding simultaneous amplification of homologous CYP11B2 gene sequences. Seventeen different mutations were detected, 6 of which are novel (p.Gln189Hisfs*70, p.Glu198Gly, p.Thr318Lys, p.Gly446Ser, IVS8+5G>C and exon 3-5 del). All of the identified mutations resulted in the classical form with severe virilization, except for the p.Gly446Ser mutation, which caused a late-onset type of 11β-hydroxylase deficiency. The c.954G>A;p.Thr318Thr mutation was the most common in our cohort, with an allele frequency of 14.6%.Of the CYP11B1 gene mutations detected, 75% were found in exons 3, 5 and 7 and the half of the mutations were nonsense, splice site, deletion or insertion mutations, causing severe virilization in female patients. The findings are important for genetic counseling and the prenatal diagnosis of Turkish patients with 11β-hydroxylase deficiency.

CAVER: Algorithms for Analyzing Dynamics of Tunnels in Macromolecules

The biological function of a macromolecule often requires that a small molecule or ion is transported through its structure. The transport pathway often leads through void spaces in the structure. The properties of transport pathways change significantly in time; therefore, the analysis of a trajectory from molecular dynamics rather than of a single static structure is needed for understanding the function of pathways. The identification and analysis of transport pathways are challenging because of the high complexity and diversity of macromolecular shapes, the thermal motion of their atoms, and the large amount of conformations needed to properly describe conformational space of protein structure. In this paper, we describe the principles of the CAVER 3.0 algorithms for the identification and analysis of properties of transport pathways both in static and dynamic structures. Moreover, we introduce the improved clustering solution for finding tunnels in macromolecules, which is included in the latest CAVER 3.02 version. Voronoi diagrams are used to identify potential pathways in each snapshot of a molecular dynamics trajectory and clustering is then used to find the correspondence between tunnels from different snapshots. Furthermore, the geometrical properties of pathways and their evolution in time are computed and visualized.

The Role of Protein-Protein and Protein-Membrane Interactions on P450 Function

This symposium summary, sponsored by the ASPET, was held at Experimental Biology 2015 on March 29, 2015, in Boston, Massachusetts. The symposium focused on: 1) the interactions of cytochrome P450s (P450s) with their redox partners; and 2) the role of the lipid membrane in their orientation and stabilization. Two presentations discussed the interactions of P450s with NADPH-P450 reductase (CPR) and cytochrome b5. First, solution nuclear magnetic resonance was used to compare the protein interactions that facilitated either the hydroxylase or lyase activities of CYP17A1. The lyase interaction was stimulated by the presence of b5 and 17α-hydroxypregnenolone, whereas the hydroxylase reaction was predominant in the absence of b5. The role of b5 was also shown in vivo by selective hepatic knockout of b5 from mice expressing CYP3A4 and CYP2D6; the lack of b5 caused a decrease in the clearance of several substrates. The role of the membrane on P450 orientation was examined using computational methods, showing that the proximal region of the P450 molecule faced the aqueous phase. The distal region, containing the substrate-access channel, was associated with the membrane. The interaction of NADPH-P450 reductase (CPR) with the membrane was also described, showing the ability of CPR to "helicopter" above the membrane. Finally, the endoplasmic reticulum (ER) was shown to be heterogeneous, having ordered membrane regions containing cholesterol and more disordered regions. Interestingly, two closely related P450s, CYP1A1 and CYP1A2, resided in different regions of the ER. The structural characteristics of their localization were examined. These studies emphasize the importance of P450 protein organization to their function.

Cytochrome P450 20A1 in zebrafish: Cloning, regulation and potential involvement in hyperactivity disorders

Cytochrome P450 (CYP) enzymes for which there is no functional information are considered "orphan" CYPs. Previous studies showed that CYP20A1, an orphan, is expressed in human hippocampus and substantia nigra, and in zebrafish (Danio rerio) CYP20A1 maternal transcript occurs in eggs, suggesting involvement in brain and in early development. Moreover, hyperactivity is reported in humans with chromosome 2 microdeletions including CYP20A1. We examined CYP20A1 in zebrafish, including impacts of chemical exposure on expression. Zebrafish CYP20A1 cDNA was cloned, sequenced, and aligned with cloned human CYP20A1 and predicted vertebrate orthologs. CYP20A1s share a highly conserved N-terminal region and unusual sequences in the I-helix and the heme-binding CYP signature motifs. CYP20A1 mRNA expression was observed in adult zebrafish organs including the liver, heart, gonads, spleen and brain, as well as the eye and optic nerve. Putative binding sites in proximal promoter regions of CYP20A1s, and response of zebrafish CYP20A1 to selected nuclear and xenobiotic receptor agonists, point to up-regulation by agents involved in steroid hormone response, cholesterol and lipid metabolism. There also was a dose-dependent reduction of CYP20A1 expression in embryos exposed to environmentally relevant levels of methylmercury. Morpholino knockdown of CYP20A1 in developing zebrafish resulted in behavioral effects, including hyperactivity and a slowing of the optomotor response in larvae. The results suggest that altered expression of CYP20A1 might be part of a mechanism linking methylmercury exposure to neurobehavioral deficits. The expanded information on CYP20A1 brings us closer to "deorphanization", that is, identifying CYP20A1 functions and its roles in health and disease.

Incorporation of charged residues in the CYP2J2 F-G loop disrupts CYP2J2-lipid bilayer interactions

CYP2J2 epoxygenase is an extrahepatic, membrane bound cytochrome P450 (CYP) that is primarily found in the heart and mediates endogenous fatty acid metabolism. CYP2J2 interacts with membranes through an N-terminal anchor and various non-contiguous hydrophobic residues. The molecular details of the motifs that mediate membrane interactions are complex and not fully understood. To gain better insights of these complex protein-lipid interactions, we employed molecular dynamics (MD) simulations using a highly mobile membrane mimetic (HMMM) model that enabled multiple independent spontaneous membrane binding events to be captured. Simulations revealed that CYP2J2 engages with the membrane at the F-G loop through hydrophobic residues Trp-235, Ille-236, and Phe-239. To explore the role of these residues, three F-G loop mutants were modeled from the truncated CYP2J2 construct (Δ34) which included Δ34-I236D, Δ34-F239H and Δ34-I236D/F239H. Using the HMMM coordinates of CYP2J2, the simulations were extended to a full POPC membrane which showed a significant decrease in the depth of insertion for each of the F-G loop mutants. The CYP2J2 F-G loop mutants were expressed in E. coli and were shown to be localized to the cytosolic fraction at a greater percentage relative to construct Δ34. Notably, the functional data demonstrated that the double mutant, Δ34-I236D/F239H, maintained native-like enzymatic activity. The membrane insertion characteristics were examined by monitoring CYP2J2 Trp-quenching fluorescence spectroscopy upon binding nanodiscs containing pyrene phospholipids. Relative to the Δ34 construct, the F-G loop mutants exhibited lower Trp quenching and membrane insertion. Taken together, the results suggest that the mutants exhibit a different membrane topology in agreement with the MD simulations and provide important evidence towards the involvement of key residues in the F-G loop of CYP2J2.

Effects of tetracycline on developmental toxicity and molecular responses in zebrafish (Danio rerio) embryos

The extensive use of pharmaceuticals has resulted in the intensive contamination of water bodies. Tetracycline is a type of antibiotic and its potential toxicity is causing environmental concern. The effects of developmental toxicity and the mechanisms of tetracycline on fish embryos are not well understood. Zebrafish embryos are used in this study to investigate the developmental toxicity of this compound. Four hour post-fertilization (hpf) zebrafish embryos are exposed to different concentrations of tetracycline until 96 hpf. The larvae display developmental delay phenotypes, including hatching delay, shorter body length, increased yolk sac area and uninflated swim bladder upon exposure to tetracycline. Delayed yolk sac absorption and swim bladder deficiency at 96 hpf are observed in the zebrafish larvae upon exposure to 20 μg/L of tetracycline. To test whether tetracycline causes oxidative damage and the resulting oxidative stress-induced apoptosis, the generation of reactive oxygen species (ROS), Acridine Orange staining and real time polymerase chain reaction have been performed in this study. The results indicate that tetracycline exposure results in significant increases in ROS production and cell apoptosis, mainly in the tail areas at 96 hpf. The gene expression pattern demonstrates that tetracycline induces ROS which causes apoptosis in the zebrafish larvae, and the results also indicate that caspase-dependent apoptotic pathways may greatly contribute to tetracycline-induced apoptosis in the early-life stages of the zebrafish. In addition, we have investigated the effects of tetracycline on marker genes related to resistance mechanisms and gene regulating drug biotransformation. The results of these gene expression studies indicate that tetracycline could induce zebrafish to resist pharmaceuticals and Cytochrome P450s that are involved in the biotransformation of tetracycline in zebrafish larvae. The overall results indicate that tetracycline can produce oxidative stress and induce apoptosis, which brings about significant developmental delay in zebrafish embryos.

Recent progress on structural bioinformatics research of cytochrome P450 and its impact on drug discovery

Cytochrome P450 is predominantly responsible for human drug metabolism, which is of critical importance for drug discovery and development. Structural bioinformatics focuses on analysis and prediction of three-dimentional structure of biological macromolecules and elucidation of structure-function relationship as well as identification of important binding interactions. Rapid advancement of structural bioinformatics has been made over the last decade. With more information available for CYP structures, the methods of structural bioinformatics may be used in the CYP field. In this review, we demonstrate three previous studies on CYP using the methods of structural bioinformatics, including the investigation of reasons for decrease of enzymatic activity of CYP1A2 caused by a peripheral mutation, the construction of a pharmacophore model specific to active site of CYP1A2 and the prediction of the functional consequences of single residue mutation in CYP. By illustrating these studies we attempt to show the potential role of structural bioinformatics in CYP research and help better understanding the importance of structural bioinformatics in drug designing.

Assessing orally bioavailable commercial silver nanoparticle product on human cytochrome P450 enzyme activity

Nanotechnology produces a wide range of medicinal compounds, including nanoparticulate silver, which are increasingly introduced in various forms for consumer use. As with all medicinal compounds, potential drug interactions are an important consideration for ingested silver nanoparticles. Nanoparticulate silver-drug interactions may be mediated through induced oxidative stress in liver tissue where the majority of systemically bioavailable silver nanoparticles is found. To investigate whether an orally ingested commercially available colloidal silver nanoproduct produces pharmacokinetic interference on select cytochrome P450 enzymes, a prospective, single-blind, controlled in vivo human study using simultaneous administration of standardized probes for P450 enzyme classes CYP1A2, CYP2C9, CYP2C19, CYP2D6 and CYP3A4 was conducted. Oral ingestion of a commercial colloidal silver nanoproduct produces detectable silver in human serum after 14 days of dosing. This silver, however, elicits no demonstrable clinically significant changes in metabolic, hematologic, urinary, physical findings or cytochrome P450 enzyme inhibition or induction activity. Given their increasingly broad, diverse human exposures, future characterization of human cytochrome P450 enzyme activity for other systemically bioavailable nanotechnology products are warranted.

Ensemble generation and the influence of protein flexibility on geometric tunnel prediction in cytochrome P450 enzymes

Computational prediction of ligand entry and egress paths in proteins has become an emerging topic in computational biology and has proven useful in fields such as protein engineering and drug design. Geometric tunnel prediction programs, such as Caver3.0 and MolAxis, are computationally efficient methods to identify potential ligand entry and egress routes in proteins. Although many geometric tunnel programs are designed to accommodate a single input structure, the increasingly recognized importance of protein flexibility in tunnel formation and behavior has led to the more widespread use of protein ensembles in tunnel prediction. However, there has not yet been an attempt to directly investigate the influence of ensemble size and composition on geometric tunnel prediction. In this study, we compared tunnels found in a single crystal structure to ensembles of various sizes generated using different methods on both the apo and holo forms of cytochrome P450 enzymes CYP119, CYP2C9, and CYP3A4. Several protein structure clustering methods were tested in an attempt to generate smaller ensembles that were capable of reproducing the data from larger ensembles. Ultimately, we found that by including members from both the apo and holo data sets, we could produce ensembles containing less than 15 members that were comparable to apo or holo ensembles containing over 100 members. Furthermore, we found that, in the absence of either apo or holo crystal structure data, pseudo-apo or -holo ensembles (e.g. adding ligand to apo protein throughout MD simulations) could be used to resemble the structural ensembles of the corresponding apo and holo ensembles, respectively. Our findings not only further highlight the importance of including protein flexibility in geometric tunnel prediction, but also suggest that smaller ensembles can be as capable as larger ensembles at capturing many of the protein motions important for tunnel prediction at a lower computational cost.

Characterisation of three novel CYP11B1 mutations in classic and non-classic 11β-hydroxylase deficiency

BACKGROUND

Congenital adrenal hyperplasia (CAH) is one of the most common autosomal recessive inherited endocrine diseases. Steroid 11β-hydroxylase (P450c11) deficiency (11OHD) is the second most common form of CAH.

AIM

The aim of the study was to study the functional consequences of three novel CYP11B1 gene mutations (p.His125Thrfs*8, p.Leu463_Leu464dup and p.Ser150Leu) detected in patients suffering from 11OHD and to correlate this data with the clinical phenotype.

METHODS

Functional analyses were done by using a HEK293 cell in vitro expression system comparing WT with mutant P450c11 activity. Mutant proteins were examined in silico to study their effect on the three-dimensional structure of the protein.

RESULTS

Two mutations (p.His125Thrfs*8 and p.Leu463_Leu464dup) detected in patients with classic 11OHD showed a complete loss of P450c11 activity. The mutation (p.Ser150Leu) detected in a patient with non-classic 11OHD showed partial functional impairment with 19% of WT activity.

CONCLUSION

Functional mutation analysis enables the correlation of novel CYP11B1 mutations to the classic and non-classic 11OHD phenotype respectively. Mutations causing a non-classic phenotype show typically partial impairment due to reduced maximum reaction velocity comparable with non-classic mutations in 21-hydroxylase deficiency. The increasing number of mutations associated with non-classic 11OHD illustrate that this disease should be considered as diagnosis in patients with otherwise unexplained hyperandrogenism.

Molecular dynamic simulations reveal the mechanism of binding between xanthine inhibitors and DPP-4

We apply molecular docking, molecular dynamics (MD) simulation, and binding free energy calculation to investigate and reveal the binding mechanism between five xanthine inhibitors and DPP-4. The electrostatic and van der Waals interactions of the five inhibitors with DPP-4 are analyzed and discussed. The computed binding free energies using MM-PBSA method are in qualitatively agreement with experimental inhibitory potency of five inhibitors. The hydrogen bonds of inhibitors with Ser630 and Asp663 can stabilize the inhibitors in binding sites. The van der Waals interactions, especially the key contacts with His740, Asn710, Trp629, and Tyr666 have larger contributions to the binding free energy and play important roles in distinguishing the variant bioactivity of five inhibitors.

Behavior of human cytochromes P450 on lipid membranes

Human cytochromes P450 (CYPs) are membrane-anchored enzymes involved in biotransformation of many marketed drugs. We constructed atomic models of six human CYPs (CYP1A2, 2A6, 2C9, 2D6, 2E1, and 3A4) anchored to a lipid bilayer to investigate the positions and orientations of CYPs on a membrane. We equilibrated the models by molecular dynamics simulations on a 100+ ns time scale. Catalytic domains of all studied CYPs were found to be partially immersed in the lipid bilayer, whereas the N-terminal part and F'/G' loop are deeply immersed. The proximal side of the enzyme faces the cytosol, whereas the distal side, where openings of substrate access and product release channels to the active site are primarily located, points toward the lipid bilayer. Access channels with openings in the vicinity of the B/C and F/G loops are typically positioned below the lipid head groups, whereas the solvent channel points toward the membrane-water interface. We found that the access channel opening positions match the preferred substrate positions, whereas the product release channel exit positions correspond closely with the positions of the products. This may indicate that membrane-anchored CYPs have evolutionarily adapted to facilitate uptake of nonpolar substrates from the membrane and uptake/release of polar substrates or products from/to the membrane-water interface.

Congenital adrenal hyperplasia due to 11-beta-hydroxylase deficiency: functional consequences of four CYP11B1 mutations

Congenital adrenal hyperplasia (CAH) is one of the most common autosomal recessive inherited endocrine disease. Steroid 11β-hydroxylase deficiency (11β-OHD) is the second most common form of CAH. The aim of the study was to study the functional consequences of three novel and one previously described CYP11B1 gene mutations (p.(Arg143Trp), p.(Ala306Val), p.(Glu310Lys) and p.(Arg332Gln)) detected in patients suffering from classical and non-classical 11β-OHD. Functional analyses were performed by using a HEK293 cell in vitro expression system comparing wild type (WT) with mutant 11β-hydroxylase activity. Mutant proteins were examined in silico to study their effect on the three-dimensional structure of the protein. Two mutations (p.(Ala306Val) and p.(Glu310Lys)) detected in patients with classical 11β-OHD showed a nearly complete loss of 11β-hydroxylase activity. The mutations p.(Arg143Trp) and p.(Arg332Gln) detected in patients with non-classical 11β-OHD showed a partial functional impairment with approximately 8% and 6% of WT activity, respectively. Functional mutation analysis allows the classification of novel CYP11B1 mutations as causes of classical and non-classical 11β-OHD. The detection of patients with non-classical phenotypes underscores the importance to screen patients with a phenotype comparable to non-classical 21-hydroxylase deficiency for mutations in the CYP11B1 gene in case of a negative analysis of the CYP21A2 gene. As CYP11B1 mutations are most often individual for a family, the in vitro analysis of novel mutations is essential for clinical and genetic counselling.

MOLE 2.0: advanced approach for analysis of biomacromolecular channels

BACKGROUND

Channels and pores in biomacromolecules (proteins, nucleic acids and their complexes) play significant biological roles, e.g., in molecular recognition and enzyme substrate specificity.

RESULTS

We present an advanced software tool entitled MOLE 2.0, which has been designed to analyze molecular channels and pores. Benchmark tests against other available software tools showed that MOLE 2.0 is by comparison quicker, more robust and more versatile. As a new feature, MOLE 2.0 estimates physicochemical properties of the identified channels, i.e., hydropathy, hydrophobicity, polarity, charge, and mutability. We also assessed the variability in physicochemical properties of eighty X-ray structures of two members of the cytochrome P450 superfamily.

CONCLUSION

Estimated physicochemical properties of the identified channels in the selected biomacromolecules corresponded well with the known functions of the respective channels. Thus, the predicted physicochemical properties may provide useful information about the potential functions of identified channels. The MOLE 2.0 software is available at http://mole.chemi.muni.cz.

Structure and function of heme proteins in non-native states: a mini-review

Heme proteins perform various biological functions ranging from electron transfer, oxygen binding and transport, catalysis, to signaling. Although adopting proper native states is very important for these functions, progresses in representative heme proteins, including cytochrome c (cyt c), cytochrome b5 (cyt b5), myoglobin (Mb), neuroglobin (Ngb), cytochrome P450 (CYP) and heme-based sensor proteins such as CO sensor CooA, showed that various native functions, or new functions evolved, are also closely associated with non-native states. The structure and function relationship of heme proteins in non-native states is thus as important as that in native states for elucidating the precise roles of heme proteins in biological systems.