Mechanism of Cationic Lipid Induced DNA Condensation: Lipid-DNA Coordination and Divalent Cation Charge Fluctuations

The condensation of nucleic acids by lipids is a widespread phenomenon in biology with crucial implications for drug delivery. However, the mechanisms of DNA assembly in lipid bilayers remain insufficiently understood due to challenges in measuring and assessing each component's contribution in the lipid-DNA-cation system. This study uses all-atom molecular dynamics simulations to investigate DNA condensation in cationic lipid bilayers. Our exhaustive exploration of the thermodynamic factors reveals unique roles for phospholipid head groups and cations. We observed that bridging cations between lipid and DNA drastically reduce charges, while mobile magnesium cations "ping-ponging" between double strands create charge fluctuations. While the first factor stabilizes the DNA-lipid complex, the latter creates attractive forces to induce the spontaneous condensation of DNAs. This novel mechanism not only sheds light on the current data regarding cationic lipid-induced DNA condensation but also provides potential design strategies for creating efficient gene delivery vectors for drug delivery.

Xeno Nucleic Acids as Functional Materials: From Biophysical Properties to Application

Xeno nucleic acid (XNA) are artificial nucleic acids, in which the chemical composition of the sugar moiety is changed. These modifications impart distinct physical and chemical properties to XNAs, leading to changes in their biological, chemical, and physical stability. Additionally, these alterations influence the binding dynamics of XNAs to their target molecules. Consequently, XNAs find expanded applications as functional materials in diverse fields. This review provides a comprehensive summary of the distinctive biophysical properties exhibited by various modified XNAs and explores their applications as innovative functional materials in expanded fields.

Counterion Release from Macroion Assemblies of Planar Geometry

Macroions such as nanoparticles, polyelectrolytes, ionic gels, and amphiphiles can form condensed, often self-assembled, phases that are embedded in a solvent region. The condensed phase contains not only the partially or fully immobile charges of their macroions but also corresponding counterions that are mobile and thus free to migrate out of their confinement into the solvent region where they benefit from high translational entropy. Based on the nonlinear Poisson-Boltzmann model for monovalent ions, we quantify the corresponding fraction of released counterions for a planar slab geometry of the macroion phase. Slab thickness, extension of the solvent phase, fixed background charge density provided by the macroions, and dielectric constants inside slab and solvent combine into three dimensionless parameters that the fraction of released counterions depends on. We calculate that fraction and analyze the limits of a thin macroion phase, a large solvent phase, and linearized theory, where simple analytic results become available. Of particular interest is the presence of a single-planar interface that separates a bulk macroion phase from an extended solvent region. We calculate the apparent surface charge density that emerges due to the released counterions. Our model yields a comprehensive description of counterion partitioning between a planar macroion phase and a solvent region on the level of mean-field electrostatics in the absence of added salt ions.

DNA functionalized programmable hybrid biomaterials for targeted multiplexed applications

With the advent of DNA nanotechnology, DNA-based biomaterials have emerged as a unique class of materials at the center of various biological advances. Owing to DNA's high modification capacity via programmable Watson-Crick base-pairing, DNA structures of desired design with increased complexity have been developed. However, the limited scalability, along with poor mechanical properties, high synthesis costs, and poor stability, reduced the adaptability of DNA-based materials to complex biological applications. DNA-based hybrid biomaterials were designed to overcome these limitations by conjugating DNA with functional materials. Today, DNA-based hybrid materials have attracted significant attention in biological engineering with broad application prospects in biomedicine, clinical diagnosis, and nanodevices. Here, we summarize the recent advances in DNA-based hybrid materials with an in-depth understanding of general molecular design principles, functionalities, and applications. Finally, the challenges and prospects associated with DNA-based hybrid materials are discussed at the end of this review.

DNA-Targeted Complexes of Tc and Re for Biomedical Applications

Due to their favorable chemical features, Re and Tc complexes have been widely used for the development of new therapeutic agents and imaging probes to solve problems of biomedical relevance. This review provides an update of the most relevant research efforts towards the development of novel cancer theranostic agents using Re and Tc-based compounds interacting with specific DNA structures. This includes a variety of homometallic complexes, namely those containing M(CO)3 (M=Re, Tc) moieties, that exhibit different modes of interaction with DNA, such as covalent binding, intercalation, groove binding or G-quadruplex DNA binding. Additionally, heterometallic complexes, designed to potentiate synergistic effects of different metal centers to improve DNA-targeting, cytotoxicity and fluorescence properties, are also reviewed. Particular attention is also given to 99m Tc- and 188 Re-labeled oligonucleotides that have been widely explored to develop imaging and therapeutic radiopharmaceuticals through the in vivo hybridization with a specific complementary DNA or RNA target sequence to provide useful molecular tools in precision medicine for cancer diagnosis and treatment. Finally, the need for further improvement of DNA-targeted Re and Tc-based compounds as potential therapeutic and diagnostic agents is highlighted, and future directions are discussed.

Base flipping mechanism and binding strength of methyl-damaged DNA during the interaction with AGT

Methyl damage to DNA bases is common in the cell nucleus. O6-alkylguanine-DNA alkyl transferase (AGT) may be a promising candidate for direct damage reversal in methylated DNA (mDNA) at the O6 point of the guanine. Indeed, atomic-level investigations in the contact region of AGT-DNA complex can provide an in-depth understanding of their binding mechanism, allowing to evaluate the silico-drug nature of AGT and its utility in removing methyl damage in DNA. In this study, molecular dynamics (MD) simulation was utilized to examine the flipping of methylated nucleotide, the binding mechanism between mDNA and AGT, and the comparison of binding strength prior and post methyl transfer to AGT. The study reveals that methylation at the O6 atom of guanine weakens the hydrogen bond (H-bond) between guanine and cytosine, permitting for the flipping of such nucleotide. The formation of a H-bond between the base pair of methylated nucleotide (i.e., cytosine) and the intercalated arginine of AGT also forces the nucleotide to rotate. Following that, electrostatics and van der Waals contacts as well as hydrogen bonding contribute to form the complex of DNA and protein. The stronger binding of AGT with DNA before methyl transfer creates the suitable condition to transfer methyl adduct from DNA to AGT.

Controlling Chiroptical Responses via Chemo-Mechanical Deformation of DNA Origami Structures

DNA origami-based templates have been widely used to fabricate chiral plasmonic metamaterials due to their precise control of the placement of nanoparticles (NPs) in a desired configuration. However, achieving various chiroptical responses inevitably requires a change in the structure of DNA origami-based templates or binding sites on them, leading to the use of significantly different sets of DNA strands. Here, we propose an approach to controlling various chiroptical responses with a single DNA origami design using its chemo-mechanical deformation induced by DNA intercalators. The chiroptical response could be finely tuned by altering the concentration of intercalators only. The silver (Ag) enhancement was used to amplify the chiroptical signal by enlarging NPs and to maintain it by stiffening the template DNA structure. Furthermore, the sensitivity in the chiroptical signal change to the concentration of intercalators could be modulated by the type of intercalator, the mixture of two intercalators, and the stiffness of DNA origami structures. This approach would be useful in a variety of optical applications that require programmed spatial modification of chiroptical responses.

Electrostatic interactions between charge regulated spherical macroions

We study the interaction between two charge regulating spherical macroions with dielectric interior and dissociable surface groups immersed in a monovalent electrolyte solution. The charge dissociation is modelled via the Frumkin-Fowler-Guggenheim isotherm, which allows for multiple adsorption equilibrium states. The interactions are derived from the solutions of the mean-field Poisson-Boltzmann type theory with charge regulation boundary conditions. For a range of conditions we find symmetry breaking transitions from symmetric to asymmetric charge distribution exhibiting annealed charge patchiness, which results in like-charge attraction even in a univalent electrolyte-thus fundamentally modifying the nature of electrostatic interactions in charge-stabilized colloidal suspensions.

Reductant-dependent DNA-templated silver nanoparticle formation kinetics

DNA molecules have been demonstrated to be good templates for producing silver nanoparticles (AgNPs), with the advantages of well-controlled sizes, shapes, and properties. Revealing the formation kinetics of DNA-templated AgNPs is crucial for their efficient synthesis. Herein, using magnetic tweezers, we studied the reduction kinetics of the Ag+-DNA structure and the subsequent nucleation kinetics by adding NaBH4, L-ascorbic acid, and sodium citrate solutions. At [Ag+] = 0.01 mM, the addition of NaBH4 solution with the same concentration resulted in the restoration of DNA. In contrast, by increasing the [NaBH4]/[Ag+] ratio (r) to 10 and 100, the DNA extension initially decreased rapidly and then increased, indicating nucleation-dissolution kinetics. With AgNO3 solutions of higher concentrations (0.1 mM and 1 mM), direct particle nucleation and growth kinetics were observed by adding a tenfold (r = 10) or a hundredfold (r = 100) amount of NaBH4, which were evidenced by a significant reduction in DNA extension. The reductant dependence of the kinetics was further investigated. Addition of L-ascorbic acid to the DNA-Ag+ solution yielded an increase-decrease kinetics that was different from that caused by NaBH4, suggesting that nucleation was not initially favored due to the lack of sufficient Ag atoms; while sodium citrate showed a weak nucleation-promoting ability to form AgNPs. We discussed the findings within the framework of classical nucleation theory, in which the supersaturation of the Ag atom is strongly influenced by multiple factors (including the reducing ability of the reductant), resulting in different kinetics.

Conformal order and Poincaré-Klein mapping underlying electrostatics-driven inhomogeneity in tethered membranes

Understanding the organization of matter under the long-range electrostatic force is a fundamental problem in multiple fields. In this work, based on the electrically charged tethered membrane model, we reveal regular structures underlying the lowest-energy states of inhomogeneously stretched planar lattices by a combination of numerical simulation and analytical geometric analysis. Specifically we show the conformal order characterized by the preserved bond angle in the lattice deformation and reveal the Poincaré-Klein mapping underlying the electrostatics-driven inhomogeneity. The discovery of the Poincaré-Klein mapping, which connects the Poincaré disk and the Klein disk for the hyperbolic plane, implies the connection of long-range electrostatic force and hyperbolic geometry. We also discuss lattices with patterned charges of opposite signs for modulating in-plane inhomogeneity and even creating 3D shapes, which may have a connection to metamaterials design. This work suggests the geometric analysis as a promising approach for elucidating the organization of matter under the long-range force.

Machine Learning-Assisted Nanoenzyme/Bioenzyme Dual-Coupled Array for Rapid Detection of Amyloids

Array-based sensing methods offer significant advantages in the simultaneous detection of multiple amyloid biomarkers and thus have great potential for diagnosing early-stage Alzheimer's disease. Yet, detecting low concentrations of amyloids remains exceptionally challenging. Here, we have developed a fluorescent sensor array based on the dual coupling of a nanoenzyme (AuNPs) and bioenzyme (horseradish peroxidase) to detect amyloids. Various ss-DNAs were bound to the nanoenzyme for regulating enzymatic activity and recognizing amyloids. A simplified sensor array was generated from a screening model via machine learning algorithms and achieved signal amplification through a two-step enzymatic reaction. As a result, our sensing system could discriminate the aggregation species and aggregation kinetics at 200 nM with 100% accuracy. Moreover, AD model mice and healthy mice were distinguished with 100% accuracy through the sensor array, providing a powerful sensing platform for diagnosing AD.

Monovalent ion-mediated charge-charge interactions drive aggregation of surface-functionalized gold nanoparticles

Monolayer-protected metal nanoparticles (NPs) are not only promising materials with a wide range of potential industrial and biological applications, but they are also a powerful tool to investigate the behaviour of matter at nanoscopic scales, including the stability of dispersions and colloidal systems. This stability is dependent on a delicate balance between attractive and repulsive interactions that occur in the solution, and it is described in quantitative terms by the classic Derjaguin-Landau-Vewey-Overbeek (DLVO) theory, that posits that aggregation between NPs is driven by van der Waals interactions and opposed by electrostatic interactions. To investigate the limits of this theory at the nanoscale, where the continuum assumptions required by the DLVO theory break down, here we investigate NP dimerization by computing the Potential of Mean Force (PMF) of this process using fully atomistic MD simulations. Serendipitously, we find that electrostatic interactions can lead to the formation of metastable NP dimers at physiological ion concentrations. These dimers are stabilized by complexes formed by negatively charged ligands belonging to distinct NPs that are bridged by positively charged monovalent ions present in solution. We validate our findings by collecting tomographic EM images of NPs in solution and by quantifying their radial distribution function, that shows a marked peak at interparticle distance comparable with that of MD simulations. Taken together, our results suggest that not only van der Waals interactions, but also electrostatic interactions mediated by monovalent ions at physiological concentrations, contribute to attraction between nano-sized charged objects at very short length scales.

Deciphering the nature of binding of dexlansoprazole with DNA: Biophysical and docking approaches

Drugs, in general, exhibit their pharmacological activity in binding with intracellular targets. Numerous anticancer and antibacterial drugs target DNA as one of their primary intracellular targets. Dexlansoprazole (DLP) is a heterocyclic compound containing benzimidazole moiety and a proton pump inhibitor used to treat gastroesophageal reflux disease. The interaction of dexlansoprazole with calf thymus DNA (ct-DNA) has been studied using biophysical methods. The UV-Visible studies revealed a binding constant of 2.15 ± 0.3 × 104 M^-1 which is close to the value of 2.44 ± 0.3 × 10⁴ M^-1 obtained from the fluorescence studies. Competitive displacement studies using the fluorescence spectroscopic method with ethidium bromide and Hoechst as DNA markers suggested the groove binding mode of DLP in ct-DNA. The groove binding mode of DLP in ct-DNA was complemented by the results of viscosity and DNA melting studies. Further studies on the effect of ionic strength and potassium iodide on DLP binding with ct-DNA supported the observed binding mode. Circular dichroism studies reflected no significant conformational variation in ct-DNA after the interaction. The binding mode obtained from the experimental studies was corroborated by the molecular docking studies that showed the position of DLP in the minor groove of ct-DNA along with the receptor interface restudies involved in the interaction.

Unraveling the Role of Charge Patterning in the Micellar Structure of Sequence-Defined Amphiphilic Peptoid Oligomers by Molecular Dynamics Simulations

Electrostatic interactions play a significant role in regulating biological systems and have received increasing attention due to their usefulness in designing advanced stimulus-responsive materials. Polypeptoids are highly tunable N-substituted peptidomimetic polymers that lack backbone hydrogen bonding and chirality. Therefore, polypeptoids are suitable systems to study the effect of noncovalent interactions of substituents without complications of backbone intramolecular and intermolecular hydrogen bonding. In this study, all-atom molecular dynamics (MD) simulations were performed on micelles formed by a series of sequence-defined ionic polypeptoid block copolymers consisting of a hydrophobic segment and a hydrophilic segment in an aqueous solution. By combining the results from MD simulations and experimental small-angle neutron scattering data, further insights were gained into the internal structure of the formed polypeptoid micelles, which is not always directly accessible from experiments. In addition, information was gained into the physics of the noncovalent interactions responsible for the self-assembly of weakly charged polypeptoids in an aqueous solution. While the aggregation number is governed by electrostatic repulsion of the negatively charged carboxylate (COO^–) substituents on the polypeptoid chain within the micelle, MD simulations indicate that the position of the charge on singly charged chains mediates the shape of the micelle through the charge–dipole interactions between the COO^– substituent and the surrounding water. Therefore, the polypeptoid micelles formed from the single-charged series offer the possibility for tailorable micelle shapes. In contrast, the polypeptoid micelles formed from the triple-charged series are characterized by more pronounced electrostatic repulsion that competes with more significant charge–sodium interactions, making it difficult to predict the shape of the micelles. This work has helped further develop design principles for the shape and structure of self-assembled micelles by controlling the position of charged moieties on the backbone of polypeptoid block copolymers.

A Bottom-Up Coarse-Grained Model for Nucleosome-Nucleosome Interactions with Explicit Ions

The nucleosome core particle (NCP) is a large complex of 145-147 base pairs of DNA and eight histone proteins and is the basic building block of chromatin that forms the chromosomes. Here, we develop a coarse-grained (CG) model of the NCP derived through a systematic bottom-up approach based on underlying all-atom MD simulations to compute the necessary CG interactions. The model produces excellent agreement with known structural features of the NCP and gives a realistic description of the nucleosome-nucleosome attraction in the presence of multivalent cations (Mg(H₂O)₆²⁺ or Co(NH₃)₆³⁺) for systems comprising 20 NCPs. The results of the simulations reveal structural details of the NCP-NCP interactions unavailable from experimental approaches, and this model opens the prospect for the rigorous modeling of chromatin fibers.

Modulation of DNA conformation in electrolytic nanodroplets

The behavior of deoxyribonucleic acid (DNA) molecules in confinement is of profound importance in various bioengineering and medical applications. In the present study, all-atom molecular dynamics simulation is utilized to investigate the transition of the double-strand DNA (dsDNA) conformation in the electrolytic nanodroplet. Three typical conformations, i.e., C-shaped, folded S-shaped, and double C-shaped, are observed for different droplet sizes and ionic concentrations. To reveal the physics underlying this phenomenon, the characteristics of the dsDNA molecules, such as the overcharging intensity, the end-to-end distance, the radius of gyration, etc. are analyzed in detail based on the numerical results. It is found that the transition can be ascribed to the buckling of the polymer molecules under the compression due to the confinement of the nanodroplet, and it can be modulated by the ionic concentration in the electrolyte. Generally, nanoscale confinement dominates dsDNA behavior over the electrostatic effects in smaller nanodroplets, while the latter becomes more important for larger nanodroplets. This competition results in the persistence length increasing with the nanodroplet radii. Based on these discussions, a non-dimensional elasto-capillary number μ is proposed to classify the dsDNA conformations into three regions.

Single Molecular Chelation Dynamics Reveals That DNA Has a Stronger Affinity toward Lead(II) than Cadmium(II)

Lead ions can bind to DNA via nonelectrostatic interactions and hence alter its structure, which may be related to their adverse effects. The dynamics of Pb²⁺-DNA interaction has not been well understood. In this study, we report the monomolecular dynamics of the Pb²⁺-DNA interaction using a magnetic tweezers (MT) setup. We found that lead cations could induce DNA compaction at ionic strengths above 1 μM, which was also confirmed by morphology characterization. The chelation behavior of the Pb²⁺-DNA and the Cd²⁺-DNA complex solutions after adding EDTA were compared. The results showed that EDTA chelated with the bound metal ions on DNA and consequently led to restoring the DNA to its original length but with different restoration speeds for the two solutions. The fast binding dynamics and the slower chelation dynamics of the Pb²⁺ scenario compared to that of Cd²⁺ suggested that Pb²⁺ was more capable to induce DNA conformational change and that the Pb²⁺-DNA complex was more stable than the Cd²⁺-DNA complex. The stronger affinities for DNA bases and the inner binding of lead cations were two possible causes of the dynamics differences. Three agents, including EDTA, sodium gluconate, and SDBS, were used to remove the bound lead ions on DNA. It was shown that EDTA was the most efficient, and sodium gluconate could not fully restore DNA from its compact state. We concluded that both EDTA and SDBS were good candidates to restore the Pb²⁺-bound DNA to its original state.

A fluorescent reporter on electrostatic DNA-ligand interactions

Among the various types of interactions between biomolecules, electrostatic interactions dominate as these are long-range interactions and are often a generic first step in the recruitment of specific ligands. DNA, being a highly charged molecule, attracts a plethora of molecules. Interactions between DNA and proteins or small molecules shape the overall function of the cell. Various processes such as DNA replication, DNA repair, synthesis of mRNA, and packaging of DNA are mediated by interactions between protein molecules and DNA that are predominantly electrostatic. Here, we present a fluorescence resonance energy transfer (FRET)-based probe which can report on the electrostatic interactions between the negatively-charged DNA and positively-charged metal ions, oligopeptides, as well as DNA groove-binding drug molecules. The simplicity, sensitivity, and versatility of the DNA-based probe makes it suited for applications where specific protein-DNA interactions can be probed, and DNA-binding drugs can be discovered in high-throughput screens of compound libraries. This is particularly relevant given that some of the most potent antitumor and antimicrobial drugs associate with DNA electrostatically.

Electrostatic co-assembly of pillar[n]pyridiniums and calix[4]arene in aqueous media

The cationic pillar[n]pyridiniums and anionic p-sulfonatocalix[4]arene co-assemble into all-organic supersalts through encaging of the supercation units within/between the capsules emerged from superanion pairs. The encapsulation occures both in the solid...

Functional DNA-Polymer Conjugates

DNA nanotechnology has seen large developments over the last 30 years through the combination of solid phase synthesis and the discovery of DNA nanostructures. Solid phase synthesis has facilitated the availability of short DNA sequences and the expansion of the DNA toolbox to increase the chemical functionalities afforded on DNA, which in turn enabled the conception and synthesis of sophisticated and complex 2D and 3D nanostructures. In parallel, polymer science has developed several polymerization approaches to build di- and triblock copolymers bearing hydrophilic, hydrophobic, and amphiphilic properties. By bringing together these two emerging technologies, complementary properties of both materials have been explored; for example, the synthesis of amphiphilic DNA-polymer conjugates has enabled the production of several nanostructures, such as spherical and rod-like micelles. Through both the DNA and polymer parts, stimuli-responsiveness can be instilled. Nanostructures have consequently been developed with responsive structural changes to physical properties, such as pH and temperature, as well as short DNA through competitive complementary binding. These responsive changes have enabled the application of DNA-polymer conjugates in biomedical applications including drug delivery. This review discusses the progress of DNA-polymer conjugates, exploring the synthetic routes and state-of-the-art applications afforded through the combination of nucleic acids and synthetic polymers.

Measuring structural parameters of crosslinked and entangled semiflexible polymer networks with single-filament tracing

Single-filament tracing has been a valuable tool to directly determine geometrical and mechanical properties of entangled polymer networks. However, systematically verifying how the stiffness of the tracer filament or its molecular interactions with the surrounding network impacts the measurement of these parameters has not been possible with the established experimental systems. Here we use mechanically programmable DNA nanotubes embedded in crosslinked and entangled F-actin networks, as well as in synthetic DNA networks, in order to measure fundamental, structural network properties like tube width and mesh size with respect to the stiffness of the tracers. While we confirm some predictions derived from models based purely on steric interactions, our results indicate that these models should be expanded to account for additional interfilament interactions, thereby describing the behavior of real polymer networks.

Anomalous Laterally Stressed Kinetically Trapped DNA Surface Conformations

DNA kinking is inevitable for the highly anisotropic 1D-1D electrostatic interaction with the one-dimensionally periodically charged surface. The double helical structure of the DNA kinetically trapped on positively charged monomolecular films comprising the lamellar templates is strongly laterally stressed and extremely perturbed at the nanometer scale. The DNA kinetic trapping is not a smooth 3D-> 2D conformational flattening but is a complex nonlinear in-plane mechanical response (bending, tensile and unzipping) driven by the physics beyond the scope of the applicability of the linear worm-like chain approximation. Up to now, the DNA molecule adsorbed on a surface was believed to always preserve its native structure. This belief implies a negligible contribution of lateral surface forces during and after DNA adsorption although their impact has never been elucidated. High-resolution atomic force microscopy was used to observe that stiff DNA molecules kinetically trapped on monomolecular films comprising one-dimensional periodically charged lamellar templates as a single layer or as a sublayer are oversaturated by sharp discontinuous kinks and can also be locally melted and supercoiled. We argue that kink/anti-kink pairs are induced by an overcritical lateral bending stress (> 30 pNnm) inevitable for the highly anisotropic 1D-1D electrostatic interaction of DNA and underlying rows of positive surface charges. In addition, the unexpected kink-inducing mechanical instability in the shape of the template-directed DNA confined between the positively charged lamellar sides is observed indicating the strong impact of helicity. The previously reported anomalously low values of the persistence length of the surface-adsorbed DNA are explained by the impact of the surface-induced low-scale bending. The sites of the local melting and supercoiling are convincingly introduced as other lateral stress-induced structural DNA anomalies by establishing a link with DNA high-force mechanics. The results open up the study in the completely unexplored area of the principally anomalous kinetically trapped DNA surface conformations in which the DNA local mechanical response to the surface-induced spatially modulated lateral electrostatic stress is essentially nonlinear. The underlying rich and complex in-plane nonlinear physics acts at the nanoscale beyond the scope of applicability of the worm-like chain approximation.

Engineering Hydrogel Adhesion for Biomedical Applications via Chemical Design of the Junction

Hydrogel adhesion inherently relies on engineering the contact surface at soft and hydrated interfaces. Upon contact, adhesion normally occurs through the formation of chemical or physical interactions between the disparate surfaces. The ability to form these adhesion junctions is challenging for hydrogels as the interfaces are wet and deformable and often contain low densities of functional groups. In this Review, we link the design of the binding chemistries or adhesion junctions, whether covalent, dynamic covalent, supramolecular, or physical, to the emergent adhesive properties of soft and hydrated interfaces. Wet adhesion is useful for bonding to or between tissues and implants for a range of biomedical applications. We highlight several recent and emerging adhesive hydrogels for use in biomedicine in the context of efficient junction design. The main focus is on engineering hydrogel adhesion through molecular design of the junctions to tailor the adhesion strength, reversibility, stability, and response to environmental stimuli.

Buckling transitions and soft-phase invasion of two-component icosahedral shells

What is the optimal distribution of two types of crystalline phases on the surface of icosahedral shells, such as of many viral capsids? We here investigate the distribution of a thin layer of soft material on a crystalline convex icosahedral shell. We demonstrate how the shapes of spherical viruses can be understood from the perspective of elasticity theory of thin two-component shells. We develop a theory of shape transformations of an icosahedral shell upon addition of a softer, but still crystalline, material onto its surface. We show how the soft component "invades" the regions with the highest elastic energy and stress imposed by the 12 topological defects on the surface. We explore the phase diagram as a function of the surface fraction of the soft material, the shell size, and the incommensurability of the elastic moduli of the rigid and soft phases. We find that, as expected, progressive filling of the rigid shell by the soft phase starts from the most deformed regions of the icosahedron. With a progressively increasing soft-phase coverage, the spherical segments of domes are filled first (12 vertices of the shell), then the cylindrical segments connecting the domes (30 edges) are invaded, and, ultimately, the 20 flat faces of the icosahedral shell tend to be occupied by the soft material. We present a detailed theoretical investigation of the first two stages of this invasion process and develop a model of morphological changes of the cone structure that permits noncircular cross sections. In conclusion, we discuss the biological relevance of some structures predicted from our calculations, in particular for the shape of viral capsids.

Recent Progress in Ionic Coassembly of Cationic Peptides and Anionic Species

Peptide assembly has been extensively exploited as a promising platform for the creation of hierarchical nanostructures and tailor-made bioactive materials. Ionic coassembly of cationic peptides and anionic species is paving the way to provide particularly important contribution to this topic. In this review, the recent progress of ionic coassembly soft materials derived from the electrostatic coupling between cationic peptides and anionic species in aqueous solution is systematically summarized. The presentation of this review starts from a brief background on the general importance and advantages of peptide-based ionic coassembly. After that, diverse combinations of cationic peptides with small anions, macro- and/or oligo-anions, anionic polymers, and inorganic polyoxometalates are described. Emphasis is placed on the hierarchical structures, value-added properties, and applications. The molecular design of cationic peptides and the general principles behind the ionic coassembled structures are discussed. It is summarized that the combination of interesting and unique characteristics that arise both from the chemical diversity of peptides and the wide range of anionic species may contribute in a variety of output, including drug delivery, tissue engineering, gene transfection, and antibacterial activity. The emergent new phenomena and findings are illustrated. Finally, the outlook for the peptide-based ionic coassembly systems is also presented.

Divalent metal ions and intermolecular interactions facilitate DNA network formation

The interactions between divalent metal ions and DNA are crucial for basic life processes. These interactions are also important in advanced technological products such as DNA-based ion sensors. Current polyelectrolyte theories cannot describe these interactions well and do not consider the corresponding dynamics. In this study, we report the single-molecule dynamics of the binding of divalent metal ions to a single DNA molecule and the morphology characterization of the complex. We found that most of the divalent metal ions (Mn²⁺, Zn²⁺, Co²⁺, Ni²⁺, and Cd²⁺), except Mg²⁺ and Ca²⁺, could cause monomolecular DNA condensation. For transition metal ions, different ionic strengths were required to induce the compaction, and different shortening speeds were displayed in the dynamics, indicating ionic specificity. Atomic force microscopy revealed that the morphologies of the metal ion-DNA complexes were affected by the ionic strength of the metal ion, DNA chain length, and DNA concentration. At low metal ion concentration, DNA tended to adopt a random coil conformation. Increasing the ionic strength led to network-like condensed structures, suggesting that divalent metal ions can induce attraction between DNA molecules. Furthermore, higher DNA concentration and longer chain length enhanced intermolecular interactions and consequently resulted in network structures with a higher degree of interconnectivity.

Critical adsorption of multiple polyelectrolytes onto a nanosphere: splitting the adsorption-desorption transition boundary

Employing extensive Monte Carlo computer simulations, we investigate in detail the properties of multichain adsorption of charged flexible polyelectrolytes (PEs) onto oppositely charged spherical nanoparticles (SNPs). We quantify the conditions of critical adsorption-the phase-separation curve between the adsorbed and desorbed states of the PEs-as a function of the SNP surface-charge density and the concentration of added salt. We study the degree of fluctuations of the PE-SNP electrostatic binding energy, which we use to quantify the emergence of the phase subtransitions, including a series of partially adsorbed PE configurations. We demonstrate how the phase-separation adsorption-desorption boundary shifts and splits into multiple subtransitions at low-salt conditions, thereby generalizing and extending the results for critical adsorption of a single PE onto the SNP. The current findings are relevant for finite concentrations of PEs around the attracting SNP, such as the conditions for PE adsorption onto globular proteins carrying opposite electric charges.

Multivalent counterions induced attraction between DNA polyelectrolytes

In this paper we study the electrostatic attraction between two parallel rodlike DNA polyelectrolytes induced by neutralizing multivalent counterions at the zero temperature limit. The counterions crystallize on the charged surfaces of DNA so that we can handle the system by using the Wigner crystal lattice model. We derived the 3D ground state configuration of counterions with minimized energy by use of the gradient descent method, and calculated the interaction between two DNA cylinders with divalent or trivalent counterions when they approach. The results show that the complex ground state configuration of counterions plays a key role in determining the caused attraction. The counterions form three-dimensional Wigner crystals on each cylinder at large separation. When the cylinders are brought together, some counterion lines will move towards the inner region and lead to strong attraction. The calculated interaction from our model is in good agreement with the simulation result, however, the single particle approximation considerably overestimates the attraction.

Counterions in the ground state partially condense in the intervening region of two approaching DNA cylinders and result in attraction.

Rational design of a DNA sequence-specific modular protein tag by tuning the alkylation kinetics

Sequence-selective chemical modification of DNA by synthetic ligands has been a long-standing challenge in the field of chemistry. Even when the ligand consists of a sequence-specific DNA binding domain and reactive group, sequence-selective reactions by these ligands are often accompanied by off-target reactions. A basic principle to design DNA modifiers that react at specific sites exclusively governed by DNA sequence recognition remains to be established. We have previously reported selective DNA modification by a self-ligating protein tag conjugated with a DNA-binding domain, termed as a modular adaptor, and orthogonal application of modular adaptors by relying on the chemoselectivity of the protein tag. The sequence-specific crosslinking reaction by the modular adaptor is thought to proceed in two steps: the first step involves the formation of a DNA-protein complex, while in the second step, a proximity-driven intermolecular crosslinking occurs. According to this scheme, the specific crosslinking reaction of a modular adaptor would be driven by the DNA recognition process only when the dissociation rate of the DNA complex is much higher than the rate constant for the alkylation reaction. In this study, as a proof of principle, a set of combinations for modular adaptors and their substrates were utilized to evaluate the reactions. Three types of modular adaptors consisting of a single type of self-ligating tag and three types of DNA binding proteins fulfill the kinetic requirements for the reaction of the self-ligating tag with a substrate and the dissociation of the DNA-protein complex. These modular adaptors actually undergo sequence-specific crosslinking reactions exclusively driven by the recognition of a specific DNA sequence. The design principle of sequence-specific modular adaptors based on the kinetic aspects of complex formation and chemical modification is applicable for developing recognition-driven selective modifiers for proteins and other biological macromolecules.

Multistage dynamics of Hg2+-DNA interactions: a single-molecule study

The metal ion-DNA interaction is key to biochemical processes and has applications in areas such as metal ion sensors and DNA nanomachines. For example, the formation of the T-Hg2+-T structure has been used in technologies such as DNA-based mercuric ion sensors. Though the interaction is widely used for practical purposes, the underlying mechanism has not been fully understood. In the present study, we used magnetic tweezers to explore the interactions between λ-DNA and two metal ions, Hg2+ and Cd2+, at the single-molecule level. Both metal ions caused considerable DNA conformational changes. The resulting DNA compaction dynamics were related to the ion concentration and the exerted force. The increase in the ion concentration promoted DNA compaction, whereas exerting greater forces inhibited this process. Application of a high force generated two-stage dynamics of the Hg2+-DNA interaction. However, at a sufficiently high Hg2+ concentration, a lower force led to a three-stage process. In contrast, the curves of the binding of Cd2+ ions to DNA had a stepwise pattern. Both the AFM scanning results and the single-molecule measurements confirmed that Hg2+ influences the DNA conformation in a more pronounced manner than Cd2+. The multistage Hg2+-DNA interaction was considered to be a result of the different binding mechanisms, including the mismatched base-pair formation. A model was then proposed to explain the peculiar dynamics.

Off-center charge model revisited: Electrical double layer with multivalent cations

The off-center charge model of ions is a relatively simple model for introducing asymmetry in Coulomb interaction while retaining the simplicity and convenience of the spherical hard core geometry. A Monte Carlo simulation analysis of the planar electric double layer formed by this ionic model for 1+:1- valence systems [S. Lamperski et al., Langmuir 33, 11554-11560 (2017)] is extended to include solutions of multivalent (2+, 3+) hard spherical cations and single valence (1-) hard spherical anions near a uniformly charged, planar electrode. The solvent is modelled as a uniform dielectric continuum with a dielectric constant equal to that of the pure solvent, viz., the primitive model. Results are reported for the ion density, the cation charge profile, and the electrostatic potential profile at 1 mol/dm³ salt concentration. Additionally, the double layer potential drop, that is, the electrode potential, and the integral and the differential capacitances are computed as functions of the electrode surface charge density. The latter two quantities show an expected asymmetry as long as the cation valence is not too great and the charge of the off-center ion cannot approach too close to the electrode surface. It is unusual that the integral and differential capacitances are negative for high valence cations and a negatively charged electrode when the off-center charge is large and can be very near the surface of the electrode. The corresponding electrode potential versus surface charge density curve becomes non-monotonic and shows a change of slope, and thus the resultant integral and differential capacitances can become negative. This nonphysical result is the result of an incipient singularity when a large positive charge is too near a negatively charged electrode. Overall, the off-center charge model suggests a useful recipe to model electrical asymmetry within the broader context of the primitive model provided that the off-center charge is not too near the surface of the electrode.

Open-Boundary Molecular Dynamics of a DNA Molecule in a Hybrid Explicit/Implicit Salt Solution

The composition and electrolyte concentration of the aqueous bathing environment have important consequences for many biological processes and can profoundly affect the behavior of biomolecules. Nevertheless, because of computational limitations, many molecular simulations of biophysical systems can be performed only at specific ionic conditions: either at nominally zero salt concentration, i.e., including only counterions enforcing the system's electroneutrality, or at excessive salt concentrations. Here, we introduce an efficient molecular dynamics simulation approach for an atomistic DNA molecule at realistic physiological ionic conditions. The simulations are performed by employing the open-boundary molecular dynamics method that allows for simulation of open systems that can exchange mass and linear momentum with the environment. In our open-boundary molecular dynamics approach, the computational burden is drastically alleviated by embedding the DNA molecule in a mixed explicit/implicit salt-bathing solution. In the explicit domain, the water molecules and ions are both overtly present in the system, whereas in the implicit water domain, only the ions are explicitly present and the water is described as a continuous dielectric medium. Water molecules are inserted and deleted into/from the system in the intermediate buffer domain that acts as a water reservoir to the explicit domain, with both water molecules and ions free to enter or leave the explicit domain. Our approach is general and allows for efficient molecular simulations of biomolecules solvated in bathing salt solutions at any ionic strength condition.

The Mixing Counterion Effect on DNA Compaction and Charge Neutralization at Low Ionic Strength

DNA compaction and charge neutralization in a mixing counterion solution involves competitive and cooperative electrostatic binding, and sometimes counterion complexation. At normal ionic strength, it has been found that the charge neutralization of DNA by the multivalent counterion is suppressed when being added extra mono- and di-valent counterions. Here, we explore the effect mixing counterion on DNA compaction and charge neutralization under the condition of low ionic strength. Being quite different from normal ionic strength, the electrophoretic mobility of DNA in multivalent counterion solution (octalysine, spermine) increases the presence of mono- and di-valent cations, such as sodium and magnesium ions. It means that the charge neutralization of DNA by the multivalent counterion is promoted rather than suppressed when introducing extra mono- and di-valent counterions into solution. This conclusion is also supported by the measurement of condensing and unraveling forces of DNA condensates under the same condition by single molecular magnetic tweezers. This mixing effect can be attributed to the cooperative electrostatic binding of counterions to DNA when the concentration of counterions in solution is below a critical concentration.

A systematic analysis of nucleosome core particle and nucleosome-nucleosome stacking structure

Chromatin condensation is driven by the energetically favourable interaction between nucleosome core particles (NCPs). The close NCP-NCP contact, stacking, is a primary structural element of all condensed states of chromatin in vitro and in vivo. However, the molecular structure of stacked nucleosomes as well as the nature of the interactions involved in its formation have not yet been systematically studied. Here we undertake an investigation of both the structural and physico-chemical features of NCP structure and the NCP-NCP stacking. We introduce an “NCP-centred” set of parameters (NCP-NCP distance, shift, rise, tilt, and others) that allows numerical characterisation of the mutual positions of the NCPs in the stacking and in any other structures formed by the NCP. NCP stacking in more than 140 published NCP crystal structures were analysed. In addition, coarse grained (CG) MD simulations modelling NCP condensation was carried out. The CG model takes into account details of the nucleosome structure and adequately describes the long range electrostatic forces as well as excluded volume effects acting in chromatin. The CG simulations showed good agreement with experimental data and revealed the importance of the H2A and H4 N-terminal tail bridging and screening as well as tail-tail correlations in the stacked nucleosomes.

Acceleration of bursty multiprotein target search kinetics on DNA by colocalisation

Proteins are capable of locating specific targets on DNA by employing a facilitated diffusion process with intermittent 1D and 3D search steps. We here uncover the implications of colocalisation of protein production and DNA binding sites via computer simulations.

Potential of mean force between oppositely charged nanoparticles: A comprehensive comparison between Poisson-Boltzmann theory and Monte Carlo simulations

Ion-mediated interactions between like-charged polyelectrolytes have been paid much attention, and the Poisson-Boltzmann (PB) theory has been shown to fail in qualitatively predicting multivalent ion-mediated like-charge attraction. However, inadequate attention has been paid to the ion-mediated interactions between oppositely charged polyelectrolytes. In this work, the potentials of mean force (PMF) between oppositely charged nanoparticles in 1:1 and 2:2 salt solutions were investigated by Monte Carlo simulations and the PB theory. Our calculations show that the PMFs between oppositely charged nanoparticles are generally attractive in 1:1 and 2:2 salt solutions and that such attractive PMFs become weaker at higher 1:1 or 2:2 salt concentrations. The comprehensive comparisons show that the PB theory can quantitatively predict the PMFs between oppositely charged nanoparticles in 1:1 salt solutions, except for the slight deviation at very high 1:1 salt concentration. However, for 2:2 salt solutions, the PB theory generally overestimates the attractive PMF between oppositely charged nanoparticles, and this overestimation becomes more pronounced for nanoparticles with higher charge density and for higher 2:2 salt concentration. Our microscopic analyses suggest that the overestimation of the PB theory on the attractive PMFs for 2:2 salt solutions is attributed to the underestimation of divalent ions bound to nanoparticles.

Divalent Ion-Mediated DNA-DNA Interactions: A Comparative Study of Triplex and Duplex

Ion-mediated interaction between DNAs is essential for DNA condensation, and it is generally believed that monovalent and nonspecifically binding divalent cations cannot induce the aggregation of double-stranded (ds) DNAs. Interestingly, recent experiments found that alkaline earth metal ions such as Mg²⁺ can induce the aggregation of triple-stranded (ts) DNAs, although there is still a lack of deep understanding of the surprising findings at the microscopic level. In this work, we employed all-atom dynamic simulations to directly calculate the potentials of mean force (PMFs) between tsDNAs, between dsDNAs, and between tsDNA and dsDNA in Mg²⁺ solutions. Our calculations show that the PMF between tsDNAs is apparently attractive and becomes more strongly attractive at higher [Mg²⁺], although the PMF between dsDNAs cannot become apparently attractive even at high [Mg²⁺]. Our analyses show that Mg²⁺ internally binds into grooves and externally binds to phosphate groups for both tsDNA and dsDNA, whereas the external binding of Mg²⁺ is much stronger for tsDNA. Such stronger external binding of Mg²⁺ for tsDNA favors more apparent ion-bridging between helices than for dsDNA. Furthermore, our analyses illustrate that bridging ions, as a special part of external binding ions, are tightly and positively coupled to ion-mediated attraction between DNAs.

Order and interactions in DNA arrays: Multiscale molecular dynamics simulation

While densely packed DNA arrays are known to exhibit hexagonal and orthorhombic local packings, the detailed mechanism governing the associated phase transition remains rather elusive. Furthermore, at high densities the atomistic resolution is paramount to properly account for fine details, encompassing the DNA molecular order, the contingent ordering of counterions and the induced molecular ordering of the bathing solvent, bringing together electrostatic, steric, thermal and direct hydrogen-bonding interactions, resulting in the observed osmotic equation of state. We perform a multiscale simulation of dense DNA arrays by enclosing a set of 16 atomistically resolved DNA molecules within a semi-permeable membrane, allowing the passage of water and salt ions, and thus mimicking the behavior of DNA arrays subjected to external osmotic stress in a bathing solution of monovalent salt and multivalent counterions. By varying the DNA density, local packing symmetry, and counterion type, we obtain osmotic equation of state together with the hexagonal-orthorhombic phase transition, and full structural characterization of the DNA subphase in terms of its positional and angular orientational fluctuations, counterion distributions, and the solvent local dielectric response profile with its order parameters that allow us to identify the hydration force as the primary interaction mechanism at high DNA densities.