Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU-), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU-, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU- inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU- is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU- to GLU0 to reduce the energy barriers. Overall, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

Study on the Assembly Mechanisms and Transport Properties of Transmembrane End-Charged Cyclic Peptide Nanotubes

In this work, two end-charged cyclic peptide nanotubes (CPNTs) embedded in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) were designed to simulate transmembrane ion channels. Density functional theory (DFT) computations at the level of M06-2X/6-31G give different assembling modes of the negatively charged ELWL-CPNT and positively charged RLWL-CPNT as (L-L)(D-L)(D-D)(L-L)(D-D)(L-L)(D-D) and (D-D)(L-L)(D-D)(L-L)(D-D)(L-L)(D-D), respectively. Molecular dynamics (MD) simulations indicate that a charge at a CPNT end obviously affects the structure of the channel water chain and the diffusion behavior of K⁺. The regions with the highest probability of H-bond defects in the channel water chains are gap5 and gap2 in ELWL/POPE-CPNT and RLWL/POPE-CPNT, respectively. K⁺ can easily enter either CPNT by desolvation, and behaves more actively in RLWL/POPE-CPNT, shuttling rapidly and frequently between an α-plane zone and an adjacent midplane region. Results of this work reveal that a charge at the end of an ionic channel may significantly alter the transport characteristics of the channel.

Molecular Dynamics Simulations of Micelle Properties and Behaviors of Sodium Lauryl Ether Sulfate Penetrating Ceramide and Phospholipid Bilayers

Molecular dynamics (MD) simulations with the umbrella sampling (US) method were used to investigate the properties of micelles formed by sodium lauryl ether sulfate with two ether groups (SLE2S) and behaviors of corresponding surfactants transferring from micelles to ceramide and DMPC bilayer surfaces. Average micelle radii based on the Einstein-Smoluchowski and Stokes-Einstein relations showed excellent agreement with those measured by dynamic light scattering, while those obtained by evaluating the gyration radius or calculating the distance between the micelle sulfur atoms and center of mass overestimate the radii. As an SLE2S micelle was pulled down to the ceramide bilayer surface in a 400 ns constant-force steered MD (cf-SMD) simulation, the micelle was partially deformed on the bilayer surface, and several SLE2S surfactants easily were partitioned from the micelle into the ceramide bilayer. In contrast, a micelle was not deformed on the DMPC bilayer surface, and SLE2S surfactants were not transferred from the micelle to the DMPC bilayer. Potential of mean force (PMF) calculations revealed that the Gibbs free energy required for an SLE2S surfactant monomer to transfer from a micelle to bulk water can be compensated by decreased Gibbs free energy when an SLE2S monomer transfers into the ceramide bilayer from bulk water. In addition, micelle deformation on the ceramide bilayer surface can reduce the Gibbs free energy barrier required for a surfactant to escape the micelle and help the surfactant partition from the micelle into the ceramide bilayer. An SLE2S surfactant partitioning into the ceramide bilayer is attributed to hydrogen bonding and favorable interactions between the hydrophilic surfactant head and ceramide molecules, which are more dominant than the dehydration penalty during bilayer insertion. Such interactions between surfactant and lipid molecule heads are considerably reduced in DMPC bilayers owing to dielectric screening by water molecules deep inside the head/tail boundary between the DMPC bilayer. This computational work demonstrates the distinct behavior of SLE2S surfactant micelles on ceramide and DMPC bilayer surfaces in terms of variation in Gibbs free energy, which offers insight into designing surfactants used in transdermal drug delivery systems and cosmetics.

Separation of chloroform from a dilute solution using a cyclic peptide nanotube: A molecular dynamics study

This work firstly explored the potential application of a cyclic peptide nanotube (CPNT) in the separation of chloroform from a dilute solution. Four hydrophobic CPNTs of 8 × (WL)_4,5 and 8 × (AL)_4,5 all exhibit excellent adsorption characteristics to chloroform. The CPNT diameter, side chain structures and the concentration of chloroform in a solution all affect the adsorption characteristics of chloroform. CHCl₃ molecules are overwhelmingly adsorbed on the surfaces of these CPNTs as a cluster, and sporadically reside inside the channels, consistent with the chloroform's potentials of mean force (PMFs) inside and outside the channels. The distribution characteristics, molecular orientations and interactions with the surroundings of chloroform inside and outside four CPNTs embedded in individual dilute CHCl₃/water solutions were analyzed in detail, providing referable information of the adsorption characteristics of a hydrophobic CPNT to chloroform.

Conformational changes during permeation of Na+ through a modified cyclic peptide nanotube promote energy landscape roughness

Using a metadynamics approach, we investigate the potential of mean force for Na⁺ permeation inside a cyclic peptide nanotube (CPN) with modified interior as a function of ion position, coordination number, and lumen chemistry. We show that functionalizing the lumen of a CPN with a methyl-benzoic acid group introduces non-periodic variations in the internal energy of the nanotube, which dictate the overall free energy roughness during the permeation of Na⁺. These non-periodic variations arise from the structural dynamics of the functional group, where changes in the dihedral angles induced by the proximity of the ion give rise to conformational changes that increase landscape roughness and thereby decrease transport rate. Our computational framework emphasizes the advantages of using the coordination number as a collective variable to investigate the available conformations during ion permeation through CPNs, and reveals new structure-function relations for chemically tunable CPNs, paving the way for rational design of nano-porous systems with tunable selectivity and flux.

Transport behavior of a single Ca(2+), K(+), and Na(+) in a water-filled transmembrane cyclic peptide nanotube

Molecular dynamics simulations have been performed to investigate the transport properties of a single Ca(2+), K(+), and Na(+) in a water-filled transmembrane cyclic peptide nanotube (CPNT). Two transmembrane CPNTs, i.e., 8×(WL)n=4,5/POPE (with uniform lengths but various radii), were applied to clarify the dependence of ionic transport properties on the channel radius. A huge energy barrier keeps Ca(2+) out of the octa-CPNT, while Na(+) and K(+) can be trapped in two CPNTs. The dominant electrostatic interaction of a cation with water molecules leads to a high distribution of channel water around the cation and D-defects in the first and last gaps, and significantly reduces the axial diffusion of channel water. Water-bridged interactions were mostly found between the artificially introduced Ca(2+) and the framework of the octa-CPNT, and direct coordinations with the tube wall mostly occur for K(+) in the octa-CPNT. A cation may drift rapidly or behave lazily in a CPNT. K(+) behaves most actively and can visit the whole deca-CPNT quickly. The first solvation shells of Ca(2+) and Na(+) are basically saturated in two CPNTs, while the hydration of K(+) is incomplete in the octa-CPNT. The solvation structure of Ca(2+) in the octa-CPNT is most stable, while that of K(+) in the deca-CPNT is most labile. Increasing the channel radius induces numerous interchange attempts between the first-shell water molecules of a cation and the ones in the outer region, especially for the K(+) system.