Solvation Structures and Transport Mechanisms of Cations in Water-in-Bisalt Electrolytes for High-Concentration Lithium-Ion Batteries Revealed by Molecular Dynamics Simulations.

The development of safe and cost-effective electrolytes for rechargeable batteries is currently underway. While water-based electrolytes hold promise, their restricted electrochemical stability window poses a challenge. Combining multiple ionic species emerges as a promising strategy to broaden this stability window and optimize Li-ion battery performance. This study focuses on dual-cation electrolytes, which blend lithium and potassium acetates to enhance the electrochemical characteristics of the solution at high concentrations. We investigated the solvation structure of each ion and its interactions on a molecular level. Our analysis reveals that ion clusters and aggregates are formed through shared acetate and water molecules at high salt concentrations. Furthermore, the residence time analyses of atom pairs indicate that cations diffuse in vehicular mode at low concentrations. In contrast, they switch to a structural mode at high concentrations due to diminishing water content. This study offers a comprehensive model for exploring diverse solvation structures of cations and gaining insights into their diffusion mechanisms within water-in-bisalt electrolytes for aqueous Li-ion batteries.

Dielectric Response and Linear Absorption Spectroscopy of Ionic Systems.

Time-dependent electric fields applied to ionic systems can induce both a dielectric and a conductive response, leading to the generation of macroscopic polarization and current, respectively. It has long been recognized that it is not possible to determine the two types of responses separately. However, this aspect is often not adequately accounted for in dielectric and absorption spectroscopies of ionic systems. To clarify this, we theoretically investigate the dielectric and conductive responses of ionic systems containing polyatomic ions based on linear response theory. We derive general expressions for the frequency-dependent dielectric functions, conductivity, and absorption coefficient, including those measured experimentally. Furthermore, we show that the dielectric and conductive responses cannot be uniquely distinguished even at the theoretical level and, therefore, cannot represent experimentally measured quantities. Instead, dielectric and absorption spectra of ionic systems should be expressed in terms of the generalized dielectric function that encompasses both dielectric and conductive responses. We propose a computational method to calculate this generalized dielectric function reliably. Model calculations on concentrated aqueous solutions of NaCl, a monatomic salt, and LiTFSI, a polyatomic salt, show that the dielectric and linear absorption spectra of the two systems based on the generalized dielectric function are significantly different from purely dielectric counterparts in the far-IR, terahertz, and lower-frequency regions. Moreover, the spectra are mainly determined by the autocorrelations of total dipole and total current, but dipole-current cross-correlation can also significantly contribute to the spectra of the LiTFSI solution. The present theoretical approach could be extended to nonlinear spectroscopy of ionic liquids and electrolyte solutions.

Revisiting Ultrafast Dynamics in Carbonate-Based Electrolytes for Li-Ion Batteries: Clarifying 2D-IR Cross-Peak Interpretation.

Understanding chemical exchange in carbonate-based electrolytes employed in Li-ion batteries (LIBs) is crucial for elucidating ion transport mechanisms. Ultrafast two-dimensional (2D) IR spectroscopy has been widely used to investigate the solvation structure and dynamics of Li-ions in organic carbonate-based electrolytes. However, the interpretation of cross-peaks observed in picosecond carbonyl stretch 2D-IR spectra has remained contentious. These cross-peaks could arise from various phenomena, including vibrational couplings between neighboring carbonyl groups in the first solvation shell around Li-ions, vibrational excitation transfers between carbonyl groups in distinct solvation environments, and local heating effects. Therefore, it is imperative to resolve the interpretation of 2D-IR cross-peaks to avoid misinterpretations regarding ultrafast dynamics found in LIB carbonate-based electrolytes. In this study, we have taken a comprehensive investigation of carbonate-based electrolytes utilizing 2D-IR spectroscopy and molecular dynamics (MD) simulations. Through meticulous analyses and interpretations, we have identified that the cross-peaks observed in the picosecond 2D-IR spectra of LIB electrolytes predominantly arise from intermolecular vibrational excitation transfer processes between the carbonyl groups of Li-bound and free carbonate molecules. We further discuss the limitations of employing a picosecond 2D-IR spectroscopic technique to study chemical exchange and intermolecular vibrational excitation transfer processes, particularly when the effects of the molecular photothermal process cannot be ignored. Our findings shed light on the dynamics of LIB electrolytes and resolve the controversy related to 2D-IR cross-peaks. By discerning the origin of these features, we could provide valuable insights for the design and optimization of next-generation Li-ion batteries.

Lithium-Ion Solvation Structure in Organic Carbonate Electrolytes at Low Temperatures.

Lithium-ion batteries face insufficient capacity at low temperatures. The lithium-ion desolvation process in the vicinity of a solid electrolyte interphase (SEI) layer is considered the major problem. Thus, an accurate determination of lithium-ion solvation structures is a prerequisite for understanding this process. Here, using a cryostat combined with an FTIR spectrometer, we found that as the temperature decreased, the number of coordinating carbonates in the first solvation shell of the lithium ion increased with a decreased population of the contact ion pair (CIP). More specifically, we found that two or more carbonate molecules replace a single PF6- anion upon CIP dissociation. This experimental finding shows that the prevailing notion that four carbonate molecules coordinate each lithium ion to form a tetrahedral structure is invalid for describing lithium-ion solvation structures. We anticipate that the present work will elucidate one of the molecular origins behind the low performance of lithium-ion batteries at low temperatures.

Two-dimensional electronic spectroscopy of bacteriochlorophyll a with synchronized dual mode-locked lasers.

How atoms and electrons in a molecule move during a chemical reaction and how rapidly energy is transferred to or from the surroundings can be studied with flashes of laser light. However, despite prolonged efforts to develop various coherent spectroscopic techniques, the lack of an all-encompassing method capable of both femtosecond time resolution and nanosecond relaxation measurement has hampered various applications of studying correlated electron dynamics and vibrational coherences in functional materials and biological systems. Here, we demonstrate that two broadband (>300 nm) synchronized mode-locked lasers enable two-dimensional electronic spectroscopy (2DES) study of chromophores such as bacteriochlorophyll a in condensed phases to measure both high-resolution coherent vibrational spectrum and nanosecond electronic relaxation. We thus anticipate that the dual mode-locked laser-based 2DES developed and demonstrated here would be of use for unveiling the correlation between the quantum coherence and exciton dynamics in light-harvesting protein complexes and semiconducting materials.

Water hydrogen-bonding structure and dynamics near lipid multibilayer surface: Molecular dynamics simulation study with direct experimental comparison.

Lipid multibilayers are excellent model systems for studying water structures and dynamics near biological membrane surfaces. In particular, the orientational distribution and rotational dynamics of water molecules near hydrophilic lipid groups are found to be sensitive to the chemical nature and charge distributions of the amphiphilic lipids. To elucidate how different parts of these lipids affect the water hydrogen-bonding structure and dynamics and to directly compare with recent experimental results, we carried out molecular dynamics (MD) simulations of lipid multibilayer systems. We found that the water molecules close to positively charged choline groups have a broad distribution of orientations due to the clathratelike shell formation around the choline groups but that those associated with phosphate groups, even in the second hydration shell, are orientationally restricted due to their strong hydrogen bonding with the phosphate group. These MD simulation results are in excellent agreement with our time-resolved infrared pump-probe anisotropy measurements, and we believe that they provide valuable insights into the role of water molecules in maintaining lipid bilayer integrity.

Theory of coherent two-dimensional vibrational spectroscopy.

Two-dimensional (2D) vibrational spectroscopy has emerged as one of the most important experimental techniques useful to study the molecular structure and dynamics in condensed phases. Theory and computation have also played essential and integral roles in its development through the nonlinear optical response theory and computational methods such as molecular dynamics (MD) simulations and electronic structure calculations. In this article, we present the fundamental theory of coherent 2D vibrational spectroscopy and describe computational approaches to simulate the 2D vibrational spectra. The classical approximation to the quantum mechanical nonlinear response function is invoked from the outset. It is shown that the third-order response function can be evaluated in that classical limit by using equilibrium or non-equilibrium MD simulation trajectories. Another simulation method is based on the assumptions that the molecular vibrations can still be described quantum mechanically and that the relevant molecular response functions are evaluated by the numerical integration of the Schrödinger equation. A few application examples are presented to help the researchers in this and related areas to understand the fundamental principles and to use these methods for their studies with 2D vibrational spectroscopic techniques. In summary, this exposition provides an overview of current theoretical efforts to understand the 2D vibrational spectra and an outlook for future developments.

Ab initio Modeling of the Vibrational Sum-Frequency Generation Spectrum of Interfacial Water.

Understanding the structural and dynamical features of interfacial water is of greatest interest in physics, chemistry, biology, and materials science. Vibrational sum-frequency generation (SFG) spectroscopy, which is sensitive to the molecular orientation and dynamics on the surfaces or at the interfaces, allows one to study a wide variety of interfacial systems. The structural and dynamical features of interfacial water at the air/water interface have been extensively investigated by SFG spectroscopy. However, the interpretations of the spectroscopic features have been under intense debate. Here, we report a simulated SFG spectrum of the air/water interface based on ab initio molecular dynamics simulations, which covers the OH stretching, bending, and libration modes of interfacial water. Quantitative agreement between our present simulations and the most recent experimental studies ensures that ab initio simulations predict unbiased structural features and electrical properties of interfacial systems. By utilizing the kinetic energy spectral density (KESD) analysis to decompose the simulated spectra, the spectroscopic features can then be assigned to specific hydrogen-bonding configurations of interfacial water molecules.

Ion aggregation in high salt solutions. VII. The effect of cations on the structures of ion aggregates and water hydrogen-bonding network.

Ions in high salt solutions have a strong propensity to form polydisperse ion aggregates with broad size and shape distributions. In a series of previous comparative investigations using femtosecond IR pump-probe spectroscopy, molecular dynamics simulation, and graph theoretical analysis, we have shown that there exists a morphological difference in the structures of ion aggregates formed in various salt solutions. As salt concentration increases, the ions in high salt solutions form either cluster-like structures excluding water molecules or network-like structures entwined with water hydrogen-bonding networks. Interestingly, such morphological characteristics of the ion aggregates have been found to be in correlation with the solubility limits of salts. An important question that still remains unexplored is why certain salts with different cations have notably different solubility limits in water. Here, carrying out a series of molecular dynamics simulations of aqueous salt solutions and analyzing the distributions and connectivity patterns of ion aggregates with a spectral graph analysis method, we establish the relationship between the salt solubility and the ion aggregate morphology with a special emphasis on the cationic effects on water structures and ion aggregation. We anticipate that the understanding of large scale ion aggregate structures revealed in this study will be critical for elucidating the specific ion effects on the solubility and conformational stability of co-solute molecules such as proteins in water.

Hydrogen bonding and vibrational energy relaxation of interfacial water: A full DFT molecular dynamics simulation.

The air-water interface has been a subject of extensive theoretical and experimental studies due to its ubiquity in nature and its importance as a model system for aqueous hydrophobic interfaces. We report on the structure and vibrational energy transfer dynamics of this interfacial water system studied with equilibrium and non-equilibrium molecular dynamics simulations employing a density functional theory -based description of the system and the kinetic energy spectral density analysis. The interfacial water molecules are found to make fewer and weaker hydrogen (H)-bonds on average compared to those in the bulk. We also find that (i) the H-bonded OH groups conjugate to the free OH exhibit rather low vibrational frequencies (3000-3500 cm-1); (ii) the presence of a significant fraction (>10%) of free and randomly oriented water molecules at the interface ("labile water"), neither of whose OH groups are strong H-bond donors; (iii) the inertial rotation of free OH groups, especially from the labile water, contribute to the population decay of excited free OH groups with comparable rate and magnitude as intramolecular energy transfer between the OH groups. These results suggest that the labile water, which might not be easily detectable by the conventional vibrational sum frequency generation method, plays an important role in the surface water dynamics.

Simultaneous spectral and temporal analyses of kinetic energies in nonequilibrium systems: theory and application to vibrational relaxation of O-D stretch mode of HOD in water.

A time series of kinetic energies (KE) from classical molecular dynamics (MD) simulation contains fundamental information on system dynamics. It can also be analyzed in the frequency domain through Fourier transformation (FT) of velocity correlation functions, providing energy content of different spectral regions. By limiting the FT time span, we have previously shown that spectral resolution of KE evolution is possible in the nonequilibrium situations [Jeon and Cho, J. Chem. Phys. 2011, 135, 214504]. In this paper, we refine the method by employing the concept of instantaneous power spectra, extending it to reflect an instantaneous time-correlation of velocities with those in the future as well as with those in the past, and present a new method to obtain the instantaneous spectral density of KE (iKESD). This approach enables the simultaneous spectral and temporal resolution of KE with unlimited time precision. We discuss the formal and novel properties of the new iKESD approaches and how to optimize computational methods and determine parameters for practical applications. The method is specifically applied to the nonequilibrium MD simulation of vibrational relaxation of the OD stretch mode in a hydrated HOD molecule by employing a hybrid quantum mechanical/molecular mechanical (QM/MM) potential. We directly compare the computational results with the OD band population relaxation time profiles extracted from the IR pump-probe measurements for 5% HOD in water. The calculated iKESD yields the OD bond relaxation time scale ∼30% larger than the experimental value, and this decay is largely frequency-independent if the classical anharmonicity is accounted for. From the integrated iKESD over intra- and intermolecular bands, the major energy transfer pathways were found to involve the HOD bending mode in the subps range, then the internal modes of the solvent until 5 ps after excitation, and eventually the solvent intermolecular modes. Also, strong hydrogen-bonding of HOD is found to significantly hinder the initial intramolecular energy transfer process.

An accurate classical simulation of a two-dimensional vibrational spectrum: OD stretch spectrum of a hydrated HOD molecule.

An accurate computational method for the classical simulation of the two-dimensional vibrational spectra is presented. The method refines our previous computational method for the third order vibrational response function in the classical limit, and it enables capturing the diagonal elongation and its waiting time (T) dependence widely observed in experimental two-dimensional infrared (2D IR) spectra of intramolecular modes. The improvement is achieved by a series of new developments including (i) a block algorithm for the stability matrix computation, (ii) new equations of motion for the position-perturbed molecular dynamics (MD) trajectory, and (iii) enhanced sampling efficiency by exploiting the time-reversal invariance of MD trajectories. The method is applied to the simulation of 2D IR spectra of the OD stretch mode in a hydrated HOD molecule, employing a hybrid quantum mechanical/molecular mechanical force field. The simulated spectra exhibit diagonal elongation of the 2D IR signal at small T, reflecting the correlation of individual transitions among the inhomogeneously broadened ensemble. The slopes of the nodal lines of the elongated signals are found to decay with a time scale of 1.6 ps as T increases, in reasonable agreement with the frequency correlation decay time of 1.2 ps. The amplitudes of the positive and negative peaks also decay as T increases, due to vibrational population relaxation and molecular rotation. The peak positions tend to blue shift with increasing T, reflecting the different relaxation rates of the strongly and weakly solvated HOD species. These results indicate that the present method can reliably predict the waiting-time-dependent changes of 2D IR spectra of a single vibrational chromophore in solution.

Induced optical activity of DNA-templated cyanine dye aggregates: exciton coupling theory and TD-DFT studies.

Certain cyanine dye molecules have been observed to self-assemble in DNA templates to form large chiral aggregates, which exhibit induced circular dichroism. The structure and circular dichroism (CD) of one such system, aggregates of a cationic DiSC2(5) cyanine dye, are investigated using the time-dependent Kohn-Sham density functional theory (TD-DFT) and exciton coupling model. A series of TD-DFT calculations on the aggregates with one, two, and four dye molecules clearly shows the onset of CD induced by the helically twisted structure compatible with the minor groove of DNA templates. More simplified exciton coupling model analysis successfully reproduces the major positive Cotton effect observed in the experiment as well as the TD-DFT calculations, but it is unable to capture minor features of the CD spectrum that are closely related to absolute configurations of constituent dyes in the complex. We assess the performance of various methods used for evaluation of the electronic coupling energies between interacting chromophores. Our results confirm that the interchromophore interactions in cyanine dye aggregates are primarily electrostatic in nature and indicate that the exciton coupling model is adequate for studying induced CD of self-assembled aggregates of cyanine dye molecules.

Infrared probing of 4-azidoproline conformations modulated by azido configurations.

4-Azidoproline (Azp) can tune the stability of the polyproline II (P(II)) conformation in collagen. The azido group in the 4R and 4S configurations stabilizes and destabilizes the P(II) conformation, respectively. To obtain insights into the dependence of the conformational stability on the azido configuration, we carried out Fourier transform (FT) IR experiments with four 4-azidoproline derivatives, Ac-(4R/S)-Azp-(NH/O)Me. We found that the amide I and azido IR spectra are different depending on the azido configuration and C-terminal structure. The origin of such spectral differences between 4R and 4S configurations and between C-terminal methylamide and ester ends was elucidated by quantum chemistry calculations in combination with (1)H NMR and time- and frequency-resolved IR pump-probe spectroscopy. We found that the azido configurations and C-terminal structures affect intramolecular interactions, which are responsible for the ensuing conformational and thereby IR spectral differences. Consequently, 4-azidoproline conformations modulated by azido configurations can be probed by IR spectroscopy. These findings suggest that 4-azidoproline can be both a structure-control and -probing element, which enables the infrared tracking of proline roles in protein structure, function, and dynamics.

Vibrational spectroscopic determination of local solvent electric field, solute-solvent electrostatic interaction energy, and their fluctuation amplitudes.

IR probes have been extensively used to monitor local electrostatic and solvation dynamics. Particularly, their vibrational frequencies are highly sensitive to local solvent electric field around an IR probe. Here, we show that the experimentally measured vibrational frequency shifts can be inversely used to determine local electric potential distribution and solute-solvent electrostatic interaction energy. In addition, the upper limits of their fluctuation amplitudes are estimated by using the vibrational bandwidths. Applying this method to fully deuterated N-methylacetamide (NMA) in D(2)O and examining the solvatochromic effects on the amide I' and II' mode frequencies, we found that the solvent electric potential difference between O(═C) and D(-N) atoms of the peptide bond is about 5.4 V, and thus, the approximate solvent electric field produced by surrounding water molecules on the NMA is 172 MV/cm on average if the molecular geometry is taken into account. The solute-solvent electrostatic interaction energy is estimated to be -137 kJ/mol, by considering electric dipole-electric field interaction. Furthermore, their root-mean-square fluctuation amplitudes are as large as 1.6 V, 52 MV/cm, and 41 kJ/mol, respectively. We found that the water electric potential on a peptide bond is spatially nonhomogeneous and that the fluctuation in the electrostatic peptide-water interaction energy is about 10 times larger than the thermal energy at room temperature. This indicates that the peptide-solvent interactions are indeed important for the activation of chemical reactions in aqueous solution.

Redistribution of carbonyl stretch mode energy in isolated and solvated N-methylacetamide: kinetic energy spectral density analyses.

The vibrational energy transfer from the excited carbonyl stretch mode in N-deuterated N-methylacetamide (NMA-d), both in isolation and in a heavy water cluster, is studied with nonequilibrium molecular dynamics (NEMD) simulations, employing a quantum mechanical/molecular mechanical (QM∕MM) force field at the semiempirical PM3 level. The nonequilibrium ensemble of vibrationally excited NMA-d is prepared by perturbing the positions and velocities of the carbonyl C and O atoms and its NEMD trajectories are obtained with a leap-frog algorithm properly modified for the initial perturbation. In addition to the time-domain analysis of the kinetic and potential energies, a novel method for the spectral analysis of the atomic kinetic energies is developed, in terms of the spectral density of kinetic energy, which provides the time-dependent changes of the frequency-resolved kinetic energies without the complications of normal mode analysis at every MD time step. Due to the QM description of the solute electronic structure, the couplings among the normal modes are captured more realistically than with classical force fields. The energy transfer in the isolated NMA-d is found to proceed first from the carbonyl bond to other modes with time scales of 3 ps or less, and then among the other modes over 3-21 ps. In the solvated NMA-d, most of the excess energy is first transferred to other intramolecular modes within 5 ps, which is subsequently dissipated to solvent with 7-19 ps time scales. The contribution of the direct energy transfer from the carbonyl bond to solvent was only 5% with ~7 ps time scale. Solvent reorganization that leads to destabilization of the electrostatic interactions is found to be crucial in the long time relaxation of the excess energy, while the water intramolecular modes do not contribute significantly. Detailed mode-specific energy transfer pathways are deduced for the isolated and solvated NMA-d and they show that the energy transfer in NMA-d is a highly cooperative process among the intramolecular modes and there is no single dominant pathway with more than 30% of transient contribution.

Computational vibrational spectroscopy of peptides and proteins in one and two dimensions.

Vibrational spectroscopy provides direct information on molecular environment and motions but, its interpretation is often hampered by band broadening. Over the past decade, two-dimensional (2D) vibrational spectroscopy has emerged as a promising technique to overcome a number of difficulties associated with linear spectroscopy and provided significantly detailed information on the structure and dynamics of complex molecules in condensed phases. This Account reviews recently developed computational methods used to simulate 1D and 2D vibrational spectra. The central quantity to calculate in computational spectroscopy is the spectroscopic response function, which is the product of many contributing factors such as vibrational transition energies, transition moments, and their modulations by fluctuating local environment around a solute. Accurate calculations of such linear and nonlinear responses thus require a concerted effort employing a wide range of methods including electronic structure calculation (ESC) and molecular dynamics (MD) simulation. The electronic structure calculation can provide fundamental quantities such as normal-mode frequencies and transition multipole strengths. However, since the treatable system size is limited with this method, classical MD simulation has also been used to account for the dynamics of the solvent environment. To achieve chemical accuracy, these two results are combined to generate time series of fluctuating transition frequencies and transition moments with the distributed multipole analysis, and this particular approach has been known as the hybrid ESC/MD method. For coupled multichromophore systems, vibrational properties of each chromophore such as a peptide are individually calculated by electronic structure methods and the Hessian matrix reconstruction scheme was used to obtain local mode frequencies and couplings of constituting anharmonic oscillators. The spectra thus obtained, especially for biomolecules including polypeptides and proteins, have proven to be reliable and in good agreement with experimental spectra. An alternative to the hybrid method has also been developed, where the classical limit of the vibrational response function was considered. Its main attraction is the capability to obtain the spectra directly from a set of MD trajectories. A novel development along this direction has been achieved by using quantum mechanical/molecular mechanical (QM/MM) force fields for the accurate description of vibrational anharmonicity and chromophore polarization effects. The latter aspects are critical in the 2D case because classical force fields employing harmonic intramolecular potential cannot produce reliable 2D signal. We anticipate that the computational methods presented here will continue to evolve along with experimental advancements and will be of use to further elucidate ultrafast dynamics of chemical and biological systems.

Molecular dynamics simulation study of a pulmonary surfactant film interacting with a carbonaceous nanoparticle.

This article reports an all-atom molecular dynamics simulation to study a model pulmonary surfactant film interacting with a carbonaceous nanoparticle. The pulmonary surfactant is modeled as a dipalmitoylphosphatidylcholine monolayer with a peptide consisting of the first 25 residues from surfactant protein B. The nanoparticle model with a chemical formula C188H53 was generated using a computational code for combustion conditions. The nanoparticle has a carbon cage structure reminiscent of the buckyballs with open ends. A series of molecular-scale structural and dynamical properties of the surfactant film in the absence and presence of nanoparticle are analyzed, including radial distribution functions, mean-square displacements of lipids and nanoparticle, chain tilt angle, and the surfactant protein B peptide helix tilt angle. The results show that the nanoparticle affects the structure and packing of the lipids and peptide in the film, and it appears that the nanoparticle and peptide repel each other. The ability of the nanoparticle to translocate the surfactant film is one of the most important predictions of this study. The potential of mean force for dragging the particle through the film provides such information. The reported potential of mean force suggests that the nanoparticle can easily penetrate the monolayer but further translocation to the water phase is energetically prohibitive. The implication is that nanoparticles can interact with the lung surfactant, as supported by recent experimental data by Bakshi et al.

Gating of the mechanosensitive channel protein MscL: the interplay of membrane and protein.

The mechanosensitive channel of large conductance (MscL) belongs to a family of transmembrane channel proteins in bacteria and functions as a safety valve that relieves the turgor pressure produced by osmotic downshock. MscL gating can be triggered solely by stretching of the membrane. This work reports an effort to understand this mechanotransduction by means of molecular dynamics (MD) simulation on the MscL of mycobacterium tuberculosis embedded in a palmitoyloleoylphosphatidylethanolamine membrane. Equilibrium MD under zero membrane tension produced a more compact protein structure, as measured by its radii of gyration, compared to the crystal structure, in agreement with previous experimental findings. Even under a large applied tension up to 1000 dyn/cm, the MscL lateral dimension largely remained unchanged after up to 20 ns of simulation. A nonequilibrium MD simulation of 3% membrane expansion showed a significant increase in membrane rigidity upon MscL inclusion, which can contribute to efficient mechanotransduction. Direct observation of channel opening was possible only when an explicit lateral bias force was applied to each of the five subunits of MscL in the radially outward direction. Using this force, open structures with a large pore of radius 10 A could be obtained. The channel opening takes place in a stepwise manner and concurrently with the water chain formation across the channel, which occurs without direct involvement of protein hydrophilic residues. The N-terminal S1 helices stabilize the open structure, and the membrane asymmetry (different lipid density on the two leaflets of membrane) promotes channel opening.

Structure of flexible and semiflexible polyelectrolyte chains in confined spaces of slit micro/nanochannels.

Understanding the behavior of a polyelectrolyte in confined spaces has direct relevance in design and manipulation of microfluidic devices, as well as transport in living organisms. In this paper, a coarse-grained model of anionic semiflexible polyelectrolyte is applied, and its structure and dynamics are fully examined with Brownian dynamics (BD) simulations both in bulk solution and under confinement between two negatively charged parallel plates. The modeling is based on the nonlinear bead-spring discretization of a continuous chain with additional long-range electrostatic, Lennard-Jones, and hydrodynamic interactions between pairs of beads. The authors also consider the steric and electrostatic interactions between the bead and the confining wall. Relevant model parameters are determined from experimental rheology data on the anionic polysaccharide xanthan reported previously. For comparison, both flexible and semiflexible models are developed accompanying zero and finite intrinsic persistence lengths, respectively. The conformational changes of the polyelectrolyte chain induced by confinements and their dependence on the screening effect of the electrolyte solution are faithfully characterized with BD simulations. Depending on the intrinsic rigidity and the medium ionic strength, the polyelectrolyte can be classified as flexible, semiflexible, or rigid. Confined flexible and semiflexible chains exhibit a nonmonotonic variation in size, as measured by the radius of gyration and end-to-end distance, with changing slit width. For the semiflexible chain, this is coupled to the variations in long-range bond vector correlation. The rigid chain, realized at low ionic strength, does not have minima in size but exhibits a sigmoidal transition. The size of confined semiflexible and rigid polyelectrolytes can be well described by the wormlike chain model once the electrostatic effects are taken into account by the persistence length measured at long length scale.