Molecular dynamics of ion transport through the open conformation of a bacterial voltage-gated sodium channel

The crystal structure of the open conformation of a bacterial voltage-gated sodium channel pore from Magnetococcus sp. (NaVMs) has provided the basis for a molecular dynamics study defining the channel’s full ion translocation pathway and conductance process, selectivity, electrophysiological characteristics, and ion-binding sites. Microsecond molecular dynamics simulations permitted a complete time-course characterization of the protein in a membrane system, capturing the plethora of conductance events and revealing a complex mixture of single and multi-ion phenomena with decoupled rapid bidirectional water transport. The simulations suggest specific localization sites for the sodium ions, which correspond with experimentally determined electron density found in the selectivity filter of the crystal structure. These studies have also allowed us to identify the ion conductance mechanism and its relation to water movement for the NavMs channel pore and to make realistic predictions of its conductance properties. The calculated single-channel conductance and selectivity ratio correspond closely with the electrophysiology measurements of the NavMs channel expressed in HEK 293 cells. The ion translocation process seen in this voltage-gated sodium channel is clearly different from that exhibited by members of the closely related family of voltage-gated potassium channels and also differs considerably from existing proposals for the conductance process in sodium channels. These studies simulate sodium channel conductance based on an experimentally determined structure of a sodium channel pore that has a completely open transmembrane pathway and activation gate

Conformational Preferences of a 14-Residue Fibrillogenic Peptide from Acetylcholinesterase†

A 14-residue fragment from near the C-terminus of the enzyme acetylcholinesterase (AChE) is believed to have a neurotoxic/neurotrophic effect acting via an unknown pathway. While the peptide is α-helical in the full-length enzyme, the structure and association mechanism of the fragment are unknown. Using multiple molecular dynamics simulations, starting from a tetrameric complex of the association domain of AChE and systematicall disassembled subsets that include the peptide fragment, we show that the fragment is incapable of retaining its helicity in solution. Extensive replica exchange Monte Carlo folding and unfolding simulations in implicit solvent with capped and uncappted termini failed to converge to any consistent cluster of structures, suggesting that the fragment remains largely unstructured in solution under the conditions considered. Furthermore, extended molecular dynamics simulations of two steric zipper models show that the peptide is likely to form a zipper with antiparallel sheets and that peptides with mutations known to prevent fibril formation likely do so by interfering with this packing. The results demonstrate how the local environment of a peptide can stabilize a particular conformation

Tracking magnetic bright point motions through the solar atmosphere

High-cadence, multiwavelength observations and simulations are employed for the analysis of solar photospheric magnetic bright points (MBPs) in the quiet Sun. The observations were obtained with the Rapid Oscillations in the Solar Atmosphere (ROSA) imager and the Interferometric Bidimensional Spectrometer at the Dunn Solar Telescope. Our analysis reveals that photospheric MBPs have an average transverse velocity of approximately 1 km s−1, whereas their chromospheric counterparts have a slightly higher average velocity of 1.4 km s−1. Additionally, chromospheric MBPs were found to be around 63 per cent larger than the equivalent photospheric MBPs. These velocity values were compared with the output of numerical simulations generated using the MURAM code. The simulated results were similar, but slightly elevated, when compared to the observed data. An average velocity of 1.3 km s−1 was found in the simulated G-band images and an average of 1.8 km s−1 seen in the velocity domain at a height of 500 km above the continuum formation layer. Delays in the change of velocities were also analysed. Average delays of ∼4 s between layers of the simulated data set were established and values of ∼29 s observed between G-band and Ca II K ROSA observations. The delays in the simulations are likely to be the result of oblique granular shock waves, whereas those found in the observations are possibly the result of a semi-rigid flux tube

Fundamental parameters, integrated RGB mass loss and dust production in the Galactic globular cluster 47 Tucanae

Fundamental parameters and time-evolution of mass loss are investigated for post-main-sequence stars in the Galactic globular cluster 47 Tucanae (NGC 104). This is accomplished by fitting spectral energy distributions (SEDs) to existing optical and infrared photometry and spectroscopy, to produce a true Hertzsprung--Russell diagram. We confirm the cluster's distance as 4611 (+213, -200) pc and age as 12 +/- 1 Gyr. Horizontal branch models appear to confirm that no more RGB mass loss occurs in 47 Tuc than in the more-metal-poor omega Centauri, though difficulties arise due to inconsistencies between the models. Using our SEDs, we identify those stars which exhibit infrared excess, finding excess only among the brightest giants: dusty mass loss begins at a luminosity of ~ 1000 Lsun, becoming ubiquitous above 2000 Lsun. Recent claims of dust production around lower-luminosity giants cannot be reproduced, despite using the same archival Spitzer imagery.Comment: 22 pages, 17 figures, accepted ApJ

Lipid membranes for membrane proteins

Andreas Kukol, ‘Lipid membranes for membrane proteins in Molecular Modeling of Proteins (Clifton: Humana Press/Sringer, 2015), ISBN: 978-1-4939-1464-7, e-BOOK ISBN: 978-1-4939-1465-4Peer reviewe

Temperature‐Dependent Re‐alignment of the Short Multifunctional Peptide BP100 in Membranes Revealed by Solid‐State NMR Spectroscopy and Molecular Dynamics Simulations

BP100 is a cationic undecamer peptide with antimicrobial and cell-penetrating activities. The orientation of this amphiphilic α-helix in lipid bilayers was examined under numerous conditions using solid-state 19F, 15N and 2H NMR. At high temperatures in saturated phosphatidylcholine lipids, BP100 lies flat on the membrane surface, as expected. Upon lowering the temperature towards the lipid phase transition, the helix is found to flip into an upright transmembrane orientation. In thin bilayers, this inserted state was stable at low peptide concentration, but thicker membranes required higher peptide concentrations. In the presence of lysolipids, the inserted state prevailed even at high temperature. Molecular dynamics simulations suggest that BP100 monomer insertion can be stabilized by snorkeling lysine side chains. These results demonstrate that even a very short helix like BP100 can span (and thereby penetrate through) a cellular membrane under suitable conditions

Free energy barrier for melittin reorientation from a membrane-bound state to a transmembrane state

An important step in a phospholipid membrane pore formation by melittin antimicrobial peptide is a reorientation of the peptide from a surface into a transmembrane conformation. In this work we perform umbrella sampling simulations to calculate the potential of mean force (PMF) for the reorientation of melittin from a surface-bound state to a transmembrane state and provide a molecular level insight into understanding peptide and lipid properties that influence the existence of the free energy barrier. The PMFs were calculated for a peptide to lipid (P/L) ratio of 1/128 and 4/128. We observe that the free energy barrier is reduced when the P/L ratio increased. In addition, we study the cooperative effect; specifically we investigate if the barrier is smaller for a second melittin reorientation, given that another neighboring melittin was already in the transmembrane state. We observe that indeed the barrier of the PMF curve is reduced in this case, thus confirming the presence of a cooperative effect

Atmospheric Heating and Wind Acceleration: Results for Cool Evolved Stars based on Proposed Processes

A chromosphere is a universal attribute of stars of spectral type later than ~F5. Evolved (K and M) giants and supergiants (including the zeta Aurigae binaries) show extended and highly turbulent chromospheres, which develop into slow massive winds. The associated continuous mass loss has a significant impact on stellar evolution, and thence on the chemical evolution of galaxies. Yet despite the fundamental importance of those winds in astrophysics, the question of their origin(s) remains unsolved. What sources heat a chromosphere? What is the role of the chromosphere in the formation of stellar winds? This chapter provides a review of the observational requirements and theoretical approaches for modeling chromospheric heating and the acceleration of winds in single cool, evolved stars and in eclipsing binary stars, including physical models that have recently been proposed. It describes the successes that have been achieved so far by invoking acoustic and MHD waves to provide a physical description of plasma heating and wind acceleration, and discusses the challenges that still remain.Comment: 46 pages, 9 figures, 1 table; modified and unedited manuscript; accepted version to appear in: Giants of Eclipse, eds. E. Griffin and T. Ake (Berlin: Springer

Wave Damping Observed in Upwardly Propagating Sausage-mode Oscillations contained within a Magnetic Pore

We present observational evidence of compressible MHD wave modes propagating from the solar photosphere through to the base of the transition region in a solar magnetic pore. High cadence images were obtained simultaneously across four wavelength bands using the Dunn Solar Telescope. Employing Fourier and wavelet techniques, sausage-mode oscillations displaying significant power were detected in both intensity and area fluctuations. The intensity and area fluctuations exhibit a range of periods from 181 to 412 s, with an average period ~290 s, consistent with the global p-mode spectrum. Intensity and area oscillations present in adjacent bandpasses were found to be out of phase with one another, displaying phase angles of 6fdg12, 5fdg82, and 15fdg97 between the 4170 Å continuum–G-band, G-band–Na i D1, and Na i D1–Ca ii K heights, respectively, reiterating the presence of upwardly propagating sausage-mode waves. A phase relationship of ~0° between same-bandpass emission and area perturbations of the pore best categorizes the waves as belonging to the "slow" regime of a dispersion diagram. Theoretical calculations reveal that the waves are surface modes, with initial photospheric energies in excess of 35,000 W m−2. The wave energetics indicate a substantial decrease in energy with atmospheric height, confirming that magnetic pores are able to transport waves that exhibit appreciable energy damping, which may release considerable energy into the local chromospheric plasma

Multiwavelength studies of MHD waves in the solar chromosphere: An overview of recent results

The chromosphere is a thin layer of the solar atmosphere that bridges the relatively cool photosphere and the intensely heated transition region and corona. Compressible and incompressible waves propagating through the chromosphere can supply significant amounts of energy to the interface region and corona. In recent years an abundance of high-resolution observations from state-of-the-art facilities have provided new and exciting ways of disentangling the characteristics of oscillatory phenomena propagating through the dynamic chromosphere. Coupled with rapid advancements in magnetohydrodynamic wave theory, we are now in an ideal position to thoroughly investigate the role waves play in supplying energy to sustain chromospheric and coronal heating. Here, we review the recent progress made in characterising, categorising and interpreting oscillations manifesting in the solar chromosphere, with an impetus placed on their intrinsic energetics.Comment: 48 pages, 25 figures, accepted into Space Science Review