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.

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.

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.

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.

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.

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.

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.

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.

The N-Terminal Peptides of the Three Human Isoforms of the Mitochondrial Voltage-Dependent Anion Channel Have Different Helical Propensities.

The voltage-dependent anion channel (VDAC) is the main mitochondrial porin allowing the exchange of ions and metabolites between the cytosol and the mitochondrion. In addition, VDAC was found to actively interact with proteins playing a fundamental role in the regulation of apoptosis and being of central interest in cancer research. VDAC is a large transmembrane ß-barrel channel, whose N-terminal helical fragment adheres to the channel interior, partially closing the pore. This fragment is considered to play a key role in protein stability and function as well as in the interaction with apoptosis-related proteins. Three VDAC isoforms are differently expressed in higher eukaryotes, for which distinct and complementary roles are proposed. In this work, the folding propensity of their N-terminal fragments has been compared. By using multiple spectroscopic techniques, and complementing the experimental results with theoretical computer-assisted approaches, we have characterized their conformational equilibrium. Significant differences were found in the intrinsic helical propensity of the three peptides, decreasing in the following order: hVDAC2 > hVDAC3 > hVDAC1. In light of the models proposed in the literature to explain voltage gating, selectivity, and permeability, as well as interactions with functionally related proteins, our results suggest that the different chemicophysical properties of the N-terminal domain are possibly correlated to different functions for the three isoforms. The overall emerging picture is that a similar transmembrane water accessible conduit has been equipped with not identical domains, whose differences can modulate the functional roles of the three VDAC isoforms.

Opening pathways of the DNA clamps proliferating cell nuclear antigen and Rad9-Rad1-Hus1.

Proliferating cell nuclear antigen and the checkpoint clamp Rad9-Rad1-Hus1 topologically encircle DNA and act as mobile platforms in the recruitment of proteins involved in DNA damage response and cell cycle regulation. To fulfill these vital cellular functions, both clamps need to be opened and loaded onto DNA by a clamp loader complex-a process, which involves disruption of the DNA clamp's subunit interfaces. Herein, we compare the relative stabilities of the interfaces using the molecular mechanics Poisson-Boltzmann solvent accessible surface method. We identify the Rad9-Rad1 interface as the weakest and, therefore, most likely to open during clamp loading. We also delineate the dominant interface disruption pathways under external forces in multiple-trajectory steered molecular dynamics runs. We show that, similar to the case of protein folding, clamp opening may not proceed through a single interface breakdown mechanism. Instead, we identify an ensemble of opening pathways, some more prevalent than others, characterized by specific groups of contacts that differentially stabilize the regions of the interface and determine the spatial and temporal patterns of breakdown. In Rad9-Rad1-Hus1, the Rad9-Rad1 and Rad9-Hus1 interfaces share the same dominant unzipping pathway, whereas the Hus1-Rad1 interface is disrupted concertedly with no preferred directionality.

Coarse-grained modeling of protein unspecifically bound to DNA.

There is now a certain consensus that transcription factors (TFs) reach their target sites, where they regulate gene transcription, via a mechanism dubbed facilitated diffusion (FD). In FD, the TF cycles between events of 3D diffusion in solution (jumps), 1D diffusion along DNA (sliding), and small jumps (hopping), achieving association rates higher than for 3D diffusion alone. We investigate the FD phenomenology through molecular dynamics simulations in the framework of coarse-grained modeling. We show that, despite the crude approximations, the model generates, upon varying the equilibrium distance of the DNA-TF interaction, a phenomenology matching a number of experimental and numerical results obtained with more refined models. In particular, focusing on the kinematics of the process, we characterize the geometrical properties of TF trajectories during sliding. We find that sliding occurs via helical paths around the DNA helix, leading to a coupling of translation along the DNA axis with rotation around it. The 1D diffusion constant measured in simulations is found to be interwoven with the geometrical properties of sliding and we develop a simple argument that can be used to quantitatively reproduce the measured values.

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.

The conformational landscape of tartrate-based inhibitors of the TACE enzyme as revealed by Hamiltonian Replica Exchange simulation.

The inhibitors of the Tumor Necrosis Factor-α Converting Enzyme represent promising tools for the treatment of Rheumatoid Arthritis, Multiple Sclerosis and other autoimmune diseases. In this work, using Hamiltonian Replica Exchange Molecular Dynamics simulations and atomistic force field we perform an accurate structural characterization of a group of tartrate-based inhibitors. The simulations highlight a correlation between the conformational landscape in bulk solvent and inhibition potency. Since the structures in bulk solvent are much more compact than the crystallographic bound state, we formulate the hypothesis of a two-step docking mechanism: (i) formation of an intermediate between the compact, hydroxyl exposing conformations in solution and the catalytic zinc ion; (ii) structural rearrangement in the active site of TACE of the zinc-tethered drug in the final binding conformation.

Chemical-physical analysis of a tartrate model compound for TACE inhibition.

We have synthesized and done an extensive chemical-physical analysis of the behavior of a new compound, named MBET306, a synthetic precursor of the recently discovered tartrate-based inhibitors of the protein Tumor Necrosis factor-α Converting Enzyme (TACE). Experimental and theoretical data have shown that in water solution MBET306 is overwhelmingly found as a monoanion at physiological pH, in a conformation that differs substantially from that detected in the known co-crystal structures of MBET306 derivatives bound to TACE. The body of collected experimental and theoretical data indicates that the monoanionic species binds Zn(ii) inducing a strong stabilization of the crystal-like arrangement of the central tartrate zinc-binding group, lending support for a two step TACE docking mechanism via a zinc-bound intermediate. The thorough chemical-physical characterization of the conformational behavior of free and zinc-bound MBET306 in water bulk solution opens new avenues for the rational drug design of tartrate-based highly specific TACE inhibitors.

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.

Histone H4 acetylation differentially modulates arginine methylation by an in Cis mechanism.

Histone H4 undergoes extensive post-translational modifications (PTMs) at its N-terminal tail. Many of these PTMs profoundly affect the on and off status of gene transcription. The molecular mechanism by which histone PTMs modulate genetic and epigenetic processes is not fully understood. In particular, how a PTM mark affects the presence and level of other histone modification marks needs to be addressed and is essential for better understanding the molecular basis of histone code hypothesis. To dissect the interplaying relationship between different histone modification marks, we investigated how individual lysine acetylations and their different combinations at the H4 tail affect Arg-3 methylation in cis. Our data reveal that the effect of lysine acetylation on arginine methylation depends on the site of acetylation and the type of methylation. Although certain acetylations present a repressive impact on PRMT1-mediated methylation (type I methylation), lysine acetylation generally is correlated with enhanced methylation by PRMT5 (type II dimethylation). In particular, Lys-5 acetylation decreases the activity of PRMT1 but increases that of PRMT5. Furthermore, circular dichroism study and computer simulation demonstrate that hyperacetylation increases the content of ordered secondary structures at the H4 tail region. These findings provide new insights into the regulatory mechanism of Arg-3 methylation by H4 acetylation and unravel the complex intercommunications that exist between different the PTM marks in cis. The divergent activities of PRMT1 and PRMT5 with respect to different acetyl-H4 substrates suggest that type I and type II protein-arginine methyltransferases use distinct molecular determinants for substrate recognition and catalysis.

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.

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.