Free energy calculations of gramicidin dimer dissociation

Design of peptides with strong binding affinity to poly(methyl methacrylate) resin by use of molecular simulation-based materials informatics

Peptides with strong binding affinities for poly(methyl methacrylate) (PMMA) resin were designed by use of materials informatics technology based on molecular dynamics simulation for the purpose of covering the resin surface with adhesive peptides, which were expected to result in eco-friendly and biocompatible biomaterials. From the results of binding affinity obtained with this molecular simulation, it was confirmed that experimental values could be predicted with errors <10%. By analyzing the simulation data with the response-surface method, we found that three peptides (RWWRPWW, EWWRPWR, and RWWRPWR), which consist of arginine (R), tryptophan (W), and proline (P), have strong binding affinity to the PMMA resin. These amino acids were effective because arginine and tryptophan have strong binding affinities for methoxycarbonyl groups and methyl groups, which are the main constituents of the PMMA resin, and proline stabilizes the flat zigzag structures of the peptides in water. The strong binding affinities of the three peptides were confirmed by experiments (surface plasmon resonance methods).

Time-resolved measurements of an ion channel conformational change driven by a membrane phase transition

Using temperature-jump infrared spectroscopy, we are able to trigger a gel-to-fluid phase transition in lipid vesicles and monitor in real time how a membrane protein responds to structural changes in the membrane. The melting of lipid domains in 1,2-dimyristoyl-sn-glycero-3-phosphocholine vesicles is observed to occur in as fast as 50 ns, with a temperature dependence characteristic of critical slowing. Gramicidin D (gD) added to the membrane responds primarily to the change in thickness of the membrane on a timescale coincident with the membrane melting. Using structure-based spectral modeling, we assign the conformational changes to compression and rotation of a partially dissociated gD dimer. Free energy calculations indicate that the high rate is a result of near-barrierless diffusion on a protein energy landscape that is radically reshaped by membrane thinning. The structural changes associated with the phase transition are similar to the fluctuation modes of fluid phase membranes, highlighting the importance of understanding the dynamic nature of the membrane environment around proteins.

Gramicidin A Channel Formation Induces Local Lipid Redistribution I: Experiment and Simulation

Integral membrane protein function can be modulated by the host bilayer. Because biological membranes are diverse and nonuniform, we explore the consequences of lipid diversity using gramicidin A channels embedded in phosphatidylcholine (PC) bilayers composed of equimolar mixtures of di-oleoyl-PC and di-erucoyl-PC (dC_18:1+dC_22:1, respectively), di-palmitoleoyl-PC and di-nervonoyl-PC (dC_16:1+dC_24:1, respectively), and di-eicosenoyl-PC (pure dC_20:1), all of which have the same average bilayer chain length. Single-channel lifetime experiments, molecular dynamics simulations, and a simple lipid compression model are used in tandem to gain insight into lipid redistribution around the channel, which partially alleviates the bilayer deformation energy associated with channel formation. The average single-channel lifetimes in the two-component bilayers (95 ± 10 ms for dC_18:1+dC_22:1 and 195 ± 20 ms for dC_16:1+dC_24:1) were increased relative to the single-component dC_20:1 control bilayer (65 ± 10 ms), implying lipid redistribution. Using a theoretical treatment of thickness-dependent changes in channel lifetimes, the effective local enrichment of lipids around the channel was estimated to be 58 ± 4% dC_18:1 and 66 ± 2% dC_16:1 in the dC_18:1+dC_22:1 and dC_16:1+dC_24:1 bilayers, respectively. 3.5-μs molecular dynamics simulations show 66 ± 2% dC_16:1 in the first lipid shell around the channel in the dC_16:1+dC_24:1 bilayer, but no significant redistribution (50 ± 4% dC_18:1) in the dC_18:1+dC_22:1 bilayer; these simulated values are within the 95% confidence intervals of the experimental averages. The strong preference for the better matching lipid (dC_16:1) near the channel in the dC_16:1+dC_24:1 mixture and lesser redistribution in the dC_18:1+dC_22:1 mixture can be explained by the energetic cost associated with compressing the lipids to match the channel's hydrophobic length.

Kinetics of peptide folding in lipid membranes

Despite our extensive understanding of water-soluble protein folding kinetics, much less is known about the folding dynamics and mechanisms of membrane proteins. However, recent studies have shown that for relatively simple systems, such as peptides that form a transmembrane α-helix, helical dimer, or helix-turn-helix, it is possible to assess the kinetics of several important steps, including peptide binding to the membrane from aqueous solution, peptide folding on the membrane surface, helix insertion into the membrane, and helix-helix association inside the membrane. Herein, we provide a brief review of these studies and also suggest new initiation and probing methods that could lead to improved temporal and structural resolution in future experiments.