Brian Novak

Molecular Dynamics Simulation Study of the Effect of DMSO on Structural and Permeation Properties of DMPC Lipid Bilayers

We present molecular dynamics simulations of dimyristoylphosphatidylcholine (DMPC) lipid bilayers in the presence of dimethyl sulfoxide (DMSO). The MD simulations focus on understanding the effect of 3 mol % DMSO on structural and permeation properties of DMPC bilayers. The potential of mean force (PMF) and the diffusivity profiles along the normal direction to the bilayer were calculated for water and DMSO molecules in systems containing 0 and 3 mol % DMSO. The simulation results indicate that while the presence of DMSO has only a small effect on diffusion coefficients of both water and DMSO molecules, it affects significantly the corresponding trans-membrane free energy profiles. Using the free energy profiles and diffusivities for water and DMSO and by employing an inhomogeneous solubility-diffusion model we calculated the permeability coefficients. Our simulations show that the increase of the concentration of DMSO in the solution to 3 mol % leads to a significant increase, by about 3 times, of the permeability of water through a DMPC bilayer; a permeability increase that might explain in part the unusual ability of DMSO, even at relatively low concentrations, to reduce the osmotic pressure imbalance present during cryopreservation protocols.

Electrophoretic Transport of Single DNA Nucleotides through Nanoslits: A Molecular Dynamics Simulation Study.

There is potential for flight time based DNA sequencing involving disassembly into individual nucleotides which would pass through a nanochannel with two or more detectors. We performed molecular dynamics simulations of electrophoretic motion of single DNA nucleotides through 3 nm wide hydrophobic slits with both smooth and rough walls. The electric field (E) varied from 0.0 to 0.6 V/nm. The nucleotides adsorb and desorb from walls multiple times during their transit through the slit. The nucleotide-wall interactions differed due to nucleotide hydrophobicities and wall roughness which determined duration and frequency of nucleotide adsorptions and their velocities while adsorbed. Transient association of nucleotides with one, two, or three sodium ions occurred, but the mean association numbers (ANs) were weak functions of nucleotide type. Nucleotide-wall interactions contributed more to separation of nucleotide flight time distributions than ion association and thus indicate that nucleotide-wall interactions play a defining role in successfully discriminating between nucleotides on the basis of their flight times through nanochannels/slits. With smooth walls, smaller nucleotides moved faster, but with rough walls larger nucleotides moved faster due to fewer favorable wall adsorption sites. This indicates that roughness, or surface patterning, might be exploited to achieve better time-of-flight based discrimination between nucleotides.

Behavior of the ATP grasp domain of biotin carboxylase monomers and dimers studied using molecular dynamics simulations

The enzyme biotin carboxylase (BC) uses adenosine triphosphate (ATP) to carboxylate biotin and is involved in fatty acid synthesis. Structural evidence suggests that the B domain of BC undergoes a large hinge motion of ∼45° when binding and releasing substrates. Escherichia coli BC can function as a natural homodimer and as a mutant monomer. Using molecular dynamics simulations, we evaluate the free energy profile along a closure angle of the B domain of E. coli BC for three cases: a monomer without bound Mg2ATP, a monomer with bound Mg2ATP, and a homodimer with bound Mg2ATP in one subunit. The simulation results show that a closed state is the most probable for the monomer with or without bound Mg2ATP. For the dimer with Mg2ATP in one of its subunits, communication between the two subunits was observed. Specifically, in the dimer, the opening of the subunit without Mg2ATP caused the other subunit to open, and hysteresis was observed upon reclosing it. The most stable state of the dimer is one in which the B domain of both subunits is closed; however, the open state for the B domain without Mg2ATP is only approximately 2kBT higher in free energy than the closed state. A simple diffusion model indicates that the mean times for opening and closing of the B domain in the monomer with and without Mg2ATP are much smaller than the overall reaction time, which is on the order of seconds. Proteins 2011.

Umbrella sampling simulations of biotin carboxylase: is a structure with an open ATP grasp domain stable in solution?

Biotin carboxylase is a homodimer that utilizes ATP to carboxylate biotin. Studies of the enzyme using X-ray crystallography revealed a prominent conformational change upon binding ATP. To determine the importance of this closing motion, the potential of mean force with the closure angle as a reaction coordinate was calculated using molecular dynamics simulations and umbrella sampling for a monomer of Escherichia coli biotin carboxylase in water with restraints to simulate attachment to a surface. The result suggests that the most stable state for the enzyme is a closed state different from both the ATP-bound and open state X-ray crystallography structures. There is also a significant motion of a region near the dimer interface not predicted by considering only open and closed configurations, which may have implications for the dynamics and activity of the dimer.

Multi-species Fluid Flow Simulations Using a Hybrid Computational Fluid Dynamics - Molecular Dynamics Approach

The constrained Lagrangian dynamics modeling in the hybrid computational fluid dynamics (CFD) molecular dynamics (MD) approach is improved for the simulation of multi-species polyatomic fluid. The primitive formulation of the classical Lagrangian dynamics equation is replaced by conservative form to account for multi-species fluid system. Also, the equation is applied on molecules instead of individual atom, to preserve the linear momentum between continuum and particle domain without encountering the unfavorable numerical break-down of molecular bonding. We verify our hybrid CFD-MD simulation package by analyzing a nano-scale transient Couette flow of a single monatomic fluid. The multi-species polyatomic Lagrangian dynamics modeling has been evaluated by analyzing two different fluid models: the mixture of two monatomic fluids and a polyatomic molecular fluid under the short-range potential. These two applications verify the accuracy of the proposed model and evaluate the hybrid CFD-MD approach as a tool to describe the complex flow field near the solid obstacle.

The role of the asymmetric bolaamphiphilic character of VECAR on the kinetic and structural aspects of its self-assembly: A molecular dynamics simulation study

VECAR are novel bolaamphiphilic molecules consisting of two hydrophilic molecular groups, a carnosine derivative and a chromanol group, covalently linked by a hydrophobic alkyl spacer of varying length. Despite the potential for application in various biomedical applications VECAR properties, including their bulk properties, are still largely unknown. The early stage of the self-assembly process of VECAR molecules in water is studied using molecular dynamics simulations. The study reveals that the length of the hydrophobic spacer in VECAR affects the aggregation kinetics as well as the size, shape, density, and atomistic structure of the self-assembled aggregates. A mechanism based on cooperative interactions between water, the hydrophilic hydroxyl group, and the hydrophobic benzene ring of the chromanol head is proposed to explain the ordered packings of chromanols in the self-assembled aggregate structures at the aggregate-water interface.