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.

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.

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.

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.

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.

Spider Toxin SNX-482 Gating Modifier Spontaneously Partitions in the Membrane Guided by Electrostatic Interactions.

Spider toxin SNX-482 is a cysteine-rich peptide that interferes with calcium channel activity by binding to voltage-sensing domains of the CaV2.3 subtype. Two mechanisms dominate the binding process of cysteine-rich peptides: direct binding from the aqueous phase or through lateral diffusion from the membrane, the so-called reduction in dimensionality mechanism. In this work, via coarse-grained and atomistic molecular dynamics simulations, we have systematically studied the spontaneous partitioning of SNX-482 with membranes of different anionic compositions and explored via diffusional analysis both binding mechanisms. Our simulations revealed a conserved protein patch that inserts in the membrane, a preference for binding towards partially negatively charged membranes, and that electrostatics guides membrane binding by incrementing and aligning the molecular dipole. Finally, diffusivity calculations showed that the toxin diffusion along the membrane plane is an order of magnitude slower than the aqueous phase suggesting that the critical factor in determining the SNX-482-CaV2.3 binding mechanism is the affinity between the membrane and SNX-482.

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.

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.

Simulations on Simple Models of Connexin Hemichannels Indicate That Ca2+ Blocking Is Not a Pure Electrostatic Effect.

Connexin hemichannels allow the unspecific but regulated interchange of molecules from ions to second messenger and ATP, between the eukariotic cell and its extracellular space. The transport of ions and water through hemichannels is important for physiological functions and also in the progression of several pathological conditions. Extracellular Ca2+ concentration is one of the regulators that drives the channel to a closed state. However the relation between their functional and structural states is far for being totally understood. In this work, we modelled connexin hemichannels using simple systems based on a fixed array of carbon atoms and assess the Ca2+ regulation using molecular dynamics simulations. The two proposed mechanism described so far for calcium action were studied combined, e.g., an electrostatic effect and a pore stretching. Our results show that the addition of positive charge density inside the channel cannot stop the flow of potassium, chloride nor water. Only a pore stretching at the center of the pore can explain the channel blocking.

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.

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.

Orientational and Folding Thermodynamics via Electric Dipole Moment Restraining.

The projection of molecular processes onto a small set of relevant descriptors, the so-called reaction coordinates or collective variables (CVs), is a technique nowadays routinely employed by the biomolecular simulation community. In this work, we implemented two CVs to manipulate the orientation (i.e., angle) (µâƒ—a) and magnitude (|µâƒ—|) of the electric dipole moment. In doing so, we studied the thermodynamics of water orientation under the application of external voltages and the folding of two polypeptides at zero-field conditions. The projection of the free-energy [potential of mean force (PMF)] along water orientation defined an upper limit of around 0.3 V for irrelevant thermodynamic effects. On the other hand, sufficiently strong µâƒ—a restraints applied on 12-alanine (Ala12) triggered structural effects because of the alignment of local dipoles; for lower restraints, a full-body rotation is achieved. The manipulation of |µâƒ—| produced strong perturbations on the secondary structure of Ala12, promoting an enhanced sampling to its configurational space. Rigorous free-energy calculations in the form of 2-D PMFs for deca-alanine showed the utility of |µâƒ—| as a reaction coordinate to study folding in small α helices. As a whole, we propose that the manipulation of both components of the dipole moment, µâƒ—a and |µâƒ—|, provides thermodynamics insights into the structural conformation and stability of biomolecules. These new CVs are implemented in the Colvars module, available for NAMD and LAMMPS.

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.

Transprotein-Electropore Characterization: A Molecular Dynamics Investigation on Human AQP4.

Electroporation characterization is a topic of intensive interest probed by extensive ongoing research efforts. Usually, these studies are carried out on lipid-bilayer electroporation. Surprisingly, the possibility of water-channel electropore formation across transmembrane proteins themselves, particularly in view of such a promising application, has not yet been elucidated. The present work examines the geometrical and kinetic aspects of electropores and their stability in such a protein milieux (as opposed through the phospholipid membranes) in depth, by means of scrutiny of such a process in human-AQP4 as a well-representative prototype. The residues forming the electropore's walls, organized in loops, reveal the formation mechanism by their dipole alignment and translational response in response to applied axial electric fields in nonequilibrium molecular dynamics simulation. The magnitude of sustaining electric fields (keeping a stable electropore open) were determined. This suggests that transmembrane proteins could play a central role in electroporation applications, e.g., in medicine and biotechnology.

RIP-MD: a tool to study residue interaction networks in protein molecular dynamics.

Protein structure is not static; residues undergo conformational rearrangements and, in doing so, create, stabilize or break non-covalent interactions. Molecular dynamics (MD) is a technique used to simulate these movements with atomic resolution. However, given the data-intensive nature of the technique, gathering relevant information from MD simulations is a complex and time consuming process requiring several computational tools to perform these analyses. Among different approaches, the study of residue interaction networks (RINs) has proven to facilitate the study of protein structures. In a RIN, nodes represent amino-acid residues and the connections between them depict non-covalent interactions. Here, we describe residue interaction networks in protein molecular dynamics (RIP-MD), a visual molecular dynamics (VMD) plugin to facilitate the study of RINs using trajectories obtained from MD simulations of proteins. Our software generates RINs from MD trajectory files. The non-covalent interactions defined by RIP-MD include H-bonds, salt bridges, VdWs, cation-π, π-π, Arginine-Arginine, and Coulomb interactions. In addition, RIP-MD also computes interactions based on distances between Cαs and disulfide bridges. The results of the analysis are shown in an user friendly interface. Moreover, the user can take advantage of the VMD visualization capacities, whereby through some effortless steps, it is possible to select and visualize interactions described for a single, several or all residues in a MD trajectory. Network and descriptive table files are also generated, allowing their further study in other specialized platforms. Our method was written in python in a parallelized fashion. This characteristic allows the analysis of large systems impossible to handle otherwise. RIP-MD is available at http://www.dlab.cl/ripmd.

Disaccharide-Based Anionic Amphiphiles as Potent Inhibitors of Lipopolysaccharide-Induced Inflammation.

Despite significant advances made in the last decade in the understanding of molecular mechanisms of sepsis and in the development of clinically relevant therapies, sepsis remains the leading cause of mortality in intensive care units with increasing incidence worldwide. Toll-like receptor 4 (TLR4)-a transmembrane pattern-recognition receptor responsible for propagating the immediate immune response to Gram-negative bacterial infection-plays a central role in the pathogenesis of sepsis and chronic inflammation-related disorders. TLR4 is complexed with the lipopolysaccharide (LPS)-sensing protein myeloid differentiation-2 (MD-2) which represents a preferred target for establishing new anti-inflammatory treatment strategies. Herein we report the development, facile synthesis, and biological evaluation of novel disaccharide-based TLR4⋅MD-2 antagonists with potent anti-endotoxic activity at micromolar concentrations. A series of synthetic anionic glycolipids entailing amide-linked ß-ketoacyl lipid residues was prepared in a straightforward manner by using a single orthogonally protected nonreducing diglucosamine scaffold. Suppression of the LPS-induced release of interleukin-6 and tumor necrosis factor was monitored and confirmed in human immune cells (MNC and THP1) and mouse macrophages. Structure-activity relationship studies and molecular dynamics simulations revealed the structural basis for the high-affinity interaction between anionic glycolipids and MD-2, and highlighted two compounds as leads for the development of potential anti-inflammatory therapeutics.

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.

Exploration of graphene oxide as an intelligent platform for cancer vaccines.

We explored an intelligent vaccine system via facile approaches using both experimental and theoretical techniques based on the two-dimensional graphene oxide (GO). Without extra addition of bio/chemical stimulators, the microsized GO imparted various immune activation tactics to improve the antigen immunogenicity. A high antigen adsorption was acquired, and the mechanism was revealed to be a combination of electrostatic, hydrophobic, and π-π stacking interactions. The "folding GO" acted as a cytokine self-producer and antigen reservoir and showed a particular autophagy, which efficiently promoted the activation of antigen presenting cells (APCs) and subsequent antigen cross-presentation. Such a "One but All" modality thus induced a high level of anti-tumor responses in a programmable way and resulted in efficient tumor regression in vivo. This work may shed light on the potential use of a new dimensional nano-platform in the development of high-performance cancer vaccines.

Protein corona mitigates the cytotoxicity of graphene oxide by reducing its physical interaction with cell membrane.

Many recent studies have shown that the way nanoparticles interact with cells and biological molecules can vary greatly in the serum-containing or serum-free culture medium. However, the underlying molecular mechanisms of how the so-called "protein corona" formed in serum medium affects nanoparticles' biological responses are still largely unresolved. Thus, it is critical to understand how absorbed proteins on the surfaces of nanoparticles alter their biological effects. In this work, we have demonstrated with both experimental and theoretical approaches that protein BSA coating can mitigate the cytotoxicity of graphene oxide (GO) by reducing its cell membrane penetration. Our cell viability and cellular uptake experiments showed that protein corona decreased cellular uptake of GO, thus significantly mitigating the potential cytotoxicity of GO. The electron microscopy images also confirmed that protein corona reduced the cellular morphological damage by limiting GO penetration into the cell membrane. Further molecular dynamics (MD) simulations validated the experimental results and revealed that the adsorbed BSA in effect weakened the interaction between the phospholipids and graphene surface due to a reduction of the available surface area plus an unfavorable steric effect, thus significantly reducing the graphene penetration and lipid bilayer damaging. These findings provide new insights into the underlying molecular mechanism of this important graphene protein corona interaction with cell membranes, and should have implications in future development of graphene-based biomedical applications.