Computing the local pressure in molecular dynamics simulations

Computer simulations of inhomogeneous soft matter systems often require accurate methods for computing the local pressure. We present a simple derivation, based on the virial relation, of two equivalent expressions for the local (atomistic) pressure in a molecular dynamics simulation. One of these expressions, previously derived by other authors via a different route, involves summation over interactions between particles within the region of interest; the other involves summation over interactions across the boundary of the region of interest. We illustrate our derivation using simulations of a simple osmotic system; both expressions produce accurate results even when the region of interest over which the pressure is measured is very small.Comment: 11 pages, 4 figure

Venom gland transcriptomes of two elapid snakes (Bungarus multicinctus and Naja atra) and evolution of toxin genes

<p>Abstract</p> <p>Background</p> <p>Kraits (genus <it>Bungarus</it>) and cobras (genus <it>Naja</it>) are two representative toxic genera of elapids in the old world. Although they are closely related genera and both of their venoms are very toxic, the compositions of their venoms are very different. To unveil their detailed venoms and their evolutionary patterns, we constructed venom gland cDNA libraries and genomic bacterial artificial chromosome (BAC) libraries for <it>Bungarus multicinctus </it>and <it>Naja atra</it>, respectively. We sequenced about 1500 cDNA clones for each of the venom cDNA libraries and screened BAC libraries of the two snakes by blot analysis using four kinds of toxin probes; <it>i.e</it>., three-finger toxin (3FTx), phospholipase A2 (PLA2), kunitz-type protease inhibitor (Kunitz), and natriuretic peptide (NP).</p> <p>Results</p> <p>In total, 1092 valid expressed sequences tags (ESTs) for <it>B. multicinctus </it>and 1166 ESTs for <it>N. atra </it>were generated. About 70% of these ESTs can be annotated as snake toxin transcripts. 3FTx (64.5%) and <it>β </it>bungarotoxin (25.1%) comprise the main toxin classes in <it>B. multicinctus</it>, while 3FTx (95.8%) is the dominant toxin in <it>N. atra</it>. We also observed several less abundant venom families in <it>B. multicinctus </it>and <it>N. atra</it>, such as PLA2, C-type lectins, and Kunitz. Peculiarly a cluster of NP precursors with tandem NPs was detected in <it>B. multicinctus</it>. A total of 71 positive toxin BAC clones in <it>B. multicinctus </it>and <it>N. atra </it>were identified using four kinds of toxin probes (3FTx, PLA2, Kunitz, and NP), among which 39 3FTx-postive BACs were sequenced to reveal gene structures of 3FTx toxin genes.</p> <p>Conclusions</p> <p>Based on the toxin ESTs and 3FTx gene sequences, the major components of <it>B. multicinctus </it>venom transcriptome are neurotoxins, including long chain alpha neurotoxins (<it>α</it>-ntx) and the recently originated <it>β </it>bungarotoxin, whereas the <it>N. atra </it>venom transcriptome mainly contains 3FTxs with cytotoxicity and neurotoxicity (short chain <it>α</it>-ntx). The data also revealed that tandem duplications contributed the most to the expansion of toxin multigene families. Analysis of nonsynonymous to synonymous nucleotide substitution rate ratios (<it>dN</it>/<it>dS</it>) indicates that not only multigene toxin families but also other less abundant toxins might have been under rapid diversifying evolution.</p

Discovery of a Distinct Superfamily of Kunitz-Type Toxin (KTT) from Tarantulas

BACKGROUND: Kuntiz-type toxins (KTTs) have been found in the venom of animals such as snake, cone snail and sea anemone. The main ancestral function of Kunitz-type proteins was the inhibition of a diverse array of serine proteases, while toxic activities (such as ion-channel blocking) were developed under a variety of Darwinian selection pressures. How new functions were grafted onto an old protein scaffold and what effect Darwinian selection pressures had on KTT evolution remains a puzzle. PRINCIPAL FINDINGS: Here we report the presence of a new superfamily of ktts in spiders (TARANTULAS: Ornithoctonus huwena and Ornithoctonus hainana), which share low sequence similarity to known KTTs and is clustered in a distinct clade in the phylogenetic tree of KTT evolution. The representative molecule of spider KTTs, HWTX-XI, purified from the venom of O. huwena, is a bi-functional protein which is a very potent trypsin inhibitor (about 30-fold more strong than BPTI) as well as a weak Kv1.1 potassium channel blocker. Structural analysis of HWTX-XI in 3-D by NMR together with comparative function analysis of 18 expressed mutants of this toxin revealed two separate sites, corresponding to these two activities, located on the two ends of the cone-shape molecule of HWTX-XI. Comparison of non-synonymous/synonymous mutation ratios (omega) for each site in spider and snake KTTs, as well as PBTI like body Kunitz proteins revealed high Darwinian selection pressure on the binding sites for Kv channels and serine proteases in snake, while only on the proteases in spider and none detected in body proteins, suggesting different rates and patterns of evolution among them. The results also revealed a series of key events in the history of spider KTT evolution, including the formation of a novel KTT family (named sub-Kuntiz-type toxins) derived from the ancestral native KTTs with the loss of the second disulfide bridge accompanied by several dramatic sequence modifications. CONCLUSIONS/SIGNIFICANCE: These finding illustrate that the two activity sites of Kunitz-type toxins are functionally and evolutionally independent and provide new insights into effects of Darwinian selection pressures on KTT evolution, and mechanisms by which new functions can be grafted onto old protein scaffolds

A Family of Diverse Kunitz Inhibitors from Echinococcus granulosus Potentially Involved in Host-Parasite Cross-Talk

The cestode Echinococcus granulosus, the agent of hydatidosis/echinococcosis, is remarkably well adapted to its definitive host. However, the molecular mechanisms underlying the successful establishment of larval worms (protoscoleces) in the dog duodenum are unknown. With the aim of identifying molecules participating in the E. granulosus-dog cross-talk, we surveyed the transcriptomes of protoscoleces and protoscoleces treated with pepsin at pH 2. This analysis identified a multigene family of secreted monodomain Kunitz proteins associated mostly with pepsin/H+-treated worms, suggesting that they play a role at the onset of infection. We present the relevant molecular features of eight members of the E. granulosus Kunitz family (EgKU-1 – EgKU-8). Although diverse, the family includes three pairs of close paralogs (EgKU-1/EgKU-4; EgKU-3/EgKU-8; EgKU-6/EgKU-7), which would be the products of recent gene duplications. In addition, we describe the purification of EgKU-1 and EgKU-8 from larval worms, and provide data indicating that some members of the family (notably, EgKU-3 and EgKU-8) are secreted by protoscoleces. Detailed kinetic studies with native EgKU-1 and EgKU-8 highlighted their functional diversity. Like most monodomain Kunitz proteins, EgKU-8 behaved as a slow, tight-binding inhibitor of serine proteases, with global inhibition constants (KI*) versus trypsins in the picomolar range. In sharp contrast, EgKU-1 did not inhibit any of the assayed peptidases. Interestingly, molecular modeling revealed structural elements associated with activity in Kunitz cation-channel blockers. We propose that this family of inhibitors has the potential to act at the E. granulosus-dog interface and interfere with host physiological processes at the initial stages of infection

Action de l’eau oxygénée sur les bases puriques et pyrimidiques, leurs déoxyribonucléotides et sur l’acide désoxyribonucléique

L’action de l’eau oxygénée sur l’ADN et ses constituants a été étudiée.L’eau oxygénée dégrade les hétérocycles des bases puriques et surtout pyrimidiques. La vitesse avec laquelle les hétérocycles des bases sont dégradées par H202 est sensiblement la même, qu’il s’agisse des bases libres ou des bases des déoxyribonucléotides; elle est beaucoup plus faible que celle que l’on observe lors de la dégradation de ces composés par le peroxyde succinique.L’eau oxygénée provoque en outre la rupture de liaisons déoxyribose-base, détachant les bases et des nucléotides libres et de l’ADN, ainsi que la rupture de liaisons déoxyribose-phosphate, qui se traduit par une diminution importante du poids moléculaire de l’ADN traité par H202.La comparaison du mode d’action de H202 sur les liaisons déoxyribose-base et déoxyribose-phosphate de l’ADN, avec celle d’autres agents qui affectent également ces liaisons, indiquerait que H202 agit suivant un mécanisme particulier