Markus S. Miettinen

Concerted diffusion of lipids in raft-like membranes

Currently, there is no comprehensive model for the dynamics of cellular membranes. The understanding of even the basic dynamic processes, such as lateral diffusion of lipids, is still quite limited. Recent studies of one-component membrane systems have shown that instead of single-particle motions, the lateral diffusion is driven by a more complex, concerted mechanism for lipid diffusion (E. Falck et al., J. Am. Chem. Soc., 2008, 130, 44-45), where a lipid and its neighbors move in unison in terms of loosely defined clusters. In this work, we extend the previous study by considering the concerted lipid diffusion phenomena in many-component raft-like membranes. This nature of diffusion phenomena emerge in all the cases we have considered, including both atom-scale simulations of lateral diffusion within rafts and coarse-grained MARTINI simulations of diffusion in membranes characterized by coexistence of raft and non-raft domains. The data allows us to identify characteristic time scales for the concerted lipid motions, which turn out to range from hundreds of nanoseconds to several microseconds. Further, we characterize typical length scales associated with the correlated lipid diffusion patterns and find them to be about 10 nm, or even larger if weak correlations are taken into account. Finally, the concerted nature of lipid motions is also found in dissipative particle dynamics simulations of lipid membranes, clarifying the role of hydrodynamics (local momentum conservation) in membrane diffusion phenomena.

Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions

Phospholipids are essential building blocks of biological membranes. Despite a vast amount of very accurate experimental data, the atomistic resolution structures sampled by the glycerol backbone and choline headgroup in phoshatidylcholine bilayers are not known. Atomistic resolution molecular dynamics simulations have the potential to resolve the structures, and to give an arrestingly intuitive interpretation of the experimental data, but only if the simulations reproduce the data within experimental accuracy. In the present work, we simulated phosphatidylcholine (PC) lipid bilayers with 13 different atomistic models, and compared simulations with NMR experiments in terms of the highly structurally sensitive C–H bond vector order parameters. Focusing on the glycerol backbone and choline headgroups, we showed that the order parameter comparison can be used to judge the atomistic resolution structural accuracy of the models. Accurate models, in turn, allow molecular dynamics simulations to be used as an interpretation tool that translates these NMR data into a dynamic three-dimensional representation of biomolecules in biologically relevant conditions. In addition to lipid bilayers in fully hydrated conditions, we reviewed previous experimental data for dehydrated bilayers and cholesterol-containing bilayers, and interpreted them with simulations. Although none of the existing models reached experimental accuracy, by critically comparing them we were able to distill relevant chemical information: (1) increase of choline order parameters indicates the P–N vector tilting more parallel to the membrane, and (2) cholesterol induces only minor changes to the PC (glycerol backbone) structure. This work has been done as a fully open collaboration, using nmrlipids.blogspot.fi as a communication platform; all the scientific contributions were made publicly on this blog. During the open research process, the repository holding our simulation trajectories and files (https://zenodo.org/collection/user-nmrlipids) has become the most extensive publicly available collection of molecular dynamics simulation trajectories of lipid bilayers.

Ion dynamics in cationic lipid bilayer systems in saline solutions.

Positively charged lipid bilayer systems are a promising class of nonviral vectors for safe and efficient gene and drug delivery. Detailed understanding of these systems is therefore not only of fundamental but also of practical biomedical interest. Here, we study bilayers comprising a binary mixture of cationic dimyristoyltrimethylammoniumpropane (DMTAP) and zwitterionic (neutral) dimyristoylphosphatidylcholine (DMPC) lipids. Using atomistic molecular dynamics simulations, we address the effects of bilayer composition (cationic to zwitterionic lipid fraction) and of NaCl electrolyte concentration on the dynamical properties of these cationic lipid bilayer systems. We find that, despite the fact that DMPCs form complexes via Na(+) ions that bind to the lipid carbonyl oxygens, NaCl concentration has a rather minute effect on lipid diffusion. We also find the dynamics of Cl(-) and Na(+) ions at the water-membrane interface to differ qualitatively. Cl(-) ions have well-defined characteristic residence times of nanosecond scale. In contrast, the binding of Na(+) ions to the carbonyl region appears to lack a characteristic time scale, as the residence time distributions displayed power-law features. As to lateral dynamics, the diffusion of Na(+) ions within the water-membrane interface consists of two qualitatively different modes of motion: very slow diffusion when ions are bound to DMPC, punctuated by fast rapid jumps when detached from the lipids. Overall, the prolonged dynamics of the Na(+) ions are concluded to be interesting for the physics of the whole membrane, especially considering its interaction dynamics with charged macromolecular surfaces.

Assessing Polyglutamine Conformation in the Nucleating Event by Molecular Dynamics Simulations

Polyglutamine (polyQ) diseases comprise a group of dominantly inherited pathology caused by an expansion of an unstable polyQ stretch which is presumed to form β-sheets. Similar to other amyloid pathologies, polyQ amyloidogenesis occurs via a nucleated polymerization mechanism, and proceeds through energetically unfavorable nucleus whose existence and structure are difficult to detect. Here, we use atomistic molecular dynamics simulations in explicit solvent to assess the conformation of the polyQ stretch in the nucleus that initiates polyQ fibrillization. Comparison of the kinetic stability of various structures of polyQ peptide with a Q-length in the pathological range (Q40) revealed that steric zipper or nanotube-like structures (β-nanotube or β-pseudohelix) are not kinetically stable enough to serve as a template to initiate polyQ fibrillization as opposed to β-hairpin-based (β-sheet and β-sheetstack) or α-helical conformations. The selection of different structures of the polyQ stretch in the aggregation-initiating event may provide an alternative explanation for polyQ aggregate polymorphism.

Unfolding and Folding Internal Friction of β-Hairpins Is Smaller than That of α-Helices

By the forced unfolding of polyglutamine and polyalanine homopeptides in competing α-helix and β-hairpin secondary structures, we disentangle equilibrium free energetics from nonequilibrium dissipative effects. We find that α-helices are characterized by larger friction or dissipation upon unfolding, regardless of whether they are free energetically preferred over β-hairpins or not. Our analysis, based on MD simulations for atomistic peptide models with explicit water, suggests that this difference is related to the internal friction and mostly caused by the different number of intrapeptide hydrogen bonds in the α-helix and β-hairpin states.

Architecture of Polyglutamine-containing Fibrils from Time-resolved Fluorescence Decay

Background: Large polyglutamine-rich inclusions are hallmarks of CAG repeat pathologies. Results: Sensitive time-resolved fluorescent decay measurements combined with Monte Carlo simulations reveal the architecture of polyglutamine amyloid fibrils. Conclusion: Polyglutamine stretches are β-stranded in monomers and are organized into layered β-sheets with alternating N termini in amyloid fibrils. Significance: This sensitive approach can be routinely used to characterize the effect of new amyloid therapeutics. The disease risk and age of onset of Huntington disease (HD) and nine other repeat disorders strongly depend on the expansion of CAG repeats encoding consecutive polyglutamines (polyQ) in the corresponding disease protein. PolyQ length-dependent misfolding and aggregation are the hallmarks of CAG pathologies. Despite intense effort, the overall structure of these aggregates remains poorly understood. Here, we used sensitive time-dependent fluorescent decay measurements to assess the architecture of mature fibrils of huntingtin (Htt) exon 1 implicated in HD pathology. Varying the position of the fluorescent labels in the Htt monomer with expanded 51Q (Htt51Q) and using structural models of putative fibril structures, we generated distance distributions between donors and acceptors covering all possible distances between the monomers or monomer dimensions within the polyQ amyloid fibril. Using Monte Carlo simulations, we systematically scanned all possible monomer conformations that fit the experimentally measured decay times. Monomers with four-stranded 51Q stretches organized into five-layered β-sheets with alternating N termini of the monomers perpendicular to the fibril axis gave the best fit to our data. Alternatively, the core structure of the polyQ fibrils might also be a zipper layer with antiparallel four-stranded stretches as this structure showed the next best fit. All other remaining arrangements are clearly excluded by the data. Furthermore, the assessed dimensions of the polyQ stretch of each monomer provide structural evidence for the observed polyQ length threshold in HD pathology. Our approach can be used to validate the effect of pharmacological substances that inhibit or alter amyloid growth and structure.

Stable Polyglutamine Dimers Can Contain β-Hairpins with Interdigitated Side Chains But Not α-Helices, β-Nanotubes, β-Pseudohelices, or Steric Zippers

A common thread connecting nine fatal neurodegenerative protein aggregation diseases is an abnormally expanded polyglutamine tract found in the respective proteins. Although the structure of this tract in the large mature aggregates is increasingly well described, its structure in the small early aggregates remains largely unknown. As experimental evidence suggests that the most toxic species along the aggregation pathway are the small early ones, developing strategies to alleviate disease pathology calls for understanding the structure of polyglutamine peptides in the early stages of aggregation. Here, we present a criterion, grounded in available experimental data, that allows for using kinetic stability of dimers to assess whether a given polyglutamine conformer can be on the aggregation path. We then demonstrate that this criterion can be assessed using present-day molecular dynamics simulations. We find that although the α-helical conformer of polyglutamine is very stable, dimers of α-helices lack the kinetic stability necessary to support further oligomerization. Dimers of steric zipper, β-nanotube, and β-pseudohelix conformers are also too short-lived to initiate aggregation. The β-hairpin-containing conformers, instead, invariably form very stable dimers when their side chains are interdigitated. Combining these findings with the implications of recent solid-state NMR data on mature fibrils, we propose a possible pathway for the initial stages of polyglutamine aggregation, in which β-hairpin-containing conformers act as templates for fibril formation.

Coiled-Coils in Phage Display Screening: Insight into Exceptional Selectivity Provided by Molecular Dynamics

Involved in numerous key biological functions, protein helix-helix interactions follow a well-defined intermolecular recognition pattern. The characteristic structure of the α-helical coiled-coil allows for the specific randomization of clearly defined interaction partners within heteromeric systems. In this work, a rationally designed heterodimeric coiled-coil was used to investigate potential factors influencing the sequence selectivity in interhelical interactions.

A flexible approach to assess fluorescence decay functions in complex energy transfer systems.

BackgroundTime-correlated Förster resonance energy transfer (FRET) probes molecular distances with greater accuracy than intensity-based calculation of FRET efficiency and provides a powerful tool to study biomolecular structure and dynamics. Moreover, time-correlated photon count measurements bear additional information on the variety of donor surroundings allowing more detailed differentiation between distinct structural geometries which are typically inaccessible to general fitting solutions.ResultsHere we develop a new approach based on Monte Carlo simulations of time-correlated FRET events to estimate the time-correlated single photon counts (TCSPC) histograms in complex systems. This simulation solution assesses the full statistics of time-correlated photon counts and distance distributions of fluorescently labeled biomolecules. The simulations are consistent with the theoretical predictions of the dye behavior in FRET systems with defined dye distances and measurements of randomly distributed dye solutions. We validate the simulation results using a highly heterogeneous aggregation system and explore the conditions to use this tool in complex systems.ConclusionThis approach is powerful in distinguishing distance distributions in a wide variety of experimental setups, thus providing a versatile tool to accurately distinguish between different structural assemblies in highly complex systems.

Concluding the Amyloid Formation Pathway of a Coiled-Coil-Based Peptide from the Size of the Critical Nucleus

The size of the critical nucleus acting as intermediate in the amyloid formation of a model peptide is calculated. The theoretical approach is based on experimentally determined amyloid formation rates and gives new insights into the amyloid formation pathway.