Evolutionary dynamics of wAu-Like Wolbachia variants in neotropical Drosophila spp.

Wolbachia bacteria are common intracellular symbionts of arthropods and have been extensively studied in Drosophila. Most research focuses on two Old Word hosts, Drosophila melanogaster and Drosophila simulans, and does not take into account that some of the Wolbachia associations in these species may have evolved only after their fast global expansion and after the exposure to Wolbachia of previously isolated habitats. Here we looked at Wolbachia of Neotropical Drosophila species. Seventy-one lines of 16 Neotropical Drosophild species sampled in different regions and at different time points were analyzed. Wolbachia is absent in lines of Drosophild willistoni collected before the 1970s, but more recent samples are infected with a strain designated wWiL Wolbachia is absent in all other species of the willistoni group. Polymorphic wWil-related strains were detected in some saltans group species, with D. septentriosaltans being coinfected with at least four variants. Based on wsp and ftsZ sequence data, wWil of D. willistoni is identical to wAu, a strain isolated from D. simulans, but can be discriminated when using a polymorphic minisatellite marker. In contrast to wAu, which infects both germ line and somatic tissues of D. simulans, wWil is found exclusively in the primordial germ line cells of D. willistoni embryos. We report on a pool of closely related Wolbachia strains in Neotropical Drosophila species as a potential source for the wAu strain in D. simulans. Possible evolutionary scenarios reconstructing the infection history of wAu-like Wolbachia in Neotropical Drosophild species and the Old World species D. simulans are discussed.

First-principles quantum dynamics in interacting Bose gases II: stochastic gauges

First principles simulations of the quantum dynamics of interacting Bose gases using the stochastic gauge representation are analysed. In a companion paper, we showed how the positive-P representation can be applied to these problems using stochastic differential equations. That method, however, is limited by increased sampling error as time evolves. Here, we show how the sampling error can be greatly reduced and the simulation time significantly extended using stochastic gauges. In particular, local stochastic gauges (a subset) are investigated. Improvements are confirmed in numerical calculations of single-, double- and multi-mode systems in the weak-mode coupling regime. Convergence issues are investigated, including the recognition of two modes by which stochastic equations produced by phase-space methods in general can diverge: movable singularities and a noise-weight relationship. The example calculated here displays wave-like behaviour in spatial correlation functions propagating in a uniform 1D gas after a sudden change in the coupling constant. This could in principle be tested experimentally using Feshbach resonance methods.

First-principles quantum dynamics in interacting Bose gases: I. The positive P representation

The performance of the positive P phase-space representation for exact many- body quantum dynamics is investigated. Gases of interacting bosons are considered, where the full quantum equations to simulate are of a Gross-Pitaevskii form with added Gaussian noise. This method gives tractable simulations of many-body systems because the number of variables scales linearly with the spatial lattice size. An expression for the useful simulation time is obtained, and checked in numerical simulations. The dynamics of first-, second- and third-order spatial correlations are calculated for a uniform interacting 1D Bose gas subjected to a change in scattering length. Propagation of correlations is seen. A comparison is made with other recent methods. The positive P method is particularly well suited to open systems as no conservation laws are hard-wired into the calculation. It also differs from most other recent approaches in that there is no truncation of any kind.

Kinematics, substructure and luminosity-weighted dynamics of six nearby galaxy groups

We have redefined group membership of six southern galaxy groups in the local universe (mean cz < 2000 km s(-1)) based on new redshift measurements from our recently acquired Anglo-Australian Telescope 2dF spectra. For each group, we investigate member galaxy kinematics, substructure, luminosity functions and luminosity-weighted dynamics. Our calculations confirm that the group sizes, virial masses and luminosities cover the range expected for galaxy groups, except that the luminosity of NGC 4038 is boosted by the central starburst merger pair. We find that a combination of kinematical, substructural and dynamical techniques can reliably distinguish loose, unvirialized groups from compact, dynamically relaxed groups. Applying these techniques, we find that Dorado, NGC 4038 and NGC 4697 are unvirialized, whereas NGC 681, NGC 1400 and NGC 5084 are dynamically relaxed.

Mimicking the action of GroEL in molecular dynamics simulations: Application to the refinement of protein structures

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results Sur port the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.

Network topology and the evolution of dynamics in an artificial genetic regulatory network model created by whole genome duplication and divergence

Topological measures of large-scale complex networks are applied to a specific artificial regulatory network model created through a whole genome duplication and divergence mechanism. This class of networks share topological features with natural transcriptional regulatory networks. Specifically, these networks display scale-free and small-world topology and possess subgraph distributions similar to those of natural networks. Thus, the topologies inherent in natural networks may be in part due to their method of creation rather than being exclusively shaped by subsequent evolution under selection. The evolvability of the dynamics of these networks is also examined by evolving networks in simulation to obtain three simple types of output dynamics. The networks obtained from this process show a wide variety of topologies and numbers of genes indicating that it is relatively easy to evolve these classes of dynamics in this model. (c) 2006 Elsevier Ireland Ltd. All rights reserved.

Many-body quantum dynamics of polarization squeezing in optical fibers

We report new experiments that test quantum dynamical predictions of polarization squeezing for ultrashort photonic pulses in a birefringent fiber, including all relevant dissipative effects. This exponentially complex many-body problem is solved by means of a stochastic phase-space method. The squeezing is calculated and compared to experimental data, resulting in excellent quantitative agreement. From the simulations, we identify the physical limits to quantum noise reduction in optical fibers. The research represents a significant experimental test of first-principles time-domain quantum dynamics in a one-dimensional interacting Bose gas coupled to dissipative reservoirs.

Molecular dynamics characterization of n-octyl-beta-D-glucopyranoside micelle structure in aqueous solution

n-Octyl-beta-D-glueopyranoside (OG) is a non-ionic glycolipid, which is used widely in biotechnical and biochemical applications. All-atom molecular dynamics simulations from two different initial coordinates and velocities in explicit solvent have been performed to characterize the structural behaviour of an OG aggregate at equilibrium conditions. Geometric packing properties determined from the simulations and small angle neutron scattering experiment state that OG micelles are more likely to exist in a non-spherical shape, even at the concentration range near to the critical micelle concentration (0.025 M). Despite few large deviations in the principal moment of inertia ratios, the average micelle shape calculated from both simulations is a prolate ellipsoid. The deviations at these time scales are presumably the temporary shape change of a micelle. However, the size of the micelle and the accessible surface areas were constant during the simulations with the micelle surface being rough and partially elongated. Radial distribution functions computed for the hydroxyl oxygen atoms of an OG show sharper peaks at a minimum van der Waals contact distance than the acetal oxygen, ring oxygen, and anomeric carbon atoms. This result indicates that these atoms are pointed outwards at the hydrophilic/hydrophobic interface, form hydrogen bonds with the water molecules, and thus hydrate the micelle surface effectively. (c) 2005 Elsevier Inc. All rights reserved.

Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface

Hydrophobins are small (similar to 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/ hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily P-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in a-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.

Molecular dynamics study of MscL interactions with a curved lipid bilayer

Mechanosensitivity is a ubiquitous sensory mechanism found in living organisms. The simplest known mechanotransducing mechanism is found in bacteria in the form of the mechanosensitive membrane channel of large conductance, MscL. This channel has been studied extensively using a variety of methods at a functional and structural level. The channel is gated by membrane tension in the lipid bilayer alone. It serves as a safety valve protecting bacterial cells against hypoosmotic shock. MscL of Escherichia coli embedded in bilayers composed of asymmetric amounts of single-tailed and double-tailed lipids has been shown to gate spontaneously, even in the absence of membrane tension. To gain insight into the effect of the lipid membrane composition and geometry on MscL structure, a fully solvated, all-atom model of MscL in a stress-free curved bilayer composed of double- and single-tailed lipids was studied using a 9.5-ns molecular dynamics simulation. The bilayer was modeled as a domed structure accommodating the asymmetric composition of the monolayers. During the course of the simulation a spontaneous restructuring of the periplasmic loops occurred, leading to interactions between one of the loops and phospholipid headgroups. Previous experimental studies of the role of the loops agree with the observation that opening starts with a restructuring of the periplasmic loop, suggesting an effect of the curved bilayer. Because of limited resources, only one simulation of the large system was performed. However, the results obtained suggest that through the geometry and composition of the bilayer the protein structure can be affected even on short timescales.

Population dynamics and gene flow of Helicoverpa armigera (Lepidoptera : Noctuidae) on cotton and grain crops in the Murrumbidgee Valley, Australia

The population dynamics of Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae) in the Murrumbidgee Valley, Australia, has been characterized using five highly variable microsatellite loci. In the 2001-2002 growing season, there were very high levels of migration into the Murrumbidgee Valley with no detectable genetic structuring, consistent with previous analyses on a national scale. By contrast, there was significant genetic structuring over the 2002-2003 growing season, with three distinct genetic types detected. The first type corresponded to the first two generations and was derived from local individuals emerging from diapause and their progeny. The second genetic type corresponded to generation 3 and resulted from substantial immigration into the region. There was another genetic shift in generation 4, which accounts for the third genetic type of the season. This genetic shift occurred despite low levels of immigration. During the third generation of the 2002-2003 growing season, different population dynamics was characterized for H. armigera on maize, Zea mays L., and cotton Gossipium hirsutum L. Populations on cotton tended to cycle independently with very little immigration from outside the region or from maize within the region. Maize acted as a major sink for immigrants from cotton and from outside the region. If resistance were to develop on cotton under these circumstances, susceptible individuals from maize or from other regions would not dilute this resistance. In addition, resistance is likely to be transferred to maize and be perpetuated until diapause, from where it may reemerge next season. If low levels of immigration were to occur on transgenic cotton, this may undermine the effectiveness of refugia, especially noncotton refugia.

Position information but not force information is used in adapting to changes in environmental dynamics

This study investigated how movement error is evaluated and used to change feedforward commands following a change in the environmental dynamics. In particular, we addressed the question of whether only position-error information is used or whether information about the force-field direction can also be used for rapid adaptation to changes in the environmental dynamics. Subjects learned to move in a position-dependent force field (PF) with a parabolic profile and the dynamics of a negative spring, which produced lateral force to the left of the target hand path. They adapted very rapidly, dramatically reducing lateral error after a single trial. Several times during training, the strength of the PF was unexpectedly doubled (PF2) for two trials. This again created a large leftward deviation, which was greatly reduced on the second PF2 trial, and an aftereffect when the force field subsequently returned to its original strength. The aftereffect was abolished if the second PF2 trial was replaced by an oppositely directed velocity-dependent force field (VF). During subsequent training in the VF, immediately after having adapted to the PF, subjects applied a force that assisted the force field for similar to 15 trials, indicating that they did not use information about the force-field direction. We concluded that the CNS uses only the position error for updating the internal model of the environmental dynamics and modifying feedforward commands. Although this strategy is not necessarily optimal, it may be the most reliable strategy for iterative improvement in performance.

Predicting population dynamics of weed biological control agents: science or gazing into crystal balls?

Various factors can influence the population dynamics of phytophages post introduction, of which climate is fundamental. Here we present an approach, using a mechanistic modelling package (CLIMEX), that at least enables one to make predictions of likely dynamics based on climate alone. As biological control programs will have minimal funding for basic work (particularly on population dynamics), we show how predictions can be made using a species geographical distribution, relative abundance across its range, seasonal phenology and laboratory rearing data. Many of these data sets are more likely to be available than long-term population data, and some can be incorporated into the exploratory phase of a biocontrol program. Although models are likely to be more robust the more information is available, useful models can be developed using information on species distribution alone. The fitted model estimates a species average response to climate, and can be used to predict likely geographical distribution if introduced, where the agent is likely to be more abundant (i.e. good locations) and more importantly for interpretation of release success, the likely variation in abundance over time due to intra- and inter-year climate variability. The latter will be useful in predicting both the seasonal and long-term impacts of the potential biocontrol agent on the target weed. We believe this tool may not only aid in the agent selection process, but also in the design of release strategies, and for interpretation of post-introduction dynamics and impacts. More importantly we are making testable predictions. If biological control is to become more of a science making and testing such hypothesis will be a key component.

Scaling and non-scaling of muscle activity, kinematics, and dynamics in sit-ups with different degrees of difficulty

The purpose of this study was to investigate how the CNS adjusts motor patterns for variants of a complex axial movement-the situp. Adjustments were induced by changing the support surface contact and mass distribution of the body. Healthy adults performed straight-legged sit-ups, 3 s in duration, with support added to or removed from the lumbar trunk, or with mass added to the head or to the legs. Each of these interventions either increased or decreased the difficulty of the task. The study addressed the extent to which changes in sit-up difficulty are compensated by scaling of muscle activity, kinematics, and dynamics versus the extent to which they are compensated by changing discretely the motor pattern. The analysis of muscle activity, kinematics, and dynamics focused on the first 30-40% of the sit-up-the trunk flexion phase-since this is the most critical part of the movement. Our results demonstrate that, in some respects, sit-up kinematics and dynamics scaled with difficulty, but in other respects, they did not. Muscle activity also scaled, in many respects, but in more difficult sit-ups, abdominal flexor activity decreased instead of increased. Non-scaling changes in these parameters suggest that complex movements, such as the sit-up, may require discrete changes in motor pattern in order to deal with large loads, which challenge the available leverage. (C) 2005 Elsevier Ltd. All rights reserved.