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.     

Distinct neurocomputational mechanisms support informational and socially normative conformity

A change of mind in response to social influence could be driven by informational conformity to increase accuracy, or by normative conformity to comply with social norms such as reciprocity. Disentangling the behavioural, cognitive, and neurobiological underpinnings of informational and normative conformity have proven elusive. Here, participants underwent fMRI while performing a perceptual task that involved both advice-taking and advice-giving to human and computer partners. The concurrent inclusion of 2 different social roles and 2 different social partners revealed distinct behavioural and neural markers for informational and normative conformity. Dorsal anterior cingulate cortex (dACC) BOLD response tracked informational conformity towards both human and computer but tracked normative conformity only when interacting with humans. A network of brain areas (dorsomedial prefrontal cortex (dmPFC) and temporoparietal junction (TPJ)) that tracked normative conformity increased their functional coupling with the dACC when interacting with humans. These findings enable differentiating the neural mechanisms by which different types of conformity shape social changes of mind.

When we change our mind in response to other people’s opinion, we may be motivated to be correct or to have a good relationship with others. This fMRI study disentangles the neurobiological underpinnings of these two different motives in the human brain, and shows that the second motive is absent when we interact with computers.     

Method for spiking soil samples with organic compounds

We examined the harmful side effects on indigenous soil microorganisms of two organic solvents, acetone and dichloromethane, that are normally used for spiking of soil with polycyclic aromatic hydrocarbons for experimental purposes. The solvents were applied in two contamination protocols to either the whole soil sample or 25% of the soil volume, which was subsequently mixed with 75% untreated soil. For dichloromethane, we included a third protocol, which involved application to 80% of the soil volume with or without phenanthrene and introduction of Pseudomonas fluorescens VKI171 SJ132 genetically tagged with luxAB::Tn5. For both solvents, application to the whole sample resulted in severe side effects on both indigenous protozoa and bacteria. Application of dichloromethane to the whole soil volume immediately reduced the number of protozoa to below the detection limit. In one of the soils, the protozoan population was able to recover to the initial level within 2 weeks, in terms of numbers of protozoa; protozoan diversity, however, remained low. In soil spiked with dichloromethane with or without phenanthrene, the introduced P. fluorescens VKI171 SJ132 was able to grow to a density 1,000-fold higher than in control soil, probably due mainly to release of predation from indigenous protozoa. In order to minimize solvent effects on indigenous soil microorganisms when spiking native soil samples with compounds having a low water solubility, we propose a common protocol in which the contaminant dissolved in acetone is added to 25% of the soil sample, followed by evaporation of the solvent and mixing with the remaining 75% of the soil sample.     

Generation,Purification, and Characterization of Cell-invasive DISC1 Protein Species

Protein aggregation is seen as a general hallmark of chronic, degenerative brain conditions like, for example, in the neurodegenerative diseases Alzheimer''s disease (Aβ, tau), Parkinson''s Disease (α-synuclein), Huntington''s disease (polyglutamine, huntingtin), and others. Protein aggregation is thought to occur due to disturbed proteostasis, i.e. the imbalance between the arising and degradation of misfolded proteins. Of note, the same proteins are found aggregated in sporadic forms of these diseases that are mutant in rare variants of familial forms.Schizophrenia is a chronic progressive brain condition that in many cases goes along with a permanent and irreversible cognitive deficit. In a candidate gene approach, we investigated whether Disrupted-in-schizophrenia 1 (DISC1), a gene cloned in a Scottish family with linkage to chronic mental disease^{1, 2}, could be found as insoluble aggregates in the brain of sporadic cases of schizophrenia³. Using the SMRI CC, we identified in approximately 20 % of cases with CMD but not normal controls or patients with neurodegenerative diseases sarkosyl-insoluble DISC1 immunoreactivity after biochemical fractionation. Subsequent studies in vitro revealed that the aggregation propensity of DISC1 was influenced by disease-associated polymorphism S704C⁴, and that DISC1 aggresomes generated in vitro were cell-invasive⁵, similar to what had been shown for Aβ⁶, tau^7-9, α-synuclein¹⁰, polyglutamine¹¹, or SOD1 aggregates¹². These findings prompted us to propose that at least a subset of cases with CMD, those with aggregated DISC1 might be protein conformational disorders. Here we describe how we generate DISC1 aggresomes in mammalian cells, purify them on a sucrose gradient and use them for cell-invasiveness studies. Similarly, we describe how we generate an exclusively multimeric C-terminal DISC1 fragment, label and purify it for cell invasiveness studies. Using the recombinant multimers of DISC1 we achieve similar cell invasiveness as for a similarly labeled synthetic α-synuclein fragment. We also show that this fragment is taken up in vivo when stereotactically injected into the brain of recipient animals.     

Generation of activated and antigen-specific T cells with cytotoxic activity after co-culture with dendritic cells

Co-culturing of immunological effector cells with antigen-pulsed DC leads to an increase of cytotoxic activity against antigen-expressing tumour cells. Using this approach, we could detect up to 2.8% antigen-specific CTLs after co-culture with antigen-pulsed DC. However, the required high effector cell numbers remain a major obstacle in immunotherapy. In this study, we show an approach for generating activated and antigen-specific effector cells that enables us to decrease effector to target cell ratios. We used an interferon-gamma secretion assay to enrich activated effector cells after co-culture with antigen-pulsed dendritic cells (DC). Purified immunological effector cells lysed 58.3% of antigen-expressing tumour cells at an effector to target ratio of 1:1. Furthermore, using MHC-IgG complexes, we enriched effector cells expressing antigen-specific T-cell receptor after co-culture with DC. Performing ELISpot, flow cytometry and TCR analysis, we could show a significant increase of activated and specific TCR-expressing effector cells after co-culture with DC.     

Cloning and heterologous expression of an enantioselective amidase from Rhodococcus erythropolis strain MP50

The gene for an enantioselective amidase was cloned from Rhodococcus erythropolis MP50, which utilizes various aromatic nitriles via a nitrile hydratase/amidase system as nitrogen sources. The gene encoded a protein of 525 amino acids which corresponded to a protein with a molecular mass of 55.5 kDa. The deduced complete amino acid sequence showed homology to other enantioselective amidases from different bacterial genera. The nucleotide sequence approximately 2.5 kb upstream and downstream of the amidase gene was determined, but no indications for a structural coupling of the amidase gene with the genes for a nitrile hydratase were found. The amidase gene was carried by an approximately 40-kb circular plasmid in R. erythropolis MP50. The amidase was heterologously expressed in Escherichia coli and shown to hydrolyze 2-phenylpropionamide, alpha-chlorophenylacetamide, and alpha-methoxyphenylacetamide with high enantioselectivity; mandeloamide and 2-methyl-3-phenylpropionamide were also converted, but only with reduced enantioselectivity. The recombinant E. coli strain which synthesized the amidase gene was shown to grow with organic amides as nitrogen sources. A comparison of the amidase activities observed with whole cells or cell extracts of the recombinant E. coli strain suggested that the transport of the amides into the cells becomes the rate-limiting step for amide hydrolysis in recombinant E. coli strains.     

The brachiopod fold: a neglected body plan hypothesis

Attention is drawn to Nielsen's radical body plan concept, here named the 'brachiopod fold hypothesis', under which brachiopods and phoronids are recognized to be transversely folded across the ontogenetic anterior–posterior axis so that, to make useful comparisons with other phyla, these organisms must be conceptually unfolded. Under the hypothesis brachiopod brachial and pedicle shell valves are respectively 'anterior' and 'posterior' rather than 'dorsal' and 'ventral' as traditionally described. The hypothesis makes sense of the symmetry axes of the brachiopod shell, is consistent with various indications from fossil and Recent brachiopods, and gives rise to predicted patterns of axis–determining gene expression that differ from those obtaining under the traditional view of the body plan, whilst the variety of folding movements in different lineages implies that superficially dissimilar morphogenetic folds may be fundamentally homologous. Convergent folding patterns are noted in some other organisms. A previous conjecture that inarticulate linguloid brachiopods were derived from halkieriid–like ancestors is elaborated with proposals that recognize possible functional continuities of coelomic and marginal sclerite functions, and it is noted that an ancestrally facultative fold could have become incorporated by genetic assimilation into the brachiopod developmental program. An experimental approach is outlined to test the possibility that some members of the 'small shelly fauna' may have been members of the halkieriid–like brachiopod stem lineage and it is also suggested that buoyancy modification may have been an important function of mineralization amongst Lower Cambrian floaters and swimmers, since negative buoyancy would facilitate access to the benthic niche.     

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.     

Freeze-quench 57Fe-Mössbauer spectroscopy: trapping reactive intermediates

⁵⁷Fe-Mössbauer spectroscopy is a method that probes transitions between the nuclear ground state (I = 1/2) and the first nuclear excited state (I = 3/2). This technique provides detailed information about the chemical environment and electronic structure of iron. Therefore, it has played an important role in studies of the numerous iron-containing proteins and enzymes. In conjunction with the freeze-quench method, ⁵⁷Fe-Mössbauer spectroscopy allows for monitoring changes of the iron site(s) during a biochemical reaction. This approach is particularly powerful for detection and characterization of reactive intermediates. Comparison of experimentally determined Mössbauer parameters to those predicted by density functional theory for hypothetical model structures can then provide detailed insight into the structures of reactive intermediates. We have recently used this methodology to study the reactions of various mononuclear non-heme-iron enzymes by trapping and characterizing several Fe(IV)-oxo reaction intermediates. In this article, we summarize these findings and demonstrate the potential of the method.