Fusing simulation and experiment: The effect of mutations on the structure and activity of the influenza fusion peptide

During the infection process, the influenza fusion peptide (FP) inserts into the host membrane, playing a crucial role in the fusion process between the viral and host membranes. In this work we used a combination of simulation and experimental techniques to analyse the molecular details of this process, which are largely unknown. Although the FP structure has been obtained by NMR in detergent micelles, there is no atomic structure information in membranes. To answer this question, we performed bias-exchange metadynamics (BE-META) simulations, which showed that the lowest energy states of the membrane-inserted FP correspond to helical-hairpin conformations similar to that observed in micelles. BE-META simulations of the G1V, W14A, G12A/G13A and G4A/G8A/G16A/G20A mutants revealed that all the mutations affect the peptide's free energy landscape. A FRET-based analysis showed that all the mutants had a reduced fusogenic activity relative to the WT, in particular the mutants G12A/G13A and G4A/G8A/G16A/G20A. According to our results, one of the major causes of the lower activity of these mutants is their lower membrane affinity, which results in a lower concentration of peptide in the bilayer. These findings contribute to a better understanding of the influenza fusion process and open new routes for future studies

Using metadynamics to explore complex free-energy landscapes

Metadynamics is an atomistic simulation technique that allows, within the same framework, acceleration of rare events and estimation of the free energy of complex molecular systems. It is based on iteratively \u2018filling\u2019 the potential energy of the system by a sum of Gaussians centred along the trajectory followed by a suitably chosen set of collective variables (CVs), thereby forcing the system to migrate from one minimum to the next. The power of metadynamics is demonstrated by the large number of extensions and variants that have been developed. The first scope of this Technical Review is to present a critical comparison of these variants, discussing their advantages and disadvantages. The effectiveness of metadynamics, and that of the numerous alternative methods, is strongly influenced by the choice of the CVs. If an important variable is neglected, the resulting estimate of the free energy is unreliable, and predicted transition mechanisms may be qualitatively wrong. The second scope of this Technical Review is to discuss how the CVs should be selected, how to verify whether the chosen CVs are sufficient or redundant, and how to iteratively improve the CVs using machine learning approaches

The inverted free energy landscape of an intrinsically disordered peptide by simulations and experiments

The free energy landscape theory has been very successful in rationalizing the folding behaviour of globular proteins, as this representation provides intuitive information on the number of states involved in the folding process, their populations and pathways of interconversion. We extend here this formalism to the case of the A\u3b240 peptide, a 40-residue intrinsically disordered protein fragment associated with Alzheimer's disease. By using an advanced sampling technique that enables free energy calculations to reach convergence also in the case of highly disordered states of proteins, we provide a precise structural characterization of the free energy landscape of this peptide. We find that such landscape has inverted features with respect to those typical of folded proteins. While the global free energy minimum consists of highly disordered structures, higher free energy regions correspond to a large variety of transiently structured conformations with secondary structure elements arranged in several different manners, and are not separated from each other by sizeable free energy barriers. From this peculiar structure of the free energy landscape we predict that this peptide should become more structured and not only more compact, with increasing temperatures, and we show that this is the case through a series of biophysical measurements

fMRI single trial discovery of spatio-temporal brain activity patterns

There is growing interest in the description of short-lived patterns in the spatiotemporal cortical activity monitored via neuroimaging. Most traditional analysis methods, designed to estimate relatively long-term brain dynamics, are not always appropriate to capture these patterns. Here we introduce a novel data-driven approach for detecting short-lived fMRI brain activity patterns. Exploiting Density Peak Clustering (Rodriguez and Laio [2014]), our approach reveals well localized clusters by identifying and grouping together voxels whose time-series are similar, irrespective of their brain location, even when very short time windows (∼10 volumes) are used. The method, which we call Coherence Density Peak Clustering (CDPC), is first tested on simulated data and compared with a standard unsupervised approach for fMRI analysis, independent component analysis (ICA). CDPC identifies activated voxels with essentially no false-positives and proves more reliable than ICA, which is troubled by a number of false positives comparable to that of true positives. The reliability of the method is demonstrated on real fMRI data from a simple motor task, containing brief iterations of the same movement. The clusters identified are found in regions expected to be involved in the task, and repeat synchronously with the paradigm. The methodology proposed is especially suitable for the study of short-time brain dynamics and single trial experiments, where the event or task of interest cannot be repeated for the same subject, as happens, for instance, in problem-solving, learning and decision-making. A GUI implementation of our method is available for download at https://github.com/micheleallegra/CDPC. Hum Brain Mapp 38:1421–1437, 2017. © 2016 Wiley Periodicals, Inc

Nucleation process of a fibril precursor in the C-terminal segment of amyloid-beta

By extended atomistic simulations in explicit solvent and bias-exchange metadynamics, we study the aggregation process of 18 chains of the C-terminal segment of amyloid-\u3b2, an intrinsically disordered protein involved in Alzheimer's disease and prone to form fibrils. Starting from a disordered aggregate, we are able to observe the formation of an ordered nucleus rich in beta sheets. The rate limiting step in the nucleation pathway involves crossing a barrier of approximately 40 kcal/mol and is associated with the formation of a very specific interdigitation of the side chains belonging to different sheets. This structural pattern is different from the one observed experimentally in a microcrystal of the same system, indicating that the structure of a "nascent" fibril may differ from the one of an "extended" fibril. \ua9 2013 American Physical Society

The inverted free energy landscape of an intrinsically disordered peptide by simulations and experiments

The free energy landscape theory has been very successful in rationalizing the folding behaviour of globular proteins, as this representation provides intuitive information on the number of states involved in the folding process, their populations and pathways of interconversion. We extend here this formalism to the case of the Aβ40 peptide, a 40-residue intrinsically disordered protein fragment associated with Alzheimer's disease. By using an advanced sampling technique that enables free energy calculations to reach convergence also in the case of highly disordered states of proteins, we provide a precise structural characterization of the free energy landscape of this peptide. We find that such landscape has inverted features with respect to those typical of folded proteins. While the global free energy minimum consists of highly disordered structures, higher free energy regions correspond to a large variety of transiently structured conformations with secondary structure elements arranged in several different manners, and are not separated from each other by sizeable free energy barriers. From this peculiar structure of the free energy landscape we predict that this peptide should become more structured and not only more compact, with increasing temperatures, and we show that this is the case through a series of biophysical measurements

Multidimensional View of Amyloid Fibril Nucleation in Atomistic Detail

Starting from a disordered aggregate, we have simulated the formation of ordered amyloid-like beta structures in a system formed by 18 polyvaline chains in explicit solvent, employing molecular dynamics accelerated by bias-exchange metadynamics. We exploited 8 different collective variables to compute the free energy of hundreds of putative aggregate structures, with variable content of parallel and antiparallel beta-sheets and different packing among the sheets. This allowed characterizing in detail a possible nucleation pathway for the formation of amyloid fibrils: first the system forms a relatively large ordered nucleus of antiparallel beta-sheets, and then a few parallel sheets start appearing. The relevant nucleation process culminates at this point: when a sufficient number of parallel sheets is formed, the free energy starts to decrease toward a new minimum in which this structure is predominant. The complex nucleation pathway we found cannot be described within classical nucleation theory, namely employing a unique simple reaction coordinate like the total content of beta-sheets. \ua9 2012 American Chemical Society