Gramicidin A channel as a test ground for molecular dynamics force fields

We use the well-known structural and functional properties of the gramicidin A channel to test the appropriateness of force fields commonly used in molecular dynamics (MD) simulations of ion channels. For this purpose, the high-resolution structure of the gramicidin A dimer is embedded in a dimyristoylphosphatidylcholine bilayer, and the potential of mean force of a K(+) ion is calculated along the channel axis using the umbrella sampling method. Calculations are performed using two of the most common force fields in MD simulations: CHARMM and GROMACS. Both force fields lead to large central barriers for K(+) ion permeation, that are substantially higher than those deduced from the physiological data by inverse methods. In long MD simulations lasting over 60 ns, several ions are observed to enter the binding site but none of them crossed the channel despite the presence of a large driving field. The present results, taken together with many earlier studies, highlights the shortcomings of the standard force fields used in MD simulations of ion channels and calls for construction of more appropriate force fields for this purpose.     

The fixation probability of a beneficial allele in a population dividing by binary fission

We derive formulae for the fixation probability, P, of a rare benefical allele segregating in a population of fixed size which reproduces by binary fission, in terms of the selection coefficient for the beneficial allele, s. We find that an earlier result P 4s does not depend on the assumption of binary fission, but depends on an assumption about the ordering of events in the life cycle. We find that P 2s for mutations occurring during chromosome replication and P 2.8s for mutations occurring at random times between replication events.     

General models of multilocus evolution

In 1991, Barton and Turelli developed recursions to describe the evolution of multilocus systems under arbitrary forms of selection. This article generalizes their approach to allow for arbitrary modes of inheritance, including diploidy, polyploidy, sex linkage, cytoplasmic inheritance, and genomic imprinting. The framework is also extended to allow for other deterministic evolutionary forces, including migration and mutation. Exact recursions that fully describe the state of the population are presented; these are implemented in a computer algebra package (available on the Web at http://helios.bto.ed.ac.uk/evolgen). Despite the generality of our framework, it can describe evolutionary dynamics exactly by just two equations. These recursions can be further simplified using a "quasi-linkage equilibrium" (QLE) approximation. We illustrate the methods by finding the effect of natural selection, sexual selection, mutation, and migration on the genetic composition of a population.     

Furin cleavage is not a requirement for Drosophila Notch function

Notch (N) is a large transmembrane protein that acts as a receptor in an evolutionarily conserved intercellular signalling pathway. Because of this conservation, it has been assumed that biochemical events mediating N function are identical in all species. For instance, intracellular maturation by furin protease and subunit assembly leading to the formation of a heterodimeric cell surface N receptor are thought to be central to its function in both mammals and flies. However, in Drosophila the majority of N appears to be full-length. It has not been determined whether this full-length N protein is on the cell surface. We describe experiments which indicate that unlike mammalian N, the majority of Drosophila N on the cell surface is full-length and that in Drosophila, in vivo, furin cleavage is not required for biological activity. We further show that the behaviour of fly and mouse N can be interchanged simply by swapping the regions in which the mammalian furin-like cleavage site is located.     

Gradiflow as a prefractionation tool for two-dimensional electrophoresis

The Gradiflow trade mark, a preparative electrophoresis instrument capable of separating proteins on the basis of their size or charge, was used to separate whole cell lysates, prepared from bakers yeast (Saccharomyces cerevisiae) and Chinese snow pea seeds (Pisum sativum macrocarpon), into protein fractions of different pH regions. Both broad and narrow range (with a difference of approximately 1 pH unit) pH fractions were obtained. Analysis of the protein fractions by isoelectric focusing gels and two-dimensional (2-D) polyacrylamide gel electrophoresis indicated minimal overlap between the pH fractions. Further, when the prefractionated acidic samples were analyzed on pH 4-7 immobilized pH gradient 2-D gels, improved resolution of the proteins within the chosen pH region was achieved compared to the unfractionated samples. This study demonstrates that the Gradiflow could be used as a preparative electrophoresis tool for the isolation of proteins into distinct pH fractions.     

Manipulation of dendritic cells as an approach to improved outcomes in transplantation

At the time of organ transplantation, a variety of non-parenchymal cells are transplanted simultaneously with the allograft. Recognition of the importance of these cells as potential immunostimulatory cells lead to the concept of 'passenger leukocytes' as the principal instigators of rejection. Passenger leukocytes include interstitial dendritic cells (DCs) and blood-derived monocytes/macrophages. As investigators have discovered the significance of DCs in influencing graft outcome, so have they begun to determine the best ways to influence DCs themselves. This review discusses the role of DCs in transplantation and then focuses on three different approaches for manipulating DCs to improve allograft survival: (1) targeting of chemokines involved in DC migration, (2) pharmacological arrest of DC maturation, and (3) genetic engineering of DCs.     

Modeling diverse range of potassium channels with Brownian dynamics

Using the experimentally determined KcsA structure as a template, we propose a plausible explanation for the diversity of potassium channels seen in nature. A simplified model of KcsA is constructed from its atomic resolution structure by smoothing out the protein-water boundary and representing the atoms forming the channel protein as a homogeneous, low dielectric medium. The properties of the simplified and atomic-detail models, deduced from electrostatic calculations and Brownian dynamics simulations, are shown to be qualitatively similar. We then study how the current flowing across the simplified model channel changes as the shape of the intrapore region is modified. This is achieved by increasing the radius of the intracellular pore systematically from 1.5 to 5 A while leaving the dimensions of the selectivity filter and inner chamber unaltered. The strengths of the dipoles located near the entrances of the channel, the carbonyl groups lining the selectivity filter, and the helix macrodipoles are kept constant. The channel conductance increases steadily as the radius of the intracellular pore is increased. The rate-limiting step for both the outward and inward current is the time it takes for an ion to cross the residual energy barrier located in the intrapore region. The current-voltage relationship obtained with symmetrical solutions is linear when the applied potential is less than approximately 100 mV but deviates slightly from Ohm's law at higher applied potentials. The nonlinearity in the current-voltage curve becomes less pronounced as the radius of the intracellular pore is increased. When the strengths of the dipoles near the intracellular entrance are reduced, the channel shows a pronounced inward rectification. Finally, the conductance exhibits the saturation property observed experimentally. We discuss the implications of these findings on the transport of ions across the potassium channels and membrane channels in general.     

Redox regulation of neuronal migration in a Down Syndrome model

Down Syndrome (DS), one of the major genetic causes of mental retardation, is characterized by disrupted corticogenesis produced, in part, by an abnormal layering of neurons in cortical laminas II and III. Because defects in the normal migration of neurons during corticogenesis can result in delayed cortical radial expansion and abnormalities in cortical layering, we have examined neuronal migration in murine trisomy 16 (Ts16), a mouse model for DS. Using an in vitro assay for chemotaxis, our data demonstrate that the number of acutely dissociated Ts16 cortical neurons migrating in response to glutamate or N-methyl-D-aspartate (NMDA), known chemotactic factors, was decreased compared to normal littermates, suggesting a defect in NMDA receptor- (NMDAR-) mediated events. Ts16 neurons did not lack NMDAR since expression of mRNA and protein for NMDAR subunits was observed in Ts16 cells. However, the number of cells that generated an observable current in response to NMDA was decreased compared to normal littermates. Similar to DS, Ts16 CNS demonstrated an inherent oxidative stress likely caused by the triplication of genes such as SOD1. To determine if the abnormal redox state was a factor in the failure of NMDAR-mediated migration in Ts16, we treated Ts16 neurons with either n-acetyl cysteine (NAC) or dithiothrietol (DTT), known antioxidants. The reduction in NMDAR-mediated migration observed in Ts16 neurons was returned to normal littermate values by NAC or DTT. Our data indicate that oxidative stress may play a key role in the abnormal glutamate-mediated responses during cortical development in the Ts16 mouse and may have an impact on neuronal migration at critical stages.