Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born

We present an implementation of generalized Born implicit solvent all-atom classical molecular dynamics (MD) within the AMBER program package that runs entirely on CUDA enabled NVIDIA graphics processing units (GPUs). We discuss the algorithms that are used to exploit the processing power of the GPUs and show the performance that can be achieved in comparison to simulations on conventional CPU clusters. The implementation supports three different precision models in which the contributions to the forces are calculated in single precision floating point arithmetic but accumulated in double precision (SPDP), or everything is computed in single precision (SPSP) or double precision (DPDP). In addition to performance, we have focused on understanding the implications of the different precision models on the outcome of implicit solvent MD simulations. We show results for a range of tests including the accuracy of single point force evaluations and energy conservation as well as structural properties pertainining to protein dynamics. The numerical noise due to rounding errors within the SPSP precision model is sufficiently large to lead to an accumulation of errors which can result in unphysical trajectories for long time scale simulations. We recommend the use of the mixed-precision SPDP model since the numerical results obtained are comparable with those of the full double precision DPDP model and the reference double precision CPU implementation but at significantly reduced computational cost. Our implementation provides performance for GB simulations on a single desktop that is on par with, and in some cases exceeds, that of traditional supercomputers

Exposure Opportunity: The Advantages of Including Men in Analyses of Female-Related Risk Factors

Abstract Intuitively, researchers do not include subjects who do not have the opportunity to be exposed, such as men in studies on oral contraceptives (OCs). We aimed to explore in which situations it is nevertheless beneficial to do so. We considered the effect of including men in case-control analyses of 8 different hypothetical data sets on the effect of OC use and venous thrombosis. In all scenarios, OC use was the exposure of interest, sex the factor that determined exposure opportunity, and air travel another risk factor. In some of these scenarios, sex and air travel were included as confounders or effect modifiers. Logistic regression was used to estimate odds ratios. Standard errors of the estimated log odds ratios, including and excluding men, were compared. We also studied the effect of including men using data from 1999–2004 from a case-control study on risk factors for venous thrombosis, conducted in the Netherlands. In all hypothetical examples, and in the real-data study, addition of men to the analysis yielded the same odds ratios when correctly adjusting for confounding. Moreover, use of additional subjects often led to more precise estimates. We suggest that subjects who do not have the opportunity to be exposed should not routinely be excluded from epidemiologic studies

Salivary and gut microbiomes play a significant role in in vitro oral bioaccessibility, biotransformation, and intestinal absorption of arsenic from food

The release of a toxicant from a food matrix during the gastrointestinal digestion is a crucial determinant of the toxicant's oral bioavailability. We present a modified setup of the human simulator of the gut microbial ecosystem (SHIME), with four sequential gastrointestinal reactors (oral, stomach, small intestine, and colon), including the salivary and colonic microbiomes. Naturally arsenic-containing rice, mussels, and nori seaweed were digested in the presence of microorganisms and in vitro oral bioaccessibility, bioavailability, and metabolism of arsenic species were evaluated following analysis by using HPLC/mass spectrometry. When food matrices were digested with salivary bacteria, the soluble arsenic in the gastric digestion stage increased for mussel and nori samples, but no coincidence impact was found in the small intestinal and colonic digestion stages. However, the simulated small intestinal absorption of arsenic was increased in all food matrices (1.2-2.7 fold higher) following digestion with salivary microorganisms. No significant transformation of the arsenic species occurred except for the arsenosugars present in mussels and nori. In those samples, conversions between the oxo arsenosugars were observed in the small intestinal digestion stage whereupon the thioxo analogs became major metabolites. These results expand our knowledge on the likely metabolism and oral bioavailabiltiy of arsenic during human digestion, and provide valuable information for future risk assessments of dietary arsenic

Transient peak-strain matching partially recovers the age-impaired mechanoadaptive cortical bone response

Mechanoadaptation maintains bone mass and architecture; its failure underlies age-related decline in bone strength. It is unclear whether this is due to failure of osteocytes to sense strain, osteoblasts to form bone or insufficient mechanical stimulus. Mechanoadaptation can be restored to aged bone by surgical neurectomy, suggesting that changes in loading history can rescue mechanoadaptation. We use non-biased, whole-bone tibial analyses, along with characterisation of surface strains and ensuing mechanoadaptive responses in mice at a range of ages, to explore whether sufficient load magnitude can activate mechanoadaptation in aged bone. We find that younger mice adapt when imposed strains are lower than in mature and aged bone. Intriguingly, imposition of short-term, high magnitude loading effectively primes cortical but not trabecular bone of aged mice to respond. This response was regionally-matched to highest strains measured by digital image correlation and to osteocytic mechanoactivation. These data indicate that aged bone’s loading response can be partially recovered, non-invasively by transient, focal high strain regions. Our results indicate that old murine bone does respond to load when the loading is of sufficient magnitude, and bones’ age-related adaptation failure may be due to insufficient mechanical stimulus to trigger mechanoadaptation

Calculations of the dynamical critical exponent using the asymptotic series summation method

We consider how the Pad'e-Borel, Pad'e-Borel-Leroy, and conformal mapping summation methods for asymptotic series can be used to calculate the dynamical critical exponent for homogeneous and disordered Ising-like systems.Comment: 21 RevTeX pages, 2 figure

Solution structure of a repeated unit of the ABA-1 nematode polyprotein allergen of ascaris reveals a novel fold and two discrete lipid-binding sites

Parasitic nematode worms cause serious health problems in humans and other animals. They can induce allergic-type immune responses, which can be harmful but may at the same time protect against the infections. Allergens are proteins that trigger allergic reactions and these parasites produce a type that is confined to nematodes, the nematode polyprotein allergens (NPAs). These are synthesized as large precursor proteins comprising repeating units of similar amino acid sequence that are subsequently cleaved into multiple copies of the allergen protein. NPAs bind small lipids such as fatty acids and retinol (Vitamin A) and probably transport these sensitive and insoluble compounds between the tissues of the worms. Nematodes cannot synthesize these lipids, so NPAs may also be crucial for extracting nutrients from their hosts. They may also be involved in altering immune responses by controlling the lipids by which the immune and inflammatory cells communicate. We describe the molecular structure of one unit of an NPA, the well-known ABA-1 allergen of Ascaris and find its structure to be of a type not previously found for lipid-binding proteins, and we describe the unusual sites where lipids bind within this structur

The transcriptional repressor protein NsrR senses nitric oxide directly via a [2Fe-2S] cluster

The regulatory protein NsrR, a member of the Rrf2 family of transcription repressors, is specifically dedicated to sensing nitric oxide (NO) in a variety of pathogenic and non-pathogenic bacteria. It has been proposed that NO directly modulates NsrR activity by interacting with a predicted [Fe-S] cluster in the NsrR protein, but no experimental evidence has been published to support this hypothesis. Here we report the purification of NsrR from the obligate aerobe Streptomyces coelicolor. We demonstrate using UV-visible, near UV CD and EPR spectroscopy that the protein contains an NO-sensitive [2Fe-2S] cluster when purified from E. coli. Upon exposure of NsrR to NO, the cluster is nitrosylated, which results in the loss of DNA binding activity as detected by bandshift assays. Removal of the [2Fe-2S] cluster to generate apo-NsrR also resulted in loss of DNA binding activity. This is the first demonstration that NsrR contains an NO-sensitive [2Fe-2S] cluster that is required for DNA binding activity

Visualizing sound emission of elephant vocalizations: evidence for two rumble production types

Recent comparative data reveal that formant frequencies are cues to body size in animals, due to a close relationship between formant frequency spacing, vocal tract length and overall body size. Accordingly, intriguing morphological adaptations to elongate the vocal tract in order to lower formants occur in several species, with the size exaggeration hypothesis being proposed to justify most of these observations. While the elephant trunk is strongly implicated to account for the low formants of elephant rumbles, it is unknown whether elephants emit these vocalizations exclusively through the trunk, or whether the mouth is also involved in rumble production. In this study we used a sound visualization method (an acoustic camera) to record rumbles of five captive African elephants during spatial separation and subsequent bonding situations. Our results showed that the female elephants in our analysis produced two distinct types of rumble vocalizations based on vocal path differences: a nasally- and an orally-emitted rumble. Interestingly, nasal rumbles predominated during contact calling, whereas oral rumbles were mainly produced in bonding situations. In addition, nasal and oral rumbles varied considerably in their acoustic structure. In particular, the values of the first two formants reflected the estimated lengths of the vocal paths, corresponding to a vocal tract length of around 2 meters for nasal, and around 0.7 meters for oral rumbles. These results suggest that African elephants may be switching vocal paths to actively vary vocal tract length (with considerable variation in formants) according to context, and call for further research investigating the function of formant modulation in elephant vocalizations. Furthermore, by confirming the use of the elephant trunk in long distance rumble production, our findings provide an explanation for the extremely low formants in these calls, and may also indicate that formant lowering functions to increase call propagation distances in this species'

Fusion of secretory vesicles isolated from rat liver

Secretory vesicles isolated from rat liver were found to fuse after exposure to Ca2+. Vescle fusion is characterized by the occurrence of twinned vesicles with a continuous cleavage plane between two vesicles in freeze-fracture electron microscopy. The number of fused vesicles increases with increasing Ca2+-concentrations and is half maximal around 10–6 m. Other divalent cations (Ba2+, Sr2+, and Mg2+) were ineffective. Mg2+ inhibits Ca2+-induced fusion. Therefore, the fusion of secretory vesiclesin vitro is Ca2+ specific and exhibits properties similar to the exocytotic process of various secretory cells. Various substances affecting secretionin vivo (microtubular inhibitors, local anethetics, ionophores) were tested for their effect on membrane fusion in our system. The fusion of isolated secretory vesicles from liver was found to differ from that of pure phospholipid membranes in its temperature dependence, in its much lower requirement for Ca2+, and in its Ca2+-specificity. Chemical and enzymatic modifications of the vesicle membrane indicate that glycoproteins may account for these differences

Dynamic Analysis of Vascular Morphogenesis Using Transgenic Quail Embryos

Background: One of the least understood and most central questions confronting biologists is how initially simple clusters or sheet-like cell collectives can assemble into highly complex three-dimensional functional tissues and organs. Due to the limits of oxygen diffusion, blood vessels are an essential and ubiquitous presence in all amniote tissues and organs. Vasculogenesis, the de novo self-assembly of endothelial cell (EC) precursors into endothelial tubes, is the first step in blood vessel formation [1]. Static imaging and in vitro models are wholly inadequate to capture many aspects of vascular pattern formation in vivo, because vasculogenesis involves dynamic changes of the endothelial cells and of the forming blood vessels, in an embryo that is changing size and shape. Methodology/Principal Findings: We have generated Tie1 transgenic quail lines Tg(tie1:H2B-eYFP) that express H2B-eYFP in all of their endothelial cells which permit investigations into early embryonic vascular morphogenesis with unprecedented clarity and insight. By combining the power of molecular genetics with the elegance of dynamic imaging, we follow the precise patterning of endothelial cells in space and time. We show that during vasculogenesis within the vascular plexus, ECs move independently to form the rudiments of blood vessels, all while collectively moving with gastrulating tissues that flow toward the embryo midline. The aortae are a composite of somatic derived ECs forming its dorsal regions and the splanchnic derived ECs forming its ventral region. The ECs in the dorsal regions of the forming aortae exhibit variable mediolateral motions as they move rostrally; those in more ventral regions show significant lateral-to-medial movement as they course rostrally. Conclusions/Significance: The present results offer a powerful approach to the major challenge of studying the relative role(s) of the mechanical, molecular, and cellular mechanisms of vascular development. In past studies, the advantages of the molecular genetic tools available in mouse were counterbalanced by the limited experimental accessibility needed for imaging and perturbation studies. Avian embryos provide the needed accessibility, but few genetic resources. The creation of transgenic quail with labeled endothelia builds upon the important roles that avian embryos have played in previous studies of vascular development