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.

A novel FIP1L1-PDGFRA mutant destabilizing the inactive conformation of the kinase domain in chronic eosinophilic leukemia/hypereosinophilic syndrome

BACKGROUND: The Fip1-like-1-platelet-derived growth factor receptor alpha (FIP1L1-PDGFRA) gene fusion is a common cause of chronic eosinophilic leukemia (CEL)/hypereosinophilic syndrome (HES), and patients suffering from this particular subgroup of CEL/HES respond to low-dose imatinib therapy. However, some patients may develop imatinib resistance because of an acquired T674I mutation, which is believed to prevent drug binding through steric hindrance. METHODS: In an imatinib resistant FIP1L1-PDGFRA positive patient, we analyzed the molecular structure of the fusion gene and analyzed the effect of several kinase inhibitors on FIP1L1-PDGFRA-mediated proliferative responses in vitro. RESULTS: Sequencing of the FIP1L1-PDGFRA fusion gene revealed the occurrence of a S601P mutation, which is located within the nucleotide binding loop. In agreement with the clinical observations, imatinib did not inhibit the proliferation of S601P mutant FIP1L1-PDGFRA-transduced Ba/F3 cells. Moreover, sorafenib, which has been described to inhibit T674I mutant FIP1L1-PDGFRA, failed to block S601P mutant FIP1L1-PDGFRA. Structural modeling revealed that the newly identified S601P mutated form of PDGFRA destabilizes the inactive conformation of the kinase domain that is necessary to bind imatinib as well as sorafenib. CONCLUSIONS: We identified a novel mutation in FIP1L1-PDGFRA resulting in both imatinib and sorafenib resistance. The identification of novel drug-resistant FIP1L1-PDGFRA variants may help to develop the next generation of target-directed compounds for CEL/HES and other leukemias.

Understanding the role of argininosuccinate lyase transcript variants in the clinical and biochemical variability of the urea cycle disorder argininosuccinic aciduria

Argininosuccinic aciduria (ASA) is an autosomal recessive urea cycle disorder caused by deficiency of argininosuccinate lyase (ASL) with a wide clinical spectrum from asymptomatic to severe hyperammonemic neonatal onset life-threatening courses. We investigated the role of ASL transcript variants in the clinical and biochemical variability of ASA. Recombinant proteins for ASL wild type, mutant p.E189G, and the frequently occurring transcript variants with exon 2 or 7 deletions were (co-)expressed in human embryonic kidney 293T cells. We found that exon 2-deleted ASL forms a stable truncated protein with no relevant activity but a dose-dependent dominant negative effect on enzymatic activity after co-expression with wild type or mutant ASL, whereas exon 7-deleted ASL is unstable but seems to have, nevertheless, a dominant negative effect on mutant ASL. These findings were supported by structural modeling predictions for ASL heterotetramer/homotetramer formation. Illustrating the physiological relevance, the predominant occurrence of exon 7-deleted ASL was found in two patients who were both heterozygous for the ASL mutant p.E189G. Our results suggest that ASL transcripts can contribute to the highly variable phenotype in ASA patients if expressed at high levels. Especially, the exon 2-deleted ASL variant may form a heterotetramer with wild type or mutant ASL, causing markedly reduced ASL activity.

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.

Social cognition as a mediator variable between neurocognition and functional outcome in schizophrenia: empirical review and new results by structural equation modeling

Cognitive impairments are currently regarded as important determinants of functional domains and are promising treatment goals in schizophrenia. Nevertheless, the exact nature of the interdependent relationship between neurocognition and social cognition as well as the relative contribution of each of these factors to adequate functioning remains unclear. The purpose of this article is to systematically review the findings and methodology of studies that have investigated social cognition as a mediator variable between neurocognitive performance and functional outcome in schizophrenia. Moreover, we carried out a study to evaluate this mediation hypothesis by the means of structural equation modeling in a large sample of 148 schizophrenia patients. The review comprised 15 studies. All but one study provided evidence for the mediating role of social cognition both in cross-sectional and in longitudinal designs. Other variables like motivation and social competence additionally mediated the relationship between social cognition and functional outcome. The mean effect size of the indirect effect was 0.20. However, social cognitive domains were differentially effective mediators. On average, 25% of the variance in functional outcome could be explained in the mediation model. The results of our own statistical analysis are in line with these conclusions: Social cognition mediated a significant indirect relationship between neurocognition and functional outcome. These results suggest that research should focus on differential mediation pathways. Future studies should also consider the interaction with other prognostic factors, additional mediators, and moderators in order to increase the predictive power and to target those factors relevant for optimizing therapy effects.

The Relationship Between Sociosexuality and Aspects of Body Image in Men and Women: A Structural Equation Modeling Approach

The present study investigated the association between individual differences in sociosexual orientation and four aspects of body image in 156 male and 136 female students. While men were characterized by a less restricted sociosexual orientation, higher self-perceived physical attractiveness, and more pronounced self-rated physical assertiveness, women placed more emphasis on accentuation of body presentation. Structural equation modeling revealed significant positive relationships between sociosexual attitudes and physical attractiveness and accentuation of body presentation as well as between sociosexual behavior and physical attractiveness for the total sample. When introducing sex as a grouping variable, the attitudinal and behavioral components of sociosexuality were reliably related to both physical attractiveness and accentuation of body presentation as two aspects of body image in men, but not in women. Furthermore, our findings suggest that accentuation of body presentation represents a goal-directed behavior in men to increase the likelihood of having uncommitted sex but serves additional functions widely unrelated to unrestrictive sociosexual behavior in women.

Modeling Features at Runtime

A feature represents a functional requirement fulfilled by a system. Since many maintenance tasks are expressed in terms of features, it is important to establish the correspondence between a feature and its implementation in source code. Traditional approaches to establish this correspondence exercise features to generate a trace of runtime events, which is then processed by post-mortem analysis. These approaches typically generate large amounts of data to analyze. Due to their static nature, these approaches do not support incremental and interactive analysis of features. We propose a radically different approach called live feature analysis, which provides a model at runtime of features. Our approach analyzes features on a running system and also makes it possible to grow feature representations by exercising different scenarios of the same feature, and identifies execution elements even to the sub-method level. We describe how live feature analysis is implemented effectively by annotating structural representations of code based on abstract syntax trees. We illustrate our live analysis with a case study where we achieve a more complete feature representation by exercising and merging variants of feature behavior and demonstrate the efficiency or our technique with benchmarks.

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.

Structural model of the CopA copper ATPase of Enterococcus hirae based on chemical cross-linking

The CopA copper ATPase of Enterococcus hirae belongs to the family of heavy metal pumping CPx-type ATPases and shares 43% sequence similarity with the human Menkes and Wilson copper ATPases. Due to a lack of suitable protein crystals, only partial three-dimensional structures have so far been obtained for this family of ion pumps. We present a structural model of CopA derived by combining topological information obtained by intramolecular cross-linking with molecular modeling. Purified CopA was cross-linked with different bivalent reagents, followed by tryptic digestion and identification of cross-linked peptides by mass spectrometry. The structural proximity of tryptic fragments provided information about the structural arrangement of the hydrophilic protein domains, which was integrated into a three-dimensional model of CopA. Comparative modeling of CopA was guided by the sequence similarity to the calcium ATPase of the sarcoplasmic reticulum, Serca1, for which detailed structures are available. In addition, known partial structures of CPx-ATPase homologous to CopA were used as modeling templates. A docking approach was used to predict the orientation of the heavy metal binding domain of CopA relative to the core structure, which was verified by distance constraints derived from cross-links. The overall structural model of CopA resembles the Serca1 structure, but reveals distinctive features of CPx-type ATPases. A prominent feature is the positioning of the heavy metal binding domain. It features an orientation of the Cu binding ligands which is appropriate for the interaction with Cu-loaded metallochaperones in solution. Moreover, a novel model of the architecture of the intramembranous Cu binding sites could be derived.

Modeling of human P450 oxidoreductase structure by in silico mutagenesis and MD simulation

P450 oxidoreductase (POR) is the obligate electron donor for microsomal cytochrome P450s and mutations in POR cause several metabolic disorders. We have modeled the structure of human P450 oxidoreductase by in silico amino acid replacements in the rat POR crystal structure. The rat POR has 94% homology with human POR and 38 amino acids were replaced to make its sequence identical to human POR. Several rounds of molecular dynamic simulations refined the model and removed structural clashes from side chain alterations of replaced amino acids. This approach has the advantage of keeping the cofactor contacts and structural features of the core enzyme intact which could not be achieved by homology based approaches. The final model from our approach was of high quality and compared well with experimentally determined structures of other PORs. This model will be used for analyzing the structural implications of mutations and polymorphisms in human POR.

Measurement of structural anisotropy in femoral trabecular bone using clinical-resolution CT images

Discrepancies in finite-element model predictions of bone strength may be attributed to the simplified modeling of bone as an isotropic structure due to the resolution limitations of clinical-level Computed Tomography (CT) data. The aim of this study is to calculate the preferential orientations of bone (the principal directions) and the extent to which bone is deposited more in one direction compared to another (degree of anisotropy). Using 100 femoral trabecular samples, the principal directions and degree of anisotropy were calculated with a Gradient Structure Tensor (GST) and a Sobel Structure Tensor (SST) using clinical-level CT. The results were compared against those calculated with the gold standard Mean-Intercept-Length (MIL) fabric tensor using micro-CT. There was no significant difference between the GST and SST in the calculation of the main principal direction (median error=28°), and the error was inversely correlated to the degree of transverse isotropy (r=−0.34, p<0.01). The degree of anisotropy measured using the structure tensors was weakly correlated with the MIL-based measurements (r=0.2, p<0.001). Combining the principal directions with the degree of anisotropy resulted in a significant increase in the correlation of the tensor distributions (r=0.79, p<0.001). Both structure tensors were robust against simulated noise, kernel sizes, and bone volume fraction. We recommend the use of the GST because of its computational efficiency and ease of implementation. This methodology has the promise to predict the structural anisotropy of bone in areas with a high degree of anisotropy, and may improve the in vivo characterization of bone.

From transient to steady state deformation and grain size: A thermodynamic approach using elasto-visco-plastic numerical modeling

Numerical simulation experiments give insight into the evolving energy partitioning during high-strain torsion experiments of calcite. Our numerical experiments are designed to derive a generic macroscopic grain size sensitive flow law capable of describing the full evolution from the transient regime to steady state. The transient regime is crucial for understanding the importance of micro structural processes that may lead to strain localization phenomena in deforming materials. This is particularly important in geological and geodynamic applications where the phenomenon of strain localization happens outside the time frame that can be observed under controlled laboratory conditions. Ourmethod is based on an extension of the paleowattmeter approach to the transient regime. We add an empirical hardening law using the Ramberg-Osgood approximation and assess the experiments by an evolution test function of stored over dissipated energy (lambda factor). Parameter studies of, strain hardening, dislocation creep parameter, strain rates, temperature, and lambda factor as well asmesh sensitivity are presented to explore the sensitivity of the newly derived transient/steady state flow law. Our analysis can be seen as one of the first steps in a hybrid computational-laboratory-field modeling workflow. The analysis could be improved through independent verifications by thermographic analysis in physical laboratory experiments to independently assess lambda factor evolution under laboratory conditions.