MOLECULAR DYNAMICS STUDIES OF SIMPLE MODEL FLUIDS AND WATER CONFINED IN CARBON NANOTUBE

Molecular Dynamics (MD) simulation is one of the most important computational techniques with broad applications in physics, chemistry, chemical engineering, materials design and biological science. Traditional computational chemistry refers to quantum calculations based on solving Schrodinger equations. Later developed Density Functional Theory (DFT) based on solving Kohn-Sham equations became the more popular ab initio calculation technique which could deal with ~1000 atoms by explicitly considering electron interactions. In contrast, MD simulation based on solving classical mechanics equations of motion is a totally different technique in the field of computational chemistry. Electron interactions were implicitly included in the empirical atom-based potential functions and the system size to be investigated can be extended to ~106 atoms. The thermodynamic properties of model fluids are mainly determined by macroscopic quantities, like temperature, pressure, density. The quantum effects on thermodynamic properties like melting point, surface tension are not dominant. In this work, we mainly investigated the melting point, surface tension (liquid-vapor and liquid-solid) of model fluids including Lennard-Jones model, Stockmayer model and a couple of water models (TIP4P/Ew, TIP5P/Ew) by means of MD simulation. In addition, some new structures of water confined in carbon nanotube were discovered and transport behaviors of water and ions through nano-channels were also revealed.

A macroscopic kinetic model for DNA polymerase elongation and high-fidelity nucleotide election

The enzymatically catalyzed template-directed extension of ssDNA/primer complex is an impor-tant reaction of extraordinary complexity. The DNA polymerase does not merely facilitate the insertion of dNMP, but it also performs rapid screening of substrates to ensure a high degree of fidelity. Several kinetic studies have determined rate constants and equilibrium constants for the elementary steps that make up the overall pathway. The information is used to develop a macro-scopic kinetic model, using an approach described by Ninio [Ninio J., 1987. Alternative to the steady-state method: derivation of reaction rates from first-passage times and pathway probabili-ties. Proc. Natl. Acad. Sci. U.S.A. 84, 663–667]. The principle idea of the Ninio approach is to track a single template/primer complex over time and to identify the expected behavior. The average time to insert a single nucleotide is a weighted sum of several terms, in-cluding the actual time to insert a nucleotide plus delays due to polymerase detachment from ei-ther the ternary (template-primer-polymerase) or quaternary (+nucleotide) complexes and time delays associated with the identification and ultimate rejection of an incorrect nucleotide from the binding site. The passage times of all events and their probability of occurrence are ex-pressed in terms of the rate constants of the elementary steps of the reaction pathway. The model accounts for variations in the average insertion time with different nucleotides as well as the in-fluence of G+C content of the sequence in the vicinity of the insertion site. Furthermore the model provides estimates of error frequencies. If nucleotide extension is recognized as a compe-tition between successful insertions and time delaying events, it can be described as a binomial process with a probability distribution. The distribution gives the probability to extend a primer/template complex with a certain number of base pairs and in general it maps annealed complexes into extension products.