Dynamical principles in biological processes: A model of charge migration in proteins and DNA

The generalized master equations (GMEs) that contain multiple time scales have been derived quantum mechanically. The GME method has then been applied to a model of charge migration in proteins that invokes the hole hopping between local amino acid sites driven by the torsional motions of the floppy backbones. This model is then applied to analyze the experimental results for sequence-dependent long-range hole transport in DNA reported by Meggers et al. [Meggers, E., Michel-Beyerle, M. E., & Giese, B. (1998) J. Am. Chem. Soc. 120, 12950–12955]. The model has also been applied to analyze the experimental results of femtosecond dynamics of DNA-mediated electron transfer reported by Zewail and co-workers [Wan, C., Fiebig, T., Kelley, S. O., Treadway, C. R., Barton, J. K. & Zewail, A. H. (1999) Proc. Natl. Acad. Sci. USA 96, 6014–6019]. The initial events in the dynamics of protein folding have begun to attract attention. The GME obtained in this paper will be applicable to this problem.

Biological impact on mineral dissolution: Application of the lichen model to understanding mineral weathering in the rhizosphere

Microorganisms modify rates and mechanisms of chemical and physical weathering and clay growth, thus playing fundamental roles in soil and sediment formation. Because processes in soils are inherently complex and difficult to study, we employ a model based on the lichen–mineral system to identify the fundamental interactions. Fixed carbon released by the photosynthetic symbiont stimulates growth of fungi and other microorganisms. These microorganisms directly or indirectly induce mineral disaggregation, hydration, dissolution, and secondary mineral formation. Model polysaccharides were used to investigate direct mediation of mineral surface reactions by extracellular polymers. Polysaccharides can suppress or enhance rates of chemical weathering by up to three orders of magnitude, depending on the pH, mineral surface structure and composition, and organic functional groups. Mg, Mn, Fe, Al, and Si are redistributed into clays that strongly adsorb ions. Microbes contribute to dissolution of insoluble secondary phosphates, possibly via release of organic acids. These reactions significantly impact soil fertility. Below fungi–mineral interfaces, mineral surfaces are exposed to dissolved metabolic byproducts. Through this indirect process, microorganisms can accelerate mineral dissolution, leading to enhanced porosity and permeability and colonization by microbial communities.

Comparing in vitro, in situ, and in vivo experimental data in a three-dimensional model of mammalian cochlear mechanics

Normal mammalian hearing is refined by amplification of the motion of the cochlear partition. This partition, comprising the organ of Corti sandwiched between the basilar and tectorial membranes, contains the outer hair cells that are thought to drive this amplification process. Force generation by outer hair cells has been studied extensively in vitro and in situ, but, to understand cochlear amplification fully, it is necessary to characterize the role played by each of the components of the cochlear partition in vivo. Observations of cochlear partition motion in vivo are severely restricted by its inaccessibility and sensitivity to surgical trauma, so, for the present study, a computer model has been used to simulate the operation of the cochlea under different experimental conditions. In this model, which uniquely retains much of the three-dimensional complexity of the real cochlea, the motions of the basilar and tectorial membranes are fundamentally different during in situ- and in vivo-like conditions. Furthermore, enhanced outer hair cell force generation in vitro leads paradoxically to a decrease in the gain of the cochlear amplifier during sound stimulation to the model in vivo. These results suggest that it is not possible to extrapolate directly from experimental observations made in vitro and in situ to the normal operation of the intact organ in vivo.

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.

A model for the structure of the P domains in the subtilisin-like prohormone convertases

The proprotein convertases are a family of at least seven calcium-dependent endoproteases that process a wide variety of precursor proteins in the secretory pathway. All members of this family possess an N-terminal proregion, a subtilisin-like catalytic module, and an additional downstream well-conserved region of ≈150 amino acid residues, the P domain, which is not found in any other subtilase. The pro and catalytic domains cannot be expressed in the absence of the P domains; their thermodynamic instability may be attributable to the presence of large numbers of negatively charged Glu and Asp side chains in the substrate binding region for recognition of multibasic residue cleavage sites. Based on secondary structure predictions, we here propose that the P domains consist of 8-stranded β-barrels with well-organized inner hydrophobic cores, and therefore are independently folded components of the proprotein convertases. We hypothesize further that the P domains are integrated through strong hydrophobic interactions with the catalytic domains, conferring structural stability and regulating the properties and activity of the convertases. A molecular model of these interdomain interactions is proposed in this report.

A yeast model for the study of Batten disease

Although the CLN3 gene for Batten disease, the most common inherited neurovisceral storage disease of childhood, was identified in 1995, the function of the corresponding protein still remains elusive. We previously cloned the Saccharomyces cerevisiae homologue to the human CLN3 gene, designated BTN1, which is not essential and whose product is 39% identical and 59% similar to Cln3p. We report that btn1-Δ deletion yeast strains are more resistant to d-(−)-threo-2-amino-1-[p-nitrophenyl]-1,3-propanediol (denoted ANP), a phenotype that is complemented in yeast by the human CLN3 gene. Furthermore, the severity of Batten disease in humans and the degree of ANP resistance in yeast are related when the equivalent amino acid replacements in Cln3p and Btn1p are compared. These results indicate that yeast can be used as a model for the study of Batten disease.

Exploring the origins of topological frustration: Design of a minimally frustrated model of fragment B of protein A

Topological frustration in an energetically unfrustrated off-lattice model of the helical protein fragment B of protein A from Staphylococcus aureus was investigated. This Gō-type model exhibited thermodynamic and kinetic signatures of a well-designed two-state folder with concurrent collapse and folding transitions and single exponential kinetics at the transition temperature. Topological frustration is determined in the absence of energetic frustration by the distribution of Fersht φ values. Topologically unfrustrated systems present a unimodal distribution sharply peaked at intermediate φ, whereas highly frustrated systems display a bimodal distribution peaked at low and high φ values. The distribution of φ values in protein A was determined both thermodynamically and kinetically. Both methods yielded a unimodal distribution centered at φ = 0.3 with tails extending to low and high φ values, indicating the presence of a small amount of topological frustration. The contacts with high φ values were located in the turn regions between helices I and II and II and III, intimating that these hairpins are in large part required in the transition state. Our results are in good agreement with all-atom simulations of protein A, as well as lattice simulations of a three- letter code 27-mer (which can be compared with a 60-residue helical protein). The relatively broad unimodal distribution of φ values obtained from the all-atom simulations and that from the minimalist model for the same native fold suggest that the structure of the transition state ensemble is determined mostly by the protein topology and not energetic frustration.

Dose-response resistance to HIV-1/MuLV pseudotype virus ex vivo in a hairpin ribozyme transgenic mouse model

We have investigated the efficacy of a hairpin ribozyme targeting the 5′ leader sequence of HIV-1 RNA in a transgenic model system. Primary spleen cells derived from transgenic or control mice were infected with HIV-1/MuLV pseudotype virus. A significantly reduced susceptibility to infection in ribozyme-expressing transgenic spleen cells (P = 0.01) was shown. Variation of transgene-expression levels between littermates revealed a dose response between ribozyme expression and viral resistance, with an estimated cut off value below 0.2 copies of hairpin ribozyme per cell. These findings open up possibilities for studies on ribozyme efficacy and anti-HIV-1 gene therapy.

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.

Vocal imitation in zebra finches is inversely related to model abundance

A juvenile male zebra finch, Taeniopygia guttata, kept singly with its father develops a fairly complete imitation of the father’s song. The imitation is less complete when other male siblings are present, possibly because as imitation commences, model abundance increases. Here we examine the consequences of allowing more or less access to a song model. Young males heard a brief song playback when they pecked at a key, but different males were allowed to hear different numbers of playbacks per day. Using an automated procedure that scored the similarity between model and pupil songs, we discovered that 40 playbacks of the song motif per day, lasting a total of 30 sec, resulted in a fairly complete imitation. More exposure led to less complete imitation. Vocal imitation often may reflect the interaction of diverse influences. Among these, we should now include the possible inhibitory effect of model overabundance, which may foster individual identity and explain the vocal diversity found in zebra finches and other songbirds.

A mouse model for X-linked adrenoleukodystrophy

X-linked adrenoleukodystrophy (X-ALD) is a peroxisomal disorder with impaired β-oxidation of very long chain fatty acids (VLCFAs) and reduced function of peroxisomal very long chain fatty acyl-CoA synthetase (VLCS) that leads to severe and progressive neurological disability. The X-ALD gene, identified by positional cloning, encodes a peroxisomal membrane protein (adrenoleukodystrophy protein; ALDP) that belongs to the ATP binding cassette transporter protein superfamily. Mutational analyses and functional studies of the X-ALD gene confirm that it and not VLCS is the gene responsible for X-ALD. Its role in the β-oxidation of VLCFAs and its effect on the function of VLCS are unclear. The complex pathology of X-ALD and the extreme variability of its clinical phenotypes are also unexplained. To facilitate understanding of X-ALD pathophysiology, we developed an X-ALD mouse model by gene targeting. The X-ALD mouse exhibits reduced β-oxidation of VLCFAs, resulting in significantly elevated levels of saturated VLCFAs in total lipids from all tissues measured and in cholesterol esters from adrenal glands. Lipid cleft inclusions were observed in adrenocortical cells of X-ALD mice under the electron microscope. No neurological involvement has been detected in X-ALD mice up to 6 months. We conclude that X-ALD mice exhibit biochemical defects equivalent to those found in human X-ALD and thus provide an experimental system for testing therapeutic intervention.

Linking topography of its potential surface with the dynamics of folding of a protein model

The “3-color, 46-bead” model of a folding polypeptide is the vehicle for adapting to proteins a mode of analysis used heretofore for atomic clusters, to relate the topography of the potential surface to the dynamics that lead to formation of selected structures. The analysis is based on sequences of stationary points—successive minima, joined by saddles—that rise monotonically in energy from basin bottoms. Like structure-seeking clusters, the potential surface of the model studied here is staircase-like, rather than sawtooth-like, with highly collective motions required for passage from one minimum to the next. The surface has several deep basins whose minima correspond to very similar structures, but which are separated by high energy barriers.

Targeted ablation of the vitamin D receptor: An animal model of vitamin D-dependent rickets type II with alopecia

Vitamin D, the major steroid hormone that controls mineral ion homeostasis, exerts its actions through the vitamin D receptor (VDR). The VDR is expressed in many tissues, including several tissues not thought to play a role in mineral metabolism. Studies in kindreds with VDR mutations (vitamin D-dependent rickets type II, VDDR II) have demonstrated hypocalcemia, hyperparathyroidism, rickets, and osteomalacia. Alopecia, which is not a feature of vitamin D deficiency, is seen in some kindreds. We have generated a mouse model of VDDR II by targeted ablation of the second zinc finger of the VDR DNA-binding domain. Despite known expression of the VDR in fetal life, homozygous mice are phenotypically normal at birth and demonstrate normal survival at least until 6 months. They become hypocalcemic at 21 days of age, at which time their parathyroid hormone (PTH) levels begin to rise. Hyperparathyroidism is accompanied by an increase in the size of the parathyroid gland as well as an increase in PTH mRNA levels. Rickets and osteomalacia are seen by day 35; however, as early as day 15, there is an expansion in the zone of hypertrophic chondrocytes in the growth plate. In contrast to animals made vitamin D deficient by dietary means, and like some patients with VDDR II, these mice develop progressive alopecia from the age of 4 weeks.

A model of ocular dominance column development by competition for trophic factor

Recent experimental evidence has shown that application of certain neurotrophic factors (NTs) to the developing primary visual cortex prevents the development of ocular dominance (OD) columns. One interpretation of this result is that afferents from the lateral geniculate nucleus compete for postsynaptic trophic factor in an activity-dependent manner. Application of excess trophic factor eliminates this competition, thereby preventing OD column formation. We present a model of OD column development, incorporating Hebbian synaptic modification and activity-driven competition for NT, which accounts for both normal OD column development as well as the prevention of that development when competition is removed. In the “control” situation, when available NT is below a critical amount, OD columns form normally. These columns form without weight normalization procedures and in the presence of positive inter-eye correlations. In the “experimental” case, OD column development is prevented in a local neighborhood in which excess NT has been added. Our model proposes a biologically plausible mechanism for competition between neural populations that is motivated by several pieces of experimental data, thereby accounting for both normal and experimentally perturbed conditions.

A functionally defined model for the M2 proton channel of influenza A virus suggests a mechanism for its ion selectivity

The M2 protein from influenza A virus forms proton-selective channels that are essential to viral function and are the target of the drug amantadine. Cys scanning was used to generate a series of mutants with successive substitutions in the transmembrane segment of the protein, and the mutants were expressed in Xenopus laevis oocytes. The effect of the mutations on reversal potential, ion currents, and amantadine resistance were measured. Fourier analysis revealed a periodicity consistent with a four-stranded coiled coil or helical bundle. A three-dimensional model of this structure suggests a possible mechanism for the proton selectivity of the M2 channel of influenza virus.