DNA groove preference shift upon phosphorylation of a protamine-like cationic peptide.

Protamines, arginine-rich DNA-binding proteins, are responsible for chromatin compaction in sperm cells, but their DNA groove preference, major or minor, is not clearly identified. We herein study the DNA groove preference of a short protamine-like cationic peptide before and after phosphorylation, using all-atom molecular dynamics and umbrella sampling simulations. According to various thermodynamic and structural analyses, a peptide in its non-phosphorylated native state prefers the minor groove over the major groove, but phosphorylation of the peptide bound to the minor groove not only reduces its binding affinity but also brings a serious deformation of the minor groove, eliminating the minor-groove preference. As protamines are heavily phosphorylated before binding to DNA, we expect that the structurally disordered phosphorylated protamines would prefer major grooves to enter into DNA during spermatogenesis.

Corrosion-driven droplet wetting on iron nanolayers.

The classical Evans' drop describes a drop of aqueous salt solution, placed on a bulk metal surface where it displays a corrosion pit that grows over time producing further oxide deposits from the metal dissolution. We focus here on the corrosion-induced droplet spreading using iron nanolayers whose semi-transparency allowed us to monitor both iron corrosion propagation and electrolyte droplet behavior by simple optical means. We thus observed that pits grow under the droplet and merge into a corrosion front. This front reached the triple contact line and drove a non radial spreading, until it propagated outside the immobile droplet. Such chemically-active wetting is only observed in the presence of a conductive substrate that provides strong adhesion of the iron nanofilm to the substrate. By revisiting the classic Evan's drop experiment on thick iron film, a weaker corrosion-driven droplet spreading is also identified. These results require further investigations, but they clearly open up new perspectives on substrate wetting by corrosion-like electrochemical reactions at the nanometer scale.

DNA-protamine condensates under low salt conditions: molecular dynamics simulation with a simple coarse-grained model focusing on electrostatic interactions.

Protamine, a small, strongly positively-charged protein, plays a key role in achieving chromatin condensation inside sperm cells and is also involved in the formulation of nanoparticles for gene therapy and packaging of mRNA-based vaccines against viral infection and cancer. The detailed mechanisms of such condensations are still poorly understood especially under low salt conditions where electrostatic interaction predominates. Our previous study, with a refined coarse-grained model in full consideration of the long-range electrostatic interactions, has demonstrated the crucial role of electrostatic interaction in protamine-controlled reversible DNA condensation. Therefore, we herein pay our attention only to the electrostatic interaction and devise a coarser-grained bead-spring model representing the right linear charge density on protamine and DNA chains but treating other short-range interactions as simply as possible, which would be suitable for real-scale simulations. Effective pair potential calculations and large-scale molecular dynamics simulations using this extremely simple model reproduce the phase behaviour of DNA in a wide range of protamine concentrations under low salt conditions, again revealing the importance of the electrostatic interaction in this process and providing a detailed nanoscale picture of bundle formation mediated by a charge disproportionation mechanism. Our simulations also show that protamine length alters DNA overcharging and in turn redissolution thresholds of DNA condensates, revealing the important role played by entropies and correlated fluctuations of condensing agents and thus offering an additional opportunity to design tailored nanoparticles for gene therapy. The control mechanism of DNA-protamine condensates will also provide a better microscopic picture of biomolecular condensates, i.e., membraneless organelles arising from liquid-liquid phase separation, that are emerging as key principles of intracellular organization. Such condensates controlled by post-translational modification of protamine, in particular phosphorylation, or by variations in protamine length from species to species may also be responsible for the chromatin-nucleoplasm patterning observed during spermatogenesis in several vertebrate and invertebrate species.

Protic Ionic Liquids for Intrinsically Stretchable Conductive Polymers.

Inspired by the classic hard-soft acid-base theory and intrigued by a theoretical prediction of spontaneous ion exchange between poly(3,4-ethylenedioxythiophene):poly(styrenesulfonate) (PEDOT:PSS) and hard-cation-soft-anion ionic liquid (IL), we treat PEDOT:PSS with a new IL composed of a protic (i.e., extremely hard) cation (3-methylimidazolium, p-MIM+) and an extremely soft anion (tetracyanoborate, TCB-). In fact, this protic IL (p-MIM:TCB) accomplishes the same levels of ion-exchange-mediated PEDOT-PSS separation, PEDOT-rich nanofibril formation, and electrical conductivity enhancement (∼2500 S/cm) as its aprotic counterpart (EMIM:TCB with 1-ethyl-3-methylimidazolium), the best IL used for this purpose so far. Furthermore, p-MIM:TCB significantly outperforms EMIM:TCB in terms of improving the stretchability (i.e., the highest tensile strain) of the PEDOT:PSS thin film. This enhancement is a result of the aromatic and protic cation p-MIM+, which acts as a molecular adhesive holding the exchanged ion pairs (PEDOT+:TCB----p-MIM+:PSS-) via ionic intercalation (at the surface of TCB--decorated PEDOT+ clusters) and hydrogen bonding (to PSS-), in which washing p-MIM+ out of the film degrades the stretchability while keeping the morphology. Our results offer molecular-level insight into the morphological, electrical, and mechanical properties of PEDOT:PSS and a molecular-interaction-based enhancement strategy that can be used for intrinsically stretchable conductive polymers.

Effect of phosphorylation of protamine-like cationic peptide on the binding affinity to DNA.

Protamines are more arginine-rich and more basic than histones and are responsible for providing a highly compacted shape to the sperm heads in the testis. Phosphorylation and dephosphorylation are two events that occur in the late phase of spermatogenesis before the maturation of sperms. In this work, we have studied the effect of phosphorylation of protamine-like cationic peptides using all-atom molecular dynamics simulations. Through thermodynamic analyses, we found that phosphorylation reduces the binding efficiency of such cationic peptides on DNA duplexes. Peptide phosphorylation leads to a less efficient DNA condensation, due to a competition between DNA-peptide and peptide-peptide interactions. We hypothesize that the decrease of peptide bonds between DNA together with peptide self-assembly might allow an optimal re-organization of chromatin and an efficient condensation through subsequent peptide dephosphorylation. Based on the globular and compact conformations of phosphorylated peptides mediated by arginine-phosphoserine H-bonding, we furthermore postulate that phosphorylated protamines could more easily intrude into chromatin and participate to histone release through disruption of histone-histone and histone-DNA binding during spermatogenesis.

Hard-Cation-Soft-Anion Ionic Liquids for PEDOT:PSS Treatment.

A promising conducting polymer poly(3,4-ethylenedioxythiophene):poly(styrenesulfonate) (PEDOT:PSS) experiences significant conductivity enhancement when treated with proper ionic liquids (ILs). Based on the hard-soft-acid-base principle, we propose a combination of a hydrophilic hard cation A+ (instead of the commonly used 1-ethyl-3-methyl imidazolium, EMIM+) and a hydrophobic soft anion X- (such as tetracyanoborate, TCB-) as the best ILs for this purpose. Such ILs would decouple hydrophilic-but-insulating PSS- from conducting-but-hydrophobic PEDOT+ most efficiently by strong interactions with hydrophilic A+ and hydrophobic X-, respectively. Such a favorable ion exchange between PEDOT+:PSS- and A+:X- ILs would allow the growth of conducting PEDOT+ domains decorated by X-, not disturbed by PSS- or A+. Using density functional theory calculations and molecular dynamics simulations, we demonstrate that a protic cation- (aliphatic N-alkyl pyrrolidinium, in particular) combined with the hydrophobic anion TCB- indeed outperforms EMIM+ by promptly leaving hydrophobic TCB- and strongly binding to hydrophilic PSS-.

DNA-assisted assembly of cationic gold nanoparticles: Monte Carlo simulation.

DNA-assisted assembly of ligand-stabilized gold nanoparticles is studied using Monte Carlo simulations with coarse-grained models for DNA and AuNP. Their interaction in a periodic simulation box is described by a combination of electrostatic and pairwise hard core potentials. We first probe the self-assembly of AuNPs resulting in an ordered distribution on a single fixed DNA strand. Subsequently, the effective force calculated between a pair of parallel DNA in the presence of AuNPs shows the attraction between them at short distance associated to a stable equilibrium position. Finally, the osmotic pressure calculated in a compact DNA-AuNP lattice with various amounts of monovalent salt ions shows that an increasing amount of salt prevents aggregate formation.

Protamine-Controlled Reversible DNA Packaging: A Molecular Glue.

Packaging paternal genome into tiny sperm nuclei during spermatogenesis requires 106-fold compaction of DNA, corresponding to a 10-20 times higher compaction than in somatic cells. While such a high level of compaction involves protamine, a small arginine-rich basic protein, the precise mechanism at play is still unclear. Effective pair potential calculations and large-scale molecular dynamics simulations using a simple idealized model incorporating solely electrostatic and steric interactions clearly demonstrate a reversible control on DNA condensates formation by varying the protamine-to-DNA ratio. Microscopic states and condensate structures occurring in semidilute solutions of short DNA fragments are in good agreement with experimental phase diagram and cryoTEM observations. The reversible microscopic mechanisms induced by protamination modulation should provide valuable information to improve a mechanistic understanding of early and intermediate stages of spermatogenesis where an interplay between condensation and liquid-liquid phase separation triggered by protamine expression and post-translational regulation might occur. Moreover, recent vaccines to prevent virus infections and cancers using protamine as a packaging and depackaging agent might be fine-tuned for improved efficiency using a protamination control.

Molecular Dynamics of PEDOT:PSS Treated with Ionic Liquids. Origin of Anion Dependence Leading to Cation Design Principles.

Conductivity enhancement of PEDOT:PSS via the morphological change of PEDOT-rich domains has been achieved by introducing a 1-ethyl-3-methylimidazolium (EMIM)-based ionic liquid (IL) into its aqueous solution, and the degree of such change varies drastically with the anion coupled to the EMIM cation constituting the IL. We carry out a series of molecular dynamics simulations on various simple model systems for the extremely complex mixtures of PEDOT:PSS and EMIM:X IL in water, varying the anion X, the IL concentration, the oligomer model of PEDOT:PSS, and the size of the model systems. The common characteristic found in all simulations is that although planar hydrophobic anions X are the most efficient for ion exchange between PEDOT:PSS and EMIM:X, they tend to bring together planar EMIM cations to PEDOT-rich domains, disrupting PEDOT π-stacks with PEDOT-X-EMIM intercalating layers. Nonplanar hydrophobic anions, which leave most of EMIM cations in water, are efficient for both ion exchange and the formation of extended PEDOT π-stacks, as observed in experiments. Based on such findings, we propose a design principle for new cations replacing EMIM; nonplanar hydrophilic cations combined with hydrophobic anions should improve IL efficiency for PEDOT:PSS treatment.

Protamine Binding Site on DNA: Molecular Dynamics Simulations and Free Energy Calculations with Full Atomistic Details.

Protamine, an arginine-rich basic protein, compacts DNAs in sperm nuclei to densities higher than those in somatic cells. The mechanism of this compaction in sperm cells is even less clear than in somatic cells. Even the preferred binding site, if any, of protamine on DNA is not clearly identified. In this work, we carry out fully atomistic (or all-atom) molecular dynamics simulations to estimate the relative stabilities of protamine binding sites on DNA. Free energy calculated with umbrella sampling on a short arginine stretch bound to the major and minor grooves suggests that a short arginine stretch would prefer the DNA major groove as its binding site. Complementary umbrella sampling simulations where an arginine stretch or a whole protamine is transferred from the major to the minor groove also lead to the same conclusion. We find that the protamine located in the major groove better utilizes the DNA backbone as the binding site and represents the best compromise between enthalpy and entropy gain.

Ionic Liquid for PEDOT:PSS Treatment. Ion Binding Free Energy in Water Revealing the Importance of Anion Hydrophobicity.

Water solubility of PEDOT:PSS conducting polymer is achieved by PSS at the expense of disturbing the crystallinity and electron mobility of PEDOT. Recently, PEDOT crystallinity and electron mobility have been improved by treating the PEDOT:PSS aqueous solution with 1-ethyl-3-methylimidazolium- (EMIM-) based ionic liquids (IL) EMIM:X. The amount of such improvement varies drastically with the anion X coupled to EMIM cation in the IL. Herein, using umbrella-sampling molecular dynamics simulations on the aqueous solutions of a minimal model of PEDOT:PSS mixed with various EMIM:X ILs, we show that the solvation of each ion in water plays a major role in the free energies of ion binding and exchange. Anions X efficient for the improvement are weakly stabilized by hydration (i.e., hydrophobic) and prefer binding to hydrophobic PEDOT than to hydrophilic EMIM, favoring the ion exchange. In order to fulfill our design principle, a quantitative criterion based on hydration free energy is proposed to select efficient hydrophobic anions X.

Liquid crystal ordering of nucleic acids.

Several analytical calculations and computer simulations propose that cylindrical monodispersive rods having an aspect ratio (ratio of length to diameter) greater than 4 can exhibit liquid crystal (LC) ordering. But, recent experiments demonstrated the signature of LC ordering in systems of 4- to 20-base pair (bp) long nucleic acids (NAs) that do not satisfy the shape anisotropy criterion. Mechanisms of end-to-end adhesion and stacking have been proposed to explain this phenomenon. In this study, using all-atom molecular dynamics (MD) simulation, we explicitly verify the end-to-end stacking of double-stranded RNA (dsRNA) and demonstrate the LC ordering at the microscopic level. Using umbrella sampling (US) calculation, we quantify the potential of mean force (PMF) between two dsRNAs for various reaction coordinates (RCs) and compare our results with previously reported PMFs for double-stranded DNA (dsDNA). The PMF profiles demonstrate the anisotropic nature of inter-NA interaction. We find that, like dsDNA, dsRNA also prefers to stack on top of each other while repelling sideways, leading to the formation of supra-molecular-columns that undergo LC ordering at high NA volume fraction (φ). We also demonstrate and quantify the nematic ordering of the RNAs using several hundred nanosecond-long MD simulations that remain almost invariant for different initial configurations and under different external physiological conditions.

Orientation Dependence of Inter-NCP Interaction: Insights into the Behavior of Liquid Crystal Phase and Chromatin Fiber Organization.

We report equilibrium and nonequilibrium molecular dynamics (MD) simulations of two nucleosome core particles (NCPs) stacked with their dyad axes oriented in parallel or antiparallel fashion. From the equilibrium trajectories, we determine the bridging behavior of different histone tails and observe that different sets of histone tails play important roles in the two orientations in stabilizing the NCP stack. While the H4 and H2A tails play important intermediary roles in the parallel stack, the H3 and H2B tails are important in the antiparallel stack. We use steered MD simulations to unstack the two NCPs and find a stark difference in their unstacking pathways. While the average rupture force was found to be higher for the parallel stack, the work done for complete unstacking was similar for both orientations. We use Jarzynski equality to determine the PMF profiles along the unstacking pathway, relate our findings to the behavior of NCP mesophases, and derive insights into the enigmatic nucleosomal organization in the chromatin fiber.

Interactions between the Molecular Components of the Cowpea Chlorotic Mottle Virus Investigated by Molecular Dynamics Simulations.

The formation of a viral particle generally involves hundreds of proteins, making the assembly process intricate. Despite its intrinsic complexity, the production of a viral particle begins through the interaction between the basic assembly components. For the cowpea chlorotic mottle virus (CCMV), the first steps of the assembly process involve dimers of the capsid protein. Here, we carried out atomistic molecular dynamics simulations to investigate the initial assembly process of CCMV to get insight into the interactions at the molecular level. We found that salinity not only affects the electrostatic interactions between dimers but also changes the conformation of the positively charged N-terminal tails and can cause a serious steric hindrance for other dimers binding to the hydrophobic domains. An RNA rod was used to mimic a long segment of a viral genome and to study its interaction with dimers. We observed that the dimer with tails prefers to bind on the RNA rod with its positively charged inner side. The dimer-RNA interaction was found to be as strong as the dimer-dimer interaction, whereas the association energies between a dimer and a pentamer or a hexamer of dimers were high but strongly depended on the presence of the tails. Upon heating, the capsid experienced a shrinkage accompanied by a loss of order in the icosahedral crystal structure.

Sub-nanometer pore formation in single-molecule-thick polyurea molecular-sieving membrane: a computational study.

A polymeric network of 1-(4-tritylphenyl)urea (TPU) built via layer-by-layer cross-linking polymerization has been proposed to be an excellent mesh equipped with single-molecule-thick pores (i.e., cyclic poly-TPU rings), which can sieve glucose (∼0.7 nm) out of its mixture with urea for hemodialysis applications. Monte Carlo search for the lowest-energy conformation of various sizes of poly-TPU rings unravels the origin of narrow pore size distribution, which is around the sizes of dimer and trimer rings (0.3-0.8 nm). Flexible rings larger than the dimer and trimer rings, in particular tetramer rings, prefer a twisted conformation in the shape of the infinity symbol (∞, which looks like two dimer rings joined together) locked by a hydrogen bond between diphenylurea linker groups facing each other. Translocation energy profiles across these TPU rings reveal their urea-versus-glucose sieving mechanism: glucose is either too large (to enter dimers and twisted tetramers) or too perfectly fit (to exit trimers), leaving only a dimer-sized free space in the ring, whereas smaller-sized urea and water pass through these effective dimer-sized rings (bare dimers, twisted tetramers, and glucose-filled trimers) without encountering a substantial energy barrier or trap.

Ionic Liquid Designed for PEDOT:PSS Conductivity Enhancement.

Poly-3,4-ethylenedioxythiophene:polystyrenesulfonate (PEDOT:PSS) is a water-processable conducting polymer with promise for use in transparent flexible electrodes and thermoelectric devices, but its conductivity is not satisfactory. Its low conductivity is attributed to the formation of hydrophilic/insulating PSS outer layers encapsulating the conducting/hydrophobic p-doped PEDOT cores. Recently a significant conductivity enhancement has been achieved by adding ionic liquid (IL). It is believed that ion exchange between PEDOT:PSS and IL components helps PEDOT to decouple from PSS and to grow into large-scale conducting domains, but the exact mechanism is still under debate. Here we show through free energy calculations using density functional theory on a minimal model that the most efficient IL pairs are the least tightly bound ones with the lowest binding energies, which would lead to the most efficient ion exchange with PEDOT:PSS. This spontaneous ion exchange followed by nanophase segregation between PEDOT and PSS, with formation of a π-stacked PEDOT aggregate decorated by IL anions, is also supported by molecular dynamics performed on larger PEDOT:PSS models in solution. We also show that the most efficient IL anions would sustain the highest amount of charge carriers uniformly distributed along the PEDOT backbone to further enhance the conductivity, providing that they remain in the PEDOT domain after the ion exchange. Hence, our design principle is that the high-performance IL should induce not only an efficient ion exchange with PEDOT:PSS to improve the PEDOT morphology (to increase mobility) but also a uniform high-level p-doping of PEDOT (to enhance intrinsic conductivity). Based on this principle, a promising (electron-withdrawing, but bulky, soft, and hydrophobic) new IL pair is proposed.

Investigating the thermal dissociation of viral capsid by lattice model.

The dissociation of icosahedral viral capsids was investigated by a homogeneous and a heterogeneous lattice model. In thermal dissociation experiments with cowpea chlorotic mottle virus and probed by small-angle neutron scattering, we observed a slight shrinkage of viral capsids, which can be related to the strengthening of the hydrophobic interaction between subunits at increasing temperature. By considering the temperature dependence of hydrophobic interaction in the homogeneous lattice model, we were able to give a better estimate of the effective charge. In the heterogeneous lattice model, two sets of lattice sites represented different capsid subunits with asymmetric interaction strengths. In that case, the dissociation of capsids was found to shift from a sharp one-step transition to a gradual two-step transition by weakening the hydrophobic interaction between AB and CC subunits. We anticipate that such lattice models will shed further light on the statistical mechanics underlying virus assembly and disassembly.

Shuttlecock-Shaped Molecular Rectifier: Asymmetric Electron Transport Coupled with Controlled Molecular Motion.

A fullerene derivative with five hydroxyphenyl groups attached around a pentagon, (4-HOC6H4)5HC60 (1), has shown an asymmetric current-voltage (I-V) curve in a conducting atomic force microscopy experiment on gold. Such molecular rectification has been ascribed to the asymmetric distribution of frontier molecular orbitals over its shuttlecock-shaped structure. Our nonequilibrium Green's function (NEGF) calculations based on density functional theory (DFT) indeed exhibit an asymmetric I-V curve for 1 standing up between two Au(111) electrodes, but the resulting rectification ratio (RR ∼ 3) is insufficient to explain the wide range of RR observed in experiments performed under a high bias voltage. Therefore, we formulate a hypothesis that high RR (>10) may come from molecular orientation switching induced by a strong electric field applied between two electrodes. Indeed, molecular dynamics simulations of a self-assembled monolayer of 1 on Au(111) show that the orientation of 1 can be switched between standing-up and lying-on-the-side configurations in a manner to align its molecular dipole moment with the direction of the applied electric field. The DFT-NEGF calculations taking into account such field-induced reorientation between up and side configurations indeed yield RR of ∼13, which agrees well with the experimental value obtained under a high bias voltage.

Understanding the origin of liquid crystal ordering of ultrashort double-stranded DNA.

Recent experiments have shown that short double-stranded DNA (dsDNA) fragments having six- to 20-base pairs exhibit various liquid crystalline phases. This violates the condition of minimum molecular shape anisotropy that analytical theories demand for liquid crystalline ordering. It has been hypothesized that the liquid crystalline ordering is the result of end-to-end stacking of dsDNA to form long supramolecular columns which satisfy the shape anisotropy criterion necessary for ordering. To probe the thermodynamic feasibility of this process, we perform molecular dynamics simulations on ultrashort (four base pair long) dsDNA fragments, quantify the strong end-to-end attraction between them, and demonstrate that the nematic ordering of the self-assembled stacked columns is retained for a large range of temperature and salt concentration.

Lithium ion solvation by ethylene carbonates in lithium-ion battery electrolytes, revisited by density functional theory with the hybrid solvation model and free energy correction in solution.

Complex formation between lithium (Li(+)) ions and electrolyte molecules would affect the ionic conductivity through the electrolyte in lithium-ion batteries (LIBs). We hence revisit the solvation number of Li(+) in the most commonly used ethylene carbonate (EC) electrolyte. The solvation number n of Li(+)(EC)n in the first solvation shell of Li(+) is estimated on the basis of the free energy calculated by the density functional theory combined with a hybrid solvation model where the explicit solvation shell of Li(+) is immersed in a free volume of an implicit bulk solvent. This new hybrid solvation (implicit and explicit) model predicts the most probable solvation number (n = 4) and solvation free energy (-91.3 kcal mol(-1)) of Li(+) in a good agreement with those predicted by calculations employing simpler solvation models (either implicit or explicit). The desolvation (n = 2) of Li(0)(EC)n upon reduction near anodes is also well described with this new hybrid model.