Solution structure of a D,L-alternating oligonorleucine as a model of double-stranded antiparallel beta-helix

Conformational characteristics of alternating D,L linear peptides are of particular interest because of their capacity to form transmembrane channels with different transport properties, as some natural antibiotics do. Single- and double-stranded beta-helical structures are common for alternating D,L peptides. The stability of the beta-helix depends on several structural factors, such as the backbone peptide length, type and position of side chains, and nature of terminal groups. The NMR and molecular dynamics solution conformation of a synthetic alternating D,L-oligopeptide with 15 norleucines (XVMe) has been used as a model to get insight in to the conformational features of double-stranded beta-helix structures. The NH chemical shift values (delta(NH)) and long-range nuclear Overhauser effects (NOE) cross peaks, in particular interstrand connectivities, clearly point to an antiparallel double-stranded beta-helix for the XVMe major conformation in solution. An extensive set of distances (from NOE cross peaks) and H-bonds (from delta(NH)) has been included in the molecular dynamics calculations. The experimental NMR data and theoretical calculations clearly indicate that the most probable conformation of XVMe in solution is a double-strand antiparallel beta(5.6) increasing decreasing-helix structure.     

Simulation of an all-beta 3-icosapeptide containing the 20 proteinogenic side chains: effect of temperature, pH, counterions, solvent, and force field on helix stability

Simulations of various beta-peptides have in the last years clarified several issues concerning peptide folding equilibria and interpretation of experimental data, especially from NMR and CD spectroscopy. These simulations involved different temperatures, pH-values, ionic strengths, solvents, and force-field parameters, but a variation of these factors for one beta-peptide has not yet been done. To investigate the influence of varying these factors, we analyze the helix stability of an all-beta3-icosapeptide bearing all 20 proteinogenic amino acid side chains, which is experimentally observed to fold into a 3(14)-helix in methanol but not in water. Structural aspects, such as hydrogen-bonded rings and salt bridges, are discussed and a comparison with NMR primary (NOE distance bounds and 3J-values) and secondary (NMR derived model structures) data is made. We further investigate the reasons for the 3(14)-helix stability/instability in methanol/water. Of all factors studied, the presence of counterions seems to be the one inducing most significant effects in the simulations.     

Conformational analysis of endomorphin-2 by molecular dynamics methods

Endomorphin-2 (EM2, H-Tyr-Pro-Phe-Phe-NH(2)) is a highly potent and selective mu-opioid receptor agonist. A conformational analysis of EM2 was carried out by simulated annealing (SA) and molecular dynamics (MD) methods. Molecular modeling was conducted on both neutral (N-terminal NH(2)) and charged (N-terminal NH(3) (+)) molecules. Based on the results of NMR investigations showing an equilibrium mixture of cis and trans Tyr(1)-Pro(2) peptide bonds for EM2 in solution, simulations were performed with restrained cis-Pro and trans-Pro peptide bonds, too. A separate SA study with unrestrained Pro peptide bonds was also conducted. Preferred conformational states are presented in Ramachandran plots. The g(+), g(-), and trans populations of the aromatic amino acid residue side chains were determined in chi(1) space. The distances between the N-terminal N atom and the other backbone N and O atoms, and the distances between the centers of the aromatic rings and the Pro(2) ring, were determined. The energy distribution of the structures obtained by SA was calculated. The preferred secondary structural elements were different kinds of beta-turns, an inverse gamma-turn located in the N-terminal region, and regular and inverse gamma-turns located in the C-terminal region. These turns were stabilized by intramolecular H-bonds and bifurcated H-bonds.     

Conformational transition states of a beta-hairpin peptide between the ordered and disordered conformations in explicit water

The conformational transition states of a beta-hairpin peptide in explicit water were identified from the free energy landscapes obtained from the multicanonical ensemble, using an enhanced conformational sampling calculation. The beta-hairpin conformations were significant at 300 K in the landscape, and the typical nuclear Overhauser effect signals were reproduced, consistent with the previously reported experiment. In contrast, the disordered conformations were predominant at higher temperatures. Among the stable conformations at 300 K, there were several free energy barriers, which were not visible in the landscapes formed with the conventional parameters. We identified the transition states around the saddle points along the putative folding and unfolding paths between the beta-hairpin and the disordered conformations in the landscape. The characteristic features of these transition states are the predominant hydrophobic contacts and the several hydrogen bonds among the side-chains, as well as some of the backbone hydrogen bonds. The unfolding simulations at high temperatures, 400 K and 500 K, and their principal component analyses also provided estimates for the transition state conformations, which agreed well with those at 400 K and 500 K deduced from the current free energy landscapes at 400 K and 500 K, respectively. However, the transition states at high temperatures were much more widely distributed on the landscape than those at 300 K, and their conformations were different.     

Structural analysis of the N- and C-termini in a peptide with consensus sequence.

We present a structural analysis of a peptide, the sequence of which includes amino acids that show preferences for specific positions near the N- and C-termini in protein helices. This peptide has the sequence ac-YMSEDELKAAEAAFKRHGVP-amide, which includes a strong version of an N-terminal Harper-Rose capping box structure as well as a Gly located close to the C-terminus designed to elucidate its role in C-terminal capping. The sequence of five residues at the middle is inserted to separate effects at the two ends via a helix-stabilizing linker. Application of a simulated annealing procedure using interproton distance constraints derived from 1H NOESY experiments in water reveals the presence of a C-terminal structure in this model. The C-terminus forms a folded back structure in a significant fraction of structures generated by the annealing, in most of which Gly assumes an alpha L conformation. This structure occurs within a highly flexible region of the molecule and hence is occupied only a fraction of the time.     

Conformational flexibility and loss of structural rigidity for a model hexapeptide,GRGDTP: 1H‐NMR and molecular dynamics studies

The NMR and molecular dynamics methods are used to study the conformations of a hexapeptide, GRGDTP, which has been shown to be accessible to various types of cell‐adhesion based cellular behaviors such as cell‐to‐matrix interactions, cell differentiation, immunogenicity development, gene expression, angiogenesis, metastasis, sex determination and gamete fusion. ¹H‐NMR results indicate the existence of weak 5→2 hydrogen bonded β‐turn type‐III. Molecular simulation studies using a mixed protocol of distance geometry, constrained minimization, restrained molecular dynamics followed by energy minimization resulted additional conformations that include about 64% of population of inverse γ‐turn (HB, 3→1) and about 35% population of γ‐turn (HB, 4→2). The inter‐proton distances observed in γ‐and inverse γ‐turns are also consistent with the NMR constraints. The variable internal hydrogen bonding due to γ‐turns initiated at Gly 1 and Arg 2 , and its tendency to inter‐convert between γ‐and inverse γ‐turn conformations imply that the peptide is flexible in nature. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 460–471, 2013.     

Structure and hydration of the amylopectin trisaccharide building blocks--Synthesis, NMR, and molecular dynamics

To gain insight into the molecular details and hydration of amylopectin, the five constituting trisaccharides have been chemically synthesized as their methyl alpha-glycosides. All five trisaccharides were subjected to 950 MHz NMR spectroscopy for complete assignment and nanosecond molecular dynamics trajectories were calculated to study the structure and dynamics of the trisaccharides in aqueous solution. Systematic analysis of the simulation data revealed several examples of bridging water molecules playing an important role in the stabilization of specific amylopectin conformations, which was also supported by the experimental NMR data such as interresidue NOE's and heteronuclear scalar couplings between nuclei from neighboring residues. Although alpha-maltotriose, alpha-iso-maltotriose, alpha-panose and alpha-isopanose are relatively well characterized structures, the study also includes one less characterized trisaccharide with the structure alphaGlcp(1-->4)alphaGlcp(1-->6)alphaGlcp. This trisaccharide, tentatively labelled alpha-forkose, is located at the branch point of amylopectin, forking the amylopectin into two strands that align into double-helical segments. The results show that the conformation of alpha-forkose takes a natural bend form which fits well into the structure of the double-helical segment of amylopectin. As the only trisaccharide in this study the structure of alpha-forkose is not significantly influenced by the hydration. In contrast, alpha-isopanose takes a restricted, but rather extended form due to an exceptionally strong localized water density. The two homo-linkage oligomers, alpha-maltotriose and alpha-iso-maltotriose, showed to be the most extended and the most flexible trimers, respectively, providing regular structure for crystalline domains and maximum linker flexibility for amorphous domains.     

Determination of alpha-helix and beta-sheet stability in the solid state: a solid-state NMR investigation of poly(L-alanine)

The relative stability of alpha-helix and beta-sheet secondary structure in the solid state was investigated using poly(L-alanine) (PLA) as a model system. Protein folding and stability has been well studied in solution, but little is known about solid-state environments, such as the core of a folded protein, where peptide packing interactions are the dominant factor in determining structural stability. (13)C cross-polarization with magic angle spinning (CPMAS) NMR spectroscopy was used to determine the backbone conformation of solid powder samples of 15-kDa and 21.4-kDa PLA before and after various sample treatments. Reprecipitation from helix-inducing solvents traps the alpha-helical conformation of PLA, although the method of reprecipitation also affects the conformational distribution. Grinding converts the secondary structure of PLA to a final steady-state mixture of 55% beta-sheet and 45% alpha-helix at room temperature regardless of the initial secondary structure. Grinding PLA at liquid nitrogen temperatures leads to a similar steady-state mixture with 60% beta-sheet and 40% alpha-helix, indicating that mechanical shear force is sufficient to induce secondary structure interconversion. Cooling the sample in liquid nitrogen or subjecting it to high pressure has no effect on secondary structure. Heating the sample without grinding results in equilibration of secondary structure to 50% alpha-helix/50% beta-sheet at 100 degrees C when starting from a mostly alpha-helical state. No change was observed upon heating a beta-sheet sample, perhaps due to kinetic effects and the different heating rate used in the experiments. These results are consistent with beta-sheet approximately 260 J/mol more stable than alpha-helix in solid-state PLA.     

A molecular dynamics study of the stoichiometric complex formed by poly (alpha, L-glutamate) and octyltrimethylammonium ions in chloroform solution.

We present a molecular dynamics simulation at 300 K in explicit solvent environment of chloroform of the stoichiometric complex formed by poly(alpha,L-glutamate) and octyltrimethylammonium ions. We observed that the alpha-helix conformation of the polypeptide chain remains stable during a 2-ns run. The surfactant ions predominantly adopted an extended conformation that is stabilized by favorable interactions with the organic solvent. Analysis of the organization of the surfactant with respect to the polypeptide chain indicated that each octyltrimethylammonium cation was preferentially bound to more than one carboxylate group. It was found that the most populated arrangement was that with the surfactant cations interacting with two carboxylate groups simultaneously.     

Conserved and nonconserved features of the folding pathway of hisactophilin, a beta-trefoil protein

Based on previous studies of interleukin-1beta (IL-1beta) and both acidic and basic fibroblast growth factors (FGFs), it has been suggested that the folding of beta-trefoil proteins is intrinsically slow and may occur via the formation of essential intermediates. Using optical and NMR-detected quenched-flow hydrogen/deuterium exchange methods, we have measured the folding kinetics of hisactophilin, another beta-trefoil protein that has < 10% sequence identity and unrelated function to IL-1beta and FGFs. We find that hisactophilin can fold rapidly and with apparently two-state kinetics, except under the most stabilizing conditions investigated where there is evidence for formation of a folding intermediate. The hisactophilin intermediate has significant structural similarities to the IL-1beta intermediate that has been observed experimentally and predicted theoretically using a simple, topology-based folding model; however, it appears to be different from the folding intermediate observed experimentally for acidic FGF. For hisactophilin and acidic FGF, intermediates are much less prominent during folding than for IL-1beta. Considering the structures of the different beta-trefoil proteins, it appears that differences in nonconserved loops and hydrophobic interactions may play an important role in differential stabilization of the intermediates for these proteins.