News
Extensive conformational heterogeneity within protein cores
Basic principles of statistical mechanics require that proteins sample an ensemble of conformations at any nonzero temperature. However, it is still common to treat the crystallographic structure of a protein as the structure of its native state, largely because high-resolution structural characterization of protein flexibility remains a profound challenge. To assess the typical degree of conformational heterogeneity within folded proteins, we construct Markov state models describing …
Read moreCloud-based simulations on Google Exacycle reveal ligand modulation of GPCR activation pathways
Simulations can provide tremendous insight into the atomistic details of biological mechanisms, but micro- to millisecond timescales are historically only accessible on dedicated supercomputers. We demonstrate that cloud computing is a viable alternative that brings long-timescale processes within reach of a broader community. We used Google’s Exacycle cloud-computing platform to simulate two milliseconds of dynamics of a major drug target, the G-protein-coupled receptor β2AR. Ma…
Read moreMarkov State Models to Elucidate Ligand Binding Mechanism
Molecular dynamics simulations can now routinely access the microsecond timescale, making feasible direct sampling of ligand association events. While Markov State Model (MSM) approaches offer a useful framework for analyzing such trajectory data to gain insight into binding mechanisms, accurate modeling of ligand association pathways and kinetics must be done…
Read moreDiscovery of multiple hidden allosteric sites by combining Markov state models and experiments
The discovery of drug-like molecules that bind pockets in proteins that are not present in crystallographic structures yet exert allosteric control over activity has generated great interest in designing pharmaceuticals that exploit allosteric effects. However, there have only been a small number of successes, so the therapeutic potential of these pockets–called hidden allosteric sites–remains unclear. One challenge for assessing their utility is that rational drug design approaches…
Read moreFluctuations within folded proteins: implications for thermodynamic and allosteric regulation
Folded protein structures are both stable and dynamic. Historically, our clearest window into these structures came from X-ray crystallography, which generally provided a static image of each protein’s singular “folded state”, highlighting its stability. Deviations away from that crystallographic structure were difficult to quantify, and as a result, their potential functional consequences were often neglected. However, several dynamical and statistical studies now highlight the struc…
Read moreAccurately modeling nanosecond protein dynamics requires at least microseconds of simulation
Advances in hardware and algorithms have greatly extended the timescales accessible to molecular simulation. This article assesses whether such long timescale simulations improve our ability to calculate order parameters that describe conformational heterogeneity on ps-ns timescales or if force fields are now a limiting factor. Order parameters from experiment are compared with order parameters calculated in three different ways from simulations ranging from 10 ns to 100 μs in leng…
Read moreChoice of Adaptive Sampling Strategy Impacts State Discovery, Transition Probabilities, and the Apparent Mechanism of Conformational Changes
Interest in atomically detailed simulations has grown significantly with recent advances in computational hardware and Markov state modeling (MSM) methods, yet outstanding questions remain that hinder their widespread adoption. Namely, how do alternative sampling strategies explore conformational space and how might this influence predictions generated from the data? Here, we seek to answer these questions for four commonly used sampling methods: (1) a single long simulation, (2) many…
Read more