Implicit solvent model studies of the interactions of the influenza hemagglutinin fusion peptide with lipid bilayers

The "fusion peptide," a segment of approximately 20 residues of the influenza hemagglutinin (HA), is necessary and sufficient for HA-induced membrane fusion. We used mean-field calculations of the free energy of peptide-membrane association (DeltaG(tot)) to deduce the most probable orientation of the fusion peptide in the membrane. The main contributions to DeltaG(tot) are probably from the electrostatic (DeltaG(el)) and nonpolar (DeltaG(np)) components of the solvation free energy; these were calculated using continuum solvent models. The peptide was described in atomic detail and was modeled as an alpha-helix based on spectroscopic data. The membrane's hydrocarbon region was described as a structureless slab of nonpolar medium embedded in water. All the helix-membrane configurations, which were lower in DeltaG(tot) than the isolated helix in the aqueous phase, were in the same (wide) basin in configurational space. In each, the helix was horizontally adsorbed at the water-bilayer interface with its principal axis parallel to the membrane plane, its hydrophobic face dissolved in the bilayer, and its polar face in the water. The associated DeltaG(tot) value was approximately -8 to -10 kcal/mol (depending on the rotameric state of one of the phenylalanine residues). In contrast, the DeltaG(tot) values associated with experimentally observed oblique orientations were found to be near zero, suggesting they are marginally stable at best. The theoretical model did not take into account the interactions of the polar headgroups with the peptide and peptide-induced membrane deformation effects. Either or both may overcompensate for the DeltaG(tot) difference between the horizontal and oblique orientations.     

Sequential unfolding of the hemolysin two‐partner secretion domain from Proteus mirabilis

Protein secretion is a major contributor to Gram‐negative bacterial virulence. Type Vb or two‐partner secretion (TPS) pathways utilize a membrane bound β‐barrel B component (TpsB) to translocate large and predominantly virulent exoproteins (TpsA) through a nucleotide independent mechanism. We focused our studies on a truncated TpsA member termed hemolysin A (HpmA265), a structurally and functionally characterized TPS domain from Proteus mirabilis. Contrary to the expectation that the TPS domain of HpmA265 would denature in a single cooperative transition, we found that the unfolding follows a sequential model with three distinct transitions linking four states. The solvent inaccessible core of HpmA265 can be divided into two different regions. The C‐proximal region contains nonpolar residues and forms a prototypical hydrophobic core as found in globular proteins. The N‐proximal region of the solvent inaccessible core, however, contains polar residues. To understand the contributions of the hydrophobic and polar interiors to overall TPS domain stability, we conducted unfolding studies on HpmA265 and site‐specific mutants of HpmA265. By correlating the effect of individual site‐specific mutations with the sequential unfolding results we were able to divide the HpmA265 TPS domain into polar core, nonpolar core, and C‐terminal subdomains. Moreover, the unfolding studies provide quantitative evidence that the folding free energy for the polar core subdomain is more favorable than for the nonpolar core and C‐terminal subdomains. This study implicates the hydrogen bonds shared among these conserved internal residues as a primary means for stabilizing the N‐proximal polar core subdomain.     

Structural basis for the unusual specificity of Escherichia coli aminopeptidase N

Aminopeptidase N from Escherichia coli is a M1 class aminopeptidase with the active-site region related to that of thermolysin. The enzyme has unusual specificity, cleaving adjacent to the large, nonpolar amino acids Phe and Tyr but also cleaving next to the polar residues Lys and Arg. To try to understand the structural basis for this pattern of hydrolysis, the structure of the enzyme was determined in complex with the amino acids L-arginine, L-lysine, L-phenylalanine, L-tryptophan, and L-tyrosine. These amino acids all bind with their backbone atoms close to the active-site zinc ion and their side chain occupying the S1 subsite. This subsite is in the form of a cylinder, about 10 A in cross-section and 12 A in length. The bottom of the cylinder includes the zinc ion and a number of polar side chains that make multiple hydrogen-bonding and other interactions with the alpha-amino group and the alpha-carboxylate of the bound amino acid. The walls of the S1 cylinder are hydrophobic and accommodate the nonpolar or largely nonpolar side chains of Phe and Tyr. The top of the cylinder is polar in character and includes bound water molecules. The epsilon-amino group of the bound lysine side chain and the guanidinium group of arginine both make multiple hydrogen bonds to this part of the S1 site. At the same time, the hydrocarbon part of the lysine and arginine side chains is accommodated within the nonpolar walls of the S1 cylinder. This combination of hydrophobic and hydrophilic binding surfaces explains the ability of ePepN to cleave Lys, Arg, Phe, and Tyr. Another favored substrate has Ala at the P1 position. The short, nonpolar side chain of this residue can clearly be bound within the hydrophobic part of the S1 cylinder, but the reason for its facile hydrolysis remains uncertain.     

Modeling of peptides and proteins in membrane environment. I. A solvation model mimicking a lipid bilayer

A theoretical solvation model of peptides and proteins that mimics the heterogeneous membrane-water system was proposed. Our approach is based on the combined use of atomic parameters of solvation for water and hydrocarbons, which approximates the hydrated polar groups and acyl chains of lipids, respectively. This model was tested in simulations of several peptides: a nonpolar 20-mer polyleucine, a hydrophobic peptide with terminal polar groups, and a strongly amphiphilic peptide. The conformational space of the peptides in the presence of the membrane was studied by the Monte Carlo method. Unlike a polar solvent and vacuum, the membrane-like environment was shown to stabilize the alpha-helical conformation: low-energy structures have a helicity index of 100% in all cases. At the same time, the energetically most favorable orientations of the peptides relative to the membrane depend on their hydrophobic properties: nonpolar polyleucine is entirely immersed in the bilayer and the hydrophobic peptide with polar groups at the termini adopts a transbilayer orientation, whereas the amphiphilic peptide lies at the interface parallel to the membrane plane. The results of the simulations agree well with the available experimental data for these systems. In the following communications of this series, we plan to describe applications of the solvation model to membrane-bound proteins and peptides with biologically important functional activities.     

Modeling of peptides and proteins in a membrane environment: I. A solvation model mimicking a lipid bilayer

A theoretical solvation model of peptides and proteins that mimics the heterogeneous membrane-water system was proposed. Our approach is based on the combined use of atomic parameters of solvation for water and hydrocarbons, which approximates the hydrated polar groups and acyl chains of lipids, respectively. This model was tested in simulations of several peptides: a nonpolar 20-mer polyleucine, a hydrophobic peptide with terminal polar groups, and a strongly amphiphilic peptide. The conformational space of the peptides in the presence of the membrane was studied by the Monte Carlo method. Unlike a polar solvent and vacuum, the membrane-like environment was shown to stabilize the α-helical conformation: low-energy structures have a helicity index of 100% in all cases. At the same time, the energetically most favorable orientations of the peptides relative to the membrane depend on their hydrophobic properties: nonpolar polyleucine is entirely immersed in the bilayer and the hydrophobic peptide with polar groups at the termini adopts a transbilayer orientation, whereas the amphiphilic peptide lies at the interface parallel to the membrane plane. The results of the simulations agree well with the available experimental data for these systems. In the following communications of this series, we plan to describe applications of the solvation model to membrane-bound proteins and peptides with biologically important functional activities.     

An open-channel blocker interacts with adjacent turns of alpha-helices in the nicotinic acetylcholine receptor.

The binding site for an open-channel blocker, QX-222, at mouse muscle nicotinic acetylcholine receptors was probed using site-directed mutagenesis, oocyte expression, and electrophysiological analysis. The proposed cytoplasmic end of the M2 transmembrane helix is termed position 1'. At position 10' (alpha S252, beta T263, gamma A261, delta A266), Ala residues yield stronger and longer binding of QX-222 than Ser or Thr residues. These effects are opposite and roughly equal (30%-50% per mutation) to previously reported effects at position 6'. The polar end of an anesthetic molecule seems to bind to the position 6' OH groups, which provide a water-like region; the nonpolar moiety is near position 10' and binds more strongly in a nonpolar environment. Interactions with adjacent OH-rich turns of an amphiphilic helix may explain the widespread blocking effects of local anesthetics at the conduction pore of ion channels.     

Semiconductor theory of ion transport in thin lipid membranes: I. Potential and field distributions

In the theory as presented in this paper and the following one, we shall attempt to apply the semiconductor principles and methods to the study of ion transport in thin lipid membranes. Detailed formulations are given on the potential energy barriers at the interfaces, voltage drops in the polar and non-polar regions, and potential and field distributions in the diffuse double layer and within a charged membrane. These results will be used mainly as the boundary conditions for the solution of ion flow as to be given in the following paper. The analysis clearly indicates that the ion transport is interface-limited and is profoundly influenced by the presence of surface charges. An explanation of Na⁺ extrusion in nerve membrane is given based on the field distribution analysis. The theory also suggests that the “membrane potential” depends mainly on surface charges but not necessarily on ion permeation through the membrane.     

A monophasic extraction strategy for the simultaneous lipidome analysis of polar and nonpolar retina lipids

Lipid extraction using a monophasic chloroform/methanol/water mixture, coupled with functional group selective derivatization and direct infusion nano-ESI-high-resolution/accurate MS, is shown to facilitate the simultaneous analysis of both highly polar and nonpolar lipids from a single retina lipid extract, including low abundance highly polar ganglioside lipids, nonpolar sphingolipids, and abundant glycerophospholipids. Quantitative comparison showed that the monophasic lipid extraction method yielded similar lipid distributions to those obtained from established “gold standard” biphasic lipid extraction methods known to enrich for either highly polar gangliosides or nonpolar lipids, respectively, with only modest relative ion suppression effects. This improved lipid extraction and analysis strategy therefore enables detailed lipidome analyses of lipid species across a broad range of polarities and abundances, from minimal amounts of biological samples and without need for multiple lipid class-specific extractions or chromatographic separation prior to analysis.     

Microscopic description of voltage effects on ion-driven cotransport systems

A microscopic model for the analysis of voltage effects on ion-driven cotransport systems is described. The model is based on the notion that the voltage dependence of a given rate constant is directly related to the amount of charge which is translocated in the corresponding reaction step. Charge translocation may result from the movement of an ion along the transport pathway, from the displacement of charged ligand groups of the ion-binding site, or from reorientation of polar residues of the protein in the course of a conformational transition. The voltage dependence of overall transport rate is described by a set of dimensionless coefficients reflecting the dielectric distances over which charge is displaced in the elementary reaction steps. The dielectric coefficients may be evaluated from the shape of the experimental flux-voltage curve if sufficient information on the rate constants of the reaction cycle is available. Examples of flux-voltage curves which are obtained by numerical simulation of the transport model are given for a number of limiting cases.     

Dynamic properties of polyelectrolyte calcium membranes

Summary Shashoua observed spontaneous oscillations in a polyelectrolyte membrane formed by interfacial precipitates of polyacid and polybase. We have here undertaken experimental and theoretical studies of polyglutamic acid-Ca⁺⁺ membrane in order to clarify the processes involved in this dynamic behavior. We find a region of distinct hysteresis in the voltage current curve for this system. A sharp transition from a state of low membrane resistance to one of high resistance occurs at a current density different from that of inverse transition.This membrane system is modeled as a two layer structure: a negatively charged layer made of ionized polyelectrolyte in series with a neutral region in which the polymeric ionic sites are masked by calcium ion. This structure results in a difference in the transference number for the mobile ions, causing salt accumulation at the interfacial region during a current flow in the to direction. This altered salt concentration induces a change of polymeric conformation, which in turn affects the membrane permeability and the rate of accumulation. Based upon nonequilibrium thermodynamic flow equations, and a two-state representation of membrane macromolecular conformation, this model displays a region of hysteresis in the current range of experimental observations.     

Dielectric properties of the polar head group region of zwitterionic lipid bilayers.

A theoretical model describing the dielectric properties of the lipid membrane-water interface region was developed. The rotating polar head groups (e.g. phosphatidylcholine) were simulated as a collection of interacting dipoles imbedded in a nonhomogeneous dielectric. The interactions between the nearest neighborhood were explicitly taken into account, while the other interactions were evaluated by means of the continuum theories. The values of the dielectric constant, its anisotropy and the spontaneous polarization of the interface were evaluated. As an application, we calculated the energy of interaction between an ion and the membrane polar head group region. The results indicate a small spontaneous polarization of the interface (1-1.7 Debyes per lipid molecule) due to the tilting angle of the choline residue with respect to the membrane surface. This dipolar field partially compensates that of opposite orientation originating from the ester group region, giving calculated overall dipolar potentials in better agreement with the experimental data. Our model suggests also a very strong dielectric anisotropy of the interface region, the component of the dielectric constant perpendicular to the membrane plane being much smaller than the parallel component.     

The noneffect of a large linear hydrocarbon, squalene, on the phosphatidylcholine packing structure.

The interaction of squalene with liposomes and monolayers of dipalmitoyl phosphatidylcholine (DPL) has been studied by differential scanning calorimetry, Raman spectroscopy, and surface potential measurements. Mole ratios of squalene to DPL up to 9 to 1 were studied. In contrast to small, nonpolar molecules, which profoundly influence the structure of lipid bilayers as detected by changes in both their thermodynamic phase transition parameters and membrane fluidity, this large, nonpolar, linear hydrocarbon is devoid of such influences. It is clear from our data that a large nonpolar molecule such as squalene, having no polar group that might anchor it to the aqueous interface, cannot intercalate between the acyl chains either below or above the phase transition of DPL. This behavior is not compatible with models that treat the bilayer interior as a bulk hydrocarbon, and suggests that great caution should be exercised in extrapolating partition coefficients based on bulk hydrocarbon measurements to lipid bilayers.     

Temperature-dependent folding pathways of Pin1 WW domain: an all-atom molecular dynamics simulation of a Gō model

We study the folding thermodynamics and kinetics of the Pin1 WW domain, a three-stranded beta-sheet protein, by using all-atom (except nonpolar hydrogens) discontinuous molecular dynamics simulations at various temperatures with a Gō model. The protein exhibits a two-state folding kinetics near the folding transition temperature. A good agreement between our simulations and the experimental measurements by the Gruebele group has been found, and the simulation sheds new insights into the structure of transition state, which is hard to be straightforwardly captured in experiments. The simulation also reveals that the folding pathways at approximately the transition temperature and at low temperatures are much different, and an intermediate state at a low temperature is predicted. The transition state of this small beta-protein at its folding transition temperature has a well-established hairpin 1 made of beta1 and beta2 strands while its low-temperature kinetic intermediate has a formed hairpin 2 composed of beta2 and beta3 strands. Theoretical results are compared with other simulation results as well as available experimental data. This study confirms that specific side-chain packing in an all-atom Gō model can yield a reasonable prediction of specific folding kinetics for a given protein. Different folding behaviors at different temperatures are interpreted in terms of the interplay of entropy and enthalpy in folding process.     

A Stochastic Model of Conductance Transitions in Voltage-Gated IonChannels

We present a statistical physics model to describe the stochastic behaviorof ion transport and channel transitions under an applied membrane voltage.To get pertinent ideas we apply our general theoretical scheme to ananalytically tractable model of the channel with a deep binding site whichinteracts with the permeant ions electrostatically. It is found that theinteraction is modulated by the average ionic occupancy in the bindingsite, which is enhanced by the membrane voltage increases. Above acritical voltage, the interaction gives rise to a emergence of a newconducting state along with shift of S4 charge residues in the channel.This exploratory study calls for further investigations to correlate thecomplex transition behaviors with a variety of ion channels, withparameters in the model, potential energy parameters, voltage, and ionicconcentration.     

Mechanism of the electric response of lipid bilayers to bitter substances.

In order to clarify by what mechanism the lipid bilayer membrane changes its potential under the stimulation of bitter substances, a microscopic model for the effects of the substances on the membrane is presented and studied theoretically. It is assumed that the substances are adsorbed on the membrane and change the partition coefficients of ions between the membrane and the stimulation solution, the dipole orientation in the polar head, and the diffusion constants of ions in the membrane. It is shown, based on the comparison of the calculated results with the experimental ones, that the response arises mainly from a change in the partition coefficients. Protons play an essential role in the membrane potential variation due to the change in their partition coefficients. The present model reproduces the following observed unique properties in the response of lipid bilayers to bitter substances, which cannot be accounted for by the usual channel model for the membrane potential: 1) the response of the membrane potential appears even under the condition that there is no ion gradient across the membrane, 2) the response remains even when the salt in the stimulating solution is replaced with the salt made of an impermeable cation, and 3) the direction of the polarization of the potential is not reversed, even when the ion gradient across the bilayer is reversed.     

Kinetics of Excitable Membranes : Voltage amplification in a diffusion regime

An understanding of the properties of excitable membranes requires the calculation of ion flow through the membrane, including the effects of nonuniformity in the transverse membrane properties (mobilities, fixed charge, electric field). Permeability is apparently controlled at the external interface. Two factors may be involved here: the statistical blocking of pores by divalent cations, and activation energy. Only the former is included in the present treatment. When the total transmembrane voltage is varied, a redistribution in ionic concentration occurs. This can cause a change in boundary (zeta) potential, large in comparison with the applied voltage change—"voltage amplification." The result is a steep change in membrane conductance. The calculated flow curves are compared with experimental results. The Appendix gives an outline of the numerical method used for solving the boundary value problem with several diffusible ions, across a nonuniform regime.     

Electrostatic modeling of dipole-ion interactions in gramicidinlike channels

Using an electrostatic model for the pore and membrane region in a gramicidinlike channel, the effect of dipoles located inside the membrane on the ion transport are analyzed. Calculated energy profiles for different orientations of dipoles show a predominant influence of their radial components. The results qualitatively agree with experimental measurements of conductance on different modified gramicidins and allow to understand the important role of polar side chains on ion permeation.     

Interactions of cholesterol with lipid bilayers: the preferred configuration and fluctuations.

The free energy difference associated with the transfer of a single cholesterol molecule from the aqueous phase into a lipid bilayer depends on its final location, namely on its insertion depth and orientation within the bilayer. We calculated desolvation and lipid bilayer perturbation contributions to the water-to-membrane transfer free energy, thus allowing us to determine the most favorable location of cholesterol in the membrane and the extent of fluctuations around it. The electrostatic and nonpolar contributions to the solvation free energy were calculated using continuum solvent models. Lipid layer perturbations, resulting from both conformational restrictions of the lipid chains in the vicinity of the (rigid) cholesterol backbone and from cholesterol-induced elastic deformations, were calculated using a simple director model and elasticity theory, respectively. As expected from the amphipathic nature of cholesterol and in agreement with the available experimental data, our results show that at the energetically favorable state, cholesterol's hydrophobic core is buried within the hydrocarbon region of the bilayer. At this state, cholesterol spans approximately one leaflet of the membrane, with its OH group protruding into the polar (headgroup) region of the bilayer, thus avoiding an electrostatic desolvation penalty. We found that the transfer of cholesterol into a membrane is mainly driven by the favorable nonpolar contributions to the solvation free energy, whereas only a small opposing contribution is caused by conformational restrictions of the lipid chains. Our calculations also predict a strong tendency of the lipid layer to elastically respond to (thermally excited) vertical fluctuations of cholesterol so as to fully match the hydrophobic height of the solute. However, orientational fluctuations of cholesterol were found to be accompanied by both an elastic adjustment of the surrounding lipids and by a partial exposure of the hydrophobic cholesterol backbone to the polar (headgroup) environment. Our calculations of the molecular order parameter, which reflects the extent of orientational fluctuations of cholesterol in the membrane, are in good agreement with available experimental data.