Search

The uneven worldwide vaccination coverage against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and emergence of variants escaping immunity call for broadly effective and easily deployable therapeutic agents. We have previously described the human single-chain scFv76 antibody, which recognizes SARS-CoV-2 Alpha, Beta, Gamma and Delta variants. We now show that scFv76 also neutralizes the infectivity and fusogenic activity of the Omicron BA.1 and BA.2 variants. Cryoelectron microscopy (cryo-EM) analysis reveals that scFv76 binds to a well-conserved SARS-CoV-2 spike epitope, providing the structural basis for its broad-spectrum activity. We demonstrate that nebulized scFv76 has therapeutic efficacy in a severe hACE2 transgenic mouse model of coronavirus disease 2019 (COVID-19) pneumonia, as shown by body weight and pulmonary viral load data. Counteraction of infection correlates with inhibition of lung inflammation, as observed by histopathology and expression of inflammatory cytokines and chemokines. Biomarkers of pulmonary endothelial damage were also significantly reduced in scFv76-treated mice. The results support use of nebulized scFv76 for COVID-19 induced by any SARS-CoV-2 variants that have emerged so far.

Protein aggregation into amyloid fibrils is the archetype of aberrant biomolecular self-assembly processes, with more than 50 associated diseases that are mostly uncurable. Understanding aggregation mechanisms is thus of fundamental importance and goes in parallel with the structural characterization of the transient oligomers formed during the process. Oligomers have been proven elusive to high-resolution structural techniques, while the large sizes and long time scales, typical of aggregation processes, have limited the use of computational methods to date. To surmount these limitations, we here present multi- e GO, an atomistic, hybrid structure-based model which, leveraging the knowledge of monomers conformational dynamics and of fibril structures, efficiently captures the essential structural and kinetics aspects of protein aggregation. Multi- e GO molecular dynamics simulations can describe the aggregation kinetics of thousands of monomers. The concentration dependence of the simulated kinetics, as well as the structural features of the resulting fibrils, are in qualitative agreement with in vitro experiments carried out on an amyloidogenic peptide from Transthyretin, a protein responsible for one of the most common cardiac amyloidoses. Multi- e GO simulations allow the formation of primary nuclei in a sea of transient lower-order oligomers to be observed over time and at atomic resolution, following their growth and the subsequent secondary nucleation events, until the maturation of multiple fibrils is achieved. Multi- e GO, combined with the many experimental techniques deployed to study protein aggregation, can provide the structural basis needed to advance the design of molecules targeting amyloidogenic diseases.

The stabilization of native states of proteins is a powerful drug discovery strategy. It is still unclear, however, whether this approach can be applied to intrinsically disordered proteins. Here, we report a small molecule that stabilizes the native state of the Aβ42 peptide, an intrinsically disordered protein fragment associated with Alzheimer's disease. We show that this stabilization takes place by a disordered binding mechanism, in which both the small molecule and the Aβ42 peptide remain disordered. This disordered binding mechanism involves enthalpically favorable local π-stacking interactions coupled with entropically advantageous global effects. These results indicate that small molecules can stabilize disordered proteins in their native states through transient non-specific interactions that provide enthalpic gain while simultaneously increasing the conformational entropy of the proteins.

Recognition motifs that mediate protein-protein interactions are usually embedded within longer intrinsically disordered regions. While binding interfaces involving the recognition motif in such interactions are well studied, less is known about the role of disordered regions flanking the motifs. The interaction between the transcriptional co-activators NCOA3 (ACTR) and CBP is mediated by coupled binding and folding of the two domains CID and NCBD. Here, we used circular dichroism and kinetics to directly quantify the contribution of the adjacent flanking regions of CID to its interaction with NCBD. Using N-and C terminal combinatorial variants we found that the flanking regions promote binding in an additive fashion while retaining a large degree of disorder in the complex. Experiments at different ionic strengths demonstrated that the increase in affinity is not mediated by electrostatic interactions from the flanking regions. Instead, site-directed mutagenesis and molecular dynamics simulations suggest that binding is promoted by short-lived non-specific hydrophobic contacts between the flanking regions and NCBD. Our findings are consistent with highly frustrated interactions outside of the canonical binding interface resulting in a slightly energetically favorable fuzzy binding. Modulation of affinity via flanking regions could represent a general mechanism for functional regulation by intrinsically disordered protein regions.(c) 2022 The Author(s). Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Protein aggregation into amyloid fibrils is the archetype of aberrant biomolecular self-assembly processes, with more than 50 diseases associated that are mostly uncurable. Understanding aggregation mechanisms is thus of fundamental importance and goes in parallel with the characterization of the structures of the transient oligomers formed in the process. Oligomers have been proven elusive to high-resolution structural techniques, while the large sizes and long-time scales typical of aggregation processes have limited, so far, the use of computational methods. To surmount these limitations, we introduce here multi-eGO, an atomistic, hybrid structure-based model, which leveraging on the knowledge of monomers conformational dynamics and of fibril structures, can efficiently capture the essential structural and kinetics aspects of protein aggregation. Multi-eGO molecular dynamics simulations can describe the aggregation kinetics of thousands of monomers. The concentration dependence of the simulated kinetics, as well as the structural features of the resulting fibrils, are in qualitative agreement with in vitro experiments on an amyloidogenic peptide of Transthyretin, a protein responsible for one of the most common cardiac amyloidosis. Multi-eGO simulations allow to observe in time and at atomic resolution the formation of primary nuclei in a sea of transient lower order oligomers, to follow their growth and the subsequent secondary nucleation events, till the maturation of multiple fibrils. Multi-eGO, combined with the many experimental techniques deployed to study protein aggregation, can provide the structural basis needed to advance the design of molecules targeting amyloidogenic diseases.Significance StatementAlzheimer’s and Parkinson’s diseases are uncurable pathologies associated to the aberrant aggregation of specific proteins into amyloid fibrils. Understanding the mechanism leading to protein aggregation, by characterizing the structures of the oligomeric species populated in the process, would have a tremendous impact on the design of therapeutic molecules. We propose that a structure-based approach to molecular dynamics simulations can allow following at high resolution the aggregation kinetics of thousands of monomers. Having shown that simulations can describe the aggregation of a Transthyretin amyloidogenic peptide, we demonstrate how their efficiency allows acquiring a wealth of structural information. We foresee that integrating the latter with the many techniques developed to study protein aggregation will support the design of molecules to modulate amyloidogenesis.

Presentation of pathogen-derived epitopes by major histocompatibility complex I (MHC-I) can lead to the activation and expansion of specific CD8

The stabilisation of native states of proteins is a powerful drug discovery strategy. It is still unclear, however, whether this approach can be applied to intrinsically disordered proteins. Here we report a small molecule that stabilises the native state of the Aβ42 peptide, an intrinsically disordered protein fragment associated with Alzheimer’s disease. We show that this stabilisation takes place by a dynamic binding mechanism, in which both the small molecule and the Aβ42 peptide remain disordered. This disordered binding mechanism involves enthalpically favourable local π-stacking interactions coupled with entropically advantageous global effects. These results indicate that small molecules can stabilise disordered proteins in their native states through transient non-specific interactions that provide enthalpic gain while simultaneously increasing the conformational entropy of the proteins.

Bacterial NADPH-dependent glutamate synthase (GltS) is a complex iron-sulfur flavoprotein that catalyzes the reductive synthesis of two L-Glu molecules from L-Gln and 2-oxo-glutarate. GltS functional unit hosts an α-subunit (αGltS) and a β-subunit (βGltS) that assemble in different αβ oligomers in solution. Here, we present the cryo-electron microscopy structures of Azospirillum brasilense GltS in four different oligomeric states (α4β3, α4β4, α6β4 and α6β6, in the 3.5- to 4.1-A resolution range). Our study provides a comprehensive GltS model that details the inter-protomeric assemblies and allows unequivocal location of the FAD cofactor and of two electron transfer [4Fe-4S]+1,+2 clusters within βGltS.

To survive in cold environments, psychrophilic organisms produce enzymes endowed with high specific activity at low temperature. The structure of these enzymes is usually flexible and mostly thermolabile. In this work, we investigate the structural basis of cold adaptation of a GH42 β-galactosidase from the psychrophilic Marinomonas ef1. This enzyme couples cold activity with astonishing robustness for a psychrophilic protein, for it retains 23% of its highest activity at 5°C and it is stable for several days at 37°C and even 50°C. Phylogenetic analyses indicate a close relationship with thermophilic β-galactosidases, suggesting that the present-day enzyme evolved from a thermostable scaffold modeled by environmental selective pressure. The crystallographic structure reveals the overall similarity with GH42 enzymes, along with a hexameric arrangement (dimer of trimers) not found in psychrophilic, mesophilic, and thermophilic homologues. In the quaternary structure, protomers form a large central cavity, whose accessibility to the substrate is promoted by the dynamic behavior of surface loops, even at low temperature. A peculiar cooperative behavior of the enzyme is likely related to the increase of the internal cavity permeability triggered by heating. Overall, our results highlight a novel strategy of enzyme cold adaptation, based on the oligomerization state of the enzyme, which effectively challenges the paradigm of cold activity coupled with intrinsic thermolability. DATABASE: Structural data are available in the Protein Data Bank database under the accession number 6Y2K.

MJ0366 from Methanocaldococcus jannaschii is the smallest topologically knotted protein known to date. 92 residues in length, MJ0366 ties a trefoil (3(1)) knot by threading its C-terminal helix through a buttonhole formed by the remainder of the secondary structure elements. By generating a library of point mutations at positions pertinent to the knot formation, we systematically evaluated the contributions of individual residues to the folding stability and kinetics of MJ0366. The experimental Φ-values were used as restraints to computationally generate an ensemble of conformations that correspond to the transition state of MJ0366, which revealed several nonnative contacts. The importance of these nonnative contacts in stabilizing the transition state of MJ0366 was confirmed by a second round of mutagenesis, which also established the pivotal role of F15 in stapling the network of hydrophobic interactions around the threading C-terminal helix. Our converging experimental and computational results show that, despite the small size, the transition state of MJ0366 is formed at a very late stage of the folding reaction coordinate, following a polarized pathway. Eventually, the formation of extensive native contacts, as well as a number of nonnative ones, leads to the threading of the C-terminal helix that defines the topological knot.

NF-Y is a transcription factor (TF) comprising three subunits (NF-YA, NF-YB, NF-YC) that binds with high specificity to the CCAAT sequence, a widespread regulatory element in gene promoters of prosurvival, cell-cycle-promoting, and metabolic genes. Tumor cells undergo &ldquo, metabolic rewiring&rdquo, through overexpression of genes involved in such pathways, many of which are under NF-Y control. In addition, NF-YA appears to be overexpressed in many tumor types. Thus, limiting NF-Y activity may represent a desirable anti-cancer strategy, which is an ongoing field of research. With virtual-screening docking simulations on a library of pharmacologically active compounds, we identified suramin as a potential NF-Y inhibitor. We focused on suramin given its high water-solubility that is an important factor for in vitro testing, since NF-Y is sensitive to DMSO. By electrophoretic mobility shift assays (EMSA), isothermal titration calorimetry (ITC), STD NMR, X-ray crystallography, and molecular dynamics (MD) simulations, we showed that suramin binds to the histone fold domains (HFDs) of NF-Y, preventing DNA-binding. Our analyses, provide atomic-level detail on the interaction between suramin and NF-Y and reveal a region of the protein, nearby the suramin-binding site and poorly conserved in other HFD-containing TFs, that may represent a promising starting point for rational design of more specific and potent inhibitors with potential therapeutic applications.

Intrinsically disordered protein domains often have multiple binding partners. It is plausible that the strength of pairing with specific partners evolves from an initial low affinity to a higher affinity. However, little is known about the molecular changes in the binding mechanism that would facilitate such a transition. We previously showed that the interaction between two intrinsically disordered domains, NCBD and CID, likely emerged in an ancestral deuterostome organism as a low-affinity interaction that subsequently evolved into a higher-affinity interaction before the radiation of modern vertebrate groups. Here we map native contacts in the transition states of the low-affinity ancestral and high-affinity human NCBD/CID interactions. We show that the coupled binding and folding mechanism is overall similar but with a higher degree of native hydrophobic contact formation in the transition state of the ancestral complex and more heterogeneous transient interactions, including electrostatic pairings, and an increased disorder for the human complex. Adaptation to new binding partners may be facilitated by this ability to exploit multiple alternative transient interactions while retaining the overall binding and folding pathway.

The use of small-angle X-ray scattering (SAXS) in combination with molecular dynamics simulation is hampered by its heavy computational cost. The calculation of SAXS from atomic structures can be speeded up by using a coarse-grain representation of the structure. Following the work of Niebling, Björling & Westenhoff [J. Appl. Cryst. (2014), 47, 1190–1198], the Martini bead form factors for nucleic acids have been derived and then implemented, together with those previously determined for proteins, in the publicly available PLUMED library. A hybrid multi-resolution strategy has also been implemented to perform SAXS restrained simulations at atomic resolution by calculating the virtual positions of the Martini beads on the fly and using them for the calculation of SAXS. The accuracy and efficiency of the method are demonstrated by refining the structure of two protein–nucleic acid complexes. Instrumental for this result is the use of metainference, which allows the consideration and alleviation of the approximations at play in the present SAXS calculations.

The discovery that disordered proteins are widespread in the human proteome has prompted the quest for methods to characterize the conformational properties that determine their functional and dysfunctional behaviour. It has become customary to describe these proteins in terms of structural ensembles and free energy landscapes, which offer conformational and thermodynamic insight. However, a current major challenge is to generalize this description to ‘kinetic ensembles’, thereby also providing information on transition rates between states. Approaches based on the theory of stochastic processes can be particularly suitable for this purpose. Here, we develop a Markov state model and illustrate its application by determining a kinetic ensemble of the 42-residue form of the amyloid-β peptide (Aβ42), whose aggregation is associated with Alzheimer’s disease. Using the Google Compute Engine, we generated 315 μs all-atom, explicit solvent molecular dynamics trajectories, validated with experimental data from nuclear magnetic resonance spectroscopy. Using a probabilistic-based definition of conformational states in a neural network approach, we found that Aβ42 is characterized by inter-state transitions no longer than the microsecond timescale, exhibiting only fully unfolded or short-lived, partially-folded states. We contextualize our findings by performing additional simulations of the oxidized form of Aβ42. Our results illustrate how the use of kinetic ensembles offers an effective means to provide information about the structure, thermodynamics, and kinetics of disordered proteins towards an understanding of these ubiquitous biomolecules.

A growing body of evidences has established that in many cases proteins may preserve most of their function and flexibility in a crystalline environment, and several techniques are today capable to detect transiently-populated states of macromolecules in tightly packed lattices. Intriguingly, in the case of amyloidogenic precursors, the presence of these conformations (hidden to conventional crystallographic studies) can be correlated to the pathological fate of the native fold.It remains unclear, however, to which extent these minor conformations reflect the protein behaviour that is more commonly studied in solution. Here, we address this question by investigating some biophysical properties of a prototypical amyloidogenic system, β2-microglobulin (β2m) in solution and in microcrystalline state.By combining NMR chemical shifts with Molecular Dynamics (MD) simulations, we confirmed that conformational dynamics of β2m native state in the crystal lattice is in keeping with what observed in solution.A comparative study of protein stability in solution andin crystallois then carried out, monitoring the change in protein secondary structure at increasing temperature by Fourier transform infrared (FTIR) spectroscopy. The increased structural order of the crystalline state contributes to provide better resolved spectral components compared to those collected in solution and crucially, the crystalline samples display thermal stabilities in good agreement with the trend observed in solution.Overall, this work shows that protein stability and occurrence of pathological hidden states in crystals parallel their solution counterpart, confirming the interest of crystals as a platform for the biophysical characterisation of processes such as unfolding and aggregation.

Poly-ubiquitin chains are flexible multidomain proteins, whose conformational dynamics enable their molecular recognition by a large number of partners in multiple biological pathways. By using alternative linkage, it is possible to obtain poly-ubiquitin molecules with different dynamical properties. This flexibility is further increased by the possibility to tune the length of poly-ubiquitin chains. Characterizing the dynamics of poly-ubiquitins as a function of their length is thus relevant to understand their biology. Structural characterization of poly-ubiquitin conformational dynamics is challenging both experimentally and computationally due to increasing system size and conformational variability. Here, by developing highly efficient and accurate small-angle X-ray scattering driven Martini coarse-grain simulations, we characterize the dynamics of linear M1-linked di-, tri- and tetra-ubiquitin chains. Our data show that the behavior of the di-ubiquitin subunits is independent of the presence of additional ubiquitin modules. We propose that the conformational space sampled by linear poly-ubiquitins, in general, may follow a simple self-avoiding polymer model. These results, combined with experimental data from small angle X-ray scattering, biophysical techniques and additional simulations show that binding of NEMO, a central regulator in the NF-κB pathway, to linear poly-ubiquitin obeys a 2:1 (NEMO:poly-ubiquitin) stoichiometry in solution, even in the context of four ubiquitin units. Eventually, we show how the conformational properties of long poly-ubiquitins may modulate the binding with their partners in a length-dependent manner.SignificanceProtein conformational dynamics plays an essential role in molecular recognition mechanisms. The characterization of conformational dynamics is hampered by the conformational averaging of observable in experimental structural biology techniques and by the limitations in the accuracy of computational methods. By developing an efficient and accurate approach to combine small-angle X-ray scattering solution experiments and coarse-grain Martini simulations, we show that the conformational dynamics of linear poly-ubiquitins can be efficiently determined and to rationalize the role of poly-ubiquitin dynamic in the molecular recognition of the UBAN domain upon binding to the signaling regulator NEMO. The analysis of the conformational ensembles allows us to propose a general model of the dynamics of linear poly-ubiquitin chains where they can be described as a self-avoiding polymer with a characteristic length associated with their specific linkage.

Significance Efficient synthesis and folding of proteins, avoiding misfolded states, are central to cell function. As folding may be initiated in parallel with translation, key experimental challenges are to map changes that occur in folding free energy landscapes as translation proceeds and to understand how these landscapes might be modulated by the ribosome and auxiliary factors. Here, we study the length-dependent folding of a domain from a tandem repeat protein and solve the structure of a stable folding intermediate. Although destabilized by the ribosome at equilibrium, modeling of the nonequilibrium folding pathway nevertheless indicates a significant role for proline isomerization during translation. We develop a simple model to explore the impact of cotranslational folding kinetics on misfolding hazards., Cotranslational folding (CTF) is a fundamental molecular process that ensures efficient protein biosynthesis and minimizes the formation of misfolded states. However, the complexity of this process makes it extremely challenging to obtain structural characterizations of CTF pathways. Here, we correlate observations of translationally arrested nascent chains with those of a systematic C-terminal truncation strategy. We create a detailed description of chain length-dependent free energy landscapes associated with folding of the FLN5 filamin domain, in isolation and on the ribosome, and thus, quantify a substantial destabilization of the native structure on the ribosome. We identify and characterize two folding intermediates formed in isolation, including a partially folded intermediate associated with the isomerization of a conserved cis proline residue. The slow folding associated with this process raises the prospect that neighboring unfolded domains might accumulate and misfold during biosynthesis. We develop a simple model to quantify the risk of misfolding in this situation and show that catalysis of folding by peptidyl-prolyl isomerases is sufficient to eliminate this hazard.

Proteins employ the information stored in the genetic code and translated into their sequences to carry out well-defined functions in the cellular environment. The possibility to encode for such functions is controlled by the balance between the amount of information supplied by the sequence and that left after that the protein has folded into its structure. We developed a computational algorithm to evaluate the amount of information necessary to specify the protein structure, keeping into account the thermodynamic properties of protein folding. We thus show that the information remaining in the protein sequence after encoding for its structure (the 'information gap') is very close to what needed to encode for its function and interactions. Then, by predicting the information gap directly from the protein sequence, we show that it may be possible to use these insights from information theory to discriminate between ordered and disordered proteins, to identify unknown functions, and to optimize designed proteins sequences.

Small molecules with rotatable bonds can occupy different conformational states in solution as a consequence of their thermal fluctuations. The accurate determination of the structures of such states, as well as of their statistical weights, has been challenging because of the technical difficulties in extracting information from experimental measurements, which are normally averaged over the conformational space available. Here, to achieve this objective, we present an approach based on a recently proposed tensor-free method for incorporating NMR residual dipolar couplings as structural restraints in replica-averaged molecular dynamics simulations. This approach enables the information provided by the experimental data to be used in the spirit of the maximum entropy principle to determine the structural ensembles of small molecules. Furthermore, in order to enhance the sampling of the conformational space we incorporated the metadynamics method in the simulations. We illustrate the method in the case of strychnine, determining the three major conformational states of this small molecule and their associated occupation probabilities.

The PLUMED consortium unifies developers and contributors to PLUMED, an open-source library for enhanced-sampling, free-energy calculations and the analysis of molecular dynamics simulations. Here, we outline our efforts to promote transparency and reproducibility by disseminating protocols for enhanced-sampling molecular simulations.

The nucleation-condensation mechanism represents a major paradigm to understand the folding process of many small globular proteins. Although substantial evidence has been acquired for this mechanism, it has remained very challenging to characterize the initial events leading to the formation of a folding nucleus. To achieve this goal, we used a combination of relaxation dispersion NMR spectroscopy and molecular dynamics simulations to determine ensembles of conformations corresponding to the denatured, transition, and native states in the folding of the activation domain of human procarboxypeptidase A2 (ADA2h). We found that the residues making up the folding nucleus tend to interact in the denatured state in a transient manner and not simultaneously, thereby forming incomplete and distorted versions of the folding nucleus. Only when all the contacts between these key residues are eventually formed can the protein reach the transition state and continue folding. Overall, our results elucidate the mechanism of formation of the folding nucleus of a protein and provide insights into how its folding rate can be modified during evolution by mutations that modulate the strength of the interactions between the residues forming the folding nucleus.

Intrinsically disordered protein (IDP) domains often have multiple binding partners. Little is known regarding molecular changes in the binding mechanism when a new interaction evolves from low to high affinity. Here we compared the degree of native contacts in the transition state of the interaction of two IDP domains, low-affinity ancestral and high-affinity human NCBD and CID. We found that the coupled binding and folding mechanism of the domains is overall similar, but with a higher degree of native hydrophobic contact formation in the transition state of the ancestral complex while more heterogenous transient interactions, including electrostatic, and an increased disorder characterize the human complex. From an evolutionary perspective, adaptation to new binding partners for IDPs may benefit from this ability to exploit multiple alternative transient interactions while retaining the overall pathway.

A small molecule binds to a disordered protein in its monomeric form, preventing its aggregation linked to Alzheimer’s disease., Disordered proteins are challenging therapeutic targets, and no drug is currently in clinical use that modifies the properties of their monomeric states. Here, we identify a small molecule (10074-G5) capable of binding and sequestering the intrinsically disordered amyloid-β (Aβ) peptide in its monomeric, soluble state. Our analysis reveals that this compound interacts with Aβ and inhibits both the primary and secondary nucleation pathways in its aggregation process. We characterize this interaction using biophysical experiments and integrative structural ensemble determination methods. We observe that this molecule increases the conformational entropy of monomeric Aβ while decreasing its hydrophobic surface area. We also show that it rescues a Caenorhabditis elegans model of Aβ-associated toxicity, consistent with the mechanism of action identified from the in silico and in vitro studies. These results illustrate the strategy of stabilizing the monomeric states of disordered proteins with small molecules to alter their behavior for therapeutic purposes.

SAXS experiments provide low-resolution but valuable information about the dynamics of biomolecular systems, which could be ideally integrated into MD simulations to accurately determine conformational ensembles of flexible proteins. The applicability of this strategy is hampered by the high computational cost required to calculate scattering intensities from three-dimensional structures. We previously presented a hybrid resolution method that makes atomistic SAXS-restrained MD simulation feasible by adopting a coarse-grained approach to efficiently back-calculate scattering intensities; here, we extend this technique, applying it in the framework of metainference with the aim to investigate the dynamical behavior of flexible biomolecules. The efficacy of the method is assessed on the K63-diubiquitin, showing that the inclusion of SAXS-restraints is effective in generating a reliable conformational ensemble, improving the agreement with independent experimental data.

Designed peptides derived from the islet amyloid polypeptide (IAPP) cross‐amyloid interaction surface with Aβ (termed interaction surface mimics or ISMs) have been shown to be highly potent inhibitors of Aβ amyloid self‐assembly. However, the molecular mechanism of their function is not well understood. Using solution‐state and solid‐state NMR spectroscopy in combination with ensemble‐averaged dynamics simulations and other biophysical methods including TEM, fluorescence spectroscopy and microscopy, and DLS, we characterize ISM structural preferences and interactions. We find that the ISM peptide R3‐GI is highly dynamic, can adopt a β‐like structure, and oligomerizes into colloid‐like assemblies in a process that is reminiscent of liquid–liquid phase separation (LLPS). Our results suggest that such assemblies yield multivalent surfaces for interactions with Aβ40. Sequestration of substrates into these colloid‐like structures provides a mechanistic basis for ISM function and the design of novel potent anti‐amyloid molecules., ISM inhibitors are highly dynamic assemblies and exchange between monomeric and high‐molecular‐weight states. In both states, the ISM peptide adopts a β‐loop‐like structure that provides a suitable surface for sequestration of Aβ40 in the colloidal state. ISM inhibitors thus exploit multivalency by self‐association to yield high substrate avidity.

Accurate protein structural ensembles can be determined with metainference, a Bayesian inference method that integrates experimental information with prior knowledge of the system and deals with all sources of uncertainty and errors as well as with system heterogeneity. Furthermore, metainference can be implemented using the metadynamics approach, which enables the computational study of complex biological systems requiring extensive conformational sampling. In this chapter, we provide a step-by-step guide to perform and analyse metadynamic metainference simulations using the ISDB module of the open-source PLUMED library, as well as a series of practical tips to avoid common mistakes. Specifically, we will guide the reader in the process of learning how to model the structural ensemble of a small disordered peptide by combining state-of-the-art molecular mechanics force fields with nuclear magnetic resonance data, including chemical shifts, scalar couplings and residual dipolar couplings.