Fuzzy supertertiary interactions within PSD-95 enable ligand binding

Experimental and computational investigation of the effect of Hsc70 structural variants on inhibiting amylin aggregation

The misfolding and aggregation of human islet amyloid polypeptide (hIAPP), also known as amylin, have been implicated in the pathogenesis of type 2 diabetes (T2D). Heat shock proteins, specifically, heat shock cognate 70 (Hsc70), are molecular chaperones that protect against hIAPP misfolding and inhibits its aggregation. Nevertheless, there is an incomplete understanding of the mechanistic interactions between Hsc70 domains and hIAPP, thus limiting their potential therapeutic role in diabetes. This study investigates the inhibitory capacities of different Hsc70 variants, aiming to identify the structural determinants that strike a balance between efficacy and cytotoxicity. Our experimental findings demonstrate that the ATPase activity of Hsc70 is not a pivotal factor for inhibiting hIAPP misfolding. We underscore the significance of the C-terminal substrate-binding domain of Hsc70 in inhibiting hIAPP aggregation, emphasizing that the removal of the lid subdomain diminishes the inhibitory effect of Hsc70. Additionally, we employed atomistic discrete molecular dynamics simulations to gain deeper insights into the interaction between Hsc70 variants and hIAPP. Integrating both experimental and computational findings, we propose a mechanism by which Hsc70's interaction with hIAPP monomers disrupts protein-protein connections, primarily by shielding the β-sheet edges of the Hsc70-β-sandwich. The distinctive conformational dynamics of the alpha helices of Hsc70 potentially enhance hIAPP binding by obstructing the exposed edges of the β-sandwich, particularly at the β5-β8 region along the alpha helix interface. This, in turn, inhibits fibril growth, and similar results were observed following hIAPP dimerization. Overall, this study elucidates the structural intricacies of Hsc70 crucial for impeding hIAPP aggregation, improving our understanding of the potential anti-aggregative properties of molecular chaperones in diabetes treatment.

Time-heterogeneity of the Förster Radius from Dipole Orientational Dynamics Explains Observed Dynamic Shift

Förster resonance energy transfer (FRET) is a quantum mechanical phenomenon involving the non-radiative transfer of energy between coupled electric dipoles. Due to the strong dependence of FRET on the distance between the dipoles, it is frequently used as a "molecular ruler" in biology, chemistry, and physics. This is done by placing dipolar molecules called dyes on molecules of interest. In time-resolved confocal single-molecule FRET (smFRET) experiments, the joint distribution of the FRET efficiency and the donor fluorescence lifetime can reveal underlying molecular conformational dynamics via deviation from their theoretical Förster relationship. This deviation is referred to as a dynamic shift. Quantifying the dynamic shift caused by the motion of the fluorescent dyes is essential to decoupling the dynamics of the studied molecules and the dyes. We develop novel Langevin models for the dye linker dynamics, including rotational dynamics, based on first physics principles and proper dye linker chemistry to match accessible volumes predicted by molecular dynamics simulations. By simulating the dyes' stochastic translational and rotational dynamics, we show that the observed dynamic shift can largely be attributed to the mutual orientational dynamics of the electric dipole moments associated with the dyes, not their accessible volume.

DNA controls the dimerization of the human FoxP1 forkhead domain

Transcription factors (TFs) regulate gene expression by binding to specific DNA sequences and gating access to genes. Even when the binding of TFs and their cofactors to DNA is reversible, indicating a reversible control of gene expression, there is little knowledge about the molecular effect DNA has on TFs. Using single-molecule multiparameter fluorescence spectroscopy, molecular dynamics simulations, and biochemical assays, we find that the monomeric form of the forkhead (FKH) domain of the human FoxP1 behaves as a disordered protein and increases its folded population when it dimerizes. Notably, DNA binding promotes a disordered FKH dimer bound to DNA, negatively controlling the stability of the dimeric FoxP1:DNA complex. The DNA-mediated reversible regulation on FKH dimers suggests that FoxP1-dependent gene suppression is unstable, and it must require the presence of other dimerization domains or cofactors to revert the negative impact exerted by the DNA.

Deciphering the influence of Y12L and N17H substitutions on the conformation and oligomerization of human calcitonin

The abnormal aggregation of human calcitonin (hCT) hormone peptides impairs their physiological function, leading to harmful immune responses and cytotoxicity, which limits their clinical utility. Interestingly, a representative hCT analog incorporating Y12L and N17H substitutions (DM-hCT) has shown reduced aggregation tendencies while maintaining bioactivity. But the molecular mechanism of Y12L and N17H substitutions on the conformational dynamics of hCT remains unclear. Here, we systematically investigated the folding and self-assembly dynamics of hCT and DM-hCT using atomistic discrete molecular dynamics (DMD) simulations. Our findings revealed that hCT monomers predominantly adopted unstructured conformations with dynamic helices. Oligomerization of hCT resulted in the formation of β-sheet-rich aggregates and β-barrel intermediates. The Y12L and N17H substitutions enhanced helical conformations and suppressed β-sheet formation in both monomers and oligomers. These substitutions stabilized the dynamic helices and disrupted aromatic interactions responsible for β-sheet formation at residue 12. Notably, DM-hCT assemblies still exhibited β-sheets in phenylalanine-rich and C-terminal hydrophobic regions, suggesting that future optimizations should focus on these areas. Our simulations provide insights into the molecular mechanisms underlying hCT aggregation and the amyloid-resistant effects of Y12L and N17H substitutions. These findings have valuable implications for the development of clinical hCT analogs.

Co-aggregation of α-synuclein with amyloid-β stabilizes β-sheet-rich oligomers and enhances the formation of β-barrels

Alzheimer's disease (AD) and Parkinson's disease (PD) are the two most common neurodegenerative diseases with markedly different pathological features of β-amyloid (Aβ) plaques and α-synuclein (αS) Lewy bodies (LBs), respectively. However, clinical overlaps in symptoms and pathologies between AD and PD are commonly observed caused by the cross-interaction between Aβ and αS. To uncover the molecular mechanisms behind their overlapping symptoms and pathologies, we computationally investigated the impact of αS on an Aβ monomer and dimerization using atomistic discrete molecular dynamics simulations (DMD). Our results revealed that αS could directly interact with Aβ monomers and dimers, thus forming β-sheet-rich oligomers, including potentially toxic β-barrel intermediates. The binding hotspot involved the second half of the N-terminal domain and NAC region in αS, along with residues 10-21 and 31-42 in Aβ. In their hetero-complex, the binding hotspot primarily assumed a β-sheet core buried inside, which was dynamically shielded by the highly charged, amyloid-resistant C-terminus of αS. Because the amyloid prion region was the same as the binding hotspot being buried, their fibrillization may be delayed, causing the toxic oligomers to increase. This study sheds light on the intricate relationship between Aβ and αS and provides insights into the overlapping pathology of AD and PD.

Unveiling Medin Folding and Dimerization Dynamics and Conformations via Atomistic Discrete Molecular Dynamics Simulations

Medin is a principal component of localized amyloid found in the vasculature of individuals over 50 years old. Its amyloid aggregation has been linked to endothelial dysfunction and vascular inflammation, contributing to the pathogenesis of various vascular diseases. Despite its significance, the structures of the medin monomer, oligomer, and fibril remain elusive, and the dynamic processes of medin aggregation are not fully understood. In this study, we comprehensively investigated the medin folding and dimerization dynamics and conformations using atomistic discrete molecular dynamics simulations. Our simulation results suggested that the folding initiation of the medin involved the formation of β-sheets around medin30-41 and medin42-50, with subsequent capping of other segments to their β-sheet edges. Medin monomers typically consisted of three or four β-strands, along with a dynamic N-terminal helix. Two isolated medin peptides readily aggregated into a β-sheet-rich dimer, displaying a strong aggregation propensity. Dimerization of medin not only enhanced the β-sheet conformations but also led to the formation of β-barrel oligomers. The aggregation tendencies of medin1-18 and medin19-29 were relatively weak. However, the segments of medin30-41 and medin42-50 played a crucial role as they primarily formed a β-sheet core and facilitated medin1-18 and medin19-29 to form intra- and interpeptide β-sheets. The findings highlight the critical role of the medin30-41 and medin42-50 regions in stabilizing the monomer structure and driving the medin amyloid aggregation. These regions could potentially serve as promising targets for designing antiamyloid inhibitors against amyloid aggregation of medin. Additionally, our study provides a full picture of the monomer conformations and dimerization dynamics for medin, which will help better understand the pathology of medin aggregation.

Direct Observation of Seeded Conformational Conversion of hIAPP In Silico Reveals the Mechanisms for Morphological Dependence and Asymmetry of Fibril Growth

Rapid growth of amyloid fibrils via a seeded conformational conversion of monomers is a critical step of fibrillization and important for disease transmission and progression. Amyloid fibrils often display diverse morphologies with distinct populations, and yet the molecular mechanisms of fibril elongation and their corresponding morphological dependence remain poorly understood. Here, we computationally investigated the single-molecular growth of two experimentally resolved human islet amyloid polypeptide fibrils of different morphologies. In both cases, the incorporation of monomers into preformed fibrils was observed. The conformational conversion dynamics was characterized by a small number of fibril growth intermediates. Fibril morphology affected monomer binding at fibril elongation and lateral surfaces as well as the seeded conformational conversion dynamics at the fibril ends, resulting in different fibril elongation rates and populations. We also observed an asymmetric fibril growth as in our prior experiments, attributing to differences of two fibril ends in terms of their local surface curvatures and exposed hydrogen-bond donors and acceptors. Together, our mechanistic findings afforded a theoretical basis for delineating different amyloid strains-entailed divergent disease progression.

Spike Protein Fragments Promote Alzheimer's Amyloidogenesis

Alzheimer's disease (AD) is a major cause of dementia inducing memory loss, cognitive decline, and mortality among the aging population. While the amyloid aggregation of peptide Aβ has long been implicated in neurodegeneration in AD, primarily through the production of toxic polymorphic aggregates and reactive oxygen species, viral infection has a less explicit role in the etiology of the brain disease. On the other hand, while the COVID-19 pandemic is known to harm human organs and function, its adverse effects on AD pathobiology and other human conditions remain unclear. Here we first identified the amyloidogenic potential of 1058HGVVFLHVTYV1068, a short fragment of the spike protein of SARS-CoV-2 coronavirus. The peptide fragment was found to be toxic and displayed a high binding propensity for the amyloidogenic segments of Aβ, thereby promoting the aggregation and toxicity of the peptide in vitro and in silico, while retarding the hatching and survival of zebrafish embryos upon exposure. Our study implicated SARS-CoV-2 viral infection as a potential contributor to AD pathogenesis, a little explored area in our quest for understanding and overcoming Long Covid.

Hybrid computational methods combining experimental information with molecular dynamics

A goal of structural biology is to understand how macromolecules carry out their biological roles by identifying their metastable states, mechanisms of action, pathways leading to conformational changes, and the thermodynamic and kinetic relationships between those states. Integrative modeling brings structural insights into systems where traditional structure determination approaches cannot help. We focus on the synergies and challenges of integrative modeling combining experimental data with molecular dynamics simulations.

Domain tethering impacts dimerization and DNA-mediated allostery in the human transcription factor FoxP1

Transcription factors are multidomain proteins with specific DNA binding and regulatory domains. In the human FoxP subfamily (FoxP1, FoxP2, FoxP3, and FoxP4) of transcription factors, a 90 residue-long disordered region links a Leucine Zipper (ZIP)-known to form coiled-coil dimers-and a Forkhead (FKH) domain-known to form domain swapping dimers. We used replica exchange discrete molecular dynamics simulations, single-molecule fluorescence experiments, and other biophysical tools to understand how domain tethering in FoxP1 impacts dimerization at ZIP and FKH domains and how DNA binding allosterically regulates their dimerization. We found that domain tethering promotes FoxP1 dimerization but inhibits a FKH domain-swapped structure. Furthermore, our findings indicate that the linker mediates the mutual organization and dynamics of ZIP and FKH domains, forming closed and open states with and without interdomain contacts, thus highlighting the role of the linkers in multidomain proteins. Finally, we found that DNA allosterically promotes structural changes that decrease the dimerization propensity of FoxP1. We postulate that, upon DNA binding, the interdomain linker plays a crucial role in the gene regulatory function of FoxP1.

Coevolution and smFRET Enhances Conformation Sampling and FRET Experimental Design in Tandem PDZ1-2 Proteins

The structural flexibility of proteins is crucial for their functions. Many experimental and computational approaches can probe protein dynamics across a range of time and length-scales. Integrative approaches synthesize the complementary outputs of these techniques and provide a comprehensive view of the dynamic conformational space of proteins, including the functionally relevant limiting conformational states and transition pathways between them. Here, we introduce an integrative paradigm to model the conformational states of multidomain proteins. As a model system, we use the first two tandem PDZ domains of postsynaptic density protein 95. First, we utilize available sequence information collected from genomic databases to identify potential amino acid interactions in the PDZ1-2 tandem that underlie modeling of the functionally relevant conformations maintained through evolution. This was accomplished through combination of coarse-grained structural modeling with outputs from direct coupling analysis measuring amino acid coevolution, a hybrid approach called SBM+DCA. We recapitulated five distinct, experimentally derived PDZ1-2 tandem conformations. In addition, SBM+DCA unveiled an unidentified, twisted conformation of the PDZ1-2 tandem. Finally, we implemented an integrative framework for the design of single-molecule Förster resonance energy transfer (smFRET) experiments incorporating the outputs of SBM+DCA with simulated FRET observables. This resulting FRET network is designed to mutually resolve the predicted limiting state conformations through global analysis. Using simulated FRET observables, we demonstrate that structural modeling with the newly designed FRET network is expected to outperform a previously used empirical FRET network at resolving all states simultaneously. Integrative approaches to experimental design have the potential to provide a new level of detail in characterizing the evolutionarily conserved conformational landscapes of proteins, and thus new insights into functional relevance of protein dynamics in biological function.

Structural Perturbation of Monomers Determines the Amyloid Aggregation Propensity of Calcitonin Variants

Human calcitonin (hCT) is a polypeptide hormone that participates in calcium-phosphorus metabolism. Irreversible aggregation of 32-amino acid hCT into β-sheet-rich amyloid fibrils impairs physiological activity and increases the risk of medullary carcinoma of the thyroid. Amyloid-resistant hCT derivatives substituting critical amyloidogenic residues are of particular interest for clinical applications as therapeutic drugs against bone-related diseases. Uncovering the aggregation mechanism of hCT at the molecular level, therefore, is important for the design of amyloid-resistant hCT analogues. Here, we investigated the aggregation dynamics of hCT, non-amyloidogenic salmon calcitonin (sCT), and two hCT analogues with reduced aggregation tendency─TL-hCT and phCT─using long timescale discrete molecular dynamics simulations. Our results showed that hCT monomers mainly adopted unstructured conformations with dynamically formed helices around the central region. hCT self-assembled into helix-rich oligomers first, followed by a conformational conversion into β-sheet-rich oligomers with β-sheets formed by residues 10-30 and stabilized by aromatic and hydrophobic interactions. Our simulations confirmed that TL-hCT and phCT oligomers featured more helices and fewer β-sheets than hCT. Substitution of central aromatic residues with leucine in TL-hCT and replacing C-terminal hydrophobic residue with hydrophilic amino acid in phCT only locally suppressed β-sheet propensities in the central region and C-terminus, respectively. Having mutations in both central and C-terminal regions, sCT monomers and dynamically formed oligomers predominantly adopted helices, confirming that both central aromatic and C-terminal hydrophobic residues played important roles in the fibrillization of hCT. We also observed the formation of β-barrel intermediates, postulated as the toxic oligomers in amyloidosis, for hCT but not for sCT. Our computational study depicts a complete picture of the aggregation dynamics of hCT and the effects of mutations. The design of next-generation amyloid-resistant hCT analogues should consider the impact on both amyloidogenic regions and also take into account the amplification of transient β-sheet population in monomers upon aggregation.

Illuminating Intrinsically Disordered Proteins with Integrative Structural Biology

Intense study of intrinsically disordered proteins (IDPs) did not begin in earnest until the late 1990s when a few groups, working independently, convinced the community that these 'weird' proteins could have important functions. Over the past two decades, it has become clear that IDPs play critical roles in a multitude of biological phenomena with prominent examples including coordination in signaling hubs, enabling gene regulation, and regulating ion channels, just to name a few. One contributing factor that delayed appreciation of IDP functional significance is the experimental difficulty in characterizing their dynamic conformations. The combined application of multiple methods, termed integrative structural biology, has emerged as an essential approach to understanding IDP phenomena. Here, we review some of the recent applications of the integrative structural biology philosophy to study IDPs.

Different Forms of Disorder in NMDA-Sensitive Glutamate Receptor Cytoplasmic Domains Are Associated with Differences in Condensate Formation

The N-methyl-D-aspartate (NMDA)-sensitive glutamate receptor (NMDAR) helps assemble downstream signaling pathways through protein interactions within the postsynaptic density (PSD), which are mediated by its intracellular C-terminal domain (CTD). The most abundant NMDAR subunits in the brain are GluN2A and GluN2B, which are associated with a developmental switch in NMDAR composition. Previously, we used single molecule fluorescence resonance energy transfer (smFRET) to show that the GluN2B CTD contained an intrinsically disordered region with slow, hop-like conformational dynamics. The CTD from GluN2B also undergoes liquid-liquid phase separation (LLPS) with synaptic proteins. Here, we extend these observations to the GluN2A CTD. Sequence analysis showed that both subunits contain a form of intrinsic disorder classified as weak polyampholytes. However, only GluN2B contained matched patterning of arginine and aromatic residues, which are linked to LLPS. To examine the conformational distribution, we used discrete molecular dynamics (DMD), which revealed that GluN2A favors extended disordered states containing secondary structures while GluN2B favors disordered globular states. In contrast to GluN2B, smFRET measurements found that GluN2A lacked slow conformational dynamics. Thus, simulation and experiments found differences in the form of disorder. To understand how this affects protein interactions, we compared the ability of these two NMDAR isoforms to undergo LLPS. We found that GluN2B readily formed condensates with PSD-95 and SynGAP, while GluN2A failed to support LLPS and instead showed a propensity for colloidal aggregation. That GluN2A fails to support this same condensate formation suggests a developmental switch in LLPS propensity.

Unmasking a two-faced protein

Single-molecule fluorescence spectroscopy and molecular dynamics simulations illuminate the structure and dynamics of PSD-95, a protein involved in neural plasticity.