Atomistic Characterization of Gramicidin Channel Formation.
Ontology highlight
ABSTRACT:
SUBMITTER: Sun D
PROVIDER: S-EPMC7808174 | biostudies-literature | 2021 Jan
REPOSITORIES: biostudies-literature
Publications
Atomistic Characterization of Gramicidin Channel Formation.
Sun Delin D He Stewart S Bennett W F Drew WFD Bilodeau Camille L CL Andersen Olaf S OS Lightstone Felice C FC Ingólfsson Helgi I HI
Journal of chemical theory and computation 20201230 1
We investigated gramicidin A (gA) subunit dimerization in lipid bilayers using microsecond-long replica-exchange umbrella sampling simulations, millisecond-long unbiased molecular dynamics simulations, and machine learning. Our simulations led to a dimer structure that is indistinguishable from the experimentally determined gA channel structures, with the two gA subunits joined by six hydrogen bonds (6HB). The simulations also uncovered two additional dimer structures, with different gA-gA stack ...[more]
Similar Datasets
Project description:Integral membrane protein function can be modulated by the host bilayer. Because biological membranes are diverse and nonuniform, we explore the consequences of lipid diversity using gramicidin A channels embedded in phosphatidylcholine (PC) bilayers composed of equimolar mixtures of di-oleoyl-PC and di-erucoyl-PC (dC18:1+dC22:1, respectively), di-palmitoleoyl-PC and di-nervonoyl-PC (dC16:1+dC24:1, respectively), and di-eicosenoyl-PC (pure dC20:1), all of which have the same average bilayer chain length. Single-channel lifetime experiments, molecular dynamics simulations, and a simple lipid compression model are used in tandem to gain insight into lipid redistribution around the channel, which partially alleviates the bilayer deformation energy associated with channel formation. The average single-channel lifetimes in the two-component bilayers (95 ± 10 ms for dC18:1+dC22:1 and 195 ± 20 ms for dC16:1+dC24:1) were increased relative to the single-component dC20:1 control bilayer (65 ± 10 ms), implying lipid redistribution. Using a theoretical treatment of thickness-dependent changes in channel lifetimes, the effective local enrichment of lipids around the channel was estimated to be 58 ± 4% dC18:1 and 66 ± 2% dC16:1 in the dC18:1+dC22:1 and dC16:1+dC24:1 bilayers, respectively. 3.5-μs molecular dynamics simulations show 66 ± 2% dC16:1 in the first lipid shell around the channel in the dC16:1+dC24:1 bilayer, but no significant redistribution (50 ± 4% dC18:1) in the dC18:1+dC22:1 bilayer; these simulated values are within the 95% confidence intervals of the experimental averages. The strong preference for the better matching lipid (dC16:1) near the channel in the dC16:1+dC24:1 mixture and lesser redistribution in the dC18:1+dC22:1 mixture can be explained by the energetic cost associated with compressing the lipids to match the channel's hydrophobic length.
| S-EPMC5375002 | biostudies-literature
Project description:To change conformation, a protein must deform the surrounding bilayer. In this work, a three-dimensional continuum elastic model for gramicidin A in a lipid bilayer is shown to describe the sensitivity to thickness, curvature stress, and the mechanical properties of the lipid bilayer. A method is demonstrated to extract the gramicidin-lipid boundary condition from all-atom simulations that can be used in the three-dimensional continuum model. The boundary condition affects the deformation dramatically, potentially much more than typical variations in the material stiffness do as lipid composition is changed. Moreover, it directly controls the sensitivity to curvature stress. The curvature stress and hydrophobic surfaces of the all-atom and continuum models are found to be in excellent agreement. The continuum model is applied to estimate the enrichment of hydrophobically matched lipids near the channel in a mixture, and the results agree with single-channel experiments and extended molecular dynamics simulations from the companion article by Beaven et al. in this issue of Biophysical Journal.
| S-EPMC5376118 | biostudies-literature
Project description:The common occurrence of Trp residues at the aqueous-lipid interface region of transmembrane channels is thought to be indicative of its importance for insertion and stabilization of the channel in membranes. To further investigate the effects of Trp-->Phe substitution on the structure and function of the gramicidin channel, four analogs of gramicidin A have been synthesized in which the tryptophan residues at positions 9, 11, 13, and 15 are sequentially replaced with phenylalanine. The three-dimensional structure of each viable analog has been determined using a combination of two-dimensional NMR techniques and distance geometry-simulated annealing structure calculations. These phenylalanine analogs adopt a homodimer motif, consisting of two beta6.3 helices joined by six hydrogen bonds at their NH2-termini. The replacement of the tryptophan residues does not have a significant effect on the backbone structure of the channels when compared to native gramicidin A, and only small effects are seen on side-chain conformations. Single-channel conductance measurements have shown that the conductance and lifetime of the channels are significantly affected by the replacement of the tryptophan residues (Wallace, 2000; Becker et al., 1991). The variation in conductance appears to be caused by the sequential removal of a tryptophan dipole, thereby removing the ion-dipole interaction at the channel entrance and at the ion binding site. Channel lifetime variations appear to be related to changing side chain-lipid interactions. This is supported by data relating to transport and incorporation kinetics.
| S-EPMC1305000 | biostudies-literature
Project description:The free energy governing K(+) conduction through gramicidin A channels is characterized by using over 0.1 micros of all-atom molecular dynamics simulations with explicit solvent and membrane. The results provide encouraging agreement with experiments and insights into the permeation mechanism. The free energy surface of K(+), as a function of both axial and radial coordinates, is calculated. Correcting for simulation artifacts due to periodicity and the lack of hydrocarbon polarizability, the calculated single-channel conductance for K(+) ions is 0.8 pS, closer to experiment than any previous calculation. In addition, the estimated single ion dissociation constants are within the range of experimental determinations. The relatively small free energy barrier to ion translocation arises from a balance of large opposing contributions from protein, single-file water, bulk electrolyte, and membrane. Mean force decomposition reveals a remarkable ability of the single-file water molecules to stabilize K(+) by -40 kcal/mol, roughly half the bulk solvation free energy. The importance of the single-file water confirms the conjecture of Mackay et al. [Mackay, D. H. J., Berens, P. H., Wilson, K. R. & Hagler, A. T. (1984) Biophys. J. 46, 229-248]. Ion association with the channel involves gradual dehydration from approximately six to seven water molecules in the first shell, to just two inside the narrow pore. Ion permeation is influenced by the orientation of the single-file water column, which can present a barrier to conduction and give rise to long-range coupling of ions on either side of the pore. Small changes in the potential function, including contributions from electronic polarization, are likely to be sufficient to obtain quantitative agreement with experiments.
| S-EPMC314148 | biostudies-literature
Project description:Using both atomistic and coarse-grained (CG) models, we compute the three-dimensional stress field around a gramicidin A (gA) dimer in lipid bilayers that feature different degrees of negative hydrophobic mismatch. The general trends in the computed stress field are similar at the atomistic and CG levels, supporting the use of the CG model for analyzing the mechanical features of protein/lipid/water interfaces. The calculations reveal that the stress field near the protein-lipid interface exhibits a layered structure with both significant repulsive and attractive regions, with the magnitude of the stress reaching 1000 bar in certain regions. Analysis of density profiles and stress field distributions helps highlight the Trp residues at the protein/membrane/water interface as mechanical anchors, suggesting that similar analysis is useful for identifying tension sensors in other membrane proteins, especially membrane proteins involved in mechanosensation. This work fosters a connection between microscopic and continuum mechanics models for proteins in complex environments and makes it possible to test the validity of assumptions commonly made in continuum mechanics models for membrane mediated processes. For example, using the calculated stress field, we estimate the free energy of membrane deformation induced by the hydrophobic mismatch, and the results for regions beyond the annular lipids are in general consistent with relevant experimental data and previous theoretical estimates using elasticity theory. On the other hand, the assumptions of homogeneous material properties for the membrane and a bilayer thickness at the protein/lipid interface being independent of lipid type (e.g., tail length) appear to be oversimplified, highlighting the importance of annular lipids of membrane proteins. Finally, the stress field analysis makes it clear that the effect of even rather severe hydrophobic mismatch propagates to only about two to three lipid layers, thus putting a limit on the range of cooperativity between membrane proteins in crowded cellular membranes.
| S-EPMC3540266 | biostudies-literature
Project description:Sodium channels are chief proteins involved in electrical signaling in the nervous system, enabling critical functions like heartbeat and brain activity. New high-resolution X-ray structures for bacterial sodium channels have created an opportunity to see how these proteins operate at the molecular level. An important challenge to overcome is establishing relationships between the structures and functions of mammalian and bacterial channels. Bacterial sodium channels are known to exhibit the main structural features of their mammalian counterparts, as well as several key functional characteristics, including selective ion conduction, voltage-dependent gating, pore-based inactivation and modulation by local anesthetic, antiarrhythmic and antiepileptic drugs. Simulations have begun to shed light on each of these features in the past few years. Despite deviations in selectivity signatures for bacterial and mammalian channels, simulations have uncovered the nature of the multiion conduction mechanism associated with Na(+) binding to a high-field strength site established by charged glutamate side chains. Simulations demonstrated a surprising level of flexibility of the protein, showing that these side chains are active participants in the permeation process. They have also uncovered changes in protein structure, leading to asymmetrical collapses of the activation gate that have been proposed to correspond to inactivated structures. These observations offer the potential to examine the mechanisms of state-dependent drug activity, focusing on pore-blocking and pore-based slow inactivation in bacterial channels, without the complexities of inactivation on multiple timescales seen in eukaryotic channels. Simulations have provided molecular views of the interactions of drugs, consistent with sites predicted in mammalian channels, as well as a wealth of other sites as potential new drug targets. In this chapter, we survey the new insights into sodium channel function that have emerged from studies of simpler bacterial channels, which provide an excellent learning platform, and promising avenues for mechanistic discovery and pharmacological development.
| S-EPMC6591151 | biostudies-literature
Project description:Understanding the mechanisms underlying ion channel function from the atomic-scale requires accurate ab initio modelling as well as careful experiments. Here, we present a density functional theory (DFT) study of the ion channel gramicidin A (gA), whose inner pore conducts only monovalent cations and whose conductance has been shown to depend on the side chains of the amino acids in the channel. We investigate the ground state geometry and electronic properties of the channel in vacuum, focusing on their dependence on the side chains of the amino acids. We find that the side chains affect the ground state geometry, while the electrostatic potential of the pore is independent of the side chains. This study is also in preparation for a full, linear scaling DFT study of gA in a lipid bilayer with surrounding water. We demonstrate that linear scaling DFT methods can accurately model the system with reasonable computational cost. Linear scaling DFT allows ab initio calculations with 10,000-100,000 atoms and beyond, and will be an important new tool for biomolecular simulations.
| S-EPMC3808544 | biostudies-literature
Project description:Channel and nonchannel forms of gramicidin A (GA) were studied by ESR in various lipid environments using new mono- and double-spin-labeled compounds. For GA channels, we demonstrate here how pulse dipolar ESR can be used to determine the orientation of the membrane-traversing molecule relative to the membrane normal and to study subtle effects of lipid environment on the interspin distance in the spin-labeled gramicidin channel. To study nonchannel forms of gramicidin, pulse dipolar ESR was used first to determine interspin distances corresponding to monomers and double-helical dimers of spin-labeled GA molecules in the organic solvents trifluoroethanol and octanol. The same distances were then observed in membranes. Since detection of nonchannel forms in the membrane is complicated by aggregation, we suppressed any dipolar spectra from intermolecular interspin distances arising from the aggregates by using double-labeled GA in a mixture with excess unlabeled GA. In hydrophobic mismatching lipids (L(β) phase of DPPC), gramicidin channels dissociate into free monomers. The backbone structure of the monomeric form is similar to a monomeric unit of the channel dimer. In addition to channels and monomers, the double-helical conformation of gramicidin is present in some membrane environments. In the gel phase of saturated phosphatidylcholines, the fraction of double helices increases in the following order: DLPC < DMPC < DSPC < DPPC. The equilibrium DHD/monomer ratio in DPPC was determined. In membranes, the double-helical form is present only in aggregates. In addition, we studied the effect of N-terminal substitution in the GA molecule upon channel formation. This work demonstrates how pulsed dipolar ESR may be utilized to study complex equilibria of peptides in membranes.
| S-EPMC3076037 | biostudies-literature
Project description:The monovalent cation selective channel formed by a dimer of the polypeptide gramicidin A has a single-stranded, right-handed helical motif with 6.5 residues per turn forming a 4-A diameter pore. The structure has been refined to high resolution against 120 orientational constraints obtained from samples in a liquid-crystalline phase lipid bilayer. These structural constraints from solid-state NMR reflect the orientation of spin interaction tensors with respect to a unique molecular axis. Because these tensors are fixed in the molecular frame and because the samples are uniformly aligned with respect to the magnetic field of the NMR spectrometer, each constraint restricts the orientation of internuclear vectors with respect to the laboratory frame of reference. The structural motif of this channel has been validated, and the high-resolution structure has led to precise models for cation binding, cation selectivity, and cation conductance efficiency. The structure is consistent with the electrophysiological data and numerous biophysical studies. Contrary to a recent claim [Burkhart, B. M., Li, N., Langs, D. A., Pangborn, W. A. & Duax, W. L. (1998) Proc. Natl. Acad. Sci. USA 95, 12950-12955], the solid-state NMR constraints for gramicidin A in a lipid bilayer are not consistent with an x-ray crystallographic structure for gramicidin having a double-stranded, right-handed helix with 7.2 residues per turn.
| S-EPMC22161 | biostudies-literature
Project description:Gramicidin A is an antimicrobial peptide that destroys gram-positive bacteria. The bactericidal mechanism of antimicrobial peptides has been linked to membrane permeation and metabolism disruption as well as interruption of DNA and protein functions. However, the exact bacterial killing mechanism of gramicidin A is not clearly understood. In the present study, we examined the antimicrobial activity of gramicidin A on Staphylococcus aureus using biochemical and biophysical methods, including hydroxyl radical and NAD+/NADH cycling assays, atomic force microscopy, and Fourier transform infrared spectroscopy. Gramicidin A induced membrane permeabilization and changed the composition of the membrane. The morphology of Staphylococcus aureus during gramicidin A destruction was divided into four stages: pore formation, water permeability, bacterial flattening, and lysis. Changes in membrane composition included the destruction of membrane lipids, proteins, and carbohydrates. Most interestingly, we demonstrated that gramicidin A not only caused membrane permeabilization but also induced the formation of hydroxyl radicals, which are a possible end product of the transient depletion of NADH from the tricarboxylic acid cycle. The latter may be the main cause of complete Staphylococcus aureus killing. This new finding may provide insight into the underlying bactericidal mechanism of gA.
| S-EPMC4306519 | biostudies-literature