Ka Yiu Wong

The Pathway of Oligomeric DNA Melting Investigated by Molecular Dynamics Simulations

Details of the reaction coordinate for DNA melting are fundamental to much of biology and biotechnology. Recently, it has been shown experimentally that there are at least three states involved. To clarify the reaction mechanism of the melting transition of DNA, we perform 100-ns molecular dynamics simulations of a homo-oligomeric, 12-basepair DNA duplex, d(A(12)).d(T(12)), with explicit salt water at 400 K. Analysis of the trajectory reveals the various biochemically important processes that occur on different timescales. Peeling (including fraying from the ends), searching for Watson-Crick complements, and dissociation are recognizable processes. However, we find that basepair searching for Watson-Crick complements along a strand is not mechanistically tied to or directly accessible from the dissociation steps of strand melting. A three-step melting mechanism is proposed where the untwisting of the duplex is determined to be the major component of the reaction coordinate at the barrier. Though the observations are limited to the characteristics of the system being studied, they provide important insight into the mechanism of melting of other more biologically relevant forms of DNA, which will certainly differ in details from those here.

Solubility and Aggregation of Gly5 in Water

Experimentally, the solubility of oligoglycines in water decreases as its length increases. Computationally, the free energy of solvation becomes more favorable with chain length for short (n = 1–5) oligoglycines. We present results of large scale simulations with over 600 pentaglycines at varying concentrations in explicit solvent to consider the mechanism of aggregation. The solubility limit of Gly5 for the force field used was calculated and compared with experimental values. We find that intermolecular interactions between pentaglycines are favored over interactions between glycine and water, leading to their aggregation. However, the interaction driving peptide associations, liquid–liquid phase separation, are not predominantly hydrogen bonding. Instead, non-hydrogen bonding interactions between partially charged atoms on the peptide backbone allow the formation of dipole–dipole and charge layering correlations that mechanistically stabilize the formation of large, stable peptide clusters.

A Non-Watson–Crick Motif of Base-pairing on Surfaces for Untethered Oligonucleotides

A structural view of DNA association/hybridization to a target oligonucleotide molecule near a surface has been developed. Recent experiments have showed a kinetically rapid hybridization between large target DNA fragments and oligonucleotides electrostatically immobilized (untethered) to a surface. Theory and computer simulations have been used to investigate the nature of the specificity and affinity in such a system. Simulations were performed for a modified silicon dioxide surface with positively charged groups at neutral pH. The dosing of a surface with unattached oligonucleotide was simulated. The oligonucleotide was found to associate with the surface in salt water in a way that some of the bases remained stacked, and most of the bases near the surface on average pointed preferentially toward the solution, away from the surface. Use of an analytic solution to the linear Poisson–Boltzmann (PB) theory of the electric double layer interaction between DNA and a hard surface predicts tight binding in this system. The simulation thus gives a mechanism for specificity and the theory a mechanism for affinity. The geometry is such that only non-helical base pairs would be accommodated with an irregular backbone.

A new boundary condition for computer simulations of interfacial systems

A new boundary condition for computer simulations of interfacial systems is presented. The simulation box used in this boundary condition is the asymmetric unit of space group Pb, and it contains only one interface. Compared to the simulation box using common periodic boundary conditions which contains two interfaces, the number of particles in the simulation is reduced by half. This boundary condition was tested against common periodic boundary conditions in molecular dynamic simulations of liquid water interacting with hydroxylated silica surfaces. It yielded results essentially identical to periodic boundary condition and consumed less CPU time for comparable statistics.

Salt Effects on Surface-Tethered Peptides in Solution

The capability to manipulate proteins/peptide fragments at liquid-solid interfaces has led to tremendous applications in detectors and biotechnology. Therefore, understanding the detailed molecular behavior of proteins and peptides tethered on a hard material surface is an interesting and important topic. The inhomogeneity presented by surfaces as well as ions in the solution plays an important role in the thermodynamics and kinetics of the tethered proteins. In this study, we perform a series of molecular dynamics simulations of a pentapeptide RHSVV, a p53 epitope, tethered on a prepared microarray surface in various salt concentrations (0, 0.14, 0.5, and 1 M NaCl), as well as free in ionic solution (0, 0.5, and 1 M). The conformational space the tethered peptide visits largely overlaps with the free peptide in solution. However, surface tethering as well as the salt concentration changes both the thermodynamics and kinetics of the peptide. Frequent conformational changes are observed during the simulations and tend to be slowed down by both increasing the salt concentration and surface tethering. The local composition of ions at different salt concentrations is also compared between the tethered and free peptide.

Transport properties of water at functionalized molecular interfaces

Understanding transport properties of solvent such as diffusion and viscosity at interfaces with biomacromolecules and hard materials is of fundamental importance to both biology and biotechnology. Our study utilizes equilibrium molecular dynamics simulations to calculate solvent transport properties at a model peptide and microarray surface. Both diffusion and selected components of viscosity are considered. Solvent diffusion is found to be affected near the peptide and surface. The stress-stress correlation function of solvent near the hard surface exhibits long time memory. Both diffusion and viscosity are shown to be closely correlated with the density distribution function of water along the microarray surface.

Peptide Solubility Limits: Backbone and Side-Chain Interactions

We calculate the solubility limit of pentapeptides in water by simulating the phase separation in an oversaturated aqueous solution. The solubility limit order followed by our model peptides (GGRGG > GGDGG > GGGGG > GGVGG > GGQGG > GGNGG > GGFGG) is found to be the same as that reported for amino acid monomers from experiment (R > D > G > V > Q > N > F). Investigation of dynamical properties of peptides shows that the higher the solubility of a peptide is, the lower the time spent by the peptide in the aggregated cluster is. We also demonstrate that fluctuations in conformation and hydration number of peptide in monomeric form are correlated with the solubility of the peptide. We considered energetic mechanisms and dynamical properties of interbackbone CO-CO and CO···HN interactions. Our results confirm that CO-CO interactions more than the interbackbone H-bonds are important in peptide self-assembly and association. Further, we find that the stability of H-bonded peptide pairs arises mainly from coexisting CO-CO and CO···HN interactions.