Abstract In biology, heterosynaptic plasticity maintains homeostasis in synaptic inputs during associative learning and memory, and initiates long-term changes in synaptic strengths that nonspecifically modulate different synapse types. In bioinspired neuromorphic circuits, heterosynaptic plasticity may be used to extend the functionality of two-terminal, biomimetic memristors. In this article, we explore how changes in the pH of droplet interface bilayer aqueous solutions modulate the memristive responses of a lipid bilayer membrane in the pH range 4.97–7.40. Surprisingly, we did not find conclusive evidence for pH-dependent shifts in the voltage thresholds ( V* ) needed for alamethicin ion channel formation in the membrane. However, we did observe a clear modulation in the dynamics of pore formation with pH in time-dependent, pulsed voltage experiments. Moreover, at the same voltage, lowering the pH resulted in higher steady-state currents because of increased numbers of conductive peptide ion channels in the membrane. This was due to increased partitioning of alamethicin monomers into the membrane at pH 4.97, which is below the pKa (~5.3–5.7) of carboxylate groups on the glutamate residues of the peptide, making the monomers more hydrophobic. Neutralization of the negative charges on these residues, under acidic conditions, increased the concentration of peptide monomers in the membrane, shifting the equilibrium concentrations of peptide aggregate assemblies in the membrane to favor greater numbers of larger, increasingly more conductive pores. It also increased the relaxation time constants for pore formation and decay, and enhanced short-term facilitation and depression of the switching characteristics of the device. Modulating these thresholds globally and independently of alamethicin concentration and applied voltage will enable the assembly of neuromorphic computational circuitry with enhanced functionality. Impact statement We describe how to use pH as a modulatory “interneuron” that changes the voltage-dependent memristance of alamethicin ion channels in lipid bilayers by changing the structure and dynamical properties of the bilayer. Having the ability to independently control the threshold levels for pore conduction from voltage or ion channel concentration enables additional levels of programmability in a neuromorphic system. In this article, we note that barriers to conduction from membrane-bound ion channels can be lowered by reducing solution pH, resulting in higher currents, and enhanced short-term learning behavior in the form of paired-pulse facilitation. Tuning threshold values with environmental variables, such as pH, provide additional training and learning algorithms that can be used to elicit complex functionality within spiking neural networks. Graphical abstract
more »
« less
Using Simulation to Understand the Role of Titration on the Stability of a Peptide–Lipid Bilayer Complex
The pH-low insertion peptide (pHLIP) is an anionic membrane-active peptide with promising potential for applications in imaging of cancer tumors and targeted delivery of chemotherapeutics. The key advantage of pHLIP lies in its acid sensitivity: in acidic cellular environments, pHLIP can insert unidirectionally into the plasma membrane. Partitioning–folding coupling is triggered by titration of the acidic residues in pHLIP, transforming pHLIP from a hydrophilic to a hydrophobic peptide. Despite this knowledge, the reverse pathway that leads to exit of the peptide from the plasma membrane is poorly understood. Our hypothesis is that sequential deprotonation of pHLIP is a prerequisite for exit of the peptide from the plasma membrane. We carried out molecular dynamics (MD) simulations to characterize the effect that deprotonation of the acidic residues of pHLIP has on the stability of the peptide when inserted into a model lipid bilayer of 1-palmitoyl-2-oleoyl-sn-3-phosphocholine (POPC). Initiation of the exit mechanism is facilitated by a complex relationship between the peptide, bulk solvent, and the membrane environment. As the N-terminal acidic residues of pHLIP are deprotonated, localized loss of helicity drives unfolding of the peptide and more pronounced interactions with the bilayer at the lipid–water interface. Deprotonation of the C-terminal acidic residues (D25, D31, D33, and E34) leads to further loss of secondary structure distal from the C-terminus, as well as formation of a water channel that stabilizes the orientation of pHLIP parallel to the membrane normal. Together, these results help explain how stabilization of intermediates between the surface-bound and inserted states of pHLIP occur and provide insights into rational design of pHLIP variants with modified abilities of insertion.
more »
« less
- Award ID(s):
- 1726534
- NSF-PAR ID:
- 10198068
- Date Published:
- Journal Name:
- Langmuir
- Volume:
- 36
- ISSN:
- 0743-7463
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Coarse-grained (CG) models have been successful in simulating the chemical properties of lipid bilayers, but accurate treatment of membrane proteins and lipid-protein molecular interactions remains a challenge. The CgProt force field, original developed with the multiscale coarse graining method, is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. Reassignment of select CG sidechain sites from the apolar to polar site type was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlates with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and Wimley-White peptides. These experimental values serve as important anchor points in choosing between alternate CG models based on their observed permeation profiles, particularly for Arg, Lys and Gln residues where the all-atom OPLS solvation energy does not agree well with experiment. Available partitioning data was also used to reparameterize the representation of the peptide backbone, which needed to be made less attractive for the bilayer hydrophobic core region. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membrane-associated peptides LS3 and WALP23. CgProt will find application in studies of lipid-protein interactions and the conformational properties of diverse membrane protein systems.more » « less
-
more » « less
Abstract Membrane tension plays an inhibitory role in clathrin-mediated endocytosis (CME) by impeding the transition of flat plasma membrane to hemispherical clathrin-coated structures (CCSs). Membrane tension also impedes the transition of hemispherical domes to omega-shaped CCSs. However, CME is not completely halted in cells under high tension conditions. Here we find that epsin, a membrane bending protein which inserts its N-terminus H0helix into lipid bilayer, supports flat-to-dome transition of a CCS and stabilizes its curvature at high tension. This discovery is supported by molecular dynamic simulation of the epsin N-terminal homology (ENTH) domain that becomes more structured when embedded in a lipid bilayer. In addition, epsin has an intrinsically disordered protein (IDP) C-terminus domain which induces membrane curvature via steric repulsion. Insertion of H0helix into lipid bilayer is not sufficient for stable epsin recruitment. Epsin’s binding to adaptor protein 2 and clathrin is critical for epsin’s association with CCSs under high tension conditions, supporting the importance of multivalent interactions in CCSs. Together, our results support a model where the ENTH and unstructured IDP region of epsin have complementary roles to ensure CME initiation and CCS maturation are unimpeded under high tension environments.
-
Activation of SARS-CoV-2 Spike deploys its fusion peptide to a membrane of the host cell to infect it. NMR in solution demonstrates that this fusion peptide transforms from intrinsic disorder in solution into a wedge-shaped structure inserted in bilayered micelles. According to NOEs and proximity to a nitroxide spin label deep in the membrane mimic, the globular fold of three helices contrasts the open, extended conformations observed in compact prefusion states. In the hydrophobic, narrow end of the wedge, helices 1 and 2 contact the fatty acyl chains of phospholipids. 50 of the resulting paramagnetic NMR relaxation enhancements and 6 lipid-protein NOEs provided ambiguous distances as collective variables (colvars) to bias and guide MD simulations. Simulations in NAMD using the CHARMM36 forcefield included colvars for 130 medium- and long-range NOEs to maintain the equilibrium structure. In the gently NMR-biased simulations, the fusion peptide maintained its insertion of helices 1 and 2 within a single leaflet while helix 3 remained exposed. A cation occasionally visited the anionic side chains in the loop joining helices 2 and 3 or at the N-terminal end of helix 1. The unoccupied leaflet is thinned and distorted opposite the fusion peptide.The thinning could be related to the fusion peptide promoting formation of the hemi-fusion intermediate in the process of viral-cell fusion. Supported by NSF Rapid award 2030473.more » « less
-
more » « less
Single amino acid mutations provide quantitative insight into the energetics that underlie the dynamics and folding of membrane proteins. Chemical denaturation is the most widely used assay and yields the change in unfolding free energy (ΔΔ
G ). It has been applied to >80 different residues of bacteriorhodopsin (bR), a model membrane protein. However, such experiments have several key limitations: 1) a nonnative lipid environment, 2) a denatured state with significant secondary structure, 3) error introduced by extrapolation to zero denaturant, and 4) the requirement of globally reversible refolding. We overcame these limitations by reversibly unfolding local regions of an individual protein with mechanical force using an atomic-force-microscope assay optimized for 2 μs time resolution and 1 pN force stability. In this assay, bR was unfolded from its native bilayer into a well-defined, stretched state. To measure ΔΔG , we introduced two alanine point mutations into an 8-amino-acid region at the C-terminal end of bR’s G helix. For each, we reversibly unfolded and refolded this region hundreds of times while the rest of the protein remained folded. Our single-molecule–derived ΔΔG for mutant L223A (−2.3 ± 0.6 kcal/mol) quantitatively agreed with past chemical denaturation results while our ΔΔG for mutant V217A was 2.2-fold larger (−2.4 ± 0.6 kcal/mol). We attribute the latter result, in part, to contact between Val217and a natively bound squalene lipid, highlighting the contribution of membrane protein–lipid contacts not present in chemical denaturation assays. More generally, we established a platform for determining ΔΔG for a fully folded membrane protein embedded in its native bilayer.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1021/acs.langmuir.0c02038
