Simulation study of a gramicidin/lipid bilayer system in excess water and lipid. II. Rates and mechanisms of water transport

Inhibition Mechanism of Anti-TB Drug SQ109: Allosteric Inhibition of TMM Translocation of Mycobacterium Tuberculosis MmpL3 Transporter

The mycolic acid transporter MmpL3 is driven by proton motive forces (PMF) and functions via an antiport mechanism. Although the crystal structures of the Mycobacterium smegmatis MmpL3 transporter alone and in complex with a trehalose monomycolate (TMM) substrate and an antituberculosis drug candidate SQ109 under Phase 2b-3 Clinical Trials are available, no water and no conformational change in MmpL3 were observed in these structures to explain SQ109's inhibition mechanism of proton and TMM transportation. In this study, molecular dynamics simulations of both apo form and inhibitor-bound MmpL3 in an explicit membrane were used to decipher the inhibition mechanism of SQ109. In the apo system, the close-open motion of the two TM domains, likely driven by the proton translocation, drives the close-open motion of the two PD domains, presumably allowing for TMM translocation. In contrast, in the holo system, the two PD domains are locked in a closed state, and the two TM domains are locked in an off pathway wider open state due to the binding of the inhibitor. Consistent with the close-open motion of the two PD domains, TMM entry size changes in the apo system, likely loading and moving the TMM, but does not vary much in the holo system and probably impair the movement of the TMM. Furthermore, we observed that water molecules passed through the central channel of the MmpL3 transporter to the cytoplasmic side in the apo system but not in the holo system, with a mean passing time of ∼135 ns. Because water wires play an essential role in transporting protons, our findings shed light on the importance of PMF in driving the close-open motion of the two TM domains. Interestingly, the key channel residues involved in water passage display considerable overlap with conserved residues within the MmpL protein family, supporting their critical function role.

Molecular mechanisms of pore formation and membrane disruption by the antimicrobial lantibiotic peptide Mutacin 1140

The emergence of antibiotic-resistance is a major concern to global human health and identification of novel antibiotics is critical to mitigate the threat.

Characterizing Residue-Bilayer Interactions Using Gramicidin A as a Scaffold and Tryptophan Substitutions as Probes

Previous experiments have shown that the lifetime of a gramicidin A dimer channel (which forms from two nonconducting monomers) in a lipid bilayer is modulated by mutations of the tryptophan (Trp) residues at the bilayer-water interface. We explore this further using extensive molecular dynamics simulations of various gA dimer and monomer mutants at the Trp positions in phosphatidylcholine bilayers with different tail lengths. gA interactions with the surrounding bilayer are strongly modulated by mutating these Trp residues. There are three principal effects: eliminating residue hydrogen bonding ability (i.e., reducing the channel-monolayer coupling strength) reduces the extent of the bilayer deformation caused by the assembled dimeric channel; a residue's size and geometry affects its orientation, leading to different hydrogen bonding partners; and increasing a residue's hydrophobicity increases the depth of gA monomer insertion relative to the bilayer center, thereby increasing the lipid bending frustration.

Water Determines the Structure and Dynamics of Proteins

Water is an essential participant in the stability, structure, dynamics, and function of proteins and other biomolecules. Thermodynamically, changes in the aqueous environment affect the stability of biomolecules. Structurally, water participates chemically in the catalytic function of proteins and nucleic acids and physically in the collapse of the protein chain during folding through hydrophobic collapse and mediates binding through the hydrogen bond in complex formation. Water is a partner that slaves the dynamics of proteins, and water interaction with proteins affect their dynamics. Here we provide a review of the experimental and computational advances over the past decade in understanding the role of water in the dynamics, structure, and function of proteins. We focus on the combination of X-ray and neutron crystallography, NMR, terahertz spectroscopy, mass spectroscopy, thermodynamics, and computer simulations to reveal how water assist proteins in their function. The recent advances in computer simulations and the enhanced sensitivity of experimental tools promise major advances in the understanding of protein dynamics, and water surely will be a protagonist.

Role of amphipathicity and hydrophobicity in the balance between hemolysis and peptide-membrane interactions of three related antimicrobial peptides

Cationic antimicrobial peptides (CAMPs) represent important self defense molecules in many organisms, including humans. These peptides have a broad spectrum of activities, killing or neutralizing many Gram-negative and Gram-positive bacteria. The emergence of multidrug resistant microbes has stimulated research on the development of alternative antibiotics. In the search for new antibiotics, cationic antimicrobial peptides (CAMPs) offer a viable alternative to conventional antibiotics, as they physically disrupt the bacterial membranes, leading to lysis of microbial membranes and eventually cell death. In particular, the group of linear α-helical cationic peptides has attracted increasing interest from clinical as well as basic research during the last decade. In this work, we studied the biophysical and microbiological characteristics of three new designed CAMPs. We modified a previously studied CAMP sequence, in order to increase or diminish the hydrophobic face, changing the position of two lysines or replacing three leucines, respectively. These mutations modified the hydrophobic moment of the resulting peptides and allowed us to study the importance of this parameter in the membrane interactions of the peptides. The structural properties of the peptides were also correlated with their membrane-disruptive abilities, antimicrobial activities and hemolysis of human red blood cells.

Singlet oxygen effects on lipid membranes: implications for the mechanism of action of broad-spectrum viral fusion inhibitors

It was reported recently that a new aryl methyldiene rhodanine derivative, LJ001, and oxazolidine-2,4-dithione, JL103, act on the viral membrane, inhibiting its fusion with a target cell membrane. The aim of the present study was to investigate the interactions of these two active compounds and an inactive analogue used as a negative control, LJ025, with biological membrane models, in order to clarify the mechanism of action at the molecular level of these new broad-spectrum enveloped virus entry inhibitors. Fluorescence spectroscopy was used to quantify the partition and determine the location of the molecules on membranes. The ability of the compounds to produce reactive oxygen molecules in the membrane was tested using 9,10-dimethylanthracene, which reacts selectively with singlet oxygen (1O2). Changes in the lipid packing and fluidity of membranes were assessed by fluorescence anisotropy and generalized polarization measurements. Finally, the ability to inhibit membrane fusion was evaluated using FRET. Our results indicate that 1O2 production by LJ001 and JL103 is able to induce several changes on membrane properties, specially related to a decrease in its fluidity, concomitant with an increase in the order of the polar headgroup region, resulting in an inhibition of the membrane fusion necessary for cell infection.

Improvement of HIV fusion inhibitor C34 efficacy by membrane anchoring and enhanced exposure

OBJECTIVES

The aim of the present work was to evaluate the interaction of two new HIV fusion inhibitors {HIVP3 [C34-polyethylene glycol (PEG)₄-cholesterol] and HIVP4 [(C34-PEG₄)₂-cholesterol]} with membrane model systems and human blood cells in order to clarify where and how the fusion inhibitors locate, allowing us to understand their mechanism of action at the molecular level, and which strategies may be followed to increase efficacy.

METHODS

Lipid vesicles with defined compositions were used for peptide partition and localization studies, based on the intrinsic fluorescence of HIVP3 and HIVP4. Lipid monolayers were employed in surface pressure studies. Finally, human erythrocytes and peripheral blood mononuclear cells (PBMCs) isolated from blood samples were used in dipole potential assays.

RESULTS

Membrane partition, dipole potential and surface pressure assays indicate that the new fusion inhibitors interact preferentially with cholesterol-rich liquid-ordered membranes, mimicking biological membrane microdomains known as lipid rafts. HIVP3 and HIVP4 are able to interact with human erythrocytes and PBMCs to a similar degree as a previously described simpler drug with monomeric C34 and lacking the PEG spacer, C34-cholesterol. However, the pocket-binding domain (PBD) of both HIVP3 and HIVP4 is more exposed to the aqueous environment than in C34-cholesterol.

CONCLUSIONS

The present data allow us to conclude that more efficient blocking of HIV entry results from the synergism between the membranotropic behaviour and the enhanced exposure of the PBD.

Density-functional theory study of gramicidin A ion channel geometry and electronic properties

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.

Biological activity of antibacterial peptides matches synergism between electrostatic and non electrostatic forces

Substitution of Ala 108 and Ala 111 in the 107-115 human lysozyme (hLz) fragment results in a 20-fold increased anti-staphylococcal activity while its hemolytic activity becomes significant (30%) at very high concentrations. This analog displays an additional positive charge near the N-terminus (108) and an extra Trp residue at the center of the molecule (111), indicating that this particular amino acid sequence improves its interaction with the bacterial plasma membrane. In order to understand the role of this arrangement in the membrane interaction, studies with model lipid membranes were carried out. The interactions of peptides, 107-115 hLz and the novel analog ([K(108)W(111)]107-115 hLz) with liposomes and lipid monolayers were evaluated by monitoring the changes in the fluorescence of the Trp residues and the variation of the monolayers surface pressure, respectively. Results obtained with both techniques revealed a significant affinity increase of [K(108)W(111)]107-115 hLz for lipids, especially when the membranes containing negatively charged lipids, such as phosphatidylglycerol. However, there is also a significant interaction with zwitterionic lipids, suggesting that other forces in addition to electrostatic interactions are involved in the binding. The analysis of adsorption isotherms and the insertion kinetics suggest that relaxation processes of the membrane structure are involved in the insertion process of novel peptide [K(108)W(111)]107-115 hLz but not in 107-115 hLz, probably by imposing a reorganization of water at the interphases. In this regard, the enhanced activity of peptide [K(108)W(111)]107-115 hLz may be explained by a synergistic effect between the increased electrostatic forces as well as the increased hydrophobic interactions.

A mechanistic paradigm for broad-spectrum antivirals that target virus-cell fusion

LJ001 is a lipophilic thiazolidine derivative that inhibits the entry of numerous enveloped viruses at non-cytotoxic concentrations (IC₅₀≤0.5 µM), and was posited to exploit the physiological difference between static viral membranes and biogenic cellular membranes. We now report on the molecular mechanism that results in LJ001's specific inhibition of virus-cell fusion.

The antiviral activity of LJ001 was light-dependent, required the presence of molecular oxygen, and was reversed by singlet oxygen (¹O₂) quenchers, qualifying LJ001 as a type II photosensitizer. Unsaturated phospholipids were the main target modified by LJ001-generated ¹O₂. Hydroxylated fatty acid species were detected in model and viral membranes treated with LJ001, but not its inactive molecular analog, LJ025. ¹O₂-mediated allylic hydroxylation of unsaturated phospholipids leads to a trans-isomerization of the double bond and concurrent formation of a hydroxyl group in the middle of the hydrophobic lipid bilayer. LJ001-induced ¹O₂-mediated lipid oxidation negatively impacts on the biophysical properties of viral membranes (membrane curvature and fluidity) critical for productive virus-cell membrane fusion. LJ001 did not mediate any apparent damage on biogenic cellular membranes, likely due to multiple endogenous cytoprotection mechanisms against phospholipid hydroperoxides.

Based on our understanding of LJ001's mechanism of action, we designed a new class of membrane-intercalating photosensitizers to overcome LJ001's limitations for use as an in vivo antiviral agent. Structure activity relationship (SAR) studies led to a novel class of compounds (oxazolidine-2,4-dithiones) with (1) 100-fold improved in vitro potency (IC₅₀<10 nM), (2) red-shifted absorption spectra (for better tissue penetration), (3) increased quantum yield (efficiency of ¹O₂ generation), and (4) 10–100-fold improved bioavailability. Candidate compounds in our new series moderately but significantly (p≤0.01) delayed the time to death in a murine lethal challenge model of Rift Valley Fever Virus (RVFV). The viral membrane may be a viable target for broad-spectrum antivirals that target virus-cell fusion.

The threat of emerging and re-emerging viruses underscores the need to develop broad-spectrum antivirals. LJ001 is a non-cytotoxic, membrane-targeted, broad-spectrum antiviral previously reported to inhibit the entry of many lipid-enveloped viruses. Here, we delineate the molecular mechanism that underlies LJ001's antiviral activity. LJ001 generates singlet oxygen (¹O₂) in the membrane bilayer; ¹O₂-mediated lipid oxidation results in changes to the biophysical properties of the viral membrane that negatively impacts its ability to undergo virus-cell fusion. These changes are not apparent on LJ001-treated cellular membranes due to their repair by cellular lipid biosynthesis. Thus, we generated a new class of membrane-targeted broad-spectrum antivirals with improved photochemical, photophysical, and pharmacokinetic properties leading to encouraging in vivo efficacy against a lethal emerging pathogen. This study provides a mechanistic paradigm for the development of membrane-targeting broad-spectrum antivirals that target the biophysical process underlying virus-cell fusion and that exploit the difference between inert viral membranes and their biogenic cellular counterparts.

Conjugation of cholesterol to HIV-1 fusion inhibitor C34 increases peptide-membrane interactions potentiating its action

Recently, the covalent binding of a cholesterol moiety to a classical HIV-1 fusion inhibitor peptide, C34, was shown to potentiate its antiviral activity. Our purpose was to evaluate the interaction of cholesterol-conjugated and native C34 with membrane model systems and human blood cells to understand the effects of this derivatization. Lipid vesicles and monolayers with defined compositions were used as model membranes. C34-cholesterol partitions more to fluid phase membranes that mimic biological membranes. Importantly, there is a preference of the conjugate for liquid ordered membranes, rich in cholesterol and/or sphingomyelin, as observed both from partition and surface pressure studies. In human erythrocytes and peripheral blood mononuclear cells (PBMC), C34-cholesterol significantly decreases the membrane dipole potential. In PBMC, the conjugate was 14- and 115-fold more membranotropic than T-1249 and enfuvirtide, respectively. C34 or cholesterol alone did not show significant membrane activity. The enhanced interaction of C34-cholesterol with biological membranes correlates with its higher antiviral potency. Higher partitions for lipid-raft like compositions direct the drug to the receptor-rich domains where membrane fusion is likely to occur. This intermediary membrane binding step may facilitate the drug delivery to gp41 in its pre-fusion state.

Design of peptide-membrane interactions to modulate single-file water transport through modified gramicidin channels

Water permeability through single-file channels is affected by intrinsic factors such as their size and polarity and by external determinants like their lipid environment in the membrane. Previous computational studies revealed that the obstruction of the channel by lipid headgroups can be long-lived, in the range of nanoseconds, and that pore-length-matching membrane mimetics could speed up water permeability. To test the hypothesis of lipid-channel interactions modulating channel permeability, we designed different gramicidin A derivatives with attached acyl chains. By combining extensive molecular-dynamics simulations and single-channel water permeation measurements, we show that by tuning lipid-channel interactions, these modifications reduce the presence of lipid headgroups in the pore, which leads to a clear and selective increase in their water permeability.

End-point targeted molecular dynamics: large-scale conformational changes in potassium channels

Large-scale conformational changes in proteins that happen often on biological time scales may be relatively rare events on the molecular dynamics time scale. We have implemented an approach to targeted molecular dynamics called end-point targeted molecular dynamics that transforms proteins between two specified conformational states through the use of nonharmonic "soft" restraints. A key feature of the method is that the protein is free to discover its own conformational pathway through the plethora of possible intermediate states. The method is applied to the Shaker K(v)1.2 potassium channel in implicit solvent. The rate of cycling between the open and closed states was varied to explore how slow the cycling rate needed to be to ensure that microscopic reversibility along the transition pathways was well approximated. Results specific to the K(+) channel include: 1), a variation in backbone torsion angles of residues near the Pro-Val-Pro motif in the inner helix during both opening and closing; 2), the identification of possible occlusion sites in the closed channel located among Pro-Val-Pro residues and downstream; 3), a difference in the opening and closing pathways of the channel; and 4), evidence of a transient intermediate structural substate. The results also show that likely intermediate conformations during the opening-closing process can be generated in computationally tractable simulation times.

The impact of peptides on lipid membranes

We review the fundamental strategies used by small peptides to associate with lipid membranes and how the different strategies impact on the structure and dynamics of the lipids. In particular we focus on the binding of amphiphilic peptides by electrostatic and hydrophobic forces, on the anchoring of peptides to the bilayer by acylation and prenylation, and on the incorporation of small peptides that form well-defined channels. The effect of lipid-peptide interactions on the lipids is characterized in terms of lipid acyl-chain order, membrane thickness, membrane elasticity, permeability, lipid-domain and annulus formation, as well as acyl-chain dynamics. The different situations are illustrated by specific cases for which experimental observations can be interpreted and supplemented by theoretical modeling and simulations. A comparison is made with the effect on lipids of trans-membrane proteins. The various cases are discussed in the context of the possible roles played by lipid-peptide interactions for the biological, physiological, and pharmacological function of peptides.

Effect of Membrane Environment on Proton Permeation through Gramicidin A Channels

Multistate empirical valence bond simulations were employed to study proton transport through gramicidin A channels embedded in two different lipid bilayers, glycerol 1-monooleate (GMO) and diphytanolphosphatidylcholine (DiPhPC). Free energy barriers to proton permeation were derived using a new internal reaction coordinate describing the proton permeation process. The large quantitative and qualitative differences between the two systems are discussed in terms of local bilayer structures, ordering of interfacial water, and channel flexibility in the two environments.

Invariance of single-file water mobility in gramicidin-like peptidic pores as function of pore length

We investigated the structural and energetic determinants underlying water permeation through peptidic nanopores, motivated by recent experimental findings that indicate that water mobility in single-file water channels displays nonlinear length dependence. To address the molecular mechanism determining the observed length dependence, we studied water permeability in a series of designed gramicidin-like channels of different length using atomistic molecular dynamics simulations. We found that within the studied range of length the osmotic water permeability is independent of pore length. This result is at variance with textbook models, where the relationship is assumed to be linear. Energetic analysis shows that loss of solvation rather than specific water binding sites in the pore form the main energetic barrier for water permeation, consistent with our dynamics results. For this situation, we propose a modified expression for osmotic permeability that fully takes into account water motion collectivity and does not depend on the pore length. Different schematic barrier profiles are discussed that explain both experimental and computational interpretations, and we propose a set of experiments aimed at validation of the presented results. Implications of the results for the design of peptidic channels with desired permeation characteristics are discussed.

Free energy simulations of single and double ion occupancy in gramicidin A

Simultaneous occupancy of the two binding sites in gramicidin A by monovalent cations is a well known property of this channel, but the energetic feasibility of this process in molecular dynamics simulations has not been established so far. Here the authors study the energetics of single and double ion occupancy in gramicidin A by constructing the potential of mean force for single and pair of cations. As representatives of small and large ions, they consider both Na+ and K+ ions in the calculations. Binding constants of ions are estimated from the free energy profiles. Comparisons with the experimental results indicate 3-4 kT discrepancy in the binding energies. They also study the coordination of the ions in their respective binding sites and the dynamic behavior of the channel water during the double ion binding process.

Electric field effects on membranes: Gramicidin A as a test ground

Electric fields due to transmembrane potential differences or ionic gradients across the membrane are presumably crucial for many reactions across membranes or close to membranes like signal transduction, control of ion channels or the generation of neural impulses. Molecular dynamics simulations have been used to study the influence of external electric fields on a mixed gramicidin/phospholipid bilayer system. At high field strengths, formation of membrane electropores occurred both close and distal to the gramicidin. Gramicidin was found to stabilize the membrane adjacent to the protein but also at larger distances of up to 2-3 nm. As a result, membrane pore formation was found to be significantly suppressed for the mixed gramicidin/DMPC system. Moderate field strengths only weakly affected the structure and dynamics of the gramicidin. Spontaneous potassium passage events in external electric fields were observed for both the head-to-head helical conformation as well as for the double helical conformation of gramicidin A. The double-helical conformation was found to facilitate ion passage compared to the head-to-head helical dimer.

Effect of quantum partial charges on the structure and dynamics of water in single-walled carbon nanotubes

In this work, using quantum partial charges, computed from 6-31G(**)B3LYP density functional theory, in molecular dynamics simulations, we found that water inside (6,6) and (10,0) single-walled carbon nanotubes with similar diameters but with different chiralities has remarkably different structural and dynamical properties. Density functional calculations indicate that tubes with different chiralities have significantly different partial charges at the ends of tubes. The partial charges at the ends of a (10,0) tube are around 4.5 times higher than those of a (6,6) tube. Molecular dynamics simulations with the partial charges show different water dipole orientations. In the (10,0) tube, dipole vectors of water molecules at the end of the tube point towards the water reservoir resulting in the formation of an L defect in the center region. This is not observed in the (6,6) tube where dipole vectors of all the water molecules inside the tube point towards either the top or the bottom water reservoir. The water diffusion coefficient is found to increase in the presence of the partial charges. Water in the partially charged (10,0) tube has a lower diffusion coefficient compared to that of in the partially charged (6,6) tube.

Chance and design—Proton transfer in water, channels and bioenergetic proteins

Proton transfer and transport in water, gramicidin and some selected channels and bioenergetic proteins are reviewed. An attempt is made to draw some conclusions about how Nature designs long distance, proton transport functionality. The prevalence of water rather than amino acid hydrogen bonded chains is noted, and the possible benefits of waters as the major component are discussed qualitatively.

A hydrophobic gate in an ion channel: the closed state of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is the prototypic member of the 'Cys-loop' superfamily of ligand-gated ion channels which mediate synaptic neurotransmission, and whose other members include receptors for glycine, gamma-aminobutyric acid and serotonin. Cryo-electron microscopy has yielded a three-dimensional structure of the nAChR in its closed state. However, the exact nature and location of the channel gate remains uncertain. Although the transmembrane pore is constricted close to its center, it is not completely occluded. Rather, the pore has a central hydrophobic zone of radius about 3 A. Model calculations suggest that such a constriction may form a hydrophobic gate, preventing movement of ions through a channel. We present a detailed and quantitative simulation study of the hydrophobic gating model of the nicotinic receptor, in order to fully evaluate this hypothesis. We demonstrate that the hydrophobic constriction of the nAChR pore indeed forms a closed gate. Potential of mean force (PMF) calculations reveal that the constriction presents a barrier of height about 10 kT to the permeation of sodium ions, placing an upper bound on the closed channel conductance of 0.3 pS. Thus, a 3 A radius hydrophobic pore can form a functional barrier to the permeation of a 1 A radius Na+ ion. Using a united-atom force field for the protein instead of an all-atom one retains the qualitative features but results in differing conductances, showing that the PMF is sensitive to the detailed molecular interactions.

Mobility of a one-dimensional confined file of water molecules as a function of file length

Confinement of water by pore geometry to a one-dimensional file of molecules interacting with the pore alters the diffusion coefficient D(W). Here we report an exponential dependence of D(W) on the number of water positions in the pore. The result is based on measurements of single channel water permeabilities of structurally similar peptidic nanopores of different length. The inconsistency with predictions from continuum or kinetic models indicates that pore occupancy is reduced in single file transport. In longer pores (e.g., in aquaporins) the presence of charged residues increases D(W).

Single-channel water permeabilities of Escherichia coli aquaporins AqpZ and GlpF

From equilibrium molecular dynamics simulations we have determined single-channel water permeabilities for Escherichia coli aquaporin Z (AqpZ) and aquaglyceroporin GlpF with the channels embedded in lipid bilayers. GlpF's osmotic water permeability constant pf exceeds by 2-3 times that of AqpZ and the diffusive permeability constant (pd) of GlpF is found to exceed that of AqpZ 2-9-fold. Achieving complete water selectivity in AqpZ consequently implies lower transport rates overall relative to the less selective, wider channel of GlpF. For AqpZ, the ratio pf/pd congruent with 12 is close to the average number of water molecules in the channel lumen, whereas for GlpF, pf/pd congruent with 4. This implies that single-file structure of the luminal water is more pronounced for AqpZ, the narrower channel of the two. Electrostatics profiles across the pore lumens reveal that AqpZ significantly reinforces water-channel interactions, and weaker water-water interactions in turn suppress water-water correlations relative to GlpF. Consequently, suppressed water-water correlations across the narrow selectivity filter become a key structural determinant for water permeation causing luminal water to permeate slower across AqpZ.

Structural Determinants of Proton Blockage in Aquaporins

Aquaporins are an important class of membrane channels selective for water and linear polyols but impermeable to ions, including protons. Recent computational studies have revealed that the relay of protons through the water-conduction pathway of aquaporin channels is opposed by a substantial free energy barrier peaking at the signature NPA motifs. Here, free-energy simulations and continuum electrostatic calculations are combined to examine the nature and the magnitude of the contribution of specific structural elements to proton blockage in the bacterial glycerol uptake facilitator, GlpF. Potential of mean-force profiles for both hop and turn steps of structural diffusion in the narrow pore are obtained for artificial variants of the GlpF channel in which coulombic interactions between the pore contents and conserved residues Asn68 and Asn203 at the NPA signature motifs, Arg206 at the selectivity filter, and the peptidic backbone of the two half-helices M3 and M7, which are arranged in head-to-head fashion around the NPA motifs, are turned off selectively. A comparison of these results with electrostatic energy profiles for the translocation of a probe cation throughout the water permeation pathway indicates that the free-energy profile for proton movement inside the narrow pore is dominated by static effects arising from the distribution of charged and polar groups of the channel, whereas dielectric effects contribute primarily to opposing the access of H+ to the pore mouths (desolvation penalty). The single most effective way to abolish the free-energy gradients opposing the movement of H+ around the NPA motif is to turn off the dipole moments of helices M3 and M7. Mutation of either of the two NPA Asn residues to Asp compensates for charge-dipole and dipole-dipole effects opposing the hop and turn steps of structural diffusion, respectively, and dramatically reduces the free energy barrier of proton translocation, suggesting that these single mutants could leak protons.

Molecular dynamics investigation of an oriented cyclic peptide nanotube in DMPC bilayers

Nanotubes resulting from the self-assembly of cyclic peptides formed by eight alpha-amino acids and inserted into lipid bilayers have been shown to function as synthetic, integral transmembrane channels. A nanotube consisting of eight cyclo[(L-Trp-D-Leu)(3)-L-Gln-D-Leu] subunits, organized in an antiparallel, beta-sheetlike channel embedded in a hydrated dimyristoylphosphatidylcholine bilayer was investigated in an 8-ns molecular dynamics trajectory. This large-scale statistical simulation brings to light not only the atomic-level structural features of the synthetic channel, but also its dynamical properties. Overall, the nanotube conserves its hollow tubular structure. The calculation reproduces the tilt of the channel with respect to the normal of the bilayer, in reasonable agreement with experiment. The results show a dislocation of the nanotube indicative of a possible disassembly process that may influence the channel conduction. The dynamics of the water in the hollow tubular structure has been characterized, and the conductance of the channel has been estimated. Transport properties of the peptide nanotube are discussed in comparison with other transporters.

Beyond the diffusion limit: Water flow through the empty bacterial potassium channel

Water molecules are constrained to move with K+ ions through the narrow part of the Streptomyces lividans K+ channel because of the single-file nature of transport. In the presence of an osmotic gradient, a water molecule requires <10 ps to cross the purified protein reconstituted into planar bilayers. Rinsing K+ out of the channel, water may be 1,000 times faster than the fastest experimentally observed K+ ion and 20 times faster than the one-dimensional bulk diffusion of water. Both the anomalously high water mobility and its inhibition observed at high K+ concentrations are consistent with the view that liquid-vapor oscillations occur because of geometrical confinements of water in the selectivity filter. These oscillations, where the chain of molecules imbedded in the channel (the "liquid") cooperatively exits the channel, leaving behind a near vacuum (the "vapor"), thus far have only been discovered in hydrophobic nanopores by molecular dynamics simulations [Hummer, G., Rasaiah, J. C. & Noworyta, J. P. (2001) Nature 414, 188-190; and Beckstein, O. & Sansom, M. S. P. (2003) Proc. Natl. Acad. Sci. USA 100, 7063-7068].

Length, time, and energy scales of photosystems

The design of photosynthetic systems reflects the length scales of the fundamental physical processes. Energy transfer is rapid at the few angstrom scale and continues to be rapid even at the 50-A scale of the membrane thickness. Electron tunneling is nearly as rapid at the shortest distances, but becomes physiologically too slow well before 20 A. Diffusion, which starts out at a relatively slow nanosecond time scale, has the most modest slowing with distance and is physiologically competent at all biologically relevant distances. Proton transfer always operates on the shortest angstrom scale. The structural consequences of these distance dependencies are that energy transfer networks can extend over large, multisubunit and multicomplex distances and take leaps of 20 A before entering the domain of charge separating centers. Electron transfer systems are effectively limited to individual distances of 15 A or less and span the 50 A dimensions of the bioenergetic membrane by use of redox chains. Diffusion processes are generally used to cover the intercomplex electron transfer distances of 50 A and greater and tend to compensate for the lack of directionality by restricting the diffusional space to the membrane or the membrane surface, and by multiplying the diffusing species through the use of pools. Proton transfer reactions act over distances larger than a few angstroms through the use of clusters or relays, which sometimes rely on water molecules and which may only be dynamically assembled. Proteins appear to place a premium on robustness of design, which is relatively easily achieved in the long-distance physical processes of energy transfer and electron tunneling. By placing cofactors close enough, the physical process is relatively rapid compared to decay processes. Thus suboptimal conditions such as cofactor orientation, energy level, or redox potential level can be tolerated and generally do not have to be finely tuned. The most fragile regions of design tend to come in areas of complex formation and catalysis involving proton management, where relatively small changes in distance or mutations can lead to a dramatic decrease in turnover, which may already be limiting the overall speed of energy conversion in these proteins. Light-activated systems also face a challenge to robust function from the ever-present dangers of high redox potential chemistry. This can turn the protein matrix and wandering oxygen molecules into unintentional redox partners, which in the case of PSII requires the frequent, costly replacement of protein subunits.

Molecular dynamics simulations of the bacterial outer membrane protein FhuA: a comparative study of the ferrichrome-free and bound states

FhuA is one of the more complex members of the superfamily of bacterial outer membrane proteins. Its primary function is to provide a binding site on the outer membrane surface for siderophores, such as ferrichrome, and subsequently to facilitate their energy-dependent transport across the membrane, presumably powered by the TonB-ExbBD protein complex that resides in the cytoplasmic membrane. Crystal structures of FhuA with and without a bound ferrichrome molecule have provided some clues as to the initial stages of the siderophore transport mechanism. In the current study, we have employed 10-ns duration molecular dynamics simulations of FhuA and of the FhuA-ferrichrome complex, both embedded in a phospholipid bilayer, to probe the short timescale dynamics of this integral membrane protein, and to explore possible mechanistic implications of this dynamic behavior. Analysis of the dynamics of the protein suggests that the extracellular loops move as relatively rigid entities relative to the transmembrane beta-barrel. Comparison of the two simulations (with and without bound ferrichrome) revealed some ligand-induced changes in loop mobility. Specifically, loop L8 appears to be involved in a mechanism whereby the binding site is gated closed upon ligand binding. Analysis of the dynamics of water molecules within the core of the FhuA protein provided no evidence for a water-permeable protopore through which the ferrichrome might pass without a major perturbation of the FhuA protein. Overall, these simulations support the proposal that binding of ferrichrome initiates a signaling mechanism that ultimately leads to the TonB-mediated partial or total removal of the core domain from the beta-barrel, thus opening up a permeable pore. These simulations are among the longest that have been performed on a complex membrane protein. However, a simple analysis of sampling reveals that the description of protein motions is far from complete.

Voltage-gated proton channels and other proton transfer pathways

Proton channels exist in a wide variety of membrane proteins where they transport protons rapidly and efficiently. Usually the proton pathway is formed mainly by water molecules present in the protein, but its function is regulated by titratable groups on critical amino acid residues in the pathway. All proton channels conduct protons by a hydrogen-bonded chain mechanism in which the proton hops from one water or titratable group to the next. Voltage-gated proton channels represent a specific subset of proton channels that have voltage- and time-dependent gating like other ion channels. However, they differ from most ion channels in their extraordinarily high selectivity, tiny conductance, strong temperature and deuterium isotope effects on conductance and gating kinetics, and insensitivity to block by steric occlusion. Gating of H(+) channels is regulated tightly by pH and voltage, ensuring that they open only when the electrochemical gradient is outward. Thus they function to extrude acid from cells. H(+) channels are expressed in many cells. During the respiratory burst in phagocytes, H(+) current compensates for electron extrusion by NADPH oxidase. Most evidence indicates that the H(+) channel is not part of the NADPH oxidase complex, but rather is a distinct and as yet unidentified molecule.

Peroxyl radicals promoted changes in water permeability through gramicidin channels in DPPC and lecithin-PC vesicles

Gramicidin incorporation to DPPC or lecithin-PC large unilamellar vesicles (LUVs) leads to pore formation that, under hyper-osmotic conditions, produces a noticeable increase in the rate of trans-membrane water flow. This pore formation is more efficient in the more fluid lecithin-PC LUVs. Exposure of these vesicles to peroxyl radicals generated in the aerobic thermolysis of 2,2'-azo-bis(2-amidinopropane) (AAPH), changes the physical properties of the bilayer (as sensed employing fluorescent probes), modifies gramicidin molecules (as sensed by the decrease in Trp fluorescence) and notably reduces the transbilayer rate of water outflow. In order to evaluate if this reduced water-transport capacity is due to changes in the membrane due to lipid-peroxidation and/or direct damage to gramicidin channels, results obtained in the oxidable vesicles (lecithin-PC) were compared to those obtained in DPPC vesicles. The data obtained show that most of the water transport efficiency loss can be ascribed to a direct disruption of gramicidin channels by AAPH derived peroxyl radicals.

Pores formed by the nicotinic receptor m2delta Peptide: a molecular dynamics simulation study

The M2delta peptide self-assembles to form a pentameric bundle of transmembrane alpha-helices that is a model of the pore-lining region of the nicotinic acetylcholine receptor. Long (>15 ns) molecular dynamics simulations of a model of the M2delta(5) bundle in a POPC bilayer have been used to explore the conformational dynamics of the channel assembly. On the timescale of the simulation, the bundle remains relatively stable, with the polar pore-lining side chains remaining exposed to the lumen of the channel. Fluctuations at the helix termini, and in the helix curvature, result in closing/opening transitions at both mouths of the channel, on a timescale of approximately 10 ns. On average, water within the pore lumen diffuses approximately 4x more slowly than water outside the channel. Examination of pore water trajectories reveals both single-file and path-crossing regimes to occur at different times within the simulation.

Ion permeation through the gramicidin channel: atomically detailed modeling by the Stochastic Difference Equation

Atomically detailed descriptions of ionic solution, membrane, and the gramicidin channel are used to compute molecular dynamics trajectories of ion permeation. The microsecond trajectories are calculated with the Stochastic Difference Equation (SDE), which provides approximate solutions to the equations of motions (with filtered high-frequency modes) of extended timescales. The relative permeations of lithium, sodium, and potassium are estimated by using a novel, kinetic cycle protocol and are compared with experiment. The transport through native gramicidin and one fluoro-valine variant is considered as well. Qualitative agreement between theory and experiment is obtained. The faster permeation rate of sodium compared to lithium is reproduced in the calculations. The calculations also reproduce the slower diffusion through a gramicidin with fluorinated valine compared to native gramicidin. The calculations are inconclusive about the relative rates of potassium and sodium. The experiment suggests that potassium permeates more quickly. We directly probe the kinetics of a biophysical process at a relevant time window without reducing the atomically detailed description of the system. The calculations were able to capture subtle balances between binding and diffusion that determine permeation rates. The same model gave the correct ordering of diffusion rates for cases in which electrostatic binding has opposite effects and must be supplemented by dynamic factors. Diffusion rates are faster when favorable electrostatic interactions of ions in the channel (compared to the solvent) are observed. Studies of a gramicidin variant suggest an opposite effect, in which permeation is faster for the less polar channel, indicating dynamic effects. Although both trends can be explained qualitatively, it is not possible to predict (before doing the SDE calculations) which factor is more important.

Large-scale molecular dynamics simulations of general anesthetic effects on the ion channel in the fully hydrated membrane: the implication of molecular mechanisms of general anesthesia

Interactions of volatile anesthetics with the central nervous system are characterized by low yet specific binding affinities. Although neurotransmitter-gated ion channels are considered the primary anesthetic targets, the mechanism of action at the molecular level remains elusive. We consider here the theoretical implications of channel dynamics on anesthetic action in a simplified membrane-channel system. Large-scale 2.2-ns all-atom molecular dynamics simulations were performed to study the effects of halothane, a clinical anesthetic, on a gramicidin A (gA) channel in a fully hydrated dimyristoyl phosphatidylcholine membrane. In agreement with experimental results, anesthetics preferentially target the anchoring residues at the channel-lipid-water interface. Although the anesthetic effect on channel structure is minimal, the presence of halothane profoundly affects channel dynamics. For 2.2-ns simulation, the rms fluctuation of gA backbone in the lipid core increases from approximately equal 1 A in the absence of anesthetics to approximately equal 1.5 A in the presence of halothane. Autocorrelation analysis reveals that halothane (i) has no effect on the subpicosecond librational motion, (ii) prolongs the backbone autocorrelation time in the 10- to 100-ps time scale, and (iii) significantly decreases the asymptotic values of generalized order parameter and correlation time of nanosecond motions for the inner but not the outer residues. The simulation results discount the viewpoint of a structure-function paradigm that overrates the importance of structural fitting between general anesthetics and yet-unidentified hydrophobic protein pockets. Instead, the results underscore the global, as opposed to local, effects of anesthetics on protein dynamics as the underlying mechanisms for the action of general anesthetics and possibly of other low-affinity drugs.

The role of Trp side chains in tuning single proton conduction through gramicidin channels

We present an extensive set of measurements of proton conduction through gramicidin A (gA), B (gB), and M (gM) homodimer channels which have 4, 3, or 0 Trp residues at each end of the channel, respectively. In gA we find a shoulder separating two domains of conductance increasing with concentration, confirming the results of Eisenman, G., B. Enos, J. Hagglund, and J. Sandblom. 1980. Ann. NY. Acad. Sci. 339:8-20. In gB, the shoulder is shifted by approximately 1/2 pH unit to higher H(+) concentrations and is very sharply defined. No shoulder appears in the gM data, but an associated transition from sublinear to superlinear I-V values occurs at a 100-fold higher [H(+)] in gM than in gA. The data in the low concentration domain are analyzed using a configuration space model of single-proton conduction, assuming that the difference in the proton potential of mean force (PMF) between gA and its analogs is constant, similar to the results of Anderson, D., R. B. Shirts, T. A. Cross, and D. D. Busath. 2001. Biophys. J. 81:1255-1264. Our results suggest that the average amplitudes of the calculated proton PMFs are nearly correct, but that the water reorientation barrier calculated for gA by molecular dynamics using the PM6 water model (Pomès, R., and B. Roux. 1997. Biophys. J. 72:246a) must be reduced in amplitude by 1.5 kcal/mol or more, and is not rate-limiting for gA.

Water permeation through gramicidin A: desformylation and the double helix: a molecular dynamics study

Multinanosecond molecular dynamics simulations of gramicidin A embedded in a dimyristoylphosphatidylcholine bilayer show a remarkable structural stability for both experimentally determined conformations: the head-to-head helical dimer and the double helix. Water permeability was found to be much higher in the double helical conformation, which is explained by lower hydrogen bond-mediated enthalpic barriers at the channel entrance and its larger pore size. Free-energy perturbation calculations show that the double helical structure is stabilized by the positive charges at the N termini introduced by the desformylation, whereas the helical dimer is destabilized. Together with the recent experimental observation that desformyl gramicidin conducts water hundredfold better than gramicidin, this suggests that desformyl gramicidin A predominantly occurs in the double helical conformation.

Molecular dynamics simulations of wild-type and mutant forms of the Mycobacterium tuberculosis MscL channel

The crystal structure of the Mycobacterium tuberculosis homolog of the bacterial mechanosensitive channel of large conductance (Tb-MscL) provides a unique opportunity to consider mechanosensitive signal transduction at the atomic level. Molecular dynamics simulations of the Tb-MscL channel embedded in an explicit lipid bilayer and of its C-terminal helical bundle alone in aqueous solvent were performed. C-terminal calculations imply that although the helix bundle structure is relatively unstable at physiological pH, it may have been stabilized under low pH conditions such as those used in the crystallization of the channel. Specific mutations to the C-terminal region, which cause a similar conservation of the crystal structure conformation, have also been identified. Full channel simulations were performed for the wild-type channel and two experimentally characterized gain-of-function mutants, V21A and Q51E. The wild-type Tb-MscL trajectory gives insight into regions of relative structural stability and instability in the channel structure. Channel mutations led to observable changes in the trajectories, such as an alteration of intersubunit interactions in the Q51E mutant. In addition, interesting patterns of protein-lipid interactions, such as hydrogen bonding, arose in the simulations. These and other observations from the simulations are relevant to previous and ongoing experimental studies focusing on characterization of the channel.

Electrostatics calculations: recent methodological advances and applications to membranes

Calculations of electrostatic energies and forces are at the heart of the theoretical modeling of biological molecules. During the past year, new methods for accurately treating electrostatic interactions have been developed for all-atom simulations and for modeling systems in which the uninteresting part of the system, for example, the solvent, is represented implicitly. Recent applications of electrostatic energy calculations have revealed new principles concerning the role of electrostatics in peptide binding to membranes and ion transport across membranes through protein channels.

Potassium and sodium ions in a potassium channel studied by molecular dynamics simulations

We have performed simulations of both a single potassium ion and a single sodium ion within the pore of the bacterial potassium channel KcsA. For both ions there is a dehydration energy barrier at the cytoplasmic mouth suggesting that the crystal structure is a closed conformation of the channel. There is a potential energy barrier for a sodium ion in the selectivity filter that is not seen for potassium. Radial distribution functions for both ions with the carbonyl oxygens of the selectivity filter indicate that sodium may interact more tightly with the filter than does potassium. This suggests that the key to the ion selectivity of KcsA is the greater dehydration energy of Na(+) ions, and helps to explain the block of KcsA by internal Na(+) ions.

Desformylgramicidin: a model channel with an extremely high water permeability

The water conductivity of desformylgramicidin exceeds the permeability of gramicidin A by two orders of magnitude. With respect to its single channel hydraulic permeability coefficient of 1.1.10(-12) cm(3) s(-1), desformylgramicidin may serve as a model for extremely permeable aquaporin water channel proteins (AQP4 and AQPZ). This osmotic permeability exceeds the conductivity that is predicted by the theory of single-file transport. It was derived from the concentration distributions of both pore-impermeable and -permeable cations that were simultaneously measured by double barreled microelectrodes in the immediate vicinity of a planar bilayer. From solvent drag experiments, approximately five water molecules were found to be transported by a single-file process along with one ion through the channel. The single channel proton, potassium, and sodium conductivities were determined to be equal to 17 pS (pH 2.5), 7 and 3 pS, respectively. Under any conditions, the desformyl-channel remains at least 10 times longer in its open state than gramicidin A.

Molecular dynamics simulations predict a tilted orientation for the helical region of dynorphin A(1-17) in dimyristoylphosphatidylcholine bilayers

The structural properties of the endogenous opioid peptide dynorphin A(1-17) (DynA), a potential analgesic, were studied with molecular dynamics simulations in dimyristoylphosphatidylcholine bilayers. Starting with the known NMR structure of the peptide in dodecylphosphocholine micelles, the N-terminal helical segment of DynA (encompassing residues 1-10) was initially inserted in the bilayer in a perpendicular orientation with respect to the membrane plane. Parallel simulations were carried out from two starting structures, systems A and B, that differ by 4 A in the vertical positioning of the peptide helix. The complex consisted of approximately 26,400 atoms (dynorphin + 86 lipids + approximately 5300 waters). After >2 ns of simulation, which included >1 ns of equilibration, the orientation of the helical segment of DynA had undergone a transition from parallel to tilted with respect to the bilayer normal in both the A and B systems. When the helix axis achieved a approximately 50 degrees angle with the bilayer normal, it remained stable for the next 1 ns of simulation. The two simulations with different starting points converged to the same final structure, with the helix inserted in the bilayer throughout the simulations. Analysis shows that the tilted orientation adopted by the N-terminal helix is due to specific interactions of residues in the DynA sequence with phospholipid headgroups, water, and the hydrocarbon chains. Key elements are the "snorkel model"-type interactions of arginine side chains, the stabilization of the N-terminal hydrophobic sequence in the lipid environment, and the specific interactions of the first residue, Tyr. Water penetration within the bilayer is facilitated by the immersed DynA, but it is not uniform around the surface of the helix. Many water molecules surround the arginine side chains, while water penetration near the helical surface formed by hydrophobic residues is negligible. A mechanism of receptor interaction is proposed for DynA, involving the tilted orientation observed from these simulations of the peptide in the lipid bilayer.

Simulations of ion channels--watching ions and water move

Ion channels mediate electrical excitability in neurons and muscle. Three-dimensional structures for model peptide channels and for a potassium (K+) channel have been combined with computer simulations to permit rigorous exploration of structure-function relations of channels. Water molecules and ions within transbilayer pores tend to diffuse more slowly than in bulk solutions. In the narrow selectivity filter of the bacterial K+ channel (i.e. the region of the channel that discriminates between different species of ions) a column of water molecules and K+ ions moves in a concerted fashion. By combining atomistic simulations (in which all atoms of the channel molecule, water and ions are treated explicitly) with continuum methods (in which the description of the channel system is considerably simplified) it is possible to simulate some of the physiological properties of channels.

Solvent drag across gramicidin channels demonstrated by microelectrodes

The competition of ion and water fluxes across gramicidin channels was assessed from the concentration distributions of both pore-impermeable and -permeable cations that were simultaneously measured by double-barreled microelectrodes in the immediate vicinity of a planar bilayer. Because water movement across the membrane led to accumulation of solutes on one side of the membrane and depletion on the other, the permeable cation was not only pushed by water across the channel (true solvent drag); it also flowed along its concentration gradient (pseudo-solvent drag). For the demonstration of true solvent drag, a difference between the bulk concentrations on the hypertonic and the hypotonic sides of the membrane was established. It was adjusted to get equal cation concentrations at both membrane/water interfaces. From the sodium and potassium fluxes measured along with membrane conductivity under these conditions, approximately five water molecules were found to be transported simultaneously with one ion through the channel. In diphytanoyl phosphatidylcholine membranes, a single-channel hydraulic permeability coefficient of 1.6 x 10(-14) cm(3) s(-1) was obtained.