Hybrid stochastic simulations of intracellular reaction-diffusion systems.

With the observation that stochasticity is important in biological systems, chemical kinetics have begun to receive wider interest. While the use of Monte Carlo discrete event simulations most accurately capture the variability of molecular species, they become computationally costly for complex reaction-diffusion systems with large populations of molecules. On the other hand, continuous time models are computationally efficient but they fail to capture any variability in the molecular species. In this study a hybrid stochastic approach is introduced for simulating reaction-diffusion systems. We developed an adaptive partitioning strategy in which processes with high frequency are simulated with deterministic rate-based equations, and those with low frequency using the exact stochastic algorithm of Gillespie. Therefore the stochastic behavior of cellular pathways is preserved while being able to apply it to large populations of molecules. We describe our method and demonstrate its accuracy and efficiency compared with the Gillespie algorithm for two different systems. First, a model of intracellular viral kinetics with two steady states and second, a compartmental model of the postsynaptic spine head for studying the dynamics of Ca+2 and NMDA receptors.

Influence of transrectal and transabdominal ultrasound examination on salivary cortisol, heart rate, and heart rate variability in mares

Pregnancy diagnostics in equine reproduction are routinely performed using transrectal ultrasonography, although it is also possible to visualize the fetus by transabdominal ultrasound examinations from the 90th day of gestation onward. We hypothesized that ultrasound examinations may stress the mare and that the gestational stage status and lactation may influence the mare's stress reaction. To investigate the stress reaction, 25 thoroughbred mares of different age, pregnancy and lactational status underwent a transrectal examination. In pregnant mares, an additional transabdominal examination was performed. Salivary cortisol concentration, mean heart rate, and heart rate variability of mares were assessed to evaluate the reactions of hypothalamic–pituitary–adrenal (HPA) axis and of the autonomic nervous system. Significant differences were observed between lactating and nonlactating mares; with a lower responsiveness to stress in lactating mares. The transrectal ultrasound examination in nonlactating mares induced a significant increase in salivary cortisol (P < 0.05), and in the heart rate variability parameter, ratio of low to high frequencies (P < 0.05). This reflects an activation of the HPA axis and a shift to more sympathetic dominance. In contrast, a transabdominally performed pregnancy check did not induce an activation of the HPA axis over basal level but increased the mean heart rate and low to high frequency ratio. The results of this study indicate that checks of advanced pregnancies can be easily performed by transabdominal ultrasonography. With regard to animal welfare, this technique should be preferred during midgestation in nonlactating mares.

Probing the Electrocatalytic Oxygen Reduction Reaction Reactivity of Immobilized Multicopper Oxidase CueO

The bioelectrocatalytic (oxygen reduction reaction, ORR) properties of the multicopper oxidase CueO immobilized on gold electrodes were investigated. Macroscopic electrochemical techniques were combined with in situ scanning tunneling microscopy (STM) and surface-enhanced Raman spectroscopy at the ensemble and at the single-molecule level. Self-assembled monolayer of mercaptopropionic acid, cysteamine, and p-aminothiophenol were chosen as redox mediators. The highest ORR activity was observed for the protein attached to amino-terminated adlayers. In situ STM experiments revealed that the presence of oxygen causes distinct structure and electronic changes in the metallic centers of the enzyme, which determine the rate of intramolecular electron transfer and, consequently, affect the rate of electron tunneling through the protein. Complementary Raman spectroscopy experiments provided access for monitoring structural changes in the redox state of the type 1 copper center of the immobilized enzyme during the CueO-catalyzed oxygen reduction cycle. These results unequivocally demonstrate the existence of a direct electronic communication between the electrode substrate and the type 1 copper center.

Monetary Policy, Exchange Rate Overshooting, and Endogenous Physical Capital

We develop an open economy macroeconomic model with real capital accumulation and microeconomic foundations. We show that expansionary monetary policy causes exchange rate overshooting, not once, but potentially twice; the secondary repercussion comes through the reaction of firms to changed asset prices and the firms' decisions to invest in real capital. The model sheds further light on the volatility of real and nominal exchange rates, and it suggests that changes in corporate sector profitability may affect exchange rates through international portfolio diversification in corporate securities.

Impact of real time polymerase chain reaction technology on the therapeutic approach to treatment of staphylococcal infections

Can the early identification of the species of staphylococcus responsible for infection by the use of Real Time PCR technology influence the approach to the treatment of these infections? ^ This study was a retrospective cohort study in which two groups of patients were compared. The first group, ‘Physician Aware’ consisted of patients in whom physicians were informed of specific staphylococcal species and antibiotic sensitivity (using RT-PCR) at the time of notification of the gram stain. The second group, ‘Physician Unaware’ consisted of patients in whom treating physicians received the same information 24–72 hours later as a result of blood culture and antibiotic sensitivity determination. ^ The approach to treatment was compared between ‘Physician Aware’ and ‘Physician Unaware’ groups for three different microbiological diagnoses—namely MRSA, MSSA and no-SA (or coagulase negative Staphylococcus). ^ For a diagnosis of MRSA, the mean time interval to the initiation of Vancomycin therapy was 1.08 hours in the ‘Physician Aware’ group as compared to 5.84 hours in the ‘Physician Unaware’ group (p=0.34). ^ For a diagnosis of MSSA, the mean time interval to the initiation of specific anti-MSSA therapy with Nafcillin was 5.18 hours in the ‘Physician Aware’ group as compared to 49.8 hours in the ‘Physician Unaware’ group (p=0.007). Also, for the same diagnosis, the mean duration of empiric therapy in the ‘Physician Aware’ group was 19.68 hours as compared to 80.75 hours in the ‘Physician Unaware’ group (p=0.003) ^ For a diagnosis of no-SA or coagulase negative staphylococcus, the mean duration of empiric therapy was 35.65 hours in the ‘Physician Aware’ group as compared to 44.38 hours in the ‘Physician Unaware’ group (p=0.07). However, when treatment was considered a categorical variable and after exclusion of all cases where anti-MRS therapy was used for unrelated conditions, only 20 of 72 cases in the ‘Physician Aware’ group received treatment as compared to 48 of 106 cases in the ‘Physician Unaware’ group. ^ Conclusions. Earlier diagnosis of MRSA may not alter final treatment outcomes. However, earlier identification may lead to the earlier institution of measures to limit the spread of infection. The early diagnosis of MSSA infection, does lead to treatment with specific antibiotic therapy at an earlier stage of treatment. Also, the duration of empiric therapy is greatly reduced by early diagnosis. The early diagnosis of coagulase negative staphylococcal infection leads to a lower rate of unnecessary treatment for these infections as they are commonly considered contaminants. ^

Thermal reaction norms of growth in Atlantic and Pacific silverside fishes

How organisms may adapt to rising global temperatures is uncertain, but concepts can emerge from studying adaptive physiological trait variations across existing spatial climate gradients. Many ectotherms, particularly fish, have evolved increasing genetic growth capacities with latitude (i.e. countergradient variation (CnGV) in growth), which are thought to be an adaptation primarily to strong gradients in seasonality. In contrast, evolutionary responses to gradients in mean temperature are often assumed to involve an alternative mode, 'thermal adaptation'. We measured thermal growth reaction norms in Pacific silverside populations (Atherinops affinis) occurring across a weak latitudinal temperature gradient with invariant seasonality along the North American Pacific coast. Instead of thermal adaptation, we found novel evidence for CnGV in growth, suggesting that CnGV is a ubiquitous mode of reaction-norm evolution in ectotherms even in response to weak spatial and, by inference, temporal climate gradients. A novel, large-scale comparison between ecologically equivalent Pacific versus Atlantic silversides (Menidia menidia) revealed how closely growth CnGV patterns reflect their respective climate gradients. While steep growth reaction norms and increasing growth plasticity with latitude in M. menidia mimicked the strong, highly seasonal Atlantic coastal gradient, shallow reaction norms and much smaller, latitude-independent growth plasticity in A. affinis resembled the weak Pacific latitudinal temperature gradient.

Swift thermal reaction norm evolution in a key coccolithopore species: Reaction norm assay data from an experiment in Kiel from December 2012 to June 2015

Temperature has a profound effect on the species composition and physiology of marine phytoplankton, a polyphyletic group of microbes responsible for half of global primary production. Here, we ask whether and how thermal reaction norms in a key calcifying species, the coccolithophore Emiliania huxleyi, change as a result of 2.5 years of experimental evolution to a temperature about 2°C below its upper thermal limit. Replicate experimental populations derived from a single genotype isolated from Norwegian coastal waters were grown at two temperatures for 2.5 years before assessing thermal responses at 6 temperatures ranging from 15 to 26°C, with pCO2 (400/1100/2200 ?atm) as a fully factorial additional factor. The two selection temperatures (15°/26.3°C) led to a marked divergence of thermal reaction norms. Optimal growth temperatures were 0.7°C higher in experimental populations selected at 26.3°C than those selected at 15.0°C. An additional negative effect of high pCO2 on maximal growth rate (8% decrease relative to lowest level) was observed. Finally, the maximum persistence temperature (Tmax) differed by 1-3°C between experimental treatments, as a result of an interaction between pCO2 and the temperature selection. Taken together, we demonstrate that several attributes of thermal reaction norms in phytoplankton may change faster than the predicted progression of ocean warming.

An Empirical Analysis of the Monetary Policy Reaction Function in India

This paper empirically analyzes India’s monetary policy reaction function by applying the Taylor (1993) rule and its open-economy version which employs dynamic OLS. The analysis uses monthly data from the period of April 1998 to December 2007. When the simple Taylor rule was estimated for India, the output gap coefficient was statistically significant, and its sign condition was found to be consistent with theoretical rationale; however, the same was not true of the inflation coefficient. When the Taylor rule with exchange rate was estimated, the coefficients of output gap and exchange rate had statistical significance with the expected signs, whereas the results of inflation remained the same as before. Therefore, the inflation rate has not played a role in the conduct of India’s monetary policy, and it is inappropriate for India to adopt an inflation-target type policy framework.

On the reaction mechanism of l-lactate oxidase: Quantitative structure-activity analysis of the reaction with para-substituted l-mandelates

The rate constants for reduction of the flavoenzyme, l-lactate oxidase, and a mutant (in which alanine 95 is replaced by glycine), by a series of para-substituted mandelates, in both the 2-1H- and 2-2H- forms, have been measured by rapid reaction spectrophotometry. In all cases, significant isotope effects (1H/2H = 3–7) on the rate constants of flavin reduction were found, indicating that flavin reduction is a direct measure of α-C-H bond breakage. The rate constants show only a small influence of the electronic characteristics of the substituents, but show a good correlation when combined with some substituent volume parameters. A surprisingly good correlation is found with the molecular mass of the substrate. The results are compatible with any mechanism in which there is little development of charge in the transition state. This could be a transfer of hydride to the flavin N(5) position or a synchronous mechanism in which the α-C-H is formally abstracted as a H+ while the resulting charge is simultaneously neutralized by another event.

Photoassembly of the manganese cluster and oxygen evolution from monomeric and dimeric CP47 reaction center photosystem II complexes

Isolated subcomplexes of photosystem II from spinach (CP47RC), composed of D1, D2, cytochrome b559, CP47, and a number of hydrophobic small subunits but devoid of CP43 and the extrinsic proteins of the oxygen-evolving complex, were shown to reconstitute the Mn4Ca1Clx cluster of the water-splitting system and to evolve oxygen. The photoactivation process in CP47RC dimers proceeds by the same two-step mechanism as observed in PSII membranes and exhibits the same stoichiometry for Mn2+, but with a 10-fold lower affinity for Ca2+ and an increased susceptibility to photodamage. After the lower Ca2+ affinity and the 10-fold smaller absorption cross-section for photons in CP47 dimers is taken into account, the intrinsic rate constant for the rate-limiting calcium-dependent dark step is indistinguishable for the two systems. The monomeric form of CP47RC also showed capacity to photoactivate and catalyze water oxidation, but with lower activity than the dimeric form and increased susceptibility to photodamage. After optimization of the various parameters affecting the photoactivation process in dimeric CP47RC subcores, 18% of the complexes were functionally reconstituted and the quantum efficiency for oxygen production by reactivated centers approached 96% of that observed for reconstituted photosystem II-enriched membranes.

Ab initio calculation of electronic coupling in the photosynthetic reaction center

We have carried out an ab initio electronic structure calculations of electron transfer couplings between chromophores in the bacterial photosynthetic reaction center. The couplings agree remarkably well with parameters obtained from recent quantum dynamical modeling of experimental data assuming an explicit intermediate mechanism. We also have computed couplings on the M-side of the reaction center and have found that the interaction of the primary donor to the M-side intermediate bacteriochlorophyll is quite small because of destructive interference of the two localized coupling matrix elements. This may explain the slow rate of electron transfer down the M-side of the reaction center.

Determination of the binding sites of the proton transfer inhibitors Cd2+ and Zn2+ in bacterial reaction centers

The reaction center (RC) from Rhodobacter sphaeroides couples light-driven electron transfer to protonation of a bound quinone acceptor molecule, QB, within the RC. The binding of Cd2+ or Zn2+ has been previously shown to inhibit the rate of reduction and protonation of QB. We report here on the metal binding site, determined by x-ray diffraction at 2.5-Å resolution, obtained from RC crystals that were soaked in the presence of the metal. The structures were refined to R factors of 23% and 24% for the Cd2+ and Zn2+ complexes, respectively. Both metals bind to the same location, coordinating to Asp-H124, His-H126, and His-H128. The rate of electron transfer from QA− to QB was measured in the Cd2+-soaked crystal and found to be the same as in solution in the presence of Cd2+. In addition to the changes in the kinetics, a structural effect of Cd2+ on Glu-H173 was observed. This residue was well resolved in the x-ray structure—i.e., ordered—with Cd2+ bound to the RC, in contrast to its disordered state in the absence of Cd2+, which suggests that the mobility of Glu-H173 plays an important role in the rate of reduction of QB. The position of the Cd2+ and Zn2+ localizes the proton entry into the RC near Asp-H124, His-H126, and His-H128. Based on the location of the metal, likely pathways of proton transfer from the aqueous surface to QB⨪ are proposed.

Identification of the proton pathway in bacterial reaction centers: Inhibition of proton transfer by binding of Zn2+ or Cd2+

The reaction center (RC) from Rhodobacter sphaeroides converts light into chemical energy through the light induced two-electron, two-proton reduction of a bound quinone molecule QB (the secondary quinone acceptor). A unique pathway for proton transfer to the QB site had so far not been determined. To study the molecular basis for proton transfer, we investigated the effects of exogenous metal ion binding on the kinetics of the proton-assisted electron transfer kAB(2) (QA−•QB−• + H+ → QA(QBH)−, where QA is the primary quinone acceptor). Zn2+ and Cd2+ bound stoichiometrically to the RC (KD ≤ 0.5 μM) and reduced the observed value of kAB(2) 10-fold and 20-fold (pH 8.0), respectively. The bound metal changed the mechanism of the kAB(2) reaction. In native RCs, kAB(2) was previously shown to be rate-limited by electron transfer based on the dependence of kAB(2) on the driving force for electron transfer. Upon addition of Zn2+ or Cd2+, kAB(2) became approximately independent of the electron driving force, implying that the rate of proton transfer was reduced (≥ 102-fold) and has become the rate-limiting step. The lack of an effect of the metal binding on the charge recombination reaction D+•QAQB−• → DQAQB suggests that the binding site is located far (>10 Å) from QB. This hypothesis is confirmed by preliminary x-ray structure analysis. The large change in the rate of proton transfer caused by the stoichiometric binding of the metal ion shows that there is one dominant site of proton entry into the RC from which proton transfer to QB−• occurs.

Identification of the proton pathway in bacterial reaction centers: Both protons associated with reduction of QB to QBH2 share a common entry point

The reaction center from Rhodobacter sphaeroides uses light energy for the reduction and protonation of a quinone molecule, QB. This process involves the transfer of two protons from the aqueous solution to the protein-bound QB molecule. The second proton, H+(2), is supplied to QB by Glu-L212, an internal residue protonated in response to formation of QA− and QB−. In this work, the pathway for H+(2) to Glu-L212 was studied by measuring the effects of divalent metal ion binding on the protonation of Glu-L212, which was assayed by two types of processes. One was proton uptake from solution after the one-electron reduction of QA (DQA→D+QA−) and QB (DQB→D+QB−), studied by using pH-sensitive dyes. The other was the electron transfer kAB(1) (QA−QB→QAQB−). At pH 8.5, binding of Zn2+, Cd2+, or Ni2+ reduced the rates of proton uptake upon QA− and QB− formation as well as kAB(1) by ≈an order of magnitude, resulting in similar final values, indicating that there is a common rate-limiting step. Because D+QA− is formed 105-fold faster than the induced proton uptake, the observed rate decrease must be caused by an inhibition of the proton transfer. The Glu-L212→Gln mutant reaction centers displayed greatly reduced amplitudes of proton uptake and exhibited no changes in rates of proton uptake or electron transfer upon Zn2+ binding. Therefore, metal binding specifically decreased the rate of proton transfer to Glu-L212, because the observed rates were decreased only when proton uptake by Glu-L212 was required. The entry point for the second proton H+(2) was thus identified to be the same as for the first proton H+(1), close to the metal binding region Asp-H124, His-H126, and His-H128.