Carlos X. Hernández

MDTraj: A Modern Open Library for the Analysis of Molecular Dynamics Trajectories

As molecular dynamics (MD) simulations continue to evolve into powerful computational tools for studying complex biomolecular systems, the necessity of flexible and easy-to-use software tools for the analysis of these simulations is growing. We have developed MDTraj, a modern, lightweight, and fast software package for analyzing MD simulations. MDTraj reads and writes trajectory data in a wide variety of commonly used formats. It provides a large number of trajectory analysis capabilities including minimal root-mean-square-deviation calculations, secondary structure assignment, and the extraction of common order parameters. The package has a strong focus on interoperability with the wider scientific Python ecosystem, bridging the gap between MD data and the rapidly growing collection of industry-standard statistical analysis and visualization tools in Python. MDTraj is a powerful and user-friendly software package that simplifies the analysis of MD data and connects these datasets with the modern interactive data science software ecosystem in Python.

Markov State Models Provide Insights into Dynamic Modulation of Protein Function

Conspectus Protein function is inextricably linked to protein dynamics. As we move from a static structural picture to a dynamic ensemble view of protein structure and function, novel computational paradigms are required for observing and understanding conformational dynamics of proteins and its functional implications. In principle, molecular dynamics simulations can provide the time evolution of atomistic models of proteins, but the long time scales associated with functional dynamics make it difficult to observe rare dynamical transitions. The issue of extracting essential functional components of protein dynamics from noisy simulation data presents another set of challenges in obtaining an unbiased understanding of protein motions. Therefore, a methodology that provides a statistical framework for efficient sampling and a human-readable view of the key aspects of functional dynamics from data analysis is required. The Markov state model (MSM), which has recently become popular worldwide for studying protein dynamics, is an example of such a framework. In this Account, we review the use of Markov state models for efficient sampling of the hierarchy of time scales associated with protein dynamics, automatic identification of key conformational states, and the degrees of freedom associated with slow dynamical processes. Applications of MSMs for studying long time scale phenomena such as activation mechanisms of cellular signaling proteins has yielded novel insights into protein function. In particular, from MSMs built using large-scale simulations of GPCRs and kinases, we have shown that complex conformational changes in proteins can be described in terms of structural changes in key structural motifs or “molecular switches” within the protein, the transitions between functionally active and inactive states of proteins proceed via multiple pathways, and ligand or substrate binding modulates the flux through these pathways. Finally, MSMs also provide a theoretical toolbox for studying the effect of nonequilibrium perturbations on conformational dynamics. Considering that protein dynamics in vivo occur under nonequilibrium conditions, MSMs coupled with nonequilibrium statistical mechanics provide a way to connect cellular components to their functional environments. Nonequilibrium perturbations of protein folding MSMs reveal the presence of dynamically frozen glass-like states in their conformational landscape. These frozen states are also observed to be rich in β-sheets, which indicates their possible role in the nucleation of β-sheet rich aggregates such as those observed in amyloid-fibril formation. Finally, we describe how MSMs have been used to understand the dynamical behavior of intrinsically disordered proteins such as amyloid-β, human islet amyloid polypeptide, and p53. While certainly not a panacea for studying functional dynamics, MSMs provide a rigorous theoretical foundation for understanding complex entropically dominated processes and a convenient lens for viewing protein motions.

MSMBuilder: Statistical Models for Biomolecular Dynamics

MSMBuilder is a software package for building statistical models of high-dimensional time-series data. It is designed with a particular focus on the analysis of atomistic simulations of biomolecular dynamics such as protein folding and conformational change. MSMBuilder is named for its ability to construct Markov state models (MSMs), a class of models that has gained favor among computational biophysicists. In addition to both well-established and newer MSM methods, the package includes complementary algorithms for understanding time-series data such as hidden Markov models and time-structure based independent component analysis. MSMBuilder boasts an easy to use command-line interface, as well as clear and consistent abstractions through its Python application programming interface. MSMBuilder was developed with careful consideration for compatibility with the broader machine learning community by following the design of scikit-learn. The package is used primarily by practitioners of molecular dynamics, but is just as applicable to other computational or experimental time-series measurements.

Osprey: Hyperparameter Optimization for Machine Learning

Osprey is a tool for hyperparameter optimization of machine learning algorithms in Python. Hyperparameter optimization can often be an onerous process for researchers, due to timeconsuming experimental replicates, non-convex objective functions, and constant tension between exploration of global parameter space and local optimization (Jones, Schonlau, and Welch 1998). We’ve designed Osprey to provide scientists with a practical, easyto-use way of finding optimal model parameters. The software works seamlessly with scikit-learn estimators (Pedregosa et al. 2011) and supports many different search strategies for choosing the next set of parameters with which to evaluate a given model, including gaussian processes (GPy 2012), tree-structured Parzen estimators (Yamins, Tax, and Bergstra 2013), as well as random and grid search. As hyperparameter optimization is an embarrassingly parallel problem, Osprey can easily scale to hundreds of concurrent processes by executing a simple command-line program multiple times. This makes it easy to exploit large resources available in high-performance computing environments.

Structure‐based network analysis of an evolved G protein‐coupled receptor homodimer interface

Crystallographic structures and experimental assays of human CXC chemokine receptor type 4 (CXCR4) provide strong evidence for the capacity to homodimerize, potentially as a means of allosteric regulation. Even so, how this homodimer forms and its biological significance has yet to be fully characterized. By applying principles from network analysis, sequence‐based approaches such as statistical coupling analysis to determine coevolutionary residues, can be used in conjunction with molecular dynamics simulations to identify residues relevant to dimerization. Here, the predominant coevolution sector lies along the observed dimer interface, suggesting functional relevance. Furthermore, coevolution scoring provides a basis for determining significant nodes, termed hubs, in the network formed by residues found along the interface of the homodimer. These node residues coincide with hotspots indicating potential druggability. Drug design efforts targeting such key residues could potentially result in modulation of binding and therapeutic benefits for disease states, such as lung cancers, lymphomas and latent HIV‐1 infection. Furthermore, this method may be applied to any protein–protein interaction.

MSMExplorer: Data Visualizations for Biomolecular Dynamics

MSMExplorer is a Python package for visualizing data generated from biomolecular dynamics. While molecular visualizations have been a large focus of the molecular dynamics (MD) community (Humphrey, Dalke, and Schulten 1996, Schrödinger, LLC (2015)), data visualizations for the analyses of MD trajectories have been less developed. MSMExplorer seeks to fill this niche by providing publication-quality statistical plots with an easy-to-use Python API that works seamlessly with commonly used Python libraries, such as numpy and scikit-learn (Walt, Colbert, and Varoquaux 2011, Pedregosa et al. (2011)). Additionally, plots are generated using already established plotting libraries, like seaborn, to provide a consistent aesthetic (Waskom et al. 2016, Hunter (2007), Hagberg, Schult, and Swart (2008), Foreman-Mackey (2016)).

MDEntropy: Information-Theoretic Analyses for Molecular Dynamics

MDEntropy is a Python package for information-theoretic (IT) analysis of data generated from molecular dynamics simulations. While correlation studies have long been of interest to the molecular dynamics (MD) community (McCammon, Gelin, and Karplus 1977, McClendon et al. (2009)), IT tools to analyze MD trajectories have been much less developed. MDEntropy seeks to fill this niche by providing an easy-to-use Python API that works seamlessly with other Python packages, such as mdtraj, msmbuilder, and numpy (R. T. McGibbon et al. 2015, Walt, Colbert, and Varoquaux (2011), M. P. Harrigan et al. (2017)).

Understanding the Origins of a Pandemic Virus

Understanding the origin of infectious diseases provides scientifically based rationales for implementing public health measures that may help to avoid or mitigate future epidemics. The recent ancestors of a pandemic virus provide invaluable information about the set of minimal genomic alterations that transformed a zoonotic agent into a full human pandemic. Since the first confirmed cases of the H1N1 pandemic virus in the spring of 2009, several hypotheses about the strains origins have been proposed. However, how, where, and when it first infected humans is still far from clear. The only way to piece together this epidemiological puzzle relies on the collective effort of the international scientific community to increase genomic sequencing of influenza isolates, especially ones collected in the months prior to the origin of the pandemic.

Modelling Intrinsically Disordered Protein Dynamics as Networks of Transient Secondary Structure

Describing the dynamics and conformational landscapes of Intrinsically Disordered Proteins (IDPs) is of paramount importance to understanding their functions. Markov State Models (MSMs) are often used to characterize the dynamics of more structured proteins, but models of IDPs built using conventional MSM modelling protocols can be difficult to interpret due to the inherent nature of IDPs, which exhibit fast transitions between disordered microstates. We propose a new method of determining MSM states from all-atom molecular dynamics simulation data of IDPs by using per-residue secondary structure assignments as input features in a MSM model. Because such secondary structure algorithms use a select set of features for assignment (dihedral angles, contact distances, etc.), they represent a knowledge-based refinement of feature sets used for model-building. This method adds interpretability to IDP conformational landscapes, which are increasingly viewed as composed of transient secondary structure, and allows us to readily use MSM analysis tools in this paradigm. We demonstrate the use of our method with the transcription factor p53 c-terminal domain (p53-CTD), a commonly-studied IDP. We are able to characterize the full secondary structure phase space observed for p53-CTD, and describe characteristics of p53-CTD as a network of transient helical and beta-hairpin structures with different network behaviors in different domains of secondary structure. This analysis provides a novel example of how IDPs can be studied and how researchers might better understand a disordered protein conformational landscape.