Using likelihood-free inference to compare evolutionary dynamics of the protein networks of H. pylori and P. falciparum

Gene duplication with subsequent interaction divergence is one of the primary driving forces in the evolution of genetic systems. Yet little is known about the precise mechanisms and the role of duplication divergence in the evolution of protein networks from the prokaryote and eukaryote domains. We developed a novel, model-based approach for Bayesian inference on biological network data that centres on approximate Bayesian computation, or likelihood-free inference. Instead of computing the intractable likelihood of the protein network topology, our method summarizes key features of the network and, based on these, uses a MCMC algorithm to approximate the posterior distribution of the model parameters. This allowed us to reliably fit a flexible mixture model that captures hallmarks of evolution by gene duplication and subfunctionalization to protein interaction network data of Helicobacter pylori and Plasmodium falciparum. The 80% credible intervals for the duplication–divergence component are [0.64, 0.98] for H. pylori and [0.87, 0.99] for P. falciparum. The remaining parameter estimates are not inconsistent with sequence data. An extensive sensitivity analysis showed that incompleteness of PIN data does not largely affect the analysis of models of protein network evolution, and that the degree sequence alone barely captures the evolutionary footprints of protein networks relative to other statistics. Our likelihood-free inference approach enables a fully Bayesian analysis of a complex and highly stochastic system that is otherwise intractable at present. Modelling the evolutionary history of PIN data, it transpires that only the simultaneous analysis of several global aspects of protein networks enables credible and consistent inference to be made from available datasets. Our results indicate that gene duplication has played a larger part in the network evolution of the eukaryote than in the prokaryote, and suggests that single gene duplications with immediate divergence alone may explain more than 60% of biological network data in both domains.     

Apicoplast lipoic acid protein ligase B is not essential for Plasmodium falciparum

Lipoic acid (LA) is an essential cofactor of alpha-keto acid dehydrogenase complexes (KADHs) and the glycine cleavage system. In Plasmodium, LA is attached to the KADHs by organelle-specific lipoylation pathways. Biosynthesis of LA exclusively occurs in the apicoplast, comprising octanoyl-[acyl carrier protein]: protein N-octanoyltransferase (LipB) and LA synthase. Salvage of LA is mitochondrial and scavenged LA is ligated to the KADHs by LA protein ligase 1 (LplA1). Both pathways are entirely independent, suggesting that both are likely to be essential for parasite survival. However, disruption of the LipB gene did not negatively affect parasite growth despite a drastic loss of LA (>90%). Surprisingly, the sole, apicoplast-located pyruvate dehydrogenase still showed lipoylation, suggesting that an alternative lipoylation pathway exists in this organelle. We provide evidence that this residual lipoylation is attributable to the dual targeted, functional lipoate protein ligase 2 (LplA2). Localisation studies show that LplA2 is present in both mitochondrion and apicoplast suggesting redundancy between the lipoic acid protein ligases in the erythrocytic stages of P. falciparum.     

GROmaρs: A GROMACS-Based Toolset to Analyze Density Maps Derived from Molecular Dynamics Simulations

We introduce a computational toolset, named GROmaρs, to obtain and compare time-averaged density maps from molecular dynamics simulations. GROmaρs efficiently computes density maps by fast multi-Gaussian spreading of atomic densities onto a three-dimensional grid. It complements existing map-based tools by enabling spatial inspection of atomic average localization during the simulations. Most importantly, it allows the comparison between computed and reference maps (e.g., experimental) through calculation of difference maps and local and time-resolved global correlation. These comparison operations proved useful to quantitatively contrast perturbed and control simulation data sets and to examine how much biomolecular systems resemble both synthetic and experimental density maps. This was especially advantageous for multimolecule systems in which standard comparisons like RMSDs are difficult to compute. In addition, GROmaρs incorporates absolute and relative spatial free-energy estimates to provide an energetic picture of atomistic localization. This is an open-source GROMACS-based toolset, thus allowing for static or dynamic selection of atoms or even coarse-grained beads for the density calculation. Furthermore, masking of regions was implemented to speed up calculations and to facilitate the comparison with experimental maps. Beyond map comparison, GROmaρs provides a straightforward method to detect solvent cavities and average charge distribution in biomolecular systems. We employed all these functionalities to inspect the localization of lipid and water molecules in aquaporin systems, the binding of cholesterol to the G protein coupled chemokine receptor type 4, and the identification of permeation pathways through the dermicidin antimicrobial channel. Based on these examples, we anticipate a high applicability of GROmaρs for the analysis of molecular dynamics simulations and their comparison with experimentally determined densities.     

High Efficiency Gene Transfer by Multiple Transfection Protocol

A high efficiency transfection protocol employing a common polycationic lipid is described. Using LipofectAMINE, a widely used transfection reagent, we transfected 293T cells with a plasmid harboring the -galactosidase (-gal) gene. The transfection efficiency was determined by direct staining for X-gal. The conventional transfection protocol achieved an efficiency of <40% while our protocol, which employs the repetition of transfection a few times, achieved an efficiency of approximately 80%. Thus, a dramatic increase in transfection efficiency can be obtained by simply repeating transfection with the use of a common polycationic lipid. This method will be useful in many molecular biological experiments.     

Ligand-mediated conformational changes of the VDR are required for gene transactivation

The central element of the molecular switch of nuclear 1alpha,25-dihydroxyvitamin D(3) (1alpha,25(OH)(2)D(3)) signaling is the ligand-binding domain (LBD) of the Vitamin D receptor (VDR), which can be stabilized by 1alpha,25(OH)(2)D(3) or its analogues in to agonistic, antagonistic or inverse agonistic conformations. The positioning of helix 12 of the LBD is of most critical importance for these conformations, because it determines the distance between the charge clamp amino acids K246 and E420 that are essential for VDR-coactivator (CoA) interaction. Most VDR ligands have been identified as agonists and only a few (e.g., ZK168281 and TEI-9647) as pure or partial antagonists. Antagonists induce corepressor (CoR) dissociation from the VDR but prevent completely or partially CoA interaction and thus transactivation. Gemini is a 1alpha,25(OH)(2)D(3) analogue with two identical side chains that despite its significantly increased volume binds to the VDR and acts under most conditions as an agonist. Interestingly, supramolar CoR concentrations shift Gemini from an agonist to an inverse agonist, which actively recruits CoR to the VDR and thus mediates repression of 1alpha,25(OH)(2)D(3) target genes. Gemini is the first described (conditional) inverse agonist to an endocrine nuclear receptor (NR) and may function as a sensor for cell-specific CoA/CoR ratios.     

Marsupial relationships and a timeline for marsupial radiation in South Gondwana

Recent marsupials include about 280 species divided into 18 families and seven orders. Approximately 200 species live in Australia/New Guinea. The remaining species inhabit South America with some of these secondarily ranging into North America. In this study, we examine marsupial relationships and estimate their divergences times using complete mitochondrial (mt) genomes. The sampling, which includes nine new mtDNAs and a total number of 19 marsupial genomes, encompasses all extant orders and 14 families. The analysis identified a basal split between Didelphimorphia and remaining orders about 69 million years before present (MYBP), while other ordinal divergences were placed in Tertiary times. The monotypic South American order Microbiotheria (Dromiciops gliroides, Monito del Monte) was solidly nested among its Australian counterparts. The results suggest that marsupials colonized Australia twice from Antarctica/South America and that the divergence between Microbiotheria and its Australian relatives coincided with the geological separation of Antarctica and Australia. Within Australia itself, several of the deepest divergences were estimated to have taken place close to the Eocene/Oligocene transition.     

Lignification and structural biomass production in tobacco with suppressed caffeic/5-hydroxy ferulic acid-O-methyl transferase activity under ambient and elevated CO concentrations

To investigate the relationship between growth, biomass partitioning and lignification we used tobacco (Nicotiana tabacum) in which O-methyl transferase (OMT) activity, an enzyme involved in the pathway of sinapyl alcohol formation for lignin synthesis, was suppressed by antisense transformation. To modulate growth, controls and transformed tobacco plants were grown under ambient (approximately 380 p.p.m) or elevated CO(2) (700 p.p.m), respectively. Lignin concentrations and composition were determined with spectrophotometric methods (thioglycolate and acetyl bromide) and Fourier transform infrared (FTIR) spectroscopy, respectively. A comparison of the thioglycolate and acetylbromide method suggested that the thioglycolate method was sensitive to changes in the syringyl/guaiacyl (S/G)-ratio in lignin and therefore not suitable for quantification in tissues with altered lignin composition. Growth under elevated CO(2) increased leaf and stem biomass of both genotypes by 40 and 20%, respectively, compared with ambient CO(2) and had no effect on root biomass. OMT suppression did not affect lignin concentrations in the leaves but caused a shift in biomass partitioning from the structural to the non-structural fraction. Elevated CO(2) caused a shift towards production of structural compounds resulting in decreased foliar lignin concentrations and indicating that the lignin/structural mass ratio is flexible in leaves. By contrast, the lignin concentrations of stems were unaffected by elevated CO(2) or OMT suppression and increased about three-fold from the apex to the base. Antisense-OMT plants produced more stem biomass than controls but showed no changes of the relative partitioning of biomass to the different pools. This indicates that the metabolic control of carbon fluxes to the production of structural versus non-structural fractions is tighter in stems than in leaves. FTIR spectroscopy indicated a relative increase in guaiacyl- as compared with syringyl-units in antisense-OMT tobacco, which was more pronounced under elevated as compared with ambient CO(2). Since there was no evidence for decreased lignin concentrations in the antisense-OMT plants but increased biomass formation we suggest that less methylated lignins are 'cheaper' in biosynthetic carbon and energy demand and, thus, may enable plants to allocate increased resources to growth.