Coupling of Conformational Switches in Calcium Sensor Unraveled with Local Markov Models and Transfer Entropy.

Proteins often have multiple switching domains that are coupled to each other and to the binding of ligands in order to realize signaling functions. Here we investigate the C2A domain of Synaptotagmin-1 (Syt-1), a calcium sensor in the neurotransmitter release machinery and a model system for the large family of C2 membrane binding domains. We combine extensive molecular dynamics (MD) simulations with Markov modeling in order to model conformational switching domains, their states, and their dependence on bound calcium ions.

Complete protein-protein association kinetics in atomic detail revealed by molecular dynamics simulations and Markov modelling.

Protein-protein association is fundamental to many life processes. However, a microscopic model describing the structures and kinetics during association and dissociation is lacking on account of the long lifetimes of associated states, which have prevented efficient sampling by direct molecular dynamics (MD) simulations. Here we demonstrate protein-protein association and dissociation in atomistic resolution for the ribonuclease barnase and its inhibitor barstar by combining adaptive high-throughput MD simulations and hidden Markov modelling.

Spatial Averaging: Sampling Enhancement for Exploring Configurational Space of Atomic Clusters and Biomolecules.

Spatial averaging Monte Carlo (SA-MC) is an efficient algorithm dedicated to the study of rare-event problems. At the heart of this method is the realization that from the equilibrium density a related, modified probability density can be constructed through a suitable transformation. This new density is more highly connected than the original density, which increases the probability for transitions between neighboring states, which in turn speeds up the sampling.

PyEMMA 2: A Software Package for Estimation, Validation, and Analysis of Markov Models.

Markov (state) models (MSMs) and related models of molecular kinetics have recently received a surge of interest as they can systematically reconcile simulation data from either a few long or many short simulations and allow us to analyze the essential metastable structures, thermodynamics, and kinetics of the molecular system under investigation. However, the estimation, validation, and analysis of such models is far from trivial and involves sophisticated and often numerically sensitive methods. In this work we present the open-source Python package PyEMMA ( http://pyemma.

Crystal structure of the dynamin tetramer.

The mechanochemical protein dynamin is the prototype of the dynamin superfamily of large GTPases, which shape and remodel membranes in diverse cellular processes. Dynamin forms predominantly tetramers in the cytosol, which oligomerize at the neck of clathrin-coated vesicles to mediate constriction and subsequent scission of the membrane. Previous studies have described the architecture of dynamin dimers, but the molecular determinants for dynamin assembly and its regulation have remained unclear.

Protein conformational plasticity and complex ligand-binding kinetics explored by atomistic simulations and Markov models.

Understanding the structural mechanisms of protein-ligand binding and their dependence on protein sequence and conformation is of fundamental importance for biomedical research. Here we investigate the interplay of conformational change and ligand-binding kinetics for the serine protease Trypsin and its competitive inhibitor Benzamidine with an extensive set of 150 μs molecular dynamics simulation data, analysed using a Markov state model. Seven metastable conformations with different binding pocket structures are found that interconvert at timescales of tens of microseconds.

The effect of classical and quantum dynamics on vibrational frequency shifts of H2 in clathrate hydrates.

Vibrational frequency shifts of H2 in clathrate hydrates are important to understand the properties and elucidate details of the clathrate structure. Experimental spectra of H2 in clathrate hydrates have been measured for different clathrate compositions, temperatures, and pressures. In order to establish reliable relationships between the clathrate structure, dynamics, and observed frequencies, calculations of vibrational frequency shifts in different clathrate environments are required.

Projected and hidden Markov models for calculating kinetics and metastable states of complex molecules.

Markov state models (MSMs) have been successful in computing metastable states, slow relaxation timescales and associated structural changes, and stationary or kinetic experimental observables of complex molecules from large amounts of molecular dynamics simulation data. However, MSMs approximate the true dynamics by assuming a Markov chain on a clusters discretization of the state space. This approximation is difficult to make for high-dimensional biomolecular systems, and the quality and reproducibility of MSMs has, therefore, been limited.

Overcoming the Rare Event Sampling Problem in Biological Systems with Infinite Swapping.

Infinite swapping (INS) is a recently developed method to address the rare event sampling problem. For INS, an expanded computational ensemble composed of a number of replicas at different temperatures is used, similar to the widely used parallel tempering (PT) method. While the basic concept of PT is to sample various replicas of the system at different temperatures and exchange information between the replicas occasionally, INS uses the symmetrized distribution of configurations in temperature space, which corresponds to the infinite swapping limit of PT.

Rare-event sampling: occupation-based performance measures for parallel tempering and infinite swapping Monte Carlo methods.

In the present paper we identify a rigorous property of a number of tempering-based Monte Carlo sampling methods, including parallel tempering as well as partial and infinite swapping. Based on this property we develop a variety of performance measures for such rare-event sampling methods that are broadly applicable, informative, and straightforward to implement. We illustrate the use of these performance measures with a series of applications involving the equilibrium properties of simple Lennard-Jones clusters, applications for which the performance levels of partial and infinite swapping approaches are found to be higher than those of conventional parallel tempering.

Structure, spectroscopy and dynamics of layered H2O and CO2 ices.

Molecular dynamics simulations of structural, spectroscopic and dynamical properties of mixed water-carbon dioxide (H(2)O-CO(2)) ices are discussed over temperature ranges relevant to atmospheric and astrophysical conditions. The simulations employ multipolar force fields to represent electrostatic interactions which are essential for spectroscopic and dynamical investigations. It is found that at the water/CO(2) interface the water surface acts as a template for the CO(2) component.

Quantifying the importance of protein conformation on ligand migration in myoglobin.

Myoglobin (Mb) is a model system for ligand binding and migration. The energy barriers (ΔG) for ligand migration in Mb have been studied in the past by experiment and theory and significant differences between different approaches were found. From experiment, it is known that Mb can assume a large number of conformational substates.

An infinite swapping approach to the rare-event sampling problem.

We describe a new approach to the rare-event Monte Carlo sampling problem. This technique utilizes a symmetrization strategy to create probability distributions that are more highly connected and, thus, more easily sampled than their original, potentially sparse counterparts. After discussing the formal outline of the approach and devising techniques for its practical implementation, we illustrate the utility of the technique with a series of numerical applications to Lennard-Jones clusters of varying complexity and rare-event character.

Structural and spectroscopic characterization of mixed planetary ices.

Mixed ices play a central role in characterizing the origin, evolution, stability and chemistry of planetary ice surfaces. Examples include the polar areas of Mars, the crust of the Jupiter moon Europa, or atmospheres of planets and their satellites, particularly in the outer solar system. Atomistic simulations using accurate representations of the interaction potentials have recently shown to be suitable to quantitatively describe both, the mid- and the far-infrared spectrum of mixed H2O/CO amorphous ices.

Spatial averaging for small molecule diffusion in condensed phase environments.

Spatial averaging is a new approach for sampling rare-event problems. The approach modifies the importance function which improves the sampling efficiency while keeping a defined relation to the original statistical distribution. In this work, spatial averaging is applied to multidimensional systems for typical problems arising in physical chemistry.

Application of multipolar charge models and molecular dynamics simulations to study stark shifts in inhomogeneous electric fields.

Atomic multipole moments are used to investigate vibrational frequency shifts of CO and H(2) in uniform and inhomogeneous electric fields using ab initio calculations and Molecular Dynamics (MD) simulations. The importance of using atomic multipole moments that can accurately represent both molecular electrostatics and the vibrational response of the molecule to changes in the local electric field is highlighted. The vibrational response of CO to applied uniform and inhomogeneous electric fields is examined using Density Functional Theory calculations for a range of test fields, and the results are used to assess the performance of different atomic multipole models.

Higher order multipole moments for molecular dynamics simulations.

In conventional force fields, the electrostatic potential is represented by atom-centred point charges. This choice is in principle arbitrary, but technically convenient. Point charges can be understood as the first term of multipole expansions, which converge with an increasing number of terms towards the accurate representation of the molecular potential given by the electron density distribution.

Atomistic simulations of CO vibrations in ices relevant to astrochemistry.

The experimental absorption band of carbon monoxide (CO) in mixed ices has been extensively studied in the past. The astrophysical interest in this band is related to its characteristic shape, which appears to depend on the surrounding ice structure. Herein, molecular dynamics simulations are carried out to analyze the relationship between the structure of the ice and the infrared (IR) spectrum of embedded CO molecules at different concentrations.

The role of higher CO-multipole moments in understanding the dynamics of photodissociated carbonmonoxide in myoglobin.

The influence of electrostatic multipole moments up to hexadecapole on the dynamics of photodissociated carbon monoxide (CO) in myoglobin is investigated. The CO electrostatic potential is expressed as an expansion into atomic multipole moments of increasing order up to octopole which are obtained from a distributed multipole analysis. Three models with increasingly accurate molecular multipoles (accurate quadrupole, octopole, and hexadecapole moments, respectively) are developed and used in molecular dynamics simulations.