Exploring Molecular Energy Landscapes by Coupling the DFTB Potential with a Tree-Based Stochastic Algorithm: Investigation of the Conformational Diversity of Phthalates.

Exploring the global energy landscape of relatively large molecules at the quantum level is a challenging problem. In this work, we report the coupling of a nonredundant conformational space exploration method, namely, the robotics-inspired iterative global exploration and local optimization (IGLOO) algorithm, with the quantum-chemical density functional tight binding (DFTB) potential. The application of this fast and efficient computational approach to three close-sized molecules of the phthalate family (DBP, BBP, and DEHP) showed that they present different conformational landscapes. These differences have been rationalized by making use of descriptors based on distances and dihedral angles. Coulomb interactions, steric hindrance, and dispersive interactions have been found to drive the geometric properties. A strong correlation has been evidenced between the two dihedral angles describing the side-chain orientation of the phthalate molecules. Our approach identifies low-energy minima without prior knowledge of the potential energy surface, paving the way for future investigations into transition paths and states.

WASCO: A Wasserstein-based Statistical Tool to Compare Conformational Ensembles of Intrinsically Disordered Proteins.

The structural investigation of intrinsically disordered proteins (IDPs) requires ensemble models describing the diversity of the conformational states of the molecule. Due to their probabilistic nature, there is a need for new paradigms that understand and treat IDPs from a purely statistical point of view, considering their conformational ensembles as well-defined probability distributions. In this work, we define a conformational ensemble as an ordered set of probability distributions and provide a suitable metric to detect differences between two given ensembles at the residue level, both locally and globally. The underlying geometry of the conformational space is properly integrated, one ensemble being characterized by a set of probability distributions supported on the three-dimensional Euclidean space (for global-scale comparisons) and on the two-dimensional flat torus (for local-scale comparisons). The inherent uncertainty of the data is also taken into account to provide finer estimations of the differences between ensembles. Additionally, an overall distance between ensembles is defined from the differences at the residue level. We illustrate the potential of the approach with several examples of applications for the comparison of conformational ensembles: (i) produced from molecular dynamics (MD) simulations using different force fields, and (ii) before and after refinement with experimental data. We also show the usefulness of the method to assess the convergence of MD simulations, and discuss other potential applications such as in machine-learning-based approaches. The numerical tool has been implemented in Python through easy-to-use Jupyter Notebooks available at https://gitlab.laas.fr/moma/WASCO.

MoMA-LoopSampler: a web server to exhaustively sample protein loop conformations.

SUMMARY: MoMA-LoopSampler is a sampling method that globally explores the conformational space of flexible protein loops. It combines a large structural library of three-residue fragments and a novel reinforcement-learning-based approach to accelerate the sampling process while maintaining diversity. The method generates a set of statistically likely loop states satisfying geometric constraints, and its ability to sample experimentally observed conformations has been demonstrated. This paper presents a web user interface to MoMA-LoopSampler through the illustration of a typical use-case. AVAILABILITY AND IMPLEMENTATION: MoMA-LoopSampler is freely available at: https://moma.laas.fr/applications/LoopSampler/. We recommend users to create an account, but anonymous access is possible. In most cases, jobs are completed within a few minutes. The waiting time may increase depending on the server load, but it very rarely exceeds an hour. For users requiring more intensive use, binaries can be provided upon request. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Predicting Secondary Structure Propensities in IDPs Using Simple Statistics from Three-Residue Fragments.

Intrinsically disordered proteins (IDPs) play key functional roles facilitated by their inherent plasticity. In most of the cases, IDPs recognize their partners through partially structured elements inserted in fully disordered chains. The identification and characterization of these elements is fundamental to understand the functional mechanisms of IDPs. Although several computational methods have been developed to identify order within disordered chains, most of the current secondary structure predictors are focused on globular proteins and are not necessarily appropriate for IDPs. Here, we present a comprehensible method, called Local Structural Propensity Predictor (LS2P), to predict secondary structure elements from IDP sequences. LS2P performs statistical analyses from a database of three-residue fragments extracted from coil regions of high-resolution protein structures. In addition to identifying scarcely populated helical and extended regions, the method pinpoints short stretches triggering ß-turn formation or promoting α-helices. The simplicity of the method enables a direct connection between experimental observations and structural features encoded in IDP sequences.

Collective information processing in human phase separation.

In our digital societies, individuals massively interact through digital interfaces whose impact on collective dynamics can be important. In particular, the combination of social media filters and recommender systems can lead to the emergence of polarized and fragmented groups. In some social contexts, such segregation processes of human groups have been shown to share similarities with phase separation phenomena in physics. Here, we study the impact of information filtering on collective segregation behaviour of human groups. We report a series of experiments where groups of 22 subjects have to perform a collective segregation task that mimics the tendency of individuals to bond with other similar individuals. More precisely, the participants are each assigned a colour (red or blue) unknown to them, and have to regroup with other subjects sharing the same colour. To assist them, they are equipped with an artificial sensory device capable of detecting the majority colour in their 'environment' (defined as their k nearest neighbours, unbeknownst to them), for which we control the perception range, k = 1, 3, 5, 7, 9, 11, 13. We study the separation dynamics (emergence of unicolour groups) and the properties of the final state, and show that the value of k controls the quality of the segregation, although the subjects are totally unaware of the precise definition of the 'environment'. We also find that there is a perception range k = 7 above which the ability of the group to segregate does not improve. We introduce a model that precisely describes the random motion of a group of pedestrians in a confined space, and which faithfully reproduces and allows interpretation of the results of the segregation experiments. Finally, we discuss the strong and precise analogy between our experiment and the phase separation of two immiscible materials at very low temperature. This article is part of the theme issue 'Multi-scale analysis and modelling of collective migration in biological systems'.

MoMA-LigPath: a web server to simulate protein-ligand unbinding.

Protein-ligand interactions taking place far away from the active site, during ligand binding or release, may determine molecular specificity and activity. However, obtaining information about these interactions with experimental or computational methods remains difficult. The computational tool presented in this article, MoMA-LigPath, is based on a mechanistic representation of the molecular system, considering partial flexibility, and on the application of a robotics-inspired algorithm to explore the conformational space. Such a purely geometric approach, together with the efficiency of the exploration algorithm, enables the simulation of ligand unbinding within short computing time. Ligand unbinding pathways generated by MoMA-LigPath are a first approximation that can provide useful information about protein-ligand interactions. When needed, this approximation can be subsequently refined and analyzed using state-of-the-art energy models and molecular modeling methods. MoMA-LigPath is available at http://moma.laas.fr. The web server is free and open to all users, with no login requirement.