Mehmet Sayar

Self-assembling multidomain peptide fibers with aromatic cores.

Self-assembling multidomain peptides have been shown to have desirable properties, such as the ability to form hydrogels that rapidly recover following shear-thinning and the potential to be tailored by amino acid selection to vary their elasticity and encapsulate and deliver proteins and cells. Here we describe the effects of substitution of aliphatic hydrophobic amino acids in the central domain of the peptide for the aromatic amino acids phenylalanine, tyrosine, and tryptophan. While the basic nanofibrous morphology is retained in all cases, selection of the particular core residues results in switching from antiparallel hydrogen bonding to parallel hydrogen bonding in addition to changes in nanofiber morphology and in hydrogel rheological properties. Peptide nanofiber assemblies are investigated by circular dichroism polarimetry, infrared spectroscopy, atomic force microscopy, transmission and scanning electron microscopy, oscillatory rheology, and molecular dynamics simulations. Results from this study will aid in designing next generation cell scaffolding materials.

Finite-size polyelectrolyte bundles at thermodynamic equilibrium

We present the results of extensive computer simulations performed on solutions of monodisperse charged rod-like polyelectrolytes in the presence of trivalent counterions. To overcome energy barriers we used a combination of parallel tempering and hybrid Monte Carlo techniques. Our results show that for small values of the electrostatic interaction the solution mostly consists of dispersed single rods. The potential of mean force between the polyelectrolyte monomers yields an attractive interaction at short distances. For a range of larger values of the Bjerrum length, we find finite-size polyelectrolyte bundles at thermodynamic equilibrium. Further increase of the Bjerrum length eventually leads to phase separation and precipitation. We discuss the origin of the observed thermodynamic stability of the finite-size aggregates.

Twist-writhe partitioning in a coarse-grained DNA minicircle model.

Here we present a systematic study of supercoil formation in DNA minicircles under varying linking number by using molecular-dynamics simulations of a two-bead coarse-grained model. Our model is designed with the purpose of simulating long chains without sacrificing the characteristic structural properties of the DNA molecule, such as its helicity, backbone directionality, and the presence of major and minor grooves. The model parameters are extracted directly from full-atomistic simulations of DNA oligomers via Boltzmann inversion; therefore, our results can be interpreted as an extrapolation of those simulations to presently inaccessible chain lengths and simulation times. Using this model, we measure the twist/writhe partitioning in DNA minicircles, in particular its dependence on the chain length and excess linking number. We observe an asymmetric supercoiling transition consistent with experiments. Our results suggest that the fraction of the linking number absorbed as twist and writhe is nontrivially dependent on chain length and excess linking number. Beyond the supercoiling transition, chains of the order of one persistence length carry equal amounts of twist and writhe. For longer chains, an increasing fraction of the linking number is absorbed by the writhe.

Buckling of stiff polymers: influence of thermal fluctuations.

Marc Emanuel, Hervé Mohrbach, Mehmet Sayar, Helmut Schiessel, and Igor M. Kulić Instituut-Lorentz, Universiteit Leiden, Postbus 9506, 2300 RA Leiden, The Netherlands Institut de Physique, Université Paul Verlaine-Metz, LPMC, CPMB1-FR CNRS 2843, 1 boulevard Arago, 57078 Metz, France Koc University, College of Engineering, 34450 Sariyer, Istanbul, Turkey School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA

A transferable coarse-grained model for diphenylalanine: How to represent an environment driven conformational transition

One of the major challenges in the development of coarse grained (CG) simulation models that aim at biomolecular structure formation processes is the correct representation of an environment-driven conformational change, for example, a folding/unfolding event upon interaction with an interface or upon aggregation. In the present study, we investigate this transferability challenge for a CG model using the example of diphenylalanine. This dipeptide displays a transition from a trans-like to a cis-like conformation upon aggregation as well as upon transfer from bulk water to the cyclohexane/water interface. Here, we show that one can construct a single CG model that can reproduce both the bulk and interface conformational behavior and the segregation between hydrophobic/hydrophilic medium. While the general strategy to obtain nonbonded interactions in the present CG model is to reproduce solvation free energies of small molecules representing the CG beads in the respective solvents, the success of the model strongly depends on nontrivial decisions one has to make to capture the delicate balance between the bonded and nonbonded interactions. In particular, we found that the peptides conformational behavior is qualitatively affected by the cyclohexane/water interaction potential, an interaction that does not directly involve the peptide at all but merely influences the properties of the hydrophobic/hydrophilic interface. Furthermore, we show that a small modification to improve the structural/conformational properties of the CG model could dramatically alter the thermodynamic properties.

Driving Forces for Adsorption of Amphiphilic Peptides to the Air−Water Interface

We have studied the partitioning of amphiphilic peptides at the air-water interface. The free energy of adsorption from bulk to interface was calculated by determining the potential of mean force via atomistic molecular dynamics simulations. To this end a method is introduced to restrain or constrain the center of mass of a group of molecules in a periodic system. The model amphiphilic peptides are composed of alternating valine and asparagine residues. The decomposition of the free energy difference between the bulk and interface is studied for different peptide block lengths. Our analysis revealed that for short amphiphilic peptides the surface driving force dominantly stems from the dehydration of hydrophobic side chains. The only opposing force is associated with the loss of orientational freedom of the peptide at the interface. For the peptides studied, the free energy difference scales linearly with the size of the molecule, since the peptides mainly adopt extended conformations both in bulk and at the interface. The free energy difference depends strongly on the water model, which can be rationalized through the hydration thermodynamics of hydrophobic solutes. Finally, we measured the reduction of the surface tension associated with complete coverage of the interface with peptides.

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

Secondary amphiphilicity is inherent to the secondary structural elements of proteins. By forming energetically favorable contacts with each other these amphiphilic building blocks give rise to the formation of a tertiary structure. Small proteins and peptides, on the other hand, are usually too short to form multiple structural elements and cannot stabilize them internally. Therefore, these molecules are often found to be structurally ambiguous up to the point of a large degree of intrinsic disorder in solution. Consequently, their conformational preference is particularly susceptible to environmental conditions such as pH, salts, or presence of interfaces. In this study we use molecular dynamics simulations to analyze the conformational behavior of two synthetic peptides, LKKLLKLLKKLLKL (LK) and EAALAEALAEALAE (EALA), with built-in secondary amphiphilicity upon forming an alpha-helix. We use these model peptides to systematically study their aggregation and the influence of macroscopic and molecular interfaces on their conformational preferences. We show that the peptides are neither random coils in bulk water nor fully formed alpha helices, but adopt multiple conformations and secondary structure elements with short lifetimes. These provide a basis for conformation-selection and population-shift upon environmental changes. Differences in these peptides’ response to macroscopic and molecular interfaces (presented by an aggregation partner) can be linked to their inherent alpha-helical tendencies in bulk water. We find that the peptides’ aggregation behavior is also strongly affected by presence or absence of an interface, and rather subtly depends on their surface charge and hydrophobicity.

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.

Adsorption, folding, and packing of an amphiphilic peptide at the air/water interface.

Peptide oligomers play an essential role as model compounds for identifying key motifs in protein structure formation and protein aggregation. Here, we present our results, based on extensive molecular dynamics simulations, on adsorption, folding, and packing within a surface monolayer of an amphiphilic peptide at the air/water interface. Experimental results suggest that these molecules spontaneously form ordered monolayers at the interface, adopting a β-hairpin-like structure within the surface layer. Our results reveal that the β-hairpin structure can be observed both in bulk and at the air/water interface. However, the presence of an interface leads to ideal partitioning of the hydrophobic and hydrophilic residues, and therefore reduces the conformational space for the molecule and increases the stability of the hairpin structure. We obtained the adsorption free energy of a single β-hairpin at the air/water interface, and analyzed the enthalpic and entropic contributions. The adsorption process is favored by two main factors: (1) Free-energy reduction due to desolvation of the hydrophobic side chains of the peptide and release of the water molecules which form a cage around these hydrophobic groups in bulk water. (2) Reduction of the total air/water contact area at the interface upon adsorption of the peptide amphiphile. By performing mutations on the original molecule, we demonstrated the relative role of key design features of the peptide. Finally, by analyzing the potential of mean force among two peptides at the interface, we investigated possible packing mechanisms for these molecules within the surface monolayer.

The introduction of hydrogen bond and hydrophobicity effects into the rotational isomeric states model for conformational analysis of unfolded peptides

Relative contributions of local and non-local interactions to the unfolded conformations of peptides are examined by using the rotational isomeric states model which is a Markov model based on pairwise interactions of torsion angles. The isomeric states of a residue are well described by the Ramachandran map of backbone torsion angles. The statistical weight matrices for the states are determined by molecular dynamics simulations applied to monopeptides and dipeptides. Conformational properties of tripeptides formed from combinations of alanine, valine, tyrosine and tryptophan are investigated based on the Markov model. Comparison with molecular dynamics simulation results on these tripeptides identifies the sequence-distant long-range interactions that are missing in the Markov model. These are essentially the hydrogen bond and hydrophobic interactions that are obtained between the first and the third residue of a tripeptide. A systematic correction is proposed for incorporating these long-range interactions into the rotational isomeric states model. Preliminary results suggest that the Markov assumption can be improved significantly by renormalizing the statistical weight matrices to include the effects of the long-range correlations.