Amyloid ß Fibril Elongation by Monomers Involves Disorder at the Tip.

The growth of amyloid fibrils from Aß1-42 peptide, one of the key pathogenic players in Alzheimer's disease, is believed to follow a nucleation-elongation mechanism. Fibril elongation is often described as a "dock-lock" procedure, where a disordered monomer adsorbs to an existing fibril in a relatively fast process (docking), followed by a slower conformational transition toward the ordered state of the template (locking). Here, we use molecular dynamics simulations of an ordered pentamer of Aß42 at fully atomistic resolution, which includes solvent, to characterize the elongation process. We construct a Markov state model from an ensemble of short trajectories generated by an advanced sampling algorithm that efficiently diversifies a subset of the system without any bias forces. This subset corresponds to selected dihedral angles of the peptide chain at the fibril tip favored to be the fast growing one experimentally. From the network model, we extract distinct locking pathways covering time scales in the high microsecond regime. Slow steps are associated with the exchange of hydrophobic contacts, between nonnative and native intermolecular contacts as well as between intra- and intermolecular ones. The N-terminal segments, which are disordered in fibrils and typically considered inert, are able to shield the lateral interfaces of the pentamer. We conclude by discussing our findings in the context of a refined dock-lock model of Aß fibril elongation, which involves structural disorder for more than one monomer at the growing tip.

A combined NMR and computational approach to investigate peptide binding to a designed Armadillo repeat protein.

The specific recognition of peptide sequences by proteins plays an important role both in biology and in diagnostic applications. Here we characterize the relatively weak binding of the peptide neurotensin (NT) to the previously developed Armadillo repeat protein VG_328 by a multidisciplinary approach based on solution NMR spectroscopy, mutational studies, and molecular dynamics (MD) simulations, totaling 20µs for all MD runs. We describe assignment challenges arising from the repetitive nature of the protein sequence, and we present novel approaches to address them. Partial assignments obtained for VG_328 in combination with chemical shift perturbations allowed us to identify the repeats not involved in binding. Their subsequent elimination resulted in a reduced-size binder with very similar affinity for NT, for which near-complete backbone assignments were achieved. A binding mode suggested by automatic docking and further validated by explicit solvent MD simulations is consistent with paramagnetic relaxation enhancement data collected using spin-labeled NT. Favorable intermolecular interactions are observed in the MD simulations for the residues that were previously shown to contribute to binding in an Ala scan of NT. We further characterized the role of residues within the N-cap for protein stability and peptide binding. Our multidisciplinary approach demonstrates that an initial low-resolution picture for a low-micromolar-peptide binder can be refined through the combination of NMR, protein design, docking, and MD simulations to establish its binding mode, even in the absence of crystallographic data, thereby providing valuable information for further design.

Spontaneous self-assembly of engineered armadillo repeat protein fragments into a folded structure.

Repeat proteins are built of modules, each of which constitutes a structural motif. We have investigated whether fragments of a designed consensus armadillo repeat protein (ArmRP) recognize each other. We examined a split ArmRP consisting of an N-capping repeat (denoted Y), three internal repeats (M), and a C-capping repeat (A). We demonstrate that the C-terminal MA fragment adopts a fold similar to the corresponding part of the entire protein. In contrast, the N-terminal YM2 fragment constitutes a molten globule. The two fragments form a 1:1 YM2:MA complex with a nanomolar dissociation constant essentially identical to the crystal structure of the continuous YM3A protein. Molecular dynamics simulations show that the complex is structurally stable over a 1 µs timescale and reveal the importance of hydrophobic contacts across the interface. We propose that the existence of a stable complex recapitulates possible intermediates in the early evolution of these repeat proteins.

Charge distribution and imperfect amphipathicity affect pore formation by antimicrobial peptides.

Antimicrobial peptides often permeabilize biological membranes via a pore mechanism. Two pore types have been proposed: toroidal, where the pore is partly lined by lipid, and barrel-stave, where a cylindrical pore is completely lined by peptides. What drives the preference of antimicrobial peptides for a certain pore type is not yet fully understood. According to neutron scattering and oriented circular dichroism, melittin and MG-H2 induce toroidal pores whereas alamethicin forms barrel-stave pores. In previous work we found that indeed melittin seems to favor toroidal pores whereas alamethicin favors cylindrical pores. Here we designed mutants of these two peptides and the magainin analog MG-H2, aimed to probe how the distribution of charges along the helix and its imperfectly amphipathic structure influence pore formation. Molecular dynamics (MD) simulations of the peptides in a pre-formed cylindrical pore have been performed. The duration of the simulations was 136ns to 216ns. We found that a melittin mutant with lysine 7 neutralized favors cylindrical pores whereas a MG-H2 mutant with lysines in the N-terminal half of these peptides neutralized and an alamethicin mutant with a positive charge at the position 7 form semitoroidal pores. These results suggest that charged residues within the N-terminal half are important for toroidal pore formation. Toroidal pores produced by MG-H2 are more disordered than the melittin pores, likely because of the charged residues located in the middle of the MG-H2 helix (K11 and K14). Imperfect amphipathicity of melittin seems to play a role in its preference for toroidal pores since the substitutions of charged residues located within the nonpolar face by hydrophobic residues suppress evolution of a toroidal pore. The mutations change the position of lysine 7 near the N-terminus, relative to the lower leaflet headgroups. The MD simulations also show that the melittin P14A mutant forms a toroidal pore, but its configuration diverges from that of melittin and it is probably metastable.

Slow dissociation of a charged ligand: analysis of the primary quinone Q(A) site of photosynthetic bacterial reaction centers.

Reaction centers (RCs) are integral membrane proteins that undergo a series of electron transfer reactions during the process of photosynthesis. In the Q(A) site of RCs from Rhodobacter sphaeroides, ubiquinone-10 is reduced, by a single electron transfer, to its semiquinone. The neutral quinone and anionic semiquinone have similar affinities, which is required for correct in situ reaction thermodynamics. A previous study showed that despite similar affinities, anionic quinones associate and dissociate from the Q(A) site at rates ≈10(4) times slower than neutral quinones indicating that anionic quinones encounter larger binding barriers (Madeo, J.; Gunner, M. R. Modeling binding kinetics at the Q(A) site in bacterial reaction centers. Biochemistry 2005, 44, 10994-11004). The present study investigates these barriers computationally, using steered molecular dynamics (SMD) to model the unbinding of neutral ground state ubiquinone (UQ) and its reduced anionic semiquinone (SQ(-)) from the Q(A) site. In agreement with experiment, the SMD unbinding barrier for SQ(-) is larger than for UQ. Multi Conformational Continuum Electrostatics (MCCE), used here to calculate the binding energy, shows that SQ(-) and UQ have comparable affinities. In the Q(A) site, there are stronger binding interactions for SQ(-) compared to UQ, especially electrostatic attraction to a bound non-heme Fe(2+). These interactions compensate for the higher SQ(-) desolvation penalty, allowing both redox states to have similar affinities. These additional interactions also increase the dissociation barrier for SQ(-) relative to UQ. Thus, the slower SQ(-) dissociation rate is a direct physical consequence of the additional binding interactions required to achieve a Q(A) site affinity similar to that of UQ. By a similar mechanism, the slower association rate is caused by stronger interactions between SQ(-) and the polar solvent. Thus, stronger interactions for both the unbound and bound states of charged and highly polar ligands can slow their binding kinetics without a conformational gate. Implications of this for other systems are discussed.

Antimicrobial peptides in toroidal and cylindrical pores.

Antimicrobial peptides (AMPs) are small, usually cationic peptides, which permeabilize biological membranes. Their mechanism of action is still not well understood. Here we investigate the preference of alamethicin and melittin for pores of different shapes, using molecular dynamics (MD) simulations of the peptides in pre-formed toroidal and cylindrical pores. When an alamethicin hexamer is initially embedded in a cylindrical pore, at the end of the simulation the pore remains cylindrical or closes if glutamines in the N-termini are not located within the pore. On the other hand, when a melittin tetramer is embedded in toroidal pore or in a cylindrical pore, at the end of the simulation the pore is lined both with peptides and lipid headgroups, and, thus, can be classified as a toroidal pore. These observations agree with the prevailing views that alamethicin forms barrel-stave pores whereas melittin forms toroidal pores. Both alamethicin and melittin form amphiphilic helices in the presence of membranes, but their net charge differs; at pH approximately 7, the net charge of alamethicin is -1 whereas that of melittin is +5. This gives rise to stronger electrostatic interactions of melittin with membranes than those of alamethicin. The melittin tetramer interacts more strongly with lipids in the toroidal pore than in the cylindrical one, due to more favorable electrostatic interactions.

Antimicrobial peptides bind more strongly to membrane pores.

Antimicrobial peptides (AMPs) are small, usually cationic peptides, which permeabilize bacterial membranes. Understanding their mechanism of action might help design better antibiotics. Using an implicit membrane model, modified to include pores of different shapes, we show that four AMPs (alamethicin, melittin, a magainin analogue, MG-H2, and piscidin 1) bind more strongly to membrane pores, consistent with the idea that they stabilize them. The effective energy of alamethicin in cylindrical pores is similar to that in toroidal pores, whereas the effective energy of the other three peptides is lower in toroidal pores. Only alamethicin intercalates into the membrane core; MG-H2, melittin and piscidin are located exclusively at the hydrophobic/hydrophilic interface. In toroidal pores, the latter three peptides often bind at the edge of the pore, and are in an oblique orientation. The calculated binding energies of the peptides are correlated with their hemolytic activities. We hypothesize that one distinguishing feature of AMPs may be the fact that they are imperfectly amphipathic which allows them to bind more strongly to toroidal pores. An initial test on a melittin-based mutant seems to support this hypothesis.

Computational prediction of the binding site of proteinase 3 to the plasma membrane.

Proteinase 3 (PR3) is a neutrophil-derived serine proteinase localized within cytoplasmic granules which can be released upon activation. PR3 is exposed at the neutrophil plasma membrane where it can mediate proinflammatory effects. Moreover, PR3 membrane expression is of special relevance in patients with Wegener's granulomatosis, a systemic vasculitis presenting anticytoplasmic neutrophil autoantibodies (ANCA) against PR3, which can bind to PR3 expressed at the surface of neutrophils and amplify their activation state. Therefore, it is of special relevance to unravel the molecular mechanisms governing its association with the membrane to be able to modulate it. To this end, we performed molecular dynamics (MD) simulations of PR3 with the implicit membrane model IMM1-GC to identify its interfacial binding site (IBS). Both the energies and structures resulting from the MD suggest that PR3 associates strongly with anionic membranes. We observe a unique IBS consisting of five basic (R177, R186A, R186B, K187, R222) and six hydrophobic (F165, F166, F224, L223, F184, W218) amino acids. The basic residues provide the driving force to orient PR3 at the membrane surface, so that the hydrophobic residues can anchor into the hydrocarbon region. Energy decomposition and in silico mutations show that only a few residues account for the membrane association. Similar calculations with HNE suggest a different membrane-binding mechanism. Our results agree with previous experimental observations and this work predicts, for the first time, the structural determinants of the binding of PR3 to membranes.

Membrane-bound structure and energetics of alpha-synuclein.

Aggregation and fibrillation of alpha-synuclein bound to membranes are believed to be involved in Parkinson's and other neurodegenerative diseases. On SDS micelles, the N-terminus of alpha-synuclein forms two curved helices linked by a short loop. However, its structure on lipid bilayers has not been experimentally resolved. Using MD simulations with an implicit membrane model we show here that, on a planar mixed membrane, the truncated alpha-synuclein (residues 1-95) forms a bent helix. Bending of the helix is not due to the protein sequence or membrane binding, but to collective motions of the long helix. The backbone of the helix is approximately 2.5 A above the membrane surface, with some residues partially inserted in the membrane core. The helix periodicity is 11/3 (11 residues complete three full turns) as opposed to 18/5 periodicity of an ideal alpha-helix, with hydrophobic residues towards the membrane, negatively charged residues towards the solvent and lysines on the polar/nonpolar interface. A series of threonines, which are characteristic for alpha-synuclein and perhaps a phosphorylation site, is also located at the hydrophobic/hydrophilic interface with their side chain often hydrogen bonded to the main-chain atom. The calculations show that the energy penalty for change in periodicity from the 18/5 to 11/3 on the anionic membrane is overcome by favorable solvation energy. The binding of truncated alpha-synuclein to membranes is weak. It prefers anionic membranes but it also binds marginally to a neutral membrane, via its C-terminus. Dimerization of helical monomers on the mixed membrane is energetically favorable. However, it slightly interferes with membrane binding. This might promote lateral diffusion of the protein on the membrane surface and facilitate assembly of oligomers that precede fibrillation.

Modeling fatty acid delivery from intestinal fatty acid binding protein to a membrane.

Intestinal fatty acid binding protein (IFABP) interacts with biological membranes and delivers fatty acid (FA) into them via a collisional mechanism. However, the membrane-bound structure of the protein and the pathway of FA transfer are not precisely known. We used molecular dynamics (MD) simulations with an implicit membrane model to determine the optimal orientation of apo- and holo-IFABP (bound with palmitate) on an anionic membrane. In this orientation, the helical portal region, delimited by the alphaII helix and the betaC-betaD and betaE-betaF turns, is oriented toward the membrane whereas the putative beta-strand portal, delimited by the betaB-betaC, betaF-betaG, betaH-betaI turns and the N terminus, is exposed to solvent. Starting from the MD structure of holo-IFABP in the optimal orientation relative to the membrane, we examined the release of palmitate via both pathways. Although the domains can widen enough to allow the passage of palmitate, fatty acid release through the helical portal region incurs smaller conformational changes and a lower energetic cost.

Calculations of pH-dependent binding of proteins to biological membranes.

Binding of proteins to membranes is often accompanied by titration of ionizable residues and is, therefore, dependent on pH. We present a theoretical treatment and computational approach for predicting absolute, pH-dependent membrane binding free energies. The standard free energy of binding, DeltaG, is defined as -RTln(P(b)/P(f)), where P(b) and P(f) are the amounts of bound and free protein. The apparent pK(a) of binding is the pH value at which P(b) and P(f) are equal. Proteins bind to the membrane in the pH range where DeltaG is negative. The components of the binding free energy are (a) the free energy cost of ionization state changes (DeltaG(ion)), (b) the effective energy of transfer from solvent to the membrane surface, (c) the translational/rotational entropy cost of binding, and (d) an ideal entropy term that depends on the relative volume of the bound and free state and therefore depends on lipid concentration. Calculation of the first term requires determination of pK(a) values in solvent and on the membrane surface. All energies required by the method are obtained from molecular dynamics trajectories on an implicit membrane (IMM1-GC). The method is tested on pentalysine and the helical peptide VEEKS, derived from the membrane-binding domain of phosphocholine cytidylyltransferase. The agreement between the measured and the calculated free energies of binding of pentalysine is good. The extent of membrane binding of VEEKS is, however, underestimated compared to experiment. Calculations of the interaction energy between two VEEKS helices on the membrane suggest that the discrepancy is mainly due to the neglect of protein-protein interactions on the membrane surface.

Interfacial slip in sheared polymer blends.

We have developed a dynamic self-consistent field theory, without any adjustable parameters, for unentangled polymer blends under shear. Our model accounts for the interaction between polymers, and enables one to compute the evolution of the local rheology, microstructure, and the conformations of the polymer chains under shear self-consistently. We use this model to study the interfacial dynamics in sheared polymer blends and make a quantitative comparison between this model and molecular dynamics simulations. We find good agreement between the two methods.

Dynamic self-consistent field theory for unentangled homopolymer fluids.

We present a lattice formulation of a dynamic self-consistent field (DSCF) theory that is capable of resolving interfacial structure, dynamics, and rheology in inhomogeneous, compressible melts and blends of unentangled homopolymer chains. The joint probability distribution of all the Kuhn segments in the fluid, interacting with adjacent segments and walls, is approximated by a product of one-body probabilities for free segments interacting solely with an external potential field that is determined self-consistently. The effect of flow on ideal chain conformations is modeled with finitely extensible, nonlinearly elastic dumbbells in the Peterlin approximation, and related to stepping probabilities in a random walk. Free segment and stepping probabilities generate statistical weights for chain conformations in a self-consistent field, and determine local volume fractions of chain segments. Flux balance across unit lattice cells yields mean field transport equations for the evolution of free segment probabilities and of momentum densities on the Kuhn length scale. Diffusive and viscous contributions to the fluxes arise from segmental hops modeled as a Markov process, with transition rates reflecting changes in segmental interaction, kinetic energy, and entropic contributions to the free energy under flow. We apply the DSCF equations to study both transient and steady-state interfacial structure, flow, and rheology in a sheared planar channel containing either a one-component melt or a phase-separated, two-component blend.