A feedback-controlled ensemble model of the stress-responsive hypothalamo-pituitary-adrenal axis

The present work develops and implements a biomathematical statement of how reciprocal connectivity drives stress-adaptive homeostasis in the corticotropic (hypothalamo-pituitary-adrenal) axis. In initial analyses with this interactive construct, we test six specific a priori hypotheses of mechanisms linking circadian (24-h) rhythmicity to pulsatile secretory output. This formulation offers a dynamic framework for later statistical estimation of unobserved in vivo neurohormone secretion and within-axis, dose-responsive interfaces in health and disease. Explication of the core dynamics of the stress-responsive corticotropic axis based on secure physiological precepts should help to unveil new biomedical hypotheses of stressor-specific system failure.

Three-dimensional model of the potyviral genome-linked protein.

The full sequence of the genome-linked viral protein (VPg) cistron located in the central part of potato virus Y (common strain) genome has been identified. The VPg gene codes for a protein of 188 amino acids, with significant homology to other known potyviral VPg polypeptides. A three-dimensional model structure of VPg is proposed on the basis of similarity of hydrophobic-hydrophilic residue distribution to the sequence of malate dehydrogenase of known crystal structure. The 5' end of the viral RNA can be fitted to interact with the protein through the exposed hydroxyl group of Tyr-64, in agreement with experimental data. The complex favors stereochemically the formation of a phosphodiester bond [5'-(O4-tyrosylphospho)adenylate] typical for representatives of picornavirus-like viruses. The chemical mechanisms of viral RNA binding to VPg are discussed on the basis of the model structure of protein-RNA complex.

DNA-mediated immunization in a transgenic mouse model of the hepatitis B surface antigen chronic carrier state.

Transgenic mice expressing the sequences coding for the envelope proteins of the hepatitis B virus (HBV) in the liver have been used as a model of the HBV chronic carrier state. We evaluated the possibility of inducing a specific immune response to the viral envelope antigens and thus potentially controlling chronic HBV infection. Using HBV-specific DNA-mediated immunization in this transgenic model, we show that the immune response induced after a single intramuscular injection of DNA resulted in the complete clearance of circulating hepatitis B surface antigen and in the long-term control of transgene expression in hepatocytes. This response does not involve a detectable cytopathic effect in the liver. Adoptive transfer of fractionated primed spleen cells from DNA-immunized mice shows that T cells are responsible for the down-regulation of HBV mRNA in the liver of transgenic mice. To our knowledge, this is the first demonstration of a potential immunotherapeutic application of DNA-mediated immunization against an infectious disease and raises the possibility of designing more effective ways of treating HBV chronic carriers.

Evolutionarily conserved Galphabetagamma binding surfaces support a model of the G protein-receptor complex.

The pivotal role of G proteins in sensory, hormonal, inflammatory, and proliferative responses has provoked intense interest in understanding how they interact with their receptors and effectors. Nonetheless, the locations of the receptors and effector binding sites remain poorly characterized, although nearly complete structures of the alphabetagamma heterotrimeric complex are available. Here we apply evolutionary trace (ET) analysis [Lichtarge, O., Bourne, H. R. & Cohen, F. E. (1996) J. Mol. Biol. 257, 342-358] to propose plausible locations for these sites. On each subunit, ET identifies evolutionarily selected surfaces composed of residues that do not vary within functional subgroups and that form spatial clusters. Four clusters correctly identify subunit interfaces, and additional clusters on Galpha point to likely receptor or effector binding sites. Our results implicate the conformationally variable region of Galpha in an effector binding role. Furthermore the range of predicted interactions between the receptor and Galphabetagamma, is sufficiently limited that we can build a low resolution and testable model of the receptor-G protein complex.

Cultural universals: Measuring the semantic structure of emotion terms in English and Japanese

Research is presented on the semantic structure of 15 emotion terms as measured by judged-similarity tasks for monolingual English-speaking and monolingual and bilingual Japanese subjects. A major question is the relative explanatory power of a single shared model for English and Japanese versus culture-specific models for each language. The data support a shared model for the semantic structure of emotion terms even though some robust and significant differences are found between English and Japanese structures. The Japanese bilingual subjects use a model more like English when performing tasks in English than when performing the same task in Japanese.

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.

A Polymer Model for the Structural Organization of Chromatin Loops and Minibands in Interphase Chromosomes

A quantitative model of interphase chromosome higher-order structure is presented based on the isochore model of the genome and results obtained in the field of copolymer research. G1 chromosomes are approximated in the model as multiblock copolymers of the 30-nm chromatin fiber, which alternately contain two types of 0.5- to 1-Mbp blocks (R and G minibands) differing in GC content and DNA-bound proteins. A G1 chromosome forms a single-chain string of loop clusters (micelles), with each loop ∼1–2 Mbp in size. The number of ∼20 loops per micelle was estimated from the dependence of geometrical versus genomic distances between two points on a G1 chromosome. The greater degree of chromatin extension in R versus G minibands and a difference in the replication time for these minibands (early S phase for R versus late S phase for G) are explained in this model as a result of the location of R minibands at micelle cores and G minibands at loop apices. The estimated number of micelles per nucleus is close to the observed number of replication clusters at the onset of S phase. A relationship between chromosomal and nuclear sizes for several types of higher eukaryotic cells (insects, plants, and mammals) is well described through the micelle structure of interphase chromosomes. For yeast cells, this relationship is described by a linear coil configuration of chromosomes.

A model of long-term memory storage in the cerebellar cortex: A possible role for plasticity at parallel fiber synapses onto stellate/basket interneurons

By evoking changes in climbing fiber activity, movement errors are thought to modify synapses from parallel fibers onto Purkinje cells (pf*Pkj) so as to improve subsequent motor performance. Theoretical arguments suggest there is an intrinsic tradeoff, however, between motor adaptation and long-term storage. Assuming a baseline rate of motor errors is always present, then repeated performance of any learned movement will generate a series of climbing fiber-mediated corrections. By reshuffling the synaptic weights responsible for any given movement, such corrections will degrade the memories for other learned movements stored in overlapping sets of synapses. The present paper shows that long-term storage can be accomplished by a second site of plasticity at synapses from parallel fibers onto stellate/basket interneurons (pf*St/Bk). Plasticity at pf*St/Bk synapses can be insulated from ongoing fluctuations in climbing fiber activity by assuming that changes in pf*St/Bk synapses occur only after changes in pf*Pkj synapses have built up to a threshold level. Although climbing fiber-dependent plasticity at pf*Pkj synapses allows for the exploration of novel motor strategies in response to changing environmental conditions, plasticity at pf*St/Bk synapses transfers successful strategies to stable long-term storage. To quantify this hypothesis, both sites of plasticity are incorporated into a dynamical model of the cerebellar cortex and its interactions with the inferior olive. When used to simulate idealized motor conditioning trials, the model predicts that plasticity develops first at pf*Pkj synapses, but with additional training is transferred to pf*St/Bk synapses for long-term storage.

The "allosteric three-site model" of elongation cannot be confirmed in a well-defined ribosome system from Escherichia coli.

For the functional role of the ribosomal tRNA exit (E) site, two different models have been proposed. It has been suggested that transient E-site binding of the tRNA leaving the peptidyl (P) site promotes elongation factor G (EF-G)-dependent translocation by lowering the energetic barrier of tRNA release [Lill, R., Robertson, J. M. & Wintermeyer, W. (1989) EMBO J. 8, 3933-3938]. The alternative "allosteric three-site model" [Nierhaus, K.H. (1990) Biochemistry 29, 4997-5008] features stable, codon-dependent tRNA binding to the E site and postulates a coupling between E and aminoacyl (A) sites that regulates the tRNA binding affinity of the two sites in an anticooperative manner. Extending our testing of the two conflicting models, we have performed translocation experiments with fully active ribosomes programmed with heteropolymeric mRNA. The results confirm that the deacylated tRNA released from the P site is bound to the E site in a kinetically labile fashion, and that the affinity of binding, i.e., the occupancy of the E site, is increased by Mg2+ or polyamines. At conditions of high E-site occupancy in the posttranslocation complex, filling the A site with aminoacyl-tRNA had no influence on the E site, i.e., there was no detectable anticooperative coupling between the two sites, provided that second-round translocation was avoided by removing EF-G. On the basis of these results, which are entirely consistent with our previous results, we consider the allosteric three-site model of elongation untenable. Rather, as proposed earlier, the E site-bound state of the leaving tRNA is a transient intermediate and, as such, is a mechanistic feature of the classic two-state model of the elongating ribosome.

Broken symmetry and strangeness of the semiconductor impurity band metal-insulator transition

The filamentary model of the metal-insulator transition in randomly doped semiconductor impurity bands is geometrically equivalent to similar models for continuous transitions in dilute antiferromagnets and even to the λ transition in liquid He, but the critical behaviors are different. The origin of these differences lies in two factors: quantum statistics and the presence of long range Coulomb forces on both sides of the transition in the electrical case. In the latter case, in addition to the main transition, there are two satellite transitions associated with disappearance of the filamentary structure in both insulating and metallic phases. These two satellite transitions were first identified by Fritzsche in 1958, and their physical origin is explained here in geometrical and topological terms that facilitate calculation of critical exponents.

Antiviral agent blocks breathing of the common cold virus

A dynamic capsid is critical to the events that shape the viral life cycle; events such as cell attachment, cell entry, and nucleic acid release demand a highly mobile viral surface. Protein mass mapping of the common cold virus, human rhinovirus 14 (HRV14), revealed both viral structural dynamics and the inhibition of such dynamics with an antiviral agent, WIN 52084. Viral capsid digestion fragments resulting from proteolytic time-course experiments provided structural information in good agreement with the HRV14 three-dimensional crystal structure. As expected, initial digestion fragments included peptides from the capsid protein VP1. This observation was expected because VP1 is the most external viral protein. Initial digestion fragments also included peptides belonging to VP4, the most internal capsid protein. The mass spectral results together with x-ray crystallography data provide information consistent with a “breathing” model of the viral capsid. Whereas the crystal structure of HRV14 shows VP4 to be the most internal capsid protein, mass spectral results show VP4 fragments to be among the first digestion fragments observed. Taken together this information demonstrates that VP4 is transiently exposed to the viral surface via viral breathing. Comparative digests of HRV14 in the presence and absence of WIN 52084 revealed a dramatic inhibition of digestion. These results indicate that the binding of the antiviral agent not only causes local conformational changes in the drug binding pocket but actually stabilizes the entire viral capsid against enzymatic degradation. Viral capsid mass mapping provides a fast and sensitive method for probing viral structural dynamics as well as providing a means for investigating antiviral drug efficacy.

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.

Multiple loop structures critical for ligand binding of the integrin α4 subunit in the upper face of the β-propeller mode 1

A non-I-domain integrin, α4β1, recognizes vascular cell adhesion molecule 1 (VCAM-1) and the IIICS portion of fibronectin. To localize regions of α4 critical for ligand binding, we swapped several predicted loops within or near the putative ligand-binding site of α4 (which spans repeats 2–5 of the seven N-terminal repeats) with the corresponding regions of α5. Swapping residues 112–131 in repeat 2, or residues 237–247 in repeat 4, completely blocked adhesion to immobilized VCAM-1 and connecting segment 1 (CS-1) peptide. However, swapping residues 40–52 in repeat 1, residues 151–164 in repeat 3, or residues 282–288 (which contain a putative cation binding motif) in repeat 5 did not affect or only slightly reduced adhesion to these ligands. The binding of several function-blocking antibodies is blocked by swapping residues 112–131, 151–164, and 186–191 (which contain previously identified residues critical for ligand binding, Tyr-187 and Gly-190). These results are consistent with the recently published β-propeller folding model of the integrin α4 subunit [Springer, T. A. (1997) Proc. Natl. Acad. Sci. USA 94, 65–72], in which seven four-stranded β-sheets are arranged in a torus around a pseudosymmetric axis. The regions of α4 critical for ligand binding are adjacent to each other and are located in the upper face, the predicted ligand-binding site, of the β-propeller model, although they are not adjacent in the primary structure.

Dual role of the nuclear factor of activated T cells insert region in DNA recognition and cooperative contacts to activator protein 1

The transcription factors nuclear factor of activated T cells (NFAT) and activator protein 1 (AP-1) coordinately regulate cytokine gene expression in activated T-cells by binding to closely juxtaposed sites in cytokine promoters. The structural basis for cooperative binding of NFAT and AP-1 to these sites, and indeed for the cooperative binding of transcription factors to composite regulatory elements in general, is not well understood. Mutagenesis studies have identified a segment of AP-1, which lies at the junction of its DNA-binding and dimerization domains (basic region and leucine zipper, respectively), as being essential for protein–protein interactions with NFAT in the ternary NFAT/AP-1/DNA complex. In a model of the ternary complex, the segment of NFAT nearest AP-1 is the Rel insert region (RIR), a feature that is notable for its hypervariability in size and in sequence amongst members of the Rel transcription factor family. Here we have used mutational analysis to study the role of the NFAT RIR in binding to DNA and AP-1. Parallel yeast one-hybrid screening assays in combination with alanine-scanning mutagenesis led to the identification of four amino acid residues in the RIR of NFAT2 (also known as NFATC1 or NFATc) that are essential for cooperativity with AP-1 (Ile-544, Glu-545, Thr-551, and Ile-553), and three residues that are involved in interactions with DNA (Lys-538, Arg-540, and Asn-541). These results were confirmed and extended through in vitro binding assays. We thus conclude that the NFAT RIR plays an essential dual role in DNA recognition and cooperative binding to AP-1 family transcription factors.

Human fatty acid synthase: Assembling recombinant halves of the fatty acid synthase subunit protein reconstitutes enzyme activity

Our model of the native fatty acid synthase (FAS) depicts it as a dimer of two identical multifunctional proteins (Mr ≈ 272,000) arranged in an antiparallel configuration so that the active Cys-SH of the β-ketoacyl synthase of one subunit (where the acyl group is attached) is juxtaposed within 2 Å of the pantetheinyl-SH of the second subunit (where the malonyl group is bound). This arrangement generates two active centers for fatty acid synthesis and predicts that if we have two appropriate halves of the monomer, we should be able to reconstitute an active fatty acid-synthesizing site. We cloned, expressed, and purified catalytically active thioredoxin (TRX) fusion proteins of the NH2-terminal half of the human FAS subunit protein (TRX-hFAS-dI; residues 1–1,297; Mr ≈ 166) and of the C-terminal half (TRX-hFAS-dII-III; residues 1,296–2,504; Mr ≈ 155). Adding equivalent amounts of TRX-hFAS-dI and TRX-hFAS-dII-III to a reaction mixture containing acetyl-CoA, malonyl-CoA, and NADPH resulted in the synthesis of long-chain fatty acids. The rate of synthesis was dependent upon the presence of both recombinant proteins and reached a constant level when they were present in equivalent amounts, indicating that the reconstitution of an active fatty acid-synthesizing site required the presence of every partial activity associated with the subunit protein. Analyses of the product acids revealed myristate to be the most abundant with small amounts of palmitate and stearate, possibly because of the way the fused recombinant proteins interacted with each other so that the thioesterase hydrolyzed the acyl group in its myristoyl state. The successful reconstitution of the human FAS activity from its domain I and domains II and III fully supports our model for the structure–function relationship of FAS in animal tissues.