Search

We report here an experimental-computational study of hydrated TiO2anatase nanoparticles interacting with glycine, where we obtain quantitative agreement of the measured adsorption free energies. Ab initio simulations are performed within the tight binding and density functional theory in combination with enhanced free-energy sampling techniques, which exploit the thermodynamic integration of the unbiased mean forces collected on-the-fly along the molecular dynamics trajectories. The experiments adopt a new and efficient setup for electrochemical impedance spectroscopy measurements based on portable screen-printed gold electrodes, which allows fast and in situ signal assessment. The measured adsorption free energy is −30 kJ/mol (both from experiment and calculation), with preferential interaction of the charged NH3+group which strongly adsorbs on the TiO2bridging oxygens. This highlights the importance of the terminal amino groups in the adsorption mechanism of amino acids on hydrated metal oxides. The excellent agreement between computation and experiment for this amino acid opens the doors to the exploration of the interaction free energies for other moderately complex bionano systems.

We propose a method of bottom-up coarse-graining, in which interactions within a coarse-grained model are determined by an artificial neural network trained on structural data obtained from multiple atomistic simulations. The method uses ideas of the inverse Monte Carlo approach, relating changes in the neural network weights with changes in average structural properties, such as radial distribution functions. As a proof of concept, we demonstrate the method on a system interacting by a Lennard–Jones potential modeled by a simple linear network and a single-site coarse-grained model of methanol–water solutions. In the latter case, we implement a nonlinear neural network with intermediate layers trained by atomistic simulations carried out at different methanol concentrations. We show that such a network acts as a transferable potential at the coarse-grained resolution for a wide range of methanol concentrations, including those not included in the training set.

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(H2O)62+or Co(NH3)63+) 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.

Understanding of interactions between inorganic nanomaterials and biomolecules, and particularly lipid bilayers, is crucial in many biotechnological and biomedical applications, as well as for the evaluation of possible toxic effects caused by nanoparticles. Here, we present a molecular dynamics study of adsorption of two important constituents of the cell membranes, 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), lipids to a number of titanium dioxide planar surfaces, and a spherical nanoparticle under physiological conditions. By constructing the number density profiles of the lipid headgroup atoms, we have identified several possible binding modes and calculated their relative prevalence in the simulated systems. Our estimates of the adsorption strength, based on the total fraction of adsorbed lipids, show that POPE binds to the selected titanium dioxide surfaces stronger than DMPC, due to the ethanolamine group forming hydrogen bonds with the surface. Moreover, while POPE shows a clear preference toward anatase surfaces over rutile, DMPC has a particularly high affinity to rutile(101) and a lower affinity to other surfaces. Finally, we study how lipid concentration, addition of cholesterol, as well as titanium dioxide surface curvature may affect overall adsorption.

Membrane curvature sensing, where the binding free energies of membrane-associated molecules depend on the local membrane curvature, is a key factor to modulate and maintain the shape and organization of cell membranes. However, the microscopic mechanisms are not well understood, partly due to absence of efficient simulation methods. Here, we describe a method to compute the curvature dependence of the binding free energy of a membrane-associated probe molecule that interacts with a buckled membrane, which has been created by lateral compression of a flat bilayer patch. This buckling approach samples a wide range of curvatures in a single simulation, and anisotropic effects can be extracted from the orientation statistics. We develop an efficient and robust algorithm to extract the motion of the probe along the buckled membrane surface, and evaluate its numerical properties by extensive sampling of three coarse-grained model systems: local lipid density in a curved environment for single-component bilayers, curvature preferences of individual lipids in two-component membranes, and curvature sensing by a homotrimeric transmembrane protein. The method can be used to complement experimental data from curvature partition assays and provides additional insight into mesoscopic theories and molecular mechanisms for curvature sensing.

Adsorption free energies of 32 small biomolecules (amino acids side chains, fragments of lipids, and sugar molecules) on 33 different nanomaterials, computed by the molecular dynamics - metadynamics methodology, have been analyzed using statistical machine learning approaches. Multiple unsupervised learning algorithms (principal component analysis, agglomerative clustering, and K-means) as well as supervised linear and nonlinear regression algorithms (linear regression, AdaBoost ensemble learning, artificial neural network) have been applied. As a result, a small set of biomolecules has been identified, knowledge of adsorption free energies of which to a specific nanomaterial can be used to predict, within the developed machine learning model, adsorption free energies of other biomolecules. Furthermore, the methodology of grouping of nanomaterials according to their interactions with biomolecules has been presented.

Due to the vast length scale inside the cell nucleus, multiscale models are required to understand chromatin folding, structure, and dynamics and how they regulate genomic activities such as DNA transcription, replication, and repair. We study the interactions and structure of condensed phases formed by the universal building block of chromatin, the nucleosome core particle (NCP), using bottom-up multiscale coarse-grained (CG) simulations with a model extracted from all-atom MD simulations. In the presence of the multivalent cations Mg(H2O)62+or CoHex3+, we analyze the internal structures of the NCP aggregates and the contributions of histone tails and ions to the aggregation patterns. We then derive a "super" coarse-grained (SCG) NCP model to study the macroscopic scale phase separation of NCPs. The SCG simulations show the formation of NCP aggregates with Mg(H2O)62+concentration-dependent densities and sizes. Variation of the CoHex3+concentrations results in highly ordered lamellocolumnar and hexagonal columnar phases in agreement with experimental data. The results give detailed insights into nucleosome interactions and for understanding chromatin folding in the cell nucleus.

The use of carbon-based nanomaterials is tremendously increasing in various areas of technological, bioengineering, and biomedical applications. The functionality of carbon-based nanomaterials can be further broadened viachemical functionalization of carbon nanomaterial surfaces. On the other hand, concern is rising on possible adverse effects when nanomaterials are taken up by biological organisms. In order to contribute into understanding of interactions of carbon-based nanomaterials with biological matter, we have investigated adsorption of small biomolecules on nanomaterials using enhanced sampling molecular dynamics. The biomolecules included amino acid side chain analogues, fragments of lipids, and sugar monomers. The adsorption behavior on unstructured amorphous carbon, pristine graphene and its derivatives (such as few-layer graphene, graphene oxide, and reduced graphene oxide) as well as pristine carbon nanotubes, and those functionalized with OH–, COOH–, COO–, NH2–, and NH3+groups was investigated with respect to surface concentration. An adsorption profile, that is, the free energy as a function of distance from the nanomaterial surfaces, was determined for each molecule and surface using the Metadynamics approach. The results were analyzed in terms of chemical specificity, surface charge, and surface concentration. It was shown that although morphology of the nanomaterial has a limited effect on the adsorption properties, functionalization of the surface by various molecular groups can drastically change the adsorption behavior that can be used in the design of nanosurfaces with highly selective adsorption properties and safe for human health and environment.

The molecular mechanics force field Slipids developed in a series of works by Jämbeck and Lyubartsev (J. Phys. Chem. B2012, 116, 3164–3179; J. Chem. Theory Comput.2012, 8, 2938–2948) generally provides a good description of various lipid bilayer systems. However, it was also found that order parameters of C–H bonds in the glycerol moiety of the phosphatidylcholine headgroup deviate significantly from NMR results. In this work, the dihedral force field parameters have been reparameterized in order to improve the agreement with experiment. For this purpose, we have computed energies for a large amount of lipid headgroup conformations using density functional theory on the B3P86/cc-pvqz level and optimized dihedral angle parameters simultaneously to provide the best fit to the quantum chemical energies. The new parameter set was validated for three lipid bilayer systems against a number of experimental properties including order parameters, area per lipid, scattering form factors, bilayer thickness, area compressibility and lateral diffusion coefficients. In addition, the order parameter dependence on cholesterol content in the POPC bilayer was investigated. It is shown that the new force field significantly improves agreement with the experimental order parameters for the lipid headgroup while keeping good agreement with other experimentally measured properties.

Accurate parametrization of force fields (FFs) is of ultimate importance for computer simulations to be reliable and to possess a predictive power. In this work, we analyzed, in multi-microsecond simulations of a 40-base-pair DNA fragment, the performance of four force fields, namely, the two recent major updates of CHARMM and two from the AMBER family. We focused on a description of double-helix DNA flexibility and dynamics both at atomistic and at mesoscale level in coarse-grained (CG) simulations. In addition to the traditional analysis of different base-pair and base-step parameters, we extended our analysis to investigate the ability of the force field to parametrize a CG DNA model by structure-based bottom-up coarse-graining, computing DNA persistence length as a function of ionic strength. Our simulations unambiguously showed that the CHARMM36 force field is unable to preserve DNA’s structural stability at over-microsecond time scale. Both versions of the AMBER FF, parmbsc0 and parmbsc1, showed good agreement with experiment, with some bias of parmbsc0 parameters for intermediate A/B form DNA structures. The CHARMM27 force field provides stable atomistic trajectories and overall (among the considered force fields) the best fit to experimentally determined DNA flexibility parameters both at atomistic and at mesoscale level.

Dissociative and molecular water adsorption on the anatase (001) surface is studied in the context of state-of-the-art density functional theory in large supercells suited for adsorption studies at various water coverage ratios. At low coverage values below 1/4 ML, water adsorption remains dissociative and a network of hydrogen bonds between the so formed hydroxyl groups favors the formation of a ridge surface structure. The hydroxyl patterned (4 × 4) surface thus undergoes a (2 × 4) reconstruction that causes the relief of the large tensile stress measured in the unreconstructed surface along the direction orthogonal to the ridge. This phenomenology is accompanied by the loss of reactivity of the reconstructed surface with respect to the dissociative water adsorption that becomes molecular above 1/4 ML. We also show that the molecular adsorption on the terrace is weaker than the one on the ridge. The present water reconstruction model is discussed and compared to the well-known ADM model of the reconstructed anatase (001) surface in dry environment.

In this work, we present results from molecular dynamics simulations of aqueous solutions of furfural and 5-hydroxymethylfurfural, which are important intermediates in the hydrothermal carbonization processes of biomass conversion. The computations were performed both at ambient and hydrothermal conditions using a two-level factorial design varying concentration, temperature, and pressure. A number of equilibrium and dynamic properties have been computed including enthalpies and free energies of vaporization, free energies of solvation, diffusion coefficients, and rotational/reorientational correlation times. Structural properties of solutions were analyzed using radial and spatial distribution functions. It was shown that the formation of hydrogen bonds among 5-hydroxymethylfurfural molecules is preferred compared to hydrogen bonding between 5-hydroxymethylfurfural and water. In addition, our results suggest that the oxygen atoms in the furan rings of furfural and 5-hydroxymethylfurfural do not participate in hydrogen bonding to the same extent as the oxygen atoms in the hydroxyl and carbonyl groups. It is also observed that furfural molecules aggregate under certain conditions, and we show how this is affected by changes in temperature, pressure, and concentration in agreement with experimental solubility data. The analysis of the computational results provides useful insight into the structure and dynamics of the considered molecules at conditions of hydrothermal carbonization, as well as at ambient conditions.

Zinc oxide nanostructures are used in an ever increasing line of applications in technology and biomedical fields. This requires a detailed understanding of the phenomena that occur at the surface particularly in aqueous environments and in contact with biomolecules. In this work, we used ab initio molecular dynamics (AIMD) simulations to determine structural details of ZnO surfaces in water and to develop a general and transferable classical force field for hydrated ZnO surfaces. AIMD simulations show that water molecules dissociate near unmodified ZnO surfaces, forming hydroxyl groups at about 65% of the surface Zn atoms and protonating 3-coordinated surface oxygen atoms, while the rest of the surface Zn atoms bind molecularly adsorbed waters. Several force field atom types for ZnO surface atoms were identified by analysis of the specific connectivities of atoms. The analysis of the electron density was then used to determine partial charges and Lennard-Jones parameters for the identified force field atom types. The obtained force field was validated by comparison with AIMD results and with available experimental data on adsorption and immersion enthalpies, as well as adsorption free energies of several amino acids in methanol. The developed force field can be used for modeling of ZnO in aqueous and other fluid environments and in interaction with biomolecules.

Adsorption profiles and adsorption free energies were determined for the side chain analogues of the 20 naturally occurring amino acids and a titanium binding peptide on the TiO2(100) surface. Microsecond simulations with umbrella sampling and metadynamics were used to sample the free energy barriers associated with desolvation of strongly bound water molecules at the TiO2surface. Polar and aromatic side chain analogues that hydrogen bond either to surface waters or directly to the metal oxide surface were found to be the strongest binders. Further, adsorption simulations of a 6-residue titanium binding peptide identified two binding modes on TiO2(100). The peptide structure with lowest free energy was shown to be stabilized by a salt bridge between the end termini. A comparison between the free energies of the side chain analogues of the peptide sequence and the peptide itself shows that the free energy contributions are not additive. The simulations emphasize that tightly bound surface waters play a key role for peptide and protein structures when bound to inorganic surfaces in biological environments.

Atomistic force field parameters were developed for the TiO2–water interface by systematic optimization with respect to experimentally determined crystal structures (lattice parameters) and surface thermodynamics (water adsorption enthalpy). Optimized force field parameters were determined for the two cases where TiO2was modeled with or without covalent bonding. The nonbonded TiO2model can be used to simulate different TiO2phases, while the bonded TiO2model is particularly useful for simulations of nanosized TiO2and biomatter, including protein–surface and nanoparticle–biomembrane simulations. The procedure is easily generalized to parametrize interactions between other inorganic surfaces and biomolecules.

Liposomes are proposed as drug delivery systems and can in principle be designed so as to cohere with specific tissue types or local environments. However, little detail is known about the exact mechanisms for drug delivery and the distributions of drug molecules inside the lipid carrier. In the current work, a coarse-grained (CG) liposome model is developed, consisting of over 2500 lipids, with varying degrees of drug loading. For the drug molecule, we chose hypericin, a natural compound proposed for use in photodynamic therapy, for which a CG model was derived and benchmarked against corresponding atomistic membrane bilayer model simulations. Liposomes with 21–84 hypericin molecules were generated and subjected to 10 microsecond simulations. Distribution of the hypericins, their orientations within the lipid bilayer, and the potential of mean force for transferring a hypericin molecule from the interior aqueous “droplet” through the liposome bilayer are reported herein.

Piezoelectric semiconductors can be polarized and used in mechanoredox systems and in photoredox catalysis. Conventional non-piezoelectric semiconductors have limitations when it comes to charge carrier recombination and slow transport rates in catalytic reactions, which can be overcome by piezoelectric polarization processes in piezoelectric semiconductors. Heterostructures based on semiconducting piezoelectrics often offer enhanced catalytic reactivities that are related to their mechanical, piezoelectric, optical, and electronic characteristics. We review how to use such heterostructures to convert mechanical energy into chemical energy, and how the related piezoelectric polarization tunes the band structures and provides advantages in piezophotocatalysis over regular photocatalysis. We discuss fundamental concepts of piezoelectricity, piezoelectric potential, and examine different piezoelectric heterostructures for piezo- and piezophotocatalysis. A review of dynamic investigations of piezo- and piezophotocatalytic processes is presented. The complementary developments in the understanding of the piezotronic and piezophototronic effects are described, which include the induced charge-transfer mechanisms for piezo- and piezophotocatalytic reactions that can occur with piezoelectric heterostructures. Finally, we derive design principles and suggest future research directions in the emerging field of piezo- and piezophotocatalysis employing semiconductive heterostructures.

We present software package MagiC, which is designed to perform systematic structure-based coarse graining of molecular models. The effective pairwise potentials between coarse-grained sites of low-resolution molecular models are constructed to reproduce structural distribution functions obtained from the modeling of the system in a high resolution (atomistic) description. The software supports coarse-grained tabulated intramolecular bond and angle interactions, as well as tabulated nonbonded interactions between different site types in the coarse-grained system, with the treatment of long-range electrostatic forces by the Ewald summation. Two methods of effective potential refinement are implemented: iterative Boltzmann inversion and inverse Monte Carlo, the latter accounting for cross-correlations between pair interactions. MagiC uses its own Metropolis Monte Carlo sampling engine, allowing parallel simulation of many copies of the system with subsequent averaging of the properties, which provides fast convergence of the method with nearly linear scaling at parallel execution.

To be able to model complex biological membranes in a more realistic manner, the force field Slipids (Stockholm lipids) has been extended to include parameters for sphingomyelin (SM), phosphatidylglycerol (PG), phosphatidylserine (PS) lipids, and cholesterol. Since the parametrization scheme was faithful to the scheme used in previous editions of Slipids, all parameters are consistent and fully compatible. The results of careful validation of a number of key structural properties for one and two component lipid bilayers are in excellent agreement with experiments. Potentials of mean force for transferring water across binary mixtures of lipids and cholesterol were also computed in order to compare water permeability rates to experiments. In agreement with experimental and simulation studies, it was found that the permeability and partitioning of water is affected by cholesterol in lipid bilayers made of saturated lipids to the largest extent. With the extensions of Slipids presented here, it is now possible to study complex systems containing many different lipids and proteins in a fully atomistic resolution in the isothermic–isobaric (NPT) ensemble, which is the proper ensemble for membrane simulations.

Biological membranes are versatile in composition and host intriguing molecular processes. In order to be able to study these systems, an accurate model Hamiltonian or force field (FF) is a necessity. Here, we report the results of our extension of earlier developed all-atomistic FF parameters for fully saturated phospholipids that complements an earlier parameter set for saturated phosphatidylcholine lipids (J. Phys. Chem. B, 2012, 116, 3164–3179). The FF, coined Slipids (Stockholm lipids), now also includes parameters for unsaturated phosphatidylcholine and phosphatidylethanolamine lipids, e.g., POPC, DOPC, SOPC, POPE, and DOPE. As the extended set of parameters is derived with the same philosophy as previously applied, the resulting FF has been developed in a fully consistent manner. The capabilities of Slipids are demonstrated by performing long simulations without applying any surface tension and using the correct isothermal–isobaric (NPT) ensemble for a range of temperatures and carefully comparing a number of properties with experimental findings. Results show that several structural properties are very well reproduced, such as scattering form factors, NMR order parameters, thicknesses, and area per lipid. Thermal dependencies of different thicknesses and area per lipid are reproduced as well. Lipid diffusion is systematically slightly underestimated, whereas the normalized lipid diffusion follows the experimental trends. This is believed to be due to the lack of collective movement in the relatively small bilayer patches used. Furthermore, the compatibility with amino acid FFs from the AMBER family is tested in explicit transmembrane complexes of the WALP23 peptide with DLPC and DOPC bilayers, and this shows that Slipids can be used to study more complex and biologically relevant systems.

The physical mechanism of the folding and unfolding of chromatin is fundamentally related to transcription but is incompletely characterized and not fully understood. We experimentally and theoretically studied chromatin compaction by investigating the salt-mediated folding of an array made of 12 positioning nucleosomes with 177 bp repeat length. Sedimentation velocity measurements were performed to monitor the folding provoked by addition of cations Na+, K+, Mg2+, Ca2+, spermidine3+, Co(NH3)63+, and spermine4+. We found typical polyelectrolyte behavior, with the critical concentration of cation needed to bring about maximal folding covering a range of almost five orders of magnitude (from 2 μM for spermine4+to 100 mM for Na+). A coarse-grained model of the nucleosome array based on a continuum dielectric description and including the explicit presence of mobile ions and charged flexible histone tails was used in computer simulations to investigate the cation-mediated compaction. The results of the simulations with explicit ions are in general agreement with the experimental data, whereas simple Debye-Hückel models are intrinsically incapable of describing chromatin array folding by multivalent cations. We conclude that the theoretical description of the salt-induced chromatin folding must incorporate explicit mobile ions that include ion correlation and ion competition effects.

Ionic liquids (ILs) are a class of new green materials that have attracted extensive attention in recent decades. Many novel properties not evident under normal conditions may appear when ionic liquids are confined to a nanometer scale. As was observed in the experiment, an anomalous phase behavior from liquid to high melting point perfect crystal occurred when 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6]) ionic liquid was confined in a carbon nanotube. In this work, we performed molecular dynamics (MD) simulations for [bmim][PF6] ionic liquid and provided direct structural evidence that the ionic crystallizes in a carbon nanotube. The ordered ionic arrangement in both the radial and the axial directions can be observed inside the channels of the CNTs to induce the form of crystallites. The ionic stacking and distributing can be determined by the sizes of the CNTs. Hydrogen bonds remain the dominant interactions between cations and anions when the ionic liquid enters into the CNT from the bulk phase. The free energies as the thermal driven forces were calculated, and it is found that it is very difficult for a single anion to enter into the channel of the CNT spontaneously. A more favorable way is through an ion-pair in which a cation “pulls” an anion to enter into the channel of the CNT together. It is predicted that other ionic liquids that possess similar structures, even including the pyridinium-based ionic liquids, can show higher melting points when confined in CNTs.

Coarse-grained Langevin molecular dynamics computer simulations were conducted for systems that mimic solutions of nucleosome core particles (NCPs). The NCP was modeled as a negatively charged spherical particle representing the complex of DNA and the globular part of the histones combined with attached strings of connected charged beads modeling the histone tails. The size, charge, and distribution of the tails relative to the core were built to match real NCPs. Three models of NCPs were constructed to represent different extents of covalent modification on the histone tails: (nonmodified) recombinant (rNCP), acetylated (aNCP), and acetylated and phosphorylated (paNCP). The simulation cell contained 10 NCPs in a dielectric continuum with explicit mobile counterions and added salt. The NCP-NCP interaction is decisively dependent on the modification state of the histone tails and on salt conditions. Increasing the monovalent salt concentration (KCl) from salt-free to physiological concentration leads to NCP aggregation in solution for rNCP, whereas NCP associates are observed only occasionally in the system of aNCPs. In the presence of divalent salt (Mg2+), rNCPs form dense stable aggregates, whereas aNCPs form aggregates less frequently. Aggregates are formed via histone-tail bridging and accumulation of counterions in the regions of NCP-NCP contacts. The paNCPs do not show NCP-NCP interaction upon addition of KCl or in the presence of Mg2+. Simulations for systems with a gradual substitution of K+for Mg2+, to mimic the Mg2+titration of an NCP solution, were performed. The rNCP system showed stronger aggregation that occurred at lower concentrations of added Mg2+, compared to the aNCP system. Additional molecular dynamics simulations performed with a single NCP in the simulation cell showed that detachment of the tails from the NCP core was modest under a wide range of salt concentrations. This implies that salt-induced tail dissociation of the histone tails from the globular NCP is not in itself a major factor in NCP-NCP aggregation. The approximation of coarse-graining, with respect to the description of the NCP as a sphere with uniform charge distribution, was tested in control simulations. A more detailed description of the NCP did not change the main features of the results. Overall, the results of this work are in agreement with experimental data reported for NCP solutions and for chromatin arrays.

The influence of the local anesthetic lidocaine on electrostatic properties of a lipid membrane bilayer was studied by molecular dynamics simulations. The electrostatic dipole potential, charge densities, and orientations of the headgroup angle have been examined in the presence of different amounts of charged or uncharged forms of lidocaine. Important changes in the membrane properties caused by the presence of both forms of lidocaine are presented and discussed. Our simulations have shown that both charged and uncharged lidocaine cause almost the same increase in the electrostatic potential in the middle of the membrane, although for different reasons. The increase, ∼90mV for 9mol % of lidocaine and 220mV for 28mol % of lidocaine, is of a size that may affect the functioning of voltage-gated ion channels.

We conducted molecular dynamics computer simulations of charged histone tail-DNA interactions in systems mimicking nucleosome core particles (NCP). In a coarse-grained model, the NCP is modeled as a negatively charged spherical particle with flexible polycationic histone tails attached to it in a dielectric continuum with explicit mobile counterions and added salt. The size, charge, and distribution of the tails relative to the core were built to mimick real NCP. In this way, we incorporate attractive ion-ion correlation effects due to fluctuations in the ion cloud and the attractive entropic and energetic tail-bridging effects. In agreement with experimental data, increase of monovalent salt content from salt-free to physiological concentration leads to the formation of NCP aggregates; likewise, in the presence of MgCl2, the NCPs form condensed systems via histone-tail bridging and accumulation of counterions. More detailed mechanisms of the histone tail-DNA interactions and dynamics have been obtained from all-atom molecular dynamics simulations (including water), comprising three DNA 22-mers and 14 short fragments of the H4 histone tail (amino acids 5–12) carrying three positive charges on lysine+ interacting with DNA. We found correlation of the DNA-DNA distance with the presence and association of the histone tail between the DNA molecules.

Summary:We extend Monte Carlo (MC) methods developed in our previous paper (J. Phys. A, Math. Gen.2004, 37, 1573) and based on entropic sampling within Wang‐Landau (WL) algorithm to the simulation of lattice and continuous models of ring polymers. For a continuous freely joined ring chain (an equilateral polygon) with hard sphere monomer units, the excess entropy of rings relative to the corresponding reference system, a phantom ring chain, is obtained. The excess specific entropy is calculated for a set of various diameters of monomer units dand chain lengths N. Its limiting values for $N\to \infty$($N^{-1}\to 0$) are estimated for each dand coincidence with those for corresponding free chains is demonstrated. We also develop a WL approach to calculate thermal properties of free and ring continuous chains with an Lennard‐Jones attraction between nonbonded beads being added to hard core repulsion at fixed d. The obtained energy distributions provide calculation of canonical properties such as conformational energy, heat capacity, entropy, and mean square radius of inertia. Thermal results for free and ring chains are being finally compared. Analogous calculations are performed for lattice‐free chains and rings.

We report a detailed comparison between calculations of inter-filament interactions based on Monte-Carlo simulations and experimental features of lateral aggregation of bacteriophages fd and M13 induced by a number of divalent metal ions. The general findings are consistent with the polyelectrolyte nature of the virus filaments and confirm that the solution electrostatics account for most of the experimental features observed. One particularly interesting discovery is resolubilization for bundles of either fd or M13 viruses when the concentration of the bundle-inducing metal ion Mg2+or Ca2+is increased to large (>100mM) values. In the range of Mg2+or Ca2+concentrations where large bundles of the virus filaments are formed, the optimal attractive interaction energy between the virus filaments is estimated to be on the order of 0.01 kT per net charge on the virus surface when a recent analytical prediction to the experimentally defined conditions of resolubilization is applied. We also observed qualitatively distinct behavior between the alkali-earth metal ions and the divalent transition metal ions in their action on the charged viruses. The understanding of metal ions-induced reversible aggregation based on solution electrostatics may lead to potential applications in molecular biology and medicine.

The interaction of DNA with the polyamine spermine4+ (Spm4+), sodium ions, and water molecules has been studied using molecular dynamics computer simulations in a system modeling a DNA crystal. The simulation model consisted of three B-DNA decamers in a periodic hexagonal cell, containing 1200 water molecules, 8 Spm4+, 32 Na+, and 4 Cl− ions. The present paper gives a more detailed account of a recently published report of this system and compares results on this mixed Spm4+/Na+-cation system with an molecular dynamics simulation carried out for the same DNA decamer under similar conditions with only sodium counterions (Korolev et al., 2001, J. Mol. Biol. 308:907). The presence of Spm4+ makes significant influence on the DNA hydration and on the interaction of the sodium ions with DNA. Spermine pushes water molecules out of the minor groove, whereas Na+ attracts and organizes water around DNA. The major binding site of the Spm4+ amino groups and the Na+ ions is the phosphate group of DNA. The flexible polyamine spermine displays a high presence in the minor groove but does not form long-lived and structurally defined complexes. Sodium ions compete with Spm4+ for binding to the DNA bases in the minor groove. Sodium ions also have several strong binding sites in the major groove. The ability of water molecules, Spm4+, and Na+ to modulate the local structure of the DNA double helix is discussed.

The association of spermine4+ (Spm4+), Mg2+ and monovalent (M+) ions with DNA in crystal form, have been studied using grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) computer simulations. GCMC calculations were used to calculate the distribution of Spm4+, Mg2+, and M+ between the equilibrating solvent and the DNA crystal under conditions mimicking the crystal-growing protocols reported in a number of recent X-ray diffraction studies of DNA oligomers. The GCMC simulations show that the composition of ions neutralizing the negative charge of DNA can vary in a broad range. The GCMC simulations were used to provide appropriate conditions for subsequent 6 ns constant pressure and temperature MD simulations of DNA in a typical crystalline environment consisting of three DNA double helix decamers in a periodic hexagonal cell, containing 1200 water molecules, eight Spm4+, 32 Na+ and four Cl− ions. Based on the simulation results, it seems possible to give an explanation why spermine molecules are usually not detected in X-ray studies in spite of their high concentration in the preparatory samples used as the crystallizing agent. It appears that this flexible polyamine molecule has several binding modes, interacting in fairly irregular manner with different sites on DNA and showing no regular ordering in the DNA crystals. Ions of Na+ and Spm4+ compete with each other and with water molecules in binding to bases in the minor groove and they influence the structure of the DNA hydration shell in different ways.

Competition of the trivalent cation, Co(NH3)³⁺6, with K⁺ and Na⁺ ions in binding to DNA was studied by equilibrating oriented DNA fibers with ethanol/water solutions (65 and 52% v/v EtOH), containing different combinations and concentrations of KCl and NaCl and constant concentration (0.8 mM) of Co(NH3)6Cl3. The degree of Co(NH3)³⁺6 binding to DNA does not depend significantly on the ethanol concentration or on the kind of univalent cation (Na⁺ or K⁺). The ion exchange selectivity coefficient of monovalent–trivalent ion competition, D¹c3, increases with the concentration of Me⁺, C^o+, and the monotonic dependence of log D¹c3 vs log C^o+ has an inflection between 100 and 300 mM that is caused by a structural transformation of DNA from A- to B-form. The ion exchange experimental data are compared with results of grand canonical Monte Carlo (GCMC) simulations of systems of parallel and hexagonally ordered, discretely charged polyions with density and spatial distribution of the charged groups modeling B- and A-forms of DNA. The GCMC method for discretely charged models of the DNA polyion produces a quantitative agreement with experimental data on trivalent–monovalent ion competition in dependence on DNA structural state and salt concentration. Based on this and previous studies it is concluded that the affinity of DNA for the cations decreases in the order Co(NH3)³⁺6 ≫ Ca²⁺ > Mg²⁺ ≫ Na⁺ ≈ K⁺ > Li⁺. DNA does not exhibit selectivity for Na⁺ or K⁺ in ethanol/water solutions either in the absence or in the presence of Co(NH3)³⁺6, Ca²⁺, and Mg²⁺. © 2001 John Wiley & Sons, Inc. Biopolymers 58: 268–278, 2001

Competitive binding of the most common cations of the cytoplasm (K+, Na+, Ca2+, and Mg2+) with DNA was studied by equilibrating oriented DNA fibers with ethanol/water solutions (65 and 52% v/v EtOH) containing different combinations and concentrations of the counterions. The affinity of DNA for the cations decreases in the order Ca>Mg≫Na≈K. The degree of Ca2+ and/or Mg2+ binding to DNA displays maximum changes just at physiological concentrations of salts (60–200mM) and does not depend significantly on the ethanol concentration or on the kind of univalent cation (Na+ or K+). Ca2+ is more tightly bound to DNA and is replaced by the monovalent cations to a lesser extent than is Mg2+. Similarly, Ca2+ is a better competitor for binding to DNA than Mg2+: the ion exchange equilibrium constant for a 1:1 mixture of Ca2+ and Mg2+ ions, KcCaMg, changes from KcCaMg≈2 in 65% EtOH (in 3–30mM NaCl and/or KCl) to KcCaMg≈1.2–1.4 in 52% EtOH (in 300mM NaCl and/or KCl). DNA does not exhibit selectivity for Na+ or K+ in ethanol/water solutions either in the absence or in the presence of Ca2+ and/or Mg2+. The ion exchange experimental data are compared with results of grand canonical Monte Carlo (GCMC) simulations of systems of parallel and hexagonally ordered, uniformly and discretely charged polyions with the density and spatial distribution of the charged groups modeling B DNA. A quantitative agreement with experimental data on divalent-monovalent competition has been obtained for discretely charged models of the DNA polyion (for the uniformly charged cylinder model, coincidence with experiment is qualitative). The GCMC method gives also a qualitative description of experimental results for DNA binding competitions of counterions of the same charge (Ca2+ with Mg2+ or K+ with Na+).

Numerical calculations, using Poisson-Boltzmann (PB) and counterion condensation (CC) polyelectrolyte theories, of the electrostatic free energy difference, ΔGel, between single-stranded (coil) and double-helical DNA have been performed for solutions of NaDNA+NaCl with and without added MgCl2. Calculations have been made for conditions relevant to systems where experimental values of helix coil transition temperature (Tm) and other thermodynamic quantities have been measured. Comparison with experimental data has been possible by invoking values of Tm for solutions containing NaCl salt only. Resulting theoretical values of enthalpy, entropy, and heat capacity (for NaCl salt-containing solutions) and of Tm as a function of NaCl concentration in NaCl+MgCl2 solutions have thus been obtained. Qualitative and, to a large extent, quantitative reproduction of the experimental Tm, ΔHm, ΔSm, and ΔCp values have been found from the results of polyelectrolyte theories. However, the quantitative resemblance of experimental data is considerably better for PB theory as compared to the CC model. Furthermore, some rather implausible qualitative conclusions are obtained within the CC results for DNA melting in NaCl+MgCl2 solutions. Our results argue in favor of the Poisson-Boltzmann theory, as compared to the counterion condensation theory.

A method to reconstruct the pair interaction potential from radial distribution functions within the Monte Carlo simulation scheme is described and some applications are briefly discussed.