Mechanism of voltage sensing in Ca2+- and voltage-activated K+ (BK) channels.

In neurosecretion, allosteric communication between voltage sensors and Ca2+ binding in BK channels is crucially involved in damping excitatory stimuli. Nevertheless, the voltage-sensing mechanism of BK channels is still under debate. Here, based on gating current measurements, we demonstrate that two arginines in the transmembrane segment S4 (R210 and R213) function as the BK gating charges. Significantly, the energy landscape of the gating particles is electrostatically tuned by a network of salt bridges contained in the voltage sensor domain (VSD). Molecular dynamics simulations and proton transport experiments in the hyperpolarization-activated R210H mutant suggest that the electric field drops off within a narrow septum whose boundaries are defined by the gating charges. Unlike Kv channels, the charge movement in BK appears to be limited to a small displacement of the guanidinium moieties of R210 and R213, without significant movement of the S4.

The voltage sensor is responsible for ΔpH dependence in Hv1 channels.

The dissipation of acute acid loads by the voltage-gated proton channel (Hv1) relies on regulating the channel's open probability by the voltage and the ΔpH across the membrane (ΔpH = pHex - pHin). Using monomeric Ciona-Hv1, we asked whether ΔpH-dependent gating is produced during the voltage sensor activation or permeation pathway opening. A leftward shift of the conductance-voltage (G-V) curve was produced at higher ΔpH values in the monomeric channel. Next, we measured the voltage sensor pH dependence in the absence of a functional permeation pathway by recording gating currents in the monomeric nonconducting D160N mutant. Increasing the ΔpH leftward shifted the gating charge-voltage (Q-V) curve, demonstrating that the ΔpH-dependent gating in Hv1 arises by modulating its voltage sensor. We fitted our data to a model that explicitly supposes the Hv1 voltage sensor free energy is a function of both the proton chemical and the electrical potential. The parameters obtained showed that around 60% of the free energy stored in the ΔpH is coupled to the Hv1 voltage sensor activation. Our results suggest that the molecular mechanism underlying the Hv1 ΔpH dependence is produced by protons, which alter the free-energy landscape around the voltage sensor domain. We propose that this alteration is produced by accessibility changes of the protons in the Hv1 voltage sensor during activation.

Size Matters: Free-Energy Calculations of Amino Acid Adsorption over Pristine Graphene.

There is still growing interest in graphene interactions with proteins, both for its possible biological applications and due to concerns over detrimental effects at the cellular level. As with any process involving proteins, an understanding of amino acid composition is desirable. In this work, we systematically studied the adsorption process of amino acids onto pristine graphene via rigorous free-energy calculations. We characterized the free energy, potential energy, and entropy of the adsorption of all proteinogenic amino acids. The energetic components were further separated into pair interaction contributions. A linear correlation was found between the free energy and the solvent accessible surface area change during adsorption (ΔSASAads) over pristine graphene and uncharged amino acids. Free energies over pristine graphene were compared with adsorption onto graphene oxide, finding an almost complete loss of the favorability of amino acid adsorption onto graphene. Finally, the correlation with ΔSASAads was used to successfully predict the free energy of adsorption of several penta-l-peptides in different structural states and sequences. Due to the relative ease of calculating the ΔSASAads compared to free-energy calculations, it could prove to be a cost-effective predictor of the free energy of adsorption for proteins onto nonpolar surfaces.

Protein-Mediated Electroporation in a Cardiac Voltage-Sensing Domain Due to an nsPEF Stimulus.

This study takes a step in understanding the physiological implications of the nanosecond pulsed electric field (nsPEF) by integrating molecular dynamics simulations and machine learning techniques. nsPEF, a state-of-the-art technology, uses high-voltage electric field pulses with a nanosecond duration to modulate cellular activity. This investigation reveals a relatively new and underexplored phenomenon: protein-mediated electroporation. Our research focused on the voltage-sensing domain (VSD) of the NaV1.5 sodium cardiac channel in response to nsPEF stimulation. We scrutinized the VSD structures that form pores and thereby contribute to the physical chemistry that governs the defibrillation effect of nsPEF. To do so, we conducted a comprehensive analysis involving the clustering of 142 replicas simulated for 50 ns under nsPEF stimuli. We subsequently pinpointed the representative structures of each cluster and computed the free energy between them. We find that the selected VSD of NaV1.5 forms pores under nsPEF stimulation, but in a way that significant differs from the traditional VSD opening. This study not only extends our understanding of nsPEF and its interaction with protein channels but also adds a new effect to further study.

On the effects of induced polarizability at the water-graphene interface via classical charge-on-spring models.

Molecular models of the water-graphene interaction are essential to describe graphene in condensed phases. Different challenges are associated with the generation of these models, in particular π-π and dispersion interactions; thus quantum and classical models have been developed and due to the numerical efficiency of the latter, they have been extensively employed. In this work, we have systematically studied, via molecular dynamics, two polarizable graphene models, denominated CCCP and CCCPD, employing the charge-on-spring model of the GROMOS forcefield, both being compatible with the polarizable water models COS/G2 and COS/D2, respectively. These models were compared with non-polarizable graphene and SPC water models. We focused the study on the water-graphene interface in two distinct systems and under the influence of an electric field: one composed of graphene immersed in water and the other composed of graphene with a water droplet above it. In the former, the orientation of water close to the graphene layer is affected by polarizable graphene in comparison to non-polarizable graphene. This effect is emphasised when an electric field is applied. In the latter, carbon polarizability reduced water contact angles, but graphene retained its hydrophobicity and the computed angles are within the experimental data. Given the significant extra computational cost, the use of polarizable models instead of the traditional fixed-charged approach for the graphene-water interaction may be justified when polarizability effects are relevant, for example, when applying relatively strong fields or in very anisotropic systems, such as the vacuum-bulk interface, as these models are more responsive to such conditions.

Aquaporin Gating: A New Twist to Unravel Permeation through Water Channels.

Aquaporins (AQPs) are small transmembrane tetrameric proteins that facilitate water, solute and gas exchange. Their presence has been extensively reported in the biological membranes of almost all living organisms. Although their discovery is much more recent than ion transport systems, different biophysical approaches have contributed to confirm that permeation through each monomer is consistent with closed and open states, introducing the term gating mechanism into the field. The study of AQPs in their native membrane or overexpressed in heterologous systems have experimentally demonstrated that water membrane permeability can be reversibly modified in response to specific modulators. For some regulation mechanisms, such as pH changes, evidence for gating is also supported by high-resolution structures of the water channel in different configurations as well as molecular dynamics simulation. Both experimental and simulation approaches sustain that the rearrangement of conserved residues contributes to occlude the cavity of the channel restricting water permeation. Interestingly, specific charged and conserved residues are present in the environment of the pore and, thus, the tetrameric structure can be subjected to alter the positions of these charges to sustain gating. Thus, is it possible to explore whether the displacement of these charges (gating current) leads to conformational changes? To our knowledge, this question has not yet been addressed at all. In this review, we intend to analyze the suitability of this proposal for the first time.

Entropy deepens loading chemical potentials of small alcohols by narrow carbon nanotubes.

Small alcohol confinement within narrow carbon nanotubes has been extensively and systematically studied via rigorous free-energy calculations. Employing molecular dynamics simulations, thermodynamic integration and thermodynamic cycling, the loading process of methanol and ethanol from aqueous solution into (6,6), (7,7) and (8,8) single-walled carbon nanotubes was computed and decomposed into its entropic and energetic terms. For all tubes and alcohols, loading is favoured from infinite dilution in water; for the same alcohol, wider tubes allow for the formation of a collective dipole which is cooperative in terms of electrostatics and reduce the rotational freedom of the loaded particles; narrow tubes only permit the formation of dipole-dipole dimers instead, with a (rotational) entropic gain that compensates for the loss of long-range dipole-dipole interactions. The latter renders deeper loading chemical potentials for narrower tubes when partitioning small alcohols from aqueous solution and it is a clear example of an entropy-energy compensation phenomenon.

Controlling ionic conductivity through transprotein electropores in human aquaporin 4: a non-equilibrium molecular-dynamics study.

Electroporation is a matter of intensive ongoing research interest, and a much-neglected topic in trans-membrane proteins, particularly in view of such promising potential applications in medicine and biotechnology. In particular, selected such novel and exciting applications are predicated on controlling ionic conductivity through electro-pores. Here, we scrutinise the mechanisms of ions' electric conductivity, by means of structural rearrangements, through quasi-stable electro-pores through human-AQP4 as a well-representative prototype of trans-membrane ionic conduction, achieving exquisite control over ionic permeability manipulated by the application of intense static electric fields.

Human aquaporin 4 gating dynamics under axially oriented electric-field impulses: A non-equilibrium molecular-dynamics study.

Human aquaporin 4 has been studied using non-equilibrium molecular dynamics simulations in the absence and presence of pulses of external electric fields. The pulses were 100 ns in duration and 0.005-0.015 V/Å in intensity acting along the pores' axes. Water diffusivity and the dipolar response of various residues of interest within the pores have been studied. Results show relatively little change in levels of water permeability per se within aquaporin channels during axially oriented field impulses, although care must be taken with regard to statistical certainty. However, the spatial variation of water permeability vis-à-vis electric-field intensity within the milieu of the channels, as revealed by heterogeneity in diffusivity-map gradients, indicates the possibility of somewhat enhanced diffusivity, owing to several residues being affected substantially by external fields, particularly for HIS 201 and 95 and ILE 93. This has the effect of increasing slightly intra-pore water diffusivity in the "pore-mouths" locale, albeit rendering it more spatially uniform overall vis-à-vis zero-field conditions (via manipulation of the selectivity filter).

Human Aquaporin 4 Gating Dynamics under Perpendicularly-Oriented Electric-Field Impulses: A Molecular Dynamics Study.

Human aquaporin 4 has been studied using molecular dynamics (MD) simulations in the absence and presence of pulses of external static electric fields. The pulses were 10 ns in duration and 0.012-0.065 V/Å in intensity acting along both directions perpendicular to the pores. Water permeability and the dipolar response of all residues of interest (including the selectivity filter) within the pores have been studied. Results showed decreased levels of water osmotic permeability within aquaporin channels during orthogonally-oriented field impulses, although care must be taken with regard to statistical certainty. This can be explained observing enhanced "dipolar flipping" of certain key residues, especially serine 211, histidine 201, arginine 216, histidine 95 and cysteine 178. These residues are placed at the extracellular end of the pore (serine 211, histidine 201, and arginine 216) and at the cytoplasm end (histidine 95 and cysteine 178), with the key role in gating mechanism, hence influencing water permeability.

Characterization of a novel water pocket inside the human Cx26 hemichannel structure.

Connexins (Cxs) are a family of vertebrate proteins constituents of gap junction channels (GJCs) that connect the cytoplasm of adjacent cells by the end-to-end docking of two Cx hemichannels. The intercellular transfer through GJCs occurs by passive diffusion allowing the exchange of water, ions, and small molecules. Despite the broad interest to understand, at the molecular level, the functional state of Cx-based channels, there are still many unanswered questions regarding structure-function relationships, perm-selectivity, and gating mechanisms. In particular, the ordering, structure, and dynamics of water inside Cx GJCs and hemichannels remains largely unexplored. In this work, we describe the identification and characterization of a believed novel water pocket-termed the IC pocket-located in-between the four transmembrane helices of each human Cx26 (hCx26) monomer at the intracellular (IC) side. Using molecular dynamics (MD) simulations to characterize hCx26 internal water structure and dynamics, six IC pockets were identified per hemichannel. A detailed characterization of the dynamics and ordering of water including conformational variability of residues forming the IC pockets, together with multiple sequence alignments, allowed us to propose a functional role for this cavity. An in vitro assessment of tracer uptake suggests that the IC pocket residue Arg-143 plays an essential role on the modulation of the hCx26 hemichannel permeability.

Entropic and enthalpic contributions to stereospecific ligand binding from enhanced sampling methods.

The stereoselective binding of R- and S-propranolol to the metabolic enzyme cytochrome P450 2D6 and its mutant F483A was studied using various computational approaches. Previously reported free-energy differences from Hamiltonian replica exchange simulations, combined with thermodynamic integration, are compared to the one-step perturbation approach, combined with local-elevation enhanced sampling, and an excellent agreement between methods was obtained. Further, the free-energy differences are interpreted in terms of enthalpic and entropic contributions where it is shown that exactly compensating contributions obscure a molecular interpretation of differences in the affinity while various reduced terms allow a more detailed analysis, which agree with heuristic observations on the interactions.

On the thermodynamics of carbon nanotube single-file water loading: free energy, energy and entropy calculations.

Single-file water chains confined in carbon nanotubes have been extensively studied using molecular dynamics simulations. Specifically, the pore loading process of periodic (6,6) and (5,5) single-walled carbon nanotubes was thermodynamically characterized by means of free-energy calculations at every loading state and compared to bulk water employing thermodynamic cycles. Long simulations of each end-state allowed for the partitioning of the free energy into its energetic and entropic components. The calculations revealed that the initial loading states are dominated by entropic (both translational and rotational) components, whereas the latter stages are energetically driven by strong dipolar interactions among the water molecules in the file.

Dewetting transition assisted clearance of (NFGAILS) amyloid fibrils from cell membranes by graphene.

Clearance of partially ordered oligomers and monomers deposited on cell membrane surfaces is believed to be an effective route to alleviate many potential protein conformational diseases (PCDs). With large-scale all-atom molecular dynamics simulations, here we show that graphene nanosheets can easily and quickly win a competitive adsorption of human islet amyloid polypeptides (hIAPP22-28) NFGAILS and associated fibrils against cell membrane, due to graphene's unique two-dimensional, highly hydrophobic surface with its all-sp(2) hybrid structure. A nanoscale dewetting transition was observed at the interfacial region between the fibril (originally deposited on the membrane) and the graphene nanosheet, which significantly assisted the adsorption of fibrils onto graphene from the membrane. The π-π stacking interaction between Phe23 and graphene played a crucial role, providing the driving force for the adsorption at the graphene surface. This study renders new insight towards the importance of water during the interactions between amyloid peptides, the phospholipidic membrane, and graphene, which might shed some light on future developments of graphene-based nanomedicine for preventing/curing PCDs like type II diabetes mellitus.

Lipid A from lipopolysaccharide recognition: structure, dynamics and cooperativity by molecular dynamics simulations.

Molecular dynamics simulations of Lipid A and its natural precursor Lipid IVA from E.coli have been carried out free in solution, bound to the myeliod differentiation protein 2 (MD2) and in the complex of MD2 with the toll like receptor 4 (TLR4). In addition, simulations of the ligand free MD2 and MD2-TLR4 complex were performed. A structural and energetic characterization of the bound and unbound states of Lipid A/IVA was generated. As the crystal structures depict, the main driving force for MD2-Lipid A/IVA are the hydrophobic interactions between the aliphatic tails and the MD2 cavity. The charged phosphate groups do strongly interact with positively charged residues, located at the surface of MD2. However, they are not essential for keeping the lipids in the cavity, indicating a more prominent role in binding recognition and ionic interactions with TLR4 at the MD2/TLR4 interface. Interestingly, in the absence of any ligand MD2 rapidly closes, blocking the binding cavity. The presence of TLR4, though changing the dynamics, was not able to impede the aforementioned closing event. We hypothesize that fluctuations of the H1 region are essential for this phenomenon, and it is plausible that an equilibrium between the open and closed states exists, although the lengths of our simulations are not sufficient to encompass the reversible process. The MD2/Lipid A-TLR4 complex simulations show that the presence of the ligand energetically stabilizes the complex relative to the ligand-free structures, indicating cooperativity in the binding process.

Free-energy differences between states with different conformational ensembles.

Multiple conformations separated by high-energy barriers represent a challenging problem in free-energy calculations due to the difficulties in achieving adequate sampling. We present an application of thermodynamic integration (TI) in conjunction with the local elevation umbrella sampling (LE/US) method to improve convergence in alchemical free-energy calculations. TI-LE/US was applied to the guanosine triphosphate (GTP) to 8-Br-GTP perturbation, molecules that present high-energy barriers between the anti and syn states and that have inverted preferences for those states. The convergence and reliability of TI-LE/US was assessed by comparing with previous results using the enhanced-sampling one-step perturbation (OSP) method. A linear interpolation of the end-state biasing potentials was sufficient to dramatically improve sampling along the chosen reaction coordinate. Conformational free-energy differences were also computed for the syn and anti states and compared to experimental and theoretical results. Additionally, a coupled OSP with LE/US was carried out, allowing the calculation of conformational and alchemical free energies of GTP and 8-substituted GTP analogs.

Human aquaporin 4 gating dynamics under and after nanosecond-scale static and alternating electric-field impulses: a molecular dynamics study of field effects and relaxation.

Water self-diffusion and the dipolar response of the selectivity filter within human aquaporin 4 have been studied using molecular dynamics (MD) simulations in the absence and presence of pulses of external static and alternating electric fields. The pulses were approximately 50 and 100 ns in duration and 0.0065 V/Å in (r.m.s.) intensity and were either static or else 2.45 or 100 GHz in frequency and applied both along and perpendicular to the channels. In addition, the relaxation of the aquaporin, water self-diffusion and gating dynamics following cessation of the impulses was studied. In previous work it was determined that switches in the dihedral angle of the selectivity filter led to boosting of water permeation events within the channels, in the presence of identical external static and alternating electric fields, although applied continuously. Here the application of field impulses (and subsequently, upon removal) has shown that it is the dipolar orientation of the histidine-201 residue in the selectivity filter which governs the dihedral angle, and hence influences water self-diffusion; this constitutes an appropriate order parameter. The dipolar response of this residue to the applied field leads to the adoption of four distinct states, which we modelled as time-homogeneous Markov jump processes, and may be distinguished in the potential of mean force (PMF) as a function of the dipolar orientation of histidine-201. The observations of enhanced "dipolar flipping" of H201 serve to explain increased levels of water self-diffusion within aquaporin channels during, and immediately following, field impulses, although the level of statistical certainty here is lower. Given the appreciable size of the energy barriers evident in PMFs computed directly from deterministic MD (whether in the absence or presence of external fields), metadynamics calculations were undertaken to explore the free-energy landscape of histidine-201 orientation with greater accuracy and precision. These indicate that electric fields do alter the free-energy profile of the H201 side-chain orientation, wherein a perturbation of the symmetric bimodal state evident in the zero-field case is observed. These effects are dependent on the field intensities.

In silico approaches to develop new phenyl-pyrimidines as glycogen synthase kinase 3 (GSK-3) inhibitors with halogen-bonding capabilities: 3D-QSAR CoMFA/CoMSIA, molecular docking and molecular dynamics studies.

Glycogen synthase kinase 3 (GSK-3) is involved in different diseases, such as manic-depressive illness, Alzheimer's disease and cancer. Studies have shown that insulin inhibits GSK-3 to keep glycogen synthase active. Inhibiting GSK-3 may have an indirect pro-insulin effect by favouring glycogen synthesis. Therefore, the development of GSK-3 inhibitors can be a useful alternative for the treatment of type II diabetes. Aminopyrimidine derivatives already proved to be interesting GSK-3 inhibitors. In the current study, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) have been performed on a series of 122 aminopyrimidine derivatives in order to generate a robust model for the rational design of new compounds with promising antidiabetic activity. The q2 values obtained for the best CoMFA and CoMSIA models have been 0.563 and 0.598, respectively. In addition, the r2 values have been 0.823 and 0.925 for CoMFA and CoMSIA, respectively. The models were statistically validated, and from the contour maps analysis, a proposal of 10 new compounds has been generated, with predicted pIC50 higher than 9. The final contribution of our work is that: (a) we provide an extensive structure-activity relationship for GSK-3 inhibitory pyrimidines; and (b) these models may speed up the discovery of GSK-3 inhibitors based on the aminopyrimidine scaffold. Finally, we carried out docking and molecular dynamics studies of the two best candidates, which were shown to establish halogen-bond interactions with the enzyme.Communicated by Ramaswamy H. Sarma.

Mechanosensitive aquaporins.

Cellular systems must deal with mechanical forces to satisfy their physiological functions. In this context, proteins with mechanosensitive properties play a crucial role in sensing and responding to environmental changes. The discovery of aquaporins (AQPs) marked a significant breakthrough in the study of water transport. Their transport capacity and regulation features make them key players in cellular processes. To date, few AQPs have been reported to be mechanosensitive. Like mechanosensitive ion channels, AQPs respond to tension changes in the same range. However, unlike ion channels, the aquaporin's transport rate decreases as tension increases, and the molecular features of the mechanism are unknown. Nevertheless, some clues from mechanosensitive ion channels shed light on the AQP-membrane interaction. The GxxxG motif may play a critical role in the water permeation process associated with structural features in AQPs. Consequently, a possible gating mechanism triggered by membrane tension changes would involve a conformational change in the cytoplasmic extreme of the single file region of the water pathway, where glycine and histidine residues from loop B play a key role. In view of their transport capacity and their involvement in relevant processes related to mechanical forces, mechanosensitive AQPs are a fundamental piece of the puzzle for understanding cellular responses.

Diffusion and interactions of carbon dioxide and oxygen in the vicinity of the active site of Rubisco: molecular dynamics and quantum chemical studies.

Molecular dynamics (MD) at the molecular mechanical level and geometry optimisation at the quantum mechanical level have been performed to investigate the transport and fixation of oxygen and carbon dioxide in the cavity of ribulose-1,5-bisphosphate carboxylase/oxygenase, or Rubisco. Multiple MD simulations have been carried out to study the diffusive behaviour of O(2) and CO(2) molecules from the Mg(2+) cation in Rubisco at 298 K and 1 bar, being one step in the overall process of carboxylation/oxygenation in Rubisco. In addition to this work, in order to gain additional perspective on the role of chemical reaction rates and thermodynamics, oxygen, and carbon dioxide uptake mechanisms have also been investigated by the aid of quantum chemical calculations. The results indicate that the activation barrier for carboxylation is slightly lower than that of oxygenation. This agrees qualitatively with experimental findings, and rationalises the observed competition between both catalytic processes in nature. Finally, the longer-lived persistence of CO(2) in the vicinity of the active centre (i.e., slower self-diffusion) may serve to explain, in part, why carboxylation is the more kinetically favoured on an overall basis compared to oxygenation.