Binding Dynamics of a Stapled Peptide Targeting the Transcription Factor NF-Y

Transcription factors (TFs) play a central role in gene regulation, and their malfunction can result in a plethora of severe diseases. TFs are therefore interesting therapeutic targets, but their involvement in protein-protein interaction networks and the frequent lack of well-defined binding pockets render them challenging targets for classical small molecules. As an alternative, peptide-based scaffolds have proven useful, in particular with an α-helical active conformation. Peptide-based strategies often require extensive structural optimization efforts, which could benefit from a more detailed understanding of the dynamics in inhibitor/protein interactions. In this study, we investigate how truncated stapled α-helical peptides interact with the transcription factor Nuclear Factor-Y (NF-Y). We identified a 13-mer minimal binding core region, for which two crystal structures with an altered C-terminal peptide conformation when bound to NF-Y were obtained. Subsequent molecular dynamics simulations confirmed that the C-terminal part of the stapled peptide is indeed relatively flexible while still showing defined interactions with NF-Y. Our findings highlight the importance of flexibility in the bound state of peptides, which can contribute to overall binding affinity.

Elucidating the role of water in collagen self-assembly by isotopically modulating collagen hydration

Water is known to play an important role in collagen self-assembly, but it is still largely unclear how water-collagen interactions influence the assembly process and determine the fibril network properties. Here, we use the H[Formula: see text]O/D[Formula: see text]O isotope effect on the hydrogen-bond strength in water to investigate the role of hydration in collagen self-assembly. We dissolve collagen in H[Formula: see text]O and D[Formula: see text]O and compare the growth kinetics and the structure of the collagen assemblies formed in these water isotopomers. Surprisingly, collagen assembly occurs ten times faster in D[Formula: see text]O than in H[Formula: see text]O, and collagen in D[Formula: see text]O self-assembles into much thinner fibrils, that form a more inhomogeneous and softer network, with a fourfold reduction in elastic modulus when compared to H[Formula: see text]O. Combining spectroscopic measurements with atomistic simulations, we show that collagen in D[Formula: see text]O is less hydrated than in H[Formula: see text]O. This partial dehydration lowers the enthalpic penalty for water removal and reorganization at the collagen-water interface, increasing the self-assembly rate and the number of nucleation centers, leading to thinner fibrils and a softer network. Coarse-grained simulations show that the acceleration in the initial nucleation rate can be reproduced by the enhancement of electrostatic interactions. These results show that water acts as a mediator between collagen monomers, by modulating their interactions so as to optimize the assembly process and, thus, the final network properties. We believe that isotopically modulating the hydration of proteins can be a valuable method to investigate the role of water in protein structural dynamics and protein self-assembly.

Unlocking novel therapies: cyclic peptide design for amyloidogenic targets through synergies of experiments, simulations, and machine learning

Existing therapies for neurodegenerative diseases like Parkinson's and Alzheimer's address only their symptoms and do not prevent disease onset. Common therapeutic agents, such as small molecules and antibodies struggle with insufficient selectivity, stability and bioavailability, leading to poor performance in clinical trials. Peptide-based therapeutics are emerging as promising candidates, with successful applications for cardiovascular diseases and cancers due to their high bioavailability, good efficacy and specificity. In particular, cyclic peptides have a long in vivo stability, while maintaining a robust antibody-like binding affinity. However, the de novo design of cyclic peptides is challenging due to the lack of long-lived druggable pockets of the target polypeptide, absence of exhaustive conformational distributions of the target and/or the binder, unknown binding site, methodological limitations, associated constraints (failed trials, time, money) and the vast combinatorial sequence space. Hence, efficient alignment and cooperation between disciplines, and synergies between experiments and simulations complemented by popular techniques like machine-learning can significantly speed up the therapeutic cyclic-peptide development for neurodegenerative diseases. We review the latest advancements in cyclic peptide design against amyloidogenic targets from a computational perspective in light of recent advancements and potential of machine learning to optimize the design process. We discuss the difficulties encountered when designing novel peptide-based inhibitors and we propose new strategies incorporating experiments, simulations and machine learning to design cyclic peptides to inhibit the toxic propagation of amyloidogenic polypeptides. Importantly, these strategies extend beyond the mere design of cyclic peptides and serve as template for the de novo generation of (bio)materials with programmable properties.

Decrypting Integrins by Mixed-Solvent Molecular Dynamics Simulations

Integrins are a family of α/β heterodimeric cell surface adhesion receptors which are capable of transmitting signals bidirectionally across membranes. They are known for their therapeutic potential in a wide range of diseases. However, the development of integrin-targeting medications has been impacted by unexpected downstream effects including unwanted agonist-like effects. Allosteric modulation of integrins is a promising approach to potentially overcome these limitations. Applying mixed-solvent molecular dynamics (MD) simulations to integrins, the current study uncovers hitherto unknown allosteric sites within the integrin α I domains of LFA-1 (α_Lβ₂; CD11a/CD18), VLA-1 (α₁β₁; CD49a/CD29), and Mac-1 (α_Mβ₂, CD11b/CD18). We show that these pockets are putatively accessible to small-molecule modulators. The findings reported here may provide opportunities for the design of novel allosteric integrin inhibitors lacking the unwanted agonism observed with earlier as well as current integrin-targeting drugs.

Antibody binding increases the flexibility of the prion protein

Prion diseases are associated with the conversion of the cellular prion protein (PrP) into a pathogenic conformer (PrP^Sc). A proposed therapeutic approach to avoid the pathogenic transformation is to develop antibodies that bind to PrP and stabilize its structure. POM1 and POM6 are two monoclonal antibodies that bind the globular domain of PrP and have different biological responses, i.e., trigger neurotoxicity mimicking prion infections (POM1) or prevent neurotoxicity (POM6). The crystal structures of PrP in complex with the two antibodies show similar epitopes which seems inconsistent with the opposite phenotypes. Here, we investigate the influence of the POM1 and POM6 antibodies on the flexibility of the mouse PrP by molecular dynamics simulations. The simulations reveal that the POM6/PrP interface is less stable than the POM1/PrP interface, ascribable to localized polar mismatches at the interface, despite the former complex having a larger epitope than the latter. In the presence of any of the two antibodies, the flexibility of the globular domain increases everywhere except for the β₁-α₁ loop in the POM1/PrP complex which suggests the involvement of this loop in the pathological conversion. The secondary structure of PrP is preserved whereas the polar interactions involving residues Glu146, Arg156 and Arg208 are modified upon antibody binding.

Interactions of Curcumin's Degradation Products with the Aβ42 Dimer: A Computational Study

Amyloid-β (Aβ) dimers are the smallest toxic species along the amyloid-aggregation pathway and among the most populated oligomeric accumulations present in the brain affected by Alzheimer's disease (AD). A proposed therapeutic strategy to avoid the aggregation of Aβ into higher-order structures is to develop molecules that inhibit the early stages of aggregation, i.e., dimerization. Under physiological conditions, the Aβ dimer is highly dynamic and does not attain a single well-defined structure but is rather characterized by an ensemble of conformations. In a recent study, a highly heterogeneous library of conformers of the Aβ dimer was generated by an efficient sampling method with constraints based on ion mobility mass spectrometry data. Here, we make use of the Aβ dimer library to study the interaction with two curcumin degradation products, ferulic aldehyde and vanillin, by molecular dynamics (MD) simulations. Ensemble docking and MD simulations are used to provide atomistic detail of the interactions between the curcumin degradation products and the Aβ dimer. The simulations show that the aromatic residues of Aβ, and in particular ¹⁹FF²⁰, interact with ferulic aldehyde and vanillin through π-π stacking. The binding of these small molecules induces significant changes on the ¹⁶KLVFF²⁰ region.

Antibody binding modulates the dynamics of the membrane-bound prion protein

Misfolding of the cellular prion protein (PrP^C) is associated with lethal neurodegeneration. PrP^C consists of a flexible tail (residues 23-123) and a globular domain (residues 124-231) whose C-terminal end is anchored to the cell membrane. The neurotoxic antibody POM1 and the innocuous antibody POM6 recognize the globular domain. Experimental evidence indicates that POM1 binding to PrP^C emulates the influence on PrP^C of the misfolded prion protein (PrP^Sc) while the binding of POM6 has the opposite biological response. Little is known about the potential interactions between flexible tail, globular domain, and the membrane. Here, we used atomistic simulations to investigate how these interactions are modulated by the binding of the Fab fragments of POM1 and POM6 to PrP^C and by interstitial sequence truncations to the flexible tail. The simulations show that the binding of the antibodies restricts the range of orientations of the globular domain with respect to the membrane and decreases the distance between tail and membrane. Five of the six sequence truncations influence only marginally this distance and the contact patterns between tail and globular domain. The only exception is a truncation coupled to a charge inversion mutation of four N-terminal residues, which increases the distance of the flexible tail from the membrane. The interactions of the flexible tail and globular domain are modulated differently by the two antibodies.

The role of the N-terminal amphipathic helix in bacterial YidC: Insights from functional studies, the crystal structure and molecular dynamics simulations

The evolutionary conserved YidC is a unique dual-function membrane protein that adopts insertase and chaperone conformations. The N-terminal helix of Escherichia coli YidC functions as an uncleaved signal sequence and is important for membrane insertion and interaction with the Sec translocon. Here, we report the first crystal structure of Thermotoga maritima YidC (TmYidC) including the N-terminal amphipathic helix (N-AH) (PDB ID: 6Y86). Molecular dynamics simulations show that N-AH lies on the periplasmic side of the membrane bilayer forming an angle of about 15° with the membrane surface. Our functional studies suggest a role of N-AH for the species-specific interaction with the Sec translocon. The reconstitution data and the superimposition of TmYidC with known YidC structures suggest an active insertase conformation for YidC. Molecular dynamics (MD) simulations of TmYidC provide evidence that N-AH acts as a membrane recognition helix for the YidC insertase and highlight the flexibility of the C1 region underlining its ability to switch between insertase and chaperone conformations. A structure-based model is proposed to rationalize how YidC performs the insertase and chaperone functions by re-positioning of N-AH and the other structural elements.

Disorder at the Tips of a Disease-Relevant Aβ42 Amyloid Fibril: A Molecular Dynamics Study

We present a simulation study of the early events of peptide dissociation from a fibril of the Alzheimer's Aβ42 peptide. The fibril consists of layers of two adjacent Aβ42 peptides each folded in an S-shaped structure which has been determined by solid state NMR spectroscopy of a monomorphic disease-relevant species. Multiple molecular dynamics runs (16 at 310 K and 15 at 370 K) were carried out starting from an 18-peptide protofibril for a cumulative sampling of about 15 μs. The simulations show structural stability of the fibrillar core and an overall increase in the twist to about 3 degrees. The N-terminal segment 1-14 is disordered in all peptides. At both ends of the fibril, the central segment 21-29, which includes part of the β₂ strand, dissociates in some of the simulations. The β₁ and β₃ strands, residues 15-20 and 35-41, respectively, are structurally stable. The transient binding of the N-terminal stretch to the β₃ strand of the adjacent peptide at the tip is likely to contribute to the arrest phase of the stop-and-go mechanism.

The attachment of α-synuclein to a fiber: A coarse-grain approach

We present simulations of the amyloidogenic core of α-synuclein, the protein causing Parkinson's disease, as a short chain of coarse-grain patchy particles. Each particle represents a sequence of about a dozen amino acids. The fluctuating secondary structure of this intrinsically disordered protein is modelled by dynamic variations of the shape and interaction characteristics of the patchy particles, ranging from spherical with weak isotropic attractions for the disordered state to spherocylindrical with strong directional interactions for a β-sheet. Flexible linkers between the particles enable sampling of the tertiary structure. This novel model is applied here to study the growth of an amyloid fibril, by calculating the free energy profile of a protein attaching to the end of a fibril. The simulation results suggest that the attaching protein readily becomes trapped in a mis-folded state, thereby inhibiting further growth of the fibril until the protein has readjusted to conform to the fibril structure, in line with experimental findings and previous simulations on small fragments of other proteins.

Intrinsic Conformational Preferences and Interactions in α-Synuclein Fibrils: Insights from Molecular Dynamics Simulations

Amyloid formation by the intrinsically disordered α-synuclein protein is the hallmark of Parkinson's disease. We present atomistic Molecular Dynamics simulations of the core of α-synuclein using enhanced sampling techniques to describe the conformational and binding free energy landscapes of fragments implicated in fibril stabilization. The theoretical framework is derived to combine the free energy profiles of the fragments into the reaction free energy of a protein binding to a fibril. Our study shows that individual fragments in solution have a propensity toward attaining non-β conformations, indicating that in a fibril β-strands are stabilized by interactions with other strands. We show that most dimers of hydrogen-bonded fragments are unstable in solution, while hydrogen bonding stabilizes the collective binding of five fragments to the end of a fibril. Hydrophobic effects make further contributions to the stability of fibrils. This study is the first of its kind where structural and binding preferences of the five major fragments of the hydrophobic core of α-synuclein have been investigated. This approach improves sampling of intrinsically disordered proteins, provides information on the binding mechanism between the core sequences of α-synuclein, and enables the parametrization of coarse grained models.