Solution structure of a 142-residue recombinant prion protein corresponding to the infectious fragment of the scrapie isoform

The scrapie prion protein (PrPSc) is the major, and possibly the only, component of the infectious prion; it is generated from the cellular isoform (PrPC) by a conformational change. N-terminal truncation of PrPSc by limited proteolysis produces a protein of ≈142 residues designated PrP 27–30, which retains infectivity. A recombinant protein (rPrP) corresponding to Syrian hamster PrP 27–30 was expressed in Escherichia coli and purified. After refolding rPrP into an α-helical form resembling PrPC, the structure was solved by multidimensional heteronuclear NMR, revealing many structural features of rPrP that were not found in two shorter PrP fragments studied previously. Extensive side-chain interactions for residues 113–125 characterize a hydrophobic cluster, which packs against an irregular β-sheet, whereas residues 90–112 exhibit little defined structure. Although identifiable secondary structure is largely lacking in the N terminus of rPrP, paradoxically this N terminus increases the amount of secondary structure in the remainder of rPrP. The surface of a long helix (residues 200–227) and a structured loop (residues 165–171) form a discontinuous epitope for binding of a protein that facilitates PrPSc formation. Polymorphic residues within this epitope seem to modulate susceptibility of sheep and humans to prion disease. Conformational heterogeneity of rPrP at the N terminus may be key to the transformation of PrPC into PrPSc, whereas the discontinuous epitope near the C terminus controls this transition.

Two functionally distinct subsites for the binding of internal blockers to the pore of voltage-activated K+ channels

Many blockers of Na+ and K+ channels act by blocking the pore from the intracellular side. For Shaker K+ channels, such intracellular blockers vary in their functional effect on slow (C-type) inactivation: Some blockers interfere with C-type inactivation, whereas others do not. These functional differences can be explained by supposing that there are two overlapping “subsites” for blocker binding, only one of which inhibits C-type inactivation through an allosteric effect. We find that the ability to bind to these subsites depends on specific structural characteristics of the blockers, and correlates with the effect of mutations in two distinct regions of the channel protein. These interactions are important because they affect the ability of blockers to produce use-dependent inhibition.

Unwinding of nucleic acids by HCV NS3 helicase is sensitive to the structure of the duplex

Hepatitis C virus (HCV) helicase, non-structural protein 3 (NS3), is proposed to aid in HCV genome replication and is considered a target for inhibition of HCV. In order to investigate the substrate requirements for nucleic acid unwinding by NS3, substrates were prepared by annealing a 30mer oligonucleotide to a 15mer. The resulting 15 bp duplex contained a single-stranded DNA overhang of 15 nt referred to as the bound strand. Other substrates were prepared in which the 15mer DNA was replaced by a strand of peptide nucleic acid (PNA). The PNA–DNA substrate was unwound by NS3, but the observed rate of strand separation was at least 25-fold slower than for the equivalent DNA–DNA substrate. Binding of NS3 to the PNA–DNA substrate was similar to the DNA–DNA substrate, due to the fact that NS3 initially binds to the single-stranded overhang, which was identical in each substrate. A PNA–RNA substrate was not unwound by NS3 under similar conditions. In contrast, morpholino–DNA and phosphorothioate–DNA substrates were utilized as efficiently by NS3 as DNA–DNA substrates. These results indicate that the PNA–DNA and PNA–RNA heteroduplexes adopt structures that are unfavorable for unwinding by NS3, suggesting that the unwinding activity of NS3 is sensitive to the structure of the duplex.

Does a neuroimmune interaction contribute to the genesis of painful peripheral neuropathies?

Painful peripheral neuropathies are precipitated by nerve injury from disease or trauma. All such injuries will be accompanied by an inflammatory reaction, a neuritis, that will mobilize the immune system. The role of the inflammation itself is difficult to determine in the presence of structural damage to the nerve. A method has been devised to produce a focal neuritis in the rat sciatic nerve that involves no more than trivial structural damage to the nerve. This experimental focal neuritis produces neuropathic pain sensations (heat- and mechano-hyperalgesia, and cold- and mechano-allodynia) in the ipsilateral hind paw. The abnormal pain sensations begin in 1–2 days and last for 4–6 days, with a subsequent return to normal. These results suggest that there is a neuroimmune interaction that occurs at the outset of nerve injury (and perhaps episodically over time in slow developing conditions like diabetic neuropathy) that produces neuropathic pain. The short duration of the phenomena suggest that they may prime the system for more slowly developing mechanisms of abnormal pain (e.g., ectopic discharge in axotomized primary afferent neurons) that underlie the chronic phase of painful neuropathy.

Characterization of a Novel Human SMC Heterodimer Homologous to the Schizosaccharomyces pombe Rad18/Spr18 Complex

The structural maintenance of chromosomes (SMC) protein encoded by the fission yeast rad18 gene is involved in several DNA repair processes and has an essential function in DNA replication and mitotic control. It has a heterodimeric partner SMC protein, Spr18, with which it forms the core of a multiprotein complex. We have now isolated the human orthologues of rad18 and spr18 and designated them hSMC6 and hSMC5. Both proteins are about 1100 amino acids in length and are 27–28% identical to their fission yeast orthologues, with much greater identity within their N- and C-terminal globular domains. The hSMC6 and hSMC5 proteins interact to form a tight complex analogous to the yeast Rad18/Spr18 heterodimer. In proliferating human cells the proteins are bound to both chromatin and the nucleoskeleton. In addition, we have detected a phosphorylated form of hSMC6 that localizes to interchromatin granule clusters. Both the total level of hSMC6 and its phosphorylated form remain constant through the cell cycle. Both hSMC5 and hSMC6 proteins are expressed at extremely high levels in the testis and associate with the sex chromosomes in the late stages of meiotic prophase, suggesting a possible role for these proteins in meiosis.

Structural studies of p21Waf1/Cip1/Sdi1 in the free and Cdk2-bound state: conformational disorder mediates binding diversity.

The cyclin-dependent kinase (Cdk) inhibitor p21Waf1/Cip1/Sdi1, important for p53-dependent cell cycle control, mediates G1/S arrest through inhibition of Cdks and possibly through inhibition of DNA replication. Cdk inhibition requires a sequence of approximately 60 amino acids within the p21 NH2 terminus. We show, using proteolytic mapping, circular dichroism spectropolarimetry, and nuclear magnetic resonance spectroscopy, that p21 and NH2-terminal fragments that are active as Cdk inhibitors lack stable secondary or tertiary structure in the free solution state. In sharp contrast to the disordered free state, however, the p21 NH2 terminus adopts an ordered stable conformation when bound to Cdk2, as shown directly by NMR spectroscopy. We have, thus, identified a striking disorder-order transition for p21 upon binding to one of its biological targets, Cdk2. This structural transition has profound implications in light of the ability of p21 to bind and inhibit a diverse family of cyclin-Cdk complexes, including cyclin A-Cdk2, cyclin E-Cdk2, and cyclin D-Cdk4. Our findings suggest that the flexibility, or disorder, of free p21 is associated with binding diversity and offer insights into the role for structural disorder in mediating binding specificity in biological systems. Further, these observations challenge the generally accepted view of proteins that stable secondary and tertiary structure are prerequisites for biological activity and suggest that a broader view of protein structure should be considered in the context of structure-activity relationships.

Three distinct structural environments of a transmembrane domain in the inwardly rectifying potassium channel ROMK1 defined by perturbation.

To probe the protein environment of an ion channel, we have perturbed the structure of a transmembrane domain by substituting side chains with those of two different sizes by using site-specific mutagenesis. We have used Trp and Ala as a high- and a low-impact perturbation probe, respectively, to replace each of 18 consecutive residues within the putative second transmembrane segment, M2, of an inwardly rectifying potassium channel, ROMK1. Our rationale is that a change in the channel function as a consequence of these mutations at a particular position will reflect the structural environment of the altered side chain. Each position can then be assigned to one of three classes of environments, as grated by different levels of perturbation: very tolerant (channel functions with both Trp and Ala substitutions), tolerant (function preserved with Ala but not with Trp substitution), and intolerant (either Ala or Trp substitution destroys function). We identify the very tolerant environment as being lipid-facing, tolerant as protein-interior-facing, and intolerant as pore-facing. We observe a strikingly ordered pattern of perturbation of all three environmental classes. This result indicates that M2 is a straight alpha-helix.

The new dysmorphology: application of insights from basic developmental biology to the understanding of human birth defects.

Information obtained from studies of developmental and cellular processes in lower organisms is beginning to make significant contributions to the understanding of the pathogenesis of human birth defects, and it is now becoming possible to treat birth defects as inborn errors of development. Mutations in genes for transcription factors, receptors, cell adhesion molecules, intercellular junctions, molecules involved in signal transduction, growth factors, structural proteins, enzymes, and transporters have been identified in genetically caused human malformations and dysplasias. The identification of these mutations and the analysis of their developmental effects have been greatly facilitated by the existence of natural or engineered models in the mouse and even of related mutations in Drosophila, and in some instances a remarkable conservation of function in development has been observed, even between widely separated species.

Mathematical analysis of activation thresholds in enzyme-catalyzed positive feedbacks: application to the feedbacks of blood coagulation.

A hierarchy of enzyme-catalyzed positive feedback loops is examined by mathematical and numerical analysis. Four systems are described, from the simplest, in which an enzyme catalyzes its own formation from an inactive precursor, to the most complex, in which two sequential feedback loops act in a cascade. In the latter we also examine the function of a long-range feedback, in which the final enzyme produced in the second loop activates the initial step in the first loop. When the enzymes generated are subject to inhibition or inactivation, all four systems exhibit threshold properties akin to excitable systems like neuron firing. For those that are amenable to mathematical analysis, expressions are derived that relate the excitation threshold to the kinetics of enzyme generation and inhibition and the initial conditions. For the most complex system, it was expedient to employ numerical simulation to demonstrate threshold behavior, and in this case long-range feedback was seen to have two distinct effects. At sufficiently high catalytic rates, this feedback is capable of exciting an otherwise subthreshold system. At lower catalytic rates, where the long-range feedback does not significantly affect the threshold, it nonetheless has a major effect in potentiating the response above the threshold. In particular, oscillatory behavior observed in simulations of sequential feedback loops is abolished when a long-range feedback is present.

Genetic variation in vulnerability to the behavioral effects of neonatal hippocampal damage in rats.

We explored how two independent variables, one genetic (i.e., specific rat strains) and another environmental (i.e., a developmental excitotoxic hippocampal lesion), contribute to phenotypic variation. Sprague-Dawley (SD), Fischer 344 (F344), and Lewis rats underwent two grades of neonatal excitotoxic damage: small and large ventral hippocampal (SVH and LVH) lesions. Locomotion was tested before puberty [postnatal day 35 (P35)] and after puberty (P56) following exposure to a novel environment or administration of amphetamine. The behavioral effects were strain- and lesion-specific. As shown previously, SD rats with LVH lesions displayed enhanced spontaneous and amphetamine-induced locomotion as compared with controls at P56, but not at P35. SVH lesions in SD rats had no effect at any age. In F344 rats with LVH lesions, enhanced spontaneous and amphetamine-induced locomotion appeared early (P35) and was exaggerated at P56. SVH lesions in F344 rats resulted in a pattern of effects analogous to LVH lesions in SD rats--i.e., postpubertal onset of hyperlocomotion (P56). In Lewis rats, LVH lesions had no significant effect on novelty- or amphetamine-induced locomotion at any age. These data show that the degree of genetic predisposition and the extent of early induced hippocampal defect contribute to the particular pattern of behavioral outcome. These results may have implications for modeling interactions of genetic and environmental factors involved in schizophrenia, a disorder characterized by phenotypic heterogeneity, genetic predisposition, a developmental hippocampal abnormality, and vulnerability to environmental stress.

Identification of a plant serine-arginine-rich protein similar to the mammalian splicing factor SF2/ASF.

We show that the higher plant Arabidopsis thaliana has a serine-arginine-rich (SR) protein family whose members contain a phosphoepitope shared by the animal SR family of splicing factors. In addition, we report the cloning and characterization of a cDNA encoding a higher-plant SR protein from Arabidopsis, SR1, which has striking sequence and structural homology to the human splicing factor SF2/ASF. Similar to SF2/ASF, the plant SR1 protein promotes splice site switching in mammalian nuclear extracts. A novel feature of the Arabidopsis SR protein is a C-terminal domain containing a high concentration of proline, serine, and lysine residues (PSK domain), a composition reminiscent of histones. This domain includes a putative phosphorylation site for the mitotic kinase cyclin/p34cdc2.

Abnormal junctions between surface membrane and sarcoplasmic reticulum in skeletal muscle with a mutation targeted to the ryanodine receptor.

Junctions that mediate excitation-contraction (e-c) coupling are formed between the sarcoplasmic reticulum (SR) and either the surface membrane or the transverse (T) tubules in normal skeletal muscle. Two structural components of the junctions, the feet of the SR and the tetrads of T tubules, have been identified respectively as ryanodine receptors (RyRs, or SR calcium-release channels), and as groups of four dihydropyridine receptors (DHPRs, or voltage sensors of e-c coupling). A targeted mutation (skrrm1) of the gene for skeletal muscle RyRs in mice results in the absence of e-c coupling in homozygous offspring of transgenic parents. The mutant gene is expected to produce no functional RyRs, and we have named the mutant mice "dyspedic" because they lack feet--the cytoplasmic domain of RyRs anchored in the SR membrane. We have examined the development of junctions in skeletal muscle fibers from normal and dyspedic embryos. Surprisingly, despite the absence of RyRs, junctions are formed in dyspedic myotubes, but the junctional gap between the SR and T tubule is narrow, presumably because the feet are missing. Tetrads are also absent from these junctions. The results confirm the identity of RyRs and feet and a major role for RyRs and tetrads in e-c coupling. Since junctions form in the absence of feet and tetrads, coupling of SR to surface membrane and T tubules appears to be mediated by additional proteins, distinct from either RyRs or DHPRs.

Finite temperature dynamics of vortices in the two dimensional anisotropic Heisenberg model

We study the effects of finite temperature on the dynamics of non-planar vortices in the classical, two-dimensional anisotropic Heisenberg model with XY- or easy-plane symmetry. To this end, we analyze a generalized Landau-Lifshitz equation including additive white noise and Gilbert damping. Using a collective variable theory with no adjustable parameters we derive an equation of motion for the vortices with stochastic forces which are shown to represent white noise with an effective diffusion constant linearly dependent on temperature. We solve these stochastic equations of motion by means of a Green's function formalism and obtain the mean vortex trajectory and its variance. We find a non-standard time dependence for the variance of the components perpendicular to the driving force. We compare the analytical results with Langevin dynamics simulations and find a good agreement up to temperatures of the order of 25% of the Kosterlitz-Thouless transition temperature. Finally, we discuss the reasons why our approach is not appropriate for higher temperatures as well as the discreteness effects observed in the numerical simulations.

Integration of different models in the design of chemical processes: Application to the design of a power plant

With advances in the synthesis and design of chemical processes there is an increasing need for more complex mathematical models with which to screen the alternatives that constitute accurate and reliable process models. Despite the wide availability of sophisticated tools for simulation, optimization and synthesis of chemical processes, the user is frequently interested in using the ‘best available model’. However, in practice, these models are usually little more than a black box with a rigid input–output structure. In this paper we propose to tackle all these models using generalized disjunctive programming to capture the numerical characteristics of each model (in equation form, modular, noisy, etc.) and to deal with each of them according to their individual characteristics. The result is a hybrid modular–equation based approach that allows synthesizing complex processes using different models in a robust and reliable way. The capabilities of the proposed approach are discussed with a case study: the design of a utility system power plant that has been decomposed into its constitutive elements, each treated differently numerically. And finally, numerical results and conclusions are presented.

A hydrological–economic model for sustainable groundwater use in sparse-data drylands: Application to the Amtoudi Oasis in southern Morocco, northern Sahara

A hydrological–economic model is introduced to describe the dynamics of groundwater-dependent economics (agriculture and tourism) for sustainable use in sparse-data drylands. The Amtoudi Oasis, a remote area in southern Morocco, in the northern Sahara attractive for tourism and with evidence of groundwater degradation, was chosen to show the model operation. Governing system variables were identified and put into action through System Dynamics (SD) modeling causal diagrams to program basic formulations into a model having two modules coupled by the nexus ‘pumping’: (1) the hydrological module represents the net groundwater balance (G) dynamics; and (2) the economic module reproduces the variation in the consumers of water, both the population and tourists. The model was operated under similar influx of tourists and different scenarios of water availability, such as the wet 2009–2010 and the average 2010–2011 hydrological years. The rise in international tourism is identified as the main driving force reducing emigration and introducing new social habits in the population, in particular concerning water consumption. Urban water allotment (PU) was doubled for less than a 100-inhabitant net increase in recent decades. The water allocation for agriculture (PI), the largest consumer of water, had remained constant for decades. Despite that the 2-year monitoring period is not long enough to draw long-term conclusions, groundwater imbalance was reflected by net aquifer recharge (R) less than PI + PU (G < 0) in the average year 2010–2011, with net lateral inflow from adjacent Cambrian formations being the largest recharge component. R is expected to be much less than PI + PU in recurrent dry spells. Some low-technology actions are tentatively proposed to mitigate groundwater degradation, such as: wastewater capture, treatment, and reuse for irrigation; storm-water harvesting for irrigation; and active maintenance of the irrigation system to improve its efficiency.