Multiscale Molecular Dynamics Simulations of an Active Self-Assembling Material.

Inspired by the adaptability observed in biological materials, self-assembly processes have attracted significant interest for their potential to yield novel materials with unique properties. However, experimental methods have often fallen short in capturing the molecular details of the assembly process. In this study, we employ a multiscale molecular dynamics simulation approach, complemented by NMR quantification, to investigate the mechanism of self-assembly in a redox-fueled bioinspired system. Contrary to conventional assumptions, we have uncovered a significant role played by the monomer precursor in the assembly process, with its presence varying with concentration and the extent of conversion of the monomer to the dimer. Experimental confirmation through NMR quantification underscores the concentration-dependent incorporation of monomers into the fibrous structures. Furthermore, our simulations also shed light on the diverse intermolecular interactions, including T-shaped and parallel π stacking, as well as hydrogen bonds, in stabilizing the aggregates. Overall, the open conformation of the dimer is preferred within these aggregates. However, inside the aggregates, the distribution of conformations shifts slightly to the closed conformation compared to on the surface. These findings contribute to the growing field of bioinspired materials science by providing valuable mechanistic and structural insights to guide the design and development of self-assembling materials with biomimetic functionalities.

Native American ataxia medicines rescue ataxia-linked mutant potassium channel activity via binding to the voltage sensing domain.

There are currently no drugs known to rescue the function of Kv1.1 voltage-gated potassium channels carrying loss-of-function sequence variants underlying the inherited movement disorder, Episodic Ataxia 1 (EA1). The Kwakwaka'wakw First Nations of the Pacific Northwest Coast used Fucus gardneri (bladderwrack kelp), Physocarpus capitatus (Pacific ninebark) and Urtica dioica (common nettle) to treat locomotor ataxia. Here, we show that extracts of these plants enhance wild-type Kv1.1 current, especially at subthreshold potentials. Screening of their constituents revealed that gallic acid and tannic acid similarly augment wild-type Kv1.1 current, with submicromolar potency. Crucially, the extracts and their constituents also enhance activity of Kv1.1 channels containing EA1-linked sequence variants. Molecular dynamics simulations reveal that gallic acid augments Kv1.1 activity via a small-molecule binding site in the extracellular S1-S2 linker. Thus, traditional Native American ataxia treatments utilize a molecular mechanistic foundation that can inform small-molecule approaches to therapeutically correcting EA1 and potentially other Kv1.1-linked channelopathies.

Insights into the solubility of γ D-crystallin from multiscale atomistic simulations.

The molecular basis underlying the rich phase behavior of globular proteins remains poorly understood. We use atomistic multiscale molecular simulations to model the solution-state conformational dynamics and interprotein interactions of γ D-crystallin and its P23T-R36S mutant, which drastically limits the protein solubility, at both infinite dilution and at a concentration where the mutant fluid phase and crystalline phase coexist. We find that while the mutant conserves the protein fold, changes to the surface exposure of residues in the neighborhood of residue-36 enhance protein-protein interactions and develop specific protein-protein contacts found in the protein crystal lattice.

Adsorption of 6-MHO on two indoor relevant surface materials: SiO2 and TiO2.

The compound 6-methyl-5-hepten-2-one (6-MHO) is a product of skin oil ozonolysis and is of significance in understanding the role of human occupants in the indoor environment. We present a joint computational and experimental study investigating the adsorption of 6-MHO on two model indoor relevant surfaces, SiO2, a model for a glass window, and TiO2, a component of paint and self-cleaning surfaces. Our classical force field-based molecular dynamics, ab initio molecular dynamics simulations, and FTIR absorption spectra indicate 6-MHO can adsorb on to both of these surfaces via hydrogen and π-hydrogen bonds and is quite stable due to the linear geometry of 6-MHO. Detailed analysis of 6-MHO on the SiO2 surface shows that relative humidity does not impact surface adsorption and adsorbed water does not displace 6-MHO from the hydroxylated SiO2 surface. Additionally, the desorption kinetics of 6-MHO from the hydroxylated SiO2 surface is compared to other compounds found in indoor environments and 6-MHO is shown to desorb with a first order rate constant that is approximately four times slower than that of limonene, but six times faster than that of carvone. In addition, our joint results indicate 6-MHO forms a stronger interaction with the TiO2 surface compared to the SiO2 surface. This study suggests that skin oil ozonolysis products can partition to indoor surfaces leading to the formation of organic films.

Network Hamiltonian Models for Unstructured Protein Aggregates, with Application to γD-Crystallin.

Network Hamiltonian models (NHMs) are a framework for topological coarse-graining of protein-protein interactions, in which each node corresponds to a protein, and edges are drawn between nodes representing proteins that are noncovalently bound. Here, this framework is applied to aggregates of γD-crystallin, a structural protein of the eye lens implicated in cataract disease. The NHMs in this study are generated from atomistic simulations of equilibrium distributions of wild-type and the cataract-causing variant W42R in solution, performed by Wong, E. K.; Prytkova, V.; Freites, J. A.; Butts, C. T.; Tobias, D. J. Molecular Mechanism of Aggregation of the Cataract-Related γD-Crystallin W42R Variant from Multiscale Atomistic Simulations. Biochemistry2019, 58 (35), 3691-3699. Network models are shown to successfully reproduce the aggregate size and structure observed in the atomistic simulation, and provide information about the transient protein-protein interactions therein. The system size is scaled from the original 375 monomers to a system of 10000 monomers, revealing a lowering of the upper tail of the aggregate size distribution of the W42R variant. Extrapolation to higher and lower concentrations is also performed. These results provide an example of the utility of NHMs for coarse-grained simulation of protein systems, as well as their ability to scale to large system sizes and high concentrations, reducing computational costs while retaining topological information about the system.

Interactions of limonene and carvone on titanium dioxide surfaces.

Limonene, a monoterpene, found in cleaning products and air fresheners can interact with a variety of surfaces in indoor environments. An oxidation product of limonene, carvone, has been reported to cause contact allergens. In this study, we have investigated the interactions of limonene and carvone with TiO2, a component of paint and self-cleaning surfaces, at 297 ± 1 K with FTIR spectroscopy and force field-based molecular dynamics and ab initio simulations. The IR absorption spectra and computational methods show that limonene forms π-hydrogen bonds with the surface O-H groups on the TiO2 surface and that carvone adsorbs on the TiO2 surface through a variety of molecular interactions including through carbonyl oxygen atoms with Ti4+ surface atoms, O-H hydrogen bonding (carbonyl O⋯HO) and π-hydrogen bonds with surface O-H groups. Furthermore, we investigated the effects of relative humidity (RH) on the adsorption of limonene and carvone on the TiO2 surface. The spectroscopic results show that the adsorbed limonene can be completely displaced by water at a relative humidity of ca. 50% RH (∼2 MLs of water) and that 25% of carvone is displaced at ca. 67% RH, which agrees with the calculated free energies of adsorption which show carvone more strongly adsorbs on the surface relative to limonene and thus would be harder to displace from the surface. Overall, this study shows how a monoterpene and its oxidation product interact with TiO2 and the impact of relative humidity on these interactions.

Ca2+ -dependent interactions between lipids and the tumor-targeting peptide pHLIP.

Cancerous tissues undergo extensive changes to their cellular environments that differentiate them from healthy tissues. These changes include changes in extracellular pH and Ca2+ concentrations, and the exposure of phosphatidylserine (PS) to the extracellular environment, which can modulate the interaction of peptides and proteins with the plasma membrane. Deciphering the molecular mechanisms of such interactions is critical for advancing the knowledge-based design of cancer-targeting molecular tools, such as pH-low insertion peptide (pHLIP). Here, we explore the effects of PS, Ca2+ , and peptide protonation state on the interactions of pHLIP with lipid membranes. Cellular studies demonstrate that exposed PS on the plasma membrane promotes pHLIP targeting. The magnitude of this effect is dependent on extracellular Ca2+ concentration, indicating that divalent cations play an important role in pHLIP targeting in vivo. The targeting mechanism is further explored with a combination of fluorescence and circular dichroism experiments in model membranes and microsecond-timescale all-atom molecular dynamics simulations. Our results demonstrate that Ca2+ is engaged in coupling peptide-lipid interactions in the unprotonated transmembrane conformation of pHLIP. The simulations reveal that while the pH-induced insertion leads to a strong depletion of PS around pHLIP, the Ca2+ -induced insertion has the opposite effect. Thus, extracellular levels of Ca2+ are crucial to linking cellular changes in membrane lipid composition with the selective targeting and insertion of pHLIP. The characterized Ca2+ -dependent coupling between pHLIP sidechains and PS provides atomistic insights into the general mechanism for lipid-coupled regulation of protein-membrane insertion by divalent cations.

Multiphase Ozonolysis of Oleic Acid-Based Lipids: Quantitation of Major Products and Kinetic Multilayer Modeling.

Commonly found in atmospheric aerosols, cooking oils, and human sebum, unsaturated lipids rapidly decay upon exposure to ozone, following the Criegee mechanism. Here, the gas-surface ozonolysis of three oleic acid-based compounds was studied in a reactor and indoors. Under dry conditions, quantitative product analyses by 1H NMR indicate up to 79% molar yield of stable secondary ozonides (SOZs) in oxidized triolein and methyl oleate coatings. Elevated relative humidity (RH) significantly suppresses the SOZ yields, enhancing the formation of condensed-phase aldehydes and volatile C9 products. Along with kinetic parameters informed by molecular dynamics simulations, these results were used as constraints in a kinetic multilayer model (KM-GAP) simulating triolein ozonolysis. Covering a wide range of coating thicknesses and ozone levels, the model predicts a much faster decay near the gas-lipid interface compared to the bulk. Although the dependence of RH on SOZ yields is well predicted, the model overestimates the production of H2O2 and aldehydes. With negligible dependence on RH, the product composition for oxidized oleic acid is substantially affected by a competitive reaction between Criegee intermediates (CIs) and carboxylic acids. The resulting α-acyloxyalkyl hydroperoxides (α-AAHPs) have much higher molar yields (29-38%) than SOZs (12-16%). Overall, the ozone-lipid chemistry could affect the indoor environment through "crust" accumulation on surfaces and volatile organic compound (VOC) emission. In the atmosphere, the peroxide formation and changes in particle hygroscopicity may have effects on climate. The related health impacts are also discussed.

Effects of Cardiolipin on the Conformational Dynamics of Membrane-Anchored Bcl-xL.

The anti-apoptotic protein Bcl-xL regulates apoptosis by preventing the permeation of the mitochondrial outer membrane by pro-apoptotic pore-forming proteins, which release apoptotic factors into the cytosol that ultimately lead to cell death. Two different membrane-integrated Bcl-xL constructs have been identified: a membrane-anchored and a membrane-inserted conformation. Here, we use molecular dynamics simulations to study the effect of the mitochondrial specific lipid cardiolipin and the protein protonation state on the conformational dynamics of membrane-anchored Bcl-xL. The analysis reveals that the protonation state of the protein and cardiolipin content of the membrane modulate the orientation of the soluble head region (helices α1 through α7) and hence the exposure of its BH3-binding groove, which is required for its interaction with pro-apoptotic proteins.

Anomalous Diffusion of Peripheral Membrane Signaling Proteins from All-Atom Molecular Dynamics Simulations.

Peripheral membrane proteins bind transiently to membrane surfaces as part of many signaling pathways. The bound proteins perform two-dimensional (2-D) diffusion on the membrane surface during the recruitment function. To better understand the interplay between the 2-D diffusion of these protein domains and their membrane binding modes, we performed multimicrosecond all-atom molecular dynamics simulations of two regulatory domains, a C2 domain and a pleckstrin homology (PH) domain, in their experimentally determined bound configuration to a lipid bilayer. The protein bound configurations are preserved throughout the simulation trajectories. Both protein domains exhibit anomalous diffusion with distinct features in their dynamics that reflect their different modes of binding. An analysis of their diffusive behavior reveals common features with the diffusion of lipid molecules in lipid bilayers, suggesting that the 2-D motion of protein domains bound to the membrane surface is modulated by the viscoelastic nature of the lipid bilayer.

HIFs: New arginine mimic inhibitors of the Hv1 channel with improved VSD-ligand interactions.

The human voltage-gated proton channel Hv1 is a drug target for cancer, ischemic stroke, and neuroinflammation. It resides on the plasma membrane and endocytic compartments of a variety of cell types, where it mediates outward proton movement and regulates the activity of NOX enzymes. Its voltage-sensing domain (VSD) contains a gated and proton-selective conduction pathway, which can be blocked by aromatic guanidine derivatives such as 2-guanidinobenzimidazole (2GBI). Mutation of Hv1 residue F150 to alanine (F150A) was previously found to increase 2GBI apparent binding affinity more than two orders of magnitude. Here, we explore the contribution of aromatic interactions between the inhibitor and the channel in the presence and absence of the F150A mutation, using a combination of electrophysiological recordings, classic mutagenesis, and site-specific incorporation of fluorinated phenylalanines via nonsense suppression methodology. Our data suggest that the increase in apparent binding affinity is due to a rearrangement of the binding site allowed by the smaller residue at position 150. We used this information to design new arginine mimics with improved affinity for the nonrearranged binding site of the wild-type channel. The new compounds, named "Hv1 Inhibitor Flexibles" (HIFs), consist of two "prongs," an aminoimidazole ring, and an aromatic group connected by extended flexible linkers. Some HIF compounds display inhibitory properties that are superior to those of 2GBI, thus providing a promising scaffold for further development of high-affinity Hv1 inhibitors.

A novel Hv1 inhibitor reveals a new mechanism of inhibition of a voltage-sensing domain.

Voltage-gated sodium, potassium, and calcium channels consist of four voltage-sensing domains (VSDs) that surround a central pore domain and transition from a down state to an up state in response to membrane depolarization. While many types of drugs bind pore domains, the number of organic molecules known to bind VSDs is limited. The Hv1 voltage-gated proton channel is made of two VSDs and does not contain a pore domain, providing a simplified model for studying how small ligands interact with VSDs. Here, we describe a ligand, named HIF, that interacts with the Hv1 VSD in the up and down states. We find that HIF rapidly inhibits proton conduction in the up state by blocking the open channel, as previously described for 2-guanidinobenzimidazole and its derivatives. HIF, however, interacts with a site slowly accessible in the down state. Functional studies and MD simulations suggest that this interaction traps the compound in a narrow pocket lined with charged residues within the VSD intracellular vestibule, which results in slow recovery from inhibition. Our findings point to a "wrench in gears" mechanism whereby side chains within the binding pocket trap the compound as the teeth of interlocking gears. We propose that the use of screening strategies designed to target binding sites with slow accessibility, similar to the one identified here, could lead to the discovery of new ligands capable of interacting with VSDs of other voltage-gated ion channels in the down state.

Heterogeneous Interactions of Prevalent Indoor Oxygenated Organic Compounds on Hydroxylated SiO2 Surfaces.

Oxygenated organic compounds (OOCs) are widely found in indoor environments and come from either the direct emissions from indoor activities or the subsequent oxidation of nonoxygenated OCs. Adsorption and partitioning of OCs on surfaces are significant processes in indoor chemistry, yet these interactions specifically involving OOCs are still poorly understood. In this study, we investigate the interactions of three prevalent indoor OOCs (dihydromyrcenol, α-terpineol, and linalool) on an indoor surface proxy (hydroxylated SiO2) by combining vibrational spectroscopy with ab initio molecular dynamics simulations. The adsorption of these compounds on the SiO2 surface is driven by π hydrogen bonding and O-H hydrogen bonding interactions, with O-H hydrogen bonding interactions being stronger. The results of kinetic measurements suggest that indoor surfaces play a significant role in the removal of these OOCs, especially under moderate and low air exchange. Additionally, indoor surfaces can also serve as a reservoir of OOCs due to their much slower desorption kinetics when compared to other indoor relevant organic compounds such as limonene. Overall, the results gleaned by experiment and theoretical simulations provide a molecular representation of the interaction of OOCs on indoor relevant surfaces as well as implications of these interactions for indoor air chemistry.

Adsorption of constitutional isomers of cyclic monoterpenes on hydroxylated silica surfaces.

We present a study of four monoterpene isomers (limonene, γ-terpinene, terpinolene, and α-pinene) that are prevalent in indoor environments and their interaction with the hydroxylated SiO2 surface, a model for the glass surface, by combining infrared spectroscopy and computational simulations. These isomers are molecularly adsorbed onto SiO2 through π-hydrogen bonds with surface hydroxyl groups. However, experimental results suggest that the strength of interaction of these compounds with the SiO2 surface varies for each isomer, with α-pinene showing the weakest interaction. This observation is supported by molecular dynamics simulations that α-pinene adsorbed on the SiO2 surface has lower free energy of desorption and a lower mass accommodation coefficient compared to other isomers. Additionally, our ab initio molecular dynamics simulations show lower π-hydrogen bonding probabilities for α-pinene compared to the other three constitutional isomers. Importantly, these interactions are most likely present for a range of other systems involving organic compounds and solid surfaces and, thus, provide a thorough framework for comparing the interactions of organic molecules on indoor relevant surfaces.

Molecular Orientation at the Squalene/Air Interface from Sum Frequency Generation Spectroscopy and Atomistic Modeling.

Human skin oils are significant scavengers of atmospheric oxidants in occupied indoor environments, and squalene is a major ozone-active constituent. Here, we present a combined spectroscopic and atomistic modeling approach to elucidate the conformational and orientational preferences of squalene at the air/oil interface and their implications for reactions with ozone. We find that squalene chains have a tendency to align with the surface normal, resulting in different concentrations of the various types of its double bonds and thus different reactivities. We also observe the presence of water at the surface of this hydrophobic compound. Both findings have possible implications for the design and outcomes of kinetic models describing this important aspect of indoor air chemistry.

Diffusivelike Motions in a Solvent-Free Protein-Polymer Hybrid.

The interaction between proteins and hydration water stabilizes protein structure and promotes functional dynamics, with water translational motions enabling protein flexibility. Engineered solvent-free protein-polymer hybrids have been shown to preserve protein structure, function, and dynamics. Here, we used neutron scattering, protein and polymer perdeuteration, and molecular dynamics simulations to explore how a polymer dynamically replaces water. Even though relaxation rates and vibrational properties are strongly modified in polymer coated compared to hydrated proteins, liquidlike polymer dynamics appear to plasticize the conjugated protein in a qualitatively similar way as do hydration-water translational motions.

Thermodynamics and Mechanism of the Membrane Permeation of Hv1 Channel Blockers.

The voltage-gated proton channel Hv1 mediates efflux of protons from the cell. Hv1 integrally contributes to various physiological processes including pH homeostasis and the respiratory burst of phagocytes. Inhibition of Hv1 may provide therapeutic avenues for the treatment of inflammatory diseases, breast cancer, and ischemic brain damage. In this work, we investigate two prototypical Hv1 inhibitors, 2-guanidinobenzimidazole (2GBI), and 5-chloro-2-guanidinobenzimidazole (GBIC), from an experimentally screened class of guanidine derivatives. Both compounds block proton conduction by binding the same site located on the intracellular side of the channel. However, when added to the extracellular medium, the compounds strongly differ in their ability to inhibit proton conduction, suggesting substantial differences in membrane permeability. Here, we compute the potential of mean force for each compound to permeate through the membrane using atomistic molecular dynamics simulations with the adaptive biasing force method. Our results rationalize the putative distinction between these two blockers with respect to their abilities to permeate the cellular membrane.

Impact of Adsorbed Water on the Interaction of Limonene with Hydroxylated SiO2: Implications of π-Hydrogen Bonding for Surfaces in Humid Environments.

The indoor environment is a dynamic one with many variables impacting indoor air quality and indoor air chemistry. These include relative humidity (RH) and the presence of different surfaces. Although it has been suggested that the indoor concentrations of gas-phase compounds increase at higher relative humidity, because of displacement of these compounds from indoor surfaces, little is known from a molecular perspective about how RH and adsorbed water impact the adsorption of indoor relevant organic compounds such as limonene with indoor relevant surfaces. Herein, we investigate the effects of RH on the adsorption of limonene, a hydrophobic molecule, on hydroxylated SiO2 surfaces, a model for glass surfaces. Experimental data using infrared spectroscopy to directly measure limonene adsorption are combined with both force field-based molecular dynamics (MD) and ab initio molecular dynamics (AIMD) simulations to understand the competitive interactions between limonene, water, and the SiO2 surface. The spectroscopic data provide evidence that adsorbed limonene is not completely displaced by adsorbed water, even at high RH (∼80%) when the water layer coverage is close to three monolayers (MLs). These experimental data are supported by AIMD and MD simulations, which indicate that limonene is present at the adsorbed water interface but displaced from direct interactions with SiO2. This study shows that although some limonene can desorb from the surface, even at the highest RH, more than half the limonene remains adsorbed on the surface that can undergo continued surface reactivity. A complex network of π-hydrogen bonds, water-water hydrogen bonds, and SiO2-water hydrogen bonds explains these interactions at the air/adsorbed water/SiO2 interface that hold the hydrophobic limonene molecule at the interface. Importantly, these interactions are most likely present for a range of other systems involving organic compounds and solid surfaces at ambient relative humidity and may be important in a range of scientific areas, from sensor development to cultural heritage science.

Insights on small molecule binding to the Hv1 proton channel from free energy calculations with molecular dynamics simulations.

Hv1 is a voltage-gated proton channel whose main function is to facilitate extrusion of protons from the cell. The development of effective channel blockers for Hv1 can lead to new therapeutics for the treatment of maladies related to Hv1 dysfunction. Although the mechanism of proton permeation in Hv1 remains to be elucidated, a series of small molecules have been discovered to inhibit Hv1. Here, we computed relative binding free energies of a prototypical Hv1 blocker on a model of human Hv1 in an open state. We used alchemical free energy perturbation techniques based on atomistic molecular dynamics simulations. The results support our proposed open state model and shed light on the preferred tautomeric state of the channel blocker. This work lays the groundwork for future studies on adapting the blocker molecule for more effective inhibition of Hv1.

Voltage-dependent structural models of the human Hv1 proton channel from long-timescale molecular dynamics simulations.

The voltage-gated Hv1 proton channel is a ubiquitous membrane protein that has roles in a variety of cellular processes, including proton extrusion, pH regulation, production of reactive oxygen species, proliferation of cancer cells, and increased brain damage during ischemic stroke. A crystal structure of an Hv1 construct in a putative closed state has been reported, and structural models for the channel open state have been proposed, but a complete characterization of the Hv1 conformational dynamics under an applied membrane potential has been elusive. We report structural models of the Hv1 voltage-sensing domain (VSD), both in a hyperpolarized state and a depolarized state resulting from voltage-dependent conformational changes during a 10-µs-timescale atomistic molecular dynamics simulation in an explicit membrane environment. In response to a depolarizing membrane potential, the S4 helix undergoes an outward displacement, leading to changes in the VSD internal salt-bridge network, resulting in a reshaping of the permeation pathway and a significant increase in hydrogen bond connectivity throughout the channel. The total gating charge displacement associated with this transition is consistent with experimental estimates. Molecular docking calculations confirm the proposed mechanism for the inhibitory action of 2-guanidinobenzimidazole (2GBI) derived from electrophysiological measurements and mutagenesis. The depolarized structural model is also consistent with the formation of a metal bridge between residues located in the core of the VSD. Taken together, our results suggest that these structural models are representative of the closed and open states of the Hv1 channel.