skip to main content


Title: Refining amino acid hydrophobicity for dynamics simulation of membrane proteins
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
Award ID(s):
1516826
NSF-PAR ID:
10108013
Author(s) / Creator(s):
Date Published:
Journal Name:
PeerJ
Volume:
6
ISSN:
2167-8359
Page Range / eLocation ID:
e4230
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Many intrinsically disordered peptides have been shown to undergo liquid–liquid phase separation and form complex coacervates, which play various regulatory roles in the cell. Recent experimental studies found that such phase separation processes may also occur at the lipid membrane surface and help organize biomolecules during signaling events; in some cases, phase separation of proteins at the membrane surface was also observed to lead to significant remodeling of the membrane morphology. The molecular mechanisms that govern the interactions between complex coacervates and lipid membranes and the impacts of such interactions on their structure and morphology, however, remain unclear. Here we study the coacervation of poly-glutamate (E 30 ) and poly-lysine (K 30 ) in the presence of lipid bilayers of different compositions. We carry out explicit-solvent coarse-grained molecular dynamics simulations by using the MARTINI (v3.0) force-field. We find that more than 20% anionic lipids are required for the coacervate to form stable contact with the bilayer. Upon wetting, the coacervate induces negative curvature to the bilayer and facilitates local lipid demixing, without any peptide insertion. The magnitude of negative curvature, extent of lipid demixing, and asphericity of the coacervate increase with the concentration of anionic lipids. Overall, we observe a decrease in the number of contacts among the polyelectrolytes as the droplet spreads over the bilayer. Therefore, unlike previous suggestions, interactions among polyelectrolytes do not constitute a driving force for the membrane bending upon wetting by the coacervate. Rather, analysis of interaction energy components suggests that bending of the membrane is favored by enhanced interactions between polyelectrolytes with lipids as well as with counterions. Kinetic studies reveal that, at the studied polyelectrolyte concentrations, the coacervate formation precedes bilayer wetting. 
    more » « less
  2. Reconstituted photosynthetic proteins which are activated upon exposure to solar energy hold enormous potential for powering future solid state devices and solar cells. The functionality and integration of these proteins into such devices has been successfully enabled by lipid-like peptides. Yet, a fundamental understanding of the organization of these peptides with respect to the photosynthetic proteins and themselves remains unknown and is critical for guiding the design of such light-activated devices. This study investigates the relative organization of one such peptide sequence V 6 K 2 (V: valine and K: lysine) within assemblies. Given the expansive spatiotemporal scales associated with this study, a hybrid coarse-grained (CG) model which captures the structure, conformation and aggregation of the peptide is adopted. The CG model uses a combination of iterative Boltzmann inversion and force matching to provide insight into the relative organization of V 6 K 2 in assemblies. The CG model reproduces the structure of a V 6 K 2 peptide sequence along with its all atom (AA) solvation structure. The relative organization of multiple peptides in an assembly, as captured by CG simulations, is in agreement with corresponding results from AA simulations. Also, a backmapping procedure reintroduces the AA details of the peptides within the aggregates captured by the CG model to demonstrate the relative organization of the peptides. Furthermore, a large number of peptides self-assemble into an elongated micelle in the CG simulation, which is consistent with experimental findings. The coarse-graining procedure is tested for transferability to longer peptide sequences, and hence can be extended to other amphiphilic peptide sequences. 
    more » « less
  3. null (Ed.)
    Abstract

    Synthetic lipid membranes are self-assembled biomolecular double layers designed to approximate the properties of living cell membranes. These membranes are employed as model systems for studying the interactions of cellular envelopes with the surrounding environment in a controlled platform. They are constructed by dispersing amphiphilic lipids into a combination of immiscible fluids enabling the biomolecules to self-assemble into ordered sheets, or monolayers at the oil-water interface. The adhesion of two opposing monolayer sheets forms the membrane, or the double layer. The mechanical properties of these synthetic membranes often differ from biological ones mainly due to the presence of residual solvent in between the leaflets. In fact, the double layer compresses in response to externally applied electrical field with an intensity that varies depending on the solvent present. While typically viewed as a drawback associated with their assembly, in this work the elasticity of the double layer is utilized to further quantify complex biophysical phenomena. The adsorption of charged molecules on the surface of a lipid bilayer is a key property to decipher biomolecule interactions at the interface of the cell membrane, as well as to develop effective antimicrobial peptides and similar membrane-active molecules. This adsorption generates a difference in the boundary potentials on either side of the membrane which may be tracked through electrophysiology. The soft synthetic membranes produced in the laboratory compress when exposed to an electric field. Tracking the minimum membrane capacitance allows for quantifying when the intrinsic electric field produced by the asymmetry is properly compensated by the supplied transmembrane voltage. The technique adopted in this work is the intramembrane field compensation (IFC). This technique focuses on the current generated by the bilayer in response to a sinusoidal voltage with a DC component, VDC. Briefly, the output sinusoidal current is divided into its harmonics and the second harmonic equals zero when VDC compensates the internal electric field. In this work, we apply the IFC technique to droplet interface bilayers (DIB) enabling the development of a biological sensor. A certain membrane elasticity is needed for accurate measurements and is tuned through the solvent selection. The asymmetric DIBs are formed, and an automated PID-controlled IFC design is implemented to rapidly track and compensate the membrane asymmetry. The closed loop system continuously reads the current and generates the corresponding voltage until the second harmonic is abated. This research describes the development and optimization of a biological sensor and examines how varying the structure of the synthetic membrane influences its capabilities for detecting membrane-environment interactions. This platform may be applied towards studying the interactions of membrane-active molecules and developing models for the associated phenomena to enhance their design.

     
    more » « less
  4. null (Ed.)
    We developed coarse-grained models of spike proteins in SARS-CoV-2 coronavirus and angiotensin-converting enzyme 2 (ACE2) receptor proteins to study the endocytosis of a whole coronavirus under physiologically relevant spatial and temporal scales. We first conducted all-atom explicit-solvent molecular dynamics simulations of the recently characterized structures of spike and ACE2 proteins. We then established coarse-grained models using the shape-based coarse-graining approach based on the protein crystal structures and extracted the force field parameters from the all-atom simulation trajectories. To further analyze the coarse-grained models, we carried out normal mode analysis of the coarse-grained models to refine the force field parameters by matching the fluctuations of the internal coordinates with the original all-atom simulations. Finally, we demonstrated the capability of these coarse-grained models by simulating the endocytosis of a whole coronavirus through the host cell membrane. We embedded the coarse-grained models of spikes on the surface of the virus envelope and anchored ACE2 receptors on the host cell membrane, which is modeled using a one-particle-thick lipid bilayer model. The coarse-grained simulations show the spike proteins adopt bent configurations due to their unique flexibility during their interaction with the ACE2 receptors, which makes it easier for them to attach to the host cell membrane than rigid spikes. 
    more » « less
  5. Aggregation of misfolded oligomeric amyloid-beta (Aβ) peptides on lipid membranes has been identified as a primary event in Alzheimer's pathogenesis. However, the structural and dynamical features of this membrane assisted Aβ aggregation have not been well characterized. The microscopic characterization of dynamic molecular-level interactions in peptide aggregation pathways has been challenging both computationally and experimentally. In this work, we explore differential patterns of membrane-induced Aβ 16–22 (K–L–V–F–F–A–E) aggregation from the microscopic perspective of molecular interactions. Physics-based coarse-grained molecular dynamics (CG-MD) simulations were employed to investigate the effect of lipid headgroup charge – zwitterionic (1-palmitoyl-2-oleoyl- sn-glycero -3-phosphocholine: POPC) and anionic (1-palmitoyl-2-oleoyl- sn-glycero -3-phospho- l -serine: POPS) – on Aβ 16–22 peptide aggregation. Our analyses present an extensive overview of multiple pathways for peptide absorption and biomechanical forces governing peptide folding and aggregation. In agreement with experimental observations, anionic POPS molecules promote extended configurations in Aβ peptides that contribute towards faster emergence of ordered β-sheet-rich peptide assemblies compared to POPC, suggesting faster fibrillation. In addition, lower cumulative rates of peptide aggregation in POPS due to higher peptide–lipid interactions and slower lipid diffusion result in multiple distinct ordered peptide aggregates that can serve as nucleation seeds for subsequent Aβ aggregation. This study provides an in-silico assessment of experimentally observed aggregation patterns, presents new morphological insights and highlights the importance of lipid headgroup chemistry in modulating the peptide absorption and aggregation process. 
    more » « less