Molecular Dynamics Simulations of Phospholipid Bilayers

Abstract

Molecular dynamics (MD) simulations at 37°C have been performed on three phospholipid bilayer systems composed of the lipids DLPE, DOPE, and DOPC. The model used included 24 explicit lipid molecules and explicit waters of solvation in the polar head group regions, together with constant-pressure periodic boundary conditions in three dimensions. Using this model, a MD simulation samples part of an infinite planar lipid bilayer. The lipid dynamics and packing behavior were characterized. Furthermore, using the results of the simulations, a number of diverse properties including bilayer structural parameters, hydrocarbon chain order parameters, dihedral conformations, electron density profile, hydration per lipid, and water distribution along the bilayer normal were calculated. Many of these properties are available for the three lipid systems chosen, making them well suited for evaluating the model and protocols used in these simulations by direct comparisons with experimental data. The calculated MD behavior, chain disorder, and lipid packing parameter, i.e. the ratio of the effective areas of hydrocarbon tails and head group per lipid (a_t/a_h), correctly predict the aggregation preferences of the three lipids observed experimentally at 37°C, namely: a gel bilayer for DLPE, a hexagonal tube for DOPE, and a liquid crystalline bilayer for DOPC. In addition, the model and conditions used in the MD simulations led to good agreement of the calculated properties of the bilayers with available experimental results, demonstrating the reliability of the simulations. The effects of the cis unsaturation in the hydrocarbon chains of DOPE and DOPC, compared to the fully saturated one in DLPE, as well as the effects of the different polar head groups of PC and PE with the same unsaturated chains on the lipid packing and bilayer structure have been investigated. The results of these studies indicate the ability of MD methods to provide molecular-level insights into the structure and dynamics of lipid assemblies.     

Sense-overlapping lncRNA as a decoy of translational repressor protein for dimorphic gene expression

Long noncoding RNAs (lncRNAs) are vastly transcribed and extensively studied but lncRNAs overlapping with the sense orientation of mRNA have been poorly studied. We analyzed the lncRNA DAPALR overlapping with the 5´ UTR of the Doublesex1 (Dsx1), the male determining gene in Daphnia magna. By affinity purification, we identified an RNA binding protein, Shep as a DAPALR binding protein. Shep also binds to Dsx1 5´ UTR by recognizing the overlapping sequence and suppresses translation of the mRNA. In vitro and in vivo analyses indicated that DAPALR increased Dsx1 translation efficiency by sequestration of Shep. This regulation was impaired when the Shep binding site in DAPALR was deleted. These results suggest that Shep suppresses the unintentional translation of Dsx1 by setting a threshold; and when the sense lncRNA DAPALR is expressed, DAPALR cancels the suppression caused by Shep. This mechanism may be important to show dimorphic gene expressions such as sex determination and it may account for the binary expression in various developmental processes.     

The Himalayas: Barrier and conduit for gene flow

The Himalayan mountain range is strategically located at the crossroads of the major cultural centers in Asia, the Middle East and Europe. Although previous Y‐chromosome studies indicate that the Himalayas served as a natural barrier for gene flow from the south to the Tibetan plateau, this region is believed to have played an important role as a corridor for human migrations between East and West Eurasia along the ancient Silk Road. To evaluate the effects of the Himalayan mountain range in shaping the maternal lineages of populations residing on either side of the cordillera, we analyzed mitochondrial DNA variation in 344 samples from three Nepalese collections (Newar, Kathmandu and Tamang) and a general population of Tibet. Our results revealed a predominantly East Asian‐specific component in Tibet and Tamang, whereas Newar and Kathmandu are both characterized by a combination of East and South Central Asian lineages. Interestingly, Newar and Kathmandu harbor several deep‐rooted Indian lineages, including M2, R5, and U2, whose coalescent times from this study (U2, >40 kya) and previous reports (M2 and R5, >50 kya) suggest that Nepal was inhabited during the initial peopling of South Central Asia. Comparisons with our previous Y‐chromosome data indicate sex‐biased migrations in Tamang and a founder effect and/or genetic drift in Tamang and Newar. Altogether, our results confirm that while the Himalayas acted as a geographic barrier for human movement from the Indian subcontinent to the Tibetan highland, it also served as a conduit for gene flow between Central and East Asia. Am J Phys Anthropol 151:169–182, 2013. © 2013 Wiley Periodicals, Inc.     

Selective processes in development: implications for the costs and benefits of phenotypic plasticity

Adaptive phenotypic plasticity, the ability of a genotype to develop a phenotype appropriate to the local environment, allows organisms to cope with environmental variation and has implications for predicting how organisms will respond to rapid, human-induced environmental change. This review focuses on the importance of developmental selection, broadly defined as a developmental process that involves the sampling of a range of phenotypes and feedback from the environment reinforcing high-performing phenotypes. I hypothesize that understanding the degree to which developmental selection underlies plasticity is key to predicting the costs, benefits, and consequences of plasticity. First, I review examples that illustrate that elements of developmental selection are common across the development of many different traits, from physiology and immunity to circulation and behavior. Second, I argue that developmental selection, relative to a fixed strategy or determinate (switch) mechanisms of plasticity, increases the probability that an individual will develop a phenotype best matched to the local environment. However, the exploration and environmental feedback associated with developmental selection is costly in terms of time, energy, and predation risk, resulting in major changes in life history such as increased duration of development and greater investment in individual offspring. Third, I discuss implications of developmental selection as a mechanism of plasticity, from predicting adaptive responses to novel environments to understanding conditions under which genetic assimilation may fuel diversification. Finally, I outline exciting areas of future research, in particular exploring costs of selective processes in the development of traits outside of behavior and modeling developmental selection and evolution in novel environments.     

Shrub cover in northern Nunavik: can herbivores limit shrub expansion?

Recent climate changes have increased the primary productivity of many Arctic and subarctic regions. Erected shrub has been shown to increase in abundance over the last decades in northern regions in response to warmer climate. At the same time, caribou herds are declining throughout the circumboreal regions. Based on observation of heavy browsing on shrubs at Deception Bay (Nunavik, Canada), we hypothesized that the densification of shrubs observed in nearby locations did not occur at our study site despite of observed warming because of a recent peak of the Rivière-aux-Feuilles caribou herd. To assess shrub cover changes, we compared a 1972 mosaic of aerial photos to a 2010 satellite image over a 5 km² area, divided into 56 grids of 100 30 m × 30 m cells. Most cells (n = 4,502) did not show any changes in the cover of shrubs but those who did were as likely to increase as to decrease. The relative cover of shrubs in cells who changed was not higher in 2010 (6.1 ± 0.2 %) than in 1972 (7.3 ± 0.4 %). More than 70 % of birch and willow had more than 50 % of their shoot browsed, suggesting that caribou may limit shrub expansion at this site. We cannot rule out that abiotic factors also contribute to the inertia in shrub cover. Increases in shrub abundance reported in Nunavik and elsewhere were located closer to the tree line or in discontinuous permafrost, whereas our site is characterized by herbaceous arctic tundra, continuous permafrost and relatively low annual precipitation.     

The role of biogeography in shaping diversity of the intestinal microbiota in house mice

The microbial communities inhabiting the mammalian intestinal tract play an important role in diverse aspects of host biology. However, little is known regarding the forces shaping variation in these communities and their influence on host fitness. To shed light on the contributions of host genetics, transmission and geography to diversity in microbial communities between individuals, we performed a survey of intestinal microbial communities in a panel of 121 house mice derived from eight locations across Western Europe using pyrosequencing of the bacterial 16S rRNA gene. The host factors studied included population structure estimated by microsatellite loci and mitochondrial DNA, genetic distance and geography. To determine whether host tissue (mucosa)‐associated communities display properties distinct from those of the lumen, both the caecal mucosa and contents were examined. We identified Bacteroides, Robinsoniella and Helicobacter as the most abundant genera in both the caecal content and mucosa‐associated communities of wild house mice. Overall, we found geography to be the most significant factor explaining patterns of diversity in the intestinal microbiota, with a comparatively weaker influence of host population structure and genetic distance. Furthermore, the influence of host genetic distance was limited to the mucosa communities, consistent with this environment being more intimately coupled to the host.     

A bedr way of genomic interval processing

Background

Next-generation sequencing is making it critical to robustly and rapidly handle genomic ranges within standard pipelines. Standard use-cases include annotating sequence ranges with gene or other genomic annotation, merging multiple experiments together and subsequently quantifying and visualizing the overlap. The most widely-used tools for these tasks work at the command-line (e.g. BEDTools) and the small number of available R packages are either slow or have distinct semantics and features from command-line interfaces.

Results

To provide a robust R-based interface to standard command-line tools for genomic coordinate manipulation, we created bedr. This open-source R package can use either BEDTools or BEDOPS as a back-end and performs data-manipulation extremely quickly, creating R data structures that can be readily interfaced with existing computational pipelines. It includes data-visualization capabilities and a number of data-access functions that interface with standard databases like UCSC and COSMIC.

Conclusions

bedr package provides an open source solution to enable genomic interval data manipulation and restructuring in R programming language which is commonly used in bioinformatics, and therefore would be useful to bioinformaticians and genomic researchers.

     

Lipidomic Analysis Links Mycobactin Synthase K to Iron Uptake and Virulence in M. tuberculosis

The prolonged survival of Mycobacterium tuberculosis (M. tb) in the host fundamentally depends on scavenging essential nutrients from host sources. M. tb scavenges non-heme iron using mycobactin and carboxymycobactin siderophores, synthesized by mycobactin synthases (Mbt). Although a general mechanism for mycobactin biosynthesis has been proposed, the biological functions of individual mbt genes remain largely untested. Through targeted gene deletion and global lipidomic profiling of intact bacteria, we identify the essential biochemical functions of two mycobactin synthases, MbtK and MbtN, in siderophore biosynthesis and their effects on bacterial growth in vitro and in vivo. The deletion mutant, ΔmbtN, produces only saturated mycobactin and carboxymycobactin, demonstrating an essential function of MbtN as the mycobactin dehydrogenase, which affects antigenicity but not iron uptake or M. tb growth. In contrast, deletion of mbtK ablated all known forms of mycobactin and its deoxy precursors, defining MbtK as the essential acyl transferase. The mbtK mutant showed markedly reduced iron scavenging and growth in vitro. Further, ΔmbtK was attenuated for growth in mice, demonstrating a non-redundant role of hydroxamate siderophores in virulence, even when other M. tb iron scavenging mechanisms are operative. The unbiased lipidomic approach also revealed unexpected consequences of perturbing mycobactin biosynthesis, including extreme depletion of mycobacterial phospholipids. Thus, lipidomic profiling highlights connections among iron acquisition, phospholipid homeostasis, and virulence, and identifies MbtK as a lynchpin at the crossroads of these phenotypes.     

Quantitative lipopolysaccharide analysis using HPLC/MS/MS and its combination with the limulus amebocyte lysate assay

Quantitation of plasma lipopolysaccharides (LPSs) might be used to document Gram-negative bacterial infection. In the present work, LPS-derived 3-hydroxymyristate was extracted from plasma samples with an organic solvent, separated by reversed phase HPLC, and quantitated by MS/MS. This mass assay was combined with the limulus amebocyte lysate (LAL) bioassay to monitor neutralization of LPS activity in biological samples. The described HPLC/MS/MS method is a reliable, practical, accurate, and sensitive tool to quantitate LPS. The combination of the LAL and HPLC/MS/MS analyses provided new evidence for the intrinsic capacity of plasma lipoproteins and phospholipid transfer protein to neutralize the activity of LPS. In a subset of patients with systemic inflammatory response syndrome, with documented infection but with a negative plasma LAL test, significant amounts of LPS were measured by the HPLC/MS/MS method. Patients with the highest plasma LPS concentration were more severely ill. HPLC/MS/MS is a relevant method to quantitate endotoxin in a sample, to assess the efficacy of LPS neutralization, and to evaluate the proinflammatory potential of LPS in vivo.