No abstract

There is increasing evidence that protein dynamics and conformational changes can play an important role in modulating biological function. As a result, experimental and computational methods are being developed, often synergistically, to study the dynamical heterogeneity of a protein or other macromolecules in solution. Thus, methods such as molecular dynamics simulations or ensemble refinement approaches have provided conformational ensembles that can be used to understand protein function and biophysics. These developments have in turn created a need for algorithms and software that can be used to compare structural ensembles in the same way as the root-mean-square-deviation is often used to compare static structures. Although a few such approaches have been proposed, these can be difficult to implement efficiently, hindering a broader applications and further developments. Here, we present an easily accessible software toolkit, called ENCORE, which can be used to compare conformational ensembles generated either from simulations alone or synergistically with experiments. ENCORE implements three previously described methods for ensemble comparison, that each can be used to quantify the similarity between conformational ensembles by estimating the overlap between the probability distributions that underlie them. We demonstrate the kinds of insights that can be obtained by providing examples of three typical use-cases: comparing ensembles generated with different molecular force fields, assessing convergence in molecular simulations, and calculating differences and similarities in structural ensembles refined with various sources of experimental data. We also demonstrate efficient computational scaling for typical analyses, and robustness against both the size and sampling of the ensembles. ENCORE is freely available and extendable, integrates with the established MDAnalysis software package, reads ensemble data in many common formats, and can work with large trajectory files.

We describe a Bayesian/Maximum entropy (BME) procedure and software to construct a conformational ensemble of a biomolecular system by integrating molecular simulations and experimental data. First, an initial conformational ensemble is constructed using for example Molecular Dynamics or Monte Carlo simulations. Due to potential inaccuracies in the model and finite sampling effects, properties predicted from simulations may not agree with experimental data. In BME we use the experimental data to refine the simulation so that the new conformational ensemble has the following properties: (i) the calculated averages are close to the experimental values taking uncertainty into account and (ii) it maximizes the relative Shannon entropy with respect to the original simulation ensemble. The output of this procedure is a set of optimized weights that can be used to calculate arbitrary properties and distributions. Here, we provide a practical guide on how to obtain and use such weights, how to choose adjustable parameters and discuss shortcomings of the method.

Nanodiscs are membrane mimetics that consist of a protein belt surrounding a lipid bilayer, and are broadly used for characterization of membrane proteins. Here, we investigate the structure, dynamics and biophysical properties of two small nanodiscs, MSP1D1ΔH5 and ΔH4H5. We combine our SAXS and SANS experiments with molecular dynamics simulations and previously obtained NMR and EPR data to derive and validate a conformational ensemble that represents the structure and dynamics of the nanodisc. We find that it displays conformational heterogeneity with various elliptical shapes, and with substantial differences in lipid ordering in the centre and rim of the discs. Together, our results reconcile previous apparently conflicting observations about the shape of nanodiscs, and paves the way for future integrative studies of larger complex systems such as membrane proteins embedded in nanodiscs.

The CorA family of proteins regulates the homeostasis of divalent metal ions in many bacteria, archaea, and eukaryotic mitochondria, making it an important target in the investigation of the mechanisms of transport and its functional regulation. Although numerous structures of open and closed channels are now available for the CorA family, the mechanism of the transport regulation remains elusive. Here, we investigated the conformational distribution and associated dynamic behaviour of the pentameric Mg2+ channel CorA at room temperature using small-angle neutron scattering (SANS) in combination with molecular dynamics (MD) simulations and solid-state nuclear magnetic resonance spectroscopy (NMR). We find that neither the Mg2+-bound closed structure nor the Mg2+-free open forms are sufficient to explain the average conformation of CorA. Our data support the presence of conformational equilibria between multiple states, and we further find a variation in the behaviour of the backbone dynamics with and without Mg2+. We propose that CorA must be in a dynamic equilibrium between different non-conducting states, both symmetric and asymmetric, regardless of bound Mg2+ but that conducting states become more populated in Mg2+-free conditions. These properties are regulated by backbone dynamics and are key to understanding the functional regulation of CorA.