Hexose-6-phosphate dehydrogenase modulates 11beta-hydroxysteroid dehydrogenase type 1-dependent metabolism of 7-keto- and 7beta-hydroxy-neurosteroids

BACKGROUND: The role of 11beta-hydroxysteroid dehydrogenase type 1 (11beta-HSD1) in the regulation of energy metabolism and immune system by locally reactivating glucocorticoids has been extensively studied. Experiments determining initial rates of enzyme activity revealed that 11beta-HSD1 can catalyze both the reductase and the dehydrogenase reaction in cell lysates, whereas it predominantly catalyzes the reduction of cortisone to cortisol in intact cells that also express hexose-6-phosphate dehydrogenase (H6PDH), which provides cofactor NADPH. Besides its role in glucocorticoid metabolism, there is evidence that 11beta-HSD1 is involved in the metabolism of 7-keto- and 7-hydroxy-steroids; however the impact of H6PDH on this alternative function of 11beta-HSD1 has not been assessed. METHODOLOGY: We investigated the 11beta-HSD1-dependent metabolism of the neurosteroids 7-keto-, 7alpha-hydroxy- and 7beta-hydroxy-dehydroepiandrosterone (DHEA) and 7-keto- and 7beta-hydroxy-pregnenolone, respectively, in the absence or presence of H6PDH in intact cells. 3D-structural modeling was applied to study the binding of ligands in 11beta-HSD1. PRINCIPAL FINDINGS: We demonstrated that 11beta-HSD1 functions in a reversible way and efficiently catalyzed the interconversion of these 7-keto- and 7-hydroxy-neurosteroids in intact cells. In the presence of H6PDH, 11beta-HSD1 predominantly converted 7-keto-DHEA and 7-ketopregnenolone into their corresponding 7beta-hydroxy metabolites, indicating a role for H6PDH and 11beta-HSD1 in the local generation of 7beta-hydroxy-neurosteroids. 3D-structural modeling offered an explanation for the preferred formation of 7beta-hydroxy-neurosteroids. CONCLUSIONS: Our results from experiments determining the steady state concentrations of glucocorticoids or 7-oxygenated neurosteroids suggested that the equilibrium between cortisone and cortisol and between 7-keto- and 7-hydroxy-neurosteroids is regulated by 11beta-HSD1 and greatly depends on the coexpression with H6PDH. Thus, the impact of H6PDH on 11beta-HSD1 activity has to be considered for understanding both glucocorticoid and neurosteroid action in different tissues.

The cytochrome P450 aromatase lacking exon 5 is associated with a phenotype of nonclassic aromatase deficiency and is also present in normal human steroidogenic tissues

OBJECTIVE: The previously described c655G>A mutation of the human cytochrome P450 aromatase gene (P450aro, CYP19) results in aberrant splicing due to disruption of a donor splice site. To explain the phenotype of partial aromatase deficiency observed in a female patient described with this mutation, molecular consequences of the c655G>A mutation were investigated. DESIGN: To investigate whether the c655G>A mutation causes an aberrant spliced mRNA lacking exon 5 (-Ex5), P450aro RNA was analysed from the patient's lymphocytes by reverse transcription polymerase chain reaction (RT-PCR) and by splicing assays performed in Y1 cells transfected with a P450aro -Ex5 expression vector. Aromatase activity of the c655G>A mutant was predicted by three dimensional (3D) protein modelling studies and analysed in transiently transfected Y1 cells. Exon 5 might be predicted as a poorly defined exon suggesting a susceptibility to both splicing mutations and physiological alternative splicing events. Therefore, expression of the -Ex5 mRNA was also assessed as a possibly naturally occurring alternative splicing transcript in normal human steroidogenic tissues. PATIENTS: An aromatase deficient girl was born with ambiguous genitalia. Elevated serum LH, FSH and androgens, as well as cystic ovaries, were found during prepuberty. At the age of 8.4 years, spontaneous breast development and a 194.6 pmol/l serum oestradiol level was observed. RESULTS: The -Ex5 mRNA was found in lymphocytes of the P450aro deficient girl and her father, who was a carrier of the mutation. Mutant minigene expression resulted in complete exon 5 skipping. As expected from 3D protein modelling, -Ex5 cDNA expression in Y1 cells resulted in loss of P450aro activity. In addition, the -Ex5 mRNA was present in placenta, prepubertal testis and adrenal tissues. CONCLUSIONS: Alternative splicing of exon 5 of the CYP19 gene occurs in the wild type (WT) as well as in the c655G>A mutant. We speculate that for the WT it might function as a regulatory mechanism for aromatization, whereas for the mutant a relative prevalence of the shorter over the full-length protein might explain the phenotype of partial aromatase deficiency.

Whole-exome sequencing reveals a C-terminal germline variant in CEBPA-associated acute myeloid leukemia: 45-year follow-up of a large family.

Familial acute myeloid leukemia is rare and linked to germline mutations in RUNX1, GATA2 or CCAAT/enhancer binding protein-α (CEBPA). We re-evaluated a large family with acute myeloid leukemia originally seen at NIH in 1969. We utilized whole-exome sequencing to study this family, and conducted in silico bioinformatics analysis, protein structural modeling and laboratory experiments to assess the impact of the identified CEBPA Q311P mutation. Unlike most previously identified germline mutations in CEBPA, which were N-terminal frameshift mutations, we identified a novel Q311P variant that was located in the C-terminal bZip domain of C/EBPα. Protein structural modeling suggested that the Q311P mutation alters the ability of the CEBPA dimer to bind DNA. Electrophoretic mobility shift assays showed that the Q311P mutant had attenuated binding to DNA, as predicted by the protein modeling. Consistent with these findings, we found that the Q311P mutation has reduced transactivation, consistent with a loss-of-function mutation. From 45 years of follow-up, we observed incomplete penetrance (46%) of CEBPA Q311P. This study of a large multi-generational pedigree reveals that a germline mutation in the C-terminal bZip domain can alter the ability of C/EBP-α to bind DNA and reduces transactivation, leading to acute myeloid leukemia.

Clinical and biochemical description of a novel CYP21A2 gene mutation 962_963insA using a new 3D model for the P450c21 protein

OBJECTIVE: A severely virilized 46, XX newborn girl was referred to our center for evaluation and treatment of congenital adrenal hyperplasia (CAH) because of highly elevated 17alpha-hydroxyprogesterone levels at newborn screening; biochemical tests confirmed the diagnosis of salt-wasting CAH. Genetic analysis revealed that the girl was compound heterozygote for a previously reported Q318X mutation in exon 8 and a novel insertion of an adenine between nucleotides 962 and 963 in exon 4 of the CYP21A2 gene. This 962_963insA mutation created a frameshift leading to a stop codon at amino acid 161 of the P450c21 protein. AIM AND METHODS: To better understand structure-function relationships of mutant P450c21 proteins, we performed multiple sequence alignments of P450c21 with three mammalian P450s (P450 2C8, 2C9 and 2B4) with known structures as well as with human P450c17. Comparative molecular modeling of human P450c21 was then performed by MODELLER using the X-ray crystal structure of rabbit P450 2B4 as a template. RESULTS: The new three dimensional model of human P450c21 and the sequence alignment were found to be helpful in predicting the role of various amino acids in P450c21, especially those involved in heme binding and interaction with P450 oxidoreductase, the obligate electron donor. CONCLUSION: Our model will help in analyzing the genotype-phenotype relationship of P450c21 mutations which have not been tested for their functional activity in an in vitro assay.

3D-modeling of Mercury's solar wind sputtered surface-exosphere environment

The efficiency of sputtered refractory elements by H+ and He++ solar wind ions from Mercury's surface and their contribution to the exosphere are studied for various solar wind conditions. A 3D solar wind-planetary interaction hybrid model is used for the evaluation of precipitation maps of the sputter agents on Mercury's surface. By assuming a global mineralogical surface composition, the related sputter yields are calculated by means of the 2013 SRIM code and are coupled with a 3D exosphere model. Because of Mercury's magnetic field, for quiet and nominal solar wind conditions the plasma can only precipitate around the polar areas, while for extreme solar events (fast solar wind, coronal mass ejections, interplanetary magnetic clouds) the solar wind plasma has access to the entire dayside. In that case the release of particles form the planet's surface can result in an exosphere density increase of more than one order of magnitude. The corresponding escape rates are also about an order of magnitude higher. Moreover, the amount of He++ ions in the precipitating solar plasma flow enhances also the release of sputtered elements from the surface in the exosphere. A comparison of our model results with MESSENGER observations of sputtered Mg and Ca elements in the exosphere shows a reasonable quantitative agreement. (C) 2015 Elsevier Ltd. All rights reserved.

Block of the hERG channel by bupivacaine: Electrophysiological and modeling insights towards stereochemical optimization

The hERG voltage-gated potassium channel mediates the cardiac I(Kr) current, which is crucial for the duration of the cardiac action potential. Undesired block of the channel by certain drugs may prolong the QT interval and increase the risk of malignant ventricular arrhythmias. Although the molecular determinants of hERG block have been intensively studied, not much is known about its stereoselectivity. Levo-(S)-bupivacaine was the first drug reported to have a higher affinity to block hERG than its enantiomer. This study strives to understand the principles underlying the stereoselectivity of bupivacaine block with the help of mutagenesis analyses and molecular modeling simulations. Electrophysiological measurements of mutated hERG channels allowed for the identification of residues involved in bupivacaine binding and stereoselectivity. Docking and molecular mechanics simulations for both enantiomers of bupivacaine and terfenadine (a non-stereoselective blocker) were performed inside an open-state model of the hERG channel. The predicted binding modes enabled a clear depiction of ligand-protein interactions. Estimated binding affinities for both enantiomers were consistent with electrophysiological measurements. A similar computational procedure was applied to bupivacaine enantiomers towards two mutated hERG channels (Tyr652Ala and Phe656Ala). This study confirmed, at the molecular level, that bupivacaine stereoselectively binds the hERG channel. These results help to lay the foundation for structural guidelines to optimize the cardiotoxic profile of drug candidates in silico.

Effect of repetitiveness on the immunogenicity and antigenicity of Trypanosoma cruzi FRA protein

Repetitive proteins (RP) of Trypanosoma cruzi are highly present in the parasite and are strongly recognized by sera from Chagas' disease patients. Flagelar Repetitive Antigen (FRA), which is expressed in all steps of the parasite life cycle, is the RP that displays the greatest number of aminoacids per repeat and has been indicated as one of the most suitable candidate for diagnostic test because of its high performance in immunoassays. Here we analyzed the influence of the number of repeats on the immunogenic and antigenic properties of the antigen. Recombinant proteins containing one, two, and four tandem repeats of FRA (FRA1, FRA2, and FRA4, respectively) were obtained and the immune response induced by an equal amount of repeats was evaluated in a mouse model. The reactivity of specific antibodies present in sera from patients naturally infected with T. cruzi was also assessed against FRA1, FRA2, and FRA4 proteins, and the relative avidity was analyzed. We determined that the number of repeats did not increase the humoral response against the antigen and this result was reproduced when the repeated motifs were alone or fused to a non-repetitive protein. By contrast, the binding affinity of specific human antibodies increases with the number of repeated motifs in FRA antigen. We then concluded that the high ability of FRA to be recognized by specific antibodies from infected individuals is mainly due to a favorable polyvalent interaction between the antigen and the antibodies. In accordance with experimental results, a 3D model was proposed and B epitope in FRA1, FRA2, and FRA4 were predicted.

A tool for the qualitative comparison of membrane-embedded and detergent-solubilized membrane protein structures in projection

The calculation of projection structures (PSs) from Protein Data Bank (PDB)-coordinate files of membrane proteins is not well-established. Reports on such attempts exist but are rare. In addition, the different procedures are barely described and thus difficult if not impossible to reproduce. Here we present a simple, fast and well-documented method for the calculation and visualization of PSs from PDB-coordinate files of membrane proteins: the projection structure visualization (PSV)-method. The PSV-method was successfully validated using the PS of aquaporin-1 (AQP1) from 2D crystals and cryo-transmission electron microscopy, and the PDB-coordinate file of AQP1 determined from 3D crystals and X-ray crystallography. Besides AQP1, which is a relatively rigid protein, we also studied a flexible membrane transport protein, i.e. the L-arginine/agmatine antiporter AdiC. Comparison of PSs calculated from the existing PDB-coordinate files of substrate-free and L-arginine-bound AdiC indicated that conformational changes are detected in projection. Importantly, structural differences were found between the PSV-method calculated PSs of the detergent-solubilized AdiC proteins and the PS from cryo-TEM of membrane-embedded AdiC. These differences are particularly exciting since they may reflect a different conformation of AdiC induced by the lateral pressure in the lipid bilayer.

Mutagenesis and computer modeling studies of a GPCR conserved residue W5.43(194) in ligand recognition and signal transduction for CB2 receptor

W5.43(194), a conserved tryptophan residue among G-protein coupled receptors (GPCRs) and cannabinoid receptors (CB), was examined in the present report for its significance in CB2 receptor ligand binding and adenylyl cyclase (AC) activity. Computer modeling postulates that this site in CB2 may be involved in the affinity of WIN55212-2 and SR144528 through aromatic contacts. In the present study, we reported that a CB2 receptor mutant, W5.43(194)Y, which had a tyrosine (Y) substitution for tryptophan (W), retained the binding affinity for CB agonist CP55940, but reduced binding affinity for CB2 agonist WIN55212-2 and inverse agonist SR144528 by 8-fold and 5-fold, respectively; the CB2 W5.43(194)F and W5.43(194)A mutations significantly affect the binding activities of CP55940, WIN55212-2 and SR144528. Furthermore, we found that agonist-mediated inhibition of the forskolin-induced cAMP production was dramatically diminished in the CB2 mutant W5.43(194)Y, whereas W5.43(194)F and W5.43(194)A mutants resulted in complete elimination of downstream signaling, suggesting that W5.43(194) was essential for the full activation of CB2. These results indicate that both aromatic interaction and hydrogen bonding are involved in ligand binding for the residue W5.43(194), and the mutations of this tryptophan site may affect the conformation of the ligand binding pocket and therefore control the active conformation of the wild type CB2 receptor. W5.43(194)Y/F/A mutations also displayed noticeable enhancement of the constitutive activation probably attributed to the receptor conformational changes resulted from the mutations.

Protein and lipid binding parameters in rainbow trout (Oncorhynchus mykiss) blood and liver fractions to extrapolate from an in vitro metabolic degradation assay to in vivo bioaccumulation potential of hydrophobic organic chemicals

Binding of hydrophobic chemicals to colloids such as proteins or lipids is difficult to measure using classical microdialysis methods due to low aqueous concentrations, adsorption to dialysis membranes and test vessels, and slow kinetics of equilibration. Here, we employed a three-phase partitioning system where silicone (polydimethylsiloxane, PDMS) serves as a third phase to determine partitioning between water and colloids and acts at the same time as a dosing device for hydrophobic chemicals. The applicability of this method was demonstrated with bovine serum albumin (BSA). Measured binding constants (K(BSAw)) for chlorpyrifos, methoxychlor, nonylphenol, and pyrene were in good agreement with an established quantitative structure-activity relationship (QSAR). A fifth compound, fluoxypyr-methyl-heptyl ester, was excluded from the analysis because of apparent abiotic degradation. The PDMS depletion method was then used to determine partition coefficients for test chemicals in rainbow trout (Oncorhynchus mykiss) liver S9 fractions (K(S9w)) and blood plasma (K(bloodw)). Measured K(S9w) and K(bloodw) values were consistent with predictions obtained using a mass-balance model that employs the octanol-water partition coefficient (K(ow)) as a surrogate for lipid partitioning and K(BSAw) to represent protein binding. For each compound, K(bloodw) was substantially greater than K(S9w), primarily because blood contains more lipid than liver S9 fractions (1.84% of wet weight vs 0.051%). Measured liver S9 and blood plasma binding parameters were subsequently implemented in an in vitro to in vivo extrapolation model to link the in vitro liver S9 metabolic degradation assay to in vivo metabolism in fish. Apparent volumes of distribution (V(d)) calculated from the experimental data were similar to literature estimates. However, the calculated binding ratios (f(u)) used to relate in vitro metabolic clearance to clearance by the intact liver were 10 to 100 times lower than values used in previous modeling efforts. Bioconcentration factors (BCF) predicted using the experimental binding data were substantially higher than the predicted values obtained in earlier studies and correlated poorly with measured BCF values in fish. One possible explanation for this finding is that chemicals bound to proteins can desorb rapidly and thus contribute to metabolic turnover of the chemicals. This hypothesis remains to be investigated in future studies, ideally with chemicals of higher hydrophobicity.

Besnoitia besnoiti protein disulfide isomerase (BbPDI): molecular characterization, expression and in silico modelling

Besnoitia besnoiti is an apicomplexan parasite responsible for bovine besnoitiosis, a disease with a high prevalence in tropical and subtropical regions and re-emerging in Europe. Despite the great economical losses associated with besnoitiosis, this disease has been underestimated and poorly studied, and neither an effective therapy nor an efficacious vaccine is available. Protein disulfide isomerase (PDI) is an essential enzyme for the acquisition of the correct three-dimensional structure of proteins. Current evidence suggests that in Neosporacaninum and Toxoplasmagondii, which are closely related to B. besnoiti, PDI play an important role in host cell invasion, is a relevant target for the host immune response, and represents a promising drug target and/or vaccine candidate. In this work, we present the nucleotide sequence of the B. besnoiti PDI gene. BbPDI belongs to the thioredoxin-like superfamily (cluster 00388) and is included in the PDI_a family (cluster defined cd02961) and the PDI_a_PDI_a'_c subfamily (cd02995). A 3D theoretical model was built by comparative homology using Swiss-Model server, using as a template the crystallographic deduced model of Tapasin-ERp57 (PDB code 3F8U chain C). Analysis of the phylogenetic tree for PDI within the phylum apicomplexa reinforces the close relationship among B. besnoiti, N. caninum and T. gondii. When subjected to a PDI-assay based on the polymerisation of reduced insulin, recombinant BbPDI expressed in E. coli exhibited enzymatic activity, which was inhibited by bacitracin. Antiserum directed against recombinant BbPDI reacted with PDI in Western blots and by immunofluorescence with B. besnoiti tachyzoites and bradyzoites.

Statistical modeling of the eye for multimodal treatment planning for external beam radiation therapy of intraocular tumors

Ocular anatomy and radiation-associated toxicities provide unique challenges for external beam radiation therapy. For treatment planning, precise modeling of organs at risk and tumor volume are crucial. Development of a precise eye model and automatic adaptation of this model to patients' anatomy remain problematic because of organ shape variability. This work introduces the application of a 3-dimensional (3D) statistical shape model as a novel method for precise eye modeling for external beam radiation therapy of intraocular tumors.

Correspondence establishment in statistical modeling of shapes with arbitrary topology

Correspondence establishment is a key step in statistical shape model building. There are several automated methods for solving this problem in 3D, but they usually can only handle objects with simple topology, like that of a sphere or a disc. We propose an extension to correspondence establishment over a population based on the optimization of the minimal description length function, allowing considering objects with arbitrary topology. Instead of using a fixed structure of kernel placement on a sphere for the systematic manipulation of point landmark positions, we rely on an adaptive, hierarchical organization of surface patches. This hierarchy can be built on surfaces of arbitrary topology and the resulting patches are used as a basis for a consistent, multi-scale modification of the surfaces' parameterization, based on point distribution models. The feasibility of the approach is demonstrated on synthetic models with different topologies.

Infantile hypophosphatasia due to a new compound heterozygous TNSALP mutation - functional evidence for a hydrophobic side-chain?

BACKGROUND: Infantile hypophosphatasia (IH) is an inherited disorder characterized by defective bone mineralization and a deficiency of alkaline phosphatase activity. OBJECTIVE/DESIGN: The aim of the study was to evaluate a new compound heterozygous TNSALP mutation for its residual enzyme activity and localization of the comprised amino acid residues in a 3D-modeling. PATIENT: We report on a 4-week old girl with craniotabes, severe defects of ossification, and failure to thrive. Typical clinical features as low serum alkaline phosphatase, high serum calcium concentration, increased urinary calcium excretion, and nephrocalcinosis were observed. Vitamin D was withdrawn and the patient was started on calcitonin and hydrochlorothiazide. Nonetheless, the girl died at the age of 5 months from respiratory failure. RESULTS: Sequence analysis of the patient's TNSALP gene revealed two heterozygous mutations [c.653T>C (I201T), c.1171C>T (R374C)]. Transfection studies of the unique I201T variant in COS-7 cells yielded a mutant TNSALP protein with only a residual enzyme activity (3.7%) compared with wild-type, whereas the R374C variant was previously shown to reduce normal activity to 10.3%. 3D-modeling of the mutated enzyme showed that I201T resides in a region that does not belong to any known functional site. CONCLUSION: We note that I201, which has been conserved during evolution, is buried in a hydrophobic pocket and, therefore, the I>T-change should affect its functional properties. Residue R374C is located in the interface between monomers and it has been previously suggested that this mutation affects dimerization. These findings explain the patient's clinical picture and severe course.