The minimal model of the hypothalamic–pituitary–adrenal axis

This paper concerns ODE modeling of the hypothalamic–pituitary– adrenal axis (HPA axis) using an analytical and numerical approach, combined with biological knowledge regarding physiological mechanisms and parameters. The three hormones, CRH, ACTH, and cortisol, which interact in the HPA axis are modeled as a system of three coupled, nonlinear differential equations. Experimental data shows the circadian as well as the ultradian rhythm. This paper focuses on the ultradian rhythm. The ultradian rhythm can mathematically be explained by oscillating solutions. Oscillating solutions to an ODE emerges from an unstable fixed point with complex eigenvalues with a positive real parts and a non-zero imaginary parts. The first part of the paper describes the general considerations to be obeyed for a mathematical model of the HPA axis. In this paper we only include the most widely accepted mechanisms that influence the dynamics of the HPA axis, i.e. a negative feedback from cortisol on CRH and ACTH. Therefore we term our model the minimal model. The minimal model, encompasses a wide class of different realizations, obeying only a few physiologically reasonable demands. The results include the existence of a trapping region guaranteeing that concentrations do not become negative or tend to infinity. Furthermore, this treatment guarantees the existence of a unique fixed point. A change in local stability of the fixed point, from stable to unstable, implies a Hopf bifurcation; thereby, oscillating solutions may emerge from the model. Sufficient criteria for local stability of the fixed point, and an easily applicable sufficient criteria guaranteeing global stability of the fixed point, is formulated. If the latter is fulfilled, ultradian rhythm is an impossible outcome of the minimal model and all realizations thereof. The second part of the paper concerns a specific realization of the minimal model in which feedback functions are built explicitly using receptor dynamics. Using physiologically reasonable parameter values, along with the results of the general case, it is demonstrated that un-physiological values of the parameters are needed in order to achieve local instability of the fixed point. Small changes in physiologically relevant parameters cause the system to be globally stable using the analytical criteria. All simulations show a globally stable fixed point, ruling out periodic solutions even when an investigation of the ‘worst case parameters’ is performed.     

The molecular dissection of the chaperone–usher pathway

Secretion systems are specialized in transport of proteins, DNA or nutrients across the cell envelope of bacteria and enable them to communicate with their environment. The chaperone–usher (CU) pathway is used for assembly and secretion of a large family of long adhesive protein polymers, termed pili, and is widespread among Gram-negative pathogens [1]. Moreover, recent evidence has indicated that CU secretion systems are also involved in sporulation  and . In this review we focus on the structural biology of the paradigmatic type 1 and P pili CU systems encoded by uropathogenic Escherichia coli (UPEC), where recent progress has provided unprecedented insights into pilus assembly and secretion mechanism. This article is part of a Special Issue entitled: Protein trafficking and secretion in bacteria. Guest Editors: Anastassios Economou and Ross Dalbey.     

The kinetics of the reversible inhibition of heart lactate dehydrogenase through the formation of the enzyme–oxidized nicotinamide–adenine dinucleotide–pyruvate compound

The inhibition of lactate dehydrogenase at high pyruvate concentration was studied in three ways. First, a rapid decrease in the rate of the enzyme reaction was observed; secondly, the rate of formation of a pyruvate-NAD(+) compound was followed by the change in E(325); thirdly, the rate of quenching of the protein fluorescence was measured. The data obtained at pH6.0 at different temperatures and ionic strengths as functions of pyruvate, NAD(+) and enzyme concentrations show that the extent of inhibition can be correlated with the reversible formation of a compound between pyruvate and enzyme-bound NAD(+). It is suggested that the detailed kinetic analysis of the formation of this abortive ternary compound will give pertinent information about properties of the enzyme-NAD(+) compound involved in the normal catalytic process.     

The mechanism of the periodate–thiobarbituric acid reaction of sialic acids

1. The chromogen formation from N-acetylneuraminic acid in the periodate-thiobarbituric acid reaction was investigated. Measurement of periodate consumption showed an uptake of approx. 3moles/mole of substrate in neutral as well as in strongly acidic solution. Therefore the chromogen beta-formylpyruvic acid is not a direct product of the periodate oxidation; it is presumed to be formed from the true oxidation product, a hexos-5-uluronic acid, by aldol splitting during the reaction in hot acidic solution with thiobarbituric acid. 2. Methyl (methyl beta-l-threo-hexos-4-enepyranosid)uronate, an analogue of the pre-chromogen, has been shown to yield with thiobarbituric acid in acidic solution a pigment exhibiting an identical absorption spectrum and showing the same behaviour on paper chromatography as the pigment obtained from N-acetylneuraminic acid in the periodate-thiobarbituric acid assay. 3. The substitution at C-2 of methoxyneuraminic acid does not inhibit the periodate-thiobarbituric acid reaction. In neutral solution methoxyneuraminic acid is oxidized by periodate to a substance that reacts readily with thiobarbituric acid in acidic solution. When periodate oxidation is attempted in acidic solution, protonation of the amino group protects this group against oxidation, rendering methoxyneuraminic acid negative in the assay systems of Warren (1959a,b) and Aminoff (1959, 1961).     

Does DNA alone determine the development of the organism? (the information aspect of the problem)

Two closely related controversial problems are discussed: whether the developmental processes can be reduced to the synthesis of polypeptides encoded in DNA, and whether the information in DNA is equivalent to that in the adult organism. Critically considered are the ideas that DNA is only responsible for the protein synthesis, whereas morphogenesis proceeds independently and according to epigenetic regularities of its own. It is stated that development is the realization of genetic information in which more elementary (molecular) processes unambiguously determine a more complex cellular level which in its turn determines morphogenesis of tissues and organs. Various mechanisms of the appearance of new information in the course of development are considered. The statement is made that new information concerns only some individual characters of the organism, whereas most of information that determines the process of development and the structure of the adult organism is created in the course of evolution, is stored in DNA and inherited.     

The Impact of the Species–Area Relationship on Estimates of Paleodiversity

Estimates of paleodiversity patterns through time have relied on datasets that lump taxonomic occurrences from geographic areas of varying size per interval of time. In essence, such estimates assume that the species–area effect, whereby more species are recorded from larger geographic areas, is negligible for fossil data. We tested this assumption by using the newly developed Miocene Mammal Mapping Project database of western North American fossil mammals and its associated analysis tools to empirically determine the geographic area that contributed to species diversity counts in successive temporal bins. The results indicate that a species–area effect markedly influences counts of fossil species, just as variable spatial sampling influences diversity counts on the modern landscape. Removing this bias suggests some traditionally recognized peaks in paleodiversity are just artifacts of the species–area effect while others stand out as meriting further attention. This discovery means that there is great potential for refining existing time-series estimates of paleodiversity, and for using species–area relationships to more reliably understand the magnitude and timing of such biotically important events as extinction, lineage diversification, and long-term trends in ecological structure.     

The role of the abnormalities in the distal pathway of cholesterol biosynthesis in the Conradi–Hünermann–Happle syndrome

Conradi–Hünermann–Happle syndrome (CDPX2, OMIM 302960) is an inherited X-linked dominant variant of chondrodysplasia punctata (CP) caused by mutations in one gene of the distal pathway of cholesterol biosynthesis. It exhibits intense phenotypic variation and primarily affects the skin, bones and eyes. The ichthyosis following Blaschko's lines, chondrodysplasia punctata and cataracts are the typical clinical findings. The cardinal biochemical features are an increase in 8(9)-cholestenol and 8-dehydrocholesterol (8DHC), which suggest a deficiency in 3β-hydroxysteroid-Δ8,Δ7-isomerase, also called emopamil binding protein (EBP). The EBP gene is located on the short arm of the X chromosome (Xp11.22–p11.23) and encodes a 230 amino acid protein with dual function. Explaining the clinical phenotype in CDPX2 implies an understanding of both the genetics and biochemical features of this disease. CDPX2 displays an X-linked dominant pattern of inheritance, which is responsible for the distribution of lesions in some tissues. The clinical phenotype in CDPX2 results directly from impairment in cholesterol biosynthesis, and indirectly from abnormalities in the hedgehog signaling protein pathways. This article is part of a Special Issue entitled The Important Role of Lipids in the Epidermis and their Role in the Formation and Maintenance of the Cutaneous Barrier. Guest Editors: Kenneth R. Feingold and Peter Elias.     

The dual nature of haemocyanin in the establishment and persistence of the squid–vibrio symbiosis

We identified and sequenced from the squid Euprymna scolopes two isoforms of haemocyanin that share the common structural/physiological characteristics of haemocyanin from a closely related cephalopod, Sepia officinalis, including a pronounced Bohr effect. We examined the potential roles for haemocyanin in the animal''s symbiosis with the luminous bacterium Vibrio fischeri. Our data demonstrate that, as in other cephalopods, the haemocyanin is primarily synthesized in the gills. It transits through the general circulation into other tissues and is exported into crypt spaces that support the bacterial partner, which requires oxygen for its bioluminescence. We showed that the gradient of pH between the circulating haemolymph and the matrix of the crypt spaces in adult squid favours offloading of oxygen from the haemocyanin to the symbionts. Haemocyanin is also localized to the apical surfaces and associated mucus of a juvenile-specific epithelium on which the symbionts gather, and where their specificity is determined during the recruitment into the association. The haemocyanin has an antimicrobial activity, which may be involved in this enrichment of V. fischeri during symbiont initiation. Taken together, these data provide evidence that the haemocyanin plays a role in shaping two stages of the squid–vibrio partnership.     

The Structural Dynamics of the Flavivirus Fusion Peptide–Membrane Interaction

Membrane fusion is a crucial step in flavivirus infections and a potential target for antiviral strategies. Lipids and proteins play cooperative roles in the fusion process, which is triggered by the acidic pH inside the endosome. This acidic environment induces many changes in glycoprotein conformation and allows the action of a highly conserved hydrophobic sequence, the fusion peptide (FP). Despite the large volume of information available on the virus-triggered fusion process, little is known regarding the mechanisms behind flavivirus–cell membrane fusion. Here, we evaluated the contribution of a natural single amino acid difference on two flavivirus FPs, FLA_G (⁹⁸DRGWGNGCGLFGK¹¹⁰) and FLA_H (⁹⁸DRGWGNHCGLFGK¹¹⁰), and investigated the role of the charge of the target membrane on the fusion process. We used an in silico approach to simulate the interaction of the FPs with a lipid bilayer in a complementary way and used spectroscopic approaches to collect conformation information. We found that both peptides interact with neutral and anionic micelles, and molecular dynamics (MD) simulations showed the interaction of the FPs with the lipid bilayer. The participation of the indole ring of Trp appeared to be important for the anchoring of both peptides in the membrane model, as indicated by MD simulations and spectroscopic analyses. Mild differences between FLA_G and FLA_H were observed according to the pH and the charge of the target membrane model. The MD simulations of the membrane showed that both peptides adopted a bend structure, and an interaction between the aromatic residues was strongly suggested, which was also observed by circular dichroism in the presence of micelles. As the FPs of viral fusion proteins play a key role in the mechanism of viral fusion, understanding the interactions between peptides and membranes is crucial for medical science and biology and may contribute to the design of new antiviral drugs.     

The making of ethology: The association for the Study of Animal Behaviour, 1936–1986

The paper reviews the history of the Association for the Study of Animal Behaviour in its 50th anniversary year. Having discussed the founding of the Institute for the Study of Animal Behaviour in 1936, the paper establishes its principal aims, follows its fortunes before and during the war, traces its role in the expansion and institutionalization of ethology in Britainbetween 1945 and 1950, and summarizes its transformation into a flourishing scientific society in the 1950s and 1960s. It is concluded that the history of the Association is, to a striking extent, a microcosm of the making of ethology as a successful scientific discipline.     

The enthalpy change for the reduction of nicotinamide–adenine dinucleotide

The heat of the reaction NAD(+)+propan-2-ol=NADH+acetone+H(+) was determined to be 42.5+/-0.6kJ/mol (10.17+/-0.15kcal/mol) from equilibrium measurements at 9-42 degrees C catalysed by yeast alcohol dehydrogenase. With the aid of thermochemical data for acetone and propan-2-ol the values of DeltaH=-29.2kJ/mol (-6.99kcal/mol) and DeltaG(0)=22.1kJ/mol (5.28kcal/mol) are derived for the reduction of NAD (NAD(+)+H(2)=NADH+H(+)). These values are consistent with analogous but less accurate data for the ethanol-acetaldehyde reaction. Thermodynamic data for the reduction of NAD and NADP are summarized.     

The Community Structure of the European Network of Interlocking Directorates 2005–2010

The boards of directors at large European companies overlap with each other to a sizable extent both within and across national borders. This could have important economic, political and management consequences. In this work we study in detail the topological structure of the networks that arise from this phenomenon. Using a comprehensive information database, we reconstruct the implicit networks of shared directorates among the top 300 European firms in 2005 and 2010, and suggest a number of novel ways to explore the trans-nationality of such business elite networks. Powerful community detection heuristics indicate that geography still plays an important role: there exist clear communities and they have a distinct national character. Nonetheless, from 2005 to 2010 we observe a densification of the boards interlocks network and a larger transnational orientation in its communities. Together with central actors and assortativity analyses, we provide statistical evidence that, at the level of corporate governance, Europe is getting closer.     

The hypothalamic–hypophyseal–gonadal regulation of hepatic glutathione S-transferases in the rat

Hepatic glutathione S-transferase activities were determined with the substrates 1,2-dichloro-4-nitrobenzene and 1-chloro-2,4-dinitrobenzene. Sexual differentiation of glutathione S-transferase activities is not evident during the prepubertal period, but glutathione conjugation with 1,2-dichloro-4-nitrobenzene is 2–3-fold greater in adult males than in females. Glutathione conjugation with 1-chloro-2,4-dinitrobenzene is slightly higher in adult males than adult females. No change in activity was observed after postpubertal gonadectomy of males or females. Neonatal castration of males results in a significant decrease in glutathione conjugation with 1,2-dichloro-4-nitrobenzene. Hypophysectomy, or hypophysectomy followed by gonadectomy did result in significantly higher glutathione S-transferase activities in both sexes. These increases can be reversed by implanting an adult male or female pituitary or four prepubertal pituitaries under the kidney capsule. Postpubertal sexual differentiation of glutathione S-transferase activities is neither dependent on pituitary sexual differentiation nor pituitary maturation. Prolactin concentrations are inversely related to glutathione S-transferase activities in hypophysectomized rats with or without ectopic pituitaries. Somatotropin exogenously administered to hypophysectomized rats results in decreased glutathione S-transferase activities, whereas prolactin has no effect. Adult male rats treated neonatally with monosodium l-glutamate to induce arcuate nucleus lesions of the hypothalamus have decreased glutathione S-transferase activities towards 1,2-dichloro-4-nitrobenzene and decreased somatotropin concentrations. Our experiments suggests that sexual differentiation of hepatic glutathione S-transferase is a result of a hypothalamic inhibiting factor in the male (absent in the female). This postpubertally expressed inhibiting factor acts on the pituitary to prevent secretion of a pituitary inhibiting factor (autonomously secreted by the female), resulting in higher glutathione S-transferase activities in the adult male than the adult female.     

Active components of frequently used β-blockers from the aspect of computational study

The aim of this study is to investigate the active components of representative drugs for blood pressure regulation by applying quantum mechanical computer codes and comparison of the same for the sake of obtaining knowledge about the properties associated with the electronic structure of given molecules. The study included three well-known, but not theoretically investigated enough, active components of β-blockers: acebutolol, metoprolol and atenolol. The results are in agreement with the experimental data and were used for initial assumptions concerning the degradation of these compounds.     

The Role of the Transport System in the Control of the Source–Sink Relations in Pinus sylvestris

Morphophysiological correlations were studied in medium-aged (20- to 60-year-old) Scots pine trees under the northern taiga conditions. Under various ecological conditions, pine trees developed a well-balanced structure, with close linear relationships between needle and root weight and their cross-section areas in all components of the continuous transport network (the coefficient of determination was between 0.88 and 0.999). When the annual cycle of soluble and insoluble carbohydrate contents was followed in various pine tissues, the total concentrations of soluble and insoluble carbohydrates were maintained at constant and tissue-specific levels, except in the growth period. The maximum level of carbohydrates was observed in all tissues at the beginning of rapid growth, and the minimum, at growth cessation. The qualitative composition and amount of carbohydrates matched the phenological phases of development and were not affected by the ecological growth conditions pertinent to the particular environment. The authors conclude that assimilate synthesis and partitioning are related to structural development, and the state of sink centers determines the attracting capacity, whereas the transport network, from roots to needles, and its conducting capacity are essential for the realization of systemic relationships and the control over growth and development in Pinus sylvestris L.     

The Variance of Identity-by-Descent Sharing in the Wright–Fisher Model

Widespread sharing of long, identical-by-descent (IBD) genetic segments is a hallmark of populations that have experienced recent genetic drift. Detection of these IBD segments has recently become feasible, enabling a wide range of applications from phasing and imputation to demographic inference. Here, we study the distribution of IBD sharing in the Wright–Fisher model. Specifically, using coalescent theory, we calculate the variance of the total sharing between random pairs of individuals. We then investigate the cohort-averaged sharing: the average total sharing between one individual and the rest of the cohort. We find that for large cohorts, the cohort-averaged sharing is distributed approximately normally. Surprisingly, the variance of this distribution does not vanish even for large cohorts, implying the existence of “hypersharing” individuals. The presence of such individuals has consequences for the design of sequencing studies, since, if they are selected for whole-genome sequencing, a larger fraction of the cohort can be subsequently imputed. We calculate the expected gain in power of imputation by IBD and subsequently in power to detect an association, when individuals are either randomly selected or specifically chosen to be the hypersharing individuals. Using our framework, we also compute the variance of an estimator of the population size that is based on the mean IBD sharing and the variance in the sharing between inbred siblings. Finally, we study IBD sharing in an admixture pulse model and show that in the Ashkenazi Jewish population the admixture fraction is correlated with the cohort-averaged sharing.IN isolated populations, even purported unrelated individuals often share genetic material that is identical-by-descent (IBD). Traditionally, the term IBD sharing referred to coancestry at a single site (or autozygosity, in the case of a diploid individual) and was widely investigated as a measure of the degree of inbreeding in a population (Hartl and Clark 2006). Recent years have brought dramatic increases in the quantity and density of available genetic data and, together with new computational tools, these data have enabled the detection of IBD sharing of entire genomic segments (see, e.g., Purcell et al. 2007; Kong et al. 2008; Albrechtsen et al. 2009; Gusev et al. 2009; Browning and Browning 2011; Carr et al. 2011; Brown et al. 2012). The availability of IBD detection tools that are efficient enough to detect shared segments in large cohorts has resulted in numerous applications, from demographic inference (Davison et al. 2009; Palamara et al. 2012) and characterization of populations (Gusev et al. 2012a) to selection detection (Albrechtsen et al. 2010), relatedness detection and pedigree reconstruction (Huff et al. 2011; Kirkpatrick et al. 2011; Stevens et al. 2011; Henn et al. 2012), prioritization of individuals for sequencing (Gusev et al. 2012b), inference of HLA type (Setty et al. 2011), detection of haplotypes associated with a disease or a trait (Akula et al. 2011; Gusev et al. 2011; Browning and Thompson 2012), imputation (Uricchio et al. 2012), and phasing (Palin et al. 2011).Recently, some of us used coalescent theory to calculate several theoretical quantities of IBD sharing under a number of demographic histories. Then, shared segments were detected in real populations, and their demographic histories were inferred (Palamara et al. 2012). Here, we expand upon Palamara et al. (2012) to investigate additional aspects of the stochastic variation in IBD sharing. Specifically, we provide a precise calculation for the variance of the total sharing in the Wright–Fisher model, either between a random pair of individuals or between one individual and all others in the cohort.Understanding the variation in IBD sharing is an important theoretical characterization of the Wright–Fisher model, and additionally, it has several practical applications. For example, it can be used to calculate the variance of an estimator of the population size that is based on the sharing between random pairs. In a different domain, the variance in IBD sharing is needed to accurately assess strategies for sequencing study design, specifically, in prioritization of individuals to be sequenced. This is because imputation strategies use IBD sharing between sequenced individuals and genotyped, not-sequenced individuals to increase the number of effective sequences analyzed in the association study (Palin et al. 2011; Gusev et al. 2012b; Uricchio et al. 2012).In the remainder of this article, we first review the derivation of the mean fraction of the genome shared between two individuals (Palamara et al. 2012). We then calculate the variance of this quantity, using coalescent theory with recombination. We provide a number of approximations, one of which results in a surprisingly simple expression, which is then generalized to a variable population size and to the sharing of segments in a length range. We also numerically investigate the pairwise sharing distribution and provide an approximate fit. We then turn to the average total sharing between each individual and the entire cohort. We show that this quantity, which we term the cohort-averaged sharing, is approximately normally distributed, but is much wider than naively expected, implying the existence of hypersharing individuals. We consider several applications: the number of individuals needed to be sequenced to achieve a certain imputation power and the implications to disease mapping, inference of the population size based on the total sharing, and the variance of the sharing between siblings. We finally calculate the mean and the variance of the sharing in an admixture pulse model and show numerically that admixture results in a broader than expected cohort-averaged sharing. Therefore, large variance of the cohort-averaged sharing can indicate admixture. In the Ashkenazi Jewish population, we show that the cohort-averaged sharing is strongly anticorrelated with the fraction of European ancestry.     

The time aspect of bioenergy – climate impacts of solid biofuels due to carbon dynamics

The climate impacts from bioenergy involve an important time aspect. Using forest residues for energy may result in high initial emissions, but net emissions are reduced over time since, if the residues were left on the ground, they would decompose and release CO₂ to the atmosphere. This article investigates the climate impacts from bioenergy with special focus on the time aspects. More specifically, we analyze the climate impacts of forest residues and stumps where combustion related emissions are compensated by avoided emissions from leaving them on the ground to decompose. These biofuels are compared with fossil gas and coal. Net emissions are defined as emissions from utilizing the fuel minus emissions from a reference case of no utilization. Climate impacts are estimated using the measures radiative forcing and global average surface temperature. We find that the climate impacts from using forest residues and stumps depend on the decomposition rates and the time perspective over which the analysis is done. Over a 100 year perspective, branches and tops have lower climate impacts than stumps which in turn have lower impacts than fossil gas and coal. Over a 20 year time perspective, branches and tops have lower climate impacts than all other fuels but the relative difference is smaller. However, stumps have slightly higher climate impacts over 20 years than fossil gas but lower impacts than coal. Regarding metrics for climate impacts, over shorter time scales, approximately 30 years or less, radiative forcing overestimates the climate impacts compared with impacts expressed by global surface temperature change, which is due to the inertia of the climate system. We also find that establishing willow on earlier crop land may reduce atmospheric CO₂, provided new land is available. However, these results are inconclusive since we haven't considered the effects of producing the agricultural crops elsewhere.     

The influence of ion concentrations on the dynamic behavior of the Hodgkin–Huxley model-based cortical network

Action potentials (APs) in the form of very short pulses arise when the cell is excited by any internal or external stimulus exceeding the critical threshold of the membrane. During AP generation, the membrane potential completes its natural cycle through typical phases that can be formatted by ion channels, gates and ion concentrations, as well as the synaptic excitation rate. On the basis of the Hodgkin–Huxley cell model, a cortical network consistent with the real anatomic structure is realized with randomly interrelated small population of neurons to simulate a cerebral cortex segment. Using this model, we investigated the effects of Na⁺ and K⁺ ion concentrations on the outcome of this network in terms of regularity, phase locking, and synchronization. The results suggested that Na⁺ concentration does slightly affect the amplitude but not considerably affects the other parameters specified by depolarization and repolarization. K⁺ concentration significantly influences the form, regularity, and synchrony of the network-generated APs. No previous study dealing directly with the effects of both Na⁺ and K⁺ ion concentrations on regularity and synchronization of the simulated cortical network-generated APs, allowing for the comparison of results obtained using our methods, was encountered in the literature. The results, however, were consistent with those obtained through studies concerning resonance and synchronization from another perspective and with the information revealed through physiological and pharmacological experiments concerning changing ion concentrations or blocking ion channels. Our results demonstrated that the regularity and reliability of brain functions have a strong relationship with cellular ion concentrations, and suggested the management of the dynamic behavior of the cellular network with ion concentrations.     

The twenty-fifth anniversary of the Limnological Institute,The Netherlands (1957–1982)

This paper marks the silver jubilee in 1982 of the Limnological Institute, The Netherlands. It describes the history, present research programme and organisation of the Institute, as well as its cooperative studies with other institutes and its training facilities.Some remarks about its future plans of work are made.