A stochastic approximation algorithm with Markov chain Monte-Carlo method for incomplete data estimation problems

We propose a general procedure for solving incomplete data estimation problems. The procedure can be used to find the maximum likelihood estimate or to solve estimating equations in difficult cases such as estimation with the censored or truncated regression model, the nonlinear structural measurement error model, and the random effects model. The procedure is based on the general principle of stochastic approximation and the Markov chain Monte-Carlo method. Applying the theory on adaptive algorithms, we derive conditions under which the proposed procedure converges. Simulation studies also indicate that the proposed procedure consistently converges to the maximum likelihood estimate for the structural measurement error logistic regression model.

Mouse model of Sanfilippo syndrome type B produced by targeted disruption of the gene encoding α-N-acetylglucosaminidase

The Sanfilippo syndrome type B is an autosomal recessive disorder caused by mutation in the gene (NAGLU) encoding α-N-acetylglucosaminidase, a lysosomal enzyme required for the stepwise degradation of heparan sulfate. The most serious manifestations are profound mental retardation, intractable behavior problems, and death in the second decade. To generate a model for studies of pathophysiology and of potential therapy, we disrupted exon 6 of Naglu, the homologous mouse gene. Naglu−/− mice were healthy and fertile while young and could survive for 8–12 mo. They were totally deficient in α-N-acetylglucosaminidase and had massive accumulation of heparan sulfate in liver and kidney as well as secondary changes in activity of several other lysosomal enzymes in liver and brain and elevation of gangliosides GM2 and GM3 in brain. Vacuolation was seen in many cells, including macrophages, epithelial cells, and neurons, and became more prominent with age. Although most vacuoles contained finely granular material characteristic of glycosaminoglycan accumulation, large pleiomorphic inclusions were seen in some neurons and pericytes in the brain. Abnormal hypoactive behavior was manifested by 4.5-mo-old Naglu−/− mice in an open field test; the hyperactivity that is characteristic of affected children was not observed even in younger mice. In a Pavlovian fear conditioning test, the 4.5-mo-old mutant mice showed normal response to context, indicating intact hippocampal-dependent learning, but reduced response to a conditioning tone, perhaps attributable to hearing impairment. The phenotype of the α-N-acetylglucosaminidase-deficient mice is sufficiently similar to that of patients with the Sanfilippo syndrome type B to make these mice a good model for study of pathophysiology and for development of therapy.

Midbrain injection of recombinant adeno-associated virus encoding rat glial cell line-derived neurotrophic factor protects nigral neurons in a progressive 6-hydroxydopamine-induced degeneration model of Parkinson’s disease in rats

A recombinant adeno-associated virus (rAAV) vector capable of infecting cells and expressing rat glial cell line-derived neurotrophic factor (rGDNF), a putative central nervous system dopaminergic survival factor, under the control of a potent cytomegalovirus (CMV) immediate/early promoter (AAV-MD-rGDNF) was constructed. Two experiments were performed to evaluate the time course of expression of rAAV-mediated GDNF protein expression and to test the vector in an animal model of Parkinson’s disease. To evaluate the ability of rAAV-rGDNF to protect nigral dopaminergic neurons in the progressive Sauer and Oertel 6-hydroxydopamine (6-OHDA) lesion model, rats received perinigral injections of either rAAV-rGDNF virus or rAAV-lacZ control virus 3 weeks prior to a striatal 6-OHDA lesion and were sacrificed 4 weeks after 6-OHDA. Cell counts of back-labeled fluorogold-positive neurons in the substantia nigra revealed that rAAV-MD-rGDNF protected a significant number of cells when compared with cell counts of rAAV-CMV-lacZ-injected rats (94% vs. 51%, respectively). In close agreement, 85% of tyrosine hydroxylase-positive cells remained in the nigral rAAV-MD-rGDNF group vs. only 49% in the lacZ group. A separate group of rats were given identical perinigral virus injections and were sacrificed at 3 and 10 weeks after surgery. Nigral GDNF protein expression remained relatively stable over the 10 weeks investigated. These data indicate that the use of rAAV, a noncytopathic viral vector, can promote delivery of functional levels of GDNF in a degenerative model of Parkinson’s disease.

Hepatocyte gene therapy in a large animal: A neonatal bovine model of citrullinemia

The development of gene-replacement therapy for inborn errors of metabolism has been hindered by the limited number of suitable large-animal models of these diseases and by inadequate methods of assessing the efficacy of treatment. Such methods should provide sensitive detection of expression in vivo and should be unaffected by concurrent pharmacologic and dietary regimens. We present the results of studies in a neonatal bovine model of citrullinemia, an inborn error of urea-cycle metabolism characterized by deficiency of argininosuccinate synthetase and consequent life-threatening hyperammonemia. Measurements of the flux of nitrogen from orally administered 15NH4 to [15N]urea were used to determine urea-cycle activity in vivo. In control animals, these isotopic measurements proved to be unaffected by pharmacologic treatments. Systemic administration of a first-generation E1-deleted adenoviral vector expressing human argininosuccinate synthetase resulted in transduction of hepatocytes and partial correction of the enzyme defect. The isotopic method showed significant restoration of urea synthesis. Moreover, the calves showed clinical improvement and normalization of plasma glutamine levels after treatment. The results show the clinical efficacy of treating a large-animal model of an inborn error of hepatocyte metabolism in conjunction with a method for sensitively measuring correction in vivo. These studies will be applicable to human trials of the treatment of this disorder and other related urea-cycle disorders.

Intrastriatal injection of an adenoviral vector expressing glial-cell-line-derived neurotrophic factor prevents dopaminergic neuron degeneration and behavioral impairment in a rat model of Parkinson disease

Glial-cell-line-derived neurotrophic factor (GDNF) is a potent neurotrophic factor for adult nigral dopamine neurons in vivo. GDNF has both protective and restorative effects on the nigro-striatal dopaminergic (DA) system in animal models of Parkinson disease. Appropriate administration of this factor is essential for the success of its clinical application. Since it cannot cross the blood–brain barrier, a gene transfer method may be appropriate for delivery of the trophic factor to DA cells. We have constructed a recombinant adenovirus (Ad) encoding GDNF and injected it into rat striatum to make use of its ability to infect neurons and to be retrogradely transported by DA neurons. Ad-GDNF was found to drive production of large amounts of GDNF, as quantified by ELISA. The GDNF produced after gene transfer was biologically active: it increased the survival and differentiation of DA neurons in vitro. To test the efficacy of the Ad-mediated GDNF gene transfer in vivo, we used a progressive lesion model of Parkinson disease. Rats received injections unilaterally into their striatum first of Ad and then 6 days later of 6-hydroxydopamine. We found that mesencephalic nigral dopamine neurons of animals treated with the Ad-GDNF were protected, whereas those of animals treated with the Ad-β-galactosidase were not. This protection was associated with a difference in motor function: amphetamine-induced turning was much lower in animals that received the Ad-GDNF than in the animals that received Ad-β-galactosidase. This finding may have implications for the development of a treatment for Parkinson disease based on the use of neurotrophic factors.

Genetic definition of a protein-splicing domain: Functional mini-inteins support structure predictions and a model for intein evolution

Inteins are protein-splicing elements, most of which contain conserved sequence blocks that define a family of homing endonucleases. Like group I introns that encode such endonucleases, inteins are mobile genetic elements. Recent crystallography and computer modeling studies suggest that inteins consist of two structural domains that correspond to the endonuclease and the protein-splicing elements. To determine whether the bipartite structure of inteins is mirrored by the functional independence of the protein-splicing domain, the entire endonuclease component was deleted from the Mycobacterium tuberculosis recA intein. Guided by computer modeling studies, and taking advantage of genetic systems designed to monitor intein function, the 440-aa Mtu recA intein was reduced to a functional mini-intein of 137 aa. The accuracy of splicing of several mini-inteins was verified. This work not only substantiates structure predictions for intein function but also supports the hypothesis that, like group I introns, mobile inteins arose by an endonuclease gene invading a sequence encoding a small, functional splicing element.

Targeted disruption of the mouse gene encoding steroidogenic acute regulatory protein provides insights into congenital lipoid adrenal hyperplasia

An essential component of regulated steroidogenesis is the translocation of cholesterol from the cytoplasm to the inner mitochondrial membrane where the cholesterol side-chain cleavage enzyme carries out the first committed step in steroidogenesis. Recent studies showed that a 30-kDa mitochondrial phosphoprotein, designated steroidogenic acute regulatory protein (StAR), is essential for this translocation. To allow us to explore the roles of StAR in a system amenable to experimental manipulation and to develop an animal model for the human disorder lipoid congenital adrenal hyperplasia (lipoid CAH), we used targeted gene disruption to produce StAR knockout mice. These StAR knockout mice were indistinguishable initially from wild-type littermates, except that males and females had female external genitalia. After birth, they failed to grow normally and died from adrenocortical insufficiency. Hormone assays confirmed severe defects in adrenal steroids—with loss of negative feedback regulation at hypothalamic–pituitary levels—whereas hormones constituting the gonadal axis did not differ significantly from levels in wild-type littermates. Histologically, the adrenal cortex of StAR knockout mice contained florid lipid deposits, with lesser deposits in the steroidogenic compartment of the testis and none in the ovary. The sex-specific differences in gonadal involvement support a two-stage model of the pathogenesis of StAR deficiency, with trophic hormone stimulation inducing progressive accumulation of lipids within the steroidogenic cells and ultimately causing their death. These StAR knockout mice provide a useful model system in which to determine the mechanisms of StAR’s essential roles in adrenocortical and gonadal steroidogenesis.

Mycoparasitic interaction relieves binding of the Cre1 carbon catabolite repressor protein to promoter sequences of the ech42 (endochitinase-encoding) gene in Trichoderma harzianum

The fungus Trichoderma harzianum is a potent mycoparasite of various plant pathogenic fungi. We have studied the molecular regulation of mycoparasitism in the host/mycoparasite system Botrytis cinerea/T. harzianum. Protein extracts, prepared from various stages of mycoparasitism, were used in electrophoretic mobility-shift assays (EMSAs) with two promoter fragments of the ech-42 (42-kDa endochitinase-encoding) gene of T. harzianum. This gene was chosen as a model because its expression is triggered during mycoparasitic interaction [Carsolio, C., Gutierrez, A., Jimenez, B., van Montagu, M. & Herrera-Estrella, A. (1994) Proc. Natl. Acad. Sci. USA 91, 10903–10907]. All cell-free extracts formed high-molecular weight protein–DNA complexes, but those obtained from mycelia activated for mycoparasitic attack formed a complex with greater mobility. Competition experiments, using oligonucleotides containing functional and nonfunctional consensus sites for binding of the carbon catabolite repressor Cre1, provided evidence that the complex from nonmycoparasitic mycelia involves the binding of Cre1 to both fragments of the ech-42 promoter. The presence of two and three consensus sites for binding of Cre1 in the two ech-42 promoter fragments used is consistent with these findings. In contrast, the formation of the protein–DNA complex from mycoparasitic mycelia is unaffected by the addition of the competing oligonucleotides and hence does not involve Cre1. Addition of equal amounts of protein of cell-free extracts from nonmycoparasitic mycelia converted the mycoparasitic DNA–protein complex into the nonmycoparasitic complex. The addition of the purified Cre1::glutathione S-transferase protein to mycoparasitic cell-free extracts produced the same effect. These findings suggest that ech-42 expression in T. harzianum is regulated by (i) binding of Cre1 to two single sites in the ech-42 promoter, (ii) binding of a “mycoparasitic” protein–protein complex to the ech-42 promoter in vicinity of the Cre1 binding sites, and (iii) functional inactivation of Cre1 upon mycoparasitic interaction to enable the formation of the mycoparasitic protein–DNA complex.

Prophylactic DNA vaccine for hepatitis C virus (HCV) infection: HCV-specific cytotoxic T lymphocyte induction and protection from HCV-recombinant vaccinia infection in an HLA-A2.1 transgenic mouse model

DNA vaccines express antigens intracellularly and effectively induce cellular immune responses. Because only chimpanzees can be used to model human hepatitis C virus (HCV) infections, we developed a small-animal model using HLA-A2.1-transgenic mice to test induction of HLA-A2.1-restricted cytotoxic T lymphocytes (CTLs) and protection against recombinant vaccinia expressing HCV-core. A plasmid encoding the HCV-core antigen induced CD8+ CTLs specific for three conserved endogenously expressed core peptides presented by human HLA-A2.1. When challenged, DNA-immunized mice showed a substantial (5–12 log10) reduction in vaccinia virus titer compared with mock-immunized controls. This protection, lasting at least 14 mo, was shown to be mediated by CD8+ cells. Thus, a DNA vaccine expressing HCV-core is a potential candidate for a prophylactic vaccine for HLA-A2.1+ humans.

Neonatal gene transfer leads to widespread correction of pathology in a murine model of lysosomal storage disease

For many inborn errors of metabolism, early treatment is critical to prevent long-term developmental sequelae. We have used a gene-therapy approach to demonstrate this concept in a murine model of mucopolysaccharidosis type VII (MPS VII). Newborn MPS VII mice received a single intravenous injection with 5.4 × 106 infectious units of recombinant adeno-associated virus encoding the human β-glucuronidase (GUSB) cDNA. Therapeutic levels of GUSB expression were achieved by 1 week of age in liver, heart, lung, spleen, kidney, brain, and retina. GUSB expression persisted in most organs for the 16-week duration of the study at levels sufficient to either reduce or prevent completely lysosomal storage. Of particular significance, neurons, microglia, and meninges of the central nervous system were virtually cleared of disease. In addition, neonatal treatment of MPS VII mice provided access to the central nervous system via an intravenous route, avoiding a more invasive procedure later in life. These data suggest that gene transfer mediated by adeno-associated virus can achieve therapeutically relevant levels of enzyme very early in life and that the rapid growth and differentiation of tissues does not limit long-term expression.

Rapid restoration of visual pigment and function with oral retinoid in a mouse model of childhood blindness

Mutations in the retinal pigment epithelium gene encoding RPE65 are a cause of the incurable early-onset recessive human retinal degenerations known as Leber congenital amaurosis. Rpe65-deficient mice, a model of Leber congenital amaurosis, have no rod photopigment and severely impaired rod physiology. We analyzed retinoid flow in this model and then intervened by using oral 9-cis-retinal, attempting to bypass the biochemical block caused by the genetic abnormality. Within 48 h, there was formation of rod photopigment and dramatic improvement in rod physiology, thus demonstrating that mechanism-based pharmacological intervention has the potential to restore vision in otherwise incurable genetic retinal degenerations.

Mouse strain differences determine severity of iron accumulation in Hfe knockout model of hereditary hemochromatosis

Hereditary hemochromatosis (HH) is a common disorder of iron metabolism caused by mutation in HFE, a gene encoding an MHC class I-like protein. Clinical studies demonstrate that the severity of iron loading is highly variable among individuals with identical HFE genotypes. To determine whether genetic factors other than Hfe genotype influence the severity of iron loading in the murine model of HH, we bred the disrupted murine Hfe allele onto three different genetically defined mouse strains (AKR, C57BL/6, and C3H), which differ in basal iron status and sensitivity to dietary iron loading. Serum transferrin saturations (percent saturation of serum transferrin with iron), hepatic and splenic iron concentrations, and hepatocellular iron distribution patterns were compared for wild-type (Hfe +/+), heterozygote (Hfe +/−), and knockout (Hfe −/−) mice from each strain. Although the Hfe −/− mice from all three strains demonstrated increased transferrin saturations and liver iron concentrations compared with Hfe +/+ mice, strain differences in severity of iron accumulation were striking. Targeted disruption of the Hfe gene led to hepatic iron levels in Hfe −/− AKR mice that were 2.5 or 3.6 times higher than those of Hfe −/− C3H or Hfe −/− C57BL/6 mice, respectively. The Hfe −/− mice also demonstrated strain-dependent differences in transferrin saturation, with the highest values in AKR mice and the lowest values in C3H mice. These observations demonstrate that heritable factors markedly influence iron homeostasis in response to Hfe disruption. Analysis of mice from crosses between C57BL/6 and AKR mice should allow the mapping and subsequent identification of genes modifying the severity of iron loading in this murine model of HH.

RNA virus error catastrophe: Direct molecular test by using ribavirin

RNA viruses evolve rapidly. One source of this ability to rapidly change is the apparently high mutation frequency in RNA virus populations. A high mutation frequency is a central tenet of the quasispecies theory. A corollary of the quasispecies theory postulates that, given their high mutation frequency, animal RNA viruses may be susceptible to error catastrophe, where they undergo a sharp drop in viability after a modest increase in mutation frequency. We recently showed that the important broad-spectrum antiviral drug ribavirin (currently used to treat hepatitis C virus infections, among others) is an RNA virus mutagen, and we proposed that ribavirin's antiviral effect is by forcing RNA viruses into error catastrophe. However, a direct demonstration of error catastrophe has not been made for ribavirin or any RNA virus mutagen. Here we describe a direct demonstration of error catastrophe by using ribavirin as the mutagen and poliovirus as a model RNA virus. We demonstrate that ribavirin's antiviral activity is exerted directly through lethal mutagenesis of the viral genetic material. A 99.3% loss in viral genome infectivity is observed after a single round of virus infection in ribavirin concentrations sufficient to cause a 9.7-fold increase in mutagenesis. Compiling data on both the mutation levels and the specific infectivities of poliovirus genomes produced in the presence of ribavirin, we have constructed a graph of error catastrophe showing that normal poliovirus indeed exists at the edge of viability. These data suggest that RNA virus mutagens may represent a promising new class of antiviral drugs.

Long-lasting reduction of inhibitory function and gamma-aminobutyric acid type A receptor subunit mRNA expression in a model of temporal lobe epilepsy.

This study evaluated hippocampal inhibitory function and the level of expression of gamma-aminobutyric acid type A (GABAA) receptor mRNA in an in vivo model of epilepsy. Chronic recurrent limbic seizures were induced in rats using injections of pilocarpine. Electrophysiological studies performed on hippocampal slices prepared from control and epileptic animals 1 to 2 months after pilocarpine injections demonstrated a significant hyperexcitability in the epileptic animals. Reduced levels of mRNA expression for the alpha 2 and alpha 5 subunits of the GABAA receptors were evident in the CA1, CA2, and CA3 regions of the hippocampus of epileptic animals. No decrease in mRNA encoding alpha 1, beta 2, or gamma 2 GABAA receptor subunits was observed. In addition, no change in the mRNA levels of alpha CaM kinase II was seen. Selective decreases in mRNA expression did not correlate with neuronal cell loss. The results indicate that selective, long-lasting reduction of GABAA subunit mRNA expression and increased excitability, possibly reflecting loss of GABAergic inhibition, occur in an in vivo model of partial complex epilepsy.

Zinc- and iron-rubredoxins from Clostridium pasteurianum at atomic resolution: a high-precision model of a ZnS4 coordination unit in a protein.

The Zn(Scys)4 unit is present in numerous proteins, where it assumes structural, regulatory, or catalytic roles. The same coordination is found naturally around iron in rubredoxins, several structures of which have been refined at resolutions of, or near to, 1 A. The fold of the small protein rubredoxin around its metal ion is an excellent model for many zinc finger proteins. Zn-substituted rubredoxin and its Fe-containing counterpart were both obtained as the products of the expression in Escherichia coli of the rubredoxin-encoding gene from Clostridium pasteurianum. The structures of both proteins have been refined with an anisotropic model at atomic resolution (1.1 A, R = 8.3% for Fe-rubredoxin, and 1.2 A, R = 9.6% for Zn-rubredoxin) and are very similar. The most significant differences are increased lengths of the M-S bonds in Zn-rubredoxin (average length, 2.345 A) as compared with Fe-rubredoxin (average length, 2.262 A). An increase of the CA-CB-SG-M dihedral angles involving Cys-6 and Cys-39, the first cysteines of each of the Cys-Xaa-Xaa-Cys metal binding motifs, has been observed. Another consequence of the replacement of iron by zinc is that the region around residues 36-46 undergoes larger displacements than the remainder of the polypeptide chain. Despite these changes, the main features of the FeS4 site, namely a local 2-fold symmetry and the characteristic network of N-H...S hydrogen bonds, are conserved in the ZnS4 site. The Zn-substituted rubredoxin provides the first precise structure of a Zn(Scys)4 unit in a protein. The nearly identical fold of rubredoxin around iron or zinc suggests that at least in some of the sites where the metal has mainly a structural role-e.g., zinc fingers-the choice of the relevant metal may be directed by its cellular availability and mobilization processes rather than by its chemical nature.