Pathogenic LRRK2 mutations do not alter gene expression in cell model systems or human brain tissue

Point mutations in LRRK2 cause autosomal dominant Parkinson's disease. Despite extensive efforts to determine the mechanism of cell death in patients with LRRK2 mutations, the aetiology of LRRK2 PD is not well understood. To examine possible alterations in gene expression linked to the presence of LRRK2 mutations, we carried out a case versus control analysis of global gene expression in three systems: fibroblasts isolated from LRRK2 mutation carriers and healthy, non-mutation carrying controls; brain tissue from G2019S mutation carriers and controls; and HEK293 inducible LRRK2 wild type and mutant cell lines. No significant alteration in gene expression was found in these systems following correction for multiple testing. These data suggest that any alterations in basal gene expression in fibroblasts or cell lines containing mutations in LRRK2 are likely to be quantitatively small. This work suggests that LRRK2 is unlikely to play a direct role in modulation of gene expression, although it remains possible that this protein can influence mRNA expression under pathogenic cicumstances.

Rescue of gene expression by modified REST decoy oligonucleotides in a cellular model of Huntington's disease

Transcriptional dysfunction is a prominent hallmark of Huntington's disease (HD). Several transcription factors have been implicated in the aetiology of HD progression and one of the most prominent is repressor element 1 (RE1) silencing transcription factor (REST). REST is a global repressor of neuronal gene expression and in the presence of mutant Huntingtin increased nuclear REST levels lead to elevated RE1 occupancy and a concomitant increase in target gene repression, including brain-derived neurotrophic factor. It is of great interest to devise strategies to reverse transcriptional dysregulation caused by increased nuclear REST and determine the consequences in HD. Thus far, such strategies have involved RNAi or mutant REST constructs. Decoys are double-stranded oligodeoxynucleotides corresponding to the DNA-binding element of a transcription factor and act to sequester it, thereby abrogating its transcriptional activity. Here, we report the use of a novel decoy strategy to rescue REST target gene expression in a cellular model of HD. We show that delivery of the decoy in cells expressing mutant Huntingtin leads to its specific interaction with REST, a reduction in REST occupancy of RE1s and rescue of target gene expression, including Bdnf. These data point to an alternative strategy for rebalancing the transcriptional dysregulation in HD.

The double Caldeira-Leggett model: Derivation and solutions of the master equations, reservoir-induced interactions and decoherence

In this paper we analyze the double Caldeira-Leggett model: the path integral approach to two interacting dissipative harmonic oscillators. Assuming a general form of the interaction between the oscillators, we consider two different situations: (i) when each oscillator is coupled to its own reservoir, and (ii) when both oscillators are coupled to a common reservoir. After deriving and solving the master equation for each case, we analyze the decoherence process of particular entanglements in the positional space of both oscillators. To analyze the decoherence mechanism we have derived a general decay function, for the off-diagonal peaks of the density matrix, which applies both to common and separate reservoirs. We have also identified the expected interaction between the two dissipative oscillators induced by their common reservoir. Such a reservoir-induced interaction, which gives rise to interesting collective damping effects, such as the emergence of relaxation- and decoherence-free subspaces, is shown to be blurred by the high-temperature regime considered in this study. However, we find that different interactions between the dissipative oscillators, described by rotating or counter-rotating terms, result in different decay rates for the interference terms of the density matrix. (C) 2010 Elsevier B.V. All rights reserved.

Growth hormone 1 gene (GH1) polymorphisms as possible markers of the production potential of beef cattle using the Brazilian Canchim breed as a model

The growth hormone 1 gene (GH1) is a candidate gene for body weight and weight gain in cattle since it plays a fundamental role in growth regulation. We investigated the GH1 gene AluI and DdeI restriction enzyme polymorphisms, located 149 bp apart in the cattle genome, as possible markers of the production potential of Canchim crossbreed cattle, a 5/8 Charolais (Bos taurus) and 3/8 Nelore (Bos indicus) breed developed in Brazil, by evaluating the birth weight, weaning weight, yearling weight and plasma insulin-like growth factor-1 (IGF-1) concentration of 7 month to 10 months old Canchim calves (n = 204) of known genealogy and which had been genotyped for the AluI and DdeI markers. Our results showed significant effect (p < 0.05) between the homozygous DdeI+/DdeI+ polymorphism and the estimated breeding value for weaning weight (ESB-WW), while the AluI leucine homozygous (L/L) and leucine/valine (L/V) heterozygous polymorphisms showed no significant effect on the traits studied. The restriction sites of the two enzymes led to the formation of haplotypes which also exerted a significant effect (p < 0.05) on the ESB-WW, with the largest difference being 8.5 kg in favor of the homozygous L plus DdeI+/L plus DdeI+ genotype over the heterozygous L plus DdeI-/V plus DdeI+ genotype.

Spatial and temporal profiles for anti-inflammatory gene expression in leukocytes during a resolving model of peritonitis

The recent appreciation of the role played by endogenous counterregulatory mechanisms in controlling the outcome of the host inflammatory response requires specific analysis of their spatial and temporal profiles. In this study, we have focused on the glucocorticoid-regulated anti-inflammatory mediator annexin 1. Induction of peritonitis in wild-type mice rapidly (4 h) produced the expected signs of inflammation, including marked activation of resident cells (e.g., mast cells), migration of blood-borne leukocytes, mirrored by blood neutrophilia. These changes subsided after 48-96 h. In annexin 1null mice, the peritonitis response was exaggerated (∼40% at 4 h), with increased granulocyte migration and cytokine production. In blood leukocytes, annexin 1 gene expression was activated at 4, but not 24, h postzymosan, whereas protein levels were increased ai both time points. Locally, endothelial and mast cell annexin 1 gene expression was not detectable in basal conditions, whereas it was switched on during the inflammatory response. The significance of annexin 1 system plasticity in the anti-inflammatory properties of dexamethasone was assessed. Clear induction of annexin 1 gene in response to dexamethasone treatment was evident in the circulating and migrated leukocytes, and in connective tissue mast cells; this was associated with the steroid failure to inhibit leukocyte trafficking, cytokine synthesis, and mast cell degranulation in the annexin 1null mouse. In conclusion, understanding how inflammation is brought under control will help clarify the complex interplay between pro- and anti-inflammatory pathways operating during the host response to injury and infection. Copyright © 2006 by The American Association of Immunologists, Inc.

Experimental glaucoma model (ischemia and reperfusion): histology, morphometry, protein and gene espression of apoptosis pathway

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Endostatin gene therapy stimulates upregulation of ICAM-1 and VCAM-1 in a metastatic renal cell carcinoma model

One of the greatest challenges in urological oncology is renal cell carcinoma (RCC), which is the third leading cause of death in genitourinary cancers. RCCs are highly vascularized and respond positively to antiangiogenic therapy. Endostatin (ES) is a fragment of collagen XVIII that possesses antiangiogenic activity. In this study, we examined the potential of ES-based antiangiogenic therapy to activate tumor-associated endothelial cells in metastatic RCC (mRCC). Balb/c-bearing Renca cells were treated with NIH/3T3-LendSN or, as a control, with NIH/3T3-LXSN cells. The T-cell subsets and lymphocyte populations of tumors, mediastinal lymph nodes and the spleen were assessed by flow cytometry. The expression of intercellular adhesion molecule-1 (ICAM-1) and vascular cell adhesion molecule-1 (VCAM-1) was assessed by real-time PCR, flow cytometry and immunohistochemistry analysis. ES gene therapy led to an increase in the percentage of infiltrating CD4-interferon (IFN)-gamma cells (P<0.05), CD8-IFN-gamma cells (P<0.01) and CD49b-tumor necrosis factor-alpha cells (P<0.01). In addition, ES therapy caused an increase at the mRNA level of ICAM-1 (1.4-fold; P<0.01) and VCAM-1 (1.5-fold) (control vs treated group; P<0.001). Through flow cytometry, we found a significant increase in the CD34/ICAM-1 cells (8.1-fold; P<0.001) and CD34/VCAM-1 cells (1.6-fold; P<0.05). ES gene therapy induced a significant increase in both T CD4 and CD8 cells in the lymph nodes and the spleen, suggesting that ES therapy may facilitate cell survival or clonal expansion. CD49b cells were also present in increased quantities in all of these organs. In this study, we demonstrate an antitumor inflammatory effect of ES in an mRCC model, and this effect is mediated by an increase in ICAM-1 and VCAM-1 expression in tumor-associated endothelial cells.

Four-gene expression model predictive of lymph node metastases in oral squamous cell carcinoma

Background. Previous knowledge of cervical lymph node compromise may be crucial to choose the best treatment strategy in oral squamous cell carcinoma (OSCC). Here we propose a set four genes, whose mRNA expression in the primary tumor predicts nodal status in OSCC, excluding tongue. Material and methods. We identified differentially expressed genes in OSCC with and without compromised lymph nodes using Differential Display RT-PCR. Known genes were chosen to be validated by means of Northern blotting or real time RT-PCR (qRT-PCR). Thereafter we constructed a Nodal Index (NI) using discriminant analysis in a learning set of 35 patients, which was further validated in a second independent group of 20 patients. Results. Of the 63 differentially expressed known genes identified comparing three lymph node positive (pN+) and three negative (pN0) primary tumors, 23 were analyzed by Northern analysis or RT-PCR in 49 primary tumors. Six genes confirmed as differentially expressed were used to construct a NI, as the best set predictive of lymph nodal status, with the final result including four genes. The NI was able to correctly classify 32 of 35 patients comprising the learning group (88.6%; p = 0.009). Casein kinase 1alpha1 and scavenger receptor class B, member 2 were found to be up regulated in pN + group in contrast to small proline-rich protein 2B and Ras-GTPase activating protein SH3 domain-binding protein 2 which were upregulated in the pN0 group. We validated further our NI in an independent set of 20 primary tumors, 11 of them pN0 and nine pN+ with an accuracy of 80.0% (p = 0.012). Conclusions. The NI was an independent predictor of compromised lymph nodes, taking into the consideration tumor size and histological grade. The genes identified here that integrate our "Nodal Index" model are predictive of lymph node metastasis in OSCC.

Bayesian model accounting for within-class biological variability in Serial Analysis of Gene Expression (SAGE)

Abstract Background An important challenge for transcript counting methods such as Serial Analysis of Gene Expression (SAGE), "Digital Northern" or Massively Parallel Signature Sequencing (MPSS), is to carry out statistical analyses that account for the within-class variability, i.e., variability due to the intrinsic biological differences among sampled individuals of the same class, and not only variability due to technical sampling error. Results We introduce a Bayesian model that accounts for the within-class variability by means of mixture distribution. We show that the previously available approaches of aggregation in pools ("pseudo-libraries") and the Beta-Binomial model, are particular cases of the mixture model. We illustrate our method with a brain tumor vs. normal comparison using SAGE data from public databases. We show examples of tags regarded as differentially expressed with high significance if the within-class variability is ignored, but clearly not so significant if one accounts for it. Conclusion Using available information about biological replicates, one can transform a list of candidate transcripts showing differential expression to a more reliable one. Our method is freely available, under GPL/GNU copyleft, through a user friendly web-based on-line tool or as R language scripts at supplemental web-site.

Modeling gene expression regulatory networks with the sparse vector autoregressive model

Abstract Background To understand the molecular mechanisms underlying important biological processes, a detailed description of the gene products networks involved is required. In order to define and understand such molecular networks, some statistical methods are proposed in the literature to estimate gene regulatory networks from time-series microarray data. However, several problems still need to be overcome. Firstly, information flow need to be inferred, in addition to the correlation between genes. Secondly, we usually try to identify large networks from a large number of genes (parameters) originating from a smaller number of microarray experiments (samples). Due to this situation, which is rather frequent in Bioinformatics, it is difficult to perform statistical tests using methods that model large gene-gene networks. In addition, most of the models are based on dimension reduction using clustering techniques, therefore, the resulting network is not a gene-gene network but a module-module network. Here, we present the Sparse Vector Autoregressive model as a solution to these problems. Results We have applied the Sparse Vector Autoregressive model to estimate gene regulatory networks based on gene expression profiles obtained from time-series microarray experiments. Through extensive simulations, by applying the SVAR method to artificial regulatory networks, we show that SVAR can infer true positive edges even under conditions in which the number of samples is smaller than the number of genes. Moreover, it is possible to control for false positives, a significant advantage when compared to other methods described in the literature, which are based on ranks or score functions. By applying SVAR to actual HeLa cell cycle gene expression data, we were able to identify well known transcription factor targets. Conclusion The proposed SVAR method is able to model gene regulatory networks in frequent situations in which the number of samples is lower than the number of genes, making it possible to naturally infer partial Granger causalities without any a priori information. In addition, we present a statistical test to control the false discovery rate, which was not previously possible using other gene regulatory network models.

Study of apoptotic deletion mediated by Bifidobacterium longum with construction of recombinant strains for Serpin encoding gene and phenotypes comparison in a pig cell model

The first part of the research project of the Co-Advisorship Ph.D Thesis was aimed to select the best Bifidobacterium longum strains suitable to set the basis of our study. We were looking for strains with the abilities to colonize the intestinal mucosa and with good adhesion capacities, so that we can test these strains to investigate their ability to induce apoptosis in “damaged” intestinal cells. Adhesion and apoptosis are the two process that we want to study to better understand the role of an adhesion protein that we have previously identified and that have top scores homologies with the recent serpin encoding gene identified in B. longum by Nestlè researchers. Bifidobacterium longum is a probiotic, known for its beneficial effects to the human gut and even for its immunomodulatory and antitumor activities. Recently, many studies have stressed out the intimate relation between probiotic bacteria and the GIT mucosa and their influence on human cellular homeostasis. We focused on the apoptotic deletion of cancer cells induced by B. longum. This has been valued in vitro, performing the incubation of three B.longum strains with enterocyte-like Caco- 2 cells, to evidence DNA fragmentation, a cornerstone of apoptosis. The three strains tested were defined for their adhesion properties using adhesion and autoaggregation assays. These features are considered necessary to select a probiotic strain. The three strains named B12, B18 and B2990 resulted respectively: “strong adherent”, “adherent” and “non adherent”. Then, bacteria were incubated with Caco-2 cells to investigate apoptotic deletion. Cocultures of Caco-2 cells with B. longum resulted positive in DNA fragmentation test, only when adherent strains were used (B12 and B18). These results indicate that the interaction with adherent B. longum can induce apoptotic deletion of Caco-2 cells, suggesting a role in cellular homeostasis of the gastrointestinal tract and in restoring the ecology of damaged colon tissues. These results were used to keep on researching and the strains tested were used as recipient of recombinant techniques aimed to originate new B.longum strains with enhanced capacity of apoptotic induction in “damaged” intestinal cells. To achieve this new goal it was decided to clone the serpin encoding gene of B. longum, so that we can understand its role in adhesion and apoptosis induction. Bifidobacterium longum has immunostimulant activity that in vitro can lead to apoptotic response of Caco-2 cell line. It secretes a hypothetical eukaryotic type serpin protein, which could be involved in this kind of deletion of damaged cells. We had previously characterised a protein that has homologies with the hypothetical serpin of B. longum (DD087853). In order to create Bifidobacterium serpin transformants, a B. longum cosmid library was screened with a PCR protocol using specific primers for serpin gene. After fragment extraction, the insert named S1 was sub-cloned into pRM2, an Escherichia coli - Bifidobacterium shuttle vector, to construct pRM3. Several protocols for B. longum transformation were performed and the best efficiency was obtained using MRS medium and raffinose. Finally bacterial cell supernatants were tested in a dotblot assay to detect antigens presence against anti-antitrypsin polyclonal antibody. The best signal was produced by one starin that has been renamed B. longum BLKS 7. Our research study was aimed to generate transformants able to over express serpin encoding gene, so that we can have the tools for a further study on bacterial apoptotic induction of Caco-2 cell line. After that we have originated new trasformants the next step to do was to test transformants abilities when exposed to an intestinal cell model. In fact, this part of the project was achieved in the Department of Biochemistry of the Medical Faculty of the University of Maribor, guest of the abroad supervisor of the Co-Advisorship Doctoral Thesis: Prof. Avrelija Cencic. In this study we examined the probiotic ability of some bacterial strains using intestinal cells from a 6 years old pig. The use of intestinal mammalian cells is essential to study this symbiosis and a functional cell model mimics a polarised epithelium in which enterocytes are separated by tight junctions. In this list of strains we have included the Bifidobacterium longum BKS7 transformant strain that we have previously originated; in order to compare its abilities. B. longum B12 wild type and B. longum BKS7 transformant and eight Lactobacillus strains of different sources were co-cultured with porcine small intestine epithelial cells (PSI C1) and porcine blood monocytes (PoM2) in Transwell filter inserts. The strains, including Lb. gasseri, Lb. fermentum, Lb. reuterii, Lb. plantarum and unidentified Lactobacillus from kenyan maasai milk and tanzanian coffee, were assayed for activation of cell lines, measuring nitric oxide by Griess reaction, H202 by tetramethylbenzidine reaction and O2 - by cytochrome C reduction. Cytotoxic effect by crystal violet staining and induction on metabolic activity by MTT cell proliferation assay were tested too. Transepithelial electrical resistance (TER) of polarised PSI C1 was measured during 48 hours co-culture. TER, used to observe epithelium permeability, decrease during pathogenesis and tissue becomes permeable to ion passive flow lowering epithelial barrier function. Probiotics can prevent or restore increased permeability. Lastly, dot-blot was achieved against Interleukin-6 of treated cells supernatants. The metabolic activity of PoM2 and PSI C1 increased slightly after co-culture not affecting mitochondrial functions. No strain was cytotoxic over PSI C1 and PoM2 and no cell activation was observed, as measured by the release of NO2, H202 and O2 - by PoM2 and PSI C1. During coculture TER of polarised PSI C1 was two-fold higher comparing with constant TER (~3000 ) of untreated cells. TER raise generated by bacteria maintains a low permeability of the epithelium. During treatment Interleukin-6 was detected in cell supernatants at several time points, confirming immunostimulant activity. All results were obtained using Lactobacillus paracasei Shirota e Carnobacterium divergens as controls. In conclusion we can state that both the list of putative probiotic bacteria and our new transformant strain of B. longum are not harmful when exposed to intestinal cells and could be selected as probiotics, because can strengthen epithelial barrier function and stimulate nonspecific immunity of intestinal cells on a pig cell model. Indeed, we have found out that none of the strains tested that have good adhesion abilities presents citotoxicity to the intestinal cells and that non of the strains tested can induce cell lines to produce high level of ROS, neither NO2. Moreover we have assayed even the capacity of producing certain citokynes that are correlated with immune response. The detection of Interleukin-6 was assayed in all our samples, including B.longum transformant BKS 7 strain, this result indicates that these bacteria can induce a non specific immune response in the intestinal cells. In fact, when we assayed the presence of Interferon-gamma in cells supernatant after bacterial exposure, we have no positive signals, that means that there is no activation of a specific immune response, thus confirming that these bacteria are not recognize as pathogen by the intestinal cells and are certainly not harmful for intestinal cells. The most important result is the measure of Trans Epithelial Electric Resistance that have shown how the intestinal barrier function get strengthen when cells are exposed to bacteria, due to a reduction of the epithelium permeability. We have now a new strain of B. longum that will be used for further studies above the mechanism of apoptotic induction to “damaged cells” and above the process of “restoring ecology”. This strain will be the basis to originate new transformant strains for Serpin encoding gene that must have better performance and shall be used one day even in clinical cases as in “gene therapy” for cancer treatment and prevention.

Gene therapy for a novel mouse model of Canavan disease

Canavan disease (CD) is a rare leukodystrophy caused by loss-of-function mutations in the gene encoding aspartoacylase (ASPA), an oligodendrocyte-enriched enzyme. It is characterised by the accumulation of the ASPA substrate N-acetylaspartate (NAA) in brain, blood and urine, leading to a spongiform vacuolisation of the brain, severe motoric and cognitive impairments and premature death. To date, no therapy is available due to the lack of a gene-transfer system allowing transgene expression in oligodendrocytes (OLs) and the restoration of the missing enzyme. Hence, the aim of this study was to establish a novel gene-transfer system and its preclinical evaluation in a CD animal model.rnIn the first part of this thesis, a novel ASPA mouse mutant was generated. A βgeo cassette (including the genes encoding β-galactosidase and neomycin) flanked by frt sites was inserted into intron 1 of the intact aspa gene. Additionally, exon 2 was flanked by loxP sites for optional conditional deletion of the targeted locus. The resulting ASPA-deficient aspalacZ/lacZ-mouse was found to be an accurate model of CD and an important tool to identify novel aspects of its complex pathology. Homozygous mutants showed a CD-like histopathology, neurological impairment, behavioural deficits as well as a reduced body weight. Additionally, MRI data revealed changes in brain metabolite composition. rnRecombinant adeno-associated viral (rAAV) vectors have become a versatile tool for gene transfer to the central nervous system because they are efficient, non-toxic and replication-deficient. Based on the natural neurotropism of AAV vectors, AAV-based gene delivery has entered the clinics for the treatment of neurodegenerative diseases. However, the lack of AAV vectors with oligodendroglial tropism has precluded gene therapy for leukodystrophies. In the second part of this work, it was shown that the transduction profile of established AAV serotypes can be targeted towards OLs in a transcriptional approach, using the oligodendrocyte-specific myelin basic protein (MBP) promoter to drive transgene expression in OLs.rnIn the last part of this work, the therapeutic efficacy of AAV-mediated aspa gene transfer to OLs of juvenile aspalacZ/lacZ mice was evaluated. AAV-aspa injections into multiple sites of the brain parenchyma resulted in transduction of OLs in the grey and white matter throughout the brain. Histological abnormalities in the brain of ASPA-deficient mice were ameliorated and accompanied by a reduction of NAA levels. Furthermore, the treatment resulted in normalisation of body weight, motor function and nest-building behaviour. These data provide a proof-of-concept for a successful gene therapy of Canavan disease. This might pave the way towards translation into clinical application and serve as the basis for the genetic treatment of other leukodystrophies.

Stocking of captive-bred fish can cause long-term population decline and gene pool replacement: predictions from a population dynamics model incorporating density-dependent mortality

Releasing captive-bred fish into natural environments (stocking) is common in fisheries worldwide. Although stocking is believed to have a positive effect on fish abundance over the short term, little is known about the long-term consequences of recurrent stocking and its influence on natural populations. In fact, there are growing concerns that genetically maladapted captive-bred fish can eventually reduce the abundance of natural population. In this study, we develop a simple model to quantitatively investigate the condition under which recurrent stocking has long-term effects on the natural population. Using a population dynamics model that takes into account a density-dependent recruitment, a gene responsible for the fitness difference between wild and captive-bred fish, and hybridization between them, we show that there is little or no contribution of recurrent stocking to the stock enhancement without a replacement of the wild gene pool by the captive-bred gene pool. The model further predicted that stocking of an intermediate level causes a reduction, rather than enhancement, of population size over the long term. The population decline due to stocking was attributed to the fitness disadvantage of captive-bred fish and strong overcompensation at recruitment stage. These results suggest that it would be difficult to simultaneously attain population size recovery and conservation of the local gene pool when captive-bred fish have fitness disadvantage in the wild, although caution is needed when applying the predictions from the simplified model to a specific species or population.

Polymorphic variants of the multidrug resistance gene Mdr1a and response to antiepileptic drug treatment in the kindling model of epilepsy

Allelic variants of the human P-glycoprotein encoding gene MDR1 (ABCB1) are discussed to be associated with different clinical conditions including pharmacoresistance of epilepsy. However, conflicting data have been reported with regard to the functional relevance of MDR1 allelic variants for the response to antiepileptic drugs. To our knowledge, it is not known whether functionally relevant genetic polymorphisms also occur in the two genes (Mdr1a/Abcb1a, Mdr1b/Abcb1b) coding for P-glycoprotein in the brain of rodents. Therefore, we have started to search for polymorphisms in the Mdr1a gene, which governs the expression of P-glycoprotein in brain capillary endothelial cells in rats. In the kindling model of temporal lobe epilepsy, subgroups of phenytoin-sensitive and phenytoin-resistant rats were selected in repeated drug trials. Sequencing of the Mdr1a gene coding sequence in the subgroups revealed no general differences between drug-resistant and drug-sensitive rats of the Wistar outbred strain. A comparison between different inbred and outbred rat strains also gave no evidence for polymorphisms in the Mdr1a coding sequence. However, in exon-flanking intron sequences, four genetic variants were identified by comparison between these rats strains. In conclusion, the finding that Wistar rats vary in their response to phenytoin, while having the same genetic background, argues against a major impact of Mdr1a genetics on pharmacosensitivity to antiepileptic drugs in the amygdala kindling model.

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.