Diego Prada-Gracia

The quest for self-consistency in hydrogen bond definitions

In the last decades several hydrogen-bond definitions were proposed by classical computer simulations. Aiming at validating their self-consistency on a wide range of conditions, here we present a comparative study of six among the most common hydrogen-bond definitions for temperatures ranging from 220 K to 400 K and six classical water models. Our results show that, in the interval of temperatures investigated, a generally weak agreement among definitions is present. Moreover, cutoff choice for geometrically based definitions depends on both temperature and water model. As such, analysis of the same water model at different temperatures as well as different water models at the same temperature would require the development of specific cutoff values. Interestingly, large discrepancies were found between two hydrogen-bond definitions which were recently introduced to improve on more conventional methods. Our results reinforce the idea that a more universal way to characterize hydrogen bonds in classical molecular systems is needed.

Towards a microscopic description of the free-energy landscape of water

Free-energy landscape theory is often used to describe complex molecular systems. Here, a microscopic description of water structure and dynamics based on configuration-space-networks and molecular dynamics simulations of the TIP4P/2005 model is applied to investigate the free-energy landscape of water. The latter is built on top of a large set of water microstates describing the kinetic stability of local hydrogen-bond arrangements up to the second solvation shell. In temperature space, the landscape displays three different regimes. At around ambient conditions, the free-energy surface is characterized by many short-lived basins of attraction which are structurally well-defined (inhomogeneous regime). At lower temperatures instead, the liquid rapidly becomes homogeneous. In this regime, the free energy is funneled-like, with fully coordinated water arrangements at the bottom of the funnel. Finally, a third regime develops below the temperature of maximal compressibility (Widom line) where the funnel becomes steeper with few interconversions between microstates other than the fully coordinated ones. Our results present a way to manage the complexity of water structure and dynamics, connecting microscopic properties to its ensemble behavior.

Accounting for the kinetics in order parameter analysis: Lessons from theoretical models and a disordered peptide

Molecular simulations as well as single molecule experiments have been widely analyzed in terms of order parameters, the latter representing candidate probes for the relevant degrees of freedom. Notwithstanding this approach is very intuitive, mounting evidence showed that such descriptions are inaccurate, leading to ambiguous definitions of states and wrong kinetics. To overcome these limitations a framework making use of order parameter fluctuations in conjunction with complex network analysis is investigated. Derived from recent advances in the analysis of single molecule time traces, this approach takes into account the fluctuations around each time point to distinguish between states that have similar values of the order parameter but different dynamics. Snapshots with similar fluctuations are used as nodes of a transition network, the clusterization of which into states provides accurate Markov-state-models of the system under study. Application of the methodology to theoretical models with a noisy order parameter as well as the dynamics of a disordered peptide illustrates the possibility to build accurate descriptions of molecular processes on the sole basis of order parameter time series without using any supplementary information.

Consensus for the Fip35 folding mechanism

Recent advances in computational power and simulation programs finally delivered the first examples of reversible folding for small proteins with an all-atom description. But having at hand the atomistic details of the process did not lead to a straightforward interpretation of the mechanism. For the case of the Fip35 WW-domain where multiple long trajectories of 100 μs are available from D. E. Shaw Research, different interpretations emerged. Some of those are in clear contradiction with each other while others are in qualitative agreement. Here, we present a network-based analysis of the same data by looking at the local fluctuations of conventional order parameters for folding. We found that folding occurs through two major pathways, one almost four times more populated than the other. Each pathway involves the formation of an intermediate with one of the two hairpins in a native configuration. The quantitative agreement of our results with a state-of-the-art reaction coordinate optimization procedure as well as qualitative agreement with other Markov-state-models and different simulation schemes provides strong evidence for a multiple folding pathways scenario with the presence of intermediates.

Structure–dynamics relationship in coherent transport through disordered systems

Quantum transport is strongly influenced by interference with phase relations that depend on the scattering medium. As even small changes in the geometry of the medium can turn constructive interference to destructive, a clear relation between structure and fast, efficient transport is difficult to identify. Here we present a complex network analysis of quantum transport through disordered systems to elucidate the relationship between transport efficiency and structural organization. Evidence is provided for the emergence of structural classes with different geometries but similar high efficiency. Specifically, a structural motif characterized by pair sites, which are not actively participating to the dynamics, renders transport properties robust against perturbations. Our results pave the way for a systematic rationalization of the design principles behind highly efficient transport, which is of paramount importance for technological applications as well as to address transport robustness in natural-light-harvesting complexes.

Water structure-forming capabilities are temperature shifted for different models.

A large number of water models exist for molecular simulations. They differ in the ability to reproduce specific features of real water instead of others, like the correct temperature for the density maximum or the diffusion coefficient. Past analysis mostly concentrated on ensemble quantities, while few data were reported on the different microscopic behavior. Here, we compare seven widely used classical water models (SPC, SPC/E, TIP3P, TIP4P, TIP4P-Ew, TIP4P/2005, and TIP5P) in terms of their local structure-forming capabilities through hydrogen bonds for temperatures ranging from 210 to 350 K by the introduction of a set of order parameters taking into account the configuration of up to the second solvation shell. We found that all models share the same structural pattern up to a temperature shift. When this shift is applied, all models overlap onto a master curve. Interestingly, increased stabilization of fully coordinated structures extending to at least two solvation shells is found for models that are able to reproduce the correct position of the density maximum. Our results provide a self-consistent atomic-level structural comparison protocol, which can be of help in elucidating the influence of different water models on protein structure and dynamics.

A comparative analysis of clustering algorithms: O2 migration in truncated hemoglobin I from transition networks

The ligand migration network for O2-diffusion in truncated Hemoglobin N is analyzed based on three different clustering schemes. For coordinate-based clustering, the conventional k-means and the kinetics-based Markov Clustering (MCL) methods are employed, whereas the locally scaled diffusion map (LSDMap) method is a collective-variable-based approach. It is found that all three methods agree well in their geometrical definition of the most important docking site, and all experimentally known docking sites are recovered by all three methods. Also, for most of the states, their population coincides quite favourably, whereas the kinetics of and between the states differs. One of the major differences between k-means and MCL clustering on the one hand and LSDMap on the other is that the latter finds one large primary cluster containing the Xe1a, IS1, and ENT states. This is related to the fact that the motion within the state occurs on similar time scales, whereas structurally the state is found to be quite diverse. In agreement with previous explicit atomistic simulations, the Xe3 pocket is found to be a highly dynamical site which points to its potential role as a hub in the network. This is also highlighted in the fact that LSDMap cannot identify this state. First passage time distributions from MCL clusterings using a one- (ligand-position) and two-dimensional (ligand-position and protein-structure) descriptor suggest that ligand- and protein-motions are coupled. The benefits and drawbacks of the three methods are discussed in a comparative fashion and highlight that depending on the questions at hand the best-performing method for a particular data set may differ.

Structure and dynamics of water in crowded environments slows down peptide conformational changes.

The concentration of macromolecules inside the cell is high with respect to conventional in vitro experiments or simulations. In an effort to characterize the effects of crowding on the thermodynamics and kinetics of disordered peptides, molecular dynamics simulations were run at different concentrations by varying the number of identical weakly interacting peptides inside the simulation box. We found that the presence of crowding does not influence very much the overall thermodynamics. On the other hand, peptide conformational dynamics was found to be strongly affected, resulting in a dramatic slowing down at larger concentrations. The observation of long lived water bridges between peptides at higher concentrations points to a nontrivial role of the solvent in the altered peptide kinetics. Our results reinforce the idea for an active role of water in molecular crowding, an effect that is expected to be relevant for problems influenced by large solvent exposure areas like in intrinsically disordered proteins.