Flow-Induced β-Hairpin Folding of the Glycoprotein Ibα β-Switch

Flow-induced shear has been identified as a regulatory driving force in blood clotting. Shear induces β-hairpin folding of the glycoprotein Ibα β-switch which increases affinity for binding to the von Willebrand factor, a key step in blood clot formation and wound healing. Through 2.1-μs molecular dynamics simulations, we investigate the kinetics of flow-induced β-hairpin folding. Simulations sampling different flow velocities reveal that under flow, β-hairpin folding is initiated by hydrophobic collapse, followed by interstrand hydrogen-bond formation and turn formation. Adaptive biasing force simulations are employed to determine the free energy required for extending the unfolded β-switch from a loop to an elongated state. Lattice and freely jointed chain models illustrate how the folding rate depends on the entropic and enthalpic energy, the latter controlled by flow. The results reveal that the free energy landscape of the β-switch has two stable conformations imprinted on it, namely, loop and hairpin—with flow inducing a transition between the two.     

Folding Simulations of a De Novo Designed Protein with a βαβ Fold

βαβ structural motifs are commonly used building blocks in protein structures containing parallel β-sheets. However, to our knowledge, no stand-alone βαβ structure has been observed in nature to date. Recently, for the first time that we know of, a small protein with an independent βαβ structure (DS119) was successfully designed in our laboratory. To understand the folding mechanism of DS119, in the study described here, we carried out all-atom molecular dynamics and coarse-grained simulations to investigate its folding pathways and energy landscape. From all-atom simulations, we successfully observed the folding event and got a stable folded structure with a minimal root mean-square deviation of 2.6 Å with respect to the NMR structure. The folding process can be described as a fast collapse phase followed by rapid formation of the central helix, and then slow formation of a parallel β-sheet. By using a native-centric Gō-like model, the cooperativity of the system was characterized in terms of the calorimetric criterion, sigmoidal transitions, conformation distribution shifts, and free-energy profiles. DS119 was found to be an incipient downhill folder that folds more cooperatively than a downhill folder, but less cooperatively than a two-state folder. This may reflect the balance between the two structural elements of DS119: the rapidly formed α-helix and the slowly formed parallel β-sheet. Folding times estimated from both the all-atom simulations and the coarse-grained model were at microsecond level, making DS119 another fast folder. Compared to fast folders reported previously, DS119 is, to the best of our knowledge, the first that exhibits a parallel β-sheet.     

Altered Backbone and Side-Chain Interactions Result in Route Heterogeneity during the Folding of Interleukin-1β (IL-1β)

Deletion of the β-bulge trigger-loop results in both a switch in the preferred folding route, from the functional loop packing folding route to barrel closure, as well as conversion of the agonist activity of IL-1β into antagonist activity. Conversely, circular permutations of IL-1β conserve the functional folding route as well as the agonist activity. These two extremes in the folding-functional interplay beg the question of whether mutations in IL-1β would result in changes in the populations of heterogeneous folding routes and the signaling activity. A series of topologically equivalent water-mediated β-strand bridging interactions within the pseudosymmetric β-trefoil fold of IL-1β highlight the backbone water interactions that stabilize the secondary and tertiary structure of the protein. Additionally, conserved aromatic residues lining the central cavity appear to be essential for both stability and folding. Here, we probe these protein backbone-water molecule and side chain-side chain interactions and the role they play in the folding mechanism of this geometrically stressed molecule. We used folding simulations with structure-based models, as well as a series of folding kinetic experiments to examine the effects of the F42W core mutation on the folding landscape of IL-1β. This mutation alters water-mediated backbone interactions essential for maintaining the trefoil fold. Our results clearly indicate that this perturbation in the primary structure alters a structural water interaction and consequently modulates the population of folding routes accessed during folding and signaling activity.     

Mechanism of Membrane Permeation Induced by Synthetic β-Hairpin Peptides

We have investigated the membrane destabilizing properties of synthetic amphiphilic cationic peptides, MAX1 and MAX35, which have the propensity to form β-hairpin structures under certain conditions, and a control non-β-hairpin-forming peptide MAX8V16E. All three peptides bind to liposomes containing a mixture of zwitterionic POPC and negatively charged POPS lipids as determined by Zeta potential measurements. Circular dichroism measurements indicated folding of MAX1 and MAX35 in the presence of the POPC/POPS liposomes, whereas no such folding was observed with MAX8V16E. There was no binding or folding of these peptides to liposomes containing only POPC. MAX1 and MAX35 induced release of contents from negatively charged liposomes, whereas MAX8V16E failed to promote solute release under identical conditions. Thus, MAX1 and MAX35 bind to, and fold at the surface of negatively charged liposomes adopting a lytic conformation. We ruled out leaky fusion as a mechanism of release by including 2 mol % PEG-PE in the liposomes, which inhibits aggregation/fusion but not folding of MAX or MAX-induced leakage. Using a concentration-dependent quenching probe (calcein), we determined that MAX-induced leakage of liposome contents was an all-or-none process. At MAX1 concentrations, which cause release of ∼50% of the liposomes that contain small (R_h <1.5 nm) markers, only ∼15% of those liposomes release a fluorescent dextran of 40 kDa. A multimeric model of the pore is presented based on these results. Atomistic molecular dynamics simulations show that barrels consisting of 10 β-hairpin MAX1 and MAX35 peptides are relatively more stable than MAX8V16E barrels in the bilayer, suggesting that barrels of this size are responsible for the peptides lytic action.     

The Structure of Aβ42 C-Terminal Fragments Probed by a Combined Experimental and Theoretical Study

The C-terminus of amyloid β-protein (Aβ) 42 plays an important role in this protein's oligomerization and may therefore be a good therapeutic target for the treatment of Alzheimer's disease. Certain C-terminal fragments (CTFs) of Aβ42 have been shown to disrupt oligomerization and to strongly inhibit Aβ42-induced neurotoxicity. Here we study the structures of selected CTFs [Aβ(x-42); x = 29-31, 39] using replica exchange molecular dynamics simulations and ion mobility mass spectrometry. Our simulations in explicit solvent reveal that the CTFs adopt a metastable β-structure: β-hairpin for Aβ(x-42) (x = 29-31) and extended β-strand for Aβ(39-42). The β-hairpin of Aβ(30-42) is converted into a turn-coil conformation when the last two hydrophobic residues are removed, suggesting that I41 and A42 are critical in stabilizing the β-hairpin in Aβ42-derived CTFs. The importance of solvent in determining the structure of the CTFs is further highlighted in ion mobility mass spectrometry experiments and solvent-free replica exchange molecular dynamics simulations. A comparison between structures with solvent and structures without solvent reveals that hydrophobic interactions are critical for the formation of β-hairpin. The possible role played by the CTFs in disrupting oligomerization is discussed.     

Synthesis and β-sheet propensity of constrained N-amino peptides

The stabilization of β-sheet secondary structure through peptide backbone modification represents an attractive approach to protein mimicry. Here, we present strategies toward stable β-hairpin folds based on peptide strand N-amination. Novel pyrazolidinone and tetrahydropyridazinone dipeptide constraints were introduced via on-resin Mitsunobu cyclization between α-hydrazino acid residues and a serine or homoserine side chain. Acyclic and cyclic N-amino peptide building blocks were then evaluated for their effect on β-hairpin stability in water using a GB1-derived model system. Our results demonstrate the strong β-sheet stabilizing effect of the peptide N-amino substituent, and provide useful insights into the impact of covalent dipeptide constraint on β-sheet folding.     

A Second Look at Mini-Protein Stability: Analysis of FSD-1 Using Circular Dichroism, Differential Scanning Calorimetry, and Simulations

Mini-proteins that contain <50 amino acids often serve as model systems for studying protein folding because their small size makes long timescale simulations possible. However, not all mini-proteins are created equal. The stability and structure of FSD-1, a 28-residue mini-protein that adopted the ββα zinc-finger motif independent of zinc binding, was investigated using circular dichroism, differential scanning calorimetry, and replica-exchange molecular dynamics. The broad melting transition of FSD-1, similar to that of a helix-to-coil transition, was observed by using circular dichroism, differential scanning calorimetry, and replica-exchange molecular dynamics. The N-terminal β-hairpin was found to be flexible. The FSD-1 apparent melting temperature of 41°C may be a reflection of the melting of its α-helical segment instead of the entire protein. Thus, despite its attractiveness due to small size and purposefully designed helix, sheet, and turn structures, the status of FSD-1 as a model system for studying protein folding should be reconsidered.     

Flow-induced structural transition in the beta-switch region of glycoprotein Ib

The impact of fluid flow on structure and dynamics of biomolecules has recently gained much attention. In this article, we present a molecular-dynamics algorithm that serves to generate stable water flow under constant temperature, for the study of flow-induced protein behavior. Flow simulations were performed on the 16-residue β-switch region of platelet glycoprotein Ibα, for which crystal structures of its N-terminal domain alone and in complex with the A1 domain of von Willebrand factor have been solved. Comparison of the two structures reveals a conformational change in this region, which, upon complex formation, switches from an unstructured loop to a β-hairpin. Interaction between glycoprotein Ibα and von Willebrand factor initiates platelet adhesion to injured vessel walls, and the adhesion is enhanced by blood flow. It has been hypothesized that the loop to β-hairpin transition in glycoprotein Ibα is induced by flow before binding to von Willebrand factor. The simulations revealed clearly a flow-induced loop→β-hairpin transition. The transition is dominated by the entropy of the protein, and is seen to occur in two steps, namely a dihedral rotation step followed by a side-group packing step.     

Nonadditivity in conformational entropy upon molecular rigidification reveals a universal mechanism affecting folding cooperativity

Previously, we employed a Maxwell counting distance constraint model (McDCM) to describe α-helix formation in polypeptides. Unlike classical helix-coil transition theories, the folding mechanism derives from nonadditivity in conformational entropy caused by rigidification of molecular structure as intramolecular cross-linking interactions form along the backbone. For example, when a hydrogen bond forms within a flexible region, both energy and conformational entropy decrease. However, no conformational entropy is lost when the region is already rigid because atomic motions are not constrained further. Unlike classical zipper models, the same mechanism also describes a coil-to-β-hairpin transition. Special topological features of the helix and hairpin structures allow the McDCM to be solved exactly. Taking full advantage of the fact that Maxwell constraint counting is a mean field approximation applied to the distribution of cross-linking interactions, we present an exact transfer matrix method that does not require any special topological feature. Upon application of the model to proteins, cooperativity within the folding transition is yet again appropriately described. Notwithstanding other contributing factors such as the hydrophobic effect, this simple model identifies a universal mechanism for cooperativity within polypeptide and protein-folding transitions, and it elucidates scaling laws describing hydrogen-bond patterns observed in secondary structure. In particular, the native state should have roughly twice as many constraints as there are degrees of freedom in the coil state to ensure high fidelity in two-state folding cooperativity, which is empirically observed.     

Free energy landscapes of two model peptides: α-helical and β-hairpin peptides explored with Brownian dynamics simulation

We applied an atomistic Brownian dynamics (BD) simulation with multiple time step method for the folding simulation of a 13-mer α-helical peptide and a 12-mer β-hairpin peptide, giving successful folding simulations. In this model, the driving energy contribution towards folding came from both electrostatic and van der Waals interactions for the α-helical peptide and from van der Waals interactions for the β-hairpin peptide. Although, many non-native structures having the same or lower energy than that of native structure were observed, the folded states formed the most populated cluster when the structures obtained by the BD simulations were subjected to the cluster analysis based on distance-based root mean square deviation of side-chains between different structures. This result indicates that we can predict the native structures from conformations sampled by BD simulation.     

Differences in β-strand Populations of Monomeric Aβ40 and Aβ42

Using homonuclear ¹H NOESY spectra, with chemical shifts, ³J_H^N_H^α scalar couplings, residual dipolar couplings, and ¹H-¹⁵N NOEs, we have optimized and validated the conformational ensembles of the amyloid-β 1–40 (Aβ40) and amyloid-β 1–42 (Aβ42) peptides generated by molecular dynamics simulations. We find that both peptides have a diverse set of secondary structure elements including turns, helices, and antiparallel and parallel β-strands. The most significant difference in the structural ensembles of the two peptides is the type of β-hairpins and β-strands they populate. We find that Aβ42 forms a major antiparallel β-hairpin involving the central hydrophobic cluster residues (16–21) with residues 29–36, compatible with known amyloid fibril forming regions, whereas Aβ40 forms an alternative but less populated antiparallel β-hairpin between the central hydrophobic cluster and residues 9–13, that sometimes forms a β-sheet by association with residues 35–37. Furthermore, we show that the two additional C-terminal residues of Aβ42, in particular Ile-41, directly control the differences in the β-strand content found between the Aβ40 and Aβ42 structural ensembles. Integrating the experimental and theoretical evidence accumulated over the last decade, it is now possible to present monomeric structural ensembles of Aβ40 and Aβ42 consistent with available information that produce a plausible molecular basis for why Aβ42 exhibits greater fibrillization rates than Aβ40.     

Stable and metastable states of human amylin in solution

Patients with type II diabetes exhibit fibrillar deposits of human amylin protein in the pancreas. It has been proposed that amylin oligomers arising along the aggregation or fibril-formation pathways are important in the genesis of the disease. In a step toward understanding these aggregation pathways, in this work we report the conformational preferences of human amylin monomer in solution using molecular simulations and infrared experiments. In particular, we identify a stable conformer that could play a key role in aggregation. We find that amylin adopts three stable conformations: one with an α-helical segment comprising residues 9-17 and a short antiparallel β-sheet comprising residues 24-28 and 31-35; one with an extended antiparallel β-hairpin with the turn region comprising residues 20-23; and one with no particular structure. Using detailed calculations, we determine the relative stability of these various conformations, finding that the β-hairpin conformation is the most stable, followed by the α-helical conformation, and then the unstructured coil. To test our predicted structure, we calculate its infrared spectrum in the amide I stretch regime, which is sensitive to secondary structure through vibrational couplings and linewidths, and compare it to experiment. We find that theoretically predicted spectra are in good agreement with the experimental line shapes presented herein. The implications of the monomer secondary structures on its aggregation pathway and on its interaction with cell membranes are discussed.     

Characterizing protein energy landscape by self-learning multiscale simulations: application to a designed β-hairpin

Characterizing the energy landscape of proteins at atomic resolution is still a very challenging problem, since it simultaneously requires high accuracy in estimating specific interactions and high efficiency in conformational sampling. Here, for these two requirements to meet, we extended the self-learning multiscale simulation (SLMS) method developed recently and applied it to the designed β-hairpin CLN025. The SLMS integrates all-atom and coarse-grained (CG) models in an iterative way such that the conformational sampling is performed by the CG model, the AA energy is used to calibrate the energy landscape, and the CG model is improved by the calibrated energy landscape. We extended the SLMS in two aspects, use of the energy decomposition for self-learning of the CG potential and a two-bead/residue CG model. The results show that the self-learning greatly improved the CG potential, and with the derived CG potential, the β-hairpin CLN025 robustly folded to the native structure. The self-learning iteration progressively enhanced the context dependence in the CG potential and increased the energy gap between the native and the denatured states of the CG model, leading to a funnel-like energy landscape. By using the SLMS method, without prior knowledge of the native structure but with the help of the AA energy, we can obtain a tailor-made CG potential specific to the target protein. The method can be useful for de novo structure prediction as well.     

Residual Structure of Unfolded Ubiquitin as Revealed by Hydrogen/Deuterium-Exchange 2D NMR

The characterization of residual structures persistent in unfolded proteins in concentrated denaturant solution is currently an important issue in studies of protein folding because the residual structure present, if any, in the unfolded state may form a folding initiation site and guide the subsequent folding reactions. Here, we studied the hydrogen/deuterium (H/D)-exchange behavior of unfolded human ubiquitin in 6 M guanidinium chloride. We employed a dimethylsulfoxide (DMSO)-quenched H/D-exchange NMR technique with the use of spin desalting columns, which allowed us to perform a quick medium exchange from 6 M guanidinium chloride to a quenching DMSO solution. Based on the backbone resonance assignment of ubiquitin in the DMSO solution, we successfully investigated the H/D-exchange kinetics of 60 identified peptide amide groups in the ubiquitin sequence. Although a majority of these amide groups were not protected, certain amide groups involved in a middle helix (residues 23–34) and an N-terminal β-hairpin (residues 2–16) were significantly protected with a protection factor of 2.1–4.2, indicating that there were residual structures in unfolded ubiquitin and that these amide groups were more than 52% hydrogen bonded in the residual structures. We show that the hydrogen-bonded residual structures in the α-helix and the β-hairpin are formed even in 6 M guanidinium chloride, suggesting that these residual structures may function as a folding initiation site to guide the subsequent folding reactions of ubiquitin.     

Smooth Functional Transition along a Mutational Pathway with an Abrupt Protein Fold Switch

Recent protein design experiments have demonstrated that proteins can migrate between folds through the accumulation of substitution mutations without visiting disordered or nonfunctional points in sequence space. To explore the biophysical mechanism underlying such transitions we use a three-letter continuous protein model with seven atoms per amino acid to provide realistic sequence-structure and sequence-function mappings through explicit simulation of the folding and interaction of model sequences. We start from two 16-amino-acid sequences folding into an α-helix and a β-hairpin, respectively, each of which has a preferred binding partner with 35 amino acids. We identify a mutational pathway between the two folds, which features a sharp fold switch. By contrast, we find that the transition in function is smooth. Moreover, the switch in preferred binding partner does not coincide with the fold switch. Discovery of new folds in evolution might therefore be facilitated by following fitness slopes in sequence space underpinned by binding-induced conformational switching.     

Antimicrobial Protegrin-1 Forms Ion Channels: Molecular Dynamic Simulation, Atomic Force Microscopy, and Electrical Conductance Studies

Antimicrobial peptides (AMPs) are an emerging class of antibiotics for controlling health effects of antibiotic-resistant microbial strains. Protegrin-1 (PG-1) is a model antibiotic among β-sheet AMPs. Antibiotic activity of AMPs involves cell membrane damage, yet their membrane interactions, their 3D membrane-associated structures and the mechanism underlying their ability to disrupt cell membrane are poorly understood. Using complementary approaches, including molecular dynamics simulations, atomic force microscopy (AFM) imaging, and planar lipid bilayer reconstitution, we provide computational and experimental evidence that PG-1, a β-hairpin peptide, forms ion channels. Simulations indicate that PG-1 forms channel-like structures with loosely attached subunits when reconstituted in anionic lipid bilayers. AFM images show the presence of channel-like structures when PG-1 is reconstituted in dioleoylphosphatidylserine/palmitoyloleoyl phosphatidylethanolamine bilayers or added to preformed bilayers. Planar lipid bilayer electrical recordings show multiple single channel conductances that are consistent with the heterogeneous oligomeric channel structures seen in AFM images. PG-1 channel formation seems to be lipid-dependent: PG-1 does not easily show ion channel electrical activity in phosphatidylcholine membranes, but readily shows channel activity in membranes rich in phosphatidylethanolamine or phosphatidylserine. The combined results support a model wherein the β-hairpin PG-1 peptide acts as an antibiotic by altering cell ionic homeostasis through ion channel formation in cell membranes.     

Effects of Topology and Sequence in Protein Folding Linked via Conformational Fluctuations

Experiments have compared the folding of proteins with different amino acid sequences but the same basic structure, or fold. Results indicate that folding is robust to sequence variations for proteins with some nonlocal folds, such as all-β, whereas the folding of more local, all-α proteins typically exhibits a stronger sequence dependence. Here, we use a coarse-grained model to systematically study how variations in sequence perturb the folding energy landscapes of three model sequences with 3α, 4β + α, and β-barrel folds, respectively. These three proteins exhibit folding features in line with experiments, including expected rank order in the cooperativity of the folding transition and stability-dependent shifts in the location of the free-energy barrier to folding. Using a generalized-ensemble simulation approach, we determine the thermodynamics of around 2000 sequence variants representing all possible hydrophobic or polar single- and double-point mutations. From an analysis of the subset of stability-neutral mutations, we find that folding is perturbed in a topology-dependent manner, with the β-barrel protein being the most robust. Our analysis shows, in particular, that the magnitude of mutational perturbations of the transition state is controlled in part by the size or “width” of the underlying conformational ensemble. This result suggests that the mutational robustness of the folding of the β-barrel protein is underpinned by its conformationally restricted transition state ensemble, revealing a link between sequence and topological effects in protein folding.     

Distinct unfolding and refolding pathways of ribonuclease a revealed by heating and cooling temperature jumps

Heating and cooling temperature jumps (T-jumps) were performed using a newly developed technique to trigger unfolding and refolding of wild-type ribonuclease A and a tryptophan-containing variant (Y115W). From the linear Arrhenius plots of the microscopic folding and unfolding rate constants, activation enthalpy (ΔH^#), and activation entropy (ΔS^#) were determined to characterize the kinetic transition states (TS) for the unfolding and refolding reactions. The single TS of the wild-type protein was split into three for the Y115W variant. Two of these transition states, TS1 and TS2, characterize a slow kinetic phase, and one, TS3, a fast phase. Heating T-jumps induced protein unfolding via TS2 and TS3; cooling T-jumps induced refolding via TS1 and TS3. The observed speed of the fast phase increased at lower temperature, due to a strongly negative ΔH^# of the folding-rate constant. The results are consistent with a path-dependent protein folding/unfolding mechanism. TS1 and TS2 are likely to reflect X-Pro¹¹⁴ isomerization in the folded and unfolded protein, respectively, and TS3 the local conformational change of the β-hairpin comprising Trp¹¹⁵. A very fast protein folding/unfolding phase appears to precede both processes. The path dependence of the observed kinetics is suggestive of a rugged energy protein folding funnel.     

Proline-glutamate chimera’s side chain conformation directs the type of β-hairpin structure

Our aim was to study the impact of two proline chimeras, containing a glutamic acid side chain in cis- or trans-configuration, on secondary structure formation. We further investigated to what extent the configuration of the side chain contributes to the overall peptide conformation. We used a 10 residue peptide (IYSNPDGTWT) that forms a β-hairpin in water. The turn-forming proline was substituted with either a cis- or trans-proline-glutamic acid chimera, resulting in the peptides IYSNP ^{cis -E} DGTWT (P1_P ^cis-E ) and IYSNP ^{trans -E} DGTWT (P1_P ^trans-E ). We studied the conformation of the modified peptides by circular dichroism (CD) and NMR-spectroscopy, and SEC/static light scattering (SLS) analysis. NMR analysis reveals that the modified peptides maintain the β-hairpin conformation in aqueous solution. At 5 °C and pH 4.3, the peptide (P1_P ^cis-E ) was found to adopt two coexisting β-hairpin conformations (2:2 β-hairpin, and 3:5 β-hairpin). In contrast to that, the peptide (P1_P ^trans-E ) adopts a 2:2 β-hairpin that exists in equilibrium with a 4:4 β-hairpin conformation. The adoption of ordered β-hairpin structures for both modified peptides could be confirmed by CD spectroscopy, while SEC/SLS analysis showed a monomeric oligomerization state for all three investigated peptides. With the combination of several NMR methods, we were able to elucidate that even small alterations in the side chain conformation of the proline-glutamate chimera (cis or trans) can significantly influence the conformation of the adopted β-hairpin.