Kinetic Analysis of Self-Replicating Peptides: Possibility of Chiral Amplification in Open Systems

A simplified kinetic model scheme is presented that addresses the main reactions of two recently reported peptide self-replicators. Experimentally observed differences in the autocatalytic efficiency between these two systems-- caused by variations in the peptide sequences--and the possible effect of chiral amplification under heterochiral reaction conditions were evaluated. Our numerical simulations indicated that differences in the catalytic performance are exclusively due to pronounced variations in the rate parameters that control the reversible and hydrophobic interactions in the reaction system but neither to alterations in the underlying reaction network nor to changes in the stoichiometry of the involved aggregation processes. Model predictions further demonstrated the possible existence of chiral amplification if peptide self-replication is performed under heterochiral reaction conditions. Pointing into the direction of a possible cause for biomolecular homochirality, it was found that in open flow reactors, keeping the system under non-equilibrium conditions, a remarkable amplification of enantiomeric excess could be achieved. According to our modeling, this is due to a chiroselective autocatalytic effect and a meso-type separation process both of which are assumed to be intrinsic for the underlying dynamics of heterochiral peptide self-replication.     

Stereochemical specificity of Alzheimer's disease beta-peptide assembly

The formation and growth of insoluble amyloid deposits composed primarily of the human beta-amyloid peptide (A beta) in brain is an essentially invariant feature of Alzheimer's disease (AD) and is widely believed to contribute to the progressive neurodegeneration of the disorder. To probe the specificity of amyloid formation and growth, we synthesized and examined the self-assembly of D- and L-stereoisomers of A beta in vitro. While both enantiomers formed insoluble aggregates at similar rates with amyloid-like fibrillar morphology, deposition of soluble A beta peptide onto preexisting A beta aggregates was stereospecific. Although the L-peptide deposited readily onto immobilized L-A beta aggregates with first-order kinetic dependence on soluble peptide concentration, essentially no association between the D-peptide and L-template was observed. Similarly, the D-peptide deposited with first-order kinetics onto a D-A beta aggregate template but did not deposit onto a similar template composed of aggregates of the L-enantiomer. Furthermore, although the L-A beta isomer deposited onto authentic AD amyloid in preparations of unfixed AD brain, no focal association between the D-peptide and brain amyloid was detected. These results establish that deposition of soluble A beta onto preexisting amyloid template is stereospecific, likely involving direct docking interactions between peptide backbone and/or side chains rather than simple hydrophobic association.     

The membrane-binding properties of a class A amphipathic peptide

The membrane-binding properties of a class A amphipathic peptide (18D) were investigated using two different immobilized model membrane systems. The first system involved the use of surface plasmon resonance (SPR) to study the binding of 18D to dimyristylphosphatidylcholine (DMPC) and dimyristylphosphatidylglycerol (DMPG), which allowed peptide binding to be monitored in real time. The SPR experiments indicated stronger binding of 18D to DMPG than DMPC, which kinetic analysis revealed was due to a faster on-rate. The second model membrane system involved immobilized membrane chromatography in which the binding of 18D to either DMPC or DMPG monolayers covalently linked to silica particles was analysed by elution chromatography. Stronger binding affinity of 18D was also obtained with the negatively charged phosphatidylglycerol (PG) monolayer compared to the phosphatidylcholine (PC) monolayer, which was consistent with the SPR results. Non-linear binding behaviour of 18D to the immobilized lipid monolayers was also observed, which suggests that the peptide undergoes conformational and orientational changes upon binding to the immobilized PC and PG ligands. Significant band broadening was also observed on both monolayers, with larger bandwidths obtained on the PC surface, indicating slower binding and orientation kinetics with the zwitterionic surface. The dependence of logk' on the percentage of methanol also demonstrated a bimodal interaction whereby hydrophobic forces predominated at higher temperatures and methanol concentrations, while at lower temperatures, electrostatic and other polar forces also made a contribution to the affinity of the peptides for the lipid monolayer particularly. Overall, these results demonstrate the complementary use of these two lipid biosensors which allows the role of hydrophobic and electrostatic forces in peptide–membrane interactions to be studied and insight gained into the kinetic factors associated with these interactions.     

Photo-crosslinking analysis of preferential interactions between a transmembrane peptide and matching lipids

In this study, a novel method is presented by which the molecular environment of a transmembrane peptide can be investigated directly. This was achieved by incorporating a photoactivatable crosslinking probe in the hydrophobic segment of a model transmembrane peptide. When this peptide was incorporated into lipid bilayers and irradiated with UV light, a covalent bond was formed between the crosslinking probe and a lipid. This crosslinking reaction could be visualized by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and the resulting product could be characterized by mass spectrometry. By use of phospholipases, it was demonstrated that the peptide crosslinks to both acyl chains of the lipids. The peptide showed a clear preference to partition into fluid lipids and was excluded from lipids in the gel phase. However, when the peptide was incorporated into bilayers containing two lipid species with different acyl chain lengths, molecular sorting of the lipids around the peptide based on hydrophobic matching was not observed. It is proposed that the size of the transmembrane part plays an important role in the dynamic interactions of membrane proteins with the surrounding lipids and hence in determining whether molecular sorting can occur.     

Deconstructing the DGAT1 Enzyme: Membrane Interactions at Substrate Binding Sites

Diacylglycerol acyltransferase 1 (DGAT1) is a key enzyme in the triacylglyceride synthesis pathway. Bovine DGAT1 is an endoplasmic reticulum membrane-bound protein associated with the regulation of fat content in milk and meat. The aim of this study was to evaluate the interaction of DGAT1 peptides corresponding to putative substrate binding sites with different types of model membranes. Whilst these peptides are predicted to be located in an extramembranous loop of the membrane-bound protein, their hydrophobic substrates are membrane-bound molecules. In this study, peptides corresponding to the binding sites of the two substrates involved in the reaction were examined in the presence of model membranes in order to probe potential interactions between them that might influence the subsequent binding of the substrates. Whilst the conformation of one of the peptides changed upon binding several types of micelles regardless of their surface charge, suggesting binding to hydrophobic domains, the other peptide bound strongly to negatively-charged model membranes. This binding was accompanied by a change in conformation, and produced leakage of the liposome-entrapped dye calcein. The different hydrophobic and electrostatic interactions observed suggest the peptides may be involved in the interactions of the enzyme with membrane surfaces, facilitating access of the catalytic histidine to the triacylglycerol substrates.     

Cooperativity of hydrophobic anchor interactions: evidence for epitope selection by MHC class II as a folding process

Peptide binding to MHC class II (MHCII) molecules is stabilized by hydrophobic anchoring and hydrogen bond formation. We view peptide binding as a process in which the peptide folds into the binding groove and to some extent the groove folds around the peptide. Our previous observation of cooperativity when analyzing binding properties of peptides modified at side chains with medium to high solvent accessibility is compatible with such a view. However, a large component of peptide binding is mediated by residues with strong hydrophobic interactions that bind to their respective pockets. If these reflect initial nucleation events they may be upstream of the folding process and not show cooperativity. To test whether the folding hypothesis extends to these anchor interactions, we measured dissociation and affinity to HLA-DR1 of an influenza hemagglutinin-derived peptide with multiple substitutions at major anchor residues. Our results show both negative and positive cooperative effects between hydrophobic pocket interactions. Cooperativity was also observed between hydrophobic pockets and positions with intermediate solvent accessibility, indicating that hydrophobic interactions participate in the overall folding process. These findings point out that predicting the binding potential of epitopes cannot assume additive and independent contributions of the interactions between major MHCII pockets and corresponding peptide side chains.     

Design and synthesis of 3alpha-helix peptides forming a cavity for a fluorescent ligand

As a model of receptor protein, a series of 3alpha-helix bundle peptides constructed on a template peptide were designed so as to possess a hydrophobic cavity. The size of cavity was modulated by simple replacements of Leu residues to Ala residues in the hydrophobic core. Binding abilities to 8-anilino-1-naphthalenesulfonic acid (ANS) were estimated by the increase of fluorescence intensity. The peptide having three or four Ala residues in the hydrophobic core remarkably increased the binding ability for ANS, though the peptide having two Ala residues gave an inefficient cavity for ANS. The peptide having six Ala residues decreased the binding ability due to crucial destabilization of the helix bundle structure. This scaffold can be utilized to a receptor model, while further tuning of the sequence is necessary.     

Control of the transmembrane orientation and interhelical interactions within membranes by hydrophobic helix length

We examined the effect of the length of the hydrophobic core of Lys-flanked poly(Leu) peptides on their behavior when inserted into model membranes. Peptide structure and membrane location were assessed by the fluorescence emission lambdamax of a Trp residue in the center of the peptide sequence, the quenching of Trp fluorescence by nitroxide-labeled lipids (parallax analysis), and circular dichroism. Peptides in which the hydrophobic core varied in length from 11 to 23 residues were found to be largely alpha-helical when inserted into the bilayer. In dioleoylphosphatidylcholine (diC18:1PC) bilayers, a peptide with a 19-residue hydrophobic core exhibited highly blue-shifted fluorescence, an indication of Trp location in a nonpolar environment, and quenching localized the Trp to the bilayer center, an indication of transmembrane structure. A peptide with an 11-residue hydrophobic core exhibited emission that was red-shifted, suggesting a more polar Trp environment, and quenching showed the Trp was significantly displaced from the bilayer center, indicating that this peptide formed a nontransmembranous structure. A peptide with a 23-residue hydrophobic core gave somewhat red-shifted fluorescence, but quenching demonstrated the Trp was still close to the bilayer center, consistent with a transmembrane structure. Analogous behavior was observed when the behavior of individual peptides was examined in model membranes with various bilayer widths. Other experiments demonstrated that in diC18:1PC bilayers the dilution of the membrane concentration of the peptide with a 23-residue hydrophobic core resulted in a blue shift of fluorescence, suggesting the red-shifted fluorescence at higher peptide concentrations was due to helix oligomerization. The intermolecular self-quenching of rhodamine observed when the peptide was rhodamine-labeled, and the concentration dependence of self-quenching, supported this conclusion. These studies indicate that the mismatch between helix length and bilayer width can control membrane location, orientation, and helix-helix interactions, and thus may mismatch control both membrane protein folding and the interactions between membrane proteins.     

Trifluoroethanol may form a solvent matrix for assisted hydrophobic interactions between peptide side chains

Several models for interactions between trifluoroethanol (TFE) and peptides and proteins have recently been proposed, but none have been able to rationalize the puzzling observations that on the one hand TFE can stabilize some hydrophobic interactions in secondary structures, but on the other can also melt the hydrophobic cores of globular proteins. The former is illustrated in this paper by the effect of TFE on a short elastin peptide, GVG(VPGVG)(3), which forms type II beta-turns stabilized by hydrophobic interactions between two intra-turn valine side chains. This folding, driven by increasing the entropy of bulk water, is stimulated in TFE-water mixtures and/or by raising the temperature. To explain these apparently contradictory observations, we propose a model in which TFE clusters locally assist the folding of secondary structures by first breaking down interfacial water molecules on the peptide and then providing a solvent matrix for further side chain--side chain interactions. This model also provides an explanation for TFE-induced transitions between secondary structures, in which the TFE clusters may redirect non-local to local interactions.     

Peptide models of folding initiation sites of bovine beta-lactoglobulin: identification of nativelike hydrophobic interactions involving G and H strands

In an attempt to characterize the early folding events in bovine beta-lactoglobulin (BLG), a set of peptides, covering the flexible N-terminal region and the stable C-terminus beta-core, was synthesized and analyzed by circular dichroism and by nuclear magnetic resonance in water, trifluoroethanol (TFE), and sodium dodecyl sulfate (SDS) below and above the critical micellar concentration. The role of local and long-range hydrophobic interactions in guiding the folding has been investigated. For the peptide fragment covering the more flexible N-terminal region of BLG (beta-strands A, B), where both theoretical predictions and kinetic refolding experiments suggested the formation of non-native alpha-helix, no native long-range contacts were identified, and a helical secondary structure was stabilized only in the presence of 25 mM SDS. At variance, in 50% (v/v) TFE, native, long-range hydrophobic interactions were observed in the peptide covering the core region comprising G and H beta-strands. The side chains involved in these interactions form a nativelike hydrophobic cluster, thus suggesting that the GH region may act as the folding initiation site for BLG. This result is reinforced by the identification, in the urea denaturated BLG, of residual structure located at the level of the GH interface, as evidenced by NMR analysis. These results, in excellent agreement with kinetic, thermodynamic, and cold denaturation folding data, once more underline the utmost importance of the GH region for the stability and folding of BLG. Severe aggregation effects prevented the structural analysis of the peptide covering the EFGH region, indicating that this larger segment does not represent an independent folding domain and that the terminal alpha-helix is necessary for stabilizing the BLG folding core.     

A theoretical model for the association of amphiphilic transmembrane peptides in lipid bilayers

A theoretical model is proposed for the association of trans-bilayer peptides in lipid bilayers. The model is based on a lattice model for the pure lipid bilayer, which accounts accurately for the most important conformational states of the lipids and their mutual interactions and statistics. Within the lattice formulation the bilayer is formed by two independent monolayers, each represented by a triangular lattice, on which sites the lipid chains are arrayed. The peptides are represented by regular objects, with no internal flexibility, and with a projected area on the bilayer plane corresponding to a hexagon with seven lattice sites. In addition, it is assumed that each peptide surface at the interface with the lipid chains is partially hydrophilic, and therefore interacts with the surrounding lipid matrix via selective anisotropic forces. The peptides would therefore assemble in order to shield their hydrophilic residues from the hydrophobic surroundings. The model describes the self-association of peptides in lipid bilayers via lateral and rotational diffusion, anisotropic lipid-peptide interactions, and peptide-peptide interactions involving the peptide hydrophilic regions. The intent of this model study is to analyse the conditions under which the association of trans-bilayer and partially hydrophilic peptides (or their dispersion in the lipid matrix) is lipid-mediated, and to what extent it is induced by direct interactions between the hydrophilic regions of the peptides. The model properties are calculated by a Monte Carlo computer simulation technique within the canonical ensemble. The results from the model study indicate that direct interactions between the hydrophilic regions of the peptides are necessary to induce peptide association in the lipid bilayer in the fluid phase. Furthermore, peptides within each aggregate are oriented in such a way as to shield their hydrophilic regions from the hydrophobic environment. The average number of peptides present in the aggregates formed depends on the degree of mismatch between the peptide hydrophobic length and the lipid bilayer hydrophobic thickness: The lower the degree of mismatch is the higher this number is. Received: 30 December 1996 / Accepted: 9 May 1997     

Interactions of cationic-hydrophobic peptides with lipid bilayers: a Monte Carlo simulation method

We present a computational model of the interaction between hydrophobic cations, such as the antimicrobial peptide, Magainin2, and membranes that include anionic lipids. The peptide's amino acids were represented as two interaction sites: one corresponds to the backbone alpha-carbon and the other to the side chain. The membrane was represented as a hydrophobic profile, and its anionic nature was represented by a surface of smeared charges. Thus, the Coulombic interactions between the peptide and the membrane were calculated using the Gouy-Chapman theory that describes the electrostatic potential in the aqueous phase near the membrane. Peptide conformations and locations near the membrane, and changes in the membrane width, were sampled at random, using the Metropolis criterion, taking into account the underlying energetics. Simulations of the interactions of heptalysine and the hydrophobic-cationic peptide, Magainin2, with acidic membranes were used to calibrate the model. The calibrated model reproduced structural data and the membrane-association free energies that were measured also for other basic and hydrophobic-cationic peptides. Interestingly, amphipathic peptides, such as Magainin2, were found to adopt two main membrane-associated states. In the first, the peptide resided mostly outside the polar headgroups region. In the second, which was energetically more favorable, the peptide assumed an amphipathic-helix conformation, where its hydrophobic face was immersed in the hydrocarbon region of the membrane and the charged residues were in contact with the surface of smeared charges. This dual behavior provides a molecular interpretation of the available experimental data.     

Interactions of hydrophobic peptides with lipid bilayers: Monte Carlo simulations with M2delta

We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.     

Conformational studies of a peptide corresponding to a region of the C-terminus of ribonuclease A: implications as a potential chain-folding initiation site

Conformational properties of the OT-16 peptide, the C-terminal 20 amino acids of RNase A, were examined by nonradiative energy transfer. A modified OT-16 peptide was prepared by solid-phase synthesis with the inclusion of diaminobutyric acid (DABA) at the C-terminus. The OT-16-DABA peptide was labeled with a fluorescent 1,5-dimethylaminonaphthalene sulfonyl (dansyl, DNS) acceptor at the N-terminal amine and a fluorescent naphthoxyacetic acid (NAA) donor at the gamma-amine of the DABA located at the C-terminus of the peptide by using an orthogonal protection scheme. Energy transfer was monitored in DNS-OT-16-DABA-NAA by using both fluorescence intensity (sensitized emission) and lifetime (donor quenching) experiments. The lifetime data indicate that the peptide system is a dynamic, flexible one. A detailed analysis, based on a dynamic model that includes a skewed Gaussian function to model the equilibrium distribution of interprobe distances and a mutual diffusion coefficient between the two probes to model conformational dynamics in the peptide [Beechem & Haas (1989) Biophys. J. 55, 1225.], identified the existence of a partially ordered structure (relatively narrow distribution of interprobe distances) at temperatures greater than or equal to 20 degrees C in the absence of denaturant. The width and the position of the average of the distributions decrease with increasing temperature, in this range; this suggests that the structure is stabilized by hydrophobic interactions. In addition, the peptide undergoes cold denaturation at around 1.5 degrees C as indicated by broadening of the distance distribution. The addition of 6 M guanidine hydrochloride (Gdn-HCl) also broadens the distance distribution significantly, presumably by eliminating the hydrophobic interactions and unfolding the peptide. The results of the analysis of the distance distribution demonstrate that (1) nonradiative energy transfer can be used to study the conformational dynamics of peptides on the nanosecond time scale, (2) a partially ordered structure of OT-16-DABA exists in solution under typical refolding conditions, and (3) structural constraints (presumably hydrophobic interactions) necessary for the formation of a chain-folding initiation site in RNase A are also present in the OT-16-DABA peptide in the absence of denaturant and are disrupted by Gdn-HCl.     

Involvement of electrostatic interactions in the mechanism of peptide folding induced by sodium dodecyl sulfate binding

Sodium dodecyl sulfate (SDS) has consistently been shown to induce secondary structure, particularly alpha-helices, in polypeptides, and is commonly used to model membrane and other hydrophobic environments. However, the precise mechanism by which SDS induces these conformational changes remains unclear. To examine the role of electrostatic interactions in this mechanism, we have designed two hydrophilic, charged amphipathic alpha-helical peptides, one basic (QAPAYKKAAKKLAES) and the other acidic (QAPAYEEAAEELAKS), and their structures were studied by CD and NMR. The design of the peptides is based on the sequence of the segment of residues 56-70 of human platelet factor 4 [PF4(56-70), QAPLYKKIIKKLLES]. Both peptides were unstructured in water, and in the presence of neutral, zwitterionic, or cationic detergents. However, in SDS at neutral pH, the basic peptide folded into an alpha-helix. By contrast, the pH needed to be lowered to 1.8 before alpha-helix formation was observed for the acidic peptide. Strong, attractive electrostatic interactions, between the anionic groups of SDS and the cationic groups of the lysines, appeared to be necessary to initiate the folding of the basic peptide. NMR analysis showed that the basic peptide was fully embedded in SDS-peptide micelles, and that its three-dimensional alpha-helical structure could be superimposed on that of the native structure of PF4(56-70). These results enabled us to propose a working model of the basic peptide-SDS complex, and a mechanism for SDS-induced alpha-helical folding. This study demonstrates that, while the folding of peptides is mostly driven by hydrophobic effects, electrostatic interactions play a significant role in the formation and the stabilization of SDS-induced structure.     

Simulation of non-enzymatic template-directed synthesis of oligonucleotides and peptides

A simulation model is proposed for the template- and sequence-directed (TSD) condensation of two trideoxyribonucleotide 3'-phosphate molecules into a hexameric template with palindromic sequence studied experimentally by von Kiedrowski (1986;Angew. Int. Ed. Engl.25, 932--935). The model simulates reasonably well the kinetics of synthesis of both the template, and the pyrophosphate product which is not directly involved in the autocatalytic reaction. It offers quantitative approximation of the different rate constants of the processes involved in the reaction. The model simulates and gives predictions for the influence of factors such as the initial concentrations of the trimers and the template, and gives predictions for the effect of temperature on the dynamics of the autocatalytic reaction. The model also simulates well the production rate of a different self-replicating system (coiled coil peptide) used in the experiments of Lee et al. (1997;Nature390, 591--594). Comparing the different rate constants, it seems that chain elongation occurs at higher rates in the peptide system (at 23 degrees C) than in the nucleotide one (at 0 degrees C), but that the relative contribution of template-directed synthesis is significantly larger with the nucleotides.     

Fundamental processes of protein folding: measuring the energetic balance between helix formation and hydrophobic interactions

Theories of protein folding often consider contributions from three fundamental elements: loops, hydrophobic interactions, and secondary structures. The pathway of protein folding, the rate of folding, and the final folded structure should be predictable if the energetic contributions to folding of these fundamental factors were properly understood. alphatalpha is a helix-turn-helix peptide that was developed by de novo design to provide a model system for the study of these important elements of protein folding. Hydrogen exchange experiments were performed on selectively 15N-labeled alphatalpha and used to calculate the stability of hydrogen bonds within the peptide. The resulting pattern of hydrogen bond stability was analyzed using a version of Lifson-Roig model that was extended to include a statistical parameter for tertiary interactions. This parameter, x, represents the additional statistical weight conferred upon a helical state by a tertiary contact. The hydrogen exchange data is most closely fit by the XHC model with an x parameter of 9.25. Thus the statistical weight of a hydrophobic tertiary contact is approximately 5.8x the statistical weight for helix formation by alanine. The value for the x parameter derived from this study should provide a basis for the understanding of the relationship between hydrophobic cluster formation and secondary structure formation during the early stages of protein folding.     

Interactions of calmodulin with death-associated protein kinase peptides: experimental and modeling studies

We have studied the interactions between calmodulin (CaM) and three target peptides from the death-associated protein kinase (DAPK) protein family using both experimental and modeling methods, aimed at determining the details of the underlying biological regulation mechanisms. Experimentally, calorimetric binding free energies were determined for the complexes of CaM with peptides representing the DAPK2 wild-type and S308D mutant, as well as DAPK1. The observed affinity of CaM was very similar for all three studied peptides. The DAPK2 and DAPK1 peptides differ significantly in sequence and total charge, while the DAPK2 S308D mutant is designed to model the effects of DAPK2 Ser308 phosphorylation. The crystal structure of the CaM-DAPK2 S308D mutant peptide is also reported. The structures of CaM-DAPK peptide complexes present a mode of CaM-kinase interaction, in which bulky hydrophobic residues at positions 10 and 14 are both bound to the same hydrophobic cleft. To explain the microscopic effects underlying these interactions, we performed free energy calculations based on the approximate MM-PBSA approach. For these highly charged systems, standard MM-PBSA calculations did not yield satisfactory results. We proposed a rational modification of the approach which led to reasonable predictions of binding free energies. All three complexes are strongly stabilized by two effects: electrostatic interactions and buried surface area. The strong favorable interactions are to a large part compensated by unfavorable entropic terms, in which vibrational entropy is the largest contributor. The electrostatic component of the binding free energy followed the trend of the overall peptide charge, with strongest interactions for DAPK1 and weakest for the DAPK2 mutant. The electrostatics was dominated by interactions of the positively charged residues of the peptide with the negatively charged residues of CaM. The nonpolar binding free energy was comparable for all three peptides, the largest contribution coming from the Trp305. About two-thirds of the buried surface area corresponds to nonpolar residues, showing that hydrophobic interactions play an important role in these CaM-peptide complexes. The simulation results agree with the experimental data in predicting a small effect of the S308D mutation on CaM interactions with DAPK2, suggesting that this mutation is not a good model for the S308 phosphorylation.     

Peptide-Templated Nucleic Acid Ligation

Abstract Short oligonucleotide and peptide replicators have been described. To determine whether cross-replication could have occurred between such systems, we have attempted to show that peptides can specifically template the ligation of nucleic acids. A complex between a 35-mer anti-Rev RNA aptamer and a 17-mer arginine-rich motif (ARM) peptide from the HIV-1 Rev protein served as a model system. Aptamer half-molecules were activated for ligation via two activation chemistries, representing two distinct kinetic possibilities for early replicators. Cyanogen bromide activation was transient relative to oligonucleotides that terminated with a 5′-iodine and a 3′phosphorothioate, respectively. The Rev ARM specifically enhanced the degree or rate of ligation by both methods: there was a 10-fold increase in the production of full-length aptamer in the presence of cyanogen bromide and a 5.9- to 7.6-fold enhancement in the rate of ligation for stably activated aptamer half-molecules. These results support the possibility that life could have originated with peptide replicators and transitioned to nucleic acid replicators or that peptide and nucleic acid replicators could have been interdependent.     

Mechanism of the cell-penetrating peptide transportan 10 permeation of lipid bilayers

The mechanism of the interaction between the cell-penetrating peptide transportan 10 (tp10) and phospholipid membranes was investigated. Tp10 induces graded release of the contents of phospholipid vesicles. The kinetics of peptide association with vesicles and peptide-induced dye efflux from the vesicle lumen were examined experimentally by stopped-flow fluorescence. The experimental kinetics were analyzed by directly fitting to the data the numerical solution of mathematical kinetic models. A very good global fit was obtained using a model in which tp10 binds to the membrane surface and perturbs it because of the mass imbalance thus created across the bilayer. The perturbed bilayer state allows peptide monomers to insert transiently into its hydrophobic core and cross the membrane, until the peptide mass imbalance is dissipated. In that transient state tp10 "catalyzes" dye efflux from the vesicle lumen. These conclusions are consistent with recent reports that used molecular dynamics simulations to study the interactions between peptide antimicrobials and phospholipid bilayers. A thermodynamic analysis of tp10 binding and insertion in the bilayer using water-membrane transfer hydrophobicity scales is entirely consistent with the model proposed. A small bilayer perturbation is both necessary and sufficient to achieve very good agreement with the model, indicating that the role of the lipids must be included to understand the mechanism of cell-penetrating and antimicrobial peptides.