Memory from the dynamics of intrinsic membrane currents

Almost all theoretical and experimental studies of the mechanisms underlying learning and memory focus on synaptic efficacy and make the implicit assumption that changes in synaptic efficacy are both necessary and sufficient to account for learning and memory. However, network dynamics depends on the complex interaction between intrinsic membrane properties and synaptic strengths and time courses. Furthermore, neuronal activity itself modifies not only synaptic efficacy but also the intrinsic membrane properties of neurons. This paper presents examples demonstrating that neurons with complex temporal dynamics can provide short-term “memory” mechanisms that rely solely on intrinsic neuronal properties. Additionally, we discuss the potential role that activity may play in long-term modification of intrinsic neuronal properties. While not replacing synaptic plasticity as a powerful learning mechanism, these examples suggest that memory in networks results from an ongoing interplay between changes in synaptic efficacy and intrinsic membrane properties.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

Evolving noisy oscillatory dynamics in genetic regulatory networks

We introduce a genetic programming (GP) approach for evolving genetic networks that demonstrate desired dynamics when simulated as a discrete stochastic process. Our representation of genetic networks is based on a biochemical reaction model including key elements such as transcription, translation and post-translational modifications. The stochastic, reaction-based GP system is similar but not identical with algorithmic chemistries. We evolved genetic networks with noisy oscillatory dynamics. The results show the practicality of evolving particular dynamics in gene regulatory networks when modelled with intrinsic noise.

Modified beam theories for bending properties of nanowires considering surface/intrinsic effects and axial extension effect

Several studies of the surface effect on bending properties of a nanowire (NW) have been conducted. However, these analyses are mainly based on theoretical predictions, and there is seldom integration study in combination between theoretical predictions and simulation results. Thus, based on the molecular dynamics (MD) simulation and different modified beam theories, a comprehensive theoretical and numerical study for bending properties of nanowires considering surface/intrinsic stress effects and axial extension effect is conducted in this work. The discussion begins from the Euler-Bernoulli beam theory and Timoshenko beam theory augmented with surface effect. It is found that when the NW possesses a relatively small cross-sectional size, these two theories cannot accurately interpret the true surface effect. The incorporation of axial extension effect into Euler-Bernoulli beam theory provides a nonlinear solution that agrees with the nonlinear-elastic experimental and MD results. However, it is still found inaccurate when the NW cross-sectional size is relatively small. Such inaccuracy is also observed for the Euler-Bernoulli beam theory augmented with both contributions from surface effect and axial extension effect. A comprehensive model for completely considering influences from surface stress, intrinsic stress, and axial extension is then proposed, which leads to good agreement with MD simulation results. It is thus concluded that, for NWs with a relatively small cross-sectional size, a simple consideration of surface stress effect is inappropriate, and a comprehensive consideration of the intrinsic stress effect is required.

Intrinsic kinetics of titania photocatalysis : simplified models for their investigation

Even though titanium dioxide photocatalysis has been promoted as a leading green technology for water purification, many issues have hindered its application on a large commercial scale. For the materials scientist the main issues have centred the synthesis of more efficient materials and the investigation of degradation mechanisms; whereas for the engineers the main issues have been the development of appropriate models and the evaluation of intrinsic kinetics parameters that allow the scale up or re-design of efficient large-scale photocatalytic reactors. In order to obtain intrinsic kinetics parameters the reaction must be analysed and modelled considering the influence of the radiation field, pollutant concentrations and fluid dynamics. In this way, the obtained kinetic parameters are independent of the reactor size and configuration and can be subsequently used for scale-up purposes or for the development of entirely new reactor designs. This work investigates the intrinsic kinetics of phenol degradation over titania film due to the practicality of a fixed film configuration over a slurry. A flat plate reactor was designed in order to be able to control reaction parameters that include the UV irradiance, flow rates, pollutant concentration and temperature. Particular attention was paid to the investigation of the radiation field over the reactive surface and to the issue of mass transfer limited reactions. The ability of different emission models to describe the radiation field was investigated and compared to actinometric measurements. The RAD-LSI model was found to give the best predictions over the conditions tested. Mass transfer issues often limit fixed film reactors. The influence of this phenomenon was investigated with specifically planned sets of benzoic acid experiments and with the adoption of the stagnant film model. The phenol mass transfer coefficient in the system was calculated to be km,phenol=8.5815x10-7Re0.65(ms-1). The data obtained from a wide range of experimental conditions, together with an appropriate model of the system, has enabled determination of intrinsic kinetic parameters. The experiments were performed in four different irradiation levels (70.7, 57.9, 37.1 and 20.4 W m-2) and combined with three different initial phenol concentrations (20, 40 and 80 ppm) to give a wide range of final pollutant conversions (from 22% to 85%). The simple model adopted was able to fit the wide range of conditions with only four kinetic parameters; two reaction rate constants (one for phenol and one for the family of intermediates) and their corresponding adsorption constants. The intrinsic kinetic parameters values were defined as kph = 0.5226 mmol m-1 s-1 W-1, kI = 0.120 mmol m-1 s-1 W-1, Kph = 8.5 x 10-4 m3 mmol-1 and KI = 2.2 x 10-3 m3 mmol-1. The flat plate reactor allowed the investigation of the reaction under two different light configurations; liquid and substrate side illumination. The latter of particular interest for real world applications where light absorption due to turbidity and pollutants contained in the water stream to be treated could represent a significant issue. The two light configurations allowed the investigation of the effects of film thickness and the determination of the catalyst optimal thickness. The experimental investigation confirmed the predictions of a porous medium model developed to investigate the influence of diffusion, advection and photocatalytic phenomena inside the porous titania film, with the optimal thickness value individuated at 5 ìm. The model used the intrinsic kinetic parameters obtained from the flat plate reactor to predict the influence of thickness and transport phenomena on the final observed phenol conversion without using any correction factor; the excellent match between predictions and experimental results provided further proof of the quality of the parameters obtained with the proposed method.

Base-induced decomposition of alkyl hydroperoxides in the gas phase. Part 3. Kinetics and dynamics in HO-+CH3OOH, C2H5OOH, and tert-C4H9OOH reactions

The E-CO(2) elimination reactions of alkyl hydroperoxides proceed via abstraction of an (x-hydrogen by a base: X- + (RRHCOOH)-R-1-H-2 -> HX + (RRC)-R-1-C-2=O + HO-. Efficiencies and product distributions for the reactions of the hydroxide anion with methyl, ethyl, and tert-butyl hydroperoxides are studied in the gas phase. On the basis of experiments using three isotopic analogues, HO- + CH3OOH, HO- + CD3OOH, and H18O- + CH3OOH. the overall intrinsic reaction efficiency is determined to be 80% or greater. The E(CO)2 decomposition is facile for these methylperoxide reactions, and predominates over competing proton transfer at the hydroperoxide moiety. The CH3CH2OOH reaction displays a similar E(CO)2 reactivity, whereas proton transfer and the formation of HOO- are the exclusive pathways observed for (CH3)(3)COOH, which has no (x-hydrogen. All results are consistent with the E-CO(2) mechanism, transition state structure, and reaction energy diagrams calculated using the hybrid density functional B3LYP approach. Isotope labeling for HO- + CH3OOH also reveals some interaction between H2O and HO- within the E(CO)2 product complex [H2O center dot center dot center dot CH2=O center dot center dot center dot HO-]. There is little evidence, however. for the formation of the most exothermic products H2O + CH2(OH)O-, which would arise from nuclephilic condensation of CH2=O and HO-. The results suggest that the product dynamics are not totally statistical but are rather direct after the E-CO(2) transition state. The larger HO- + CH3CH2OOH system displays more statistical behavior during complex dissociation.

Complex dynamics in polymer nanocomposites

Polymer nanocomposites offer the potential to create a new type of hybrid material with unique thermal, optical, or electrical properties. Understanding their structure, phase behavior, and dynamics is crucial for realizing such potentials. In this work we provide an experimental insight into the dynamics of such composites in terms of the temperature, wave vector, and volume fraction of nanoparticles, using multispeckle synchrotron x-ray photon correlation spectroscopy measurements on gold nanoparticles embedded in polymethylmethacrylate. Detailed analysis of the intermediate scattering functions reveals possible existence of an intrinsic length scale for dynamic heterogeneity in polymer nanocomposites similar to that seen in other soft materials like colloidal gels and glasses.

Stability of Dimeric Interface in Banana Lectin: Insight from Molecular Dynamics Simulations

Banana lectin (Banlec) is a homodimeric non-glycosylated protein. It exhibits the b-prism I structure. High-temperature molecular dynamics simulations have been utilized to monitor and understand early stages of thermally induced unfolding of Banlec. The present study elucidates the behavior of the dimeric protein at four different temperatures and compares the structural and conformational changes to that of the minimized crystal structure. The process of unfolding was monitored by following the radius of gyration, the rms deviation of each residue, change in relative solvent accessibility and the pattern of inter- and intra-subunit interactions. The overall study demonstrates that the Banlec dimer is a highly stable structure, and the stability is mostly contributed by interfacial interactions. It maintains its overall conformation during high-temperature (400–500 K) simulations, with only the unstructured loop regions acquiring greater momentum under such condition. Nevertheless, at still higher temperatures (600 K) the tertiary structure is gradually lost which later extends to loss of secondary structural elements. The pattern of hydrogen bonding within the subunit and at the interface across different stages has been analyzed and has provided rationale for its intrinsic high stability.

Role of conformational dynamics in kinetics of an enzymatic cycle in a nonequilibrium steady state

Enzyme is a dynamic entity with diverse time scales, ranging from picoseconds to seconds or even longer. Here we develop a rate theory for enzyme catalysis that includes conformational dynamics as cycling on a two-dimensional (2D) reaction free energy surface involving an intrinsic reaction coordinate (X) and an enzyme conformational coordinate (Q). The validity of Michaelis-Menten (MM) equation, i.e., substrate concentration dependence of enzymatic velocity, is examined under a nonequilibrium steady state. Under certain conditions, the classic MM equation holds but with generalized microscopic interpretations of kinetic parameters. However, under other conditions, our rate theory predicts either positive (sigmoidal-like) or negative (biphasic-like) kinetic cooperativity due to the modified effective 2D reaction pathway on X-Q surface, which can explain non-MM dependence previously observed on many monomeric enzymes that involve slow or hysteretic conformational transitions. Furthermore, we find that a slow conformational relaxation during product release could retain the enzyme in a favorable configuration, such that enzymatic turnover is dynamically accelerated at high substrate concentrations. The effect of such conformation retainment in a nonequilibrium steady state is evaluated.

Conceptual and statistical modelling of environmental effects in population dynamics

Population dynamics are generally viewed as the result of intrinsic (purely density dependent) and extrinsic (environmental) processes. Both components, and potential interactions between those two, have to be modelled in order to understand and predict dynamics of natural populations; a topic that is of great importance in population management and conservation. This thesis focuses on modelling environmental effects in population dynamics and how effects of potentially relevant environmental variables can be statistically identified and quantified from time series data. Chapter I presents some useful models of multiplicative environmental effects for unstructured density dependent populations. The presented models can be written as standard multiple regression models that are easy to fit to data. Chapters II IV constitute empirical studies that statistically model environmental effects on population dynamics of several migratory bird species with different life history characteristics and migration strategies. In Chapter II, spruce cone crops are found to have a strong positive effect on the population growth of the great spotted woodpecker (Dendrocopos major), while cone crops of pine another important food resource for the species do not effectively explain population growth. The study compares rate- and ratio-dependent effects of cone availability, using state-space models that distinguish between process and observation error in the time series data. Chapter III shows how drought, in combination with settling behaviour during migration, produces asymmetric spatially synchronous patterns of population dynamics in North American ducks (genus Anas). Chapter IV investigates the dynamics of a Finnish population of skylark (Alauda arvensis), and point out effects of rainfall and habitat quality on population growth. Because the skylark time series and some of the environmental variables included show strong positive autocorrelation, the statistical significances are calculated using a Monte Carlo method, where random autocorrelated time series are generated. Chapter V is a simulation-based study, showing that ignoring observation error in analyses of population time series data can bias the estimated effects and measures of uncertainty, if the environmental variables are autocorrelated. It is concluded that the use of state-space models is an effective way to reach more accurate results. In summary, there are several biological assumptions and methodological issues that can affect the inferential outcome when estimating environmental effects from time series data, and that therefore need special attention. The functional form of the environmental effects and potential interactions between environment and population density are important to deal with. Other issues that should be considered are assumptions about density dependent regulation, modelling potential observation error, and when needed, accounting for spatial and/or temporal autocorrelation.

A multistep transversal linearization (MTL) method in non-linear dynamics through a Magnus characterization

A new form of a multi-step transversal linearization (MTL) method is developed and numerically explored in this study for a numeric-analytical integration of non-linear dynamical systems under deterministic excitations. As with other transversal linearization methods, the present version also requires that the linearized solution manifold transversally intersects the non-linear solution manifold at a chosen set of points or cross-section in the state space. However, a major point of departure of the present method is that it has the flexibility of treating non-linear damping and stiffness terms of the original system as damping and stiffness terms in the transversally linearized system, even though these linearized terms become explicit functions of time. From this perspective, the present development is closely related to the popular practice of tangent-space linearization adopted in finite element (FE) based solutions of non-linear problems in structural dynamics. The only difference is that the MTL method would require construction of transversal system matrices in lieu of the tangent system matrices needed within an FE framework. The resulting time-varying linearized system matrix is then treated as a Lie element using Magnus’ characterization [W. Magnus, On the exponential solution of differential equations for a linear operator, Commun. Pure Appl. Math., VII (1954) 649–673] and the associated fundamental solution matrix (FSM) is obtained through repeated Lie-bracket operations (or nested commutators). An advantage of this approach is that the underlying exponential transformation could preserve certain intrinsic structural properties of the solution of the non-linear problem. Yet another advantage of the transversal linearization lies in the non-unique representation of the linearized vector field – an aspect that has been specifically exploited in this study to enhance the spectral stability of the proposed family of methods and thus contain the temporal propagation of local errors. A simple analysis of the formal orders of accuracy is provided within a finite dimensional framework. Only a limited numerical exploration of the method is presently provided for a couple of popularly known non-linear oscillators, viz. a hardening Duffing oscillator, which has a non-linear stiffness term, and the van der Pol oscillator, which is self-excited and has a non-linear damping term.

Nonlinear intrinsic instability of solid propellant combustion including gas-phase thermal inertia

The problem of homogeneous solid propellant combustion instability is studied with a one-dimensional flame model, including the effects of gas-phase thermal inertia and nonlinearity. Computational results presented in this paper show nonlinear instabilities inherent in the equations, due to which periodic burning is found even under steady ambient conditions such as pressure. The stability boundary is obtained in terms of Denison-Baum parameters. It is found that inclusion of gas-phase thermal inertia stabilizes the combustion. Also, the effect of a distributed heat release in the gas phase, compared to the flame sheet model, is to destabilize the burning. Direct calculations for finite amplitude pressure disturbances show that two distinct resonant modes exist, the first one near the natural frequency as obtained from intrinsic instability analysis and a second mode occurring at a much higher driving frequency. It is found that er rn in the low frequency region, the response of the propellant is significantly affected by the specific type of gas-phase chemical heat-release model employed. Examination of frequency response function reveals that the role of gas-phase thermal inertia is to stabilize the burning near the first resonant mode. Calculations made for different amplitudes of driving pressure show that the mean burning rate decreases with increasing amplitude. Also, with an increase in the driving amplitude, higher harmonics are generated in the burning rate.

Dynamics of Fractional Brownian Walks

We investigate the dynamics of polymers whose solution configurations are represented by fractional Brownian walks. The calculation of the two dynamical quantities considered here, the longest relaxation time tau(r) and the intrinsic viscosity [eta], is formulated in terms of Langevin equations and is carried out within the continuum approach developed in an earlier paper. Our results for tau(r) and [eta] reproduce known scaling relations and provide reasonable numerical estimates of scaling amplitudes. The possible relevance of the work to the study of globular proteins and other compact polymeric phases is discussed.

Reaction dynamics under confinement: an exact path integral treatment of a two-stage model of stochastic gene expression

Gene expression in living systems is inherently stochastic, and tends to produce varying numbers of proteins over repeated cycles of transcription and translation. In this paper, an expression is derived for the steady-state protein number distribution starting from a two-stage kinetic model of the gene expression process involving p proteins and r mRNAs. The derivation is based on an exact path integral evaluation of the joint distribution, P(p, r, t), of p and r at time t, which can be expressed in terms of the coupled Langevin equations for p and r that represent the two-stage model in continuum form. The steady-state distribution of p alone, P(p), is obtained from P(p, r, t) (a bivariate Gaussian) by integrating out the r degrees of freedom and taking the limit t -> infinity. P(p) is found to be proportional to the product of a Gaussian and a complementary error function. It provides a generally satisfactory fit to simulation data on the same two-stage process when the translational efficiency (a measure of intrinsic noise levels in the system) is relatively low; it is less successful as a model of the data when the translational efficiency (and noise levels) are high.

Zinc oxide based photocatalysis: tailoring surface-bulk structure and related interfacial charge carrier dynamics for better environmental applications

As an alternative to the gold standard TiO2 photocatalyst, the use of zinc oxide (ZnO) as a robust candidate for wastewater treatment is widespread due to its similarity in charge carrier dynamics upon bandgap excitation and the generation of reactive oxygen species in aqueous suspensions with TiO2. However, the large bandgap of ZnO, the massive charge carrier recombination, and the photoinduced corrosion-dissolution at extreme pH conditions, together with the formation of inert Zn(OH)(2) during photocatalytic reactions act as barriers for its extensive applicability. To this end, research has been intensified to improve the performance of ZnO by tailoring its surface-bulk structure and by altering its photogenerated charge transfer pathways with an intention to inhibit the surface-bulk charge carrier recombination. For the first time, the several strategies, such as tailoring the intrinsic defects, surface modification with organic compounds, doping with foreign ions, noble metal deposition, heterostructuring with other semiconductors and modification with carbon nanostructures, which have been successfully employed to improve the photoactivity and stability of ZnO are critically reviewed. Such modifications enhance the charge separation and facilitate the generation of reactive oxygenated free radicals, and also the interaction with the pollutant molecules. The synthetic route to obtain hierarchical nanostructured morphologies and study their impact on the photocatalytic performance is explained by considering the morphological influence and the defect-rich chemistry of ZnO. Finally, the crystal facet engineering of polar and non-polar facets and their relevance in photocatalysis is outlined. It is with this intention that the present review directs the further design, tailoring and tuning of the physico-chemical and optoelectronic properties of ZnO for better applications, ranging from photocatalysis to photovoltaics.