An official website of the United States government
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
Interaction of Camptothecin with Model Cellular Membranes
Accurate and efficient prediction of drug partitioning in model membranes is of significant interest to the pharmaceutical industry. Herein, we utilize advanced sampling methods, specifically, the adaptive biasing force methodology to calculate the potential of mean force for a model hydrophobic anticancer drug, camptothecin (CPT), across three model interfaces. We consider an octanol bilayer, a thick octanol/water interface, and a model 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/water interface. We characterize the enthalpic and entropic contributions of the drug to the potential of mean force. We show that the rotational entropy of the drug is inversely related to the probability of hydrogen bond formation of the drug with the POPC membrane. In addition, in long-time microsecond simulations of a high concentration of CPT above the POPC membrane, we show that strong drug–drug aromatic interactions shift the spatial orientation of the drug with the membrane. Stacks of hydrophobic drugs form, allowing penetration of the drug just under the POPC head groups. These results imply that inhomogeneous membrane models need to take into account the effect of drug aggregation on the membrane environment.
more »
« less
- Award ID(s):
- 1750694
- PAR ID:
- 10145459
- Date Published:
- Journal Name:
- Journal of Chemical Theory and Computation
- ISSN:
- 1549-9618
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Nano-formulated, combinatory therapeutics that control the spatiotemporal aspects of drug release have potential to overcome many of the challenges faced in cancer therapy. Herein, we describe a peptide nanotube functionalized with two anticancer drugs, 5-fluoruracil (5-FU) and camptothecin (CPT). The nanotube was formed via peptide self-assembly, which positioned 5-FU on the surface at the aqueous interface; whereas, CPT was sequestered within the hydrophobic walls. Thus, two different release profiles were observed: rapid release of 5-FU, followed by slower, sustained production of CPT. This profile emerged from the rapid hydrolytic cleavage of 5-FU at the aqueous/nanotube interface, which produced a smaller nanotube comprised of the peptide fragment.more » « less
-
Nghiem, Long (Ed.)Porous membranes having a particular wetting characteristic, hydrophobic or hydrophilic, are used for nondispersive membrane solvent extraction (MSX) where two immiscible phases flow on two sides of the membrane. The aqueous-organic phase interface across which solvent extraction/back extraction occurs remains immobilized on one surface of the membrane. This process requires the pressure of the phase not present in membrane pores to be either equal to or higher than that of the phase present in membrane pores; the excess phase pressure should not exceed a breakthrough pressure. In countercurrent MSX with significant flow pressure drop in each phase, this often poses a problem. To overcome this problem, flat porous Janus membranes were developed using either a base polypropylene (PP) or polyvinylidene fluoride (PVDF) or polyamide (Nylon) membrane, one side of which is hydrophobic and the other being hydrophilic. Such membranes were characterized using the contact angle, liquid entry pressure (LEP) and the droplet breakthrough pressure from each side of the membrane along with characterizations via scanning electron microscopy (SEM), Fourier transform infrared spectroscopy (FTIR). Nondispersive solvent extractions were carried out successfully for two systems, octanol-phenol (solute)-water, toluene-acetone (solute)-water, with either flowing phase at a pressure higher than that of the other phase. The phenol extraction system had a high solute distribution coefficient whereas acetone prefers both phases almost identically. The potential practical utility of the MSX technique will be substantially enhanced via Janus MSX membranes.more » « less
-
Physical interactions between polypeptide chains and lipid membranes underlie critical cellular processes. Yet, despite fundamental importance, key mechanistic aspects of these interactions remain elusive. Bulk experiments have revealed a linear relationship between free energy and peptide chain length in a model system, but does this linearity extend to the interaction strength and to the kinetics of lipid binding? To address these questions, we utilized a combination of coarse-grained molecular dynamics (CG MD) simulations, analytical modeling, and atomic force microscopy (AFM)-based single molecule force spectroscopy. Following previous bulk experiments, we focused on interactions between short hydrophobic peptides (WLn, n = 1, ..., 5) with 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) bilayers, a simple system that probes peptide primary structure effects. Potentials of mean force extracted from CG MD recapitulated the linearity of free energy with the chain length. Simulation results were quantitatively connected to bulk biochemical experiments via a single scaling factor of order unity, corroborating the methodology. Additionally, CG MD revealed an increase in the distance to the transition state, a result that weakens the dependence of the dissociation force on the peptide chain length. AFM experiments elucidated rupture force distributions and, through modeling, intrinsic dissociation rates. Taken together, the analysis indicates a rupture force plateau in the WLn−POPC system, suggesting that the final rupture event involves the last 2 or 3 residues. In contrast, the linear dependence on chain length was preserved in the intrinsic dissociation rate. This study advances the understanding of peptide−lipid interactions and provides potentially useful insights for the design of peptides with tailored membrane-interacting properties.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1021/acs.jctc.9b00541
