Genome-wide meta-analysis for rheumatoid arthritis

Meta-analysis is being increasingly used as a tool for integrating data from different studies of complex phenotypes, because the power of any one study to identify causal loci is limited. We applied a novel meta-analytical approach (Loesgen et al. in Genet Epidemiol 21(Suppl 1):S142–S147, 2001) in compiling results from four studies of rheumatoid arthritis in Caucasians including two studies from NARAC (Jawaheer et al. in Am J Hum Genet 68:927–936, 2001; Jawaheer et al. in Arthritis Rheum 48:906–916, 2003), one study from the UK (MacKay et al. in Arthritis Rheum 46:632–639, 2001) and one from France (Cornelis et al. in Proc Natl Acad Sci USA 95:10746–10750, 1998). For each study, we obtained NPL scores by performing interval mapping (2 cM intervals) using GeneHunter2 (Kruglyak et al. in Am J Hum Genet 58:1347–1363, 1996; Markianos et al. in Am J Hum Genet 68:963–977, 2001). The marker maps differed among the three consortium groups, therefore, the marker maps were aligned after the interval mapping was completed and the NPL scores that were within 1 cM of each other were combined using the method of Loesgen et al. (Genet Epidemiol 21(Suppl 1):S142–S147, 2001) by calculating the weighted average of the NPL score. This approach avoids some problems in analysis encountered by using GeneHunter2 when some markers in the sample are not genotyped. This procedure provided marginal evidence (P<0.05) of linkage on chromosome 1, 2, 5 and 18, strong evidence (P<0.01) on chromosomes 8 and 16, and overwhelming evidence in the HLA region of chromosome 6.     

Bone tissue,lyophilized and stored at room temperature for 15 days or more,is not capable of transmitting HIV,HCV or HBV

Over the past 57 years, 17 recipients of frozen bone have been infected with: HIV (Centers for Disease Control and Prevention in Morb Mortal Wkly Rep MMWR 37(39):597–599, 1988; Li et al. in J Formos Med Assoc 100(5):350–351, 2001; Simonds et al. in NEJM 326(11):726–732, 1992; Schratt et al. in Unfallchirurg 99(9):679–684, 1996); HCV (Eggen and Nordbo in NEJM 326(6):411, 1992; Conrad et al. in J Bone Joint Surg Am 77:214–224, 1995; Trotter in J Bone Joint Surg Am 851(11):2215–2217, 2003; Tugwell et al. in Ann of Internal Med 143(9):648–654, 2005); or HBV (Shutkin in J Bone Joint Surg Am 36:160–162, 1954). However, bone, lyophilized and stored at room temperature, has never transmitted these viral diseases. A literature review was undertaken to determine whether there is any evidence that lyophilized bone is capable of transmitting HIV, HCV and HBV.     

Chaperone receptors: guiding proteins to intracellular compartments

Despite mitochondria and chloroplasts having their own genome, 99% of mitochondrial proteins (Rehling et al., Nat Rev Mol Cell Biol 5:519–530, 2004) and more than 95% of chloroplast proteins (Soll, Curr Opin Plant Biol 5:529–535, 2002) are encoded by nuclear DNA, synthesised in the cytosol and imported post-translationally. Protein targeting to these organelles depends on cytosolic targeting factors, which bind to the precursor, and then interact with membrane receptors to deliver the precursor into a translocase. The molecular chaperones Hsp70 and Hsp90 have been widely implicated in protein targeting to mitochondria and chloroplasts, and receptors capable of recognising these chaperones have been identified at the surface of both these organelles (Schlegel et al., Mol Biol Evol 24:2763–2774, 2007). The role of these chaperone receptors is not fully understood, but they have been shown to increase the efficiency of protein targeting (Young et al., Cell 112:41–50, 2003; Qbadou et al., EMBO J 25:1836–1847, 2006). Whether these receptors contribute to the specificity of targeting is less clear. A class of chaperone receptors bearing tetratricopeptide repeat domains is able to specifically bind the highly conserved C terminus of Hsp70 and/or Hsp90. Interestingly, at least of one these chaperone receptors can be found on each organelle (Schlegel et al., Mol Biol Evol 24:2763–2774, 2007), which suggests a universal role in protein targeting for these chaperone receptors. This review will investigate the role that chaperone receptors play in targeting efficiency and specificity, as well as examining recent in silico approaches to find novel chaperone receptors.     

The NMR restraints grid at BMRB for 5,266 protein and nucleic acid PDB entries

Several pilot experiments have indicated that improvements in older NMR structures can be expected by applying modern software and new protocols (Nabuurs et al. in Proteins 55:483–186, 2004; Nederveen et al. in Proteins 59:662–672, 2005; Saccenti and Rosato in J Biomol NMR 40:251–261, 2008). A recent large scale X-ray study also has shown that modern software can significantly improve the quality of X-ray structures that were deposited more than a few years ago (Joosten et al. in J. Appl Crystallogr 42:376–384, 2009; Sanderson in Nature 459:1038–1039, 2009). Recalculation of three-dimensional coordinates requires that the original experimental data are available and complete, and are semantically and syntactically correct, or are at least correct enough to be reconstructed. For multiple reasons, including a lack of standards, the heterogeneity of the experimental data and the many NMR experiment types, it has not been practical to parse a large proportion of the originally deposited NMR experimental data files related to protein NMR structures. This has made impractical the automatic recalculation, and thus improvement, of the three dimensional coordinates of these structures. We here describe a large-scale international collaborative effort to make all deposited experimental NMR data semantically and syntactically homogeneous, and thus useful for further research. A total of 4,014 out of 5,266 entries were ‘cleaned’ in this process. For 1,387 entries, human intervention was needed. Continuous efforts in automating the parsing of both old, and newly deposited files is steadily decreasing this fraction. The cleaned data files are available from the NMR restraints grid at .     

1H, 13C and 15N NMR assignments of the C1A and C1B subdomains of PKC-delta

The Protein Kinase C family of enzymes is a group of serine/threonine kinases that play central roles in cell-cycle regulation, development and cancer. A key step in the activation of PKC is translocation to membranes and binding of membrane-associated activators including diacylglycerol (DAG). Interaction of novel and conventional isotypes of PKC with DAG and phorbol esters occurs through the two C1 regulatory domains (C1A and C1B), which exhibit distinct ligand binding selectivity that likely controls enzyme activation by different co-activators. PKC has also been implicated in physiological responses to alcohol consumption and it has been proposed that PKCα (Slater et al. J Biol Chem 272(10):6167–6173, 1997; Slater et al. Biochemistry 43(23):7601–7609, 2004), PKCε (Das et al. Biochem J 421(3):405–413, 2009) and PKCδ (Das et al. J Biol Chem 279(36):37964–37972, 2004; Das et al. Protein Sci 15(9):2107–2119, 2006) contain specific alcohol-binding sites in their C1 domains. We are interested in understanding how ethanol affects signal transduction processes through its affects on the structure and function of the C1 domains of PKC. Here we present the ¹H, ¹⁵N and ¹³C NMR chemical shift assignments for the Rattus norvegicus PKCδ C1A and C1B proteins.     

The synaptonemal complex and meiotic recombination in humans: new approaches to old questions

Meiotic prophase serves as an arena for the interplay of two important cellular activities, meiotic recombination and synapsis of homologous chromosomes. Synapsis is mediated by the synaptonemal complex (SC), originally characterized as a structure linked to pairing of meiotic chromosomes (Moses (1958) J Biophys Biochem Cytol 4:633–638). In 1975, the first electron micrographs of human pachytene stage SCs were presented (Moses et al. (1975) Science 187:363–365) and over the next 15 years the importance of the SC to normal meiotic progression in human males and females was established (Jhanwar and Chaganti (1980) Hum Genet 54:405–408; Pathak and Elder (1980) Hum Genet 54:171–175; Solari (1980) Chromosoma 81:315–337; Speed (1984) Hum Genet 66:176–180; Wallace and Hulten (1985) Ann Hum Genet 49(Pt 3):215–226). Further, these studies made it clear that abnormalities in the assembly or maintenance of the SC were an important contributor to human infertility (Chaganti et al. (1980) Am J Hum Genet 32:833–848; Vidal et al. (1982) Hum Genet 60:301–304; Bojko (1983) Carlsberg Res Commun 48:285–305; Bojko (1985) Carlsberg Res Commun 50:43–72; Templado et al. (1984) Hum Genet 67:162–165; Navarro et al. (1986) Hum Reprod 1:523–527; Garcia et al. (1989) Hum Genet 2:147–53). However, the utility of these early studies was limited by lack of information on the structural composition of the SC and the identity of other SC-associated proteins. Fortunately, studies of the past 15 years have gone a long way toward remedying this problem. In this minireview, we highlight the most important of these advances as they pertain to human meiosis, focusing on temporal aspects of SC assembly, the relationship between the SC and meiotic recombination, and the contribution of SC abnormalities to human infertility.The synaptonemal complex–50 years     

ROS-Mediated ABA Signaling

In plants, reactive oxygen species (ROS) are short-lived molecules produced through various cellular mechanisms in response to biotic and abiotic stimuli. ROS function as second messengers for hormone signaling, development, oxygen deprivation, programmed cell death, and plant–pathogen interactions. Recent research on ROS-mediated responses has produced stimulating findings such as the specific sources of ROS production, molecular elements that work in ROS-mediated signaling and homeostasis, and a ROS-regulated gene network (Neill et al., Curr Opin Plant Biol 5:388–395, 2002a; Apel and Hirt, Annu Rev Plant Biol 55:373–399, 2004; Mittler et al., Trends Plant Sci 9:490–498, 2004; Mori and Schroeder, Plant Physiol 135:702–708, 2004; Kwak et al., Plant Physiol 141:323–329, 2006; Torres et al., Plant Physiol 141:373–378, 2006; Miller et al., Physiol Plant 133:481–489, 2008). In this review, we highlight new discoveries in ROS-mediated abscisic acid (ABA) signaling. Drs. Daeshik Cho and June M. Kwak are the corresponding authors for this paper.     

Adding Adhesion to a Chemical Signaling Model for Somite Formation

Somites are condensations of mesodermal cells that form along the two sides of the neural tube during early vertebrate development. They are one of the first instances of a periodic pattern, and give rise to repeated structures such as the vertebrae. A number of theories for the mechanisms underpinning somite formation have been proposed. For example, in the “clock and wavefront” model (Cooke and Zeeman in J. Theor. Biol. 58:455–476, 1976), a cellular oscillator coupled to a determination wave progressing along the anterior-posterior axis serves to group cells into a presumptive somite. More recently, a chemical signaling model has been developed and analyzed by Maini and coworkers (Collier et al. in J. Theor. Biol. 207:305–316, 2000; Schnell et al. in C. R. Biol. 325:179–189, 2002; McInerney et al. in Math. Med. Biol. 21:85–113, 2004), with equations for two chemical regulators with entrained dynamics. One of the chemicals is identified as a somitic factor, which is assumed to translate into a pattern of cellular aggregations via its effect on cell–cell adhesion. Here, the authors propose an extension to this model that includes an explicit equation for an adhesive cell population. They represent cell adhesion via an integral over the sensing region of the cell, based on a model developed previously for adhesion driven cell sorting (Armstrong et al. in J. Theor. Biol. 243:98–113, 2006). The expanded model is able to reproduce the observed pattern of cellular aggregates, but only under certain parameter restrictions. This provides a fuller understanding of the conditions required for the chemical model to be applicable. Moreover, a further extension of the model to include separate subpopulations of cells is able to reproduce the observed differentiation of the somite into separate anterior and posterior halves. N.J. Armstrong was supported by a Doctoral Training Account Studentship from EPSRC. K.J. Painter and J.A. Sherratt were supported in part by Integrative Cancer Biology Program Grant CA113004 from the US National Institute of Health and in part by BBSRC grant BB/D019621/1 for the Centre for Systems Biology at Edinburgh.     

A Computational Analysis of Localized Ca<Superscript>2+</Superscript>-Dynamics Generated by Heterogeneous Release Sites

We investigate the role of heterogeneous expression of IP₃R and RyR in generating diverse elementary Ca²⁺ signals. It has been shown empirically (Wojcikiewicz and Luo in Mol. Pharmacol. 53(4):656–662, 1998; Newton et al. in J. Biol. Chem. 269(46):28613–28619, 1994; Smedt et al. in Biochem. J. 322(Pt. 2):575–583, 1997) that tissues express various proportions of IP₃ and RyR isoforms and this expression is dynamically regulated (Parrington et al. in Dev. Biol. 203(2):451–461, 1998; Fissore et al. in Biol. Reprod. 60(1):49–57, 1999; Tovey et al. in J. Cell Sci. 114(Pt. 22):3979–3989, 2001). Although many previous theoretical studies have investigated the dynamics of localized calcium release sites (Swillens et al. in Proc. Natl. Acad. Sci. U.S.A. 96(24):13750–13755, 1999; Shuai and Jung in Proc. Natl. Acad. Sci. U.S.A. 100(2):506–510, 2003a; Shuai and Jung in Phys. Rev. E, Stat. Nonlinear Soft Matter Phys. 67(3 Pt. 1):031905, 2003b; Thul and Falcke in Biophys. J. 86(5):2660–2673, 2004; DeRemigio and Smith in Cell Calcium 38(2):73–86, 2005; Nguyen et al. in Bull. Math. Biol. 67(3):393–432, 2005), so far all such studies focused on release sites consisting of identical channel types. We have extended an existing mathematical model (Nguyen et al. in Bull. Math. Biol. 67(3):393–432, 2005) to release sites with two (or more) receptor types, each with its distinct channel kinetics. Mathematically, the release site is represented by a transition probability matrix for a collection of nonidentical stochastically gating channels coupled through a shared Ca²⁺ domain. We demonstrate that under certain conditions a previously defined mean-field approximation of the coupling strength does not accurately reproduce the release site dynamics. We develop a novel approximation and establish that its performance in these instances is superior. We use this mathematical framework to study the effect of heterogeneity in the Ca²⁺-regulation of two colocalized channel types on the release site dynamics. We consider release sites consisting of channels with both Ca²⁺-activation and inactivation (“four-state channels”) and channels with Ca²⁺-activation only (“two-state channels”) and show that for the appropriate parameter values, synchronous channel openings within a release site with any proportion of two-state to four-state channels are possible, however, the larger the proportion of two-state channels, the more sensitive the dynamics are to the exact spatial positioning of the channels and the distance between channels. Specifically, the clustering of even a small number of two-state channels interferes with puff/spark termination and increases puff durations or leads to a tonic response.     

The heme environment of mouse neuroglobin: histidine imidazole plane orientations obtained from solution NMR and EPR spectroscopy as compared with X-ray crystallography

The ¹H NMR chemical shifts of the heme methyl groups of the ferriheme complex of metneuroglobin (Du et al. in J. Am. Chem. Soc. 125:8080–8081, 2003) predict orientations of the axial histidine ligands (Shokhirev and Walker in J. Biol. Inorg. Chem. 3:581–594, 1998) that are not consistent with the X-ray data (Vallone et al. in Proteins Struct. Funct. Bioinf. 56:85–94, 2004), and the EPR spectrum (Vinck et al. in J. Am. Chem. Soc. 126:4516–4517, 2004) is only marginally consistent with these data. The reasons for these inconsistencies appear to be rooted in the high degree of aqueous solution exposure of the heme group and the fact that there are no strong hydrogen-bond acceptors for the histidine imidazole N–H protons provided by the protein. Similar inconsistencies may exist for other water-soluble heme proteins, and ¹H NMR spectroscopy provides a simple means to verify whether the solution structure of the heme center is the same as or different from that in the crystalline state.     

Temperature controls a latitudinal gradient in the proportion of watershed nitrogen exported to coastal ecosystems

Increased export of biologically available nitrogen (N) to the coastal zone is strongly linked to eutrophication, which is a major problem in coastal marine ecosystems (NRC (2000) Clean Coastal Waters: Understanding and Reducing the Effects of Nutrient Pollution. National Academy Press, Washington, DC; Bricker et al. (1999) National Estuarine Eutrophication Assessment. Effects of nutrient enrichment in the nation’s estuaries. NOAA-NOS Special Projects Office, Silver Spring, MD). However, not all of the nitrogen input to a watershed is exported to the coast (Howarth et al. (1996) Biogeochemistry 35:75–139; Jordan and Weller (1996) Bioscience 46:655–664). Global estimates of nitrogen export to coasts have been taken to be 25% of watershed input, based largely on northeastern U.S. observations (Galloway et al. (2004) Biogeochemistry 70:153–226; Boyer et al. (2006) Global Biogeochem Cycle 20:Art. No. GB1S91). We applied the N budgeting methodology developed for the International SCOPE Nitrogen project (Howarth et al. (1996) Biogeochemistry 35:75–139; Boyer et al. (2002) Biogeochemistry 57:137–169) to 12 watersheds in the southeastern U.S., and compared them with estimates of N export for 16 watersheds in the northeastern U.S. (Boyer et al. (2002) Biogeochemistry 57:137–169). In southeastern watersheds, average N export was only 9% of input, suggesting the need for downward revision of global estimates. The difference between northern and southern watersheds is not a function of the absolute value of N inputs, which spanned a comparable range and were positively related to export in both cases. Rather, the proportion of N exported was significantly related to average watershed temperature (% N export = 58.41 e^{−0.11 * temperature}; R ² = 0.76), with lower proportionate nitrogen export in warmer watersheds. In addition, we identified a threshold in proportionate N export at 38°N latitude that corresponds to a reported breakpoint in the rate of denitrification at 10–12°C. We hypothesize that temperature, by regulating denitrification, results in increased proportionate N export at higher latitudes. Regardless of the mechanism, these observations suggest that temperature increases associated with future climate change may well reduce the amount of nitrogen that reaches estuaries, which will have implications for coastal eutrophication.     

Endocrine modulation of a pheromone-responsive gene in the honey bee brain

Pheromones cause dramatic changes in behavior and physiology, and are critical for honey bee colony organization. Queen mandibular pheromone (QMP) regulates multiple behaviors in worker bees (Slessor et al. in J Chem Ecol 31(11):2731–2745, 2005). We also identified genes whose brain expression levels were altered by exposure to QMP (Grozinger et al. in Proc Natl Acad Sci USA 100(Suppl 2):14519–14525, 2003). Krüppel-homolog 1 (Kr-h1) RNA levels were significantly downregulated by QMP, and were higher in foragers than in nurses (Whitfield et al. in Science 302(5643):296–299, 2003). Here we report on results of behavioral and pharmacological experiments that characterize factors regulating expression of Kr-h1. Foragers have higher brain levels of Kr-h1 than in-hive bees, regardless of age and pheromone exposure. Furthermore, forager Kr-h1 levels were not affected by QMP. Since the onset of foraging is caused, in part, by increasing juvenile hormone blood titers and brain octopamine levels, we investigated the effects of octopamine and methoprene (a juvenile hormone analog) on Kr-h1 expression. Methoprene produced a marginal (not significant) increase in Kr-h1 expression, but Kr-h1 brain levels in methoprene-treated bees were no longer downregulated by QMP. Octopamine did not modulate Kr-h1 expression. Our results demonstrate that the gene expression response to QMP is not hard-wired in the brain but is instead dependent on worker behavioral state.     

Use of weighted reference panels based on empirical estimates of ancestry for capturing untyped variation

Many association methods use a subset of genotyped single nucleotide polymorphisms (SNPs) to capture or infer genotypes at other untyped SNPs. We and others previously showed that tag SNPs selected to capture common variation using data from The International HapMap Consortium (Nature 437:1299–1320, 2005), The International HapMap Consortium (Nature 449:851–861, 2007) could also capture variation in populations of similar ancestry to HapMap reference populations (de Bakker et al. in Nat Genet 38:1298–1303, 2006; González-Neira et al. in Genome Res 16:323–330, 2006; Montpetit et al. in PLoS Genet 2:282–290, 2006; Mueller et al. in Am J Hum Genet 76:387–398, 2005). To capture variation in admixed populations or populations less similar to HapMap panels, a “cosmopolitan approach,” in which all samples from HapMap are used as a single reference panel, was proposed. Here we refine this suggestion and show that use of a “weighted reference panel,” constructed based on empirical estimates of ancestry in the target population (relative to available reference panels), is more efficient than the cosmopolitan approach. Weighted reference panels capture, on average, only slightly fewer common variants (minor allele frequency > 5%) than the cosmopolitan approach (mean r ² = 0.977 vs. 0.989, 94.5% variation captured vs. 96.8% at r ² > 0.8), across the five populations of the Multiethnic Cohort, but entail approximately 25% fewer tag SNPs per panel (average 538 vs. 718). These results extend a recent study in two Indian populations (Pemberton et al. in Ann Hum Genet 72:535–546, 2008). Weighted reference panels are potentially useful for both the selection of tag SNPs in diverse populations and perhaps in the design of reference panels for imputation of untyped genotypes in genome-wide association studies in admixed populations. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comment arising from a paper by Wittmer et al.: hypothesis testing for top-down and bottom-up effects in woodland caribou population dynamics

Conservation strategies for populations of woodland caribou Rangifer tarandus caribou frequently emphasize the importance of predator–prey relationships and the availability of lichen-rich late seral forests, yet the importance of summer diet and forage availability to woodland caribou survival is poorly understood. In a recent article, Wittmer et al. (Can J Zool 83:407–418, 2005b) concluded that woodland caribou in British Columbia were declining as a consequence of increased predation that was facilitated by habitat alteration. Their conclusion is consistent with the findings of other authors who have suggested that predation is the most important proximal factor limiting woodland caribou populations (Bergerud and Elliot in Can J Zool 64:1515–1529, 1986; Edmonds in Can J Zool 66:817–826, 1988; Rettie and Messier in Can J Zool 76:251–259, 1998; Hayes et al. in Wildl Monogr 152:1–35, 2003). Wittmer et al. (Can J Zool 83:407–418, 2005b) presented three alternative, contrasting hypotheses for caribou decline that differed in terms of predicted differences in instantaneous rates of increase, pregnancy rates, causes of mortality, and seasonal vulnerability to mortality (Table 1, p 258). These authors rejected the hypotheses that food or an interaction between food and predation was responsible for observed declines in caribou populations; however, the use of pregnancy rate, mortality season and cause of mortality to contrast the alternative hypotheses is problematic. We argue here that the data employed in their study were insufficient to properly evaluate a predation-sensitive foraging hypothesis for caribou decline. Empirical data on seasonal forage availability and quality and plane of nutrition of caribou would be required to test the competing hypotheses. We suggest that methodological limitations in studies of woodland caribou population dynamics prohibit proper evaluation of the mechanism of caribou population declines and fail to elucidate potential interactions between top-down and bottom-up effects on populations. An erratum to this article can be found at     

Mechanisms of age-specific regulation of dopamine metabolism by juvenile hormone and 20-hydroxyecdysone in <Emphasis Type="Italic">Drosophila</Emphasis> females

20-hydroxyecdysone (20E) and the juvenile hormone (JH) have an age-specific effect on total dopamine (DA) content in Drosophila (Gruntenko and Rauschenbach 2008). Earlier we studied the mechanism of influence of 20E and JH on DA metabolism in young females (Rauschenbach et al. in J Insect Physiol 53:587–591, 2007a: Arch Insect Biochem Physiol 65:95–102, 2008a; Gruntenko et al. in Arch Insect Biochem Physiol 72:263–269, 2009). Here we investigate the effects of 20E and JH on the activities of the alkaline phosphatase (ALP), tyrosine hydroxylase (TH) and DA-dependent arylalkylamine N-acetyltransferase (AANAT) in mature females of wild type D. virilis under normal conditions and under heat stress (38°C). 20E feeding of the flies led to a substantial decrease in ALP and TH activities and to an increase in AANAT activity in mature females. JH application resulted in an increasing of ALP and TH activities, but did not influence AANAT activity in mature females. A rise in JH and 20E levels was found to change ALP and TH stress reactivities. Mechanisms of age-specific regulation of DA level by 20E and JH in Drosophila females are discussed.     

Structure and dynamics of the antifungal molecules Syringotoxin-B and Syringopeptin-25A from molecular dynamics simulation

The phytopathogen Pseudomonas syringae pv. syringae produces toxic cyclic lipodepsipeptides (CLPs): nona-peptides and syringopeptins. All CLPs inhibit the growth of many fungal species, including human pathogens, although different fungi display different degrees of sensitivity. The best studied CLPs are Syringomycin-E (SR-E), Syringotoxin-B (ST-B) and Syringopeptin-25A (SP-25A). Their biological activity is affected by membrane composition and their structural differences. We previously (Matyus et al. in Eur Biophys J 35:459–467, 2006) reported the molecular features and structural preferences of SR-E in water and octane environments. Here we investigate in atomic detail the molecular features of the two other main CLP components, ST-B and SP-25A, in water and octane by 200 ns molecular dynamics simulations (MD), using distance restraints derived from NMR NOE data (Ballio et al. in Eur J Biochem 234:747–758, 1995). We have obtained three-dimensional models of ST-B and SP-25A CLPs in different environments. These models can now be used as a basis to investigate the interactions of ST-B and SP-25A with lipid membranes an important further step towards a better understanding of the antifungal and antibacterial activity of these peptides. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.