beta-Hairpins, alpha-helices, and the intermediates among the secondary structures in the energy landscape of a peptide from a distal beta-hairpin of SH3 domain

Multiple Antimicrobial Effects of Hybrid Peptides Synthesized Based on the Sequence of Ribosomal S1 Protein from Staphylococcus aureus

The need to develop new antimicrobial peptides is due to the high resistance of pathogenic bacteria to traditional antibiotics now and in the future. The creation of synthetic peptide constructs is a common and successful approach to the development of new antimicrobial peptides. In this work, we use a simple, flexible, and scalable technique to create hybrid antimicrobial peptides containing amyloidogenic regions of the ribosomal S1 protein from Staphylococcus aureus. While the cell-penetrating peptide allows the peptide to enter the bacterial cell, the amyloidogenic site provides an antimicrobial effect by coaggregating with functional bacterial proteins. We have demonstrated the antimicrobial effects of the R23F, R23DI, and R23EI hybrid peptides against Staphylococcus aureus, methicillin-resistant S. aureus (MRSA), Pseudomonas aeruginosa, Escherichia coli, and Bacillus cereus. R23F, R23DI, and R23EI can be used as antimicrobial peptides against Gram-positive and Gram-negative bacteria resistant to traditional antibiotics.

Enhanced and effective conformational sampling of protein molecular systems for their free energy landscapes

Protein folding and protein-ligand docking have long persisted as important subjects in biophysics. Using multicanonical molecular dynamics (McMD) simulations with realistic expressions, i.e., all-atom protein models and an explicit solvent, free-energy landscapes have been computed for several systems, such as the folding of peptides/proteins composed of a few amino acids up to nearly 60 amino-acid residues, protein-ligand interactions, and coupled folding and binding of intrinsically disordered proteins. Recent progress in conformational sampling and its applications to biophysical systems are reviewed in this report, including descriptions of several outstanding studies. In addition, an algorithm and detailed procedures used for multicanonical sampling are presented along with the methodology of adaptive umbrella sampling. Both methods control the simulation so that low-probability regions along a reaction coordinate are sampled frequently. The reaction coordinate is the potential energy for multicanonical sampling and is a structural identifier for adaptive umbrella sampling. One might imagine that this probability control invariably enhances conformational transitions among distinct stable states, but this study examines the enhanced conformational sampling of a simple system and shows that reasonably well-controlled sampling slows the transitions. This slowing is induced by a rapid change of entropy along the reaction coordinate. We then provide a recipe to speed up the sampling by loosening the rapid change of entropy. Finally, we report all-atom McMD simulation results of various biophysical systems in an explicit solvent.

Ab initio simulation of a 57-residue protein in explicit solvent reproduces the native conformation in the lowest free-energy cluster

An enhanced conformational sampling method, multicanonical molecular dynamics (McMD), was applied to the ab intio folding of the 57-residue first repeat of human glutamyl- prolyl-tRNA synthetase (EPRS-R1) in explicit solvent. The simulation started from a fully extended structure of EPRS-R1 and did not utilize prior structural knowledge. A canonical ensemble, which is a conformational ensemble thermodynamically probable at an arbitrary temperature, was constructed by reweighting the sampled structures. Conformational clusters were obtained from the canonical ensemble at 300 K, and the largest cluster (i.e., the lowest free-energy cluster), which contained 34% of the structures in the ensemble, was characterized by the highest similarity to the NMR structure relative to all alternative clusters. This lowest free-energy cluster included native-like structures composed of two anti-parallel α-helices. The canonical ensemble at 300 K also showed that a short Gly-containing segment, which adopts an α-helix in the native structure, has a tendency to be structurally disordered. Atomic-level analyses demonstrated clearly that inter-residue hydrophobic interactions drive the helix formation of the Gly-containing segment, and that increasing the hydrophobic contacts accompanies exclusion of water molecules from the vicinity of this segment. This study has shown, for the first time, that the free-energy landscape of a structurally well-ordered protein of about 60 residues is obtainable with an all atom model in explicit water without prior structural knowledge.

Conformational requirement on peptides to exert laminin's activities and search for protein segments with laminin's activities

The human laminin alpha3 chain LG4 module has biological activities of cell adhesion, heparin binding, migration, and neurite outgrowth. The authors had previously identified that the active site of this protein is in residues 1411-1429 (amino-acid sequence = KNSFMALYLSKGRLVFALG called A3G756) and that a three-amino-acid sequence KGR in A3G756 is crucial for exerting the activities. An experiment has shown that a cyclo-hEF3A peptide (a cyclic analog of A3G756) exhibits stronger activities than a linear-hEF3A peptide (a linearized peptide of the cyclo-hEF3A peptide). This experiment implies that adopting a loop conformation may be important for exerting the activities. In this study, the authors first computed the solution structures of the cyclo-hEF3A and linear-hEF3A peptides by molecular dynamics simulations. The obtained conformational ensembles consisted of a variety of conformations, which is a usual property of short peptides in solution. The ensembles involved a fraction where the peptide adopted beta-hairpins and KGR was located at the hairpin head. If there are protein segments that adopt beta-hairpins similar to those sampled from the simulation and have the KGR sequence at the hairpin head, these segments may have some activities. Then, the authors searched a database for segments satisfying these requirements and detected six functional segments. Three of them had laminin's activity, and the remaining three had activities similar to laminin's activities. Analyses on the conformational ensembles of cyclo- and linear-hEF3A peptides suggest that not only the KGR position in the hairpin but also the inter-strand packing is important for exerting laminin's activities.

Simulation study on the disordered state of an Alzheimer's beta amyloid peptide Abeta(12 36) in water consisting of random-structural, beta-structural, and helical clusters

The monomeric Alzheimer's beta amyloid peptide, Abeta, is known to adopt a disordered state in water at room temperature, and a circular dichroism (CD) spectroscopy experiment has provided the secondary-structure contents for the disordered state: 70% random, 25% beta-structural, and 5% helical. We performed an enhanced conformational sampling (multicanonical molecular dynamics simulation) of a 25-residue segment (residues 12-36) of Abeta in explicit water and obtained the conformational ensemble over a wide temperature range. The secondary-structure contents calculated from the conformational ensemble at 300 degrees K reproduced the experimental secondary-structure contents. The constructed free-energy landscape at 300 degrees K was not plain but rugged with five clearly distinguishable clusters, and each cluster had its own characteristic tertiary structure: a helix-structural cluster, two beta-structural clusters, and two random-structural clusters. This indicates that the contribution from the five individual clusters determines the secondary-structure contents experimentally measured. The helical cluster had a similarity with a stable helical structure for monomeric Abeta in 2,2,2-trifluoroethanol (TFE)/water determined by an NMR experiment: The positions of helices in the helical cluster were the same as those in the NMR structure, and the residue-residue contact patterns were also similar with those of the NMR structure. The cluster-cluster separation in the conformational space indicates that free-energy barriers separate the clusters at 300 degrees K. The two beta-structural clusters were characterized by different strand-strand hydrogen-bond (H-bond) patterns, suggesting that the free-energy barrier between the two clusters is due to the H-bond rearrangements.

Folding of the 25 Residue Aβ(12−36) Peptide in TFE/Water: Temperature-Dependent Transition from a Funneled Free-Energy Landscape to a Rugged One

The free-energy landscape of the Alzheimer beta-amyloid peptide Abeta(12-36) in a 40% (v/v) 2,2,2-trifluoroethanol (TFE)/water solution was determined by using multicanonical molecular dynamics simulations. Simulations using this enhanced conformational sampling technique were initiated from a random unfolded polypeptide conformation. Our simulations reliably folded the peptide to the experimental NMR structure, which consists of two linked helices. The shape of the free energy landscape for folding was found to be strongly dependent on temperature: Above 325 K, the overall shape was funnel-like, with the bottom of the funnel coinciding exactly with the NMR structure. Below 325 K, on the other hand, the landscape became increasingly rugged, with the emergence of new conformational clusters connected by low free-energy pathways. Finally, our simulations reveal that water and TFE solvate the polypeptide in different ways: The hydrogen bond formation between TFE and Abeta was enhanced with decreasing temperature, while that between water and Abeta was depressed.

Secondary structure provides a template for the folding of nearby polypeptides

Although protein structures are primarily encoded by their sequences, they are also critically dependent on environmental factors such as solvents and interactions with other molecules. Here we investigate how the folding-energy landscape of a short peptide is altered by interactions with another peptide, by performing atomistic replica-exchange molecular dynamics simulations of polyalanines in various environments. We analyzed the free-energy landscapes of Ala7 and Ala8 in isolation, near an alpha-helix template, and near a beta-strand template. The isolated Ala7 and Ala8 at 270 K were mainly in polyproline II helix conformations and in equilibrium between the alpha-helix and polyproline II helix, respectively, in harmony with the experiment. Interestingly, we found remarkably strong secondary-structure "templating"; namely, the alpha-helix template enhanced alpha-helix conformation and the beta-strand template induced beta-strand conformation in the simulated Ala8. The alpha-helix template lowered the nearby dielectric constant, which strengthened hydrogen bonds in the simulated Ala8, leading to alpha-helix stabilization. The beta-strand template provided hydrogen bond positions to the simulated Ala8, sharply inducing beta-strand structure. With or without templates, the energy landscape of Ala8 is always funnel-like and centered at the alpha-helix conformation, whereas entropic contribution disfavors the alpha-helix, leading to subtle competition. Secondary-structure templating may play a critical role in protein conformation dynamics in the cellular environment.

Calculation of the entropy and free energy of peptides by molecular dynamics simulations using the hypothetical scanning molecular dynamics method

Hypothetical scanning (HS) is a method for calculating the absolute entropy S and free energy F from a sample generated by any simulation technique. With this approach each sample configuration is reconstructed with the help of transition probabilities (TPs) and their product leads to the configuration's probability, hence to the entropy. Recently a new way for calculating the TPs by Monte Carlo (MC) simulations has been suggested, where all system interactions are taken into account. Therefore, this method--called HSMC--is in principle exact where the only approximation is due to insufficient sampling. HSMC has been applied very successfully to liquid argon, TIP3P water, self-avoiding walks on a lattice, and peptides. Because molecular dynamics (MD) is considered to be significantly more efficient than MC for a compact polymer chain, in this paper HSMC is extended to MD simulations as applied to peptides. Like before, we study decaglycine in vacuum but for the first time also a peptide with side chains, (Val)(2)(Gly)(6)(Val)(2). The transition from MC to MD requires implementing essential changes in the reconstruction process of HSMD. Results are calculated for three microstates, helix, extended, and hairpin. HSMD leads to very stable differences in entropy TDeltaS between these microstates with small errors of 0.1-0.2 kcal/mol (T=100 K) for a wide range of calculation parameters with extremely high efficiency. Various aspects of HSMD and plans for future work are discussed.

Evaluation of efficiency of reference potential spatial warping algorithm in conformational sampling of peptides

Molecular dynamics (MD) simulation using the reference potential spatial warping algorithm (REPSWA) is a promising method for computing the conformational equilibrium of a system with a rugged energy surface. Its effectiveness has previously been demonstrated using only a simplified model system of a hydrocarbon chain omitting nonbonded interaction terms from the potential energy function. To evaluate the applicability of REPSWA MD simulation to more realistic problems, we applied it to small peptide systems in an aqueous environment. Difficulties were encountered, however, forcing us to devise several modifications. We investigated their effectiveness in comparison to conventional constant-temperature and multicanonical MD simulations. We found that the sampling efficiency of the modified REPSWA MD after a careful optimization of its parameters was better than that of the constant-temperature MD and comparable to that of the multicanonical MD in several cases.

Molecular dynamics sampling scheme realizing multiple distributions

We present a molecular-dynamics sampling scheme in which any summation of multiple arbitrarily given distributions can be realized deterministically by a certain theoretical guide. Our scheme thus provides a static view of the composition of the distributions, rather than a dynamic view in which some parameter is switched at a certain time in the simulation process. The proposed method induces the orbit to jump automatically among different phase space regions, without the use of any artificial timing parameter. In addition, the proposed method does not require an assumption that the density of states is smooth. We applied it to multiple Tsallis distributions and established a suitable series of parameter values for which the sum of the distributions allows broad sampling. Numerical simulations applied to fundamental models with multi-Tsallis distributions showed efficient sampling, characterized by an energy trajectory that was totally different from that associated with each single distribution.

Visualization of conformational distribution of short to medium size segments in globular proteins and identification of local structural motifs

Analysis of the conformational distribution of polypeptide segments in a conformational space is the first step for understanding a principle of structural diversity of proteins. Here, we present a statistical analysis of protein local structures based on interatomic C(alpha) distances. Using principal component analysis (PCA) on the intrasegment C(alpha)-C(alpha) atomic distances, the conformational space of protein segments, which we call the protein segment universe, has been visualized, and three essential coordinate axes, suitable for describing the universe, have been identified. Three essential axes specified radius of gyration, structural symmetry, and separation of hairpin structures from other structures. Among the segments of arbitrary length, 6-22 residues long, the conservation of those axes was uncovered. Further application of PCA to the two largest clusters in the universe revealed local structural motifs. Although some of motifs have already been reported, we identified a possibly novel strand motif. We also showed that a capping box, which is one of the helix capping motifs, was separated into independent subclusters based on the C(alpha) geometry. Implications of the strand motif, which may play a role for protein-protein interaction, are discussed. The currently proposed method is useful for not only mapping the immense universe of protein structures but also identification of structural motifs.

Free-energy landscape of a chameleon sequence in explicit water and its inherent alpha/beta bifacial property

A sequence in yeast MATalpha2/MCM1/DNA complex that folds into alpha-helix or beta-hairpin depending on the surroundings has been known as "chameleon" sequence. We obtained the free-energy landscape of this sequence by using a generalized-ensemble method, multicanonical molecular dynamics simulation, to sample the conformational space. The system was expressed with an all-atom model in explicit water, and the initial conformation for the simulation was a random one. The free-energy landscape demonstrated that this sequence inherently has an ability to form either alpha or beta structure: The conformational distribution in the landscape consisted of two alpha-helical clusters with different packing patterns of hydrophobic residues, and four beta-hairpin clusters with different strand-strand interaction patterns. Narrow pathways connecting the clusters were found, and analysis on the pathways showed that a compact structure formed at the N-terminal root of the chameleon sequence controls the cluster-cluster transitions. The free-energy landscape indicates that a small conditional change induces alpha-beta transitions. Additional unfolding simulations done with replacing amino acids showed that the chameleon sequence has an advantage to form an alpha-helix. Current study may be useful to understand the mechanism of diseases resulting from abnormal chain folding, such as amyloid disease.

The Solution Structure of a Chimeric LEKTI Domain Reveals a Chameleon Sequence

The conversion of an alpha-helical to a beta-strand conformation and the presence of chameleon sequences are fascinating from the perspective that such structural features are implicated in the induction of amyloid-related fatal diseases. In this study, we have determined the solution structure of a chimeric domain (Dom1PI) from the multidomain Kazal-type serine proteinase inhibitor LEKTI using multidimensional NMR spectroscopy. This chimeric protein was constructed to investigate the reasons for differences in the folds of the homologous LEKTI domains 1 and 6 [Lauber, T., et al. (2003) J. Mol. Biol. 328, 205-219]. In Dom1PI, two adjacent phenylalanine residues (F28 and F29) of domain 1 were substituted with proline and isoleucine, respectively, as found in the corresponding P4' and P5' positions of domain 6. The three-dimensional structure of Dom1PI is significantly different from the structure of domain 1 and closely resembles the structure of domain 6, despite the sequence being identical to that of domain 1 except for the two substituted phenylalanine residues and being only 31% identical to the sequence of domain 6. The mutation converted a short 3(10)-helix into an extended loop conformation and parts of the long COOH-terminal alpha-helix of domain 1 into a beta-hairpin structure. The latter conformational change occurs in a sequence stretch distinct from the region containing the substituted residues. Therefore, this switch from an alpha-helical structure to a beta-hairpin structure indicates a chameleon sequence of seven residues. We conclude that the secondary structure of Dom1PI is determined not only by the local protein sequence but also by nonlocal interactions.

Transition states for protein folding have native topologies despite high structural variability

We present a structural analysis of the folding transition states of three SH3 domains. Our results reveal that the secondary structure is not yet fully formed at this stage of folding and that the solvent is only partially excluded from the interior of the protein. Comparison of the members of the transition state ensemble with a database of native folds shows that, despite substantial local variability, the transition state structures can all be classified as having the topology characteristic of an SH3 domain. Our results suggest a mechanism for folding in which the formation of a network of interactions among a subset of hydrophobic residues ensures that the native topology is generated. Such a mechanism enables high fidelity in folding while minimizing the need to establish a large number of specific interactions in the conformational search.