Modeling Two-Oscillator Circadian Systems Entrained by Two Environmental Cycles

Several experimental studies have altered the phase relationship between photic and non-photic environmental, 24 h cycles (zeitgebers) in order to assess their role in the synchronization of circadian rhythms. To assist in the interpretation of the complex activity patterns that emerge from these ""conflicting zeitgeber'' protocols, we present computer simulations of coupled circadian oscillators forced by two independent zeitgebers. This circadian system configuration was first employed by Pittendrigh and Bruce (1959), to model their studies of the light and temperature entrainment of the eclosion oscillator in Drosophila. Whereas most of the recent experiments have restricted conflicting zeitgeber experiments to two experimental conditions, by comparing circadian oscillator phases under two distinct phase relationships between zeitgebers (usually 0 and 12 h), Pittendrigh and Bruce compared eclosion phase under 12 distinct phase relationships, spanning the 24 h interval. Our simulations using non-linear differential equations replicated complex non-linear phenomena, such as ""phase jumps'' and sudden switches in zeitgeber preferences, which had previously been difficult to interpret. Our simulations reveal that these phenomena generally arise when inter-oscillator coupling is high in relation to the zeitgeber strength. Manipulations in the structural symmetry of the model indicated that these results can be expected to apply to a wide range of system configurations. Finally, our studies recommend the use of the complete protocol employed by Pittendrigh and Bruce, because different system configurations can generate similar results when a ""conflicting zeitgeber experiment'' incorporates only two phase relationships between zeitgebers.

Aerosol number fluxes over the Amazon rain forest during the wet season

Number fluxes of particles with diameter larger than 10 nm were measured with the eddy covariance method over the Amazon rain forest during the wet season as part of the LBA (The Large Scale Biosphere Atmosphere Experiment in Amazonia) campaign 2008. The primary goal was to investigate whether sources or sinks dominate the aerosol number flux in the tropical rain forest-atmosphere system. During the measurement campaign, from 12 March to 18 May, 60% of the particle fluxes pointed downward, which is a similar fraction to what has been observed over boreal forests. The net deposition flux prevailed even in the absolute cleanest atmospheric conditions during the campaign and therefore cannot be explained only by deposition of anthropogenic particles. The particle transfer velocity v(t) increased with increasing friction velocity and the relation is described by the equation v(t) = 2.4x10(-3)xu(*) where u(*) is the friction velocity. Upward particle fluxes often appeared in the morning hours and seem to a large extent to be an effect of entrainment fluxes into a growing mixed layer rather than primary aerosol emission. In general, the number source of primary aerosol particles within the footprint area of the measurements was small, possibly because the measured particle number fluxes reflect mostly particles less than approximately 200 nm. This is an indication that the contribution of primary biogenic aerosol particles to the aerosol population in the Amazon boundary layer may be low in terms of number concentrations. However, the possibility of horizontal variations in primary aerosol emission over the Amazon rain forest cannot be ruled out.

Nonaxisymmetric magnetorotational instability in ideal and viscous plasmas

The excitation of magnetorotational instability (MRI) in rotating laboratory plasmas is investigated. In contrast to astrophysical plasmas, in which gravitation plays an important role, in laboratory plasmas it can be neglected and the plasma rotation is equilibrated by the pressure gradient. The analysis is restricted to the simple model of a magnetic confinement configuration with cylindrical symmetry, in which nonaxisymmetric perturbations are investigated using the local approximation. Starting from the simplest case of an ideal plasma, the corresponding dispersion relations are derived for more complicated models including the physical effects of parallel and perpendicular viscosities. The Friemann-Rotenberg approach used for ideal plasmas is generalized for the viscous model and an analytical expression for the instability boundary is obtained. It is shown that, in addition to the standard effect of radial derivative of the rotation frequency (the Velikhov effect), which can be destabilizing or stabilizing depending on the sign of this derivative in the ideal plasma, there is a destabilizing effect proportional to the fourth power of the rotation frequency, or, what is the same, to the square of the plasma pressure gradient, and to the square of the azimuthal mode number of the perturbations. It is shown that the instability boundary also depends on the product of the plasma pressure and density gradients, which has a destabilizing effect when it is negative. In the case of parallel viscosity, the MRI looks like an ideal instability independent of viscosity, while, in the case of strong perpendicular viscosity, it is a dissipative instability with the growth rate inversely proportional to the characteristic viscous decay rate. We point out, however, that the modes of the continuous range of the magnetohydrodynamics spectrum are not taken into account in this paper, and they can be more dangerous than those that are considered. (c) 2008 American Institute of Physics.

Nonlocal magnetorotational instability

An analytical theory of the nonlocal magnetorotational instability (MRI) is developed for the simplest astrophysical plasma model. It is assumed that the rotation frequency profile has a steplike character, so that there are two regions in which it has constant different values, separated by a narrow transition layer. The surface wave approach is employed to investigate the MRI in this configuration. It is shown that the main regularities of the nonlocal MRI are similar to those of the local instability and that driving the nonaxisymmetric MRI is less effective than the axisymmetric one, also for the case of the nonlocal instability. The existence of nonlocal instabilities in nonmagnetized plasma is predicted. (c) 2008 American Institute of Physics.

The Importance of Protein-Protein Interactions on the pH-Induced Conformational Changes of Bovine Serum Albumin: A Small-Angle X-Ray Scattering Study

The combined effects of concentration and pH on the conformational states of bovine serum albumin (BSA) are investigated by small-angle x-ray scattering. Serum albumins, at physiological conditions, are found at concentrations of similar to 35-45 mg/mL (42 mg/mL in the case of humans). In this work, BSA at three different concentrations (10, 25, and 50 mg/mL) and pH values (2.0-9.0) have been studied. Data were analyzed by means of the Global Fitting procedure, with the protein form factor calculated from human serum albumin (HSA) crystallographic structure and the interference function described, considering repulsive and attractive interaction potentials within a random phase approximation. Small-angle x-ray scattering data show that BSA maintains its native state from pH 4.0 up to 9.0 at all investigated concentrations. A pH-dependence of the absolute net protein charge is shown and the charge number per BSA is quantified to 10(2), 8(l), 13(2), 20(2), and 26(2) for pH values 4.0, 5.4, 7.0, 8.0, and 9.0, respectively. The attractive potential diminishes as BSA concentration increases. The coexistence of monomers and dimers is observed at 50 mg/mL and pH 5.4, near the BSA isoelectric point. Samples at pH 2.0 show a different behavior, because BSA overall shape changes as a function of concentration. At 10 mg/mL, BSA is partially unfolded and a strong repulsive protein-protein interaction occurs due to the high amount of exposed charge. At 25 and 50 mg/mL, BSA undergoes some refolding, which likely results in a molten-globule state. This work concludes by confirming that the protein concentration plays an important role on the pH-unfolded BSA state, due to a delicate compromise between interaction forces and crowding effects.

Thermodynamic origin of cooperativity in actomyosin interactions: The coupling of short-range interactions with actin bending stiffness in an Ising-like model

We present Monte Carlo simulations for a molecular motor system found in virtually all eukaryotic cells, the acto-myosin motor system, composed of a group of organic macromolecules. Cell motors were mapped to an Ising-like model, where the interaction field is transmitted through a tropomyosin polymer chain. The presence of Ca(2+) induces tropomyosin to block or unblock binding sites of the myosin motor leading to its activation or deactivation. We used the Metropolis algorithm to find the transient and the equilibrium states of the acto-myosin system composed of solvent, actin, tropomyosin, troponin, Ca(2+), and myosin-S1 at a given temperature, including the spatial configuration of tropomyosin on the actin filament surface. Our model describes the short- and long-range cooperativity during actin-myosin binding which emerges from the bending stiffness of the tropomyosin complex. We found all transition rates between the states only using the interaction energy of the constituents. The agreement between our model and experimental data also supports the recent theory of flexible tropomyosin.

Combined Monte Carlo and quantum mechanics study of the solvatochromism of phenol in water. The origin of the blue shift of the lowest pi-pi* transition

A combined and sequential use of Monte Carlo simulations and quantum mechanical calculations is made to analyze the spectral shift of the lowest pi-pi* transition of phenol in water. The solute polarization is included using electrostatic embedded calculations at the MP2/aug-cc-pVDZ level giving a dipole moment of 2.25 D, corresponding to an increase of 76% compared to the calculated gas-phase value. Using statistically uncorrelated configurations sampled from the MC simulation,first-principle size-extensive calculations are performed to obtain the solvatochromic shift. Analysis is then made of the origin of the blue shift. Results both at the optimized geometry and in room-temperature liquid water show that hydrogen bonds of water with phenol promote a red shift when phenol is the proton-donor and a blue shift when phenol is the proton-acceptor. In the case of the optimized clusters the calculated shifts are in very good agreement with results obtained from mass-selected free jet expansion experiments. In the liquid case the contribution of the solute-solvent hydrogen bonds partially cancels and the total shift obtained is dominated by the contribution of the outer solvent water molecules. Our best result, including both inner and outer water molecules, is 570 +/- 35 cm(-1), in very good agreement with the small experimental shift of 460 cm(-1) for the absorption maximum.

Finding Quasi-Optimal Network Topologies for Information Transmission in Active Networks

This work clarifies the relation between network circuit (topology) and behaviour (information transmission and synchronization) in active networks, e.g. neural networks. As an application, we show how one can find network topologies that are able to transmit a large amount of information, possess a large number of communication channels, and are robust under large variations of the network coupling configuration. This theoretical approach is general and does not depend on the particular dynamic of the elements forming the network, since the network topology can be determined by finding a Laplacian matrix (the matrix that describes the connections and the coupling strengths among the elements) whose eigenvalues satisfy some special conditions. To illustrate our ideas and theoretical approaches, we use neural networks of electrically connected chaotic Hindmarsh-Rose neurons.