Tipping the Scale from Disorder to Alpha-helix: Folding of Amphiphilic Peptides in the Presence of Macroscopic and Molecular Interfaces

Double Stranded DNA Binding Stapled Peptides: An Emerging Tool for Transcriptional Regulation

Stapled peptides have rapidly established themselves as a powerful technique to mimic α-helical interactions with a short peptide sequence. There are many examples of stapled peptides that successfully disrupt α-helix-mediated protein-protein interactions, with an example currently in clinical trials. DNA-protein interactions are also often mediated by α-helices and are involved in all transcriptional regulation processes. Unlike DNA-binding small molecules, which typically lack DNA sequence selectivity, DNA-binding proteins bind with high affinity and high selectivity. These are ideal candidates for the design DNA-binding stapled peptides. Despite the parallel to protein-protein interaction disrupting stapled peptides and the need for sequence specific DNA binders, there are very few DNA-binding stapled peptides. In this review we examine all the known DNA-binding stapled peptides. Their design concepts are compared to stapled peptides that disrupt protein-protein interactions and based on the few examples in the literature, DNA-binding stapled peptide trends are discussed.

A class of peptides designed to replicate and enhance the Receptor for Hyaluronic Acid Mediated Motility binding domain

The extra-cellular matrix (ECM) is a complex and rich microenvironment that is exposed and over-expressed across several injury or disease pathologies. Biomaterial therapeutics are often enriched with peptide binders to target the ECM with greater specificity. Hyaluronic acid (HA) is a major component of the ECM, yet to date, few HA adherent peptides have been discovered. A class of HA binding peptides was designed using B(X₇)B hyaluronic acid binding domains inspired from the helical face of the Receptor for Hyaluronic Acid Mediated Motility (RHAMM). These peptides were bioengineered using a custom alpha helical net method, allowing for the enrichment of multiple B(X₇)B domains and the optimisation of contiguous and non-contiguous domain orientations. Unexpectedly, the molecules also exhibited the behaviour of nanofiber forming self-assembling peptides and were investigated for this characteristic. Ten 23-27 amino acid residue peptides were assessed. Simple molecular modelling was used to depict helical secondary structures. Binding assays were performed with varying concentrations (1-10 mg/mL) and extra-cellular matrices (HA, collagens I-IV, elastin, and Geltrex). Concentration mediated secondary structures were assessed using circular dichroism (CD), and higher order nanostructures were visualized using transmission electron microscopy (TEM). All peptides formed the initial apparent 3₁₀/alpha-helices, yet peptides 17x-3, 4, BHP3 and BHP4 were HA specific and potent (i.e., a significant effect) binders at increasing concentrations. These peptides shifted from apparent 3₁₀/alpha-helical structures at low concentration to beta-sheets at increasing concentration and also formed nanofibers which are noted as self-assembling structures. Several of the HA binding peptides outperformed our positive control (mPEP35) at 3-4 times higher concentrations, and were enhanced by self-assembly as each of these groups had observable nanofibers. STATEMENT OF SIGNIFICANCE: Specific biomolecules or peptides have played a crucial role in developing materials or systems to deliver key drugs and therapeutics to a broad spectrum of diseases and disorders. In these diseased tissues, cells build protein/sugar networks, which are uniquely exposed and great targets to deliver drugs to. Hyaluronic acid (HA) is involved in every stage of injury and is abundant in cancer. To date, only two HA specific peptides have been discovered. In our work, we have designed a way to model and trace binding regions as they appear on the face of a helical peptide. Using this method we have created a family of peptides enriched with HA binding domains that stick with 3-4 higher affinity than those previously discovered.

Understanding the helical stability of charged peptides

Cationic helical peptides play a crucial role in applications such as anti-microbial and anticancer activity. The activity of these peptides directly correlates with their helicity. In this study, we have performed extensive all-atom molecular dynamics simulations of 25 Lysine-Leucine co-polypeptide sequences of varying charge density ( λ ) and patterns. Our findings showed that, an increase in the charge density on the peptide leads to a gradual decrease in the helicity up to a critical charge density λ c . Beyond λ c _, a complete helix to coil transition was observed. The decrease in the helicity is correlated with the increased number of water molecules in first solvation shell, solvent-exposed surface area, and a higher value of the radius of gyration of the peptide.

Functional and Conformational Plasticity of an Animal Group 1 LEA Protein

Group 1 (Dur-19, PF00477, LEA_5) Late Embryogenesis Abundant (LEA) proteins are present in organisms from all three domains of life, Archaea, Bacteria, and Eukarya. Surprisingly, Artemia is the only genus known to include animals that express group 1 LEA proteins in their desiccation-tolerant life-history stages. Bioinformatics analysis of circular dichroism data indicates that the group 1 LEA protein AfLEA1 is surprisingly ordered in the hydrated state and undergoes during desiccation one of the most pronounced disorder-to-order transitions described for LEA proteins from A. franciscana. The secondary structure in the hydrated state is dominated by random coils (42%) and β-sheets (35%) but converts to predominately α-helices (85%) when desiccated. Interestingly, AfLEA1 interacts with other proteins and nucleic acids, and RNA promotes liquid–liquid phase separation (LLPS) of the protein from the solvent during dehydration in vitro. Furthermore, AfLEA1 protects the enzyme lactate dehydrogenase (LDH) during desiccation but does not aid in restoring LDH activity after desiccation-induced inactivation. Ectopically expressed in D. melanogaster Kc167 cells, AfLEA1 localizes predominantly to the cytosol and increases the cytosolic viscosity during desiccation compared to untransfected control cells. Furthermore, the protein formed small biomolecular condensates in the cytoplasm of about 38% of Kc167 cells. These findings provide additional evidence for the hypothesis that the formation of biomolecular condensates to promote water stress tolerance during anhydrobiosis may be a shared feature across several groups of LEA proteins that display LLPS behaviors.

Inside and Out of the Pore: Comparing Interactions and Molecular Dynamics of Influenza A M2 Viroporin Complexes in Standard Lipid Bilayers

Ion channels located at viral envelopes (viroporins) have a critical function for the replication of infectious viruses and are important drug targets. Over the last decade, the number and duration of molecular dynamics (MD) simulations of the influenza A M2 ion channel owing to the increased computational efficiency. Here, we aimed to define the system setup and simulation conditions for the correct description of the protein-pore and the protein-lipid interactions for influenza A M2 in comparison with experimental data. We performed numerous MD simulations of the influenza A M2 protein in complex with adamantane blockers in standard lipid bilayers using OPLS2005 and CHARMM36 (C36) force fields. We explored the effect of varying the M2 construct (M2(22-46) and M2(22-62)), the lipid buffer size and type (stiffer DMPC or softer POPC with or without 20% cholesterol), the simulation time, the H37 protonation site (N_δ or Ν_ε), the conformational state of the W41 channel gate, and M2's cholesterol binding sites (BSs). We report that the 200 ns MD with M2(22-62) (having N_ε Η37) in the 20 Å lipid buffer with the C36 force field accurately describe: (a) the M2 pore structure and interactions inside the pore, that is, adamantane channel blocker location, water clathrate structure, and water or chloride anion blockage/passage from the M2 pore in the presence of a channel blocker and (b) interactions between M2 and the membrane environment as reflected by the calculation of the M2 bundle tilt, folding of amphipathic helices, and cholesterol BSs. Strikingly, we also observed that the C36 1 μs MD simulations using M2(22-62) embedded in a 20 Å POPC:cholesterol (5:1) scrambled membrane produced frequent interactions with cholesterol, which when combined with computational kinetic analysis, revealed the experimentally observed BSs of cholesterol and suggested three similarly long-interacting positions in the top leaflet that have previously not been observed experimentally. These findings promise to be useful for other viroporin systems.

Liposome fusion with orthogonal coiled coil peptides as fusogens: the efficacy of roleplaying peptides

Biological membrane fusion is a highly specific and coordinated process as a multitude of vesicular fusion events proceed simultaneously in a complex environment with minimal off-target delivery. In this study, we develop a liposomal fusion model system with specific recognition using lipidated derivatives of a set of four de novo designed heterodimeric coiled coil (CC) peptide pairs. Content mixing was only obtained between liposomes functionalized with complementary peptides, demonstrating both fusogenic activity of CC peptides and the specificity of this model system. The diverse peptide fusogens revealed important relationships between the fusogenic efficacy and the peptide characteristics. The fusion efficiency increased from 20% to 70% as affinity between complementary peptides decreased, (from K_F ≈ 10⁸ to 10⁴ M⁻¹), and fusion efficiency also increased due to more pronounced asymmetric role-playing of membrane interacting ‘K’ peptides and homodimer-forming ‘E’ peptides. Furthermore, a new and highly fusogenic CC pair (E₃/P1_K) was discovered, providing an orthogonal peptide triad with the fusogenic CC pairs P2_E/P2_K and P3_E/P3_K. This E₃/P1_k pair was revealed, via molecular dynamics simulations, to have a shifted heptad repeat that can accommodate mismatched asparagine residues. These results will have broad implications not only for the fundamental understanding of CC design and how asparagine residues can be accommodated within the hydrophobic core, but also for drug delivery systems by revealing the necessary interplay of efficient peptide fusogens and enabling the targeted delivery of different carrier vesicles at various peptide-functionalized locations.

We developed a liposomal fusion model system with specific recognition using a set of heterodimeric coiled coil peptide pairs. This study unravels important structure–fusogenic efficacy relationships of peptide fusogens.

Intramolecular interactions play key role in stabilization of pHLIP at acidic conditions

The pH-Low Insertion Peptide (pHLIP) is a membrane-active peptide that spontaneously folds into a transmembrane α-helix upon acidification. This activity enables pHLIP to potentially act as a vector for drugs related to diseases characterized by acidosis such as cancer or heart ischemia. Presently, due to aggregation-based effects, formulations of pHLIP are only viable at near-μM concentrations. In addition, since most of pHLIP's measurable qualities involve a membrane, probing the details of pHLIP in the interstitial region is difficult. In attempts to shed light on these issues, we performed constant pH molecular dynamics simulations on pHLIP as well as P20G, a variant with increased helicity, in solution at 0 and 150 mM NaCl over a broad range of pHs. In general, the addition of ions reduced the effective pK_a of the acidic residues in pHLIP. P20G exhibits a higher helicity than pHLIP in general and is more compact than pHLIP at pH values under 4. In terms of charge effects, sodium cations localized predominantly to the C-terminus of the peptide with a high density of acidic residues. Additionally, the salt bridge between R11 and D14 is by far the most favored and particularly so with pHLIP at 150 mM NaCl. We expect that this approach will be a valuable tool to screen variants of pHLIP for favorable properties in solution, an aspect of pHLIP design that to this point has largely been neglected.

Location and Conformation of the LKα14 Peptide in Water/Ethanol Mixtures

It is widely recognized that solvation is one of the major factors determining structure and functionality of proteins and long peptides, however it is a formidable challenge to address it both experimentally and computationally. For this reason, simple peptides are used to study fundamental aspects of solvation. It is well established that alcohols can change the peptide conformation and tuning of the alcohol content in solution can dramatically affect folding and, as a consequence, the function of the peptide. In this work, we focus on the leucine and lysine based LKα14 peptide designed to adopt an α-helical conformation at an apolar-polar interface. We investigate LKα14 peptide's bulk and interfacial behavior in water/ethanol mixtures combining a suite of experimental techniques (namely, circular dichroism and nuclear magnetic resonance spectroscopy for the bulk solution, surface pressure measurements and vibrational sum frequency generation spectroscopy for the air-solution interface) with molecular dynamics simulations. We observe that ethanol highly affects both the peptide location and conformation. At low ethanol content LKα14 lacks a clear secondary structure in bulk and shows a clear preference to reside at the air-solution interface. When the ethanol content in solution increases, the peptide's interfacial affinity is markedly reduced and the peptide approaches a stable α-helical conformation in bulk facilitated by the amphiphilic nature of the ethanol molecules.

Computational studies of membrane pore formation induced by Pin2

Understanding, at the molecular level, the effect of AMPs on biological membranes is of crucial importance given the increasing number of multidrug-resistant bacteria. Being part of an ancient type of innate immunity system, AMPs have emerged as a potential solution for which bacteria have not developed resistance. Traditional antibiotics specifically act on biosynthetic pathways, while AMPs may directly destabilize the lipid membrane, but it is unclear how AMPs affect the membrane's stability. We performed multiscale molecular dynamics simulations to investigate the structural features leading to membrane pores formation on zwitterionic and anionic membranes by the antimicrobial peptide (AMP) Pandinin 2 (Pin2). Some experimental reports propose that Pin2 could form barrel-stave pores, while others suggest that it could form toroidal pores. Since there is no conclusive evidence of which type of pore is formed by Pin2 on bilayers, performing molecular dynamics simulations on these systems could shed some light on whether or not or what type of pore Pin2 forms on model membranes. Our results are focused on a detailed description of the pore formation by Pin2 in POPC and POPE:POPG membranes., which strongly suggest that Pin2 forms a toroidal pore and not a barrel-shaped pore; this type of pore also affects the membrane properties. In the process, a phospholipid remodeling in the POPE:POPG membrane takes place. Moreover, the pores formed by Pin2 indicate that they are selective for the chlorine ion. There are no previous ion selectivity reports for other AMPs with similar physicochemical properties, such as melittin and magainin.Communicated by Ramaswamy H. Sarma.

Effect of the air-water interface on the conformation of amyloid beta

It has long been recognized that liquid interfaces, such as the air-water interface (AWI), can enhance the formation of protein fibrils. This makes liquid interfaces attractive templates for fibril formation but fully realizing this requires knowledge of protein behavior at interfaces, which is currently lacking. To address this, molecular dynamics simulation is used to investigate fragments of amyloid beta, a model fibril forming protein, at the air-water interface. At the air-water interface, the enrichment of aggregation-prone helical conformations provides a mechanism for the enhancement of fibrillation at interfaces. The conformational ensemble at the air-water interface was also considerably reduced compared to bulk solution due to the tendency of hydrophobic side chains partitioning into the air restricting the range of conformations. Little overlap between the conformational ensembles at the AWI and in the bulk solution was found, suggesting that AWI induces the formation of a different set of structures compared to bulk solution. The smaller Aβ(16-22) and Aβ(25-35) fragments show an increase in the propensity for an ordered secondary structure at the air-water interface but with a increased propensity for turn over other motifs, illustrating the importance of intra-protein interactions for stabilizing helical and extended conformations.

Effect of surface chemistry on islet amyloid polypeptide conformation

The formation of dense, linear arrays (fibrils) by biomolecules is the hallmark of a number of degenerative diseases, such as Alzheimer's and type-2 diabetes. Protein fibrils have also attracted interest as building blocks for new materials. It has long been recognized that surfaces can affect the fibrillation process. Recent work on the model fibril forming protein human islet amyloid polypeptide (hIAPP) has shown that while the protein concentration is highest at hydrophobic surfaces, the rate of fibril formation is lower than on other surfaces. To understand this, replica exchange molecular dynamics simulations were used to investigate the conformations that hIAPP adopts on surfaces of different hydrophobicities. The hydrophobic surface stabilizes α-helical structures which are significantly different to those found on the hydrophilic surface and in bulk solution. There is also a greatly reduced conformational ensemble on the hydrophobic surface due to long-lived contacts between hydrophobic residues on the protein and the surface. This new microscopic information will help us determine the mechanism of the enhancement of fibril formation on surfaces and provides new insight into the effect of nanointerfaces and protein conformation.

Membrane Lipids Assist Catalysis by CTP: Phosphocholine Cytidylyltransferase

While most of the articles in this issue review the workings of integral membrane enzymes, in this review, we describe the catalytic mechanism of an enzyme that contains a soluble catalytic domain but appears to catalyze its reaction on the membrane surface, anchored and assisted by a separate regulatory amphipathic helical domain and inter-domain linker. Membrane partitioning of CTP: phosphocholine cytidylyltransferase (CCT), a key regulatory enzyme of phosphatidylcholine metabolism, is regulated chiefly by changes in membrane phospholipid composition, and boosts the enzyme's catalytic efficiency >200-fold. Catalytic enhancement by membrane binding involves the displacement of an auto-inhibitory helix from the active site entrance-way and promotion of a new conformational ensemble for the inter-domain, allosteric linker that has an active role in the catalytic cycle. We describe the evidence for close contact between membrane lipid, a compact allosteric linker, and the CCT active site, and discuss potential ways that this interaction enhances catalysis.

Combining Molecular and Spin Dynamics Simulations with Solid-State NMR: A Case Study of Amphiphilic Lysine-Leucine Repeat Peptide Aggregates

Interpreting dynamics in solid-state molecular systems requires characterization of the potentially heterogeneous environmental contexts of molecules. In particular, the analysis of solid-state nuclear magnetic resonance (ssNMR) data to elucidate molecular dynamics (MD) involves modeling the restriction to overall tumbling by neighbors, as well as the concentrations of water and buffer. In this exploration of the factors that influence motion, we utilize atomistic MD trajectories of peptide aggregates with varying hydration to mimic an amorphous solid-state environment and predict ssNMR relaxation rates. We also account for spin diffusion in multiply spin-labeled (up to 19 nuclei) residues, with several models of dipolar-coupling networks. The framework serves as a general approach to determine essential spin couplings affecting relaxation, benchmark MD force fields, and reveal the hydration dependence of dynamics in a crowded environment. We demonstrate the methodology on a previously characterized amphiphilic 14-residue lysine-leucine repeat peptide, LKα14 (Ac-LKKLLKLLKKLLKL-c), which has an α-helical secondary structure and putatively forms leucine-burying tetramers in the solid state. We measure the R₁ relaxation rates of uniformly ¹³C-labeled and site-specific ²H-labeled leucines in the hydrophobic core of LKα14 at multiple hydration levels. Studies of 9 and 18 tetramer bundles reveal the following: (a) for the incoherent component of ¹³C relaxation, the nearest-neighbor spin interactions dominate, while the ¹H-¹H interactions have minimal impact; (b) the AMBER ff14SB dihedral barriers for the leucine C_γ-C_δ bond ("methyl rotation barriers") must be lowered by a factor of 0.7 to better match the ²H data; (c) proton-driven spin diffusion explains some of the discrepancy between experimental and simulated rates for the C_β and C_α nuclei; and (d) ¹³C relaxation rates are mostly underestimated in the MD simulations at all hydrations, and the discrepancies identify likely motions missing in the 50 ns MD trajectories.

Diverse Bilayer Morphologies Achieved via α-Helix-to-β-Sheet Transitions in a Short Amphiphilic Peptide

Transmembrane proteins are functional macromolecules that direct the flow of small molecules and ions across a lipid bilayer. Here, we propose the development of helical peptide amphiphiles that will serve as both the bilayer and the functional unit of a self-assembled peptide bilayer membrane. The peptide, K₃L₁₂, was designed not only to possess dimensions similar to that of a lipid bilayer but also to yield a structurally robust, α-helical bilayer. The formation of α-helices is pH-dependent, and upon annealing the sample, a transition from α-helices to β-sheets can be controlled, as indicated by optical and vibrational spectroscopies. Imaging the materials confirms morphologies similar to that of a lipid bilayer but rich in α-helices. Annealing the samples yields a shift in the morphology from bilayers to curled disks, fibers, and sheets. The structural robustness of the material can facilitate the incorporation of many functions into the bilayer assembly.

Oligomer Formation Propensities of Dimeric Bundle Peptides Correlate with Cell Penetration Abilities

LK-3, an amphipathic dimeric peptide linked by two disulfide bonds, and related isomeric bundles were synthesized, and their cell penetrating abilities were investigated. The measurements using size exclusion chromatography and dynamic light scattering show that LK-3 and its isomers form cell penetrating oligomers. Calculations, performed for various types of peptide isomers, elucidate a strong correlation between the amphipathic character of dimers and cell penetration ability. The results suggest that the amphipathicities of LK-3 and related bundle dimers are responsible for their oligomerization propensities which in turn determine their cell penetrating abilities. The observations made in this study provide detailed information about the mechanism of cell uptake of LK-3 and suggest a plausible insight of the early stage of nanoparticle formation of the cell penetrating amphipathic peptides.

Adsorption and conformations of lysozyme and α-lactalbumin at a water-octane interface

As proteins contain both hydrophobic and hydrophilic amino acids, they will readily adsorb onto interfaces between water and hydrophobic fluids such as oil. This adsorption normally causes changes in the protein structure, which can result in loss of protein function and irreversible adsorption, leading to the formation of protein interfacial films. While this can be advantageous in some applications (e.g., food technology), in most cases it limits our ability to exploit protein functionality at interfaces. To understand and control protein interfacial adsorption and function, it is necessary to understand the microscopic conformation of proteins at liquid interfaces. In this paper, molecular dynamics simulations are used to investigate the adsorption and conformation of two similar proteins, lysozyme and α-lactalbumin, at a water-octane interface. While they both adsorb onto the interface, α-lactalbumin does so in a specific orientation, mediated by two amphipathic helices, while lysozyme adsorbs in a non-specific manner. Using replica exchange simulations, both proteins are found to possess a number of distinct interfacial conformations, with compact states similar to the solution conformation being most common for both proteins. Decomposing the different contributions to the protein energy at oil-water interfaces suggests that conformational change for α-lactalbumin, unlike lysozyme, is driven by favourable protein-oil interactions. Revealing these differences between the factors that govern the conformational change at interfaces in otherwise similar proteins can give insight into the control of protein interfacial adsorption, aggregation, and function.

Clustered versus Uniform Display of GALA-Peptides on Carrier Nanoparticles: Enhancing the Permeation of Noncharged Fluid Lipid Membranes

GALA-peptide is a random coil in neutral pH; in acidic pH, it becomes an amphipathic α-helix that aggregates in solution, possibly via its hydrophobic facet that runs along the helix's long axis. In the presence of fluid lipid membranes, the GALA-helix exhibits membrane-active properties that originate from the same hydrophobic facet; these properties make GALA a candidate for inclusion in drug delivery systems requiring permeation of the endosomal membrane to enable drug escape into the cytoplasm. Previous work has shown that the uniform functionalization of carrier nanoparticles with GALA-peptides improved their membrane activity and enhanced the endosomal escape of delivered therapeutics. The present study aims to evaluate the potential role of altering membrane activity via cluster-displayed GALA-peptides (for higher local valency) on the surface of carrier nanoparticles. The presentation of GALA-peptides on carrier nanoparticles was also designed to be pH-dependent. The peptide display on the surface of the carrier nanoparticles was uniform in neutral pH; in the acidic endosomal pH, the surface of nanocarriers formed domains (patches) with high local densities of GALA-peptides. The interactions between GALA-functionalized carrier nanoparticles and target lipid vesicles, utilized as endosome membrane surrogates that were used to primarily capture the fluid nature of these membranes, were studied as a function of pH. At endosomal pH values, ranging from 5.5 to 5.0, the greatest permeability of target membranes was induced by nanocarriers with clustered rather than uniformly displayed GALA. This enhancing effect had an optimum; at even more acidic pH values, too close an approximation of GALA-peptides residing within the same patches resulted in preferential intrapatch peptide interactions rather than interactions with the apposing target lipid membranes. This behavior could have the same physicochemical origin as the aforementioned GALA-peptide aggregation, observed in solution with decreasing pH at increasing peptide concentrations. The findings of this study support the potential of utilizing the clustered display of GALA-peptides on carrier nanoparticles to increase the permeation of fluid membranes used herein to capture this critical physical property of endosomal membranes and therefore to ultimately improve the endosomal escape of delivered therapeutics, enhancing therapeutic efficacy.

Assembly of Huntingtin headpiece into α-helical bundles

Protein aggregation is a hallmark of neurodegenerative disorders. In this group of brain-related disorders, a disease-specific "host" protein or fragment misfolds and adopts a metastatic, aggregate-prone conformation. Often, this misfolded conformation is structurally and thermodynamically different from its native state. Intermolecular contacts, which arise in this non-native state, promote aggregation. In this regard, understanding the molecular principles and mechanisms that lead to the formation of such a non-native state and further promote the formation of the critical nucleus for fiber growth is essential. In this study, the authors analyze the aggregation propensity of Huntingtin headpiece (htt^NT), which is known to facilitate the polyQ aggregation, in relation to the helix mediated aggregation mechanism proposed by the Wetzel group. The authors demonstrate that even though htt^NT displays a degenerate conformational spectrum on its own, interfaces of macroscopic or molecular origin can promote the α-helix conformation, eliminating all other alternatives in the conformational phase space. Our findings indicate that htt^NT molecules do not have a strong orientational preference for parallel or antiparallel orientation of the helices within the aggregate. However, a parallel packed bundle of helices would support the idea of increased polyglutamine concentration, to pave the way for cross-β structures.

Assembly of Triblock Amphiphilic Peptides into One-Dimensional Aggregates and Network Formation

Peptide assembly plays a key role in both neurological diseases and development of novel biomaterials with well-defined nanostructures. Synthetic model peptides provide a unique platform to explore the role of intermolecular interactions in the assembly process. A triblock peptide architecture designed by the Hartgerink group is a versatile system which relies on Coulomb interactions, hydrogen bonding, and hydrophobicity to guide these peptides' assembly at three different length scales: β-sheets, double-wall ribbon-like aggregates, and finally a highly porous network structure which can support gels with ≤1% by weight peptide concentration. In this study, by using molecular dynamics simulations of a structure based implicit solvent coarse grained model, we analyzed this hierarchical assembly process. Parametrization of our CG model is based on multiple-state points from atomistic simulations, which enables this model to represent the conformational adaptability of the triblock peptide molecule based on the surrounding medium. Our results indicate that emergence of the double-wall β-sheet packing mechanism, proposed in light of the experimental evidence, strongly depends on the subtle balance of the intermolecular forces. We demonstrate that, even though backbone hydrogen bonding dominates the early nucleation stages, depending on the strength of the hydrophobic and Coulomb forces, alternative structures such as zero-dimensional aggregates with two β-sheets oriented orthogonally (which we refer to as a cross-packed structure) and β-sheets with misoriented hydrophobic side chains are also feasible. We discuss the implications of these competing structures for the three different length scales of assembly by systematically investigating the influence of density, counterion valency, and hydrophobicity.

Conformations of Myoglobin-Derived Peptides at the Air-Water Interface

The conformational change exhibited by proteins at liquid interfaces, such as the air-water and oil-water interfaces, has long been of interest both for understanding protein structure outside of native environments and for applications in areas including food technology and pharmaceuticals. Using molecular simulation, this article studies the conformations of two peptides derived from myoglobin, for which the emulsification behavior has been studied. Both peptides were found to readily adsorb onto the air-water interface, with one of these (experimentally, the more effective stabilizer) adopting a flat, extended conformation and the other peptide remaining close to its solution conformation.

Conformation and Aggregation of LKα14 Peptide in Bulk Water and at the Air/Water Interface

Historically, the protein folding problem has mainly been associated with understanding the relationship between amino acid sequence and structure. However, it is known that both the conformation of individual molecules and their aggregation strongly depend on the environmental conditions. Here, we study the aggregation behavior of the model peptide LKα14 (with amino acid sequence LKKLLKLLKKLLKL) in bulk water and at the air/water interface. We start by a quantitative analysis of the conformational space of a single LKα14 in bulk water. Next, in order to analyze the aggregation tendency of LKα14, by using the umbrella sampling technique we calculate the potential of mean force for pulling a single peptide from an n-molecule aggregate. In agreement with the experimental results, our calculations yield the optimal aggregate size as four. This equilibrium state is achieved by two opposing forces: Coulomb repulsion between the lysine side chains and the reduction of solvent accessible hydrophobic surface area upon aggregation. At the vacuum/water interface, however, even dimers of LKα14 become marginally stable, and any larger aggregate falls apart instantaneously. Our results indicate that even though the interface is highly influential in stabilizing the α-helix conformation for a single molecule, it significantly reduces the attraction between two LKα14 peptides, along with their aggregation tendency.