The beta-strand-loop-beta-strand conformation is marginally populated in beta2-microglobulin (20-41) peptide in solution as revealed by replica exchange molecular dynamics simulations

Single-molecule visualization determines conformational substate ensembles in β-sheet-rich peptide fibrils

An understanding of protein conformational ensembles is essential for revealing the underlying mechanisms of interpeptide recognition and association. However, experimentally resolving multiple simultaneously existing conformational substates remains challenging. Here, we report the use of scanning tunneling microscopy (STM) to analyze the conformational substate ensembles of β sheet peptides with a submolecular resolution (in-plane <2.6 Å). We observed ensembles of more than 10 conformational substates (with free energy fluctuations between several k_BTs) in peptide homoassemblies of keratin (KRT) and amyloidal peptides (-5Aβ42 and TDP-43 341-357). Furthermore, STM reveals a change in the conformational ensemble of peptide mutants, which is correlated with the macroscopic properties of peptide assemblies. Our results demonstrate that the STM-based single-molecule imaging can capture a thorough picture of the conformational substates with which to build an energetic landscape of interconformational interactions and can rapidly screen conformational ensembles, which can complement conventional characterization techniques.

The Early Phase of β2-Microglobulin Aggregation: Perspectives From Molecular Simulations

Protein β2-microglobulin is the causing agent of two amyloidosis, dialysis related amyloidosis (DRA), affecting the bones and cartilages of individuals with chronic renal failure undergoing long-term hemodialysis, and a systemic amyloidosis, found in one French family, which impairs visceral organs. The protein’s small size and its biomedical significance attracted the attention of theoretical scientists, and there are now several studies addressing its aggregation mechanism in the context of molecular simulations. Here, we review the early phase of β2-microglobulin aggregation, by focusing on the identification and structural characterization of monomers with the ability to trigger aggregation, and initial small oligomers (dimers, tetramers, hexamers etc.) formed in the so-called nucleation phase. We focus our analysis on results from molecular simulations and integrate our views with those coming from in vitro experiments to provide a broader perspective of this interesting field of research. We also outline directions for future computer simulation studies.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Molecular dynamics simulation of β₂-microglobulin in denaturing and stabilizing conditions

β₂-Microglobulin has been a model system for the study of fibril formation for 20 years. The experimental study of β₂-microglobulin structure, dynamics, and thermodynamics in solution, at atomic detail, along the pathway leading to fibril formation is difficult because the onset of disorder and aggregation prevents signal resolution in Nuclear Magnetic Resonance experiments. Moreover, it is difficult to characterize conformers in exchange equilibrium. To gain insight (at atomic level) on processes for which experimental information is available at molecular or supramolecular level, molecular dynamics simulations have been widely used in the last decade. Here, we use molecular dynamics to address three key aspects of β₂-microglobulin, which are known to be relevant to amyloid formation: (1) 60 ns molecular dynamics simulations of β₂-microglobulin in trifluoroethanol and in conditions mimicking low pH are used to study the behavior of the protein in environmental conditions that are able to trigger amyloid formation; (2) adaptive biasing force molecular dynamics simulation is used to force cis-trans isomerization at Proline 32 and to calculate the relative free energy in the folded and unfolded state. The native-like trans-conformer (known as intermediate 2 and determining the slow phase of refolding), is simulated for 10 ns, detailing the possible link between cis-trans isomerization and conformational disorder; (3) molecular dynamics simulation of highly concentrated doxycycline (a molecule able to suppress fibril formation) in the presence of β₂-microglobulin provides details of the binding modes of the drug and a rationale for its effect.

K3 fragment of amyloidogenic beta(2)-microglobulin forms ion channels: implication for dialysis related amyloidosis

Beta(2)-microglobulin (beta(2)m) amyloid deposits are linked to dialysis-related amyloidosis (DRA) in hemodialysis patients. The mechanism by which beta(2)m causes DRA is not understood. It is also unclear whether only the full-length beta(2)m induces pathophysiology or if proteolytic fragments are sufficient for inducing this effect. Ser20-Lys41 (K3) is a digestion fragment of full-length beta(2)m. Solid state NMR (ssNMR) combined with X-ray diffraction and atomic force microscopy (AFM) revealed the characteristic oligomeric amyloid conformation of the U-turn beta-strand-turn-beta-strand motif stacked in parallel and stabilized by intermolecular interactions also shown by Abeta(9-40)/Abeta(17-42) and the CA150 WW domain. Here we use the K3 U-turn atomic coordinates and molecular dynamic (MD) simulations to model K3 channels in the membrane. Consistent with previous AFM imaging of other amyloids that show channel-like structures in the membrane, in the simulations K3 also forms ion channels with 3-6 loosely attached mobile subunits. We carry out AFM, single channel electrical recording, and fluorescence imaging experiments. AFM images display 3D ion channel topography with shapes, morphologies, and dimensions consistent with the theoretical model. Electrical conductance measurements indicate multiple single channel conductances, suggesting that various K3 oligomer sizes can constitute the channel structure. Fluorescence measurements in kidney cells show channel-mediated cell calcium uptake. These results suggest that the beta(2)m-induced DRA can be mediated by ion channels formed by its K3 fragment. Because the beta-strand-turn-beta-strand motif appears to be a universal amyloid feature, its ability to form ion channels further suggests that the motif may play a generic role in toxicity.

Structures and thermodynamics of Alzheimer's amyloid-beta Abeta(16-35) monomer and dimer by replica exchange molecular dynamics simulations: implication for full-length Abeta fibrillation

Many proteins display a strand-loop-strand motif in their amyloid fibrillar states. For instance, the amyloid beta-protein, Abeta1-40, associated with Alzheimer's disease, displays a loop at positions 22-28 in its amyloid fibril state. It has been suggested that this loop could appear early in the aggregation process, but quantitative information regarding its presence in small oligomers remains scant. Because residues 1-15 are disordered in Abeta1-42 fibrils and Abeta10-35 forms fibrils in vitro, we select the peptide Abeta16-35, centered on residues 22-28 and determine the structures and thermodynamics of the monomer and dimer using coarse-grained implicit solvent replica exchange molecular dynamics simulations. Our simulations totalling 5 mus for the monomer and 12 micros for the dimer show no sign of strong secondary structure signals in both instances and the significant impact of dimerization on the global structure of Abeta16-35. They reveal however that the loop 22-28 acts as a quasi-independent unit in both species. The loop structure ensemble we report in Abeta16-35 monomer and dimer has high similarity to the loop formed by the Abeta21-30 peptide in solution and, to a lesser extent, to the loop found in Abeta1-40 fibrils. We discuss the implications of our findings on the assembly of full-length Abeta.