Search

The weighted ensemble (WE) strategy has been demonstrated to be highly efficient in generating pathways and rate constants for rare events such as protein folding and protein binding using atomistic molecular dynamics simulations. Here we present two sets of tutorials instructing users in the best practices for preparing, carrying out, and analyzing WE simulations for various applications using the WESTPA software. The first set of more basic tutorials describes a range of simulation types, from a molecular association process in explicit solvent to more complex processes such as host-guest association, peptide conformational sampling, and protein folding. The second set ecompasses six advanced tutorials instructing users in the best practices of using key new features and plugins/extensions of the WESTPA 2.0 software package, which consists of major upgrades for larger systems and/or slower processes. The advanced tutorials demonstrate the use of the following key features: (i) a generalized resampler module for the creation of binless schemes, (ii) a minimal adaptive binning scheme for more efficient surmounting of free energy barriers, (iii) streamlined handling of large simulation datasets using an HDF5 framework, (iv) two different schemes for more efficient rate-constant estimation, (v) a Python API for simplified analysis of WE simulations, and (vi) plugins/extensions for Markovian Weighted Ensemble Milestoning and WE rule-based modeling for systems biology models. Applications of the advanced tutorials include atomistic and non-spatial models, and consist of complex processes such as protein folding and the membrane permeability of a drug-like molecule. Users are expected to already have significant experience with running conventional molecular dynamics or systems biology simulations.

Vascular tone is dependent on smooth muscle KATP channels comprising pore-forming Kir6.1 and regulatory SUR2B subunits, in which mutations cause Cantú syndrome. Unique among KATP isoforms, they lack spontaneous activity and require Mg-nucleotides for activation. Structural mechanisms underlying these properties are unknown. Here, we determined cryogenic electron microscopy structures of vascular KATP channels bound to inhibitory ATP and glibenclamide, which differ informatively from similarly determined pancreatic KATP channel isoform (Kir6.2/SUR1). Unlike SUR1, SUR2B subunits adopt distinct rotational “propeller” and “quatrefoil” geometries surrounding their Kir6.1 core. The glutamate/aspartate-rich linker connecting the two halves of the SUR-ABC core is observed in a quatrefoil-like conformation. Molecular dynamics simulations reveal MgADP-dependent dynamic tripartite interactions between this linker, SUR2B, and Kir6.1. The structures captured implicate a progression of intermediate states between MgADP-free inactivated, and MgADP-bound activated conformations wherein the glutamate/aspartate-rich linker participates as mobile autoinhibitory domain, suggesting a conformational pathway toward KATP channel activation. [ABSTRACT FROM AUTHOR]

Carotenoids have been found to be important in improving the integrity of biomembranes in eukaryotes. However, the molecular details of how carotenoids modulate the physical properties of biomembranes are unknown. To this end, we have conducted a series of molecular dynamics simulations of different biologically-relevant membranes in the presence of carotenoids. The carotenoid effect on the membrane was found to be specific to the identity of the carotenoid and the composition of the membrane itself. Therefore, different classes of carotenoids produce a different effect on the membrane, and different membrane phases are affected differently by carotenoids. It is apparent from our data that carotenoids do trigger the bilayer to become thinner. The mechanism by which this occurs depends on two competing factors, the ability of the lipid tails of opposing monolayers to either (1) compress or (2) interdigitate as the bilayer condenses. Indeed, carotenoids directly influence the physical properties via these two mechanisms, thus compacting the bilayer. However, the degree to which these competing mechanisms are utilized depends on the bilayer phase and the carotenoid identity. [ABSTRACT FROM AUTHOR]

In liquids, the ability of neighboring molecules to rearrange and jostle past each other is directly related to viscosity, the property which describes the propensity to flow. The presence of hydrogen bonds (H-bonds) complicates the molecular scale picture of viscosity. H-Bonds are attractive, directional interactions between molecules which, in some cases, result in transient network structures. In this work, we use experimental and computational methods to demonstrate that the timescale of H-bond network reorganization is the dominant dynamical timescale associated with viscosity for the case of the model H-bonding liquid n-methylacetamide (NMA). This molecule is a peptide analog which forms a transient linear H-bond network. Individual H-bond lifetimes and dynamical fluctuations were observed on the timescale of 1.5 ps, while collective motions and the longest lived population of H-bond partner lifetimes were observed on the order of 20 ps, in agreement with the Maxwell relaxation time. This identifies a mechanism which may aid in understanding the emergence of various complex phenomena arising from transient molecular structures, with implications ranging from the internal dynamics of proteins, to the glass transition, to better understanding the origins of the unique properties of H-bonding liquids. [ABSTRACT FROM AUTHOR]