Lipid-mediated hydrophobic gating in the BK potassium channel.

The large-conductance, calcium-activated potassium (BK) channel lacks the typical intracellular bundle-crossing gate present in most ion channels of the 6TM family. This observation, initially inferred from Ca$^{2+}$-free-pore accessibility experiments and recently corroborated by a CryoEM structure of the non-conductive state, raises a puzzling question: how can gating occur in absence of steric hindrance? To answer this question, we carried out molecular simulations and accurate free energy calculations to obtain a microscopic picture of the sequence of events that, starting from a Ca$^{2+}$-free state leads to ion conduction upon Ca$^{2+}$ binding. Our results highlight an unexpected role for annular lipids, which turn out to be an integral part of the gating machinery. Due to the presence of fenestrations, the "closed" Ca$^{2+}$-free pore can be occupied by the methyl groups from the lipid alkyl chains. This dynamic occupancy triggers and stabilizes the nucleation of a vapor bubble into the inner pore cavity, thus hindering ion conduction. By contrast, Ca$^{2+}$ binding results into a displacement of these lipids outside the inner cavity, lowering the hydrophobicity of this region and thus allowing for pore hydration and conduction. This lipid-mediated hydrophobic gating rationalizes several seemingly problematic experimental observations, including the state-dependent pore accessibility of blockers.

Hydrophobically gated memristive nanopores for neuromorphic applications.

Signal transmission in the brain relies on voltage-gated ion channels, which exhibit the electrical behaviour of memristors, resistors with memory. State-of-the-art technologies currently employ semiconductor-based neuromorphic approaches, which have already demonstrated their efficacy in machine learning systems. However, these approaches still cannot match performance achieved by biological neurons in terms of energy efficiency and size. In this study, we utilise molecular dynamics simulations, continuum models, and electrophysiological experiments to propose and realise a bioinspired hydrophobically gated memristive nanopore. Our findings indicate that hydrophobic gating enables memory through an electrowetting mechanism, and we establish simple design rules accordingly. Through the engineering of a biological nanopore, we successfully replicate the characteristic hysteresis cycles of a memristor and construct a synaptic device capable of learning and forgetting. This advancement offers a promising pathway for the realization of nanoscale, cost- and energy-effective, and adaptable bioinspired memristors.

Charting Nanocluster Structures via Convolutional Neural Networks.

A general method to obtain a representation of the structural landscape of nanoparticles in terms of a limited number of variables is proposed. The method is applied to a large data set of parallel tempering molecular dynamics simulations of gold clusters of 90 and 147 atoms, silver clusters of 147 atoms, and copper clusters of 147 atoms, covering a plethora of structures and temperatures. The method leverages convolutional neural networks to learn the radial distribution functions of the nanoclusters and distills a low-dimensional chart of the structural landscape. This strategy is found to give rise to a physically meaningful and differentiable mapping of the atom positions to a low-dimensional manifold in which the main structural motifs are clearly discriminated and meaningfully ordered. Furthermore, unsupervised clustering on the low-dimensional data proved effective at further splitting the motifs into structural subfamilies characterized by very fine and physically relevant differences such as the presence of specific punctual or planar defects or of atoms with particular coordination features. Owing to these peculiarities, the chart also enabled tracking of the complex structural evolution in a reactive trajectory. In addition to visualization and analysis of complex structural landscapes, the presented approach offers a general, low-dimensional set of differentiable variables that has the potential to be used for exploration and enhanced sampling purposes.

What keeps nanopores boiling.

The liquid-to-vapor transition can occur under unexpected conditions in nanopores, opening the door to fundamental questions and new technologies. The physics of boiling in confinement is progressively introduced, starting from classical nucleation theory, passing through nanoscale effects, and terminating with the material and external parameters that affect the boiling conditions. The relevance of boiling in specific nanoconfined systems is discussed, focusing on heterogeneous lyophobic systems, chromatographic columns, and ion channels. The current level of control of boiling in nanopores enabled by microporous materials such as metal organic frameworks and biological nanopores paves the way to thrilling theoretical challenges and to new technological opportunities in the fields of energy, neuromorphic computing, and sensing.

Structural transformations in Cu, Ag, and Au metal nanoclusters.

Finite-temperature structures of Cu, Ag, and Au metal nanoclusters are calculated in the entire temperature range from 0 K to melting using a computational methodology that we proposed recently [M. Settem et al., Nanoscale 14, 939 (2022)]. In this method, Harmonic Superposition Approximation (HSA) and Parallel Tempering Molecular Dynamics (PTMD) are combined in a complementary manner. HSA is accurate at low temperatures and fails at higher temperatures. PTMD, on the other hand, effectively samples the high temperature region and melts. This method is used to study the size- and system-dependent competition between various structural motifs of Cu, Ag, and Au nanoclusters in the size range 1-2 nm. Results show that there are mainly three types of structural changes in metal nanoclusters, depending on whether a solid-solid transformation occurs. In the first type, the global minimum is the dominant motif in the entire temperature range. In contrast, when a solid-solid transformation occurs, the global minimum transforms either completely to a different motif or partially, resulting in the co-existence of multiple motifs. Finally, nanocluster structures are analyzed to highlight the system-specific differences across the three metals.

Integrated Approach Including Docking, MD Simulations, and Network Analysis Highlights the Action Mechanism of the Cardiac hERG Activator RPR260243.

hERG is a voltage-gated potassium channel involved in the heart contraction whose defections are associated with the cardiac arrhythmia Long QT Syndrome type 2. The activator RPR260243 (RPR) represents a possible candidate to pharmacologically treat LQTS2 because it enhances the opening of the channel. However, the molecular detail of its action mechanism remains quite elusive. Here, we address the problem using a combination of docking, molecular dynamics simulations, and network analysis. We show that the drug preferably binds at the interface between the voltage sensor and the pore, enhancing the canonical activation path and determining a whole-structure rearrangement of the channel that slightly impairs inactivation.

Mechanism of Water Intrusion into Flexible ZIF-8: Liquid Is Not Vapor.

Zeolitic Imidazolate Frameworks (ZIF) find application in storage and dissipation of mechanical energy. Their distinctive properties linked to their (sub)nanometer size and hydrophobicity allow for water intrusion only under high hydrostatic pressure. Here we focus on the popular ZIF-8 material investigating the intrusion mechanism in its nanoscale cages, which is the key to its rational exploitation in target applications. In this work, we used a joint experimental/theoretical approach combining in operando synchrotron experiments during high-pressure intrusion experiments, molecular dynamics simulations, and stochastic models to reveal that water intrusion into ZIF-8 occurs by a cascade filling of connected cages rather than a condensation process as previously assumed. The reported results allowed us to establish structure/function relations in this prototypical microporous material, representing an important step to devise design rules to synthesize porous media.

Optimization of the wetting-drying characteristics of hydrophobic metal organic frameworks via crystallite size: The role of hydrogen bonding between intruded and bulk liquid.

HYPOTHESIS: The behavior of Heterogeneous Lyophobic Systems (HLSs) comprised of a lyophobic porous material and a corresponding non-wetting liquid is affected by a variety of different structural parameters of the porous material. Dependence on exogenic properties such as crystallite size is desirable for system tuning as they are much more facilely modified. We explore the dependence of intrusion pressure and intruded volume on crystallite size, testing the hypothesis that the connection between internal cavities and bulk water facilitates intrusion via hydrogen bonding, a phenomenon that is magnified in smaller crystallites with a larger surface/volume ratio. EXPERIMENTS: Water intrusion/extrusion pressures and intrusion volume were experimentally measured for ZIF-8 samples of various crystallite sizes and compared to previously reported values. Alongside the practical research, molecular dynamics simulations and stochastic modeling were performed to illustrate the effect of crystallite size on the properties of the HLSs and uncover the important role of hydrogen bonding within this phenomenon. FINDINGS: A reduction in crystallite size led to a significant decrease of intrusion and extrusion pressures below 100 nm. Simulations indicate that this behavior is due to a greater number of cages being in proximity to bulk water for smaller crystallites, allowing cross-cage hydrogen bonds to stabilize the intruded state and lower the threshold pressure of intrusion and extrusion. This is accompanied by a reduction in the overall intruded volume. Simulations demonstrate that this phenomenon is linked to ZIF-8 surface half-cages exposed to water being occupied by water due to non-trivial termination of the crystallites, even at atmospheric pressure.

An atomistically informed multiscale approach to the intrusion and extrusion of water in hydrophobic nanopores.

Understanding intrusion and extrusion in nanoporous materials is a challenging multiscale problem of utmost importance for applications ranging from energy storage and dissipation to water desalination and hydrophobic gating in ion channels. Including atomistic details in simulations is required to predict the overall behavior of such systems because the statics and dynamics of these processes depend sensitively on microscopic features of the pore, such as the surface hydrophobicity, geometry, and charge distribution, and on the composition of the liquid. On the other hand, the transitions between the filled (intruded) and empty (extruded) states are rare events that often require long simulation times, which are difficult to achieve with standard atomistic simulations. In this work, we explored the intrusion and extrusion processes using a multiscale approach in which the atomistic details of the system, extracted from molecular dynamics simulations, informed a simple Langevin model of water intrusion/extrusion in the pore. We then used the Langevin simulations to compute the transition times at different pressures, validating our coarse-grained model by comparing it with nonequilibrium molecular dynamics simulations. The proposed approach reproduces experimentally relevant features such as the time and temperature dependence of the intrusion/extrusion cycles, as well as specific details about the shape of the cycle. This approach also drastically increases the timescales that can be simulated, reducing the gap between simulations and experiments and showing promise for more complex systems.

Classical nucleation of vapor between hydrophobic plates.

In this work, an extended classical nucleation theory (CNT), including line tension, is used to disentangle classical and non-classical effects in the nucleation of vapor from a liquid confined between two hydrophobic plates at a nanometer distance. The proposed approach allowed us to gauge, from the available simulation work, the importance of elusive nanoscale effects, such as line tension and non-classical modifications of the nucleation mechanism. Surprisingly, the purely macroscopic theory is found to be in quantitative accord with the microscopic data, even for plate distances as small as 2 nm, whereas in extreme confinement (<1.5 nm), the CNT approximations proved to be unsatisfactory. These results suggest how classical nucleation theory still offers a computationally inexpensive and predictive tool useful in all domains where nanoconfined evaporation occurs-including nanotechnology, surface science, and biology.

Correction: Tempering of Au nanoclusters: capturing the temperature-dependent competition among structural motifs.

Correction for 'Tempering of Au nanoclusters: capturing the temperature-dependent competition among structural motifs' by Manoj Settem et al., Nanoscale, 2022, 14, 939-952, https://doi.org/10.1039/D1NR05078H.

Noncanonical electromechanical coupling paths in cardiac hERG potassium channel.

Voltage-gated potassium channels are involved in many physiological processes such as nerve impulse transmission, the heartbeat, and muscle contraction. However, for many of them the molecular determinants of the gating mechanism remain elusive. Here, using a combination of theoretical and experimental approaches, we address this problem focusing on the cardiac hERG potassium channel. Network analysis of molecular dynamics trajectories reveals the presence of a kinematic chain of residues that couples the voltage sensor domain to the pore domain and involves the S4/S1 and S1/S5 subunit interfaces. Mutagenesis experiments confirm the role of these residues and interfaces in the activation and inactivation mechanisms. Our findings demonstrate the presence of an electromechanical transduction path crucial for the non-domain-swapped hERG channel gating that resembles the noncanonical path identified in domain-swapped K+ channels.

Two-Steps Versus One-Step Solidification Pathways of Binary Metallic Nanodroplets.

The solidification of AgCo, AgNi, and AgCu nanodroplets is studied by molecular dynamics simulations in the size range of 2-8 nm. All these systems tend to phase separate in the bulk solid with surface segregation of Ag. Despite these similarities, the simulations reveal clear differences in the solidification pathways. AgCo and AgNi already separate in the liquid phase, and they solidify in configurations close to equilibrium. They can show a two-step solidification process in which Co-/Ni-rich parts solidify at higher temperatures than the Ag-rich part. AgCu does not separate in the liquid and solidifies in one step, thereby remaining in a kinetically trapped state down to room temperature. The solidification mechanisms and the size dependence of the solidification temperatures are analyzed, finding qualitatively different behaviors in AgCo/AgNi compared to AgCu. These differences are rationalized by an analytical model.

Molecular dynamics simulations suggest possible activation and deactivation pathways in the hERG channel.

The elusive activation/deactivation mechanism of hERG is investigated, a voltage-gated potassium channel involved in severe inherited and drug-induced cardiac channelopathies, including the Long QT Syndrome. Firstly, the available structural data are integrated by providing a homology model for the closed state of the channel. Secondly, molecular dynamics combined with a network analysis revealed two distinct pathways coupling the voltage sensor domain with the pore domain. Interestingly, some LQTS-related mutations known to impair the activation/deactivation mechanism are distributed along the identified pathways, which thus suggests a microscopic interpretation of their role. Split channels simulations clarify a surprising feature of this channel, which is still able to gate when a cut is introduced between the voltage sensor domain and the neighboring helix S5. In summary, the presented results suggest possible activation/deactivation mechanisms of non-domain-swapped potassium channels that may aid in biomedical applications.

Tempering of Au nanoclusters: capturing the temperature-dependent competition among structural motifs.

A computational approach to determine the equilibrium structures of nanoclusters in the whole temperature range from 0 K to melting is developed. Our approach relies on Parallel Tempering Molecular Dynamics (PTMD) simulations complemented by Harmonic Superposition Approximation (HSA) calculations and global optimization searches, thus combining the accuracy of global optimization and HSA in describing the low-energy part of configuration space, together with the PTMD thorough sampling of high-energy configurations. This combined methodology is shown to be instrumental towards revealing the temperature-dependent structural motifs in Au nanoclusters of sizes 90, 147, and 201 atoms. The reported phenomenology is particularly rich, displaying a size- and temperature-dependent competition between the global energy minimum and other structural motifs. In the case of Au90 and Au147, the global minimum is also the dominant structure at finite temperatures. In contrast, the Au201 cluster undergoes a solid-solid transformation at low temperature (<200 K). Results indicate that PTMD and HSA very well agree at intermediate temperatures, between 300 and 400 K. For higher temperatures, PTMD gives an accurate description of equilibrium, while HSA fails in describing the melting range. On the other hand, HSA is more efficient in catching low-temperature structural transitions. Finally, we describe the elusive structures close to the melting region which can present complex and defective geometries, that are otherwise difficult to characterize through experimental imaging.

Unveiling the Gating Mechanism of CRAC Channel: A Computational Study.

CRAC channel is ubiquitous and its importance in the regulation of the immune system is testified by the severe immunodeficiencies caused by its mutations. In this work we took advantage of the availability of open and closed structures of this channel to run for the first time simulations of the whole gating process reaching the relevant time-scale with an enhanced sampling technique, Targeted Molecular Dynamics. Our simulations highlighted a complex allosteric propagation of the conformational change from peripheral helices, where the activator STIM1 binds, to the central pore helices. In agreement with mutagenesis data, our simulations revealed the key role of residue H206 whose displacement creates an empty space behind the hydrophobic region of the pore, thus releasing a steric brake and allowing the opening of the channel. Conversely, the process of pore closing culminates with the formation of a bubble that occludes the pore even in the absence of steric block. This mechanism, known as "hydrophobic gating", has been observed in an increasing number of biological ion channels and also in artificial nanopores. Our study therefore shows promise not only to better understand the molecular origin of diseases caused by disrupted calcium signaling, but also to clarify the mode of action of hydrophobically gated ion channels, possibly even suggesting strategies for the biomimetic design of synthetic nanopores.

Structure and dynamics of water confined in cylindrical nanopores with varying hydrophobicity.

We report a detailed study of the main structural and dynamical features of water confined in model Lennard-Jones nanopores with tunable hydrophobicity and finite length ([Formula: see text] Å). The generic model of cylindrical confinement used is able to reproduce the wetting features of a large class of technologically and biologically relevant systems spanning from crystalline nanoporous materials, to mesoporous silica and ion channels. The aim of this work is to discuss the influence of parameters such as wall hydrophobicity, temperature, and pore size on the structural and dynamical features of confined water. Our simulation campaign confirmed the existence of a core domain in which water displays bulk-like structural features even in extreme ([Formula: see text] Å) confinement, while dynamical properties were shown to depend non-trivially on the size and hydrophobicity of the pores. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

Giant Negative Compressibility by Liquid Intrusion into Superhydrophobic Flexible Nanoporous Frameworks.

Materials or systems demonstrating negative linear compressibility (NLC), whose size increases (decreases) in at least one of their dimensions upon compression (decompression) are very rare. Materials demonstrating this effect in all their dimensions, negative volumetric compressibility (NVC), are exceptional. Here, by liquid porosimetry and in situ neutron diffraction, we show that one can achieve exceptional NLC and NVC values by nonwetting liquid intrusion in flexible porous media, namely in the ZIF-8 metal-organic framework (MOF). Atomistic simulations show that the volumetric expansion is due to the presence of liquid in the windows connecting the cavities of ZIF-8. This discovery paves the way for designing novel materials with exceptional NLC and NVC at reasonable pressures suitable for a wide range of applications.

Exploring K v 1.2 Channel Inactivation Through MD Simulations and Network Analysis.

The KCNA2 gene encodes the K v 1.2 channel, a mammalian Shaker-like voltage-gated K+ channel, whose defections are linked to neuronal deficiency and childhood epilepsy. Despite the important role in the kinetic behavior of the channel, the inactivation remained hereby elusive. Here, we studied the K v 1.2 inactivation via a combined simulation/network theoretical approach that revealed two distinct pathways coupling the Voltage Sensor Domain and the Pore Domain to the Selectivity Filter. Additionally, we mutated some residues implicated in these paths and we explained microscopically their function in the inactivation mechanism by computing a contact map. Interestingly, some pathological residues shown to impair the inactivation lay on the paths. In summary, the presented results suggest two pathways as the possible molecular basis of the inactivation mechanism in the K v 1.2 channel. These pathways are consistent with earlier mutational studies and known mutations involved in neuronal channelopathies.