Timothy W. Sirk

Molecular binding of catechins to biomembranes: relationship to biological activity.

Molecular dynamics simulations were used to study the interactions of four green tea catechin compounds with lipid bilayers. Reported studies have shown that catechins are linked to beneficial health effects, specifically those related to interactions with the cell membrane. To better understand the molecular interaction of catechins with membranes, simulations were carried out of interactions of four catechin molecules [epicatechin (EC), epigallocatechin (EGC), epicatechin gallate (ECG), and epigallocatechin gallate (EGCG)] with a 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) lipid bilayer. The simulations show that catechins possess a strong affinity for the lipid bilayer. Some are absorbed into the bilayer. The molecular structure and aggregated condition of the catechins significantly influences their absorption, as well as their ability to form hydrogen bonds with the lipid headgroups. Insight into these molecular interactions helps to distinguish the structure-function relationship of the catechins with lipid bilayers and provides a foundation for a better understanding of the role of catechins in biological processes.

Molecular Dynamics Study on the Biophysical Interactions of Seven Green Tea Catechins with Lipid Bilayers of Cell Membranes

Molecular dynamics simulations were performed to study the interactions of bioactive catechins (flavonoids) commonly found in green tea with lipid bilayers, as a model for cell membranes. Previously, multiple experimental studies rationalized catechins anticarcinogenic, antibacterial, and other beneficial effects in terms of physicochemical molecular interactions with the cell membranes. To contribute toward understanding the molecular role of catechins on the structure of cell membranes, we present simulation results for seven green tea catechins in lipid bilayer systems representative of HepG2 cancer cells. Our simulations show that the seven tea catechins evaluated have a strong affinity for the lipid bilayer via hydrogen bonding to the bilayer surface, with some of the smaller catechins able to penetrate underneath the surface. Epigallocatechin-gallate (EGCG) showed the strongest interaction with the lipid bilayer based on the number of hydrogen bonds formed with lipid headgroups. The simulations also provide insight into the functional characteristics of the catechins that distinguish them as effective compounds to potentially alter the lipid bilayer properties. The results on the hydrogen-bonding effects, described here for the first time, may contribute to a better understanding of proposed multiple molecular mechanisms of the action of catechins in microorganisms, cancer cells, and tissues.

Characteristics of thermal conductivity in classical water models

The thermal conductivities of common water models are compared using equilibrium (EMD) and non-equilibrium molecular dynamics (NEMD) simulation. A complete accounting for electrostatic contributions to the heat flux was found to resolve the previously reported differing results of NEMD and EMD Green-Kubo measurements for the extended simple point-charge (SPC/E) model. Accordingly, we demonstrate the influence of long-range electrostatics on the thermal conductivity with a simple coulomb cutoff, Ewald summation, and by an extended particle-particle particle-mesh method. For each water model, the thermal conductivity is computed and decomposed in terms of frequency-dependent thermodynamic and topological contributions. The rigid, three-site SPC, SPC/E, and transferable intermolecular potential (TIP3P-Ew) water models are shown to have similar thermal conductivity values at standard conditions, whereas models that include bond stretching and angle bending have higher thermal conductivities.

Molecular binding of black tea theaflavins to biological membranes: relationship to bioactivities.

Molecular dynamics simulations were used to study the interactions of three theaflavin compounds with lipid bilayers. Experimental studies have linked theaflavins to beneficial health effects, some of which are related to interactions with the cell membrane. The molecular interaction of theaflavins with membranes was explored by simulating the interactions of three theaflavin molecules (theaflavin, theaflavin-3-gallate, and theaflavin-3,3-digallate) with a mixed bilayer composed of 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC) and 1-palmitoyl-2-oleoyl phosphatidylethanolamine (POPE). The simulations show that the theaflavins evaluated have an affinity for the lipid bilayer surface via hydrogen bonding. The molecular structure of theaflavins influenced their configuration when binding to the bilayer surface, as well as their ability to form hydrogen bonds with the lipid headgroups. The theaflavin-bilayer interactions studied here help to define structure-function relationships of the theaflavins and provide a better understanding of the role of theaflavins in biological processes. The significance of the results are discussed in the context of black tea composition and bioactivity.

Entangled triblock copolymer gel: morphological and mechanical properties.

The morphological and mechanical properties of entangled ABA triblock copolymer gels, where solvent were selective to the midblock, were studied as a function of polymer concentration using a novel dissipative particle dynamics model which includes a modified segmental repulsive potential that restricts chain crossing. Morphological properties, such as micelle size, distance between micelles, and the bridge fraction, were calculated as a function of concentration. Although the micelle size was shown to have a strong dependence on concentration, the bridge fraction and distance between micelles were shown to plateau at moderate concentrations. Deformation under uni-axial tension was also performed to extract the cross-link and entanglement contribution to the elastic modulus. Scaling results qualitatively agree with other theoretical predications.

A Study of Nanoparticle-Enhanced Heat Transfer

The study of the heat transfer characteristics of nanofluids, i.e. fluids that are suspensions of nanometer size particles, has gained significant attention in the search for new coolants that can effectively service a variety of needs ranging from the increasing heat transfer demands of ever smaller microelectronic devices to mitigating the effects of loss of coolant accidents in nuclear power plants. Experimental data has shown large increases in thermal conductivity and associated increases in the level of critical heat flux in nuclear reactors; however, in some cases the range of the applicability of the experimental results is uncertain and there is a lack of a theory by which this can be resolved. Complicating the theoretical description of heat transfer in nanofluids is the fact that fluids in the vicinity of the nanoparticles are a complex combination of phase transition, interfacial, and transport phenomena. This paper describes a study in which molecular dynamics simulations were used to enhance the understanding of the effect of nanoparticles on heat transfer. The molecular dynamics (MD) simulations presented here model a Lennard-Jones fluid in a channel where the walls are maintained at different temperatures. The heat flux is calculated for a variety of nanoparticle sizes and concentrations. The results are compared to experimental data in order to provide information that will more confidently bound the data and provide information that will guide the development of more comprehensive theories. We also anticipate that this work could contribute to the design of biosensors where suspended molecules are transported through micro- and nano-channels in the presence of heat transfer.Copyright

A Diffusion-Based Approach for Nanobead Drag Forces

A central process of the Bead ARray Counter (BARC) biosensor involves propelling a dilute solution of streptavidin-coated micron-scale beads through water. It is desirable to downsize these beads, leading to a basic engineering problem: What is the drag force acting on a nanoscale bead having an unknown hydrodynamic boundary condition? We use molecular dynamics simulation to estimate the drag force acting on a hypothetical biosensor “nanobead” without assuming a no-slip condition or a user-supplied slip parameter.Copyright