Biomedical subjects

Eric Jakobsson

Publications and source records attributed to Eric Jakobsson.

15 recordsLinked to original sources

PubMed2006

The cPLA2 C2alpha domain in solution: structure and dynamics of its Ca2+-activated and cation-free states.

Cytosolic phospholipase A2 is involved in several signal transduction pathways where it catalyses release of arachidonic acid from intracellular lipid membranes. Its membrane insertion is facilitated by its independently folding C2alpha domain, which is activated by the binding of two intracellular Ca2+ ions. However, the details of its membrane-insertion mechanism, including its Ca2+-activation mechanism, are not understood. There are several unresolved issues, including the following. There are two experimentally resolved structures of the Ca2+-activated state of its isolated C2alpha domain, one determined using x-ray crystallography and the other determined using NMR spectroscopy, which differ from each other significantly in the spatial region that inserts into the membrane. This by itself adds to ambiguities associated with investigations targeting its mechanism of membrane insertion. Furthermore, there is no experimentally determined structure of its cation-free state, which hinders investigations associated with its cation-activation mechanism. In this work, we generate several unrestrained molecular dynamics trajectories of its isolated C2alpha domain in solution (equivalent to approximately 60 ns) and investigate these issues. Our main results are as follows: a), the Ca2+ coordination scheme of the domain is consistent with the x-ray structure and with previous mutagenesis studies; b), the helical segment of the Ca2+-binding loop, CBL-I, undergoes nanosecond timescale flexing (but not an unwinding), as can be inferred from physiological temperature NMR data and in contrast to low temperature x-ray data; and c), removal of the two activating Ca2+ ions from their binding pockets does not alter the backbone structure of the domain, a result consistent with electron paramagnetic resonance data.

Binding Sites↗

PubMed2006

Lateral organization in lipid-cholesterol mixed bilayers.

Interactions between lipid and cholesterol molecules in membranes play an important role in the structural and functional properties of cell membranes. Although structural properties of lipid-cholesterol mixtures have been extensively studied, an understanding of the role of cholesterol in the lateral organization of bilayers has been elusive. In this article, we propose a simple yet powerful model, based on self-consistent mean-field theory and molecular dynamics simulations, for lipid bilayers containing cholesterol. Properties predicted by our model are shown to be in excellent agreement with experimental data. Our model predicts that cholesterol induces structural changes in the bilayer through the formation of regions of ordered lipids surrounding each cholesterol molecule. We find that the "smooth" and "rough" sides of cholesterol play crucial roles in formation and distribution of the ordered regions. Our model is predictive in that input parameters are obtained from independent atomistic molecular dynamics simulations. The model and method are general enough to describe other heterogeneous lipid bilayers, including lipid rafts.

Cholesterol↗

PubMed2006

Cholesterol surrogates: a comparison of cholesterol and 16:0 ceramide in POPC bilayers.

Experimental evidence indicates that, under some circumstances, "surrogate" molecules may play the same role as cholesterol in ordering membrane lipids. The simplest molecule in this class is Ceramide. In this article, we describe atomic-level molecular dynamics simulations designed to shed light on this phenomenon. We run simulations of hydrated phosphoryl-oleoyl phosphatidylcholine (POPC) bilayers containing cholesterol, and containing ceramide, in concentrations ranging from 5% to 33%. We also perform a simulation of a pure POPC bilayer to verify the simulation force fields against experimental structural data for POPC. Our simulation data are in good agreement with experimental data for the partial molecular volumes, areas, form factors, and order parameters. These simulations suggest that ceramide and cholesterol have a very similar effect on the POPC bilayer, although ceramide is less effective in inducing order in the bilayer compared with cholesterol at the same concentrations.

Ceramides↗

PubMed2006

Calcium channel auxiliary subunits.

Many channels and carriers associate with auxiliary subunits which modify their activities and facilitate biogenesis. Advances in genome sequencing as well as biochemical, molecular genetic, and physiological experimentation have allowed for the discovery of many transport auxiliary subunits. Recent interests in the pharmacology of the calcium auxiliary subunits prompted a large amount of effort in deciphering their specific role in the conductance of calcium ions. In this review, we evaluate the functions of the 'extra' subunits of the voltage-gated calcium channels in animals as an example of auxiliary subunits of transporters in general. We discuss the functional data available for each of these subunits, present phylogenetic analyses, and discuss their potential evolutionary origins. Our analyses also reveal novel homologues of these subunits which might be of interest to the community.

Amino Acid Sequence↗

PubMed2005

The influence of amino acid protonation states on molecular dynamics simulations of the bacterial porin OmpF.

Several groups, including our own, have found molecular dynamics (MD) calculations to result in the size of the pore of an outer membrane bacterial porin, OmpF, to be reduced relative to its size in the x-ray crystal structure. At the narrowest portion of its pore, loop L3 was found to move toward the opposite face of the pore, resulting in decreasing the cross-section area by a factor of approximately 2. In an earlier work, we computed the protonation states of titratable residues for this system and obtained values different from those that had been used in previous MD simulations. Here, we show that MD simulations carried out with these recently computed protonation states accurately reproduce the cross-sectional area profile of the channel lumen in agreement with the x-ray structure. Our calculations include the investigation of the effect of assigning different protonation states to the one residue, D(127), whose protonation state could not be modeled in our earlier calculation. We found that both assumptions of charge states for D(127) reproduced the lumen size profile of the x-ray structure. We also found that the charged state of D(127) had a higher degree of hydration and it induced greater mobility of polar side chains in its vicinity, indicating that the apparent polarizability of the D(127) microenvironment is a function of the D(127) protonation state.

Algorithms↗

PubMed2005

Digital biology: an emerging and promising discipline.

This article examines the role of computation and quantitative methods in modern biomedical research to identify emerging scientific, technical, policy and organizational trends. It identifies common concerns and practices in the emerging community of computationally-oriented bio-scientists by reviewing a national symposium, Digital Biology: the Emerging Paradigm, held at the National Institutes of Health in Bethesda, Maryland, November 6th and 7th 2003. This meeting showed how biomedical computing promises scientific breakthroughs that will yield significant health benefits. Three key areas that define the emerging discipline of digital biology are: scientific data integration, multi-scale modeling and networked science. Each area faces unique technical challenges and information policy issues that must be addressed as the field matures. Here we summarize the emergent challenges and offer suggestions to academia, industry and government on how best to expand the role of computation in their scientific activities.

Computational Biology↗

PubMed2004

Ca2+ activation of the cPLA2 C2 domain: ordered binding of two Ca2+ ions with positive cooperativity.

During Ca(2+) activation, the Ca(2+)-binding sites of C2 domains typically bind multiple Ca(2+) ions in close proximity. These binding events exhibit positive cooperativity, despite the strong charge repulsion between the adjacent divalent cations. Using both experimental and computational approaches, the present study probes the detailed mechanisms of Ca(2+) activation and positive cooperativity for the C2 domain of cytosolic phospholipase A(2), which binds two Ca(2+) ions in sites I and II, separated by only 4.1 A. First, each of the five coordinating side chains in the Ca(2+)-binding cleft is individually mutated and the effect on Ca(2+)-binding affinity and cooperativity is measured. The results identify Asp 43 as the major contributor to Ca(2+) affinity, while the two coordinating side chains that provide bridging coordination to both Ca(2+) ions, Asp 43 and Asp 40, are observed to make the largest contributions to positive cooperativity. Electrostatic calculations reveal that Asp 43 possesses the highest pseudo-pK(a) of the coordinating acidic residues, as well as the highest general cation affinity, due to its relatively buried location within 3.5 A of seven protein oxygens with full or partial negative charges. These calculations therefore explain the greater importance of Asp 43 in defining the Ca(2+) affinity. Overall, the experimental and computational results support an activation model in which the first Ca(2+) ion binds usually to site I, thereby preordering both bridging side chains Asp 40 and 43, and partially or fully deprotonating the three coordinating Asp residues. This initial binding event prepares the conformation and protonation state of the remaining site for Ca(2+) binding, enabling the second Ca(2+) ion to bind with higher affinity than the first as required for positive cooperativity.

Amino Acid Substitution↗

PubMed2004

Identification of the prokaryotic ligand-gated ion channels and their implications for the mechanisms and origins of animal Cys-loop ion channels.

BACKGROUND: Acetylcholine receptor type ligand-gated ion channels (ART-LGIC; also known as Cys-loop receptors) are a superfamily of proteins that include the receptors for major neurotransmitters such as acetylcholine, serotonin, glycine, GABA, glutamate and histamine, and for Zn2+ ions. They play a central role in fast synaptic signaling in animal nervous systems and so far have not been found outside of the Metazoa. RESULTS: Using sensitive sequence-profile searches we have identified homologs of ART-LGICs in several bacteria and a single archaeal genus, Methanosarcina. The homology between the animal receptors and the prokaryotic homologs spans the entire length of the former, including both the ligand-binding and channel-forming transmembrane domains. A sequence-structure analysis using the structure of Lymnaea stagnalis acetylcholine-binding protein and the newly detected prokaryotic versions indicates the presence of at least one aromatic residue in the ligand-binding boxes of almost all representatives of the superfamily. Investigation of the domain architectures of the bacterial forms shows that they may often show fusions with other small-molecule-binding domains, such as the periplasmic binding protein superfamily I (PBP-I), Cache and MCP-N domains. Some of the bacterial forms also occur in predicted operons with the genes of the PBP-II superfamily and the Cache domains. Analysis of phyletic patterns suggests that the ART-LGICs are currently absent in all other eukaryotic lineages except animals. Moreover, phylogenetic analysis and conserved sequence motifs also suggest that a subset of the bacterial forms is closer to the metazoan forms. CONCLUSIONS: From the information from the bacterial forms we infer that cation-pi or hydrophobic interactions with the ligand are likely to be a pervasive feature of the entire superfamily, even though the individual residues involved in the process may vary. The conservation pattern in the channel-forming transmembrane domains also suggests similar channel-gating mechanisms in the prokaryotic versions. From the distribution of charged residues in the prokaryotic M2 transmembrane segments, we expect that there will be examples of both cation and anion selectivity within the prokaryotic members. Contextual connections suggest that the prokaryotic forms may function as chemotactic receptors for low molecular weight solutes. The phyletic patterns and phylogenetic relationships suggest the possibility that the metazoan receptors emerged through an early lateral transfer from a prokaryotic source, before the divergence of extant metazoan lineages.

Animals↗

PubMed2004

Simulation of the early stages of nano-domain formation in mixed bilayers of sphingomyelin, cholesterol, and dioleylphosphatidylcholine.

It is known from experimental studies that lipid bilayers composed of unsaturated phospholipids, sphingomyelin, and cholesterol contain microdomains rich in sphingomyelin and cholesterol. These domains are similar to "rafts" isolated from cell membranes, although the latter are much smaller in lateral size. Such domain formation can be a result of very specific and subtle lipid-lipid interactions. To identify and study these interactions, we have performed two molecular dynamics simulations, of 200-ns duration, of dioleylphosphatidylcholine (DOPC), sphingomyelin (SM), and cholesterol (Chol) systems, a 1:1:1 mixture of DOPC/SM/Chol, and a 1:1 mixture of DOPC/SM. The simulations show initial stages of the onset of spontaneous phase-separated domains in the systems. On the simulation timescale cholesterol favors a position at the interface between the ordered SM region and the disordered DOPC region in the ternary system and accelerates the process of domain formation. We find that the smooth alpha-face of Chol preferentially packs next to SM molecules. Based on a comparative analysis of interaction energies, we find that Chol molecules do not show a preference for SM or DOPC. We conclude that Chol molecules assist in the process of domain formation and the process is driven by entropic factors rather than differences in interaction energies.

Cholesterol↗

PubMed2004

Ionization states of residues in OmpF and mutants: effects of dielectric constant and interactions between residues.

To understand ion permeation, one must assign correct ionization states to titratable amino acid residues in protein channels. We report on the effects of physical and methodological assumptions in calculating the protonation states at neutral bulk pH of titratable residues lining the lumen of the native Escherichia coli OmpF channel, and five mutants. We systematically considered a wide range of assumed protein dielectric constants and all plausible combinations of protonation states for electrostatically interacting side chains, and three different levels of accounting for solute shielding: 1), full nonlinear Poisson-Boltzmann; 2), linearized Poisson-Boltzmann; and 3), neglect of solute shielding. For this system we found it acceptable to neglect solute shielding, a result we postulate to be generalizable to narrow lumens of other protein channels. For the large majority of residues, the protonation state at neutral bulk pH was found to be independent of the assumed dielectric constant of the protein, and unambiguously determined by the calculation; for native OmpF only Asp-127 has a protonation state that is sensitive to the assumed protein dielectric constant. Our results are significant for understanding two published experimental observations: the structure of the narrow part of the channel, and the ionic selectivity of OmpF mutants.

Amino Acid Substitution↗

PubMed2004

Sphingomyelin-cholesterol domains in phospholipid membranes: atomistic simulation.

We have carried out an atomic-level molecular dynamics simulation of a system of nanoscopic size containing a domain of 18:0 sphingomyelin and cholesterol embedded in a fully hydrated dioleylposphatidylcholine (DOPC) bilayer. To analyze the interaction between the domain and the surrounding phospholipid, we calculate order parameters and area per molecule as a function of molecule type and proximity to the domain. We propose an algorithm based on Voronoi tessellation for the calculation of the area per molecule of various constituents in this ternary mixture. The calculated areas per sphingomyelin and cholesterol are in agreement with previous simulations. The simulation reveals that the presence of the liquid-ordered domain changes the packing properties of DOPC bilayer at a distance as large as approximately 8 nm. We calculate electron density profiles and also calculate the difference in the thickness between the domain and the surrounding DOPC bilayer. The calculated difference in thickness is consistent with data obtained in atomic force microscopy experiments.

Algorithms↗

PubMed2003

Sequence-function analysis of the K+-selective family of ion channels using a comprehensive alignment and the KcsA channel structure.

Sequence-function analysis of K(+)-selective channels was carried out in the context of the 3.2 A crystal structure of a K(+) channel (KcsA) from Streptomyces lividans (Doyle et al., 1998). The first step was the construction of an alignment of a comprehensive set of K(+)-selective channel sequences forming the putative permeation path. This pathway consists of two transmembrane segments plus an extracellular linker. Included in the alignment are channels from the eight major classes of K(+)-selective channels from a wide variety of species, displaying varied rectification, gating, and activation properties. Segments of the alignment were assigned to structural motifs based on the KcsA structure. The alignment's accuracy was verified by two observations on these motifs: 1), the most variability is shown in the turret region, which functionally is strongly implicated in susceptibility to toxin binding; and 2), the selectivity filter and pore helix are the most highly conserved regions. This alignment combined with the KcsA structure was used to assess whether clusters of contiguous residues linked by hydrophobic or electrostatic interactions in KcsA are conserved in the K(+)-selective channel family. Analysis of sequence conservation patterns in the alignment suggests that a cluster of conserved residues is critical for determining the degree of K(+) selectivity. The alignment also supports the near-universality of the "glycine hinge" mechanism at the center of the inner helix for opening K channels. This mechanism has been suggested by the recent crystallization of a K channel in the open state. Further, the alignment reveals a second highly conserved glycine near the extracellular end of the inner helix, which may be important in minimizing deformation of the extracellular vestibule as the channel opens. These and other sequence-function relationships found in this analysis suggest that much of the permeation path architecture in KcsA is present in most K(+)-selective channels. Because of this finding, the alignment provides a robust starting point for homology modeling of the permeation paths of other K(+)-selective channel classes and elucidation of sequence-function relationships therein. To assay these applications, a homology model of the Shaker A channel permeation path was constructed using the alignment and KcsA as the template, and its structure evaluated in light of established structural criteria.

Bacterial Proteins↗

PubMed2003

Structure of sphingomyelin bilayers: a simulation study.

We have carried out a molecular dynamics simulation of a hydrated 18:0 sphingomyelin lipid bilayer. The bilayer contained 1600 sphingomyelin (SM) molecules, and 50,592 water molecules. After construction and initial equilibration, the simulation was run for 3.8 ns at a constant temperature of 50 degrees C and a constant pressure of 1 atm. We present properties of the bilayer calculated from the simulation, and compare with experimental data and with properties of dipalmitoyl phosphatidylcholine (DPPC) bilayers. The SM bilayers are significantly more ordered and compact than DPPC bilayers at the same temperature. SM bilayers also exhibit significant intramolecular hydrogen bonding between phosphate ester oxygen and hydroxyl hydrogen atoms. This results in a decreased hydration in the polar region of the SM bilayer compared with DPPC. Since our simulation system is very large we have calculated the power spectrum of bilayer undulation and peristaltic modes, and we compare these data with similar calculations for DPPC bilayers. We find that the SM bilayer has significantly larger bending modulus and area compressibility compared to DPPC.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2003

Phylogenetic analyses of potassium channel auxiliary subunits.

Results of recent genome-sequencing projects together with advances in biochemical, molecular genetic and physiological experimentation have allowed discovery of many transport auxiliary subunits. These subunits facilitate the proper movement of substrates across cell membranes. Mutations of any of these subunits can cause catastrophic effects to the transport mechanism and cause certain genetic diseases. Auxiliary subunits of ion channels are of particular interest because of their potential to diversify the transport properties of the principal subunits. Furthermore, ion channel auxiliary subunits may function in the capacity of enhancing surface expression, allowing gating, and providing chaperone-like activities. As a result of their evolutionary histories, these proteins can be grouped exclusively by phylogenetic techniques. Many of these families are found to be restricted to a single kingdom of life while others seem to be ubiquitous. Here we report the results of systematic analyses of three families of ion channel auxiliary subunits. Some subunits contain one or more transmembrane segments while others exist only in the cytoplasm. We have also observed potential horizontal transfer across kingdoms with these auxiliary subunits. In this report, we present tabulated results of homology searches, partial multiple alignments, secondary structure analyses, and phylogenetic trees for these families.

Animals↗

PubMed2002

Cholesterol-induced modifications in lipid bilayers: a simulation study.

We present analysis of new configurational bias Monte Carlo and molecular dynamics simulation data for bilayers of dipalmitoyl phosphatidyl choline and cholesterol for dipalmitoyl phosphatidyl choline:cholesterol ratios of 24:1, 47:3, 11.5:1, 8:1, 7:1, 4:1, 3:1, 2:1, and 1:1, using long molecular dynamics runs and interspersed configurational bias Monte Carlo to ensure equilibration and enhance sampling. In all cases with cholesterol concentrations above 12.5% the area per molecule of the heterogeneous membrane varied linearly with cholesterol fraction. By extrapolation to pure cholesterol, we find the cross-sectional area of cholesterol in these mixtures is approximately 22.3 A(2). From the slope of the area/molecule relationship, we also find that the phospholipid in these mixtures is in a liquid ordered state with an average cross-sectional area per lipid of 50.7 A(2), slightly above the molecular area of a pure phospholipid gel. For lower concentrations of cholesterol, the molecular area rises above the straight line, indicating the "melting" of at least some of the phospholipid into a fluid state. Analysis of the lateral distribution of cholesterol molecules in the leaflets reveals peaks in radial distributions of cholesterols at multiples of approximately 5 A. These peaks grow in size as the simulation progresses, suggesting a tendency for small subunits of one lipid plus one cholesterol, hydrogen bonded together, to act as one composite particle, and perhaps to aggregate with other composites. Our results are consistent with experimentally observed effects of cholesterol, including the condensation effect of cholesterol in phospholipid monolayers and the tendency of cholesterol-rich domains to form in cholesterol-lipid bilayers. We are continuing to analyze this tendency on longer timescales and for larger bilayer patches.

Biophysical Phenomena↗