Dataset Information

The impact of interchain hydrogen bonding on ?-hairpin stability is readily predicted by molecular dynamics simulation.

ABSTRACT: Peptides are frequently used model systems for protein folding. They are also gaining increased importance as therapeutics. Here, the ability of molecular dynamics (MD) simulation for describing the structure and dynamics of ?-hairpin peptides was investigated, with special attention given to the impact of a single interstrand sidechain to sidechain interaction. The MD trajectories were compared to structural information gained from solution NMR. By assigning frames from restraint-free MD simulations to an intuitive hydrogen bond on/off pattern, folding ratios and folding pathways were predicted. The computed molecular model successfully reproduces the folding ratios determined by NMR, indicating that MD simulation may be straightforwardly used as a screening tool in ?-hairpin design.

SUBMITTER: Niebling S

PROVIDER: S-EPMC4744756 | biostudies-literature | 2015 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

The impact of interchain hydrogen bonding on β-hairpin stability is readily predicted by molecular dynamics simulation.

Niebling Stephan S Danelius Emma E Brath Ulrika U Westenhoff Sebastian S Erdélyi Máté M

Biopolymers 20151101 6

Peptides are frequently used model systems for protein folding. They are also gaining increased importance as therapeutics. Here, the ability of molecular dynamics (MD) simulation for describing the structure and dynamics of β-hairpin peptides was investigated, with special attention given to the impact of a single interstrand sidechain to sidechain interaction. The MD trajectories were compared to structural information gained from solution NMR. By assigning frames from restraint-free MD simula ...[more]

PMID: 25968880

Similar Datasets

Project description:Hydrogen exchange (HX) studies have provided critical insight into our understanding of protein folding, structure, and dynamics. More recently, hydrogen exchange mass spectrometry (HX-MS) has become a widely applicable tool for HX studies. The interpretation of the wealth of data generated by HX-MS experiments as well as other HX methods would greatly benefit from the availability of exchange predictions derived from structures or models for comparison with experiment. Most reported computational HX modeling studies have employed solvent-accessible-surface-area based metrics in attempts to interpret HX data on the basis of structures or models. In this study, a computational HX-MS prediction method based on classification of the amide hydrogen bonding modes mimicking the local unfolding model is demonstrated. Analysis of the NH bonding configurations from molecular dynamics (MD) simulation snapshots is used to determine partitioning over bonded and nonbonded NH states and is directly mapped into a protection factor (PF) using a logistics growth function. Predicted PFs are then used for calculating deuteration values of peptides and compared with experimental data. Hydrogen exchange MS data for fatty acid synthase thioesterase (FAS-TE) collected for a range of pHs and temperatures was used for detailed evaluation of the approach. High correlation between prediction and experiment for observable fragment peptides is observed in the FAS-TE and additional benchmarking systems that included various apo/holo proteins for which literature data were available. In addition, it is shown that HX modeling can improve experimental resolution through decomposition of in-exchange curves into rate classes, which correlate with prediction from MD. Successful rate class decompositions provide further evidence that the presented approach captures the underlying physical processes correctly at the single residue level. This assessment is further strengthened in a comparison of residue resolved protection factor predictions for staphylococcal nuclease with NMR data, which was also used to compare prediction performance with other algorithms described in the literature. The demonstrated transferable and scalable MD based HX prediction approach adds significantly to the available tools for HX-MS data interpretation based on available structures and models.

Project description:All atom molecular dynamics simulations of the 18-residue beta-hairpin antimicrobial peptide protegrin-1 (PG-1, RGGRLCYCRRRFCVCVGR-NH(2)) in a fully hydrated dilauroylphosphatidylcholine (DLPC) lipid bilayer have been implemented. The goal of the reported work is to investigate the structure of the peptide in a membrane environment (previously solved only in solution [R.L. Fahrner, T. Dieckmann, S.S.L. Harwig, R.I. Lehrer, D. Eisenberg, J. Feigon, Solution structure of protegrin-1, a broad-spectrum antimicrobial peptide from porcine leukocytes. Chemistry and Biology, 3 (1996) 543-550]), and to delineate specific peptide-membrane interactions which are responsible for the peptide's membrane binding properties. A novel, previously unknown, "kick" shaped conformation of the peptide was detected, where a bend at the C-terminal beta-strand of the peptide caused the peptide backbone at residues 16-18 to extend perpendicular to the beta-hairpin plane. This bend was driven by a highly persistent hydrogen-bond between the polar peptide side-chain of TYR7 and the unshielded backbone carbonyl oxygen atom of GLY17. The H-bond formation relieves the unfavorable free energy of insertion of polar groups into the hydrophobic membrane core. PG-1 was anchored to the membrane by strong electrostatic binding of the protonated N-terminus of the peptide to the lipid head group phosphate anions. The orientation of the peptide in the membrane, and its influence on bilayer structural and dynamic properties are in excellent agreement with solid state NMR measurements [S. Yamaguchi, T. Hong, A. Waring, R.I. Lehrer, M. Hong, Solid-State NMR Investigations of Peptide-Lipid Interaction and Orientation of a b-Sheet Antimicrobial Peptide, Protegrin, Biochemistry, 41 (2002) 9852-9862]. Importantly, two simulations which started from different initial orientations of the peptide converged to the same final equilibrium orientation of the peptide relative to the bilayer. The kick-shaped conformation was observed only in one of the two simulations.

Dataset Information

The impact of interchain hydrogen bonding on ?-hairpin stability is readily predicted by molecular dynamics simulation.

Publications

The impact of interchain hydrogen bonding on β-hairpin stability is readily predicted by molecular dynamics simulation.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets