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.