Comparison of the structural characteristics of Cu2+-bound and unbound α-syn12 peptide obtained in simulations using different force fields

The effects of Cu²⁺ binding and the utilization of different force fields when modeling the structural characteristics of α-syn12 peptide were investigated. To this end, we performed extensive temperature replica exchange molecular dynamics (T-REMD) simulations on Cu²⁺-bound and unbound α-syn12 peptide using the GROMOS 43A1, OPLS-AA, and AMBER03 force fields. Each replica was run for 300 ns. The structural characteristics of α-syn12 peptide were studied based on backbone dihedral angle distributions, free-energy surfaces obtained with different reaction coordinates, favored conformations, the formation of different Turn structures, and the solvent exposure of the hydrophobic residues. The findings show that AMBER03 prefers to sample helical structures for the unbound α-syn12 peptide and does not sample any β-hairpin structure for the Cu²⁺-bound α-syn12 peptide. In contrast, the central structure of the major conformational clusters for the Cu²⁺-bound and unbound α-syn12 peptide according to simulations performed using the GROMOS 43A1 and OPLS-AA force fields is a β-hairpin with Turn_9-6. Cu²⁺ can also promote the formation of the β-hairpin and increase the solvent exposure of hydrophobic residues, which promotes the aggregation of α-syn12 peptide. This study can help us to understand the mechanisms through which Cu²⁺ participates in the fibrillation of α-syn12 peptide at the atomic level, which in turn represents a step towards elucidating the nosogenesis of Parkinson’s disease.

Effects of pH and temperature on the structural and thermodynamic character of alpha-syn12 peptide in aqueous solution

The structural and thermodynamic characters of alpha-syn12 peptide in aqueous solution at different pH and temperatures have been investigated through temperature replica exchange molecular dynamics (T-REMD) simulations with GROMOS 43A1 force field. The two independent T-REMD simulations were completed at pH = 7.0 and 10.0, respectively. Each replica was run for 300 ns. The structural and thermodynamic characters of alpha syn12 peptide were studied based on the distributions of backbone dihedral angles, the free energy surface, and the stability of different type structure and the favorite conformations of the peptide. The results showed that the simulation at pH = 10.0 produced more sampling in alpha region than the simulation at pH = 7.0. The temperature changes from 283 K to 308 K result in negligible effects on the distributions of backbone dihedral angle. The beta hairpin conformation with Turn(9-6) and four hydrogen bonds (HB(4-11), HB(6-9), HB(9-6) and HB(11-4)) is the lowest free energy state in the simulation at pH = 7.0. However, for the simulation at pH = 10.0, the lowest free energy state corresponds to a structure with Turn(9-6) and two hydrogen bonds (HB(6-10) and HB(10-6)) induced by an overly strong residue-residue interaction effect between lysine residues. For the simulation at pH = 7.0, the free energy change of the alpha-syn12 peptide from the unfolded state to the beta hairpin state was in good agreement with the experiments and molecular dynamics simulation results for the other beta-peptides, the beta hairpin state of the alpha-syn12 peptide included the conformations that not only the Turn(9-6) is formed, but also the terminus are closed together in space. However, the subtle balances between lysine-lysine interactions and lysine-solvent interaction are disrupted in the simulation at pH = 10.0, which induced the assembly of lysine residues, the beta hairpin conformation is destabilized by the deprotonation of the Lys side chain. This study can help us to understand the conformation changes and the thermodynamic character of alpha;-syn12 peptide at atomic level induced by changing pH and temperature, which is propitious to reveal the nosogenesis of Parkinson disease. In our knowledge, this is the first report to study the influence of pH and temperature on isolated alpha-syn12 peptide in water by T-REMD.     

Why the OPLS-AA force field cannot produce the β-hairpin structure of H1 peptide in solution when comparing with the GROMOS 43A1 force field?

The optimal combination of force field and water model is an essential problem that is able to increase molecular dynamics simulation quality for different types of proteins and peptides. In this work, an attempt has been made to explore the problem by studying H1 peptide using four different models based on different force fields, water models and electrostatic schemes. The driving force for H1 peptide conformation transition and the reason why the OPLS-AA force field cannot produce the β-hairpin structure of H1 peptide in solution while the GROMOS 43A1 force field can do were investigated by temperature replica exchange molecular dynamics simulation (T-REMD). The simulation using the GROMOS 43A1 force field preferred to adopt a β-hairpin structure, which was in good agreement with the several other simulations and the experimental evidences. However, the simulation using the OPLS-AA force field has a significant difference from the simulations with the GROMOS 43A1 force field simulation. The results show that the driving force in H1 peptide conformation transition is solvent exposure of its hydrophobic residues. However, the subtle balances between residue-residue interactions and residue-solvent interaction are disrupted by using the OPLS-AA force field, which induced the reduction in the number of residue-residue contact. Similar solvent exposure of the hydrophobic residues is observed for all the conformations sampled using the OPLS-AA force field. For H1 peptide which exhibits large solvent exposure of the hydrophobic residues, the GROMOS 43A1 force field with the SPC water model can provide more accurate results.     

Effects of different force fields on the structural character of α synuclein β-hairpin peptide (35–56) in aqueous environment

The hallmark of Parkinson’s disease (PD) is the intracellular protein aggregation forming Lewy Bodies (LB) and Lewy neuritis which comprise mostly of a protein, alpha synuclein (α-syn). Molecular dynamics (MD) simulation methods can augment experimental techniques to understand misfolding and aggregation pathways with atomistic resolution. The quality of MD simulations for proteins and peptides depends greatly on the accuracy of empirical force fields. The aim of this work is to investigate the effects of different force fields on the structural character of β hairpin fragment of α-syn (residues 35–56) peptide in aqueous solution. Six independent MD simulations are done in explicit solvent using, AMBER03, AMBER99SB, GROMOS96 43A1, GROMOS96 53A6, OPLS-AA, and CHARMM27 force fields with CMAP corrections. The performance of each force field is assessed from several structural parameters such as root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent accessible surface area (SASA), formation of β-turn, the stability of folded β-hairpin structure, and the favourable conformations obtained for different force fields. In this study, CMAP correction of CHARMM27 force field is found to overestimate the helical conformation, while GROMOS96 53A6 is found to most successfully capture the conformational dynamics of α-syn β-hairpin fragment as elicited from NMR.     

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.     

The role of alpha-, 3(10)-, and pi-helix in helix-->coil transitions

The conformational equilibrium between 3(10)- and alpha-helical structure has been studied via high-resolution NMR spectroscopy by Millhauser and coworkers using the MW peptide Ac-AMAAKAWAAKA AAARA-NH2. Their 750-MHz nuclear Overhauser effect spectroscopy (NOESY) spectra were interpreted to reflect appreciable populations of 3(10)-helix throughout the peptide, with the greatest contribution at the N and C termini. The presence of simultaneous alphaN(i,i + 2) and alphaN(i,i + 4) NOE cross-peaks was proposed to represent conformational averaging between 3(10)- and alpha-helical structures. In this study, we describe 25-nsec molecular dynamics simulations of the MW peptide at 298 K, using both an 8 A and a 10 A force-shifted nonbonded cutoff. The ensemble averages of both simulations are in reasonable agreement with the experimental helical content from circular dichroism (CD), the (3)J(HNalpha) coupling constants, and the 57 observed NOEs. Analysis of the structures from both simulations revealed very little formation of contiguous i --> i + 3 hydrogen bonds (3(10)-helix); however, there was a large population of bifurcated i --> i + 3 and i --> i + 4 alpha-helical hydrogen bonds. In addition, both simulations contained considerable populations of pi-helix (i --> i + 5 hydrogen bonds). Individual turns formed over residues 1-9, which we predict contribute to the intensities of the experimentally observed alphaN(i,i + 2) NOEs. Here we show how sampling of both folded and unfolded structures can provide a structural framework for deconvolution of the conformational contributions to experimental ensemble averages.     

Solution structure of kurtoxin: a gating modifier selective for Cav3 voltage-gated Ca(2+) channels

Kurtoxin is a 63-amino acid polypeptide isolated from the venom of the South African scorpion Parabuthus transvaalicus. It is the first and only peptide ligand known to interact with Cav3 (T-type) voltage-gated Ca(2+) channels with high affinity and to modify the voltage-dependent gating of these channels. Here we describe the nuclear magnetic resonance (NMR) solution structure of kurtoxin determined using two- and three-dimensional NMR spectroscopy with dynamical simulated annealing calculations. The molecular structure of the toxin was highly similar to those of scorpion α-toxins and contained an α-helix, three β-strands, and several turns stabilized by four disulfide bonds. This so-called "cysteine-stabilized α-helix and β-sheet (CSαβ)" motif is found in a number of functionally varied small proteins. A detailed comparison of the backbone structure of kurtoxin with those of the scorpion α-toxins revealed that three regions [first long loop (Asp(8)-Ile(15)), β-hairpin loop (Gly(39)-Leu(42)), and C-terminal segment (Arg(57)-Ala(63))] in kurtoxin significantly differ from the corresponding regions in scorpion α-toxins, suggesting that these regions may be important for interacting with Cav3 (T-type) Ca(2+) channels. In addition, the surface profile of kurtoxin shows a larger and more focused electropositive patch along with a larger hydrophobic surface compared to those seen on scorpion α-toxins. These distinct surface properties of kurtoxin could explain its binding to Cav3 (T-type) voltage-gated Ca(2+) channels.     

A mechanism by which binding of the broadly neutralizing antibody b12 unfolds the inner domain α1 helix in an engineered HIV-1 gp120

Using all-atom simulations, we examine the role of the I109C/Q428C disulfide "stitch" in altering the conformational distribution of engineered HIV-1 gp120 core relevant for binding of the broadly neutralizing recombinant antibody b12. In particular, we propose that the I109C/Q428C stitch results in a conformational distribution favoring an unfolded inner-domain α1-helix upon binding of b12. Using targeted molecular dynamics, we show that folded α1 in the b12-bound conformation of gp120 is stable both with and without the stitch, but that with folded α1, the stitch requires an orientation of the β20/β21 sheet that is sterically incompatible with b12 binding. Forcing β20/β21 into the orientation displayed by the b12-bound conformation after folding α1 with the stitch intact results in partial unfolding of α1, whereas without the stitch, β20/β21 reorientation does not affect the conformation of α1. These findings collectively support the hypothesis that the disulfide stitch shifts the conformational distribution of α1 to the unfolded state, meaning an unfolded α1 is not a strict requirement of the b12-bound conformational ensemble of gp120's lacking the I109C/Q428C stitch.     

Amyloid beta-protein monomer folding: free-energy surfaces reveal alloform-specific differences

Alloform-specific differences in structural dynamics between amyloid β-protein (Aβ) 40 and Aβ42 appear to underlie the pathogenesis of Alzheimer's disease. To elucidate these differences, we performed microsecond timescale replica-exchange molecular dynamics simulations to sample the conformational space of the Aβ monomer and constructed its free-energy surface. We find that neither peptide monomer is unstructured, but rather that each may be described as a unique statistical coil in which five relatively independent folding units exist, comprising residues 1-5, 10-13, 17-22, 28-37, and 39-42, which are connected by four turn structures. The free-energy surfaces of both peptides are characterized by two large basins, comprising conformers with either substantial α-helix or β-sheet content. Conformational transitions within and between these basins are rapid. The two additional hydrophobic residues at the Aβ42 C-terminus, Ile41 and Ala42, significantly increase contacts within the C-terminus, and between the C-terminus and the central hydrophobic cluster (Leu17-Ala21). As a result, the β-structure of Aβ42 is more stable than that of Aβ40, and the conformational equilibrium in Aβ42 shifts towards β-structure. These results suggest that drugs stabilizing α-helical Aβ conformers (or destabilizing the β-sheet state) would block formation of neurotoxic oligomers. The atomic-resolution conformer structures determined in our simulations may serve as useful targets for this purpose. The conformers also provide starting points for simulations of Aβ oligomerization—a process postulated to be the key pathogenetic event in Alzheimer's disease.     

Role of β-Hairpin Formation in Aggregation: The Self-Assembly of the Amyloid-β(25-35) Peptide

The amyloid-β(25-35) peptide plays a key role in the etiology of Alzheimer's disease due to its extreme toxicity even in the absence of aging. Because of its high tendency to aggregate and its low solubility in water, the structure of this peptide is still unknown. In this work, we sought to understand the early stages of aggregation of the amyloid-β(25-35) peptide by conducting simulations of oligomers ranging from monomers to tetramers. Our simulations show that although the monomer preferentially adopts a β-hairpin conformation, larger aggregates have extended structures, and a clear transition from compact β-hairpin conformations to extended β-strand structures occurs between dimers and trimers. Even though β-hairpins are not present in the final architecture of the fibril, our simulations indicate that they play a critical role in fibril growth. Our simulations also show that β-sheet structures are stabilized when a β-hairpin is present at the edge of the sheet. The binding of the hairpin to the sheet leads to a subsequent destabilization of the hairpin, with part of the hairpin backbone dangling in solution. This free section of the peptide can then recruit an extra monomer from solution, leading to further sheet extension. Our simulations indicate that the peptide must possess sufficient conformational flexibility to switch between a hairpin and an extended conformation in order for β-sheet extension to occur, and offer a rationalization for the experimental observation that overstabilizing a hairpin conformation in the monomeric state (for example, through chemical cross-linking) significantly hampers the fibrillization process.     

Role of cholesterol in functional association between K(+)-Cl(-) cotransporter-3a and Na+,K(+)-ATPase

The conformation of amyloid-beta peptide (Aβ) determines if toxic aggregates are formed. The peptide structure by its turn depends on the environment and molecule-molecule interactions. We characterized the secondary structure of Aβ-(1-40) in surfactant solutions and interacting with monolayers. The peptide adopts β-sheet structure in solutions of ionic surfactants at sub-micelle concentrations and α-helix in the presence of ionic micelles. Uncharged micelles induce β-sheets. Aβ-(1-40) alters the critical micelle concentration value of the non-ionic surfactant, underlining hydrophobic interactions. At ionic monolayers the peptide forms β-sheets when its concentration at the surface is high enough. These results suggest that only electrostatic interactions of charged micelles that surround completely the peptide are able to induce non-aggregated α-helix structure.     

Molecular dynamics simulations of a bacterial autotransporter: NalP from Neisseria meningitidis

NalP is an autotransporter secretory protein found in the outer membrane of Neisseria meningitidis. The crystal structure of the NalP translocator domain revealed a transmembrane β-barrel containing a central α-helix. The role of this α-helix, and of the conformational dynamics of the β-barrel pore have been studied via atomistic molecular dynamics simulations. Three simulations, each of 10 ns duration, of NalP embedded within a solvated DMPC bilayer were performed. The helix was removed from the barrel interior in one simulation. The conformational stability of the protein is similar to that of other outer membrane proteins, e.g., OmpA, in comparable simulations. The transmembrane β-barrel is stable even in the absence of the α-helix. Removal of the helix results in an influx of water into the pore region, suggesting the helix acts as a ‘plug’. Water molecules entering the resultant pore form hydrogen bonds with the barrel lining that compensate for the loss of helix-barrel hydrogen bonds. The dimensions of the pore fluctuate over the course of the simulation revealing it to be flexible, but only wide enough to allow transport of the passenger domain in an unfolded or extended conformation. The simulations help us to understand the role of the central helix in plugging the pore and in maintaining the width of the barrel, and show that the NalP monomer is sufficient for the transport of the passenger domain in an unfolded or extended conformation.     

Design,synthesis, and conformation of a model peptide of endothelin with cystine-stabilized α-helix motif

A model 16-peptide of endothelin-1 (MET-1), which has the minimized sequence homology to the corresponding pan of endothelin-1 (ET-1), was designed to confirm the cystine-stabilized α-helix motif. The model structure consists of an extended structure, a β-turn part, and an α-helix structure that is stabilized by two disulfide bonds. The α-helix segment was designed to emphasize the amphiphilic nature. In order to combine the extended structure and the α-helix segment, a D -Ala-Pro sequence was selected to fix the β-turn. The model endothelin 16-peptide amide was synthesized by solid-phase synthesis on a 4-methylbenzhydrylamine resin. Its conformation was examined by CD and two-dimensional (2D) ¹H-nmr measurements. MET-1 showed similar CD patterns to ET-1 in both buffer and 50% aqueous trifluoroethanol solution. The 2D nmr experiments in 50% aqueous ethylene glycol revealed that MET-1 closely resembles the conformation of ET-1 with an extended structure, an α-helix, and a β-turn unit in the same position of the sequence. Furthermore, model peptides without disulfide bond(s) could not assume a stable structure in aqueous solution, while they did have similar α-helical content in 50% trifluoroethanol with MET-1. When the two disulfide bridges were simultaneously formed, the peptide with the correct disulfide bonds (MET-1) was obtained in threefold excess to the isomer (apamin type. MET-2). These findings obtained by the modeling of ET-1 showed an important role for the stabilization of peptide conformation with disulfide bonds. © 1994 John Wiley & Sons, Inc.     

Ultrafast folding and molecular dynamics of a linear hydrophobic β-hairpin

The first computational study of the folding and dynamics of a hydrophobic β-hairpin containing a central heterochiral diproline segment is reported. Linear hydrophobic sequences containing centrally positioned diproline motifs, heterochiral (DL/LD) and homochiral (LL/DD)), are investigated for their ability to form β-hairpins. Heterochiral diproline motifs (LD/DL) reveal the formation of stable β-hairpins with the backbone adopting β-turn conformation and the formation of backbone hydrogen bonds with antiparallel cross-strand registry, whereas the homochiral diproline (LL/DD) containing sequences tend to adopt PP_II helix conformation. The competition between the β-turn formation and the backbone H-bond ladder of the antiparallel β-strands in heterochiral diproline containing sequences is employed to validate the hypothesis that β-turn formation precedes inter-strand registry in the folding of a β-hairpin (“zipper” mechanism). The observation of noncanonical hydrogen bonds leads to a folded β-hairpin-like conformation and points to the existence of relatively stable transition state intermediates, between the unfolded (extended) and folded (β-hairpin) states. The MD simulations are in excellent agreement with the experimental studies on the model system and constitute the very first computational investigation of the folding and dynamics of a completely hydrophobic synthetic β-hairpin containing heterogeneous residues of mixed chirality.     

Molecular Origin of Gerstmann-Sträussler-Scheinker Syndrome: Insight from Computer Simulation of an Amyloidogenic Prion Peptide

Prion proteins become pathogenic through misfolding. Here, we characterize the folding of a peptide consisting of residues 109–122 of the Syrian hamster prion protein (the H1 peptide) and of a more amyloidogenic A117V point mutant that leads in humans to an inheritable form of the Gerstmann-Sträussler-Scheinker syndrome. Atomistic molecular dynamics simulations are performed for 2.5 μs. Both peptides lose their α-helical starting conformations and assume a β-hairpin that is structurally similar in both systems. In each simulation several unfolding/refolding events occur, leading to convergence of the thermodynamics of the conformational states to within 1 kJ/mol. The similar stability of the β-hairpin relative to the unfolded state is observed in the two peptides. However, substantial differences are found between the two unfolded states. A local minimum is found within the free energy unfolded basin of the A117V mutant populated by misfolded collapsed conformations of comparable stability to the β-hairpin state, consistent with increased amyloidogenicity. This population, in which V117 stabilizes a hydrophobic core, is absent in the wild-type peptide. These results are supported by simulations of oligomers showing a slightly higher stability of the associated structures and a lower barrier to association for the mutated peptide. Hence, a single point mutation carrying only two additional methyl groups is here shown to be responsible for rather dramatic differences of structuring within the unfolded (misfolded) state.     

Simulation Studies on the Stabilities of Aggregates Formed by Fibril-Forming Segments of α-Synuclein

Abstract

We performed molecular dynamics simulations for various oligomers with different β-sheet conformations consisting of α-Synuclein 71–82 residues using an all atom force field and explicit water model. Tetramers of antiparallel β-sheet are shown to be stable, whereas parallel sheets are highly unstable due to the repulsive interactions between bulky and polar side chains as well as the weaker backbone hydrogen bonds. We also investigated the stabilities of double antiparallel β-sheets stacked with asymmetric and symmetric geometries. Our results show that this 12 amino acid residue peptide can form stable β-sheet conformers at 320K and higher temperatures. The backbone hydrogen bonds in β-sheet and the steric packing between hydrophobic side chains between β-sheets are shown to give conformational stabilities.     

Smoothing of the GB1 Hairpin Folding Landscape by Interfacial Confinement

We study the effects of confinement between planar walls on the folding thermodynamics of a β-hairpin, using large-scale replica-exchange molecular-dynamics simulations with an all-atom model and explicit solvent. We find that the folding free-energy landscape of this peptide observed in bulk is significantly modified when the peptide is confined between the walls. Most notably, the propensity of the peptide to form a misfolded state observed in the bulk solution becomes negligible under confinement. The absence of the misfolded state under confinement can be explained by an increased tendency of hydrophobic aromatic side chains to stay near the walls, because the misfolded state is characterized by a nonnative arrangement of aromatic side chains. These results from a simple confinement model may provide clues about the role of chaperonin confinement in smoothing folding landscapes by avoiding trapped intermediates.     

Flow-induced conformational changes and phase behavior of aqueous poly-L-lysine solutions

Poly-L -lysine exists as an α-helix at high pH and a random coil at neutral pH. When the α-helix is heated above 27°C, the macromolecule undergoes a conformational transition to a β-sheet. In this study, the stability of the secondary structure of poly-L -lysine in solutions subjected to shear flow, at temperatures below the α-helix to β-sheet transition temperature, were examined using Raman spectroscopy and CD. Solutions initially in the α-helical state showed time-dependent increases in viscosity with shearing, rising as much as an order of magnitude. Visual observation and turbidity measurements showed the formation of a gel-like phase under flow. Laser Raman measurements demonstrated the presence of small amounts of β-sheet structure evidenced by the amide I band at 1666 cm⁻¹. CD measurements indicated that solutions of predominantly α-helical conformation at 20°C transformed into 85% α-helix and 15% β-sheet after being sheared for 20 min. However, on continued shearing the content of β-sheet conformation decreased. The observed phenomena were explained in terms of a “zipping-up” molecular model based on flow enhanced hydrophobic interactions similar to that observed in gel-forming flexible polymers. © 1998 John Wiley & Sons, Inc. Biopoly 45: 239–246, 1998     

Conformational choreography of a molecular switch region in myelin basic protein--molecular dynamics shows induced folding and secondary structure type conversion upon threonyl phosphorylation in both aqueous and membrane-associated environments

The 18.5 kDa isoform of myelin basic protein is essential to maintaining the close apposition of myelin membranes in central nervous system myelin, but its intrinsic disorder (conformational dependence on environment), a variety of post-translational modifications, and a diversity of protein ligands (e.g., actin and tubulin) all indicate it to be multifunctional. We have performed molecular dynamics simulations of a conserved central segment of 18.5 kDa myelin basic protein (residues Glu80-Gly103, murine sequence numbering) in aqueous and membrane-associated environments to ascertain the stability of constituent secondary structure elements (α-helix from Glu80-Val91 and extended poly-proline type II from Thr92-Gly103) and the effects of phosphorylation of residues Thr92 and Thr95, individually and together. In aqueous solution, all four forms of the peptide bent in the middle to form a hydrophobic cluster. The phosphorylated variants were stabilized further by electrostatic interactions and formation of β-structures, in agreement with previous spectroscopic data. In simulations performed with the peptide in association with a dimyristoylphosphatidylcholine bilayer, the amphipathic α-helical segment remained stable and membrane-associated, although the degree of penetration was less in the phosphorylated variants, and the tilt of the α-helix with respect to the plane of the membrane also changed significantly with the modifications. The extended segment adjacent to this α-helix represents a putative SH3-ligand and remained exposed to the cytoplasm (and thus accessible to binding partners). The results of these simulations demonstrate how this segment of the protein can act as a molecular switch: an amphipathic α-helical segment of the protein is membrane-associated and presents a subsequent proline-rich segment to the cytoplasm for interaction with other proteins. Phosphorylation of threonyl residues alters the degree of membrane penetration of the α-helix and the accessibility of the proline-rich ligand and can stabilize a β-bend. A bend in this region of 18.5 kDa myelin basic protein suggests that the N- and C-termini of the proteins can interact with different leaflets of the myelin membrane and explain how a single protein can bring them close together.