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 are pigment molecules that protect biomembranes against degradation and may be involved in the formation of functional bacterial membrane microdomains. Little is known on whether different types of carotenoids have different effects on the membrane or if there is any concentration dependence of these effects. In this work, we present results from molecular dynamics simulations of phospholipid bilayers containing different amounts of either β-carotene or zeaxanthin. Both β-carotene and zeaxanthin show the ability to laterally condense the membrane lipids and reduce their inter-leaflet interactions. With increasing concentrations, both carotenoids increase the bilayer thickness and rigidity. The results reveal that carotenoids have similar effects to cholesterol on regulating the behavior of fluid-phase membranes, suggesting that they could function as sterol substitutes and confirming their potential role in the formation of functional membrane domains. [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]

Using temperature replica-exchange molecular dynamics, we characterize a globule-to-coil transition for a softwood-like lignin biopolymer in a tetrahydrofuran (THF)–water cosolvent system at temperatures at which the cosolvent undergoes a de-mixing transition. The lignin is found to be in a coil state, similar to that in the high-temperature miscible region. Analysis of the transition kinetics indicates that THF acts in a surfactant-like fashion. The present study thus suggests that THF–water based pretreatments may efficiently remove lignin from biomass even at relatively low (non-water boiling) temperatures. [ABSTRACT FROM AUTHOR]

The tetrahydrofuran-water binary system exhibits an unusual closed-loop miscibility gap (transitions from a miscible regime to an immiscible regime back to another miscible regime as the temperature increases). Here, using all-atom molecular dynamics simulations, we probe the structural and dynamical behavior of the binary system in the temperature regime of this gap at four different mass ratios, and we compare the behavior of bulk water and tetrahydrofuran. The changes in structure and dynamics observed in the simulations indicate that the temperature region associated with the miscibility gap is distinctive. Within the miscibility-gap temperature region, the self-diffusion of water is significantly altered and the second virial coefficients (pair-interaction strengths) show parabolic-like behavior. Overall, the results suggest that the gap is the result of differing trends with temperature of minor structural changes, which produces interaction virials with parabolic temperature dependence near the miscibility gap. [ABSTRACT FROM AUTHOR]

[Display omitted] • Semi-ordered multivalent IDPs are highly prevalent complexes in the cell. • These difficult targets studied by novel single-particle EM distribution analysis. • Conformational and compositional heterogeneity disentangled and quantified. • Validated and applied to multivalent IDPs bound to the universal hub protein, LC8. • This automated approach may be applied to other IDP systems. Multivalent intrinsically disordered protein (IDP) complexes are prevalent in biology and act in regulation of diverse processes, including transcription, signaling events, and the assembly and disassembly of complex macromolecular architectures. These systems pose significant challenges to structural investigation, due to continuum dynamics imparted by the IDP and compositional heterogeneity resulting from characteristic low-affinity interactions. Here, we developed a modular pipeline for automated single-particle electron microscopy (EM) distribution analysis of common but relatively understudied semi-ordered systems: 'beads-on-a-string' assemblies, composed of IDPs bound at multivalent sites to the ubiquitous ∼20 kDa cross-linking hub protein LC8. This approach quantifies conformational geometries and compositional heterogeneity on a single-particle basis, and statistically corrects spurious observations arising from random proximity of bound and unbound LC8. The statistical correction is generically applicable to oligomer characterization and not specific to our pipeline. Following validation, the approach was applied to the nuclear pore IDP Nup159 and the transcription factor ASCIZ. This analysis unveiled significant compositional and conformational diversity in both systems that could not be obtained from ensemble single particle EM class-averaging strategies, and new insights for exploring how these architectural properties might contribute to their physiological roles in supramolecular assembly and transcriptional regulation. We expect that this approach may be adopted to many other intrinsically disordered systems that have evaded traditional methods of structural characterization. [ABSTRACT FROM AUTHOR]

Ionic liquids have become a popularsolvent for cellulose pretreatmentin biorefineries due to their efficiency in dissolution and theirreusability. Understanding the interactions between cations, anions,and cellulose is key to the development of better solvents and theimprovement of pretreatment conditions. While previous studies describedthe interactions between ionic liquids and cellulose fibers, sheddinglight on the initial stages of the cellulose dissolution process,we study the end state of that process by exploring the structureand dynamics of a single cellulose decamer solvated in 1-butyl-3-methyl-imidazoliumchloride (BmimCl) and in water using replica-exchange molecular dynamics.In both solvents, global structural features of the cellulose chainare similar. However, analyses of local structural properties showthat cellulose explores greater conformational variability in theionic liquid than in water. For instance, in BmimCl the celluloseintramolecular hydrogen bond O3H′···O5 is disruptedmore often resulting in greater flexibility of the solute. Our resultsindicate that the cellulose chain is more dynamic in BmimCl than inwater, which may play a role in the favorable dissolution of cellulosein the ionic liquid. Calculation of the configurational entropy ofthe cellulose decamer confirms its higher conformational flexibilityin BmimCl than in water at elevated temperatures. [ABSTRACT FROM AUTHOR]

[Display omitted] • K ATP channel activity depends on the interplay of its ligands and its subunits. • Inhibitory ligands bias Kir6.2-cytoplasmic domain towards the membrane. • SUR1 cooperates with Kir6.2 to stabilize inhibitor ATP and activator PIP 2 binding. • ATP and PIP 2 compete by having overlapping but non-identical binding residues. • Ligands shift channel conformational dynamics to regulate channel activity. Regulation of pancreatic K ATP channels involves orchestrated interactions of their subunits, Kir6.2 and SUR1, and ligands. Previously we reported K ATP channel cryo-EM structures in the presence and absence of pharmacological inhibitors and ATP, focusing on the mechanisms by which inhibitors act as pharmacological chaperones of K ATP channels (Martin et al., 2019). Here we analyzed the same cryo-EM datasets with a focus on channel conformational dynamics to elucidate structural correlates pertinent to ligand interactions and channel gating. We found pharmacological inhibitors and ATP enrich a channel conformation in which the Kir6.2 cytoplasmic domain is closely associated with the transmembrane domain, while depleting one where the Kir6.2 cytoplasmic domain is extended away into the cytoplasm. This conformational change remodels a network of intra- and inter-subunit interactions as well as the ATP and PIP 2 binding pockets. The structures resolved key contacts between the distal N-terminus of Kir6.2 and SUR1′s ABC module involving residues implicated in channel function and showed a SUR1 residue, K134, participates in PIP 2 binding. Molecular dynamics simulations revealed two Kir6.2 residues, K39 and R54, that mediate both ATP and PIP 2 binding, suggesting a mechanism for competitive gating by ATP and PIP 2. [ABSTRACT FROM AUTHOR]

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., Competing Interests: 10Potentially Conflicting Interests The authors declare the following competing financial interests): L.T.C. is a current member of the Scientific Advisory Board of OpenEye Scientific and an Open Science Fellow with Psivant Sciences. S Zhang, JP Thompson, and DN LeBard are employees of OpenEye Scientific.

Membrane nanodomains have been implicated in Ras signaling, but what these domains are and how they interact with Ras remain obscure. Here, using single particle tracking with photoactivated localization microscopy (spt-PALM) and detailed trajectory analysis, we show that distinct membrane domains dictate KRas G12D (an active KRas mutant) diffusion and trafficking in U2OS cells. KRas G12D exhibits an immobile state in ~70 nm domains, each embedded in a larger domain (~200 nm) that confers intermediate mobility, while the rest of the membrane supports fast diffusion. Moreover, KRas G12D is continuously removed from the membrane via the immobile state and replenished to the fast state, reminiscent of Ras internalization and recycling. Importantly, both the diffusion and trafficking properties of KRas G12D remain invariant over a broad range of protein expression levels. Our results reveal how membrane organization dictates membrane diffusion and trafficking of Ras and offer new insight into the spatial regulation of Ras signaling., Competing Interests: YL, CP, TH, BM, LW, YZ, KT, YC, PS, JG, DZ, XN No competing interests declared, (© 2019, Lee et al.)

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 five tutorials instructing users in the best practices for preparing, carrying out, and analyzing WE simulations for various applications using the WESTPA software. Users are expected to already have significant experience with running standard molecular dynamics simulations using the underlying dynamics engine of interest (e.g. Amber, Gromacs, OpenMM). The tutorials range from a molecular association process in explicit solvent to more complex processes such as host-guest association, peptide conformational sampling, and protein folding., Competing Interests: 9Potentially Conflicting Interests The authors declare no conflicting interests.

Deconstruction of cellulose is crucial for the chemical conversion of lignocellulose into fuel/bioproduct precursors. Recently, a water-organosolv cosolvent system (THF-water) has been shown to both phase-separate on cellulose surfaces and partially deconstruct Avicel  (cellulose) in the absence of acid. Here we employ molecular dynamics simulations to determine whether other common water-organosolv cosolvent systems (acetone, ethanol, and γ-valerolactone) exhibit phase separation at cellulose surface and whether this alters a purely physical cellulose dissociation pathway. Despite finding varied degrees of phase-separation of organosolv on cellulose surfaces, physical dissociation is not enhanced. Interestingly, however, the total amount the median water-cellulose contact lifetimes increases for the cosolvent systems in the order of THF > acetone > ethanol > γ-valerolactone. Together our results indicate two points: a purely physical process for deconstruction of cellulose is unlikely for these cosolvents, and in THF-water, unlike γ-valerolactone- (and some concentrations of acetone and ethanol) water cosolvents, a significant fraction of surface water is slowed. This slowing may be of importance in enhancing chemical deconstruction of cellulose, as it permits an increase in potential THF-water-cellulose reactions, even while the amount of water near cellulose is decreased.

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.

Lipid extracts are an excellent choice of model biomembrane; however at present, there are no commercially available lipid extracts or computational models that mimic microbial membranes containing the branched-chain fatty acids found in many pathogenic and industrially relevant bacteria. We advance the extract of Bacillus subtilis as a standard model for these diverse systems, providing a detailed experimental description and equilibrated atomistic bilayer model included as Supporting Information to this Letter and at ( http://cmb.ornl.gov/members/cheng ). The development and validation of this model represents an advance that enables more realistic simulations and experiments on bacterial membranes and reconstituted bacterial membrane proteins.

Pretreatment facilitates more complete deconstruction of plant biomass to enable more economic production of lignocellulosic biofuels and byproducts. Various co-solvent pretreatments have demonstrated advantages relative to aqueous-only methods by enhancing lignin removal to allow unfettered access to cellulose. However, there is a limited mechanistic understanding of the interactions between the co-solvents and cellulose that impedes further improvement of such pretreatment methods. Recently, tetrahydrofuran (THF) has been identified as a highly effective co-solvent for the pretreatment and fractionation of biomass. To elucidate the mechanism of the THF-water interactions with cellulose, we pair simulation and experimental data demonstrating that enhanced solubilization of cellulose can be achieved by the THF-water co-solvent system at equivolume mixtures and moderate temperatures (≤445 K). The simulations show that THF and water spontaneously phase separate on the local surface of a cellulose fiber, owing to hydrogen bonding of water molecules with the hydrophilic cellulose faces and stacking of THF molecules on the hydrophobic faces. Furthermore, a single fully solvated cellulose chain is shown to be preferentially bound by water molecules in the THF-water mixture. In light of these findings, co-solvent reactions were performed on microcrystalline cellulose and maple wood to show that THF significantly enhanced cellulose deconstruction and lignocellulose solubilization at simulation conditions, enabling a highly versatile and efficient biomass pretreatment and fractionation method.

Cholesterol and ether lipids are ubiquitous in mammalian cell membranes, and their interactions are crucial in ether lipid mediated cholesterol trafficking. We report on cholesterol's molecular interactions with ether lipids as determined using a combination of small-angle neutron and X-ray scattering, and all-atom molecular dynamics (MD) simulations. A scattering density profile model for an ether lipid bilayer was developed using MD simulations, which was then used to simultaneously fit the different experimental scattering data. From analysis of the data the various bilayer structural parameters were obtained. Surface area constrained MD simulations were also performed to reproduce the experimental data. This iterative analysis approach resulted in good agreement between the experimental and simulated form factors. The molecular interactions taking place between cholesterol and ether lipids were then determined from the validated MD simulations. We found that in ether membranes cholesterol primarily hydrogen bonds with the lipid headgroup phosphate oxygen, while in their ester membrane counterparts cholesterol hydrogen bonds with the backbone ester carbonyls. This different mode of interaction between ether lipids and cholesterol induces cholesterol to reside closer to the bilayer surface, dehydrating the headgroup's phosphate moiety. Moreover, the three-dimensional lipid chain spatial density distribution around cholesterol indicates anisotropic chain packing, causing cholesterol to tilt. These insights lend a better understanding of ether lipid-mediated cholesterol trafficking and the roles that the different lipid species have in determining the structural and dynamical properties of membrane associated biomolecules.