Computational Models of Claudin Assembly in Tight Junctions and Strand Properties.

Claudins are one of the major components of tight junctions (TJs) that polymerize within the cell membrane and form interactions between cells. Some claudins seal the paracellular space, limiting paracellular flux, while others form selectively permeable ion channels that control the paracellular permeability of small ions. Claudin strands are known to be dynamic and reshape within TJs to accommodate large-scale movements and rearrangements of epithelial tissues. Here, we summarize the recent computational and modeling studies on claudin assembly into tetrameric ion channels and their polymerization into µm long strands within the membrane. Computational studies ranging from all-atom molecular dynamics, coarse-grained simulations, and hybrid-resolution simulations elucidate the molecular nature of claudin assembly and function and provide a framework that describes the lateral flexibility of claudin strands.

Phase Separation and Ion Diffusion in Ionic Liquid, Organic Solvent, and Lithium Salt Electrolyte Mixtures.

The highly desirable characteristics of ternary mixtures of ionic liquids, organic solvents, and metal salts make them a promising candidate for use in various electrothermal energy storage and conversion systems. In this study, using large-scale classical molecular dynamics simulations, we looked into 10 different ternary electrolyte mixtures using combinations of [EMIM]+, [BMIM]+, and [OMIM]+ cations with [NO3]-, [BF4]-, [PF6]-, [ClO4]-, [TFO]-, and [NTf2]- anions, tetraglyme, and Li salt to study the effect of ionic liquid composition on the phase behavior of ternary electrolyte mixtures. We uncovered that in these electrolytes, phase separation is mainly a function of pairwise binding energy of the constituents of the mixture. To corroborate this theory, several simulations are performed at various temperatures ranging from 260 to 500 K for each mixture, followed by calculating the binding energy of ionic liquid pairs using density functional theory. Our results verify that the transition temperature for the phase separation of each system is indeed a function of the pairwise binding energy of its ionic liquid pairs. It is also found that in some cases, the diffusion coefficient of the Li+ ions decreased even with the increase in the temperature, an effect that is attributed to the presence of condensed ionic domains in the electrolyte. This study provides a new insight for the design of multicomponent electrolyte mixtures for a wide range of energy applications.

Thermodynamics and Kinetics in Anisotropic Growth of One-Dimensional Midentropy Nanoribbons.

One-dimensional (1D) materials demonstrate anisotropic in-plane physical properties that enable a wide range of functionalities in electronics, photonics, valleytronics, optoelectronics, and catalysis. Here, we undertake an in-depth study of the growth mechanism for equimolar midentropy alloy of (NbTaTi)0.33S3 nanoribbons as a model system for 1D transition metal trichalcogenide structures. To understand the thermodynamic and kinetic effects in the growth process, the energetically preferred phases at different synthesis temperatures and times are investigated, and the phase evolution is inspected at a sequence of growth steps. It is uncovered that the dynamics of the growth process occurs at four different stages via preferential incorporation of chemical species at high-surface-energy facets. Also, a sequence of temperature and time dependent nonuniform to uniform phase evolutions has emerged in the composition and structure of (NbTaTi)0.33S3 which is described based on an anisotropic vapor-solid (V-S) mechanism. Furthermore, direct evidence for the 3D structure of the charge density wave (CDW) phase (width less than 100 nm) is provided by three-dimensional electron diffraction (3DED) in individual nanoribbons at cryogenic temperature, and detailed comparisons are made between the phases obtained before and after CDW transformation. This study provides important fundamental information for the design and synthesis of future 1D alloy structures.

Molecular mechanism of claudin-15 strand flexibility: A computational study.

Claudins are one of the major components of tight junctions that play a key role in the formation and maintenance of the epithelial barrier function. Tight junction strands are dynamic and capable of adapting their structure in response to large-scale tissue rearrangement and cellular movement. Here, we present molecular dynamics simulations of claudin-15 strands of up to 225 nm in length in two parallel lipid membranes and characterize their mechanical properties. The persistence length of claudin-15 strands is comparable with those obtained from analyses of freeze-fracture electron microscopy. Our results indicate that lateral flexibility of claudin strands is due to an interplay of three sets of interfacial interaction networks between two antiparallel double rows of claudins in the membranes. In this model, claudins are assembled into interlocking tetrameric ion channels along the strand that slide with respect to each other as the strands curve over submicrometer-length scales. These results suggest a novel molecular mechanism underlying claudin-15 strand flexibility. It also sheds light on intermolecular interactions and their role in maintaining epithelial barrier function.

Role of TM3 in claudin-15 strand flexibility: A molecular dynamics study.

Claudins are cell-cell adhesion proteins within tight junctions that connect epithelial cells together. Claudins polymerize into a network of strand-like structures within the membrane of adjoining cells and create ion channels that control paracellular permeability to water and small molecules. Tight junction morphology and barrier function is tissue specific and regulated by claudin subtypes. Here, we present a molecular dynamics study of claudin-15 strands within lipid membranes and the role of a single-point mutation (A134P) on the third transmembrane helix (TM3) of claudin-15 in determining the morphology of the strand. Our results indicate that the A134P mutation significantly affects the lateral flexibility of the strands, increasing the persistence length of claudin-15 strands by a factor of three. Analyses of claudin-claudin contact in our µsecond-long trajectories show that the mutation does not alter the intermolecular contacts (interfaces) between claudins. However, the dynamics and frequency of interfacial contacts are significantly affected. The A134P mutation introduces a kink in TM3 of claudin-15 similar to the one observed in claudin-3 crystal structure. The kink on TM3 skews the rotational flexibility of the claudins in the strands and limits their fluctuation in one direction. This asymmetric movement in the context of the double rows reduces the lateral flexibility of the strand and leads to higher persistence lengths of the mutant.

Multicomponent Phase Separation in Ternary Mixture Ionic Liquid Electrolytes.

We investigate the phase behavior of ternary mixtures of ionic liquid, organic solvent, and lithium salt by molecular dynamics simulations. We find that at room temperature, the electrolyte separates into distinct phases with specific compositions; an ion-rich domain that contains a fraction of solvent molecules and a second domain of pure solvent. The phase separation is shown to be entropy-driven and is independent of lithium salt concentration. Phase separation is only observed at microsecond time scales and greatly affects the transport properties of the electrolyte.

Computational Modeling of Claudin Structure and Function.

Tight junctions form a barrier to control passive transport of ions and small molecules across epithelia and endothelia. In addition to forming a barrier, some of claudins control transport properties of tight junctions by forming charge- and size-selective ion channels. It has been suggested claudin monomers can form or incorporate into tight junction strands to form channels. Resolving the crystallographic structure of several claudins in recent years has provided an opportunity to examine structural basis of claudins in tight junctions. Computational and theoretical modeling relying on atomic description of the pore have contributed significantly to our understanding of claudin pores and paracellular transport. In this paper, we review recent computational and mathematical modeling of claudin barrier function. We focus on dynamic modeling of global epithelial barrier function as a function of claudin pores and molecular dynamics studies of claudins leading to a functional model of claudin channels.

New Class of Electrocatalysts Based on 2D Transition Metal Dichalcogenides in Ionic Liquid.

The optimization of traditional electrocatalysts has reached a point where progress is impeded by fundamental physical factors including inherent scaling relations among thermokinetic characteristics of different elementary reaction steps, non-Nernstian behavior, and electronic structure of the catalyst. This indicates that the currently utilized classes of electrocatalysts may not be adequate for future needs. This study reports on synthesis and characterization of a new class of materials based on 2D transition metal dichalcogenides including sulfides, selenides, and tellurides of group V and VI transition metals that exhibit excellent catalytic performance for both oxygen reduction and evolution reactions in an aprotic medium with Li salts. The reaction rates are much higher for these materials than previously reported catalysts for these reactions. The reasons for the high activity are found to be the metal edges with adiabatic electron transfer capability and a cocatalyst effect involving an ionic-liquid electrolyte. These new materials are expected to have high activity for other core electrocatalytic reactions and open the way for advances in energy storage and catalysis.

Molecular determination of claudin-15 organization and channel selectivity.

Tight junctions are macromolecular structures that traverse the space between adjacent cells in epithelia and endothelia. Members of the claudin family are known to determine tight junction permeability in a charge- and size-selective manner. Here, we use molecular dynamics simulations to build and refine an atomic model of claudin-15 channels and study its transport properties. Our simulations indicate that claudin-15 forms well-defined channels for ions and molecules and otherwise "seals" the paracellular space through hydrophobic interactions. Ionic currents, calculated from simulation trajectories of wild-type as well as mutant channels, reflect in vitro measurements. The simulations suggest that the selectivity filter is formed by a cage of four aspartic acid residues (D55), contributed by four claudin-15 molecules, which creates a negative electrostatic potential to favor cation flux over anion flux. Charge reversal or charge ablation mutations of D55 significantly reduce cation permeability in silico and in vitro, whereas mutations of other negatively charged pore amino acid residues have a significantly smaller impact on channel permeability and selectivity. The simulations also indicate that water and small ions can pass through the channel, but larger cations, such as tetramethylammonium, do not traverse the pore. Thus, our model provides an atomic view of claudin channels, their transport function, and a potential three-dimensional organization of its selectivity filter.

A lithium-oxygen battery with a long cycle life in an air-like atmosphere.

Lithium-air batteries are considered to be a potential alternative to lithium-ion batteries for transportation applications, owing to their high theoretical specific energy. So far, however, such systems have been largely restricted to pure oxygen environments (lithium-oxygen batteries) and have a limited cycle life owing to side reactions involving the cathode, anode and electrolyte. In the presence of nitrogen, carbon dioxide and water vapour, these side reactions can become even more complex. Moreover, because of the need to store oxygen, the volumetric energy densities of lithium-oxygen systems may be too small for practical applications. Here we report a system comprising a lithium carbonate-based protected anode, a molybdenum disulfide cathode and an ionic liquid/dimethyl sulfoxide electrolyte that operates as a lithium-air battery in a simulated air atmosphere with a long cycle life of up to 700 cycles. We perform computational studies to provide insight into the operation of the system in this environment. This demonstration of a lithium-oxygen battery with a long cycle life in an air-like atmosphere is an important step towards the development of this field beyond lithium-ion technology, with a possibility to obtain much higher specific energy densities than for conventional lithium-ion batteries.

Modeling and Analysis of Intercalant Effects on Circular DNA Conformation.

The large-scale conformation of DNA molecules plays a critical role in many basic elements of cellular functionality and viability. By targeting the structural properties of DNA, many cancer drugs, such as anthracyclines, effectively inhibit tumor growth but can also produce dangerous side effects. To enhance the development of innovative medications, rapid screening of structural changes to DNA can provide important insight into their mechanism of interaction. In this study, we report changes to circular DNA conformation from intercalation with ethidium bromide using all-atom molecular dynamics simulations and characterized experimentally by translocation through a silicon nitride solid-state nanopore. Our measurements reveal three distinct current blockade levels and a 6-fold increase in translocation times for ethidium bromide-treated circular DNA as compared to untreated circular DNA. We attribute these increases to changes in the supercoiled configuration hypothesized to be branched or looped structures formed in the circular DNA molecule. Further evidence of the conformational changes is demonstrated by qualitative atomic force microscopy analysis. These results expand the current methodology for predicting and characterizing DNA tertiary structure and advance nanopore technology as a platform for deciphering structural changes of other important biomolecules.

Stable and Selective Humidity Sensing Using Stacked Black Phosphorus Flakes.

Black phosphorus (BP) atomic layers are known to undergo chemical degradation in humid air. Yet in more robust configurations such as films, composites, and embedded structures, BP can potentially be utilized in a large number of practical applications. In this study, we explored the sensing characteristics of BP films and observed an ultrasensitive and selective response toward humid air with a trace-level detection capability and a very minor drift over time. Our experiments show that the drain current of the BP sensor increases by ∼4 orders of magnitude as the relative humidity (RH) varies from 10% to 85%, which ranks it among the highest ever reported values for humidity detection. The mechanistic studies indicate that the operation principle of the BP film sensors is based on the modulation in the leakage ionic current caused by autoionization of water molecules and ionic solvation of the phosphorus oxoacids produced on moist BP surfaces. Our stability tests reveal that the response of the BP film sensors remains nearly unchanged after prolonged exposures (up to 3 months) to ambient conditions. This study opens up the route for utilizing BP stacked films in many potential applications such as energy generation/storage systems, electrocatalysis, and chemical/biosensing.

Bimodal Phonon Scattering in Graphene Grain Boundaries.

Graphene has served as the model 2D system for over a decade, and the effects of grain boundaries (GBs) on its electrical and mechanical properties are very well investigated. However, no direct measurement of the correlation between thermal transport and graphene GBs has been reported. Here, we report a simultaneous comparison of thermal transport in supported single crystalline graphene to thermal transport across an individual graphene GB. Our experiments show that thermal conductance (per unit area) through an isolated GB can be up to an order of magnitude lower than the theoretically anticipated values. Our measurements are supported by Boltzmann transport modeling which uncovers a new bimodal phonon scattering phenomenon initiated by the GB structure. In this novel scattering mechanism, boundary roughness scattering dominates the phonon transport in low-mismatch GBs, while for higher mismatch angles there is an additional resistance caused by the formation of a disordered region at the GB. Nonequilibrium molecular dynamics simulations verify that the amount of disorder in the GB region is the determining factor in impeding thermal transport across GBs.

A structural rearrangement of the Na+/K+-ATPase traps ouabain within the external ion permeation pathway.

With the use of the energy of ATP hydrolysis, the Na+/K+-ATPase is able to transport across the cell membrane Na+ and K+ against their electrochemical gradients. The enzyme is strongly inhibited by ouabain and its derivatives, some that are therapeutically used for patients with heart failure (cardiotonic steroids). Using lanthanide resonance energy transfer, we trace here the conformational changes occurring on the external side of functional Na+/K+-ATPases induced by the binding of ouabain. Changes in donor/acceptor pair distances are mainly observed within the α subunit of the enzyme. To derive a structural model matching the experimental lanthanide resonance energy transfer distances measured with bound ouabain, we carried out molecular dynamics simulations with energy restraints applied simultaneously using a novel methodology with multiple non-interacting fragments. The restrained simulation, initiated from the X-ray structure of the E2(2K+) state, became strikingly similar to the X-ray structure of the sodium-bound state. The final model shows that ouabain is trapped within the external ion permeation pathway of the pump.

Molecular dynamics simulations of membrane proteins under asymmetric ionic concentrations.

A computational method is developed to allow molecular dynamics simulations of biomembrane systems under realistic ionic gradients and asymmetric salt concentrations while maintaining the conventional periodic boundary conditions required to minimize finite-size effects in an all-atom explicit solvent representation. The method, which consists of introducing a nonperiodic energy step acting on the ionic species at the edge of the simulation cell, is first tested with illustrative applications to a simple membrane slab model and a phospholipid membrane bilayer. The nonperiodic energy-step method is then used to calculate the reversal potential of the bacterial porin OmpF, a large cation-specific ß-barrel channel, by simulating the I-V curve under an asymmetric 10:1 KCl concentration gradient. The calculated reversal potential of 28.6 mV is found to be in excellent agreement with the values of 26-27 mV measured from lipid bilayer experiments, thereby demonstrating that the method allows realistic simulations of nonequilibrium membrane transport with quantitative accuracy. As a final example, the pore domain of Kv1.2, a highly selective voltage-activated K(+) channel, is simulated in a lipid bilayer under conditions that recreate, for the first time, the physiological K(+) and Na(+) concentration gradients and the electrostatic potential difference of living cells.

An emerging consensus on voltage-dependent gating from computational modeling and molecular dynamics simulations.

Developing an understanding of the mechanism of voltage-gated ion channels in molecular terms requires knowledge of the structure of the active and resting conformations. Although the active-state conformation is known from x-ray structures, an atomic resolution structure of a voltage-dependent ion channel in the resting state is not currently available. This has motivated various efforts at using computational modeling methods and molecular dynamics (MD) simulations to provide the missing information. A comparison of recent computational results reveals an emerging consensus on voltage-dependent gating from computational modeling and MD simulations. This progress is highlighted in the broad context of preexisting work about voltage-gated channels.

Molecular dynamics investigation of the ω-current in the Kv1.2 voltage sensor domains.

Voltage sensor domains (VSD) are transmembrane proteins that respond to changes in membrane voltage and modulate the activity of ion channels, enzymes, or in the case of proton channels allow permeation of protons across the cell membrane. VSDs consist of four transmembrane segments, S1-S4, forming an antiparallel helical bundle. The S4 segment contains several positively charged residues, mainly arginines, located at every third position along the helix. In the voltage-gated Shaker K(+) channel, the mutation of the first arginine of S4 to a smaller uncharged amino acid allows permeation of cations through the VSD. These currents, known as ω-currents, pass through the VSD and are distinct from K(+) currents passing through the main ion conduction pore. Here we report molecular dynamics simulations of the ω-current in the resting-state conformation for Kv1.2 and for four of its mutants. The four tested mutants exhibit various degrees of conductivity for K(+) and Cl(-) ions, with a slight selectivity for K(+) over Cl(-). Analysis of the ion permeation pathway, in the case of a highly conductive mutant, reveals a negatively charged constriction region near the center of the membrane that might act as a selectivity filter to prevent permeation of anions through the pore. The residues R1 in S4 and E1 in S2 are located at the narrowest region of the ω-pore for the resting state conformation of the VSD, in agreement with experiments showing that the largest increase in current is produced by the double mutation E1D and R1S.

Constant electric field simulations of the membrane potential illustrated with simple systems.

Advances in modern computational methods and technology make it possible to carry out extensive molecular dynamics simulations of complex membrane proteins based on detailed atomic models. The ultimate goal of such detailed simulations is to produce trajectories in which the behavior of the system is as realistic as possible. A critical aspect that requires consideration in the case of biological membrane systems is the existence of a net electric potential difference across the membrane. For meaningful computations, it is important to have well validated methodologies for incorporating the latter in molecular dynamics simulations. A widely used treatment of the membrane potential in molecular dynamics consists of applying an external uniform electric field E perpendicular to the membrane. The field acts on all charged particles throughout the simulated system, and the resulting applied membrane potential V is equal to the applied electric field times the length of the periodic cell in the direction perpendicular to the membrane. A series of test simulations based on simple membrane-slab models are carried out to clarify the consequences of the applied field. These illustrative tests demonstrate that the constant-field method is a simple and valid approach for accounting for the membrane potential in molecular dynamics studies of biomolecular systems. This article is part of a Special Issue entitled: Membrane protein structure and function.

Ouabain binding site in a functioning Na+/K+ ATPase.

The Na(+)/K(+) ATPase is an almost ubiquitous integral membrane protein within the animal kingdom. It is also the selective target for cardiotonic derivatives, widely prescribed inhibitors for patients with heart failure. Functional studies revealed that ouabain-sensitive residues distributed widely throughout the primary sequence of the protein. Recently, structural work has brought some consensus to the functional observations. Here, we use a spectroscopic approach to estimate distances between a fluorescent ouabain and a lanthanide binding tag (LBT), which was introduced at five different positions in the Na(+)/K(+) ATPase sequence. These five normally functional LBT-Na(+)/K(+) ATPase constructs were expressed in the cell membrane of Xenopus laevis oocytes, operating under physiological internal and external ion conditions. The spectroscopic data suggest two mutually exclusive distances between the LBT and the fluorescent ouabain. From the estimated distances and using homology models of the LBT-Na(+)/K(+) ATPase constructs, approximate ouabain positions could be determined. Our results suggest that ouabain binds at two sites along the ion permeation pathway of the Na(+)/K(+) ATPase. The external site (low apparent affinity) occupies the same region as previous structural findings. The high apparent affinity site is, however, slightly deeper toward the intracellular end of the protein. Interestingly, in both cases the lactone ring faces outward. We propose a sequential ouabain binding mechanism that is consistent with all functional and structural studies.

On the sensing mechanism in carbon nanotube chemiresistors.

There has been recent controversy whether the response seen in carbon nanotube (CNT) chemiresistors is associated with a change in the resistance of the individual nanotubes or changes in the resistance of the junctions. In this study, we carry out a network analysis to understand the relative contributions of the nanotubes and the junctions to the change in resistance of the nanotube network. We find that the dominant mode of detection in nanotube networks changes according to the conductance level (defect level) in the nanotubes. In networks with perfect nanotubes, changes in the junctions between adjacent nanotubes and junctions between the contacts and the CNTs can cause a detectable change in the resistance of the nanotube networks, while adsorption on the nanotubes has a smaller effect. In contrast, in networks with highly defective nanotubes, the changes in the resistance of the individual nanotubes cause a detectable change in the overall resistance of a chemiresistor network, while changes in the junctions have smaller effects. The combinational effect is also observed for the case in between. The results show that the sensing mechanism of a nanotube network can change according to the defect levels of the nanotubes, which may explain the apparently contradictory results in the literature.