Monte Carlo simulations of the hydrophobic effect in aqueous electrolyte solutions

The hydrophobic interaction between two methane molecules in salt-free and high salt-containing aqueous solutions and the structure in such solutions have been investigated using an atomistic model solved by Monte Carlo simulations. Monovalent salt representing NaCl and divalent salt with the same nonelectrostatic properties as the monovalent salt have been used to examine the influence of the valence of the salt species. In salt-free solution the effective interaction between the two methane molecules displayed a global minimum at close contact of the two methane molecules and a solvent-separated secondary minimum. In 3 and 5 M monovalent salt solution the potential of mean force became slightly more attractive, and in a 3 M divalent salt solution the attraction became considerably stronger. The structure of the aqueous solutions was determined by radial distribution functions and angular probability functions. The distortion of the native water structure increased with ion valence. The increase of the hydrophobic attraction was associated with (i) a breakdown of the tetrahedral structure formed by neighboring water molecules and of the hydrogen bonds between them and (i) the concomitant increase of the solution density.     

Stretching,Tearing, and Dissecting Single Molecules of DNA

Optical tweezers, bendable microneedles, and scanning force microscope probes make it possible to play with individual molecules of DNA, to stretch them beyond their natural length, to unzip and pull apart their strands (see schematic diagram), and to dissect them to create new molecules in situ. Depending on the method of measurement, the mechanical force necessary to separate the strands was in the range of 10–50 pN per base pair.     

"Like-charge attraction" between anionic polyelectrolytes: molecular dynamics simulations

"Like-charge attraction" is a phenomenon found in many biological systems containing DNA or proteins, as well as in polyelectrolyte systems of industrial importance. "Like-charge attraction" between polyanions is observed in the presence of mobile multivalent cations. At a certain limiting concentration of cations, the negatively charged macroions cease to repel each other and even an attractive force between the anions is found. With classical molecular dynamics simulations it is possible to elucidate the processes that govern the attractive behavior with atomistic resolution. As an industrially relevant example we study the interaction of negatively charged carboxylate groups of sodium polyacrylate molecules with divalent cationic Ca2+ counterions. Here we show that Ca2+ ions initially associate with single chains of polyacrylates and strongly influence sodium ion distribution; shielded polyanions approach each other and eventually "stick" together (precipitate), contrary to the assumption that precipitation is initially induced by intermolecular Ca2+ bridging.     

DNA Structural Changes Under Different Stretching Methods Studied By Molecular Dynamics Simulations

We present a molecular dynamics simulation study of 22‐mer DNA conformational variations obtained by stretching both 3′‐termini and both 5′‐termini. Stretching 3′‐termini by 3.5 nm required 142 kJ mol^?1 and the force plateau was ～80 pN, whereas stretching 5′‐termini by the same length required 190 kJ mol^?1 and the force plateau was ～100 pN. Stretching 3′‐termini led to a larger untwisting of the double helix and the successive base pairs rolled to the side of the DNA minor groove, while stretching 5′‐termini resulted in the base pairs rolling to the major groove side and reducing of the diameter of DNA molecule. The most distinctive difference between stretching 3′‐termini and 5′‐termini was that at the force plateau region stretching the 5′‐termini resulted in breakage of the base pairs, which considerably disturbed the structure of the DNA double helix. All of the variations of base rotation and translation for both stretching methods took place when the relative length of DNA l was longer than 1.2, which was the point the force plateau appeared.     

Effects of electrostatic screening on the conformation of single DNA molecules confined in a nanochannel

Single T4-DNA molecules were confined in rectangular-shaped channels with a depth of 300 nm and a width in the range of 150-300 nm casted in a poly(dimethylsiloxane) nanofluidic chip. The extensions of the DNA molecules were measured with fluorescence microscopy as a function of the ionic strength and composition of the buffer as well as the DNA intercalation level by the YOYO-1 dye. The data were interpreted with the scaling theory for a wormlike polymer in good solvent, including the effects of confinement, charge, and self-avoidance. It was found that the elongation of the DNA molecules with decreasing ionic strength can be interpreted in terms of an increase of the persistence length. Self-avoidance effects on the extension are moderate, due to the small correlation length imposed by the channel cross-sectional diameter. Intercalation of the dye results in an increase of the DNA contour length and a partial neutralization of the DNA charge, but besides effects of electrostatic origin it has no significant effect on the bare bending rigidity. In the presence of divalent cations, the DNA molecules were observed to contract, but they do not collapse into a condensed structure. It is proposed that this contraction results from a divalent counterion mediated attractive force between the segments of the DNA molecule.     

Direct detection of the formation of V-amylose helix by single molecule force spectroscopy

An important polysaccharide, amylose crystallizes as a regular single left-handed helix from a propanol, butanol, or iodine solution. However, its solution structure remains elusive because amylose does not form molecular solutions in these solvents, and standard spectroscopic techniques cannot be exploited to determine its structure. Using AFM, we forced individual amylose chains adsorbed to a surface to enter these poor solvents and carried out stretch-release measurements on them in solution. In this manner, we directly captured the formation of individual amylose helices induced by butanol and iodine. With an accuracy approaching that of X-ray diffraction on amylose crystals, we determined that the pitch of the helix in solution is 1.3 angstroms/ring. We also directly measured the force driving the formation of the helix in solution to be 50 pN. SMD simulations in explicit butanol reproduced the AFM-measured force-extension curves and revealed that the long plateau feature is caused by the rupture of O(2)n-O(6)(n+6) and O(3)n-O(6)(n+6) hydrogen bonds and by the unwinding of the helix. We also found that amylose helices formed in iodine solution are more compliant and hysteretic as compared to helices in butanol, which extend/relax reversibly. In iodine solution, the formation of the helix is inhibited by force and limited by the slow kinetics of the amylose-iodine complex. By forcing individual molecules into poor solvents and performing force spectroscopy measurements in solution, our AFM approach uniquely supplements X-ray diffraction and NMR methods for investigating solution conformations of insoluble biopolymers.     

Probing multivalent interactions in a synthetic host-guest complex by dynamic force spectroscopy

Multivalency is present in many biological and synthetic systems. Successful application of multivalency depends on a correct understanding of the thermodynamics and kinetics of this phenomenon. In this Article, we address the stability and strength of multivalent bonds with force spectroscopy techniques employing a synthetic adamantane/β-cyclodextrin model system. Comparing the experimental findings to theoretical predictions for the rupture force and the kinetic off-rate, we find that when the valency of the complex is increased from mono- to di- to trivalent, there is a transition from quasi-equilibrium, with a constant rupture force of 99 pN, to a kinetically dependent state, with loading-rate-dependent rupture forces from 140 to 184 pN (divalent) and 175 to 210 pN (trivalent). Additional binding geometries, parallel monovalent ruptures, single-bound divalent ruptures, and single- and double-bound trivalent ruptures are identified. The experimental kinetic off-rates of the multivalent complexes show that the stability of the complexes is significantly enhanced with the number of bonds, in agreement with the predictions of a noncooperative multivalent model.     

Single‐Molecule Force Spectroscopy on Carrageenan by Means of AFM

A marked difference in force‐extension curves is observed for carrageenan before and after adding NaI buffer in single‐molecule force spectroscopy by means of atomic force microscopy (AFM). The salt‐induced helix conformation in carrageenan treated with an 0.1 M NaI solution was unfolded under the external force, and a long plateau about 300 pN high could be observed in the force‐extension curves.     

ONIOM approach for non-adiabatic on-the-fly molecular dynamics demonstrated for the backbone controlled Dewar valence isomerization

Non-adiabatic on-the-fly molecular dynamics (NA-O-MD) simulations require the electronic wavefunction, energy gradients, and derivative coupling vectors in every timestep. Thus, they are commonly restricted to the excited state dynamics of molecules with up to ≈20 atoms. We discuss an approximation that combines the ONIOM(QM:QM) method with NA-O-MD simulations to allow calculations for larger molecules. As a proof of principle we present the excited state dynamics of a (6-4)-lesion containing dinucleotide (63 atoms), and especially the importance to include the confinement effects of the DNA backbone. The method is able to include electron correlation on a high level of theory and offers an attractive alternative to QM:MM approaches for moderate sized systems with unknown force fields.     

Electrolyte distribution around two like-charged rods: their effective attractive interaction and angular dependent charge reversal

A simple model for two like-charged parallel rods immersed in an electrolyte solution is considered. We derived the three point extension (TPE) of the hypernetted chain/mean spherical approximation (TPE-HNC/MSA) and Poisson-Boltzmann (TPE-PB) integral equations. We numerically solve these equations and compare them to our results of Monte Carlo (MC) simulations. The effective interaction force, F(T), the charge distribution profiles, rho(el)(x,y), and the angular dependent integrated charge function, P(theta), are calculated for this system. The analysis of F(T) is carried out in terms of the electrostatic and entropic (depletion) contributions, F(E) and F(C). We studied several cases of monovalent and divalent electrolytes, for which the ionic size and concentration are varied. We find good qualitative agreement between TPE-HNC/MSA and MC in all the cases studied. The rod-rod force is found to be attractive when immersed in large size, monovalent or divalent electrolytes. In general, the TPE-PB has poor agreement with the MC. For large monovalent and divalent electrolytes, we find angular dependent charge reversal charge inversion and polarizability. We discuss the intimate relationship between this angular dependent charge reversal and rod-rod attraction.     

Effects of contact force and salt concentration on the unbinding of a DNA duplex by force spectroscopy

We report that varying the contact force in force spectroscopy results in a significant shift in DNA unbinding forces, measured from short oligonucleotides using a PicoForce microscope. The contact force between a 30-mer complementary DNA-coated probe and surface was varied from 100 pN to 10 nN, resulting in a significant shift in the most abundant unbinding force measured between the duplex. When contact forces were set at 200 pN or less, which is generally considered to be a low contact force region for biomolecular force spectroscopy studies, the shift in DNA unbinding forces was significant with changes in contact force. The effect of the salt concentration on the DNA unbinding forces was also examined for a range of salt concentrations from 5 to 500 mM because the presence of salt ions is necessary to facilitate the hybridization process. Although an increase in salt concentration resulted in the facilitation of DNA multiple binding events during force spectroscopy measurements, no significant shift in unbinding forces was observed. Our experiment demonstrates that the wide variation in DNA unbinding forces reported in the literature (50-600 pN) for short oligonucleotides can be accounted for by the different contact forces used and shows little or no effect of the salt concentration used in those studies. Furthermore, this study demonstrates the importance of reporting contact forces in force spectroscopy measurements for quantitative comparisons between different biomolecular systems, especially for noncovalent-type interactions.     

How does huperzine A enter and leave the binding gorge of acetylcholinesterase? Steered molecular dynamics simulations

The entering and leaving processes of Huperzine A (HupA) binding with the long active-site gorge of Torpedo californica acetylcholinesterase (TcAChE) have been investigated by using steered molecular dynamics simulations. The analysis of the force required along the pathway shows that it is easier for HupA to bind to the active site of AChE than to disassociate from it, which for the first time interprets at the atomic level the previous experimental result that unbinding process of HupA is much slower than its binding process to AChE. The direct hydrogen bonds, water bridges, and hydrophobic interactions were analyzed during two steered molecular dynamics (SMD) simulations. Break of the direct hydrogen bond needs a great pulling force. The steric hindrance of bottleneck might be the most important factor to produce the maximal rupture force for HupA to leave the binding site but it has a little effect on the binding process of HupA with AChE. Residue Asp72 forms a lot of water bridges with HupA leaving and entering the AChE binding gorge, acting as a clamp to take out HupA from or put HupA into the active site. The flip of the peptide bond between Gly117 and Gly118 has been detected during both the conventional MD and SMD simulations. The simulation results indicate that this flip phenomenon could be an intrinsic property of AChE and the Gly117-Gly118 peptide bond in both HupA bound and unbound AChE structures tends to adopt the native enzyme structure. At last, in a vacuum the rupture force is increased up to 1500 pN while in water solution the greatest rupture force is about 800 pN, which means water molecules in the binding gorge act as lubricant to facilitate HupA entering or leaving the binding gorge.     

Two-dimensional LNA/DNA arrays: estimating the helicity of LNA/DNA hybrid duplex

We measured the helical repeats of a non-natural nucleic acid, locked nucleic acid (LNA), by incorporating LNA strands into the outer arms of a DNA double crossover (DX) molecule; atomic force microscopy (AFM) imaging of the two-dimensional (2D) arrays self-assembled from these DX molecules allows us to derive the helical repeat of the LNA/DNA hetero-duplex to be 13.2 +/- 0.9 base pairs per turn.     

Theoretical study of intramolecular interaction energies during dynamics simulations of oligopeptides by the fragment molecular orbital-Hamiltonian algorithm method

We analyzed the interaction energies between residues (fragments) in an oligopeptide occurring during dynamic simulations by using the fragment molecular orbital-Hamiltonian algorithm (FMO-HA) method, an ab initio MO-molecular dynamics technique. The FMO method enables not only calculation of large molecules based on ab initio MO but also easy evaluation of interfragment interaction energies. The glycine pentamer [(Gly)(5)] and decamer [(Gly)(10)] were divided into five and ten fragments, respectively. alpha-helix structures of (Gly)(5) and (Gly)(10) were stabilized by attractive interaction energies owing to intramolecular hydrogen bonds between fragments n and n+3 (and n+4), and beta-strand structures were characterized by repulsive interaction energies between fragments n and n+2. We analyzed interfragment interaction energies during dynamics simulations as the peptides' geometries changed from alpha helix to beta strand. Intramolecular hydrogen bonds between fragments 2-4 and 2-5 control the geometrical preference of (Gly)(5) for the beta-strand structure. The pitch of one turn of the alpha helix of (Gly)(10) elongated and thus weakened during dynamics due to a shifting of the intramolecular hydrogen bonds, and enabled the beta-strand structure to form. Changes in interaction energies due to the intramolecular hydrogen bonds controlled the tendency toward alpha-helix or beta-strand structure of (Gly)(5) and (Gly)(10). Evaluation of interfragment interaction energies during dynamics simulations thus enabled detailed analysis of the process of the geometrical changes occurring in oligopeptides.     

Effects of degassing and ionic strength on AFM force measurements in octadecyltrimethylammonium chloride solutions

Sakamoto et al. (Langmuir 2002, 18, 5713) conducted AFM force measurements between silica sphere and fused-silica plate in aqueous octadecyltrimethylammonium chloride (C18TACl) solutions and concluded that long-range attractive force is not observed in carefully degassed solutions. In the present work, AFM force measurements were conducted by following the procedures described by Sakamoto et al. The results showed the presence of an attractive force that was much stronger than the van der Waals force both in air-saturated and degassed solutions. The force was most attractive at 5 x 10(-6) M C18TACl, where contact angle was maximum. At this concentration, which is close to the charge compensation point (ccp) of the glass sphere, the long-range decay lengths (D) were 34 and 38 nm in air-saturated and degassed solutions, respectively. At 10(-5) M, the decay length decreased from 30 to 4 nm upon degassing. This decrease in decay length can be explained by a pH increase (from 5.7 to 6.6), which in turn causes additional surfactant molecules to adsorb on the surface with inverse orientation. The attractive force was screened by an added electrolyte (NaCl), indicating that the attractive force may be of electrostatic origin. Therefore, the very long decay lengths observed in the absence of electrolyte may be ascribed to the fact that the ccp occurs at a very low surfactant concentration.     

Toward a designed, functioning genetic system with expanded-size base pairs: solution structure of the eight-base xDNA double helix

We describe the NMR-derived solution structure of the double-helical form of a designed eight-base genetic pairing system, termed xDNA. The benzo-homologous xDNA design contains base pairs that are wider than natural DNA pairs by ca. 2.4 A (the width of a benzene ring). The eight component bases of this xDNA helix are A, C, G, T, xA, xT, xC, and xG. The structure was solved in aqueous buffer using 1D and 2D NMR methods combined with restrained molecular dynamics. The data show that the decamer duplex is right-handed and antiparallel, and hydrogen-bonded in a way analogous to that of Watson-Crick DNA. The sugar-phosphate backbone adopts a regular conformation similar to that of B-form DNA, with small dihedral adjustments due to the larger circumference of the helix. The grooves are much wider and more shallow than those of B-form DNA, and the helix turn is slower, with ca. 12 base pairs per 360 degrees turn. There is an extensive intra- and interstrand base stacking surface area, providing an explanation for the greater stability of xDNA relative to natural DNA. There is also evidence for greater motion in this structure compared to a previous two-base-expanded helix; possible chemical and structural reasons for this are discussed. The results confirm paired self-assembly of the designed xDNA system. This suggests the possibility that other genetic system structures besides the natural one might be functional in encoding information and transferring it to new complementary strands.     

Interaction between dendrons directly studied by single-molecule force spectroscopy

In this article, we have investigated the interaction between two poly(benzyl ether) dendrons directly by single-molecule force spectroscopy. For this purpose, one dendron was immobilized on an AFM tip through a poly(ethylene glycol) (PEG) spacer, and the other dendron was anchored on a gold substrate as a self-assembled monolayer. Two dendrons approached and then interacted with each other when the AFM tip and the substrate moved close together. The rupture force between dendrons was measured while the AFM tip and the substrate separated. PEG as a flexible spacer can function as a length window for recognizing the force signals and avoiding the disturbance of the interaction between the AFM tip and the substrate. The interaction between two first-generation dendrons is measured to be about 224 pN at a force loading rate of 40 nN/s. The interaction between second- and first-generation dendrons rises to 315 pN at the same loading rate. Such interactions depend on the force loading rate in the range of several to hundreds of nanonewtons per second, indicating that the rupture between dendrons is a dynamic process. The study of the interaction between surface-bound dendrons of different generations provides a model system for understanding the surface adhesion of molecules with multiple branches. In addition, this multiple-branch molecule may be used to mimic the sticky feet of geckos as a man-made adhesive.     

Anionic polyelectrolyte adsorption on mica mediated by multivalent cations: a solution to DNA imaging by atomic force microscopy under high ionic strengths

Adsorption of DNA molecules on mica, a highly negatively charged surface, mediated by divalent or trivalent cations is considered. By analyzing atomic force microscope (AFM) images of DNA molecules adsorbed on mica, phase diagrams of DNA molecules interacting with a mica surface are established in terms of concentrations of monovalent salt (NaCl) and divalent (MgCl2) or multivalent (spermidine, cobalt hexamine) salts. These diagrams show two transitions between nonadsorption and adsorption. The first one arises when the concentration of multivalent counterions is larger than a limit value, which is not sensitive to the monovalent salt concentration. The second transition is due to the binding competition between monovalent and multivalent counterions. In addition, we develop a model of polyelectrolyte adsorption on like-charged surfaces with multivalent counterions. This model shows that the correlations of the multivalent counterions at the interface between DNA and mica play a critical role. Furthermore, it appears that DNA adsorption takes place when the energy gain in counterion correlations overcomes an energy barrier. This barrier is induced by the entropy loss in confining DNA in a thin adsorbed layer, the entropy loss in the interpenetration of the clouds of mica and DNA counterions, and the electrostatic repulsion between DNA and mica. The analysis of the experimental results provides an estimation of this energy barrier. We then discuss some important issues, including DNA adsorption under physiological conditions.     

Nonequilibrium molecular dynamics of the rheological and structural properties of linear and branched molecules. Simple shear and poiseuille flows; instabilities and slip

Nonequilibrium molecular-dynamics simulations are performed for linear and branched chain molecules to study their rheological and structural properties under simple shear and Poiseuille flows. Molecules are described by a spring-monomer model with a given intermolecular potential. The equations of motion are solved for shear and Poiseuille flows with Lees and Edward's [A. W. Lees and S. F. Edwards, J. Phys. C 5, 1921 (1972)] periodic boundary conditions. A multiple time-scale algorithm extended to nonequilibrium situations is used as the integration method, and the simulations are performed at constant temperature using Nose-Hoover [S. Nose, J. Chem. Phys. 81, 511 (1984)] dynamics. In simple shear, molecules with flow-induced ellipsoidal shape, having significant segment concentrations along the gradient and neutral directions, exhibit substantial flow resistance. Linear molecules have larger zero-shear-rate viscosity than that of branched molecules, however, this behavior reverses as the shear rate is increased. The relaxation time of the molecules is associated with segment concentrations directed along the gradient and neutral directions, and hence it depends on structure and molecular weight. The results of this study are in qualitative agreement with other simulation studies and with experimental data. The pressure (Poiseuille) flow is induced by an external force F(e) simulated by confining the molecules in the region between surfaces which have attractive forces. Conditions at the boundary strongly influence the type of the slip flow predicted. A parabolic velocity profile with apparent slip on the wall is predicted under weakly attractive wall conditions, independent of molecular structure. In the case of strongly attractive walls, a layer of adhered molecules to the wall produces an abrupt distortion of the velocity profile which leads to slip between fluid layers with magnitude that depends on the molecular structure. Finally, the molecular deformation under flow depends on the attractive force of the wall, in such a way that molecules are highly deformed in the case of strong attracting walls.