High level of divergence of male-reproductive-tract proteins, between Drosophila melanogaster and its sibling species, D. simulans

We compared male-reproductive-tract polypeptides of Drosophila melanogaster and D. simulans by using two-dimensional gel electrophoresis. Approximately 64% of male-reproductive-tract polypeptides were identical between two randomly chosen isofemale lines from these two species, compared with 83% identity for third-instar imaginal wing-disc polypeptides. Qualitatively similar differences were found between reproductive tracts and imaginal discs when D. sechellia was compared with D. melanogaster and with D. simulans. When genic polymorphism was taken into account, approximately 10% of male- reproductive-tract polypeptides were apparently fixed for different alleles between D. melanogaster and D. simulans; this proportion is the same as that found for soluble enzymes by one-dimensional gel electrophoresis. Strikingly, approximately 20% of male-reproductive- tract polypeptides of either D. melanogaster or D. simulans had no detectable homologue in the other species. We propose that proteins of the Drosophila male reproductive tract may have diverged more extensively between species than have other types of proteins and that much of this divergence may involve large changes in levels of polypeptide expression.      

Microelectrode studies of the active Na transport pathway of frog skin

When the outer surface of short-circuited frog skin was penetrated with microelectrodes, stable negative potentials that averaged near -100 mV were recorded consistently, confirming the results of Nagel (W. Nagel. 1975. Abstracts of the 5th International Biophysics Congress, Copenhagen. P-147.). The appearance of these stable potentials, V(O), concurrent with the observations that (a) a high resistance outer barrier R(O) accounting for approximately 75 percent or more of the transcellular resistance of control skins had been penetrated and that (b) 10(-5) M amiloride and reduced [Na] outside caused the values of V(O) to increase towards means value near -130 mV while the values of percent R(O) increased to more than 90 percent. It was of relationships were the same as the values of E(1) observed in studies of the current-voltage relationships were the same as the values of E’(1) defined as the values of voltage at the inner barrier when the V(O) of the outer barrier was reduced to zero by voltage clamping of the skins. Accordingly, these data are interpreted to mean that the values of E(1), approximately 130 mV, represent the E(Na) of the sodium pump at the inner barrier. 2,4-DNP was observed to decrease the values of transepithelial voltage less than E(1) the V(O) was negative. These data can be interpreted with a simple electrical equivalent circuit of the active sodium transport pathway of the frog skin that includes the idea that the outer membrane behaves as an electrical rectifier for ion transport.     

Generalized Verlet Algorithm for Efficient Molecular Dynamics Simulations with Long-range Interactions

For the purpose of molecular dynamics simulations of large biopolymers we have developed a new method to accelerate the calculation of long-range pair interactions (e.g. Coulomb interaction). The algorithm introduces distance classes to schedule updates of non-bonding interactions and to avoid unnecessary computations of interactions between particles which are far apart. To minimize the error caused by the updating schedule, the Verlet integration scheme has been modified. The results of the method are compared to those of other approximation schemes as well as to results obtained by numerical integration without approximation. For simulation of a protein with 12 637 atoms our approximation scheme yields a reduction of computer time by a factor of seven. The approximation suggested can be implemented on sequential as well as on parallel computers. We describe an implementation on a (Transputer-based) MIMD machine with a systolic ring architecture.     

Fitting genetic mapping functions based on sperm typing: results for three chromosomal segments in cattle

Genetic mapping functions translate the observed recombination rate between two loci into the corresponding map distance in Morgan units. Different mapping functions give different weights to multiple crossing overs and therefore lead to different results. This points out that not every function is best suited to fit a data set. The data used in this study originated from 2214 sperm from 37 Norwegian bulls, which were genotyped for 11 markers. The optimal functions for the chromosomes 6, 23 and the sex chromosome of cattle were derived using the maximum likelihood method, the likelihood ratio test and empirical discriminant analysis. It became apparent that for each chromosome a different function fitted the data best. These were the function of Rao et al. (Human Heredity 1977, 27, 99–104) with p = 0·63 for chromosome 6, the function of Goldgar & Fain (American Journal of Human Genetics 1988, 43, 38–45) with c₀ = 0·42, c₁ = 0·47, c₂ = 0·07 and c₃ = 0·04 for chromosome 23 and the function of Felsenstein (Genetics 1979, 91, 769–75) with K = 0·23 for the sex chromosome. The well known functions of Haldane (Journal of Genetics 1919, 8, 299–309) and Kosambi (Annals of Eugenics 1944, 12, 172–5) were shown to be suboptimal in most cases. A function is said to be multilocus feasible if the evaluation of the probability of all possible recombination events does not lead to negative values. The optimal function for chromosome 23 turned out to be multilocus feasible, whereas the functions for chromosome 6 and the sex chromosome were not. The choice of the correct mapping function is shown to have a considerable impact in mapping studies, when double recombinations have to be taken into account. Since there is no unique best mapping function, it is argued that it might be useful to use a simple parametric mapping function (like the one of Felsenstein 1979) and to estimate the respective parameter specifically for a given data set.     

Direct intracardiac injection of umbilical cord-derived stromal cells and umbilical cord blood-derived endothelial cells for the treatment of ischemic cardiomyopathy

The development of new therapeutic strategies is necessary to reduce the worldwide social and economic impact of cardiovascular disease, which produces high rates of morbidity and mortality. A therapeutic option that has emerged in the last decade is cell therapy. The aim of this study was to compare the effect of transplanting human umbilical cord-derived stromal cells (UCSCs), human umbilical cord blood-derived endothelial cells (UCBECs) or a combination of these two cell types for the treatment of ischemic cardiomyopathy (IC) in a Wistar rat model. IC was induced by left coronary artery ligation, and baseline echocardiography was performed seven days later. Animals with a left ventricular ejection fraction (LVEF) of ≤40% were selected for the study. On the ninth day after IC was induced, the animals were randomized into the following experimental groups: UCSCs, UCBECs, UCSCs plus UCBECs, or vehicle (control). Thirty days after treatment, an echocardiographic analysis was performed, followed by euthanasia. The animals in all of the cell therapy groups, regardless of the cell type transplanted, had less collagen deposition in their heart tissue and demonstrated a significant improvement in myocardial function after IC. Furthermore, there was a trend of increasing numbers of blood vessels in the infarcted area. The median value of LVEF increased by 7.19% to 11.77%, whereas the control group decreased by 0.24%. These results suggest that UCSCs and UCBECs are promising cells for cellular cardiomyoplasty and can be an effective therapy for improving cardiac function following IC.     

Most nuclear systemic autoantigens are extremely disordered proteins: implications for the etiology of systemic autoimmunity

Patients with systemic autoimmune diseases usually produce high levels of antibodies to self-antigens (autoantigens). The repertoire of common autoantigens is remarkably limited, yet no readily understandable shared thread links these apparently diverse proteins. Using computer prediction algorithms, we have found that most nuclear systemic autoantigens are predicted to contain long regions of extreme structural disorder. Such disordered regions would generally make poor B cell epitopes and are predicted to be under-represented as potential T cell epitopes. Consideration of the potential role of protein disorder may give novel insights into the possible role of molecular mimicry in the pathogenesis of autoimmunity. The recognition of extreme autoantigen protein disorder has led us to an explicit model of epitope spreading that explains many of the paradoxical aspects of autoimmunity – in particular, the difficulty in identifying autoantigen-specific helper T cells that might collaborate with the B cells activated in systemic autoimmunity. The model also explains the experimentally observed breakdown of major histocompatibility complex (MHC) class specificity in peptides associated with the MHC II proteins of activated autoimmune B cells, and sheds light on the selection of particular T cell epitopes in autoimmunity. Finally, the model helps to rationalize the relative rarity of clinically significant autoimmunity despite the prevalence of low specificity/low avidity autoantibodies in normal individuals.     

EFFECTS OF HUMAN ACTIVITIES ON STRUCTURE AND COMPOSITION OF WOODY SPECIES OF THE NOKREK BIOSPHERE RESERVE OF MEGHALAYA,NORTHEAST INDIA

Aims Our study was conducted in the Nokrek Biosphere Reserve (NBR) in the Garo hills districts of Meghalaya, Northeast India. Our aim was to assess the effects of human activities on plant diversity,population structure and regeneration.Methods We selected a representative 1.2 hm2 stand in both the core and buffer zones of NBR. Structure and composition were determined by randomly sampling square quadrats, population structure was assessed by determining age structure, and regeneration was assessed by measuring densities of seedling, sapling and adult trees.Important findings More woody species were recorded from the core zone than the buffer zone (87 vs. 81 species), and there were a large number of tropical, temperate, and Sino-Himalayan, Burma-Malaysian and Malayan elements, primitive families and primitive genera. The trees were distributed in three distinct strata,canopy, subcanopy and sapling. Subcanopy and sapling layers had the highest species richness (81％ -88％ ). Lauraceae and Euphorbiaceae were the dominant families in terms of the number of species, and a large number of families were represented by single species. Most woody species (57 ％ - 79 ％ ) were contagiously distributed and had low frequency ( ＜ 20％ ). Although stand density was high in the buffer zone, its basal area was low compared to the stand in the core zone. Low similarity and high β-diversity indicate marked differences in species composition of the stands. Shannon diversity index was high in both the stands, while Simpson dominance index was low. The diameter-class distribution for dominant species revealed that the most had a large number of young individuals in their populations. Preponderance of tree seedlings, followed by a steep decline in population density of saplings and adult trees, indicated that the seedling to sapling stage was the most critical in the life cycle of the tree populations. Most species (42 ％ - 48 ％ ) had no regeneration,25 ％ - 35 ％ had good/fair regeneration, and the rest had poor regeneration or reoccurred as immigrants.     

On the calculation of binding free energies using continuum methods: Application to MHC class I protein-peptide interactions

This paper describes a methodology to calculate the binding free energy (ΔG) of a protein-ligand complex using a continuum model of the solvent. A formal thermodynamic cycle is used to decompose the binding free energy into electrostatic and non-electrostatic contributions. In this cycle, the reactants are discharged in water, associated as purely nonpolar entities, and the final complex is then recharged. The total electrostatic free energies of the protein, the ligand, and the complex in water are calculated with the finite difference Poisson-Boltzmann (FDPB) method. The nonpolar (hydrophobic) binding free energy is calculated using a free energy-surface area relationship, with a single alkane/water surface tension coefficient (γ_aw). The loss in backbone and side-chain configurational entropy upon binding is estimated and added to the electrostatic and the nonpolar components of ΔG. The methodology is applied to the binding of the murine MHC class I protein H-2K^b with three distinct peptides, and to the human MHC class I protein HLA-A2 in complex with five different peptides. Despite significant differences in the amino acid sequences of the different peptides, the experimental binding free energy differences (ΔΔG_exp) are quite small (<0.3 and <2.7 kcal/mol for the H-2K^b and HLA-A2 complexes, respectively). For each protein, the calculations are successful in reproducing a fairly small range of values for ΔΔG_calc (<4.4 and <5.2 kcal/mol, respectively) although the relative peptide binding affinities of H-2K^b and HLA-A2 are not reproduced. For all protein-peptide complexes that were treated, it was found that electrostatic interactions oppose binding whereas nonpolar interactions drive complex formation. The two types of interactions appear to be correlated in that larger nonpolar contributions to binding are generally opposed by increased electrostatic contributions favoring dissociation. The factors that drive the binding of peptides to MHC proteins are discussed in light of our results.     

In vivo uptake of a haem analogue Zn protoporphyrin IX by the human malaria parasite P. falciparum‐infected red blood cells

The cellular traffic of haem during the development of the human malaria parasite Plasmodium falciparum, through the stages R (ring), T (trophozoite) and S (schizonts), was investigated within RBC (red blood cells). When Plasmodium cultures were incubated with a fluorescent haem analogue, ZnPPIX (Zn protoporphyrin IX) the probe was seen at the cytoplasm (R stage), and the vesicle‐like structure distribution pattern was more evident at T and S stages. The temporal sequence of ZnPPIX uptake byP. falciparum‐infected erythrocytes shows that at R and S stages, a time‐increase acquisition of the porphyrin reaches the maximum fluorescence distribution after 60 min; in contrast, at the T stage, the maximum occurs after 120 min of ZnPPIX uptake. The difference in time‐increase acquisition of the porphyrin is in agreement with a maximum activity of haem uptake at the T stage. To gain insights into haem metabolism, recombinant PfHO (P. falciparum haem oxygenase) was expressed, and the conversion of haem into BV (biliverdin) was detected. These findings point out that, in addition to haemozoin formation, the malaria parasite P. falciparum has evolved two distinct mechanisms for dealing with haem toxicity, namely, the uptake of haem into a cellular compartment where haemozoin is formed and HO activity. However, the low Plasmodium HO activity detected reveals that the enzyme appears to be a very inefficient way to scavenge the haem compared with the Plasmodium ability to uptake the haem analogue ZnPPIX and delivering it to the food vacuole.