Structure of gramicidin a in a lipid bilayer environment determined using molecular dynamics simulations and solid-state NMR data

The gramicidin channel ion permeation free-energy profile: direct and indirect effects of CHARMM force field improvements

A revised CHARMM force field for tryptophan residues is studied as well as a new grid-based correction algorithm, called CMAP, using molecular dynamics simulations of gramicidin A (1JNO) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl or KCl saline solution. The conformational stability of the interfacial side chains is studied, which shows good stability on the 10 ns time scale. The revised force field for the tryptophan side chain produces, in the decomposition, a Na(+) PMF(Trp) profile that is consonant with the prediction from the experimental results, analyzed with rate theory by Durrant et al. (2006), but in stark contrast to the prediction of the original CHARMM force field, version 22. However, the effect is diluted in the PMF profile due to indirect effects mediated by other components of the system (polypeptide, lipid molecules, ions, and water molecules). CMAP corrections to the L-amino acids help reduce the excessive translocation barrier. Decomposition demonstrates that this effect is due to effects on the K(+) PMF(H(2)O) profile rather than on the K(+) PMF(gA) profile. The results have been confirmed to be robust using an alternative umbrella-potential method. Further force field balancing efforts (direct and indirect) are required for future studies to evaluate whether these effects give rise to predictions that are consistent with those observables extracted from real experiments.

Sensitivity enhancement of separated local field experiments: application to membrane proteins

Separated local field (SLF) experiments have been used for almost three decades to obtain structural information in solid-state NMR. These experiments resolve chemical shift anisotropy (CSA) from dipole-dipole interactions (dipolar couplings, DC) in isolated spin systems. Both CSA and DC data can be converted into orientational constraints to elucidate the secondary structure and topology of membrane proteins in oriented lipid bilayers. Here, we propose a new suite of sensitivity enhanced SLF pulse sequences to measure CSA and DC for aligned membrane proteins and liquid crystalline molecules that will decrease the time needed for data acquisition. We demonstrate the efficacy of these new sensitivity enhanced experiments using both a single crystal of N-acetyl leucine and a single pass membrane protein sarcolipin reconstituted in aligned lipid bicelles. These results lay the groundwork for the routine application of this methodology for studying the structure and topology of membrane proteins.

Systematic study of anharmonic features in a principal component analysis of gramicidin A

We use principal component analysis (PCA) to detect functionally interesting collective motions in molecular-dynamics simulations of membrane-bound gramicidin A. We examine the statistical and structural properties of all PCA eigenvectors and eigenvalues for the backbone and side-chain atoms. All eigenvalue spectra show two distinct power-law scaling regimes, quantitatively separating large from small covariance motions. Time trajectories of the largest PCs converge to Gaussian distributions at long timescales, but groups of small-covariance PCs, which are usually ignored as noise, have subdiffusive distributions. These non-Gaussian distributions imply anharmonic motions on the free-energy surface. We characterize the anharmonic components of motion by analyzing the mean-square displacement for all PCs. The subdiffusive components reveal picosecond-scale oscillations in the mean-square displacement at frequencies consistent with infrared measurements. In this regime, the slowest backbone mode exhibits tilting of the peptide planes, which allows carbonyl oxygen atoms to provide surrogate solvation for water and cation transport in the channel lumen. Higher-frequency modes are also apparent, and we describe their vibrational spectra. Our findings expand the utility of PCA for quantifying the essential features of motion on the anharmonic free-energy surface made accessible by atomistic molecular-dynamics simulations.

Lipid bilayer regulation of membrane protein function: gramicidin channels as molecular force probes

Membrane protein function is regulated by the host lipid bilayer composition. This regulation may depend on specific chemical interactions between proteins and individual molecules in the bilayer, as well as on non-specific interactions between proteins and the bilayer behaving as a physical entity with collective physical properties (e.g. thickness, intrinsic monolayer curvature or elastic moduli). Studies in physico-chemical model systems have demonstrated that changes in bilayer physical properties can regulate membrane protein function by altering the energetic cost of the bilayer deformation associated with a protein conformational change. This type of regulation is well characterized, and its mechanistic elucidation is an interdisciplinary field bordering on physics, chemistry and biology. Changes in lipid composition that alter bilayer physical properties (including cholesterol, polyunsaturated fatty acids, other lipid metabolites and amphiphiles) regulate a wide range of membrane proteins in a seemingly non-specific manner. The commonality of the changes in protein function suggests an underlying physical mechanism, and recent studies show that at least some of the changes are caused by altered bilayer physical properties. This advance is because of the introduction of new tools for studying lipid bilayer regulation of protein function. The present review provides an introduction to the regulation of membrane protein function by the bilayer physical properties. We further describe the use of gramicidin channels as molecular force probes for studying this mechanism, with a unique ability to discriminate between consequences of changes in monolayer curvature and bilayer elastic moduli.

CHARMM: the biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

Molecular resolution imaging of an antibiotic peptide in a lipid matrix

In this work, we show molecular resolution scanning tunneling microscopy (STM) images of gramicidin, a model antibacterial peptide, inserted into a phospholipid matrix. The resolution of the images is superior to that obtained in previous attempts to image gramicidin in a lipid environment using atomic force microscopy (AFM). This breakthrough has allowed visualization of individual peptide molecules surrounded by lipid molecules. We have observed several important features: the peptide molecules do not aggregate, the peptide molecules adopt a single conformation corresponding to a specific ion channel form, and the lipid molecules adjacent to the peptide molecules are systematically longer than those in the lipid matrix. These results constitute a new approach to obtain structural characteristics of antibiotic peptides in lipid assemblies that is necessary for the understanding of their biological activity.

Tilt and azimuthal angles of a transmembrane peptide: a comparison between molecular dynamics calculations and solid-state NMR data of sarcolipin in lipid membranes

We report molecular dynamics simulations in the explicit membrane environment of a small membrane-embedded protein, sarcolipin, which regulates the sarcoplasmic reticulum Ca-ATPase activity in both cardiac and skeletal muscle. In its monomeric form, we found that sarcolipin adopts a helical conformation, with a computed average tilt angle of 28 +/- 6 degrees and azymuthal angles of 66 +/- 22 degrees, in reasonable accord with angles determined experimentally (23 +/- 2 degrees and 50 +/- 4 degrees, respectively) using solid-state NMR with separated-local-field experiments. The effects of time and spatial averaging on both (15)N chemical shift anisotropy and (1)H/(15)N dipolar couplings have been analyzed using short-time averages of fast amide out-of-plane motions and following principal component dynamic trajectories. We found that it is possible to reproduce the regular oscillatory patterns observed for the anisotropic NMR parameters (i.e., PISA wheels) employing average amide vectors. This work highlights the role of molecular dynamics simulations as a tool for the analysis and interpretation of solid-state NMR data.

Computational studies of gramicidin permeation: an entry way sulfonate enhances cation occupancy at entry sites

The impact on the cation-transport free-energy profile of replacing the C-terminal ethanolamine in the gramicidin A channel with a taurine residue is studied using molecular dynamics simulations of gramicidin A (1JNO) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl saline solution. The potential of mean force for ion transport is obtained by umbrella sampling. The presence of a negatively charged sulfonate group at the entrance of the gramicidin channel affects the depth and the location of the binding sites, producing a strong attraction for the cations in the bulk. The potential of mean force by the sulfonate acting directly through electrostatics and van der Waals interactions on the test ion is highly modulated by indirect effects (i.e., sulfonate effects on other components of the system that, in turn, affect the ion free-energy profile in the channel). Because the "entry" sites are located symmetrically at both entry and exit of the channel, the deeper free-energy wells should inhibit exit. Given that the channel has increased conductance experimentally, the simulation results suggest that the channel conductance is normally entry limited.

Understanding the structural and functional differences between mouse thyrotropin-releasing hormone receptors 1 and 2

Multiple computational methods have been employed in a comparative study of thyrotropin-releasing hormone receptors 1 and 2 (TRH-R1 and TRH-R2) to explore the structural bases for the different functional properties of these G protein-coupled receptors. Three-dimensional models of both murine TRH receptors have been built and optimized by means of homology modeling based on the crystal structure of bovine rhodopsin, molecular dynamics simulations, and energy minimizations in a membrane-aqueous environment. The comparison between the two models showed a correlation between the higher flexibility and higher basal activity of TRH-R2 versus the lesser flexibility and lower basal activity of TRH-R1 and supported the involvement of the highly conserved W6.48 in the signaling process. A correlation between the level of basal activity and conformational changes of TM5 was detected also. Comparison between models of the wild type receptors and their W6.48A mutants, which have reversed basal activities compared with their respective wild types, further supported these correlations. A flexible molecular docking procedure revealed that TRH establishes a direct interaction with W6.48 in TRH-R2 but not in TRH-R1. We designed and performed new mutagenesis experiments that strongly supported these observations.

Role of tryptophan residues in gramicidin channel organization and function

The linear peptide gramicidin forms prototypical ion channels specific for monovalent cations and has been used extensively to study the organization, dynamics, and function of membrane-spanning channels. The tryptophan residues in gramicidin channels are crucial for maintaining the structure and function of the channel. We explored the structural basis for the reduction in channel conductance in the case of single-tryptophan analogs of gramicidin with three Trp --> hydrophobic substitutions using a combination of fluorescence approaches, which include red edge excitation shift and membrane penetration depth analysis, size-exclusion chromatography, and circular dichroism spectroscopy. We show here that the gramicidin analogs containing single-tryptophan residues adopt a mixture of nonchannel and channel conformations, as evident from analysis of membrane penetration depth, size-exclusion chromatography, and backbone circular dichroism data. These results are potentially useful in analyzing the effect of tryptophan substitution on the functioning of other ion channels and membrane proteins.

The gramicidin ion channel: A model membrane protein

The linear peptide gramicidin forms prototypical ion channels specific for monovalent cations and has been extensively used to study the organization, dynamics and function of membrane-spanning channels. In recent times, the availability of crystal structures of complex ion channels has challenged the role of gramicidin as a model membrane protein and ion channel. This review focuses on the suitability of gramicidin as a model membrane protein in general, and the information gained from gramicidin to understand lipid-protein interactions in particular. Special emphasis is given to the role and orientation of tryptophan residues in channel structure and function and recent spectroscopic approaches that have highlighted the organization and dynamics of the channel in membrane and membrane-mimetic media.

Bilayer thickness and membrane protein function: an energetic perspective

The lipid bilayer component of biological membranes is important for the distribution, organization, and function of bilayer-spanning proteins. This regulation is due to both specific lipid-protein interactions and general bilayer-protein interactions, which modulate the energetics and kinetics of protein conformational transitions, as well as the protein distribution between different membrane compartments. The bilayer regulation of membrane protein function arises from the hydrophobic coupling between the protein's hydrophobic domains and the bilayer hydrophobic core, which causes protein conformational changes that involve the protein/bilayer boundary to perturb the adjacent bilayer. Such bilayer perturbations, or deformations, incur an energetic cost, which for a given conformational change varies as a function of the bilayer material properties (bilayer thickness, intrinsic lipid curvature, and the elastic compression and bending moduli). Protein function therefore is regulated by changes in bilayer material properties, which determine the free-energy changes caused by the protein-induced bilayer deformation. The lipid bilayer thus becomes an allosteric regulator of membrane function.

Fluorescence quenching of gramicidin D in model membranes by halothane

Inhaled anesthetics were introduced in surgery over a century ago. To this day, the molecular mechanism of anesthetic action remains largely unknown. However, ion-channels of neuronal membranes are believed to be the most- likely molecular targets of inhaled anesthetics. In the study presented here, we investigated the interaction of a simplified ion-channel system, gramicidin, with halothane, a small haloalkane inhaled anesthetic in various environments. Fluorescence-quenching experiments of gramicidin D in dioleoylphosphatidylcholine (DOPC) large unilamellar vesicles (LUVS) have shown that halothane can directly interact with the ion channel (KSV = 66 M–1). Halothane quenched the fluorescence from tryptophan residues located at the lipid bilayer – aqueous interfaces as well as those tryptophans located deeper in the bilayer. Quenching data from gramicidin D in sodium dodecyl sulfide (SDS) micelles revealed that the tryptophan residues located at the micelle–solvent interface were preferentially quenched by halothane (KSV = 22 M–1). In 1-octanol, fluorescence quenching was observed, but with a lower KSV value (KSV = 6 M–1) than in DOPC LUVS and SDS micelles. Taken together, these results indicate that halothane interactions with gramicidin, mediated by a lipid bilayer, are the strongest, and that the mechanism of anesthetic action may also be lipid-mediated.

Biological and structural characterization of new linear gomesin analogues with improved therapeutic indices

Gomesin (Gm) is a potent antimicrobial peptide isolated from the spider Acanthoscurria gomesiana. The two disulfide bridges Cys(2,15) and Cys(6,11) facilitate the folding of the molecule in a beta-hairpin structure, conferring on the peptide a high stability in human plasma. We report herein biological and structural features of new linear Gm analogues, obtained by combining the removal of both disulfide bridges and the incorporation of a D- or L-proline. Regarding their biological properties, two analogues, namely, [D-Thr(2,6,11,15), Pro(9)]-D-Gm and [Thr(2,6,11,15), D-Pro(9)]-Gm, are as potent as Gm against Candida albicans and only fourfold less against Staphylococcus aureus and Escherichia coli. In addition, at 100 microM they are approximately threefold less hemolytic than Gm. The best therapeutic indices were found for [D-Thr(2,6,11,15), Pro(9)]-D-Gm and for [(Des-pGlu(1), -Thr(2), -Arg(3)), Thr(6,11,15), D-Pro(9)]-Gm with a 32-fold increase of their activity against bacteria, and from 128- to 512-fold against yeast when compared with Gm. Regarding the stability, [D-Thr(2,6,11,15), Pro(9)]-D-Gm appeared to be the most resistant in human serum, along with [D-Thr(2,6,11,15), Pro(8)]-D-Gm and [Thr(2,6,11,15), D-Arg(4,16), D-Pro(9)]-Gm. When evaluating their conformation by CD spectroscopy in sodium dodecyl sulfate (SDS), most linear analogues display beta-conformation characteristics. Moreover, considering its high therapeutic index and stability in serum, [D-Thr(2,6,11,15), Pro(9)]-D-Gm was further analyzed by NMR spectroscopy. (1)H NMR experiments in SDS micelles demonstrated that [D-Thr(2,6,11,15), Pro(9)]-D-Gm presents a conformation very similar to that of Gm. In our search for Gm analogues with enhanced potential for drug development, we demonstrated that designing cysteine-free analogues can improve the therapeutic index of Gm derivatives.

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.

Molecular dynamics simulation of the P2Y14 receptor. Ligand docking and identification of a putative binding site of the distal hexose moiety

A rhodopsin-based homology model of the P2Y14 receptor was inserted into a phospholipid bilayer and refined by molecular dynamics (MD) simulation. The binding modes of several known agonists, namely UDP-glucose and its analogues, were proposed using automatic molecular docking combined with Monte Carlo Multiple Minimum calculations. Compared to other P2Y receptors, the P2Y14 receptor has an atypical binding mode of the nucleobase, ribose, and phosphate moieties. The diphosphate moiety interacts with only one cationic residue, namely Lys171 of EL2, while in other P2Y receptor subtypes three Arg or Lys residues interact with the phosphate chain. Two other conserved cationic residues, namely Arg253 (6.55) and Lys277 (7.35) of the P2Y14 receptor together with two anionic residues (Glu166 and Glu174, located in EL2), are likely involved in interactions with the distal hexose moiety.

Structure, topology, and tilt of cell-signaling peptides containing nuclear localization sequences in membrane bilayers determined by solid-state NMR and molecular dynamics simulation studies

Cell-signaling peptides have been extensively used to transport functional molecules across the plasma membrane into living cells. These peptides consist of a hydrophobic sequence and a cationic nuclear localization sequence (NLS). It has been assumed that the hydrophobic region penetrates the hydrophobic lipid bilayer and delivers the NLS inside the cell. To better understand the transport mechanism of these peptides, in this study, we investigated the structure, orientation, tilt of the peptide relative to the bilayer normal, and the membrane interaction of two cell-signaling peptides, SA and SKP. Results from CD and solid-state NMR experiments combined with molecular dynamics simulations suggest that the hydrophobic region is helical and has a transmembrane orientation with the helical axis tilted away from the bilayer normal. The influence of the hydrophobic mismatch, between the hydrophobic length of the peptide and the hydrophobic thickness of the bilayer, on the tilt angle of the peptides was investigated using thicker POPC and thinner DMPC bilayers. NMR experiments showed that the hydrophobic domain of each peptide has a tilt angle of 15 +/- 3 degrees in POPC, whereas in DMPC, 25 +/- 3 degree and 30 +/- 3 degree tilts were observed for SA and SKP peptides, respectively. These results are in good agreement with molecular dynamics simulations, which predict a tilt angle of 13.3 degrees (SA in POPC), 16.4 degrees (SKP in POPC), 22.3 degrees (SA in DMPC), and 31.7 degrees (SKP in DMPC). These results and simulations on the hydrophobic fragment of SA or SKP suggest that the tilt of helices increases with a decrease in bilayer thickness without changing the phase, order, and structure of the lipid bilayers.

The open state gating mechanism of gramicidin a requires relative opposed monomer rotation and simultaneous lateral displacement

The gating mechanism of the open state of the gramicidin A (gA) channel is studied by using a new Monte Carlo Normal Mode Following (MC-NMF) technique, one applicable even without a target structure. The results demonstrate that the lowest-frequency normal mode (NM) at approximately 6.5 cm(-1) is the crucial mode that initiates dissociation. Perturbing the gA dimer in either direction along this NM leads to opposed, nearly rigid-body rotations of the gA monomers around the central pore axis. Tracking this NM by using the eigenvector-following technique reveals the channel's gating mechanism: dissociation via relative opposed monomer rotation and simultaneous lateral displacement. System evolution along the lowest-frequency eigenvector shows that the large-amplitude motions required for gating (dissociation) are not simple relative rigid-body motions of the monomers. Gating involves coupling intermonomer hydrogen bond breaking, backbone realignment, and relative monomer tilt with complex side chain reorganization at the intermonomer junction.

Membrane protein simulations with a united-atom lipid and all-atom protein model: lipid-protein interactions, side chain transfer free energies and model proteins

We have reparameterized the dihedral parameters in a commonly used united-atom lipid force field so that they can be used with the all-atom OPLS force field for proteins implemented in the molecular dynamics simulation software GROMACS. Simulations with this new combination give stable trajectories and sensible behaviour of both lipids and protein. We have calculated the free energy of transfer of amino acid side chains between water and 'lipid-cyclohexane', made of lipid force field methylene groups, as a hydrophobic mimic of the membrane interior, for both the OPLS-AA and a modified OPLS-AA force field which gives better hydration free energies under simulation conditions close to those preferred for the lipid force field. The average error is 4.3 kJ mol(-1) for water-'lipid-cyclohexane' compared to 3.2 kJ mol(-1) for OPLS-AA cyclohexane and 2.4 kJ mol(-1) for the modified OPLS-AA water-'lipid-cyclohexane'. We have also investigated the effect of different methods to combine parameters between the united-atom lipid force field and the united-atom protein force field ffgmx. In a widely used combination, the strength of interactions between hydrocarbon lipid tails and proteins is significantly overestimated, causing a decrease in the area per lipid and an increase in lipid ordering. Using straight combination rules improves the results. Combined, we suggest that using OPLS-AA together with the united-atom lipid force field implemented in GROMACS is a reasonable approach to membrane protein simulations. We also suggest that using partial volume information and free energies of transfer may help to improve the parameterization of lipid-protein interactions and point out the need for accurate experimental data to validate and improve force field descriptions of such interactions.

Regulation of membrane protein function by lipid bilayer elasticity-a single molecule technology to measure the bilayer properties experienced by an embedded protein

Membrane protein function is generally regulated by the molecular composition of the host lipid bilayer. The underlying mechanisms have long remained enigmatic. Some cases involve specific molecular interactions, but very often lipids and other amphiphiles, which are adsorbed to lipid bilayers, regulate a number of structurally unrelated proteins in an apparently non-specific manner. It is well known that changes in the physical properties of a lipid bilayer (e.g., thickness or monolayer spontaneous curvature) can affect the function of an embedded protein. However, the role of such changes, in the general regulation of membrane protein function, is unclear. This is to a large extent due to lack of a generally accepted framework in which to understand the many observations. The present review summarizes studies which have demonstrated that the hydrophobic interactions between a membrane protein and the host lipid bilayer provide an energetic coupling, whereby protein function can be regulated by the bilayer elasticity. The feasibility of this 'hydrophobic coupling mechanism' has been demonstrated using the gramicidin channel, a model membrane protein, in planar lipid bilayers. Using voltage-dependent sodium channels, N-type calcium channels and GABA(A) receptors, it has been shown that membrane protein function in living cells can be regulated by amphiphile induced changes in bilayer elasticity. Using the gramicidin channel as a molecular force transducer, a nanotechnology to measure the elastic properties experienced by an embedded protein has been developed. A theoretical and technological framework, to study the regulation of membrane protein function by lipid bilayer elasticity, has been established.

Molecular dynamics simulations of gramicidin A in a lipid bilayer: From structure–function relations to force fields

Molecular dynamics simulations of membrane proteins have become a popular tool for studying their dynamic features, which are not easily accessible by experiments. Whether the force fields developed for globular proteins are adequate this purpose is an important question that is often glossed over. Here we determine the permeation properties of potassium ions in the gramicidin A channel in a lipid bilayer from free energy simulations, and compare the results to experimental data. In particular, we check the dependence of the free energy barriers ions face at the channel center on the membrane size. The results indicate that there is a serious problem with the current rigid force fields independent of the membrane size, and new, possibly polarizable, force fields need to be developed to resolve this problem.

Ion permeation through a narrow channel: using gramicidin to ascertain all-atom molecular dynamics potential of mean force methodology and biomolecular force fields

We investigate methods for extracting the potential of mean force (PMF) governing ion permeation from molecular dynamics simulations (MD) using gramicidin A as a prototypical narrow ion channel. It is possible to obtain well-converged meaningful PMFs using all-atom MD, which predict experimental observables within order-of-magnitude agreement with experimental results. This was possible by careful attention to issues of statistical convergence of the PMF, finite size effects, and lipid hydrocarbon chain polarizability. When comparing the modern all-atom force fields of CHARMM27 and AMBER94, we found that a fairly consistent picture emerges, and that both AMBER94 and CHARMM27 predict observables that are in semiquantitative agreement with both the experimental conductance and dissociation coefficient. Even small changes in the force field, however, result in significant changes in permeation energetics. Furthermore, the full two-dimensional free-energy surface describing permeation reveals the location and magnitude of the central barrier and the location of two binding sites for K(+) ion permeation near the channel entrance--i.e., an inner site on-axis and an outer site off-axis. We conclude that the MD-PMF approach is a powerful tool for understanding and predicting the function of narrow ion channels in a manner that is consistent with the atomic and thermally fluctuating nature of proteins.

Molecular dynamics - potential of mean force calculations as a tool for understanding ion permeation and selectivity in narrow channels

Ion channels catalyze the permeation of charged molecules across cell membranes and are essential for many vital physiological functions, including nerve and muscle activity. To understand better the mechanisms underlying ion conduction and valence selectivity of narrow ion channels, we have employed free energy techniques to calculate the potential of mean force (PMF) for ion movement through the prototypical gramicidin A channel. Employing modern all-atom molecular dynamics (MD) force fields with umbrella sampling methods that incorporate one hundred 1-2 ns trajectories, we find that it is possible to achieve semi-quantitative agreement with experimental binding and conductance measurements. We also examine the sensitivity of the MD-PMF results to the choice of MD force field and compare PMFs for potassium, calcium and chloride ions to explore the basis for the valence selectivity of this narrow and uncharged ion channel. A large central barrier is observed for both anions and divalent ions, consistent with lack of experimental conductance. Neither anion or divalent cation is seen to be stabilized inside the channel relative to the bulk electrolyte and each leads to large disruptions to the protein and membrane structure when held deep inside the channel. Weak binding of calcium ions outside the channel corresponds to a free energy well that is too shallow to demonstrate channel blocking. Our findings emphasize the success of the MD-PMF approach and the sensitivity of ion energetics to the choice of biomolecular force field.

Conformational sampling and dynamics of membrane proteins from 10-nanosecond computer simulations

In the current report, we provide a quantitative analysis of the convergence of the sampling of conformational space accomplished in molecular dynamics simulations of membrane proteins of duration in the order of 10 nanoseconds. A set of proteins of diverse size and topology is considered, ranging from helical pores such as gramicidin and small beta-barrels such as OmpT, to larger and more complex structures such as rhodopsin and FepA. Principal component analysis of the C(alpha)-atom trajectories was employed to assess the convergence of the conformational sampling in both the transmembrane domains and the whole proteins, while the time-dependence of the average structure was analyzed to obtain single-domain information. The membrane-embedded regions, particularly those of small or structurally simple proteins, were found to achieve reasonable convergence. By contrast, extra-membranous domains lacking secondary structure are often markedly under-sampled, exhibiting a continuous structural drift. This drift results in a significant imprecision in the calculated B-factors, which detracts from any quantitative comparison to experimental data. In view of such limitations, we suggest that similar analyses may be valuable in simulation studies of membrane protein dynamics, in order to attach a level of confidence to any biologically relevant observations.

Importance of tensor asymmetry for the analysis of 2H NMR spectra from deuterated aromatic rings

We have used ab initio calculations to compute all of the tensor elements of the electric field gradient for each carbon-deuterium bond in the ring of deuterated 3-methyl-indole. Previous analyses have ignored the smaller tensor elements perpendicular to principal component Vzz which is aligned with the C-2H bond (local bond z-axis). At each ring position, the smallest element Vxx is in the molecular plane and Vyy is normal to the plane of the ring. The asymmetry parameter = (Vyy - Vxx)/Vzz ranges from 0.07 at C4 to 0.11 at C2. We used the perpendicular (off-bond) tensor elements, in concert with an improved understanding of the indole ring geometry, to analyze prototype 2H NMR spectra from well-oriented, hydrated peptide/lipid samples. For each of the four tryptophans of membrane-spanning gramicidin A (gA) channels, the inclusion of the perpendicular elements changes the deduced ring tilt by nearly 10 and increases the ring principal order parameter Szz for overall "wobble" with respect to the membrane normal (molecular z-axis). With the improved analysis, the magnitude of Szz for the outermost indole rings of Trp13 and Trp15 is indistinguishable from that observed previously for backbone atoms (0.93 +/- 0.03). For the Trp9 and Trp11 rings, which are slightly more buried within the membrane, Szz is slightly lower (0.86 +/- 0.03). The results show that the perpendicular elements are important for the detailed analysis of 2H NMR spectra from aromatic ring systems.

The orientation effect of gramicidin A on bicelles and Eu3+-doped bicelles as studied by solid-state NMR and FT-IR spectroscopy

We have explored the effect of gramicidin A (gA) on bicelle (Bic) orientation in the absence and presence of Eu(3+) by (31)P and (2)H NMR at different DMPC/gA ratios. FT-IR spectroscopy was used to assess the lipid chain ordering and verify the transmembrane peptide conformation. Our results show a time-dependent flipping of the bilayer normal alignment at high temperatures and high proportion of gA. The results are explained by both the diamagnetic susceptibility anisotropy of the beta(6.3) helical peptides and viscosity of the lipid mixture. The concentration effect of gramicidin on Bic/Eu(3+) is compared to that on Eu(3+)-doped DMPC liposomes. The Bic/Eu(3+) system is no longer oriented in the presence of gA and adopts a vesicular morphology while the peptide incorporation induces the formation of ellipsoidal DMPC/Eu(3+) assemblies aligned with their normal parallel to the magnetic field. The difference is explained in terms of lipid chain disorder and size of the bilayers.

Role of protein flexibility in ion permeation: a case study in gramicidin A

Proteins have a flexible structure, and their atoms exhibit considerable fluctuations under normal operating conditions. However, apart from some enzyme reactions involving ligand binding, our understanding of the role of flexibility in protein function remains mostly incomplete. Here we investigate this question in the realm of membrane proteins that form ion channels. Specifically, we consider ion permeation in the gramicidin A channel, and study how the energetics of ion conduction changes as the channel structure is progressively changed from completely flexible to a fixed one. For each channel structure, the potential of mean force for a permeating potassium ion is determined from molecular dynamics (MD) simulations. Using the same molecular dynamics data for completely flexible gramicidin A, we also calculate the average densities and fluctuations of the peptide atoms and investigate the correlations between these fluctuations and the motion of a permeating ion. Our results show conclusively that peptide flexibility plays an important role in ion permeation in the gramicidin A channel, thus providing another reason--besides the well-known problem with the description of single file pore water--why this channel cannot be modeled using continuum electrostatics with a fixed structure. The new method developed here for studying the role of protein flexibility on its function clarifies the contributions of the fluctuations to energy and entropy, and places limits on the level of detail required in a coarse-grained model.

Interfacial tryptophan residues: a role for the cation-pi effect?

Integral membrane proteins are characterized by having a preference for aromatic residues, e.g., tryptophan (W), at the interface between the lipid bilayer core and the aqueous phase. The reason for this is not clear, but it seems that the preference is related to a complex interplay between steric and electrostatic forces. The flat rigid paddle-like structure of tryptophan, associated with a quadrupolar moment (aromaticity) arising from the pi-electron cloud of the indole, interacts primarily with moieties in the lipid headgroup region hardly penetrating into the bilayer core. We have studied the interaction between the nitrogen moiety of lipid molecule headgroups and the pi-electron distribution of gramicidin (gA) tryptophan residues (W9, W11, W13, and W15) using molecular dynamics (MD) simulations of gA embedded in two hydrated lipid bilayers composed of 1-palmitoyl-2-oleoylphosphatidylethanolamine (POPE) and 1-palmitoyl-2-oleoylphosphatidyl-choline (POPC), respectively. We use a force field model for tryptophan in which polarizability is only implicit, but we believe that classical molecular dynamics force fields are sufficient to capture the most prominent features of the cation-pi interaction. Our criteria for cation-pi interactions are based on distance and angular requirements, and the results from our model suggest that cation-pi interactions are relevant for W(PE)1), W(PE)13, W(PE)15, and, to some extent, W(PC)11 and W(PC)13. In our model, W9 does not seem to engage in cation-pi interactions with lipids, neither in POPE nor POPC. The criteria for the cation-pi effect are satisfied more often in POPE than in POPC, whereas the H-bonding ability between the indole donor and the carbonyl acceptor is similar in POPE and POPC. This suggests an increased affinity for lipids with ethanolamine headgroups to transmembrane proteins enriched in interfacial tryptophans.

Molecular dynamics simulations of the influenza hemagglutinin fusion peptide in micelles and bilayers: conformational analysis of peptide and lipids

Molecular dynamics simulations of the influenza hemagglutinin fusion peptide in two differently sized dodecylphosphocholine micelles and a palmitoyl oleoyl phosphatidylcholine bilayer were generated to analyze the influence of the environment. Four independent trajectories (5 ns each for the bilayer, and 2 ns each for the micelles) were generated for each system. The peptide lies at the surface of the micelles, while its N-terminal region inserts deeply in the bilayer. This leads to a substantial increase of the solvation and rigidity of the peptide in micelles as compared to the bilayer. The average structures, nevertheless, are similar in all three systems and agree reasonably with micelle-based NMR structures. When in the bilayer, the peptide increases the chain gauche population and area of adjacent lipids in the same binding leaflet, while it has the opposite effect for the nearby lipids of the other leaflet. These changes, which occur spontaneously to fill voids and defects, cause a decrease in the thickness of the membrane in the neighborhood of the peptide. They would be expected to promote positive curvature, as consistent with the formation of the convex bulge, or "nipple", in the initial stage of membrane fusion. An extension of the classical surfactant theory of Israelachvili based on shapes is proposed to introduce the concept of a "dynamically induced shape" of the membrane lipids by the peptide.

Effect of structural transition of the host assembly on dynamics of an ion channel peptide: a fluorescence approach

Structural transition can be induced in charged micelles by increasing the ionic strength of the medium. We have monitored the organization and dynamics of the functionally important tryptophan residues of gramicidin in spherical and rod-shaped sodium dodecyl sulfate micelles utilizing a combination of wavelength-selective fluorescence and related fluorescence approaches. Our results show that tryptophans in gramicidin, present in the single-stranded beta(6.3) conformation, experience slow solvent relaxation giving rise to red edge excitation shift in spherical and rod-shaped micelles. In addition, changes in fluorescence polarization with increasing excitation or emission wavelength reinforce that the gramicidin tryptophans are localized in motionally restricted regions of these micelles. Fluorescence quenching experiments using acrylamide as a quencher of tryptophan fluorescence show that there is reduced water penetration in rod-shaped micelles. Taken together, we show that gramicidin conformation and dynamics is sensitive to the salt-induced structural transition in charged micelles. In addition, these results demonstrate that deformation of the host assembly could modulate protein conformation and dynamics.

Gramicidin channels

Gramicidin channels are mini-proteins composed of two tryptophan-rich subunits. The conducting channels are formed by the transbilayer dimerization of nonconducting subunits, which are tied to the bilayer/solution interface through hydrogen bonds between the indole NH groups and the phospholipid backbone and water. The channel structure is known at atomic resolution and the channel's permeability characteristics are particularly well defined: gramicidin channels are selective for monovalent cations, with no measurable permeability to anions or polyvalent cations; ions and water move through a pore whose wall is formed by the peptide backbone; and the single-channel conductance and cation selectivity vary when the amino acid sequence is varied, even though the permeating ions make no contact with the amino acid side chains. Given the amount of experimental information that is available--for both the wild-type channels and for channels formed by amino acid-substituted gramicidin analogues--gramicidin channels provide important insights into the microphysics of ion permeation through bilayer-spanning channels. For the same reason, gramicidin channels constitute the system of choice for evaluating computational strategies for obtaining mechanistic insights into ion permeation through the complex channels formed by integral membrane proteins.

On the importance of atomic fluctuations, protein flexibility, and solvent in ion permeation

Proteins, including ion channels, often are described in terms of some average structure and pictured as rigid entities immersed in a featureless solvent continuum. This simplified view, which provides for a convenient representation of the protein's overall structure, incurs the risk of deemphasizing important features underlying protein function, such as thermal fluctuations in the atom positions and the discreteness of the solvent molecules. These factors become particularly important in the case of ion movement through narrow pores, where the magnitude of the thermal fluctuations may be comparable to the ion pore atom separations, such that the strength of the ion channel interactions may vary dramatically as a function of the instantaneous configuration of the ion and the surrounding protein and pore water. Descriptions of ion permeation through narrow pores, which employ static protein structures and a macroscopic continuum dielectric solvent, thus face fundamental difficulties. We illustrate this using simple model calculations based on the gramicidin A and KcsA potassium channels, which show that thermal atomic fluctuations lead to energy profiles that vary by tens of kcal/mol. Consequently, within the framework of a rigid pore model, ion-channel energetics is extremely sensitive to the choice of experimental structure and how the space-dependent dielectric constant is assigned. Given these observations, the significance of any description based on a rigid structure appears limited. Creating a conducting channel model from one single structure requires substantial and arbitrary engineering of the model parameters, making it difficult for such approaches to contribute to our understanding of ion permeation at a microscopic level.

Role of ions on structure and stability of a synthetic gramicidin ion channel in solution. A molecular dynamics study

We performed a molecular dynamics (MD) simulation to the investigate structure and stability of a synthetic gramicidin-like peptide in solution with and without ions. The starting structures of the MD simulations were taken from two recently solved NMR structures of this peptide in isotropic solution, which forms stable monomers or dimers in the presence or absence of ions, respectively. The monomeric structure is channel-like and is assumed to be stabilized by the presence of two Cs(+) ions bound in the channel, each one close to one channel entrance. In our MD simulations, we observed how the Cs(+) ions bind in the channel formed by the monomeric gramicidin-like peptide using implicit solvent and explicit ions with a concentration of 2 M. MD simulations were performed with and without explicit ions but with an implicit solvent model defined by the generalized Born approximation, which was used to mimic the dielectric properties of the solvent and to speed up the computations.

Computing numerically the access resistance of a pore

The access resistance (AR) of a channel is an important component of the conductance of ion channels, particularly in wide and short channels, where it accounts for a substantial fraction of the total resistance to the movement of ions. The AR is usually calculated by using a classical and simple expression derived by Hall from electrostatics (J.E. Hall 1975 J. Gen. Phys. 66:531-532), though other expressions, both analytical and numerical, have been proposed. Here we report some numerical results for the AR of a channel obtained by solving the Poisson-Nernst-Planck equations at the entrance of a circular pore. Agreement is found between numerical calculations and analytical results from Hall's equation for uncharged pores in neutral membranes. However, for channels embedded in charged membranes, Hall's expression overestimates the AR, which is much lower and can even be neglected in some cases. The weak dependence of AR on the pore radius for charged membranes at low salt concentration can be exploited to separate the channel and the access contributions to the measured conductance.

Test of molecular dynamics force fields in gramicidin A

The force fields commonly used in molecular dynamics simulations of proteins are optimized under bulk conditions. Whether the same force fields can be used in simulations of membrane proteins is not well established, although they are increasingly being used for such purposes. Here we consider ion permeation in the gramicidin A channel as a test of the AMBER force field in a membrane environment. The potentials of mean force for potassium ions are calculated along the channel axis and compared with the one deduced from the experimental conductance data. The calculated result indicates a rather large central barrier similar to those obtained from other force fields, which are incompatible with the conductance data. We suggest that lack of polarizability is the most likely cause of this problem, and, therefore, urge development of polarizable force fields for simulations of membrane proteins.

Ionic permeation free energy in gramicidin: a semimicroscopic perspective

Ion permeation through the gramicidin channel is studied using a model that circumvents two major difficulties inherent to standard simulational methods. It exploits the timescale separation between electronic and structural contributions to dielectric stabilization, accounting for the influence of electronic polarization by embedding the channel in a dielectric milieu that describes this polarization in a mean sense. The explicit mobile moieties are the ion, multipolar waters, and the carbonyls and amides of the peptide backbone. The model treats the influence of aromatic residues and the membrane dipole potential. A new electrical geometry is introduced that treats long-range electrostatics exactly and avoids problems related to periodic boundary conditions. It permits the translocating ion to make a seamless transition from nearby electrolyte to the channel interior. Other degrees of freedom (more distant bulk electrolyte and nonpolar lipid) are treated as dielectric continua. Reasonable permeation free energy profiles are obtained for potassium, rubidium, and cesium; binding wells are shallow and the central barrier is small. Estimated cationic single-channel conductances are smaller than experiment, but only by factors between 2 (rubidium) and 50 (potassium). When applied to chloride the internal barrier is large, with a corresponding miniscule single-channel conductance. The estimated relative single-channel conductances of gramicidin A, B, and C agree well with experiment.

Gramicidin-perforated patch recording revealed the oscillatory nature of secretory Cl- movements in salivary acinar cells

Elevations of cytoplasmic free calcium concentrations ([Ca²⁺]_i) evoked by cholinergic agonists stimulate isotonic fluid secretion in salivary acinar cells. This process is driven by the apical exit of Cl⁻ through Ca²⁺-activated Cl⁻ channels, while Cl⁻ enters the cytoplasm against its electrochemical gradient via a loop diuretic-sensitive Na⁺-K⁺-2Cl⁻ cotransporter (NKCC) and/or parallel operations of Cl⁻-HCO₃⁻ and Na⁺-H⁺ exchangers, located in the basolateral membrane. To characterize the contributions of those activities to net Cl⁻ secretion, we analyzed carbachol (CCh)-activated Cl⁻ currents in submandibular acinar cells using the “gramicidin-perforated patch recording configuration.” Since the linear polypeptide antibiotic gramicidin creates monovalent cation-selective pores, CCh-activated Cl⁻ currents in the gramicidin-perforated patch recording were carried by Cl⁻ efflux via Cl⁻ channels, dependent upon Cl⁻ entry through Cl⁻ transporters expressed in the acinar cells. CCh-evoked oscillatory Cl⁻ currents were associated with oscillations of membrane potential. Bumetanide, a loop diuretic, decreased the CCh-activated Cl⁻ currents and hyperpolarized the membrane potential. In contrast, neither methazolamide, a carbonic anhydrase inhibitor, nor elimination of external HCO₃⁻ had significant effects, suggesting that the cotransporter rather than parallel operations of Cl⁻-HCO₃⁻ and Na⁺-H⁺ exchangers is the primary Cl⁻ uptake pathway. Pharmacological manipulation of the activities of the Ca²⁺-activated Cl⁻ channel and the NKCC revealed that the NKCC plays a substantial role in determining the amplitude of oscillatory Cl⁻ currents, while adjusting to the rate imposed by the Ca²⁺-activated Cl⁻ channel, in the gramicidin-perforated patch configuration. By concerting with and being controlled by the cation steps, the oscillatory form of secretory Cl⁻ movements may effectively provide a driving force for fluid secretion in intact acinar cells.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

Biosynthesis of nonribosomal peptides1

Bacteria and fungi use large multifunctional enzymes, the so-called nonribosomal peptide synthetases (NRPSs), to produce peptides of broad structural and biological activity. Biochemical studies have contributed substantially to the understanding of the key principles of these modular enzymes that can draw on a much larger number of catalytic tools for the incorporation of unusual features compared with the ribosomal system. Several crystal structures of NRPS-domains have yielded deep insight into the catalytic mechanisms involved and have led to a better prediction of the products assembled and to the construction of hybrid enzymes. In addition to the structure-function relationship of the core- and tailoring-domains of NRPSs, which is the main focus of this review, different biosynthetic strategies and essential enzymes for posttranslational modification and editing are discussed.

Regulation of sodium channel function by bilayer elasticity: the importance of hydrophobic coupling. Effects of Micelle-forming amphiphiles and cholesterol

Membrane proteins are regulated by the lipid bilayer composition. Specific lipid–protein interactions rarely are involved, which suggests that the regulation is due to changes in some general bilayer property (or properties). The hydrophobic coupling between a membrane-spanning protein and the surrounding bilayer means that protein conformational changes may be associated with a reversible, local bilayer deformation. Lipid bilayers are elastic bodies, and the energetic cost of the bilayer deformation contributes to the total energetic cost of the protein conformational change. The energetics and kinetics of the protein conformational changes therefore will be regulated by the bilayer elasticity, which is determined by the lipid composition. This hydrophobic coupling mechanism has been studied extensively in gramicidin channels, where the channel–bilayer hydrophobic interactions link a “conformational” change (the monomer↔dimer transition) to an elastic bilayer deformation. Gramicidin channels thus are regulated by the lipid bilayer elastic properties (thickness, monolayer equilibrium curvature, and compression and bending moduli). To investigate whether this hydrophobic coupling mechanism could be a general mechanism regulating membrane protein function, we examined whether voltage-dependent skeletal-muscle sodium channels, expressed in HEK293 cells, are regulated by bilayer elasticity, as monitored using gramicidin A (gA) channels. Nonphysiological amphiphiles (β-octyl-glucoside, Genapol X-100, Triton X-100, and reduced Triton X-100) that make lipid bilayers less “stiff”, as measured using gA channels, shift the voltage dependence of sodium channel inactivation toward more hyperpolarized potentials. At low amphiphile concentration, the magnitude of the shift is linearly correlated to the change in gA channel lifetime. Cholesterol-depletion, which also reduces bilayer stiffness, causes a similar shift in sodium channel inactivation. These results provide strong support for the notion that bilayer–protein hydrophobic coupling allows the bilayer elastic properties to regulate membrane protein function.

Homology modeling and molecular dynamics simulations of transmembrane domain structure of human neuronal nicotinic acetylcholine receptor

A three-dimensional model of the transmembrane domain of a neuronal-type nicotinic acetylcholine receptor (nAChR), (alpha4)2(beta2)3, was constructed from a homology structure of the muscle-type nAChR recently determined by cryo-electron microscopy. The neuronal channel model was embedded in a fully hydrated DMPC lipid bilayer, and molecular-dynamics simulations were performed for 5 ns. A comparative analysis of the neuronal- versus muscle-type nAChR models revealed many conserved pore-lining residues, but an important difference was found near the periplasmic mouth of the pore. A flickering salt-bridge of alpha4-E266 with its adjacent beta2-K260 was observed in the neuronal-type channel during the course of the molecular-dynamics simulations. The narrowest region, with a pore radius of approximately 2 A formed by the salt-bridges, does not seem to be the restriction site for a continuous water passage. Instead, two hydrophobic rings, formed by alpha4-V259, alpha4-L263, and the homologous residues in the beta2-subunits, act as the gates for water flow, even though the region has a slightly larger pore radius. The model offers new insight into the water transport across the (alpha4)2(beta2)3 nAChR channel, and may lead to a better understanding of the structures, dynamics, and functions of this family of ion channels.

Effect of graded hydration on the dynamics of an ion channel peptide: a fluorescence approach

Water plays an important role in determining the folding, structure, dynamics, and, in turn, the function of proteins. We have utilized a combination of fluorescence approaches such as the wavelength-selective fluorescence approach to monitor the effect of varying degrees of hydration on the organization and dynamics of the functionally important tryptophan residues of gramicidin in reverse micelles formed by sodium bis(2-ethylhexyl) sulfosuccinate. Our results show that tryptophans in gramicidin, present in the single-stranded beta6.3 conformation, experience slow solvent relaxation giving rise to red-edge excitation shift (REES). In addition, changes in fluorescence polarization with increasing excitation or emission wavelength reinforce that the gramicidin tryptophans are localized in motionally restricted regions of the reverse micelle. Interestingly, the extent of REES is found to be independent of the [water]/[surfactant] molar ratio (w(o)). We attribute this to heterogeneity in gramicidin tryptophan localization. Fluorescence intensity and mean fluorescence lifetime of the gramicidin tryptophans show significant reductions with increasing w(o) indicating sensitivity to increased polarity. Since the dynamics of hydration is related to folding, structure, and eventually function of proteins, we conclude that REES could prove to be a potentially sensitive tool to explore the dynamics of proteins under conditions of changing hydration.

Influence of protein flexibility on the electrostatic energy landscape in gramicidin A

We describe an electrostatic model of the gramicidin A channel that allows protein atoms to move in response to the presence of a permeating ion. To do this, molecular dynamics simulations are carried out with a permeating ion at various positions within the channel. Then an ensemble of atomic coordinates taken from the simulations are used to construct energy profiles using macroscopic electrostatic calculations. The energy profiles constructed are compared to experimentally-determined conductance data by inserting them into Brownian dynamics simulations. We find that the energy landscape seen by a permeating ion changes significantly when we allow the protein atoms to move rather than using a rigid protein structure. However, the model developed cannot satisfactorily reproduce all of the experimental data. Thus, even when protein atoms are allowed to move, the dielectric model used in our electrostatic calculations breaks down when modeling the gramicidin channel.

Lipid-protein interactions of integral membrane proteins: a comparative simulation study

The interactions between membrane proteins and their lipid bilayer environment play important roles in the stability and function of such proteins. Extended (15-20 ns) molecular dynamics simulations have been used to explore the interactions of two membrane proteins with phosphatidylcholine bilayers. One protein (KcsA) is an alpha-helix bundle and embedded in a palmitoyl oleoyl phosphatidylcholine bilayer; the other (OmpA) is a beta-barrel outer-membrane protein and is in a dimyristoyl phosphatidylcholine bilayer. The simulations enable analysis in detail of a number of aspects of lipid-protein interactions. In particular, the interactions of aromatic amphipathic side chains (i.e., Trp, Tyr) with lipid headgroups, and "snorkeling" interactions of basic side chains (i.e., Lys, Arg) with phosphate groups are explored. Analysis of the number of contacts and of H-bonds reveal fluctuations on an approximately 1- to 5-ns timescale. There are two clear bands of interacting residues on the surface of KcsA, whereas there are three such bands on OmpA. A large number of Arg-phosphate interactions are seen for KcsA; for OmpA, the number of basic-phosphate interactions is smaller and shows more marked fluctuations with respect to time. Both classes of interaction occur in clearly defined interfacial regions of width approximately 1 nm. Analysis of lateral diffusion of lipid molecules reveals that "boundary" lipid molecules diffuse at about half the rate of bulk lipid. Overall, these simulations present a dynamic picture of lipid-protein interactions: there are a number of more specific interactions but even these fluctuate on an approximately 1- to 5-ns timescale.

De novo folding of membrane proteins: an exploration of the structure and NMR properties of the fd coat protein

De novo folding simulations of the major pVIII coat protein from filamentous fd bacteriophage, using a newly developed implicit membrane generalized Born model and replica-exchange molecular dynamics, are presented and discussed. The quality of the predicted structures, judged by comparison of the root-mean-square deviations of a room temperature ensemble of conformations from the replica-exchange simulations and experimental structures from both solid-state NMR in lipid bilayers and solution-phase NMR on the protein in micelles, was quite good, reinforcing the general quality of the folding simulations. The transmembrane helical segment of the protein was well defined in comparison with experiment and the amphipathic helical fragment remained at the membrane/aqueous phase boundary while undergoing significant conformational flexibility due to the loop connecting the two helical segments of the protein. Additional comparisons of computed solid-state NMR properties, the 15N chemical shift and 15N-1H dipolar coupling constants, showed semi-quantitative agreement with the corresponding measurements. These findings suggest an emerging potential for the de novo investigation of integral membrane peptides and proteins and a mechanism to assist experimental approaches to the characterization and structure determination of these important systems.

Energetics of ion conduction through the gramicidin channel

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.