Binding of SV40's Viral Capsid Protein VP1 to Its Glycosphingolipid Receptor GM1 Induces Negative Membrane Curvature: A Molecular Dynamics Study.

The binding of the pentameric capsid protein VP1 of simian virus 40 to its glycosphingolipid receptor GM1 is a key step for the entry of the virus into the host cell. Recent experimental studies have shown that the interaction of variants of soluble VP1 pentamers with giant unilamellar vesicles composed of GM1, DOPC, and cholesterol leads to the formation of tubular membrane invaginations to the inside of the vesicles, mimicking the initial steps of endocytosis. We have used coarse-grained and atomistic molecular dynamics (MD) simulations to study the interaction of VP1 with GM1/DOPC/cholesterol bilayers. In the presence of one VP1 protein, we monitor the formation of small local negative curvature and membrane thinning at the protein binding site as well as reduction of area per lipid. These membrane deformations are also observed under cholesterol-free conditions. However, here, the number of GM1 molecules attached to the VP1 binding pockets increases. The membrane curvature is slightly increased for asymmetric GM1 distribution that mimics conditions in vivo, compared to symmetric GM1 distributions which are often applied in experiments. Slightly smaller inward curvature was observed in atomistic control simulations. Binding of four VP1 proteins leads to an increase of the average intrinsic area per lipid in the protein binding leaflet. Membrane fluctuations appear to be the driving force of VP1 aggregation, as was previously shown for membrane-adhering particles because no VP1 aggregation is observed in the absence of a lipid membrane.

Thermodynamic Origin of Multilayer Structures in Langmuir Polymer Films.

The emergence of polymer-free water surface in a Langmuir polymer film at conditions where a homogeneous coverage has been expected previously is explained on the basis of the surface tensions of polymer and water, γpv and γwv, respectively, as well as the interfacial tension between the two materials, γpw. The polymer molecules considered are 22-residue poly(γ-benzyl-l-glutamate) (PBLG) peptides in α-helical conformation. Values for γpv and γpw derived from MD simulations are consistent with values inferred from experiments considering the emergence of polymer-free surface area for ultrathin films studied using the surface forces apparatus in earlier work. Based on these surface properties, the behavior of individual PBLG peptides at the air-water interface, the dimerization of PBLG peptides, the equilibrium height and width of fibers with given cross section, and the lateral fusion of fibers are described. We show that a prerequisite for the emergence of multilayer structures, which appear locally in domains of sizes of tens to hundreds of micrometers in the considered Langmuir polymer film, is that the condition γpv + γpw - γwv > 0 holds true.

Structure Formation in Langmuir Peptide Films As Revealed from Coarse-Grained Molecular Dynamics Simulations.

Molecular dynamics simulations in conjunction with the Martini coarse-grained model have been used to investigate the (nonequilibrium) behavior of helical 22-residue poly(γ-benzyl-l-glutamate) (PBLG) peptides at the water/vapor interface. Preformed PBLG mono- or bilayers homogeneously covering the water surface laterally collapse in tens of nanoseconds, exposing significant proportions of empty water surface. This behavior was also observed in recent AFM experiments at similar areas per monomer, where a complete coverage had been assumed in earlier work. In the simulations, depending on the area per monomer, either elongated clusters or fibrils form, whose heights (together with the portion of empty water surface) increase over time. Peptides tend to align with respect to the fiber axis or with the major principal axis of the cluster, respectively. The aspect ratio of the cluster observed is 1.7 and, hence, comparable to though somewhat smaller than the aspect ratio of the peptides in α-helical conformation, which is 2.2. The heights of the fibrils is 3 nm after 20 ns and increases to 4.5 nm if the relaxation time is increased by 2 orders of magnitude, in agreement with the experiment. Aggregates with heights of about 3 or 4.5 nm are found to correspond to local bi- or trilayer structures, respectively.

Mechanisms for allosteric activation of protease DegS by ligand binding and oligomerization as revealed from molecular dynamics simulations.

A local perturbation of a protein may lead to functional changes at some distal site, a phenomenon denoted as allostery. Here, we study the allosteric control of a protease using molecular dynamics simulations. The system considered is the bacterial protein DegS which includes a protease domain activated on ligand binding to an adjacent PDZ domain. Starting from crystallographic structures of DegS homo-trimers, we perform simulations of the ligand-free and -bound state of DegS at equilibrium. Considering a single protomer only, the trimeric state was mimicked by applying restraints on the residues in contact with other protomers in the DegS trimer. In addition, the bound state was also simulated without any restraints to mimic the monomer. Our results suggest that not only ligand release but also disassembly of a DegS trimer inhibits proteolytic activity. Considering various observables for structural changes, we infer allosteric pathways from the interface with other protomers to the active site. Moreover, we study how ligand release leads to (i) catalytically relevant changes involving residues 199-201 and (ii) a transition from a stretched to a bent conformation for residues 217-219 (which prohibits proper substrate binding). Finally, based on ligand-induced Cα shifts we identify residues in contact with other protomers in the DegS trimer that likely transduce the perturbation from ligand release from a given protomer to adjacent protomers. These residues likely play a key role in the experimentally known effect of ligand release from a protomer on the proteolytic activity of the other protomers. Proteins 2016; 84:1690-1705. © 2016 Wiley Periodicals, Inc.

Stable polyglutamine dimers can contain ß-hairpins with interdigitated side chains-but not α-helices, ß-nanotubes, ß-pseudohelices, or steric zippers.

A common thread connecting nine fatal neurodegenerative protein aggregation diseases is an abnormally expanded polyglutamine tract found in the respective proteins. Although the structure of this tract in the large mature aggregates is increasingly well described, its structure in the small early aggregates remains largely unknown. As experimental evidence suggests that the most toxic species along the aggregation pathway are the small early ones, developing strategies to alleviate disease pathology calls for understanding the structure of polyglutamine peptides in the early stages of aggregation. Here, we present a criterion, grounded in available experimental data, that allows for using kinetic stability of dimers to assess whether a given polyglutamine conformer can be on the aggregation path. We then demonstrate that this criterion can be assessed using present-day molecular dynamics simulations. We find that although the α-helical conformer of polyglutamine is very stable, dimers of α-helices lack the kinetic stability necessary to support further oligomerization. Dimers of steric zipper, ß-nanotube, and ß-pseudohelix conformers are also too short-lived to initiate aggregation. The ß-hairpin-containing conformers, instead, invariably form very stable dimers when their side chains are interdigitated. Combining these findings with the implications of recent solid-state NMR data on mature fibrils, we propose a possible pathway for the initial stages of polyglutamine aggregation, in which ß-hairpin-containing conformers act as templates for fibril formation.

The free energy of nanopores in tense membranes.

Membrane nanopores are central players for a range of important cellular membrane remodeling processes as well as membrane rupture. Understanding pore formation in tense membranes requires comprehension of the molecular mechanism of pore formation and the associated free energy change as a function of the membrane tension. Here we propose a scheme to calculate the free energy change associated with the formation of a nanometer sized pore in molecular dynamics simulations as a function of membrane tension, which requires the calculation of only one computationally expensive potential of mean force. We show that membrane elastic theory can be used to estimate the pore formation free energy at different tension values from the free energy change in a relaxed membrane and the area expansion curves of the membranes. We have computed the pore formation free energy for a dipalmitoyl-phosphatidylcholine (DPPC) membrane at two different lateral pressure values, 1 bar and -40 bar, by calculating the potential of mean force acting on the head group of a single lipid molecule. Unrestrained simulations of the closing process confirm that the intermediate states along this reaction coordinate are reasonable and show that hydrophilic indentations spanning half the bilayer connected by a hydrophobic pore segment represent the corresponding high energy transition state. A comparison of the stability of simulated membranes to experiment at high loading rates show that, contrary to expectation, pores form too easily in small simulated membrane patches. This discrepancy originates from a combination of the absence of ions in the simulations and the small membrane size.

Allosteric control of kinesin's motor domain by tubulin: a molecular dynamics study.

Molecular motors such as kinesin are essential for many biological processes. These motors have two motor domains, which bind to tubulin filaments, hydrolyze ATP, and transduce the released chemical energy into directed movements. The general principles of this chemomechanical coupling are now well-established but the underlying molecular mechanisms remain elusive because small conformational changes within large proteins are difficult to detect experimentally. Here, we use atomistic molecular dynamics simulations to monitor such changes within a single motor domain of KIF1A, which belongs to the kinesin-3 motor family. The nucleotide binding pocket of this domain can be empty or occupied by ATP or ADP. For these three nucleotide states, we determine the mobility of the backbone of the protein, both in solution and attached to tubulin. Only one subdomain of the motor domain is found to exhibit a strongly increased mobility upon binding to tubulin: the neck linker that presumably acts as a mechanical transmitter to the other motor domain in dimeric kinesin-3 motors. Furthermore, upon binding to tubulin, the neck linker mobility becomes sensitive to the bound nucleotide and is highly increased after phosphate release, which implies undocking of this linker from the core of the motor domain. These simulation results are consistent with experimental data from EPR spectroscopy, FRET, and cryo-electron microscopy. A detailed analysis of our simulation data also reveals that the undocking of the neck linker in the ADP-kinesin-tubulin state arises from allosteric interactions between the nucleotide and tubulin and that the ß-sheet core undergoes a twist both during phosphate release and ATP binding. The computational approach used here can be applied to other motor domains and mechanoenzymes in order to identify allosteric interactions between the subdomains of these proteins.

Specific binding of chloride ions to lipid vesicles and implications at molecular scale.

Biological membranes composed of lipids and proteins are in contact with electrolytes like aqueous NaCl solutions. Based on molecular dynamics studies it is widely believed that Na(+) ions specifically bind to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membranes, whereas Cl(-) ions stay in solution. Here, we present a careful comparison of recent data from electrophoresis and isothermal titration calorimetry experiments as well as molecular dynamics simulations suggesting that in fact both ions show very similar affinities. The corresponding binding constants are 0.44(±0.05) M(-1) for Na(+) and 0.40(±0.04) M(-1) for Cl(-) ions. This is highlighted by our observation that a widely used simulation setup showing asymmetric affinities of Na(+) and Cl(-) for POPC bilayers overestimates the effect of NaCl on the electrophoretic mobility of a POPC membrane by an order of magnitude. Implications for previous simulation results on the effect of NaCl on polarization of interfacial water, transmembrane potentials, and mechanisms for ion transport through bilayers are discussed. Our findings suggest that a range of published simulations results on the interaction of NaCl with phosphocholine bilayers have to be reconsidered and revised and that force field refinements are necessary for reliable simulation studies of membranes at physiological conditions on a molecular level.

Antimicrobial selectivity based on zwitterionic lipids and underlying balance of interactions.

An important feature of antimicrobial peptides is their ability to distinguish pro- from eukaryotic membranes. In vitro experiments on the antimicrobial peptide NK-2 indicate that the discrimination between zwitterionic phosphatidylethanolamine lipids exposed by prokaryotes and phosphatidylcholine lipids exposed by eukaryotes plays an important role. The underlying mechanism is not understood. Here we present molecular dynamics simulations in conjunction with a coarse grained model and thermodynamic integration showing that NK-2 binds more strongly to palmitoyloleoylphosphatidylethanolamine (POPE) than to palmitoyloleoylphosphatidylcholine (POPC) bilayers. Finite size effects on the relative free energy have been corrected for with a method that may also be useful in future studies of the affinities of macromolecules for lipid membranes. Our results support the previous hypothesis that the stronger binding to PE compared to PC arises from a better accessibility of the phosphates of the lipids to the cationic peptide in a sense that a similar number of peptide-lipid salt bridges requires to break more favorable electrostatic headgroup-headgroup interactions for PC relative to PE. The transfer of NK-2 from POPC to POPE is found to lead to a decrease in electrostatic peptide-lipid but an increase in lipid-lipid and ion-lipid interactions, correlating with a dehydration of the lipids and the ions but an increased hydration of the peptide. The increase in affinity of NK-2 for POPE compared to POPC hence arises from a complex interplay of competing interactions. This work opens the perspective to study how the affinity of antimicrobial peptides changes with amino acid sequence and lipid composition.

Macro- versus microscopic view on the electrokinetics of a water-membrane interface.

Electrophoresis is an experimental method widely used to study electrostatic properties of interfaces. Here, we question the validity of the macroscopic theory for the planar geometry by Helmholtz and Smoluchowski by considering a POPC bilayer in an aqueous solution with 500 mM NaCl, using molecular dynamics simulations. We find that POPC shows positive electrophoretic mobility due to adsorption of sodium ions at the lipid headgroups. The theory assumes that the region in which the water density undergoes a transition from the bulk value to zero (interfacial width) is small compared to the Debye screening length. This separation of length scale is not fullfilled in the present case. Hence, contrasting the theory, we observe that the surface is not sharply defined, continuum hydrodynamics is not applicable, the effective viscosity in the double layer is increased compared to the bulk, and the zeta potential is dominated by the dipole potential. Our results might have widespread implications for interpretation of electrokinetic studies in general.

Effect of tension and curvature on the chemical potential of lipids in lipid aggregates.

Understanding the factors that influence the free energy of lipids in bilayer membranes is an essential step toward understanding exchange processes of lipids between membranes. In general, both lipid composition and membrane geometry can affect lipid exchange rates between bilayer membranes. Here, the free energy change ΔG(des) for the desorption of dipalmitoyl-phosphatidylcholine (DPPC) lipids from different lipid aggregates has been computed using molecular dynamics simulations and umbrella sampling. The value of ΔG(des) is found to depend strongly on the local properties of the aggregate, in that both tension and curvature lead to an increase in ΔG(des). A detailed analysis shows that the increased desorption free energy for tense bilayers arises from the increased conformational entropy of the lipid tails, which reduces the favorable component -TΔS(L) of the desorption free energy.

Negligible water surface charge determined using Kelvin probe and total reflection X-ray fluorescence techniques.

The water surface charge has been extensively debated in recent decades. Electrophoretic mobilities of air bubbles in water and disjoining pressures between the surfaces of aqueous films suggest that the surface of water exhibits a significant negative charge. This is commonly attributed to a strong adsorption of hydroxide ions at the interface, though spectroscopic measurements and simulation studies suggest surface depletion of hydroxide ions. Alternatively, the negative surface charge could arise from surface contamination with trace charged surfactants. We have probed the variation in the surface charge of water with pH by measuring surface potentials using the Kelvin probe technique. Independently, the abundance in the interfacial layer of "reporter ions" (Rb(+) and Br(-)), which must be affected by a charged surface, has been monitored using the total reflection X-ray fluorescence (TRXF) technique. Special care was taken to prove the high sensitivity of this technique as well as to avoid surface contaminants. The magnitude of the surface charge was found to be below 1 e per 500 nm(2) (TRXF). No evidence of variations in the surface potential between pH 2-3 and pH 9-12 was detected within the accuracies of the methods (5 mV for Kelvin probe and 2 mV for TRXF). Hence, our findings suggest that the clean water surface exhibits negligible charge in a wide pH range.

Common force field thermodynamics of cholesterol.

Four different force fields are examined for dynamic characteristics using cholesterol as a case study. The extent to which various types of internal degrees of freedom become thermodynamically relevant is evaluated by means of principal component analysis. More complex degrees of freedom (angle bending, dihedral rotations) show a trend towards force field independence. Moreover, charge assignments for membrane-embedded compounds are revealed to be critical with significant impact on biological reasoning.

Competing interactions for antimicrobial selectivity based on charge complementarity.

Antimicrobial peptides (AMPs) are an evolutionary conserved component of the innate immune system and possible templates for the development of new antibiotics. An important property of antimicrobial peptides is their ability to discriminate bacterial from eucaryotic cells which is attributed to the difference in lipid composition of the outer leaflet of the plasma membrane between the two types of cells. Whereas eucaryotic cells usually expose zwitterionic lipids, procaryotic cells expose also anionic lipids which bind the cationic antimicrobial peptides electrostatically. An example is the antimicrobial peptide NK-2 which is highly cationic and favors binding to anionic membranes. In the present study, the difference in binding affinity of NK-2 for palmitoyl-oleoyl-phosphatidyl-glycerol (POPG) and palmitoyl-oleoyl-phosphatidyl-choline (POPC) is studied using molecular dynamics simulations in conjunction with a coarse grained model and thermodynamic integration, by computing the change in free energy and its components upon the transfer of NK-2 from POPC to POPG. The transfer is indeed found to be highly favorable. Interestingly, the favorable contribution from the electrostatic interaction between the peptide and the anionic lipids is overcompensated by an unfavorable contribution from the change in lipid-cation interactions due to the release of counterions from the lipids. The increase in entropy due to the release of the cations is compensated by other entropic components. The largest favorable contribution arises from the solvation of the counterions. Overall the interaction between NK-2 and POPG is not determined by a single driving force but a subtle balance of competing interactions.

The glomerular filtration barrier function: new concepts.

PURPOSE OF REVIEW: Each day, the human kidneys filter about 140 l of primary urine from plasma. Although this ultrafiltrate is virtually free of plasma protein, the glomerular filter never clogs under physiological conditions. Upto today it is still not entirely resolved as to how the kidney accomplishes this extraordinary task. Most of the proposed models for glomerular filtration have not considered electrical effects. RECENT FINDINGS: In micropuncture studies, we have directly measured an electrical field across the glomerular filtration barrier. This potential difference is most likely generated by forced passage of the ionic solution of the plasma across the charged glomerular filter ('electrokinetic potential'). As all plasma proteins are negatively charged, the electrical field across the glomerular filtration barrier is predicted to drive plasma proteins from the filter toward the capillary lumen by electrophoresis. SUMMARY: In this review, we examine our novel model for glomerular filtration in more detail. We outline the physical mechanisms by which electrokinetic effects (streaming potentials) are generated. We investigate the potential impact of the electrical field on the passage of albumin across the glomerular filtration barrier. We review the mathematical heteroporous model including electrical effects and analyse a selection of experimental studies for indications that electrical effects influence glomerular permeability significantly.

Energetic basis for drug resistance of HIV-1 protease mutants against amprenavir.

Amprenavir (APV) is a high affinity (0.15 nM) HIV-1 protease (PR) inhibitor. However, the affinities of the drug resistant protease variants V32I, I50V, I54V, I54M, I84V and L90M to amprenavir are decreased 3 to 30-fold compared to the wild-type. In this work, the popular molecular mechanics Poisson-Boltzmann surface area method has been used to investigate the effectiveness of amprenavir against the wild-type and these mutated protease variants. Our results reveal that the protonation state of Asp25/Asp25' strongly affects the dynamics, the overall affinity and the interactions of the inhibitor with individual residues. We emphasize that, in contrast to what is often assumed, the protonation state may not be inferred from the affinities but requires pK(a) calculations. At neutral pH, Asp25 and Asp25' are ionized or protonated, respectively, as suggested from pK(a) calculations. This protonation state was thus mainly considered in our study. Mutation induced changes in binding affinities are in agreement with the experimental findings. The decomposition of the binding free energy reveals the mechanisms underlying binding and drug resistance. Drug resistance arises from an increase in the energetic contribution from the van der Waals interactions between APV and PR (V32I, I50V, and I84V mutant) or a rise in the energetic contribution from the electrostatic interactions between the inhibitor and its target (I54M and I54V mutant). For the V32I mutant, also an increased free energy for the polar solvation contributes to the drug resistance. For the L90M mutant, a rise in the van der Waals energy for APV-PR interactions is compensated by a decrease in the polar solvation free energy such that the net binding affinity remains unchanged. Detailed understanding of the molecular forces governing binding and drug resistance might assist in the design of new inhibitors against HIV-1 PR variants that are resistant against current drugs.

Validating affinities for ion-lipid association from simulation against experiment.

Understanding biological membranes at physiological conditions requires comprehension of the interaction of lipid bilayers with sodium and potassium ions. These cations are adsorbed at palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers as indicated from previous studies. Here we compare the affinity of Na(+) and K(+) for POPC in molecular dynamics (MD) simulations with recent data from electrophoresis experiments and isothermal calorimetry (ITC) at neutral pH. NaCl and KCl were described using GROMOS or parameters matching solution activities on the basis of Kirkwood-Buff theory (KBFF), and K(+) was also described using parameters by Dang et al., all in conjunction with the Berger parameters for the lipids and the SPC water model. Apparent binding constants of GROMOS-Na(+) and KBFF-K(+) are the same within error and in good agreement with values from ITC. Although these force fields yield the same number of bound ions per number of lipids for Na(+) and K(+), they give a larger number of Na(+) ions per surface area compared to K(+), in agreement with the electrophoresis experiments, because Na(+) causes a stronger reduction in the area per lipid than K(+). The intrinsic binding constants, on the other hand, are reproduced by Dang-K(+) but overestimated by GROMOS-Na(+) and KBFF-K(+). That no ion force field reproduces the intrinsic and the apparent binding constant simultaneously arises from the fact that in MD simulations, implicitly meant to mimic neutral pH, pure PC is usually modeled with zero surface charge. In contrast, POPC at neutral conditions in experiment carries a low but significant negative surface charge and is uncharged only at acidic pH as indicated from electrophoretic mobilities. Implications for future simulation and experimental studies are discussed.

Electrical forces determine glomerular permeability.

There is ongoing controversy about the mechanisms that determine the characteristics of the glomerular filter. Here, we tested whether flow across the glomerular filter generates extracellular electrical potential differences, which could be an important determinant of glomerular filtration. In micropuncture experiments in Necturus maculosus, we measured a potential difference across the glomerular filtration barrier that was proportional to filtration pressure (-0.045 mV/10 cm H2O). The filtration-dependent potential was generated without temporal delay and was negative within Bowman's space. Perfusion with the cationic polymer protamine abolished the potential difference. We propose a mathematical model that considers the relative contributions of diffusion, convection, and electrophoretic effects on the total flux of albumin across the filter. According to this model, potential differences of -0.02 to -0.05 mV can induce electrophoretic effects that significantly influence the glomerular sieving coefficient of albumin. This model of glomerular filtration has the potential to provide a mechanistic theory, based on experimental data, about the filtration characteristics of the glomerular filtration barrier. It provides a unique approach to the microanatomy of the glomerulus, renal autoregulation, and the pathogenesis of proteinuria.

Multiscale Molecular Dynamics Studies Reveal Different Modes of Receptor Clustering by Gb3-Binding Lectins.

The recognition of carbohydrate receptors on host cell membranes by pathogenic lectins is a crucial step in the microbial invasion. Two bacterial lectins, the B-subunit of Shiga toxin from Shigella dysenteria (StxB) and lectin I from Pseudomonas aeruginosa (LecA), are specific to the same galactolipid-globotriaosylceramide (Gb3). In this study we present a coarse-grained (cg) model of Gb3, which we further apply to unravel the molecular details of glycolipid binding by two lectins on the surface of a DOPC/cholesterol/Gb3 bilayer. In cg molecular dynamics simulations with time scales of dozens of microseconds, Gb3 was randomly distributed. The binding of both StxB or LecA is accompanied by Gb3 clustering in a cholesterol environment and with exclusion of DOPC in protein vicinity. StxB being bound by all 15 binding sites induced membrane bending, while LecA interacted with two out of four binding sites for most of the time causing a smaller inward curvature of the model membrane. Stable interactions occurred preferably when LecA was normal to the membrane surface. Furthermore, all-atom simulations revealed that LecA bound Gb3's headgroup at only one out of two possible conformations of the carbohydrate moiety observed at protein-free conditions. The results shed light on the mechanism of interactions between two lectins and Gb3 on the membrane surface and offer a coarse-grained model to study more complex systems at large spatiotemporal scales.

Solvent-exposed tails as prestalk transition states for membrane fusion at low hydration.

Membrane fusion is a key step in intracellular trafficking and viral infection. The underlying molecular mechanism is poorly understood. We have used molecular dynamics simulations in conjunction with a coarse grained model to study early metastable and transition states during the fusion of two planar palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers separated by five waters per lipid in the cis leaflets at zero tension. This system mimics the contact area between two vesicles with large diameters compared to the membrane thickness at conditions where fusion may start in the core of the contact area. At elevated temperatures, the two proximal leaflets become connected via multiple lipid molecules and form a stalklike structure. At room temperature, this structure has a free energy of 3k(B)T and is separated from the unconnected state by a significant free energy barrier of 20k(B)T. Stalk formation is initiated by the establishment of a localized hydrophobic contact between the bilayers. This contact is either formed by two partially splayed lipids or a single fully splayed one leading to the formation of a (metastable) splayed lipid bond intermediate. These findings indicate that, for low hydration, early membrane fusion kinetics is not determined by the stalk energy but by the energy of prestalk transition states involving solvent-exposed lipid tails.