Water-Tryptophan Interactions: Lone-pair⋅⋅⋅π or O-H⋅⋅⋅π? Molecular Dynamics Simulations of β-Galactosidase Suggest that Both Modes Can Co-exist.

In proteins, the indole side chain of tryptophan can interact with water molecules either in-plane, forming hydrogen bonds, or out-of-plane, with the water molecule contacting the aromatic π face. The latter interaction can be either of the lone pair⋅⋅⋅π (lp⋅⋅⋅π) type or corresponds to the O-H⋅⋅⋅π binding mode, an ambiguity that X-ray structures usually do not resolve. Here, we report molecular dynamics (MD) simulations of a solvated β-galactosidase monomer, which illustrate how a water molecule located at the π face of an indole side chain of tryptophan can adapt to the position of proximate residues and "select" its binding mode. In one such site, the water molecule is predicted to rapidly oscillate between the O-H⋅⋅⋅π and lp⋅⋅⋅π binding modes, thus gaining entropic advantage. Our MD simulations provide support for the role of lp⋅⋅⋅π interactions in the stabilization of protein structures.