Nanoconfinement of Carbon Dioxide within Interfacial Aqueous/Ionic Liquid Systems.

Nanoporous, gas-selective membranes have shown encouraging results for the removal of CO2 from flue gas, yet the optimal design for such membranes is often unknown. Therefore, we used molecular dynamics simulations to elucidate the behavior of CO2 within aqueous and ionic liquid (IL) systems ([EMIM][TFSI] and [OMIM][TFSI]), both confined individually and as an interfacial aqueous/IL system. We found that within aqueous systems the mobility of CO2 is reduced due to interactions between the CO2 oxygens and hydroxyl groups on the pore surface. Within the IL systems, we found that confinement has a greater effect on the [EMIM][TFSI] system as opposed to the [OMIM][TFSI] system. Paradoxically, the larger and more asymmetrical [OMIM]+ molecule undergoes less efficient packing, resulting in fewer confinement effects. Free energy surfaces of the nanoconfined aqueous/IL interface demonstrate that CO2 will transfer spontaneously from the aqueous to the IL phase.

Energetics of high temperature degradation of fentanyl into primary and secondary products.

Fentanyl is a synthetic opioid used for managing chronic pain. Due to its higher potency (50-100×) than morphine, fentanyl is also an abused drug. A sensor that could detect illicit fentanyl by identifying its thermally degraded fragments would be helpful to law enforcement. While experimental studies have probed the thermal degradation of fentanyl, little theoretical work has been done to understand the mechanism. Here, we studied the thermal degradation pathways of fentanyl using extensive ab initio molecular dynamics simulations combined with enhanced sampling via multiple-walker metadynamics. We calculated the free energy profile for each bond suggested earlier as a potential degradation point to map the thermodynamic driving forces. We also estimated the forward attempt rate of each bond degradation reaction to gain information about degradation kinetics.

Mechanical Activation of MscL Revealed by a Locally Distributed Tension Molecular Dynamics Approach.

Membrane tension perceived by mechanosensitive (MS) proteins mediates cellular responses to mechanical stimuli and osmotic stresses, and it also guides multiple biological functions including cardiovascular control and development. In bacteria, MS channels function as tension-activated pores limiting excessive turgor pressure, with MS channel of large conductance (MscL) acting as an emergency release valve preventing cell lysis. Previous attempts to simulate gating transitions in MscL by either directly applying steering forces to the protein or by increasing the whole-system tension were not fully successful and often disrupted the integrity of the system. We present a novel, to our knowledge, locally distributed tension molecular dynamics (LDT-MD) simulation method that allows application of forces continuously distributed among lipids surrounding the channel using a specially constructed collective variable. We report reproducible and reversible transitions of MscL to the open state with measured parameters of lateral expansion and conductivity that exactly satisfy experimental values. The LDT-MD method enables exploration of the MscL-gating process with different pulling velocities and variable tension asymmetry between the inner and outer membrane leaflets. We use LDT-MD in combination with well-tempered metadynamics to reconstruct the tension-dependent free-energy landscape for the opening transition in MscL. The flexible definition of the LDT collective variable allows general application of our method to study mechanical activation of any membrane-embedded protein.

Hydration Mimicry by Membrane Ion Channels.

Ions transiting biomembranes might pass readily from water through ion-specific membrane proteins if these protein channels provide environments similar to the aqueous solution hydration environment. Indeed, bulk aqueous solution is an important reference condition for the ion permeation process. Assessment of this hydration mimicry concept depends on understanding the hydration structure and free energies of metal ions in water in order to provide a comparison for the membrane channel environment. To refine these considerations, we review local hydration structures of ions in bulk water and the molecular quasi-chemical theory that provides hydration free energies. In doing so, we note some current views of ion binding to membrane channels and suggest new physical chemical calculations and experiments that might further clarify the hydration mimicry concept.

Combined molecular/continuum modeling reveals the role of friction during fast unfolding of coiled-coil proteins.

Coiled-coils are filamentous proteins that form the basic building block of important force-bearing cellular elements, such as intermediate filaments and myosin motors. In addition to their biological importance, coiled-coil proteins are increasingly used in new biomaterials including fibers, nanotubes, or hydrogels. Coiled-coils undergo a structural transition from an α-helical coil to an unfolded state upon extension, which allows them to sustain large strains and is critical for their biological function. By performing equilibrium and out-of-equilibrium all-atom molecular dynamics (MD) simulations of coiled-coils in explicit solvent, we show that two-state models based on Kramers' or Bell's theories fail to predict the rate of unfolding at high pulling rates. We further show that an atomistically informed continuum rod model accounting for phase transformations and for the hydrodynamic interactions with the solvent can reconcile two-state models with our MD results. Our results show that frictional forces, usually neglected in theories of fibrous protein unfolding, reduce the thermodynamic force acting on the interface, and thus control the dynamics of unfolding at different pulling rates. Our results may help interpret MD simulations at high pulling rates, and could be pertinent to cytoskeletal networks or protein-based artificial materials subjected to shocks or blasts.

Insertion of Dengue E into lipid bilayers studied by neutron reflectivity and molecular dynamics simulations.

The envelope (E) protein of Dengue virus rearranges to a trimeric hairpin to mediate fusion of the viral and target membranes, which is essential for infectivity. Insertion of E into the target membrane serves to anchor E and possibly also to disrupt local order within the membrane. Both aspects are likely to be affected by the depth of insertion, orientation of the trimer with respect to the membrane normal, and the interactions that form between trimer and membrane. In the present work, we resolved the depth of insertion, the tilt angle, and the fundamental interactions for the soluble portion of Dengue E trimers (sE) associated with planar lipid bilayer membranes of various combinations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-rac-glycerol (POPG), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and cholesterol (CHOL) by neutron reflectivity (NR) and by molecular dynamics (MD) simulations. The results show that the tip of E containing the fusion loop (FL) is located at the interface of the headgroups and acyl chains of the outer leaflet of the lipid bilayers, in good agreement with prior predictions. The results also indicate that E tilts with respect to the membrane normal upon insertion, promoted by either the anionic lipid POPG or CHOL. The simulations show that tilting of the protein correlates with hydrogen bond formation between lysines and arginines located on the sides of the trimer close to the tip (K246, K247, and R73) and nearby lipid headgroups. These hydrogen bonds provide a major contribution to the membrane anchoring and may help to destabilize the target membrane.

Examining the Mechanical Equilibrium of Microscopic Stresses in Molecular Simulations.

The microscopic stress field provides a unique connection between atomistic simulations and mechanics at the nanoscale. However, its definition remains ambiguous. Rather than a mere theoretical preoccupation, we show that this fact acutely manifests itself in local stress calculations of defective graphene, lipid bilayers, and fibrous proteins. We find that popular definitions of the microscopic stress violate the continuum statements of mechanical equilibrium, and we propose an unambiguous and physically sound definition.

Structural Rearrangement of the AT1 Receptor Modulated by Membrane Thickness and Tension.

Membrane-embedded mechanosensitive (MS) proteins, including ion channels and G-protein coupled receptors (GPCRs), are essential for the transduction of external mechanical stimuli into biological signals. The angiotensin II type 1 (AT1) receptor plays many important roles in cardiovascular regulation and is associated with diseases such as hypertension and congestive heart failure. The membrane-mediated activation of the AT1 receptor is not well understood, despite this being one of the most widely studied GPCRs within the context of biased agonism. Here, we use extensive molecular dynamics (MD) simulations to characterize the effect of the local membrane environment on the activation of the AT1 receptor. We show that membrane thickness plays an important role in the stability of active and inactive states of the receptor, as well as the dynamic interchange between states. Furthermore, our simulation results show that membrane tension is effective in driving large-scale structural changes in the inactive state such as the outward movement of transmembrane helix 6 to stabilize intermediate active-like conformations. We conclude by comparing our simulation observations with AlphaFold 2 predictions, as a proxy to experimental structures, to provide a framework for how membrane mediated stimuli can facilitate activation of the AT1 receptor through the ß-arrestin signaling pathway.

Metadynamics simulations reveal mechanisms of Na+ and Ca2+ transport in two open states of the channelrhodopsin chimera, C1C2.

Cation conducting channelrhodopsins (ChRs) are a popular tool used in optogenetics to control the activity of excitable cells and tissues using light. ChRs with altered ion selectivity are in high demand for use in different cell types and for other specialized applications. However, a detailed mechanism of ion permeation in ChRs is not fully resolved. Here, we use complementary experimental and computational methods to uncover the mechanisms of cation transport and valence selectivity through the channelrhodopsin chimera, C1C2, in the high- and low-conducting open states. Electrophysiology measurements identified a single-residue substitution within the central gate, N297D, that increased Ca2+ permeability vs. Na+ by nearly two-fold at peak current, but less so at stationary current. We then developed molecular models of dimeric wild-type C1C2 and N297D mutant channels in both open states and calculated the PMF profiles for Na+ and Ca2+ permeation through each protein using well-tempered/multiple-walker metadynamics. Results of these studies agree well with experimental measurements and demonstrate that the pore entrance on the extracellular side differs from original predictions and is actually located in a gap between helices I and II. Cation transport occurs via a relay mechanism where cations are passed between flexible carboxylate sidechains lining the full length of the pore by sidechain swinging, like a monkey swinging on vines. In the mutant channel, residue D297 enhances Ca2+ permeability by mediating the handoff between the central and cytosolic binding sites via direct coordination and sidechain swinging. We also found that altered cation binding affinities at both the extracellular entrance and central binding sites underly the distinct transport properties of the low-conducting open state. This work significantly advances our understanding of ion selectivity and permeation in cation channelrhodopsins and provides the insights needed for successful development of new ion-selective optogenetic tools.

Erratum: Addition to "Carbonic Anhydrase Robustness for Use in Nanoscale CO2 Capture Technologies".

[This corrects the article DOI: 10.1021/acsomega.3c02630.].

Membrane mediated mechanical stimuli produces distinct active-like states in the AT1 receptor.

The Angiotensin II Type 1 (AT1) receptor is one of the most widely studied GPCRs within the context of biased signaling. While the AT1 receptor is activated by agonists such as the peptide AngII, it can also be activated by mechanical stimuli such as membrane stretch or shear in the absence of a ligand. Despite the importance of mechanical activation of the AT1 receptor in biological processes such as vasoconstriction, little is known about the structural changes induced by external physical stimuli mediated by the surrounding lipid membrane. Here, we present a systematic simulation study that characterizes the activation of the AT1 receptor under various membrane environments and mechanical stimuli. We show that stability of the active state is highly sensitive to membrane thickness and tension. Structural comparison of membrane-mediated vs. agonist-induced activation shows that the AT1 receptor has distinct active conformations. This is supported by multi-microsecond free energy calculations that show unique landscapes for the inactive and various active states. Our modeling results provide structural insights into the mechanical activation of the AT1 receptor and how it may produce different functional outcomes within the framework of biased agonism.

Tight hydrophobic core and flexible helices yield MscL with a high tension gating threshold and a membrane area mechanical strain buffer.

The mechanosensitive (MS) channel of large conductance, MscL, is the high-tension threshold osmolyte release valve that limits turgor pressure in bacterial cells in the event of drastic hypoosmotic shock. Despite MscL from Mycobacterium tuberculosis (TbMscL) being the first structurally characterized MS channel, its protective mechanism of activation at nearly-lytic tensions has not been fully understood. Here, we describe atomistic simulations of expansion and opening of wild-type (WT) TbMscL in comparison with five of its gain-of-function (GOF) mutants. We show that under far-field membrane tension applied to the edge of the periodic simulation cell, WT TbMscL expands into a funnel-like structure with trans-membrane helices bent by nearly 70°, but does not break its 'hydrophobic seal' within extended 20 µs simulations. GOF mutants carrying hydrophilic substitutions in the hydrophobic gate of increasing severity (A20N, V21A, V21N, V21T and V21D) also quickly transition into funnel-shaped conformations but subsequently fully open within 1-8 µs. This shows that solvation of the de-wetted (vapor-locked) constriction is the rate-limiting step in the gating of TbMscL preceded by area-buffering silent expansion. Pre-solvated gates in these GOF mutants reduce this transition barrier according to hydrophilicity and the most severe V21D eliminates it. We predict that the asymmetric shape-change of the periplasmic side of the channel during the silent expansion provides strain-buffering to the outer leaflet thus re-distributing the tension to the inner leaflet, where the gate resides.

Effect of Local Stress on Accurate Modeling of Bacterial Outer Membranes Using All-Atom Molecular Dynamics.

Biological membranes are fundamental components of living organisms that play an undeniable role in their survival. Molecular dynamics (MD) serves as an essential computational tool for studying biomembranes on molecular and atomistic scales. The status quo of MD simulations of biomembranes studies a nanometer-sized membrane patch periodically extended under periodic boundary conditions (PBCs). In nature, membranes are usually composed of different lipids in their two layers (referred to as leaflets). This compositional asymmetry imposes a fixed ratio of lipid numbers between the two leaflets in a periodically constrained membrane, which needs to be set appropriately. The widely adopted methods of defining a leaflet lipid ratio suffer from the lack of control over the mechanical tension of each leaflet, which could significantly influence research findings. In this study, we investigate the role of membrane-building protocol and the resulting initial stress state on the interaction between small molecules and asymmetric membranes. We model the outer membrane of Pseudomonas aeruginosa bacteria using two different building protocols and probe their interactions with the Pseudomonas quinolone signal (PQS). Our results show that differential stress could shift the position of free energy minimum for the PQS molecule between the two leaflets of the asymmetric membrane. This work provides critical insights into the relationship between the initial per-leaflet tension and the spontaneous intercalation of PQS.

Carbonic Anhydrase Robustness for Use in Nanoscale CO2 Capture Technologies.

Continued dependence on crude oil and natural gas resources for fossil fuels has caused global atmospheric carbon dioxide (CO2) emissions to increase to record-setting proportions. There is an urgent need for efficient and inexpensive carbon sequestration systems to mitigate large-scale emissions of CO2 from industrial flue gas. Carbonic anhydrase (CA) has shown high potential for enhanced CO2 capture applications compared to conventional absorption-based methods currently utilized in various industrial settings. This study aims to understand structural aspects that contribute to the stability of CA enzymes critical for their applications in industrial processes, which require the ability to withstand conditions different from those in their native environments. Here, we evaluated the thermostability and enzyme activity of mesophilic and thermophilic CA variants at different temperature conditions and in the presence of atmospheric gas pollutants like nitrogen oxides and sulfur oxides. Based on our enzyme activity assays and molecular dynamics simulations, we see increased conformational stability and CA activity levels in thermostable CA variants incubated week-long at different temperature conditions. The thermostable CA variants also retained high levels of CA activity despite changes in solution pH due to increasing NO and SO2 concentrations. A loss of CA activity was observed only at high concentrations of NO/SO2 that possibly can be minimized with the appropriate buffered solutions.

Identification, Heterologous Expression, and Characterization of the Tolypodiol Biosynthetic Gene Cluster through an Integrated Approach.

Cyanobacteria are tremendous producers of biologically active natural products, including the potent anti-inflammatory compound tolypodiol. However, linking biosynthetic gene clusters with compound production in cyanobacteria has lagged behind that in other bacterial genera. Tolypodiol is a meroterpenoid originally isolated from the cyanobacterium HT-58-2. Here we describe the identification of the tolypodiol biosynthetic gene cluster through heterologous expression in Anabaena and in vitro protein assays of a methyltransferase found in the tolypodiol biosynthetic gene cluster. We have also identified similar biosynthetic gene clusters in cyanobacterial and actinobacterial genomes, suggesting that meroterpenoids with structural similarity to the tolypodiols may be synthesized by other microbes. We also report the identification of two new analogs of tolypodiol that we have identified in both the original and heterologous producer. This work further illustrates the usefulness of Anabaena as a heterologous expression host for cyanobacterial compounds and how integrated approaches can help to link natural product compounds with their producing biosynthetic gene clusters.

Transport and Energetics of Carbon Dioxide in Ionic Liquids at Aqueous Interfaces.

A major hurdle in utilizing carbon dioxide (CO2) lies in separating it from industrial flue gas mixtures and finding suitable storage methods that enable its application in various industries. To address this issue, we utilized a combination of molecular dynamics simulations and experiments to investigate the behavior of CO2 in common room-temperature ionic liquids (RTIL) when in contact with aqueous interfaces. Our investigation of RTILs, [EMIM][TFSI] and [OMIM][TFSI], and their interaction with a pure water layer mimics the environment of a previously developed ultrathin enzymatic liquid membrane for CO2 separation. We analyzed diffusion constants and viscosity, which reveals that CO2 molecules exhibit faster mobility within the selected ILs compared to what would be predicted solely based on the viscosity of the liquids using the standard Einstein-Stokes relation. Moreover, we calculated the free energy of translocation for various species across the aqueous-IL interface, including CO2 and HCO3-. Free energy profiles demonstrate that CO2 exhibits a more favorable partitioning behavior in the RTILs compared to that in pure water, while a significant barrier hinders the movement of HCO3- from the aqueous layer. Experimental measurement of the CO2 transport in the RTILs corroborates the model. These findings strongly suggest that hydrophobic RTILs could serve as a promising option for selectively transporting CO2 from aqueous media and concentrating it as a preliminary step toward storage.

Role of unsaturated lipid and ergosterol in ethanol tolerance of model yeast biomembranes.

We present a combined atomic force microscopy and fluorescence microscopy study of the behavior of a ternary supported lipid bilayer system containing a saturated lipid (DPPC), an unsaturated lipid (DOPC), and ergosterol in the presence of high ethanol (20 vol %). We find that the fluorescent probe Texas Red DHPE preferentially partitions into the ethanol-induced interdigitated phase, which allows the use of fluorescence imaging to investigate the phase behavior of the system. Atomic force microscopy and fluorescence images of samples with the same lipid mixture show good agreement in sample morphology and area fractions of the observed phases. Using area fractions obtained from fluorescence images over a broad range of compositions, we constructed a phase diagram of the DPPC/DOPC/ergosterol system at 20 vol % ethanol. The phase diagram clearly shows that increasing unsaturated lipid and/or ergosterol protects the membrane by preventing the formation of the interdigitated phase. This result supports the hypothesis that yeast cells increase ergosterol and unsaturated lipid content to prevent interdigitation and maintain an optimal membrane thickness as ethanol concentration increases during anaerobic fermentations. Changes in plasma membrane composition provide an important survival factor for yeast cells to deter ethanol toxicity.

Spinodal de-wetting of light liquids on graphene.

We demonstrate theoretically the possibility of spinodal de-wetting in heterostructures made of light-atom liquids (hydrogen, helium, and nitrogen) deposited on suspended graphene. Extending our theory of film growth on two-dimensional (2D) materials to include analysis of surface instabilities via the hydrodynamic Cahn-Hilliard-type equation, we characterize in detail the spatial and temporal scales of the resulting spinodal de-wetting patterns. Both linear stability analysis and direct numerical simulations of the surface hydrodynamics show micron-sized (generally material dependent) patterns of 'dry' regions. The physical reason for the development of such instabilities on graphene can be traced back to the inherently weak van der Waals interactions between atomically thin materials and atoms in the liquid. Thus 2D materials could represent a new theoretical and technological platform for studies of spinodal de-wetting.

Crystalline, ordered and disordered lipid membranes: convergence of stress profiles due to ergosterol.

We present a simulation study focusing on modulations of the stress, or lateral pressure, profiles of lipid bilayer phases by addition of a sterol, ergosterol, at multiple temperatures. A major redistribution of lateral and normal pressures across the gel-phase bilayer required 10 mol % sterol in comparison to a gradual redistribution beginning at 20 mol % for the liquid phase. Stress profiles across all temperatures converged at 30 mol % ergosterol. Redistribution and merging of stress profiles, associated with structural alterations, are coincident with experimentally observed modulations in mechanical properties and therefore are suggested as the mechanism of action for this biologically necessary role of sterols.

AFM for structure and dynamics of biomembranes.

We review structure and dynamic measurements of biomembranes by atomic force microscopy (AFM). We focus mainly on studies involving supported lipid bilayers (SLBs), particularly formation by vesicle rupture on flat and corrugated surfaces, nucleation and growth of domains in phase-separated systems, anesthetic-lipid interactions, and protein/peptide interactions in multicomponent systems. We show that carefully designed experiments along with real-time AFM imaging with superior lateral and z resolution (0.1 nm) have revealed quantitative details of the mechanisms and factors controlling vesicle rupture, domain shape and size, phase transformations, and some model biological interactions. The AFM tip can also be used as a mechanical transducer and incorporated in electrochemical measurements of membrane components; therefore, we touch on these important applications in both model and cell membranes.