Search

Consistent buckling distortions of a large membrane patch (200 × 200 Å) are observed during molecular dynamics (MD) simulations using the Monte-Carlo (MC) barostat in combination with a hard Lennard-Jones (LJ) cutoff. The buckling behavior is independent of both the simulation engine and the force field but requires the MC barostat-hard LJ cutoff combination. Similar simulations of a smaller patch (90 × 90 Å) do not show buckling, but do show a small, systematic reduction in the surface area accompanied by ~1 Å thickening suggestive of compression. We show that a mismatch in the way potentials and forces are handled in the dynamical equations versus the MC barostat results in a compressive load on the membrane. Moreover, a straightforward application of elasticity theory reveals that a minimal compression of the linear dimensions of the membrane, inversely proportional to the edge length, is required for buckling, explaining this differential behavior. We recommend always using LJ force or potential-switching when the MC barostat is employed to avoid undesirable membrane deformations.

The influence of the membrane on transmembrane proteins is central to a number of biological phenomena, notably the gating of stretch activated ion channels. Conversely, membrane proteins can influence the bilayer, leading to the stabilization of particular membrane shapes, topological changes that occur during vesicle fission and fusion, and shape-dependent protein aggregation. Continuum elastic models of the membrane have been widely used to study protein-membrane interactions. These mathematical approaches produce physically interpretable membrane shapes, energy estimates for the cost of deformation, and a snapshot of the equilibrium configuration. Moreover, elastic models are much less computationally demanding than fully atomistic and coarse-grained simulation methodologies; however, it has been argued that continuum models cannot reproduce the distortions observed in fully atomistic molecular dynamics simulations. We suggest that this failure can be overcome by using chemically and geometrically accurate representations of the protein. Here, we present a fast and reliable hybrid continuum-atomistic model that couples the protein to the membrane. We show that the model is in excellent agreement with fully atomistic simulations of the ion channel gramicidin embedded in a POPC membrane. Our continuum calculations not only reproduce the membrane distortions produced by the channel but also accurately determine the channel's orientation. Finally, we use our method to investigate the role of membrane bending around the charged voltage sensors of the transient receptor potential cation channel TRPV1. We find that membrane deformation significantly stabilizes the energy of insertion of TRPV1 by exposing charged residues on the S4 segment to solution.

The dynamics of suspensions of particles has been an active area of research since Einstein first calculated the leading-order correction to the viscosity of a suspension of spherical particles (Einstein, Proc. R. Soc., vol. A102, 1906, pp. 161-179). Since then, researchers have strived to develop an accurate description of the behaviours of suspensions that goes beyond just leading order in the particle volume fraction. Here, we consider the low-Reynolds-number behaviour of a suspension of spherical particles. Working from the Green's functions for the flow due to a single particle, we derive a continuum-level description of the dynamics of suspensions. Our analysis corrects an error in the derivation of these equations in the work of Jackson (Chem. Engng Sci., vol. 52, 1997, pp. 2457-2469) and leads to stable equations of motion for the particles and fluid. In addition, our resulting equations naturally give the sedimentation speed for suspended particles and correct a separate error in the calculation by Batchelor (J. Fluid Mech., vol. 52, 1972, pp. 245-268). Using the pair-correlation function for hard spheres, we are able to compute the sedimentation speed out to seventh order in the volume fraction, which agrees with experimental data up to 30%-35 %, and also get higher-order corrections to the suspension viscosity, which agree with experiments up to ~15 %. Then, using the pair distribution for spheres in shear flow, we find alterations to both the first and second normal stresses. [ABSTRACT FROM AUTHOR]

Lyme disease is caused by transmission of the spirochete Borrelia burgdorferi from ticks to humans. Although much is known about B. burgdorferi replication, the routes and mechanisms by which it disseminates within the tick remain unclear. To better understand this process, we imaged live, infectious B. burgdorferi expressing a stably integrated, constitutively expressed GFP reporter. Using isolated tick midguts and salivary glands, we observed B. burgdorferi progress through the feeding tick via what we believe to be a novel, biphasic mode of dissemination. In the first phase, replicating spirochetes, positioned at varying depths throughout the midgut at the onset of feeding, formed networks of nonmotile organisms that advanced toward the basolateral surface of the epithelium while adhering to differentiating, hypertrophying, and detaching epithelial cells. In the second phase of dissemination, the nonmotile spirochetes transitioned into motile organisms that penetrated the basement membrane and entered the hemocoel, then migrated to and entered the salivary glands. We designated the first phase of dissemination 'adherence-mediated migration' and provided evidence that it involves the inhibition of spirochete motility by one or more diffusible factors elaborated by the feeding tick midgut. Our studies, which we believe are the first to relate the transmission dynamics of spirochetes to the complex morphological and developmental changes that the midgut and salivary glands undergo during engorgement, challenge the conventional viewpoint that dissemination of Lyme disease-causing spirochetes within ticks is exclusively motility driven., Introduction Almost half of the world's population is at risk for acquiring a potentially life-threatening vector-borne disease, such as malaria, dengue, Chagas disease, and plague (1). Conventional approaches for decreasing [...]

Electron cryotomography was used to analyze the structure of the Lyme disease spirochete, Borrelia burgdorferi. This methodology offers a new means for studying the native architecture of bacteria by eliminating the chemical fixing, dehydration, and staining steps of conventional electron microscopy. Using electron cryotomography, we noted that membrane blebs formed at the ends of the cells. These blebs may be precursors to vesicles that are released from cells grown in vivo and in vitro. We found that the periplasmic space of B. burgdorferi was quite narrow (16.0 nm) compared to those of Escherichia coli and Pseudomonas aeruginosa. However, in the vicinity of the periplasmic flagella, this space was considerably wider (42.3 nm). In contrast to previous results, the periplasmic flagella did not form a bundle but rather formed a tight-fitting ribbon that wraps around the protoplasmic cell cylinder in a right-handed sense. We show how the ribbon configuration of the assembled periplasmic flagella is more advantageous than a bundle for both swimming and forming the flat-wave morphology. Previous results indicate that B. burgdorferi motility is dependent on the rotation of the periplasmic flagella in generating backward-moving waves along the length of the cell. This swimming requires that the rotation of the flagella exerts force on the cell cylinder. Accordingly, a ribbon is more beneficial than a bundle, as this configuration allows each periplasmic flagellum to have direct contact with the cell cylinder in order to exert that force, and it minimizes interference between the rotating filaments.

The filaments of spirochete periplasmic flagella (PFs) have a unique structure and protein composition. In most spirochetes, the PFs consist of a core of at least three related proteins (FlaB1, FlaB2, and FlaB3) and a sheath of FlaA protein. The functions of these filament proteins remain unknown. In this study, we used a multidisciplinary approach to examine the role of these proteins in determining the composition, shape, and stiffness of the PFs and how these proteins impact motility by using the spirochete Brachyspira (formerly Treponema, Serpulina) hyodysenteriae as a genetic model. A series of double mutants lacking combinations of these PF proteins was constructed and analyzed. The results show the following. First, the diameters of PFs are primarily determined by the sheath protein FlaA, and that FlaA can form a sheath in the absence of an intact PF core. Although the sheath is important to the PF structure and motility, it is not essential. Second, the three core proteins play unequal roles in determining PF structure and swimming speed. The functions of the core proteins FlaB1 and FlaB2 overlap such that either one of these proteins is essential for the spirochete to maintain the intact PF structure and for cell motility. Finally, linear elasticity theory indicates that flagellar stiffness directly affects the spirochete's swimming speed.

The swimming movement of a helical bacterium called Spiroplasma is described. Spiroplasma propels itself by generating two temporally distinct kinks that travel the length of the bacterium, thus, indicating that the existence of a contractile apparatus that drives cell movement.