Computer simulation study of collective dynamics in the glass former Ca(NO3)(2)center dot 4H(2)O

Time correlation functions of current fluctuations were calculated by molecular dynamics (MD) simulations in order to investigate sound waves of high wavevectors in the glass-forming liquid Ca(NO3)(2)center dot 4H(2)O. Dispersion curves, omega(k), were obtained for longitudinal (LA) and transverse acoustic (TA) modes, and also for longitudinal optic (LO) modes. Spectra of LA modes calculated by MD simulations were modeled by a viscoelastic model within the memory function framework. The viscoelastic model is used to rationalize the change of slope taking place at k similar to 0.3 angstrom(-1) in the omega(k) curve of acoustic modes. For still larger wavevectors, mixing of acoustic and optic modes is observed. Partial time correlation functions of longitudinal mass currents were calculated separately for the ions and the water molecules. The wavevector dependence of excitation energies of the corresponding partial LA modes indicates the coexistence of a relatively stiff subsystem made of cations and anions, and a softer subsystem made of water molecules. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4751548]

A general framework for multi-compartmental analysis of drug chemotherapy dynamics in human immunodeficiency virus type-1 infected individuals

An optimal control strategy for the highly active antiretroviral therapy associated to the acquired immunodeficiency syndrome should be designed regarding a comprehensive analysis of the drug chemotherapy behavior in the host tissues, from major viral replication sites to viral sanctuary compartments. Such approach is critical in order to efficiently explore synergistic, competitive and prohibitive relationships among drugs and, hence, therapy costs and side-effect minimization. In this paper, a novel mathematical model for HIV-1 drug chemotherapy dynamics in distinct host anatomic compartments is proposed and theoretically evaluated on fifteen conventional anti-retroviral drugs. Rather than interdependence between drug type and its concentration profile in a host tissue, simulated results suggest that such profile is importantly correlated with the host tissue under consideration. Furthermore, the drug accumulative dynamics are drastically affected by low patient compliance with pharmacotherapy, even when a single dose lacks. (C) 2012 Elsevier Inc. All rights reserved.

Extracellular matrix in airway smooth muscle is associated with dynamics of airway function in asthma

Background: Altered deposition of extracellular matrix (ECM) in the airway smooth muscle (ASM) layer as observed in asthma may influence ASM mechanical properties. We hypothesized that ECM in ASM is associated with airway function in asthma. First, we investigated the difference in ECM expression in ASM between asthma and controls. Second, we examined whether ECM expression is associated with bronchoconstriction and bronchodilation in vivo. Methods: Our cross-sectional study comprised 19 atopic mild asthma patients, 15 atopic and 12 nonatopic healthy subjects. Spirometry, methacholine responsiveness, deep-breath-induced bronchodilation (Delta R-rs) and bronchoscopy with endobronchial biopsies were performed. Positive staining of elastin, collagen I, III and IV, decorin, versican, fibronectin, laminin and tenascin in ASM was quantified as fractional area and mean density. Data were analysed using Pearson's or Spearman's correlation coefficient. Results: Extracellular matrix expression in ASM was not different between asthma and controls. In asthmatics, fractional area and mean density of collagen I and III were correlated with methacholine dose-response slope and DRrs, respectively (r = 0.71, P < 0.01; r = 0.60, P = 0.02). Furthermore, ASM collagen III and laminin in asthma were correlated with FEV1 reversibility (r = -0.65, P = 0.01; r = -0.54, P = 0.04). Conclusion: In asthma, ECM in ASM is related to the dynamics of airway function in the absence of differences in ECM expression between asthma and controls. This indicates that the ASM layer in its full composition is a major structural component in determining variable airways obstruction in asthma.

Effect of the thermostat in the molecular dynamics simulation on the folding of the model protein chignolin

Molecular dynamics simulations of the model protein chignolin with explicit solvent were carried out, in order to analyze the influence of the Berendsen thermostat on the evolution and folding of the peptide. The dependence of the peptide behavior on temperature was tested with the commonly employed thermostat scheme consisting of one thermostat for the protein and another for the solvent. The thermostat coupling time of the protein was increased to infinity, when the protein is not in direct contact with the thermal bath, a situation known as minimally invasive thermostat. In agreement with other works, it was observed that only in the last situation the instantaneous temperature of the model protein obeys a canonical distribution. As for the folding studies, it was shown that, in the applications of the commonly utilized thermostat schemes, the systems are trapped in local minima regions from which it has difficulty escaping. With the minimally invasive thermostat the time that the protein needs to fold was reduced by two to three times. These results show that the obstacles to the evolution of the extended peptide to the folded structure can be overcome when the temperature of the peptide is not directly controlled.

Exact dynamics of single-qubit-gate fidelities under the measurement-based quantum computation scheme

Measurement-based quantum computation is an efficient model to perform universal computation. Nevertheless, theoretical questions have been raised, mainly with respect to realistic noise conditions. In order to shed some light on this issue, we evaluate the exact dynamics of some single-qubit-gate fidelities using the measurement-based quantum computation scheme when the qubits which are used as a resource interact with a common dephasing environment. We report a necessary condition for the fidelity dynamics of a general pure N-qubit state, interacting with this type of error channel, to present an oscillatory behavior, and we show that for the initial canonical cluster state, the fidelity oscillates as a function of time. This state fidelity oscillatory behavior brings significant variations to the values of the computational results of a generic gate acting on that state depending on the instants we choose to apply our set of projective measurements. As we shall see, considering some specific gates that are frequently found in the literature, the fast application of the set of projective measurements does not necessarily imply high gate fidelity, and likewise the slow application thereof does not necessarily imply low gate fidelity. Our condition for the occurrence of the fidelity oscillatory behavior shows that the oscillation presented by the cluster state is due exclusively to its initial geometry. Other states that can be used as resources for measurement-based quantum computation can present the same initial geometrical condition. Therefore, it is very important for the present scheme to know when the fidelity of a particular resource state will oscillate in time and, if this is the case, what are the best times to perform the measurements.

Records of sedimentary dynamics in the continental shelf and upper slope between Aveiro-Espinho (N Portugal)

The sedimentary unconsolidated cover of the Aveiro-Espinho continental shelf and upper slope (NW Portugal) records a complex interplay of processes including wave energy and currents, fluvial input, sediment transport alongshore and cross-shelf, geological and oceanographic processes and sediment sources and sinks. In order to study this record, a set of surface sediment samples was studied. Sediment grain size and composition, as well as the mineralogical composition (by XRD) of the fine (<63 mu m) and clay (<2 mu m) fractions and benthic microfaunal (foraminifera) data were analysed. Cluster analysis applied to the sedimentological data (grain size, sediment composition and mineralogy) allowed the establishment of three main zones corresponding to the: inner-, mid- and outer-shelf/upper slope. On the inner-shelf, the sedimentary coverture is composed of siliciclastic fine to very fine sand, essentially comprising modern (immature) terrigenous particles. The sediment grain size, as well as mineralogical and microfaunal composition, denote the high energetic conditions of this sector in which the alongshore transport of sand is predominantly southward and occurs mostly during the spring-summer oceanographic regime, when the main river providing sediments to this area, the River Douro, undergoes periods of drought. This effect may emphasize the erosive character of this coastal sector at present, since the Ria de Aveiro provides the shelf with few sediments. On the mid-shelf, an alongshore siliciclastic band of coarse sand and gravel can be found between the 40 m and 60 m isobaths. This gravelly deposit includes relic sediments deposited during lower sea-level stands. This structure stays on the surface due to the high bottom energy, which promotes the remobilization of the fine-grained sediments, and/or events of sediments bypassing. Benthic foraminifera density and "Benthic Foraminifera High Productivity" (BFHP) proxy values are in general low, which is consistent with the overall small supply of organic matter to the oceanic bottom in the inner- and mid-shelf. However, the Ria de Aveiro outflow, which delivers organic matter to the shelf, leaves its imprint mainly on the mid-shelf, identifiable by the increase in foraminifera density and BFHP values in front of the lagoon mouth. The higher values of BFHP along the 100 m isobath trace the present position of an oceanic thermal front whose situation may have changed in the last 3/5 ka BP. This zone marks a clear difference in the density, diversity and composition of benthic foraminifera assemblages. Here, in addition, sediment composition changes significantly, giving rise to carbonate-rich fine to medium sand in the deeper sector. The low bottom energy and the small sedimentation rate of the outer-shelf contributed to the preservation of a discontinuous carbonate-rich gravel band, between the 100 m and 140 m isobaths, also related to paleo-littorals, following the transgression that has occurred since the Last Glacial Maximum. The winter oceanographic regime favours the transport of fine grained sediments to the outer-shelf and upper slope. The inner- and mid-shelf, however, have low amounts of this kind of sediment and the Cretacic carbonated complexes Pontal da Galega and Pontal da Cartola, rocky outcrops located at the mid- and outer-shelf, act as morphological barriers to the cross-shelf transport of sediments. Thus a reduced sedimentation rate occurs in these deeper sectors, as indicated by the lower abundance of detrital minerals, which is compensated for the high sedimentary content of biogenic carbonates. The relatively high BFHP and Shannon Index values indicate water column stratification, high supply of organic matter and environmental stability, which provide favourable conditions for a diversified benthic fauna to flourish. These conditions also encourage authigenic chemical changes, favourable to glauconite formation, as well as illite and kaolinite degradation. Benthic foraminifera and clay mineral assemblages also reveal the effect of the internal waves pushing upward, and downslope losses of the sediments on the outer-shelf and upper slope.

Assessing the effects of multiple infections and long latency in the dynamics of tuberculosis

In order to achieve a better understanding of multiple infections and long latency in the dynamics of Mycobacterium tuberculosis infection, we analyze a simple model. Since backward bifurcation is well documented in the literature with respect to the model we are considering, our aim is to illustrate this behavior in terms of the range of variations of the model's parameters. We show that backward bifurcation disappears (and forward bifurcation occurs) if: (a) the latent period is shortened below a critical value; and (b) the rates of super-infection and re-infection are decreased. This result shows that among immunosuppressed individuals, super-infection and/or changes in the latent period could act to facilitate the onset of tuberculosis. When we decrease the incubation period below the critical value, we obtain the curve of the incidence of tuberculosis following forward bifurcation; however, this curve envelops that obtained from the backward bifurcation diagram.

Applying many-body expansion aiming to make ab in initio molecular dynamics more efficient

In this present work we present a methodology that aims to apply the many-body expansion to decrease the computational cost of ab initio molecular dynamics, keeping acceptable accuracy on the results. We implemented this methodology in a program which we called ManBo. In the many-body expansion approach, we partitioned the total energy E of the system in contributions of one body, two bodies, three bodies, etc., until the contribution of the Nth body [1-3]: E = E1 + E2 + E3 + …EN. The E1 term is the sum of the internal energy of the molecules; the term E2 is the energy due to interaction between all pairs of molecules; E3 is the energy due to interaction between all trios of molecules; and so on. In Manbo we chose to truncate the expansion in the contribution of two or three bodies, both for the calculation of the energy and for the calculation of the atomic forces. In order to partially include the many-body interactions neglected when we truncate the expansion, we can include an electrostatic embedding in the electronic structure calculations, instead of considering the monomers, pairs and trios as isolated molecules in space. In simulations we made we chose to simulate water molecules, and use the Gaussian 09 as external program to calculate the atomic forces and energy of the system, as well as reference program for analyzing the accuracy of the results obtained with the ManBo. The results show that the use of the many-body expansion seems to be an interesting approach for reducing the still prohibitive computational cost of ab initio molecular dynamics. The errors introduced on atomic forces in applying such methodology are very small. The inclusion of an embedding electrostatic seems to be a good solution for improving the results with only a small increase in simulation time. As we increase the level of calculation, the simulation time of ManBo tends to largely decrease in relation to a conventional BOMD simulation of Gaussian, due to better scalability of the methodology presented. References [1] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 3, 46 (2007). [2] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 4, 1 (2008). [3] R. Rivelino, P. Chaudhuri and S. Canuto; J. Chem. Phys., 118, 10593 (2003).

Model for the Dynamics of the Cytoskeleton.

Particle tracking of microbeads attached to the cytoskeleton (CSK) reveals an intermittent dynamic. The mean squared displacement (MSD) is subdiffusive for small Δt and superdiffusive for large Δt, which are associated with periods of traps and periods of jumps respectively. The analysis of the displacements has shown a non-Gaussian behavior, what is indicative of an active motion, classifying the cells as a far from equilibrium material. Using Langevin dynamics, we reconstruct the dynamic of the CSK. The model is based on the bundles of actin filaments that link themself with the bead RGD coating, trapping it in an harmonic potential. We consider a one- dimensional motion of a particle, neglecting inertial effects (over-damped Langevin dynamics). The resultant force is decomposed in friction force, elastic force and random force, which is used as white noise representing the effect due to molecular agitation. These description until now shows a static situation where the bead performed a random walk in an elastic potential. In order to modeling the active remodeling of the CSK, we vary the equilibrium position of the potential. Inserting a motion in the well center, we change the equilibrium position linearly with time with constant velocity. The result found exhibits a MSD versus time ’tau’ with three regimes. The first regime is when ‘tau’ < ‘tau IND 0’, where ‘tau IND 0’ is the relaxation time, representing the thermal motion. At this regime the particle can diffuse freely. The second regime is a plateau, ‘tau IND 0’ < ‘tau’ < ‘tau IND 1’, representing the particle caged in the potential. Here, ‘tau IND 1’ is a characteristic time that limit the confinement period. And the third regime, ‘tau’ > ‘tau IND 1’, is when the particles are in the superdiffusive behavior. This is where most of the experiments are performed, under 20 frames per second (FPS), thus there is no experimental evidence that support the first regime. We are currently performing experiments with high frequency, up to 100 FPS, attempting to visualize this diffusive behavior. Beside the first regime, our simple model can reproduce MSD curves similar to what has been found experimentally, which can be helpful to understanding CSK structure and properties.