Multi-eGO: Model Improvements toward the Study of Complex Self-Assembly Processes.

Structure-based models have been instrumental in simulating protein folding and suggesting hypotheses about the mechanisms involved. Nowadays, at least for fast-folding proteins, folding can be simulated in explicit solvent using classical molecular dynamics. However, other self-assembly processes, such as protein aggregation, are still far from being accessible. Recently, we proposed that a hybrid multistate structure-based model, multi-eGO, could help to bridge the gap toward the simulation of out-of-equilibrium, concentration-dependent self-assembly processes. Here, we further improve the model and show how multi-eGO can effectively and accurately learn the conformational ensemble of the amyloid ß42 intrinsically disordered peptide, reproduce the well-established folding mechanism of the B1 immunoglobulin-binding domain of streptococcal protein G, and reproduce the aggregation as a function of the concentration of the transthyretin 105-115 amyloidogenic peptide. We envision that by learning from the dynamics of a few minima, multi-eGO can become a platform for simulating processes inaccessible to other simulation techniques.

Multi-eGO: An in silico lens to look into protein aggregation kinetics at atomic resolution.

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-eGO, 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-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 carried out on an amyloidogenic peptide from Transthyretin, a protein responsible for one of the most common cardiac amyloidoses. Multi-eGO 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-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.

Inactivation of the Pyrimidine Biosynthesis pyrD Gene Negatively Affects Biofilm Formation and Virulence Determinants in the Crohn's Disease-Associated Adherent Invasive Escherichia coli LF82 Strain.

In Crohn's disease (CD) patients, the adherent-invasive Escherichia coli (AIEC) pathovar contributes to the chronic inflammation typical of the disease via its ability to invade gut epithelial cells and to survive in macrophages. We show that, in the AIEC strain LF82, inactivation of the pyrD gene, encoding dihydroorotate dehydrogenase (DHOD), an enzyme of the de novo pyrimidine biosynthetic pathway, completely abolished its ability of to grow in a macrophage environment-mimicking culture medium. In addition, pyrD inactivation reduced flagellar motility and strongly affected biofilm formation by downregulating transcription of both type 1 fimbriae and curli subunit genes. Thus, the pyrD gene appears to be essential for several cellular processes involved in AIEC virulence. Interestingly, vidofludimus (VF), a DHOD inhibitor, has been proposed as an effective drug in CD treatment. Despite displaying a potentially similar binding mode for both human and E. coli DHOD in computational molecular docking experiments, VF showed no activity on either growth or virulence-related processes in LF82. Altogether, our results suggest that the crucial role played by the pyrD gene in AIEC virulence, and the presence of structural differences between E. coli and human DHOD allowing for the design of specific inhibitors, make E. coli DHOD a promising target for therapeutical strategies aiming at counteracting chronic inflammation in CD by acting selectively on its bacterial triggers.

Retromer stabilization results in neuroprotection in a model of Amyotrophic Lateral Sclerosis.

Amyotrophic Lateral Sclerosis (ALS) is a fatal disease characterized by the degeneration of upper and lower motor neurons (MNs). We find a significant reduction of the retromer complex subunit VPS35 in iPSCs-derived MNs from ALS patients, in MNs from ALS post mortem explants and in MNs from SOD1G93A mice. Being the retromer involved in trafficking of hydrolases, a pathological hallmark in ALS, we design, synthesize and characterize an array of retromer stabilizers based on bis-guanylhydrazones connected by a 1,3-phenyl ring linker. We select compound 2a as a potent and bioavailable interactor of VPS35-VPS29. Indeed, while increasing retromer stability in ALS mice, compound 2a attenuates locomotion impairment and increases MNs survival. Moreover, compound 2a increases VPS35 in iPSCs-derived MNs and shows brain bioavailability. Our results clearly suggest the retromer as a valuable druggable target in ALS.

The structure of N184K amyloidogenic variant of gelsolin highlights the role of the H-bond network for protein stability and aggregation properties.

Mutations in the gelsolin protein are responsible for a rare conformational disease known as AGel amyloidosis. Four of these mutations are hosted by the second domain of the protein (G2): D187N/Y, G167R and N184K. The impact of the latter has been so far evaluated only by studies on the isolated G2. Here we report the characterization of full-length gelsolin carrying the N184K mutation and compare the findings with those obtained on the wild type and the other variants. The crystallographic structure of the N184K variant in the Ca2+-free conformation shows remarkable similarities with the wild type protein. Only minimal local rearrangements can be observed and the mutant is as efficient as the wild type in severing filamentous actin. However, the thermal stability of the pathological variant is compromised in the Ca2+-free conditions. These data suggest that the N to K substitution causes a local disruption of the H-bond network in the core of the G2 domain. Such a subtle rearrangement of the connections does not lead to significant conformational changes but severely affects the stability of the protein.

High-resolution crystal structure of gelsolin domain 2 in complex with the physiological calcium ion.

The second domain of gelsolin (G2) hosts mutations responsible for a hereditary form of amyloidosis. The active form of gelsolin is Ca2+-bound; it is also a dynamic protein, hence structural biologists often rely on the study of the isolated G2. However, the wild type G2 structure that have been used so far in comparative studies is bound to a crystallographic Cd2+, in lieu of the physiological calcium. Here, we report the wild type structure of G2 in complex with Ca2+ highlighting subtle ion-dependent differences. Previous findings on different G2 mutations are also briefly revised in light of these results.