Free-Energy Analysis of Peptide Binding in Lipid Membrane Using All-Atom Molecular Dynamics Simulation Combined with Theory of Solutions

Solvation dynamics on the diffusion timescale elucidated using energy-represented dynamics theory

Photoexcitation of a solute alters the solute-solvent interaction, resulting in the nonequilibrium relaxation of the solvation structure, often called a dynamic Stokes shift or solvation dynamics. Thanks to the local nature of the solute-solvent interaction, the characteristics of the local solvent environment dissolving the solute can be captured by the observation of this process. Recently, we derived the energy-represented Smoluchowski-Vlasov (ERSV) equation, a diffusion equation for molecular liquids, which can be used to analyze the solvation dynamics on the diffusion timescale. This equation expresses the time development for the solvent distribution on the solute-solvent pair interaction energy (energy coordinate). Since the energy coordinate can effectively treat the solvent flexibility in addition to the position and orientation, the ERSV equation can be utilized in various solvent systems. Here, we apply the ERSV equation to the solvation dynamics of 6-propionyl-2-dimethylamino naphthalene (Prodan) in water and different alcohol solvents (methanol, ethanol, and 1-propanol) for clarifying the differences of the relaxation processes among these solvents. Prodan is a solvent-sensitive fluorescent probe and is thus widely utilized for investigating heterogeneous environments. On the long timescale, the ERSV equation satisfactorily reproduces the relaxation time correlation functions obtained from the molecular dynamics (MD) simulations for these solvents. We reveal that the relaxation time coefficient on the diffusion timescale linearly correlates with the inverse of the translational diffusion coefficients for the alcohol solvents because of the Prodan-solvent energy distributions among the alcohols. In the case of water, the time coefficient deviates from the linear relationship for the alcohols due to the difference in the extent of importance of the collective motion between the water and alcohol solvents.

Diffusion theory of molecular liquids in the energy representation and application to solvation dynamics

The generalized Langevin equation (GLE) formalism is a useful theoretical fundament for analyzing dynamical phenomena rigorously. Despite the systematic formulation of dynamics theories with practical approximations, however, the applicability of GLE-based methods is still limited to simple polyatomic liquids due to the approximate treatment of molecular orientations involved in the static molecular liquid theory. Here, we propose an exact framework of dynamics based on the GLE formalism incorporating the energy representation theory of solution, an alternative static molecular liquid theory. A fundamental idea is the projection of the relative positions and orientations of solvents around a solute onto the solute-solvent interaction, namely the energy coordinate, enabling us to describe the dynamics on a one-dimensional coordinate. Introducing systematic approximations, such as the overdamped limit, leads to the molecular diffusion equation in the energy representation that is described in terms of the distribution function of solvents on the energy coordinate and the diffusion coefficients. The present theory is applied to the solvation dynamics triggered by the photoexcitation of benzonitrile. The long-time behavior of the solvation time correlation function is in good agreement with that obtained by the molecular dynamics simulation.

Adsorption Energetics of Amino Acid Analogs on Polymer/Water Interfaces Studied by All-Atom Molecular Dynamics Simulation and a Theory of Solutions

Energetics of adsorption was addressed with all-atom molecular dynamics simulation on the interfaces of poly(2-methoxyethyl acrylate) (PMEA), poly(methyl methacrylate) (PMMA), and poly(butyl acrylate) (PBA) with water. A wide variety of adsorbate solutes were examined, and the free energy of adsorption was computed with the method of energy representation. It was found that the adsorption free energy was favorable (negative) for all the combinations of solute and polymer, and among PMEA, PMMA, and PBA, the strongest adsorption was observed on PMMA for the hydrophobic solutes and on PMEA for the hydrophilic ones. According to the decomposition of the adsorption free energy into the contributions from polymer and water, it was seen that the polymer contribution is larger in magnitude with the solute size. The total free energy of adsorption was correlated well with the solvation free energy in bulk water only for hydrophobic solutes. The roles of the intermolecular interaction components such as electrostatic, van der Waals, and excluded-volume were further studied, and the electrostatic component was influential only in determining the polymer dependences of the adsorption propensities of hydrophilic solutes. The extent of adsorption was shown to be ranked by the van der Waals component in the solute-polymer interaction separately over the hydrophilic and hydrophobic solutes, with the excluded-volume effect from water pointed out to also drive the adsorption.

The effects of molecular weight and orientation on the membrane permeation and partitioning of polycyclic aromatic hydrocarbons: a computational study

Membrane permeation and the partitioning of polycyclic aromatic hydrocarbons (PAHs) are crucial aspects affecting their carcinogenicity and mutagenicity. However, a clear understanding of these processes is still rare due to the difficulty of determining the details experimentally. Here, the interactions between PAHs and lipid bilayers were studied by molecular simulations, mainly to check the influence of molecular weight and orientation. The liposome-water partition coefficient (K_LW), transmembrane time (τ), and permeability coefficient (P) of the PAHs were calculated by integrating free energy profiles from umbrella sampling. For selected PAHs, the membrane adsorption is a spontaneous process. The preferred location is near the CC bond and the orientation is related to the molecular structure. The P values of all the PAHs are basically the same order of magnitude, which means that the molecular weight contributes little to the process. As for K_LW and τ, they show obvious increases with different molecular weights. Unconstrained simulations showed that a flat orientation on the membrane surface would prevent PAHs from being transported through the membrane. Highly hydrophobic driving forces are not always good for the absorption of PAHs, especially the formation of aggregates. In addition, the orientations and energetic barriers of PAHs near the midplane of the lipid bilayer explain the different transitions of high- and low-weight PAHs. This work provides molecular level details relating to the interactions of PAHs with lipid membranes, with significance for understanding the health effects of PAHs.

Solvation energetics of proteins and their aggregates analyzed by all-atom molecular dynamics simulations and the energy-representation theory of solvation

Solvation is a controlling factor for the structure and function of proteins. This article addresses the effects of solvation from an energetic perspective for the fluctuations and cosolvent-induced changes in protein structures and the equilibrium of aggregate formation for a peptide. A theoretical framework to analyze the solvation effects with an explicit solvent is introduced by adopting the energy-representation theory of solvation, and the connection of the solvation free energy to the protein structure and the aggregation tendency is quantitatively described in combination with all-atom molecular dynamics simulations. The interaction components that govern the solvation effects on the structural variations of proteins are further identified through correlation analysis, and a computational scheme to assess the shift of an aggregation equilibrium due to the addition of a cosolvent is provided.

Advances in Molecular Understanding of α-Helical Membrane-Active Peptides

Biological membranes separate the interior of cells or cellular compartments from their outer environments. This barrier function of membranes can be disrupted by membrane-active peptides, some of which can spontaneously penetrate through the membranes or open leaky transmembrane pores. However, the origin of their activity/toxicity is not sufficiently understood for the development of more potent peptides. To this day, there are no design rules that would be generally valid, and the role of individual amino acids tends to be sequence-specific.In this Account, we describe recent progress in understanding the design principles that govern the activity of membrane-active peptides. We focus on α-helical amphiphilic peptides and their ability to (1) translocate across phospholipid bilayers, (2) form transmembrane pores, or (3) act synergistically, i.e., to produce a significantly more potent effect in a mixture than the individual components.We refined the description of peptide translocation using computer simulations and demonstrated the effect of selected residues. Our simulations showed the necessity to explicitly include charged residues in the translocation description to correctly sample the membrane perturbations they can cause. Using this description, we calculated the translocation of helical peptides with and without the kink induced by the proline/glycine residue. The presence of the kink had no effect on the translocation barrier, but it decreased the peptide affinity to the membrane and reduced the peptide stability inside the membrane. Interestingly, the effects were mainly caused by the peptide's increased polarity, not the higher flexibility of the kink.Flexibility plays a crucial role in pore formation and affects distinct pore structures in different ways. The presence of a kink destabilizes barrel-stave pores, because the kink prevents the tight packing of peptides in the bundle, which is characteristic of the barrel-stave structure. In contrast, the kink facilitates the formation of toroidal pores, where the peptides are only loosely arranged and do not need to closely assemble. The exact position of the kink in the sequence further determines the preferred arrangement of peptides in the pore, i.e., an hourglass or U-shaped structure. In addition, we demonstrated that two self-associated (via termini) helical peptides could mimic the behavior of peptides with a helix-kink-helix motif.Finally, we review the recent findings on the peptide synergism of the archetypal mixture of Magainin 2 and PGLa peptides. We focused on a bacterial plasma membrane mimic that contains negatively charged lipids and lipids with negative intrinsic curvature. We showed that the synergistic action of peptides was highly dependent on the lipid composition. When the lipid composition and peptide/lipid ratios were changed, the systems exhibited more complex behavior than just the previously reported pore formation. We observed membrane adhesion, fusion, and even the formation of the sponge phase in this regime. Furthermore, enhanced adhesion/partitioning to the membrane was reported to be caused by lipid-induced peptide aggregation.In conclusion, the provided molecular insight into the complex behavior of membrane-active peptides provides clues for the design and modification of antimicrobial peptides or toxins.

Structure Changes of a Membrane Polypeptide under an Applied Voltage Observed with Surface-Enhanced 2D IR Spectroscopy

The structures of many membrane-bound proteins and polypeptides depend on the membrane potential. However, spectroscopically studying their structures under an applied field is challenging, because a potential is difficult to generate across more than a few bilayers. We study the voltage-dependent structures of the membrane-bound polypeptide, alamethicin, using a spectroelectrochemical cell coated with a rough, gold film to create surface plasmons. The plasmons sufficiently enhance the 2D IR signal to measure a single bilayer. The film is also thick enough to conduct current and thereby apply a potential. The 2D IR spectra resolve features from both 3₁₀- and α-helical structures and cross-peaks connecting the two. We observe changes in the peak intensity, not their frequencies, upon applying a voltage. A similar change occurs with pH, which is known to alter the angle of alamethicin relative to the surface normal. The spectra are modeled using a vibrational exciton Hamiltonian, and the voltage-dependent spectra are consistent with a change in angle of the 3₁₀- and α-helices in the membrane from 55 to 44°and from 31 to 60°, respectively. The 3₁₀- and α-helices are coupled by approximately 10 cm^-1. These experiments provide new structural information about alamethicin under a potential difference and demonstrate a technique that might be applied to voltage-gated membrane proteins and compared to molecular dynamics structures.

Identification and antioxidant activity of bovine bone collagen-derived novel peptides prepared by recombinant collagenase from Bacillus cereus

Millions of tons of collagen-rich bovine bone are produced as byproducts of the consumption of beef. Hydrolyzing bovine bone collagen (BBC) is an effective measure for both increasing its added value and protecting the environment. In this study, a kind of recombinant bacterial collagenase mining from Bacillus cereus was successfully performed and applied to hydrolyze BBC to collagen-soluble peptides (CPP). Response surface methodology (RSM) was applied to optimize the processing conditions of antioxidant CPP, attaining a distinguished ABTS free radical scavenging activity of 99.21 ± 0.35% while keeping DPPH free radical scavenging activity and reducing power at high levels under the optimal condition. Furthermore, we identified five new antioxidant peptides by LC-MS/MS with typical collagen repeated Gly-Xaa-Yaa sequence units within the CPP. These results suggest that our recombinant collagenase is a powerful tool for degrading collagen and the CPP are promising candidates for antioxidant and related functional food applications.

Methodology for Further Thermostabilization of an Intrinsically Thermostable Membrane Protein Using Amino Acid Mutations with Its Original Function Being Retained

We develop a new methodology best suited to the identification of thermostabilizing mutations for an intrinsically stable membrane protein. The recently discovered thermophilic rhodopsin, whose apparent midpoint temperature of thermal denaturation T_m is measured to be ∼91.8 °C, is chosen as a paradigmatic target. In the methodology, we first regard the residues whose side chains are missing in the crystal structure of the wild type (WT) as the "residues with disordered side chains," which make no significant contributions to the stability, unlike the other essential residues. We then undertake mutating each of the residues with disordered side chains to another residue except Ala and Pro, and the resultant mutant structure is constructed by modifying only the local structure around the mutated residue. This construction is based on the postulation that the structure formed by the other essential residues, which is nearly optimized in such a highly stable protein, should not be modified. The stability changes arising from the mutations are then evaluated using our physics-based free-energy function (FEF). We choose the mutations for which the FEF is much lower than for the WT and test them by experiments. We successfully find three mutants that are significantly more stable than the WT. A double mutant whose T_m reaches ∼100 °C is also discovered.

Spatially-Decomposed Free Energy of Solvation Based on the Endpoint Density-Functional Method

A spatially resolved version of the density-functional method for solvation thermodynamics is presented by extending the free-energy functional previously established in the one-dimensional, energy representation and formulating a new expression in a mixed four-dimensional representation (three dimensions for position and one dimension for energy). The space was further divided into a set of discrete regions with respect to the relative position of a solvent molecule from the solute, and the spatially decomposed energetics of solvation were analyzed for small molecules with a methyl, amine, or hydroxyl group and alanine dipeptide in solvent water. It was observed that the density of the solvation free energy is weakly dependent on the solute site in the excluded-volume region and is distinctively favorable in the first shells of the solute atoms that can readily form hydrogen bonds with water. The solvent-reorganization term reduces faster with the separation from the solute than the direct interaction between the solute and solvent, and the latter governs the energetics in the second shell and outer regions. The sum of the contributions to the free energy from the excluded volume and first shell was found to deviate significantly from the total sum over all the regions, implying that the solvation free energy is not spatially localized near the solute in a quantitative sense. Still, a local description was shown to be valid as confirmed by the correlation of the total value of free energy with the corresponding value obtained by integrating the free-energy density to the second shell. The theoretical framework developed in the present work to spatially decompose the solvation free energy can thus be useful to identify stabilizing or destabilizing regions of solvent proximate to a solute and to analyze the role that the displacement of interfacial water plays in the thermodynamics of molecular association.

Energetics of cosolvent effect on peptide aggregation

The cosolvent effect on the equilibrium of peptide aggregation is reviewed from the energetic perspective. It is shown that the excess chemical potential is stationary against the variation of the distribution function for the configuration of a flexible solute species and that the derivative of the excess chemical potential with respect to the cosolvent concentration is determined by the corresponding derivative of the solvation free energy averaged over the solute configurations. The effect of a cosolvent at low concentrations on a chemical equilibrium can then be addressed in terms of the difference in the solvation free energy between pure-water solvent and the mixed solvent with the cosolvent, and illustrative analyses with all-atom model are presented for the aggregation of an 11-residue peptide by employing the energy-representation method to compute the solvation free energy. The solvation becomes more favorable with addition of the urea or DMSO cosolvent, and the extent of stabilization is smaller for larger aggregate. This implies that these cosolvents inhibit the formation of an aggregate, and the roles of such interaction components as the electrostatic, van der Waals, and excluded-volume are discussed.

Antimicrobial cell penetrating peptides with bacterial cell specificity: pharmacophore modelling, quantitative structure activity relationship and molecular dynamics simulation

Current research has shown cell-penetrating peptides and antimicrobial peptides (AMPs) as probable vectors for use in drug delivery and as novel antibiotics. It has been reported that the higher the therapeutic index (TI) the higher would be the bacterial cell penetrating ability. To the best of our knowledge, no in-silico study has been performed to determine bacterial cell specificity of the antimicrobial cell penetrating peptides (aCPP's) based on their TI. The aim of this study was to develop a quantitative structure activity relationship (QSAR) model, which can estimate antimicrobial potential and cell-penetrating ability of aCPPs against S. aureus, to confirm the relationship between the TI and aCPPs and to identify specific descriptors responsible for aCPPs penetrating ability. Molecular dynamics (MD) simulation was also performed to confirm the membrane insertion of the most active aCPPs obtained from the QSAR study. The most appropriate pharmacophore was identified to predict the aCPP's activity. The statistical results confirmed the validity of the model. The QSAR model was successful in identifying the optimal aCPP with high activity prediction and provided insights into the structural requirements to correlate their TI to cell penetrating ability. MD simulation of the best aCPP with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer confirmed its interaction with the membrane and the C-terminal residues of the aCPP played a key role in membrane penetration. The strategy of combining QSAR and molecular dynamics, allowed for optimal estimation of ligand-target interaction and confirmed the importance of Trp and Lys in interacting with the POPC bilayer. Communicated by Ramaswamy H. Sarma.

Optimal Hydrophobicity and Reorientation of Amphiphilic Peptides Translocating through Membrane

Cell-penetrating and some antimicrobial peptides can translocate across lipid bilayers without disrupting the membrane structure. However, the molecular properties required for efficient translocation are not fully understood. We employed the Metropolis Monte Carlo method together with coarse-grained models to systematically investigate free-energy landscapes associated with the translocation of secondary amphiphilic peptides. We studied α-helical peptides with different length, amphiphilicity, and distribution of hydrophobic content and found a common translocation path consisting of adsorption, tilting, and insertion. In the adsorbed state, the peptides are parallel to the membrane plane, whereas, in the inserted state, the peptides are perpendicular to the membrane. Our simulations demonstrate that, for all tested peptides, there is an optimal ratio of hydrophilic/hydrophobic content at which the peptides cross the membrane the easiest. Moreover, we show that the hydrophobicity of peptide termini has an important effect on the translocation barrier. These results provide general guidance to optimize peptides for use as carriers of molecular cargos or as therapeutics themselves.

Free-energy analysis of physisorption on solid-liquid interface with the solution theory in the energy representation

Physisorption of urea on its crystal in contact with water was subject to energetics analysis with all-atom molecular dynamics simulation. The transfer free energy of urea to an adsorption site was treated in the framework of the energy-representation theory of solutions, which allows a fast computation of the free energy in an inhomogeneous environment with solid-liquid interface. The preference of adsorption was then compared between the (001) and (110) faces, and it was found that the physisorption is more favorable on (001) than on (110) in correspondence to the hydrogen bonding between the adsorbed urea and the crystal urea. Among the terrace configurations of adsorption, the attractive interaction governs the preferable site with a minor role of the repulsive interaction. The effect of an edge was also treated by examining the terrace and step and was shown to be strongly operative on the (110) face when the CO group of the adsorbed urea points toward the edge. The present work demonstrates that the solution theory can be a framework for analyzing the energetics of physisorption and addressing the roles of the crystal and liquid at the interface through the systematic decomposition of free energy.