Evolving better solvate electrolytes for lithium secondary batteries.

The overall performance of lithium batteries remains unmatched to this date. Decades of optimisation have resulted in long-lasting batteries with high energy density suitable for mobile applications. However, the electrolytes used at present suffer from low lithium transference numbers, which induces concentration polarisation and reduces efficiency of charging and discharging. Here we show how targeted modifications can be used to systematically evolve anion structural motifs which can yield electrolytes with high transference numbers. Using a multidisciplinary combination of theoretical and experimental approaches, we screened a large number of anions. Thus, we identified anions which reach lithium transference numbers around 0.9, surpassing conventional electrolytes. Specifically, we find that nitrile groups have a coordination tendency similar to SO2 and are capable of inducing the formation of Li+ rich clusters. In the bigger picture, we identified a balanced anion/solvent coordination tendency as one of the key design parameters.

Understanding the effects of ethanol on the liposome bilayer structure using microfluidic-based time-resolved small-angle X-ray scattering and molecular dynamics simulations.

Lipid nanoparticles (LNPs) are essential carrier particles in drug delivery systems, particularly in ribonucleic acid delivery. In preparing lipid-based nanoparticles, microfluidic-based ethanol injection may produce precisely size-controlled nanoparticles. Ethanol is critical in LNP formation and post-treatment processes and affects liposome size, structure, lamellarity, and drug-loading efficiency. However, the effects of time-dependent changes in the ethanol concentration on the structural dynamics of liposomes are not clearly understood. Herein, we investigated ethanol-induced lipid bilayer changes in liposomes on a time scale from microseconds to tens of seconds using a microfluidic-based small-angle X-ray scattering (SAXS) measurement system coupled with molecular dynamics (MD) simulations. The time-resolved SAXS measurement system revealed that single unilamellar liposomes were converted to multilamellar liposomes within 0.8 s of contact with ethanol, and the d-spacing was decreased from 6.1 (w/o ethanol) to 4.4 nm (80% ethanol) with increasing ethanol concentration. We conducted 1 µs MD simulations to understand the molecular-level structural changes in the liposomes. The MD simulations revealed that the changes in the lamellar structure caused by ethanol at the molecular level could explain the structural changes in the liposomes observed via time-resolved SAXS. Therefore, the post-treatment process to remove residual ethanol is critical in liposome formation.

Light-controllable cell-membrane disturbance for intracellular delivery.

Highly polar and charged molecules, such as oligonucleotides, face significant barriers in crossing the cell membrane to access the cytoplasm. To address this problem, we developed a light-triggered twistable tetraphenylethene (TPE) derivative, TPE-C-N, to facilitate the intracellular delivery of charged molecules through an endocytosis-independent pathway. The central double bond of TPE in TPE-C-N is planar in the ground state but becomes twisted in the excited state. Under light irradiation, this planar-to-twisted structural change induces continuous cell membrane disturbances. Such disturbance does not lead to permanent damage to the cell membrane. TPE-C-N significantly enhanced the intracellular delivery of negatively charged molecules under light irradiation when endocytosis was inhibited through low-temperature treatment, confirming the endocytosis-independent nature of this delivery method. We have successfully demonstrated that the TPE-C-N-mediated light-controllable method can efficiently promote the intracellular delivery of charged molecules, such as peptides and oligonucleotides, with molecular weights ranging from 1000 to 5000 Da.

Protein-lipid acyl chain interactions: Depth-dependent changes of segmental mobility of phospholipid in contact with bacteriorhodopsin.

Boundary lipids surrounding membrane proteins play an essential role in protein function and structure. These protein-lipid interactions are mainly divided into electrostatic interactions between the polar amino acids of proteins and polar heads of phospholipids, and hydrophobic interactions between protein transmembrane sites and phospholipid acyl chains. Our previous report (Kawatake et al., Biochim. Biophys. Acta 1858 [2016] 2106-2115) covered a method for selectively analyzing boundary lipid interactions and showed differences in membrane protein-peripheral lipid interactions due to differences in their head group. Interactions in the hydrophobic acyl chains of phospholipids are relatively consistent among proteins, but the details of these interactions have not been elucidated. In this study, we reconstituted bacteriorhodopsin as a model protein into phospholipid membranes labeled with 2H and 13C for solid-state NMR measurement to investigate the depth-dependent effect of the head group structure on the lipid bilayer. The results showed that the position of the phospholipid near the carbonyl carbon was affected by the head group in terms of selectivity for protein surfaces, whereas in the deep interior of the bilayer near the leaflet interface, there was little difference between the head groups, indicating that the dependence of their interactions on the head group was much reduced.

pSPICA Force Field Extended for Proteins and Peptides.

Many coarse-grained (CG) molecular dynamics (MD) studies have been performed to investigate biological processes involving proteins and lipids. CG force fields (FFs) in these MD studies often use implicit or nonpolar water models to reduce computational costs. CG-MD using water models cannot properly describe electrostatic screening effects owing to the hydration of ionic segments and thus cannot appropriately describe molecular events involving water channels and pores through lipid membranes. To overcome this issue, we developed a protein model in the pSPICA FF, in which a polar CG water model showing the proper dielectric response was adopted. The developed CG model greatly improved the transfer free energy profiles of charged side chain analogues across the lipid membrane. Application studies on melittin-induced membrane pores and mechanosensitive channels in lipid membranes demonstrated that CG-MDs using the pSPICA FF correctly reproduced the structure and stability of the pores and channels. Furthermore, the adsorption behavior of the highly charged nona-arginine peptides on lipid membranes changed with salt concentration, indicating the pSPICA FF is also useful for simulating protein adsorption on membrane surfaces.

Molecular Level Origin of Ion Dynamics in Highly Concentrated Electrolytes.

Single-ion conducting liquid electrolytes are key to achieving rapid charge/discharge in Li secondary batteries. The Li+ transference (or transport) numbers are the defining properties of such electrolytes and have been discussed in the framework of concentrated solution theories. However, the connection between macroscopic transference and microscopic ion dynamics remains unclear. Molecular dynamics simulations were performed to obtain direct information regarding the microscopic behaviors in highly concentrated electrolytes, and the relationships between these behaviors and the transference number were determined under anion-blocking conditions. Various solvents with different donor numbers (DNs) were used along with a Li salt of the weakly Lewis basic bis(fluorosulfonyl)amide anion for electrolyte preparation. Favorable ordered Li+ structuring and a continuous Li+ conduction pathway were observed for the fluoroethylene carbonate-based electrolyte due to its low DN. The properties were less pronounced at higher DNs, e.g., for the dimethyl sulfoxide-based electrolyte. The τLi-solventlife/τdipolerelax ratio was introduced as a factor for ion dynamics, and the two mechanisms of ion transport were considered an exchange mechanism (τLi-solventlife/τdipolerelax < 1) and a vehicle mechanism (translational motion of solvated Li+) (τLi-solventlife/τdipolerelax ≥ 1). Vehicle-type transport was dominant with high DNs, while exchangeable transport was preferable at lower DNs. These findings should aid the further selection of solvents and Li salts to prepare single-ion conducting electrolytes.

Improved Protein Model in SPICA Force Field.

The previous version of the SPICA coarse-grained (CG) force field (FF) protein model focused primarily on membrane proteins and successfully reproduced the dimerization free energies of several transmembrane helices and the stable structures of various membrane protein assemblies. However, that model had limited accuracy when applied to other proteins, such as intrinsically disordered proteins (IDPs) and peripheral proteins, because the dimensions of the IDPs in an aqueous solution were too compact, and protein binding on the lipid membrane surface was overstabilized. To improve the accuracy of the SPICA FF model for the simulation of such systems, in this study, we introduce protein secondary structure-dependent nonbonded interaction parameters to the backbone segments and reoptimize almost all nonbonded parameters for amino acids. The improved FF proposed here successfully reproduces the radii of gyration of various IDPs, the binding sensitivity of several peripheral membrane proteins, and the dimerization free energies of several transmembrane helices. The new model also shows improved agreement with experiments on the free energy of peptide association in water. In addition, an extensive library of nonbonded interactions between proteins and lipids, including various glycerophospholipids, sphingolipids, and cholesterol, allows the study of specific interactions between lipids and peripheral and transmembrane proteins. Hence, the new SPICA FF (version 2) proposed herein is applicable with high accuracy for simulating a wide range of protein systems.

Pragmatic Coarse-Graining of Proteins: Models and Applications.

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

HER-2-Targeted Boron Neutron Capture Therapy with Carborane-integrated Immunoliposomes Prepared via an Exchanging Reaction.

Boron neutron capture therapy (BNCT) is a promising modality for cancer treatment because of its minimal invasiveness. To maximize the therapeutic benefits of BNCT, the development of efficient platforms for the delivery of boron agents is indispensable. Here, carborane-integrated immunoliposomes were prepared via an exchanging reaction to achieve HER-2-targeted BNCT. The conjugation of an anti-HER-2 antibody to carborane-integrated liposomes successfully endowed these liposomes with targeting properties toward HER-2-overexpressing human ovarian cancer cells (SK-OV3); the resulting BNCT activity toward SK-OV3 cells obtained using the current immunoliposomal system was 14-fold that of the l-BPA/fructose complex, which is a clinically available boron agent. Moreover, the growth of spheroids treated with this system followed by thermal neutron irradiation was significantly suppressed compared with treatment with the l-BPA/fructose complex.

Complex Energy Landscapes of Self-Assembled Vesicles.

The field of supramolecular chemistry has witnessed tremendous progress in bringing the system away from equilibrium for traditionally inaccessible structures and functions. Vesicular assemblies with complex energy landscapes and pathways, which are reminiscent of diverse cellular vesicles like exosomes, remain exceedingly rare. Here, relying on the activation of oligo(ethylene glycol) (OEG) interdigitation and the encoded conformational freedom in monodisperse Janus dendrimers, we reveal a rich landscape and a pathway selection of distinct vesicles. The interdigitation can be selectively switched on and off using temperature ramps, and the critical temperatures can be further determined by molecular design. Our findings suggest that synthetic vesicles, with different energy states and unexpected transition pathways, emulate dynamic cellular vesicles in nature. We anticipate that vesicles with an activated OEG corona conformation will open new routes for nanomedicine and advanced materials.

Molecular Dynamics Simulations of High-Concentration Li[TFSA] Sulfone Solution: Effect of Easy Conformation Change of Sulfolane on Fast Diffusion of Li Ion.

The parameters of the polarizable force field used for molecular dynamics simulations of Li diffusion in high-concentration lithium bis(trifluoromethanesulfonyl)amide (Li[TFSA]) sulfone (sulfolane, dimethylsulfone, ethylmethylsulfone, and ethyl-i-propylsulfone) solutions were refined. The densities of the solutions obtained by molecular dynamics simulations reproduced well the experimental values. The calculated concentration, temperature, and solvent dependencies of self-diffusion coefficients of ions and solvents in the mixtures well reproduce the experimentally observed dependencies. Ab initio calculations show that the intermolecular interactions between Li ions and four sulfones are not largely different. Conformational analyses show that sulfolane can change the conformation more easily owing to lower barrier height for pseudorotation compared to the rotational barrier heights of diethylsulfone and ethylmethylsulfone. Molecular dynamics simulations indicate that the easy conformation change of solvent affects the rotational relaxation of the solvent and the diffusion of Li ion in the mixture. The easy conformation change of sulfolane is one of the causes of faster diffusion of Li ion in the mixture of Li[TFSA] and sulfolane compared to the mixtures of smaller dimethylsulfone and ethylmethylsulfone.

Mutation detection of urinary cell-free DNA via catch-and-release isolation on nanowires for liquid biopsy.

Cell-free DNA (cfDNA) and extracellular vesicles (EVs) are molecular biomarkers in liquid biopsies that can be applied for cancer detection, which are known to carry information on the necessary conditions for oncogenesis and cancer cell-specific activities after oncogenesis, respectively. Analyses for both cfDNA and EVs from the same body fluid can provide insights into screening and identifying the molecular subtypes of cancer; however, a major bottleneck is the lack of efficient and standardized techniques for the isolation of cfDNA and EVs from clinical specimens. Here, we achieved catch-and-release isolation by hydrogen bond-mediated binding of cfDNA in urine to zinc oxide (ZnO) nanowires, which also capture EVs by surface charge, and subsequently we identified genetic mutations in urinary cfDNA. The binding strength of hydrogen bonds between single-crystal ZnO nanowires and DNA was found to be equal to or larger than that of conventional hydrophobic interactions, suggesting the possibility of isolating trace amounts of cfDNA. Our results demonstrated that nanowire-based cancer screening assay can screen cancer and can identify the molecular subtypes of cancer in urine from brain tumor patients through EV analysis and cfDNA mutation analysis. We anticipate our method to be a starting point for more sophisticated diagnostic models of cancer screening and identification.

Liquid Structures and Ion Dynamics of Ionic Liquids Viewed from Intermolecular Interactions.

The elucidation of the factors determining liquid structures and transport properties of ionic liquids is important for the design and development of ionic liquid electrolytes. This personal account introduces the importance of computational methods for studying ionic liquids. Molecular dynamics simulations provide detailed information on liquid structures of ionic liquid such as the structures of solvated cation complexes in equimolar mixtures of glymes and Li[TFSA] and the effects of the charges of electrode on liquid structure near the electrode. Ab initio calculations reveal that the magnitude of the attraction between ions and conformational flexibility ions play important roles in determining transport properties of ionic liquids. First principle molecular dynamics simulations elucidate why solvated cation complex is stable in the equimolar mixtures, although the Li+ -[TFSA]- interaction is greater than Li+ -glyme interaction.

Self-Assembly of Glycerol-Amphiphilic Janus Dendrimers Amplifies and Indicates Principles for the Selection of Stereochemistry by Biological Membranes.

The principles for the selection of the stereochemistry of phospholipids of biological membranes remain unclear and continue to be debated. Therefore, any new experiments on this topic may help progress in this field. To address this question, three libraries of constitutional isomeric glycerol-amphiphilic Janus dendrimers (JDs) with nonsymmetric homochiral, racemic, and symmetric achiral branching points were synthesized by an orthogonal-modular-convergent methodology. These JDs amplify self-assembly, and therefore, monodisperse vesicles known as dendrimersomes (DSs) with predictable dimensions programmed by JD concentration were assembled by rapid injection of their ethanol solution into water. DSs of homochiral JD enantiomers, racemic, including mixtures of different enantiomers, and achiral exhibited similar DS size-concentration dependence. However, the number of bilayers of DSs assembled from homochiral, achiral, and racemic JDs determined by cryo-TEM were different. Statistical analysis of the number of bilayers and coarse-grained molecular dynamics simulations demonstrated that homochiral JDs formed predominantly unilamellar DSs. Symmetric achiral JDs assembled only unilamellar DSs while racemic JDs favored multilamellar DSs. Since cell membranes are unilamellar, these results indicate a new rationale for nonsymmetric homochiral vs racemic selection. Simultaneously, these experiments imply that the symmetric achiral lipids forming more stable membrane, probably had been the preferable assemblies of prebiotic cell membranes.

Fluorescence Turn-on of Tetraphenylethylene Derivative by Transfer from Cyclodextrin to Liposomes, HeLa Cells, and E. coli.

Herein, trimethyl-ß-cyclodextrin (TMe-ß-CDx) and γ-cyclodextrin (γ-CDx) could dissolve a tetraphenylethylene derivative (TPE-OH4 ) in water through high-speed vibration milling. The fluorescence intensity of the TMe-ß-CDx-TPE-OH4 complex was much higher than that of the γ-CDx-TPE-OH4 complex, as the rotation of the central C=C double bond of TPE-OH4 after photoactivation was inhibited in a smaller TMe-ß-CDx cavity in comparison with the γ-CDx cavity. In contrast, the fluorescence intensity of the γ-CDx-TPE-OH4 complex was very weak; nevertheless, it increased after the addition of liposomes due to the transfer of TPE-OH4 from the γ-CDx cavity to the lipid membrane as a "turn-on" phenomenon. Furthermore, to apply temperature sensor, it was demonstrated that the fluorescence intensity in the liposomes depended on the phase-transition temperature. By using the fluorescence turn-on phenomenon, TPE-OH4 could detect the presence of HeLa cells and E. coli by fluorescence.

Amphotericin B assembles into seven-molecule ion channels: An NMR and molecular dynamics study.

Amphotericin B, an antifungal drug with a long history of use, forms fungicidal ion-permeable channels across cell membranes. Using solid-state nuclear magnetic resonance spectroscopy and molecular dynamics simulations, we experimentally elucidated the three-dimensional structure of the molecular assemblies formed by this drug in membranes in the presence of the fungal sterol ergosterol. A stable assembly consisting of seven drug molecules was observed to form an ion conductive channel. The structure is somewhat similar to the upper half of the barrel-stave model proposed in the 1970s but substantially different in the number of molecules and in their arrangement. The present structure explains many previous findings, including structure-activity relationships of the drug, which will be useful for improving drug efficacy and reducing adverse effects.

Adsorption characteristics of peptides on ω-functionalized self-assembled monolayers: a molecular dynamics study.

Molecular dynamics simulations were employed to investigate the adsorption behavior of a variety of amino-acid side-chain analogs (SCAs) and a ß-hairpin (HP7) peptide on a series of liquid-like self-assembled monolayers (SAMs) with terminal functional groups of -OH, -OCH3, -CH3, and -CF3. The relationships between the adsorption free energy of the SCAs and the interfacial properties of water on the SAMs were examined to determine the acute predictors of protein adsorption on the SAM surfaces. The structural changes of HP7 on the SAM surfaces were also investigated to understand the relationship between the surface nature and protein denaturation. It was found that the adsorption free energy of the SCAs was linearly related to the surface hydrophobicity, which was computed as the free energy of cavity formation near the SAM-water interfaces. In addition, the hydrophobic -CH3 and -CF3 SAMs produced substantial conformational changes in HP7 because of the strong hydrophobic attractions to the nonpolar side chains. The hydrophilic surface terminated by -OH also promoted structural changes in HP7 resulting from the formation of hydrogen bonds between the hydrophilic tail and HP7. Consequently, the moderate amphiphilic surface terminated by -OCH3 avoided the denaturation of HP7 most efficiently, thus improving the biocompatibility of the surface. In conclusion, these results provide a deep understanding of protein adsorption for a wide range of polymeric surfaces, and they can potentially aid the design of appropriate biocompatible coatings for medical applications.

Cooperative antimicrobial action of melittin on lipid membranes: A coarse-grained molecular dynamics study.

We conducted a series of coarse-grained molecular dynamics (CG-MD) simulations to investigate the complicated actions of melittin, which is an antimicrobial peptide (AMP) derived from honey bee venom, on a lipid membrane. To accurately simulate the AMP action, we developed and used a protein CG model as an extension of the pSPICA force field (FF), which was designed to reproduce several thermodynamic quantities and structural properties. At a low peptide-to-lipid (P/L) ratio (1/102), no defect was detected. At P/L = 1/51, toroidal pore formation was observed due to collective insertion of multiple melittin peptides from the N-termini. The pore formation was initiated by a local increase in membrane curvature in the vicinity of the peptide aggregate. At a higher P/L ratio (1/26), two more modes were detected, seemingly not controlled by the P/L ratio but by a local arrangement of melittin peptides: 1. Pore formation accompanied by lipid extraction by melittin peptides:a detergent-like mechanism. 2. A rapidly formed large pore in a significantly curved membrane: bursting. Thus, we observed three pore formation modes (toroidal pore formation, lipid extraction, and bursting) depending on the peptide concentration and local arrangement. These observations were consistent with experimental observations and hypothesized melittin modes. Through this study, we found that the local arrangements and population of melittin peptides and the area expansion rate by membrane deformation were key to the initiation of and competition among the multiple pore formation mechanisms.

SPICA Force Field for Proteins and Peptides.

A coarse-grained (CG) model for peptides and proteins was developed as an extension of the Surface Property fItting Coarse grAined (SPICA) force field (FF). The model was designed to examine membrane proteins that are fully compatible with the lipid membranes of the SPICA FF. A preliminary version of this protein model was created using thermodynamic properties, including the surface tension and density in the SPICA (formerly called SDK) FF. In this study, we improved the CG protein model to facilitate molecular dynamics (MD) simulations with a reproduction of multiple properties from both experiments and all-atom (AA) simulations. An elastic network model was adopted to maintain the secondary structure within a single chain. The side-chain analogues reproduced the transfer free energy profiles across the lipid membrane and demonstrated reasonable association free energy (potential of mean force) in water compared to those from AA MD. A series of peptides/proteins adsorbed onto or penetrated into the membrane simulated by the CG MD correctly predicted the penetration depths and tilt angles of peripheral and transmembrane peptides/proteins as comparable to those in the orientations of proteins in membranes (OPM) database. In addition, the dimerization free energies of several transmembrane helices within a lipid bilayer were comparable to those from experimental estimation. Application studies on a series of membrane protein assemblies, scramblases, and poliovirus capsids demonstrated the good performance of the SPICA FF.

All-atom molecular dynamics study of hepatitis B virus containing pregenome RNA in solution.

Immature hepatitis B virus (HBV) captures nucleotides in its capsid for reverse transcription. The nucleotides and nucleotide analog drugs, which are triphosphorylated and negatively charged in the cell, approach the capsid via diffusion and are absorbed into it. In this study, we performed a long-time molecular dynamics calculation of the entire HBV capsid containing pregenome RNA to investigate the interactions between the capsid and negatively charged substances. Electric field analysis demonstrated that negatively charged substances can approach the HBV capsid by thermal motion, avoiding spikes. The substances then migrate all over the floor of the HBV capsid. Finally, they find pores through which they can pass through the HBV capsid shell. Free energy profiles were calculated along these pores for small ions to understand their permeability through the pores. Anions (Cl-) showed higher free energy barriers than cations (Na+ and K+) through all pores, and the permeation rate of Cl- was eight times slower than that of K+ or Na+. Furthermore, the ions were more stable in the capsid than in the bulk water. Thus, the HBV capsid exerts ion selectivity for uptake and provides an environment for ions, such as nucleotides and nucleotide analog drugs, to be stabilized within the capsid.