Bubbles enable volumetric negative compressibility in metastable elastocapillary systems.

Although coveted in applications, few materials expand when subject to compression or contract under decompression, i.e., exhibit negative compressibility. A key step to achieve such counterintuitive behaviour is the destabilisations of (meta)stable equilibria of the constituents. Here, we propose a simple strategy to obtain negative compressibility exploiting capillary forces both to precompress the elastic material and to release such precompression by a threshold phenomenon - the reversible formation of a bubble in a hydrophobic flexible cavity. We demonstrate that the solid part of such metastable elastocapillary systems displays negative compressibility across different scales: hydrophobic microporous materials, proteins, and millimetre-sized laminae. This concept is applicable to fields such as porous materials, biomolecules, sensors and may be easily extended to create unexpected material susceptibilities.

Optimization of Surface Functionalizations for Ring Resonator-Based Biosensors.

Liquid biopsy is expected to become widespread in the coming years thanks to point of care devices, which can include label-free biosensors. The surface functionalization of biosensors is a crucial aspect that influences their overall performance, resulting in the accurate, sensitive, and specific detection of target molecules. Here, the surface of a microring resonator (MRR)-based biosensor was functionalized for the detection of protein biomarkers. Among the several existing functionalization methods, a strategy based on aptamers and mercaptosilanes was selected as the most highly performing approach. All steps of the functionalization protocol were carefully characterized and optimized to obtain a suitable protocol to be transferred to the final biosensor. The functionalization protocol comprised a preliminary plasma treatment aimed at cleaning and activating the surface for the subsequent silanization step. Different plasma treatments as well as different silanes were tested in order to covalently bind aptamers specific to different biomarker targets, i.e., C-reactive protein, SARS-CoV-2 spike protein, and thrombin. Argon plasma and 1% v/v mercaptosilane were found as the most suitable for obtaining a homogeneous layer apt to aptamer conjugation. The aptamer concentration and time for immobilization were optimized, resulting in 1 µM and 3 h, respectively. A final passivation step based on mercaptohexanol was also implemented. The functionalization protocol was then evaluated for the detection of thrombin with a photonic biosensor based on microring resonators. The preliminary results identified the successful recognition of the correct target as well as some limitations of the developed protocol in real measurement conditions.

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.

Correlation, response and entropy approaches to allosteric behaviors: a critical comparison on the ubiquitin case.

Correlation analysis and its close variant principal component analysis are tools widely applied to predict the biological functions of macromolecules in terms of the relationship between fluctuation dynamics and structural properties. However, since this kind of analysis does not necessarily imply causation links among the elements of the system, its results run the risk of being biologically misinterpreted. By using as a benchmark the structure of ubiquitin, we report a critical comparison of correlation-based analysis with the analysis performed using two other indicators, response function and transfer entropy, that quantify the causal dependence. The use of ubiquitin stems from its simple structure and from recent experimental evidence of an allosteric control of its binding to target substrates. We discuss the ability of correlation, response and transfer-entropy analysis in detecting the role of the residues involved in the allosteric mechanism of ubiquitin as deduced by experiments. To maintain the comparison as much as free from the complexity of the modeling approach and the quality of time series, we describe the fluctuations of ubiquitin native state by the Gaussian network model which, being fully solvable, allows one to derive analytical expressions of the observables of interest. Our comparison suggests that a good strategy consists in combining correlation, response and transfer entropy, such that the preliminary information extracted from correlation analysis is validated by the two other indicators in order to discard those spurious correlations not associated with true causal dependencies.

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.

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.

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.

Field-Dependent Dehydration and Optimal Ionic Escape Paths for C2N Membranes.

Most analytic theories describing electrostatically driven ion transport through water-filled nanopores assume that the corresponding permeation barriers are bias-independent. While this assumption may hold for sufficiently wide pores under infinitely small bias, transport through subnanometer pores under finite bias is difficult to interpret analytically. Given recent advances in subnanometer pore fabrication and the rapid progress in detailed computer simulations, it is important to identify and understand the specific field-induced phenomena arising during ion transport. Here we consider an atomistic model of electrostatically driven ion permeation through subnanoporous C2N membranes. We analyze probability distributions of ionic escape trajectories and show that the optimal escape path switches between two different configurations depending on the bias magnitude. We identify two distinct mechanisms contributing to field-induced changes in transport-opposing barriers: a weak one arising from field-induced ion dehydration and a strong one due to the field-induced asymmetry of the hydration shells. The simulated current-voltage characteristics are compared with the solution of the 1D Nernst-Planck model. Finally, we show that the deviation of simulated currents from analytic estimates for large fields is consistent with the field-induced barriers and the observed changes in the optimal ion escape path.

Application of a Statistical and Linear Response Theory to Multi-Ion Na+ Conduction in NaChBac.

Biological ion channels are fundamental to maintaining life. In this manuscript we apply our recently developed statistical and linear response theory to investigate Na+ conduction through the prokaryotic Na+ channel NaChBac. This work is extended theoretically by the derivation of ionic conductivity and current in an electrochemical gradient, thus enabling us to compare to a range of whole-cell data sets performed on this channel. Furthermore, we also compare the magnitudes of the currents and populations at each binding site to previously published single-channel recordings and molecular dynamics simulations respectively. In doing so, we find excellent agreement between theory and data, with predicted energy barriers at each of the four binding sites of ∼4,2.9,3.6, and 4kT.

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.

Changes in Ion Selectivity Following the Asymmetrical Addition of Charge to the Selectivity Filter of Bacterial Sodium Channels.

Voltage-gated sodium channels (NaVs) play fundamental roles in eukaryotes, but their exceptional size hinders their structural resolution. Bacterial NaVs are simplified homologues of their eukaryotic counterparts, but their use as models of eukaryotic Na+ channels is limited by their homotetrameric structure at odds with the asymmetric Selectivity Filter (SF) of eukaryotic NaVs. This work aims at mimicking the SF of eukaryotic NaVs by engineering radial asymmetry into the SF of bacterial channels. This goal was pursued with two approaches: the co-expression of different monomers of the NaChBac bacterial channel to induce the random assembly of heterotetramers, and the concatenation of four bacterial monomers to form a concatemer that can be targeted by site-specific mutagenesis. Patch-clamp measurements and Molecular Dynamics simulations showed that an additional gating charge in the SF leads to a significant increase in Na+ and a modest increase in the Ca2+ conductance in the NavMs concatemer in agreement with the behavior of the population of random heterotetramers with the highest proportion of channels with charge -5e. We thus showed that charge, despite being important, is not the only determinant of conduction and selectivity, and we created new tools extending the use of bacterial channels as models of eukaryotic counterparts.

Exploring the pore charge dependence of K+ and Cl- permeation across a graphene monolayer: a molecular dynamics study.

Selective permeation through graphene nanopores is attracting increasing interest as an efficient and cost-effective technique for water desalination and purification. In this work, using umbrella sampling and molecular dynamics simulations with constant electric field, we analyze the influence of pore charge on potassium and chloride ion permeation. As pore charge is increased, the barrier of the potential of mean force (PMF) gradually decreases until it turns into a well split in two subminima. While in the case of K+ this pattern can be explained as an increasing electrostatic compensation of the desolvation cost, in the case of Cl- the pattern can be attributed to the accumulation of a concentration polarization layer of potassium ions screening pore charge. The analysis of potassium PMFs in terms of forces revealed a conflicting influence on permeation of van der Waals and electrostatic forces that both undergo an inversion of their direction as pore charge is increased. Even if the most important transition involves the interplay between the electrostatic forces exerted by graphene and water, the simulations also revealed an important role of the changing distribution of potassium and chloride ions. The influence of pore charge on the orientation of water molecules was also found to affect the van der Waals forces they exert on potassium.

Different roles for aspartates and glutamates for cation permeation in bacterial sodium channels.

A key driving force for ion channel selectivity is represented by the negative charge of the Selectivity Filter carried by aspartate (D) and glutamate (E) residues. However, the structural effects and specific properties of D and E residues have not been extensively studied. In order to investigate this issue we studied the mutants of NaChBac channel with all possible combinations of D and E in the charged rings in position 191 and 192. Electrophysiological measurements showed significant Ca2+ currents only when position 191 was occupied by E. Equilibrium Molecular Dynamics simulations revealed the existence of two binding sites, corresponding to the charged rings and another one, more internal, at the level of L190. The simulations showed that the ion in the innermost site can interact with the residue in position 191 only when this is glutamate. Based on the MD simulations, we suggest that a D in position 191 leads to a high affinity Ca2+ block site resulting from a significant drop in the free energy of binding for an ion moving between the binding sites; in contrast, the free energy change is more gradual when an E residue occupies position 191, resulting in Ca2+ permeability. This scenario is consistent with the model of ion channel selectivity through stepwise changes in binding affinity proposed by Dang and McCleskey. Our study also highlights the importance of the structure of the selectivity filter which should contribute to the development of more detailed physical models for ion channel selectivity.

yVDAC2, the second mitochondrial porin isoform of Saccharomyces cerevisiae.

The yeast Saccharomyces cerevisiae genome is endowed with two distinct isoforms of Voltage-Dependent Anion Channel (VDAC). The isoform yVDAC2 is currently understudied with respect to the best known yVDAC1. Yet, since the discovery, the function of yVDAC2 was unclear, leading to the hypothesis that it might be devoid of a channel function. In this work we have elucidated, by bioinformatics modeling and electrophysiological analysis, the functional activity of yVDAC2. The conformation of yVDAC2 and, for comparison, of yVDAC1 were modeled using a multiple template approach involving mouse, human and zebrafish structures and both showed to arrange the sequences as the typical 19-stranded VDAC ß-barrel. Molecular dynamics simulations showed that yVDAC2, in comparison with yVDAC1, has a different number of permeation paths of potassium and chloride ions. yVDAC2 protein was over-expressed in the S. cerevisiae cells depleted of functional yVDAC1 (Δpor1 mutant) and, after purification, it was reconstituted in artificial membranes (planar lipid bilayer (PLB) system). The protein displayed channel-forming activity and the calculated conductance, voltage-dependence and ion selectivity values were similar to those of yVDAC1 and other members of VDAC family. This is the first time that yVDAC2 channel features are detected and characterized.

On the selectivity of the NaChBac channel: an integrated computational and experimental analysis of sodium and calcium permeation.

Ion channel selectivity is essential for their function, yet the molecular basis of a channel's ability to select between ions is still rather controversial. In this work, using a combination of molecular dynamics simulations and electrophysiological current measurements we analyze the ability of the NaChBac channel to discriminate between calcium and sodium. Our simulations show that a single calcium ion can access the Selectivity Filter (SF) interacting so strongly with the glutamate ring so as to remain blocked inside. This is consistent with the tiny calcium currents recorded in our patch-clamp experiments. Two reasons explain this scenario. The first is the higher free energy of ion/SF binding of Ca2+ with respect to Na+. The second is the strong electrostatic repulsion exerted by the resident ion that turns back a second potentially incoming Ca2+, preventing the knock-on permeation mechanism. Finally, we analyzed the possibility of the Anomalous Mole Fraction Effect (AMFE), i.e. the ability of micromolar Ca2+ concentrations to block Na+ currents. Current measurements in Na+/Ca2+ mixed solutions excluded the AMFE, in agreement with metadynamics simulations showing the ability of a sodium ion to by-pass and partially displace the resident calcium. Our work supports a new scenario for Na+/Ca2+ selectivity in the bacterial sodium channel, challenging the traditional notion of an exclusion mechanism strictly confining Ca2+ ions outside the channel.

Sodium Binding Sites and Permeation Mechanism in the NaChBac Channel: A Molecular Dynamics Study.

NaChBac was the first discovered bacterial sodium voltage-dependent channel, yet computational studies are still limited due to the lack of a crystal structure. In this work, a pore-only construct built using the NavMs template was investigated using unbiased molecular dynamics and metadynamics. The potential of mean force (PMF) from the unbiased run features four minima, three of which correspond to sites IN, CEN, and HFS discovered in NavAb. During the run, the selectivity filter (SF) is spontaneously occupied by two ions, and frequent access of a third one is often observed. In the innermost sites IN and CEN, Na+ is fully hydrated by six water molecules and occupies an on-axis position. In site HFS sodium interacts with a glutamate and a serine from the same subunit and is forced to adopt an off-axis placement. Metadynamics simulations biasing one and two ions show an energy barrier in the SF that prevents single-ion permeation. An analysis of the permeation mechanism was performed both computing minimum energy paths in the axial-axial PMF and through a combination of Markov state modeling and transition path theory. Both approaches reveal a knock-on mechanism involving at least two but possibly three ions. The currents predicted from the unbiased simulation using linear response theory are in excellent agreement with single-channel patch-clamp recordings.

VDAC3 as a sensor of oxidative state of the intermembrane space of mitochondria: the putative role of cysteine residue modifications.

Voltage-Dependent Anion selective Channels (VDAC) are pore-forming mitochondrial outer membrane proteins. In mammals VDAC3, the least characterized isoform, presents a set of cysteines predicted to be exposed toward the intermembrane space. We find that cysteines in VDAC3 can stay in different oxidation states. This was preliminary observed when, in our experimental conditions, completely lacking any reducing agent, VDAC3 presented a pattern of slightly different electrophoretic mobilities. This observation holds true both for rat liver mitochondrial VDAC3 and for recombinant and refolded human VDAC3. Mass spectroscopy revealed that cysteines 2 and 8 can form a disulfide bridge in native VDAC3. Single or combined site-directed mutagenesis of cysteines 2, 8 and 122 showed that the protein mobility in SDS-PAGE is influenced by the presence of cysteine and by the redox status. In addition, cysteines 2, 8 and 122 are involved in the stability control of the pore as shown by electrophysiology, complementation assays and chemico-physical characterization. Furthermore, a positive correlation between the pore conductance of the mutants and their ability to complement the growth of porin-less yeast mutant cells was found. Our work provides evidence for a complex oxidation pattern of a mitochondrial protein not directly involved in electron transport. The most likely biological meaning of this behavior is to buffer the ROS load and keep track of the redox level in the inter-membrane space, eventually signaling it through conformational changes.

A computational study of ion current modulation in hVDAC3 induced by disulfide bonds.

The human VDAC channel exists in three isoforms characterized by high sequence homology and structural similarity. Yet the function and mode of action of hVDAC3 are still elusive. The presence of six surface cysteines exposed to the oxidizing environment of the mitochondrial inter-membrane space suggests the possible establishment of intramolecular disulfide bonds. Two natural candidates for disulfide bridge formation are Cys2 and Cys8 that, located on the flexible N-terminal domain, can easily come in contact. A third potentially important residue is Cys122 that is close to Cys2 in the homology model of VDAC3. Here we analyzed the impact of SS bonds through molecular dynamics simulations of derivatives of hVDAC3 (dubbed SS-2-8, SS-2-122, SS-8-122) including a single disulfide bond. Simulations showed that in SS-8-122, the fragment 1-7 crosses the top part of the barrel partially occluding the pore and causing a 20% drop of conductance. In order to identify other potential channel-occluding disulfide bonds, we used a set of neural networks and structural bioinformatics algorithms, after filtering with the steric constraints imposed by the 3D-structure. We identified other three species, namely SS-8-65, SS-2-36 and SS-8-36. While the conductance of SS-8-65 and SS-2-36 is about 30% lower than that of the species without disulfide bonds, the conductance of SS-8-36 was 40-50% lower. The results show how VDAC3 is able to modulate its pore size and current by exploiting the mobility of the N-terminal and forming, upon external stimuli, disulfide bridges with cysteine residues located on the barrel and exposed to the inter-membrane space.

Efficient nonequilibrium method for binding free energy calculations in molecular dynamics simulations.

We introduce an effective technique for the calculation of the binding free energy in drug-receptor systems using nonequilibrium molecular dynamics and application of the Jarzynski theorem. In essence, this novel methodology constitutes the nonequilibrium adaptation of an ancient free energy perturbation technique called Double Annihilation Method, invented more than 25 years ago [J. Chem. Phys. 1988, 89, 3742-3746] and upon which modern approaches of binding free energy computation in drug-receptor systems are heavily based. The proposed computational approach, termed Fast Switching Double Annihilation Methods (FS-DAM) in honor of its ancient ancestor, is applied to a prototypical example system with multiple binding sites, proving its computational potential and versatility in unraveling multiple site or allosteric binding processes.