Residue Folding Degree—Relationship to Secondary Structure Categories and Use as Collective Variable

Recently, we have shown that the residue folding degree, a network-based measure of folded content in proteins, is able to capture backbone conformational transitions related to the formation of secondary structures in molecular dynamics (MD) simulations. In this work, we focus primarily on developing a collective variable (CV) for MD based on this residue-bound parameter to be able to trace the evolution of secondary structure in segments of the protein. We show that this CV can do just that and that the related energy profiles (potentials of mean force, PMF) and transition barriers are comparable to those found by others for particular events in the folding process of the model mini protein Trp-cage. Hence, we conclude that the relative segment folding degree (the newly proposed CV) is a computationally viable option to gain insight into the formation of secondary structures in protein dynamics. We also show that this CV can be directly used as a measure of the amount of α-helical content in a selected segment.

Structure of benzothiadiazine at zwitterionic phospholipid cell membranes

The use of drugs derived from benzothiadiazine, which is a bicyclic heterocyclic benzene derivative, has become a widespread treatment for diseases such as hypertension (treated with diuretics such as bendroflumethiazide or chlorothiazide), low blood sugar (treated with non-diuretic diazoxide) or the human immunodeficiency virus, among others. In this work we have investigated the interactions of benzothiadiazine with the basic components of cell membranes and solvents such as phospholipids, cholesterol, ions and water. The analysis of the mutual microscopic interactions is of central importance to elucidate the local structure of benzothiadiazine as well as the mechanisms responsible for the access of benzothiadiazine to the interior of the cell. We have performed molecular dynamics simulations of benzothiadiazine embedded in three different model zwitterionic bilayer membranes made by dimyristoilphosphatidylcholine, dioleoylphosphatidylcholine, 1,2- dioleoyl-sn-glycero-3-phosphoserine and cholesterol inside aqueous sodium-chloride solution in order to systematically examine microscopic interactions of benzothiadiazine with the cell membrane at liquid-crystalline phase conditions. From data obtained through radial distribution functions, hydrogen-bonding lengths and potentials of mean force based on reversible work calculations, we have observed that benzothiadiazine has a strong affinity to stay at the cell membrane interface although it can be fully solvated by water in short periods of time. Furthermore, benzothiadiazine is able to bind lipids and cholesterol chains by means of single and double hydrogen-bonds of different characteristic lengths.

Multi-Scale Modeling of Aqueous-Phase Methane Diffusion in Silicate Channels

Silica-based materials including zeolites are commonly used for wide ranging applications including separations and catalysis.<br>Substrate transport rates in these materials often significantly influence the efficiency of such applications.<br>Two factors that contribute to transport rates include<br>1) the porosity of the silicate matrix and<br>2) non-bonding interactions between the diffusing species and the silicate surface.<br>Here, we utilize computer simulation to resolve the relative contribution of these factors to effective methane transport rates in a silicate channel.<br>Specifically, we develop a `homogenized' model of methane transport valid at micron and longer length scales that incorporates atomistic-scale kinetic information.<br>The atomistic-scale data are obtained from extensive molecular dynamics simulations that yield local diffusion coefficients and potentials of mean force.<br>With this model, we demonstrate how nuances in silicate hydration and silica/methane interactions impact 'macroscale' methane diffusion rates in bulk silicate materials.<br>This hybrid homogenization/molecular dynamics approach will be of general use for describing small molecule transport in materials with detailed molecular interactions.<br><br>

Multi-Scale Modeling of Aqueous-Phase Methane Diffusion in Silicate Channels

Silica-based materials including zeolites are commonly used for wide ranging applications including separations and catalysis.<br>Substrate transport rates in these materials often significantly influence the efficiency of such applications.<br>Two factors that contribute to transport rates include<br>1) the porosity of the silicate matrix and<br>2) non-bonding interactions between the diffusing species and the silicate surface.<br>Here, we utilize computer simulation to resolve the relative contribution of these factors to effective methane transport rates in a silicate channel.<br>Specifically, we develop a `homogenized' model of methane transport valid at micron and longer length scales that incorporates atomistic-scale kinetic information.<br>The atomistic-scale data are obtained from extensive molecular dynamics simulations that yield local diffusion coefficients and potentials of mean force.<br>With this model, we demonstrate how nuances in silicate hydration and silica/methane interactions impact 'macroscale' methane diffusion rates in bulk silicate materials.<br>This hybrid homogenization/molecular dynamics approach will be of general use for describing small molecule transport in materials with detailed molecular interactions.<br><br>

Metal Cation-Binding Mechanisms of Q-Proline Peptoid Macrocycles in Solution

<div>The rational design of foldable and functionalizable peptidomimetic scaffolds requires the concerted application of both computational and experimental methods. Recently, a new class of designed peptoid macrocycle incorporating spiroligomer proline mimics (Q-prolines) has been found to pre-organize when bound by monovalent metal cations. To determine the solution-state structure of these cation-bound macrocycles, we employ a Bayesian inference method (BICePs) to reconcile enhanced-sampling molecular simulations with sparse ROESY correlations from experimental NMR studies. The BICePs approach circumvents the need for bespoke force field parameterization, instead relying on experimental restraints to help narrow the possible set of <i>cis</i>/<i>trans</i> amide isomers in solution. Conformations predicted to be most populated in solution were then simulated in the presence of explicit cations to yield trajectories with observed binding events, revealing a highly-preorganized all-<i>trans</i> amide conformation, whose formation is likely limited by the slow rate of <i>cis</i>/<i>trans</i> isomerization. Interestingly, this conformation differs from a racemic crystal structure solved in the absence of cation. Free energies of cation binding computed from distance-dependent potentials of mean force suggest Na⁺ has higher affinity to the macrocycle than K⁺, with both cations binding much more strongly in acetonitrile than water. The simulated affinities are able to correctly rank the extent to which different macrocycle sequences exhibit preorganization in the presence of different metal cations and solvents, suggesting our approach is suitable for solution-state computational design.</div>

Metal Cation-Binding Mechanisms of Q-Proline Peptoid Macrocycles in Solution

<div>The rational design of foldable and functionalizable peptidomimetic scaffolds requires the concerted application of both computational and experimental methods. Recently, a new class of designed peptoid macrocycle incorporating spiroligomer proline mimics (Q-prolines) has been found to pre-organize when bound by monovalent metal cations. To determine the solution-state structure of these cation-bound macrocycles, we employ a Bayesian inference method (BICePs) to reconcile enhanced-sampling molecular simulations with sparse ROESY correlations from experimental NMR studies. The BICePs approach circumvents the need for bespoke force field parameterization, instead relying on experimental restraints to help narrow the possible set of <i>cis</i>/<i>trans</i> amide isomers in solution. Conformations predicted to be most populated in solution were then simulated in the presence of explicit cations to yield trajectories with observed binding events, revealing a highly-preorganized all-<i>trans</i> amide conformation, whose formation is likely limited by the slow rate of <i>cis</i>/<i>trans</i> isomerization. Interestingly, this conformation differs from a racemic crystal structure solved in the absence of cation. Free energies of cation binding computed from distance-dependent potentials of mean force suggest Na⁺ has higher affinity to the macrocycle than K⁺, with both cations binding much more strongly in acetonitrile than water. The simulated affinities are able to correctly rank the extent to which different macrocycle sequences exhibit preorganization in the presence of different metal cations and solvents, suggesting our approach is suitable for solution-state computational design.</div>

Ion-mediated interactions between like-charged polyelectrolytes with bending flexibility

Ion-mediated interactions between polyelectrolytes (PEs) are crucial to the properties of flexible biopolymers such as nucleic acids and proteins but the effect of PE flexibility on such interactions has not been explicitly addressed until now. In this work, the potentials of mean force (PMFs) between like-charged PEs with different bending flexibility have been investigated by Monte Carlo simulations and a cylindrical confinement around each PE was involved to model two PEs in an array. We found that in the absence of trivalent salt, the PMFs between like-charged PEs in an array are apparently repulsive while the bending flexibility can visibly decrease the repulsive PMFs. With the addition of high trivalent salt, the PMFs become significantly attractive whereas the attractive PMFs can be apparently weakened by the bending flexibility. Our analyses reveal that the effect of bending flexibility is attributed to the increased PE conformational space, which allows the PEs to fluctuate away to decrease the monovalent ion-mediated repulsion or to weaken the trivalent ion-mediated attraction through disrupting trivalent ion-bridging configuration. Additionally, our further calculations show that the effect of bending flexibility on the ion-mediated interactions is less apparent for PEs without cylindrical confinement.

Bulk Ionic Screening Lengths from Extremely Large-Scale Molecular Dynamics Simulations

Recent experiments have reported anomalously large screening lengths of interactions between charged surfaces confining concentrated electrolytes and ionic liquids. Termed underscreening, this effect was ascribed to bulk properties of dense ionic systems. Herein, we study bulk ionic screening with extremely large-scale molecular dynamics simulations, allowing us to assess the range of distances relevant to the experiments. Our results yield two screening lengths satisfying distinct scaling relations. However, with an accuracy of 10^-5 kT in interionic potentials of mean force, we find no signs of underscreening, suggesting that other than bulk effects might be at play in the experiments.

Bulk Ionic Screening Lengths from Extremely Large-Scale Molecular Dynamics Simulations

Recent experiments have reported anomalously large screening lengths of interactions between charged surfaces confining concentrated electrolytes and ionic liquids. Termed underscreening, this effect was ascribed to bulk properties of dense ionic systems. Herein, we study bulk ionic screening with extremely large-scale molecular dynamics simulations, allowing us to assess the range of distances relevant to the experiments. Our results yield two screening lengths satisfying distinct scaling relations. However, with an accuracy of 10^-5 kT in interionic potentials of mean force, we find no signs of underscreening, suggesting that other than bulk effects might be at play in the experiments.