Although water molecules are small and only consist of two different atom types, they play various roles in cellular systems. This review discusses their influence on the binding process between biomacromolecular targets and small molecule ligands and how this influence can be modeled in computational drug design approaches. Both the structure and the thermodynamics of active site waters will be discussed as these influence the binding process significantly. Structurally conserved waters cannot always be determined experimentally and if observed, it is not clear if they will be replaced upon ligand binding, even if sufficient space is available. Methods to predict the presence of water in protein-ligand complexes will be reviewed. Subsequently, we will discuss methods to include water in computational drug research. Either as an additional factor in automated docking experiments, or explicitly in detailed molecular dynamics simulations, the effect of water on the quality of the simulations is significant, but not easily predicted. The most detailed calculations involve estimates of the free energy contribution of water molecules to protein-ligand complexes. These calculations are computationally demanding, but give insight in the versatility and importance of water in ligand binding.

Inversion of stereoselectivity: screening of a minimal mutant library revealed a cytochrome P450 BM3 variant M01 A82W S72I capable of producing 16 α-OH-testosterone. Remarkably, a single active site mutation S72I in M01 A82W inverted the stereoselectivity of hydroxylation from 16 β to 16 α. Introduction of S72I mutation in another 16 β-OH-selective variant M11 V87I, also resulted in similar inversion of stereoselectivity.

Abstract:The periplasmic oligopeptide binding protein A (OppA) represents a well-known example of water-mediated protein-ligand interactions. Here, we perform free-energy calculations for three different ligands binding to OppA, using a thermodynamic integration approach. The tripeptide ligands share a high structural similarity (all have the sequence KXK), but their experimentally-determined binding free energies differ remarkably. Thermodynamic cycles were constructed for the ligands, and simulations conducted in the bound and (freely solvated) unbound states. In the unbound state, it was observed that the difference in conformational freedom between alanine and glycine leads to a surprisingly slow convergence, despite their chemical similarity. This could be overcome by increasing the softness parameter during alchemical transformations. Discrepancies remained in the bound state however, when comparing independent simulations of the three ligands. These difficulties could be traced to a slow relaxation of the water network within the active site. Fluctuations in the number of water molecules residing in the binding cavity occur mostly on a timescale larger than the simulation time along the alchemical path. After extensive simulations, relative binding free energies that were converged to within thermal noise could be obtained, which agree well with available experimental data.

Abstract:The selective hydroxylation of an unactivated C À H bond is a crucial step in the synthesis of fine chemicals such as hydroxylated terpenoids. In the present study, the ability of 40 cytochrome P450 BM3 mutants to perform the regio-and stereoselective hydroxylation of a-ionone has been investigated. Based on their activity and selectivity to produce 3-hydroxy-a-ionone from racemic a-ionone, 6 BM3 mutants were selected. Out of these, 3 mutants (M01 A82W, M11 A82W and M11 V87I) showed high selectivity for trans-3-hydroxy-a-ionone formation while 3 other mutants (M11 L437N, M11 L437S and M11 L437T) formed almost equal amounts of both cis-3-hydroxy-and trans-3-hydroxy-a-ionone. Incubation with individual enantiomers showed that M11 L437N, M11 L437S and M11 L437T exhibited opposite stereoselectivity producing (3S,6S)-hydroxya-ionone with the (6S)-enantiomer and (3S,6R)-hydroxy-a-ionone with the (6R)-enantiomer. Thus for the first time, BM3 mutants that can selectively produce diastereomers of 3-hydroxy-a-ionone (> 90% de), with high turnover numbers and minimal secondary metabolism, have been identified. Docking studies have been performed to rationalize the basis of the experimentally observed selectivity. In conclusion, engineered P450 BM3s are promising biocatalysts for regio-and stereoselective production of hydroxylated a-ionones for industrial applications.

Three newly discovered drug metabolizing mutants of cytochrome P450 BM3 (van Vugt-Lussenburg et al., Identification of critical residues in novel drug metabolizing mutants of Cytochrome P450 BM3 using random mutagenesis, J Med Chem 2007;50:455-461) have been studied at an atomistic level to provide structural explanations for a number of their characteristics. In this study, computational methods are combined with experimental techniques. Molecular dynamics simulations, resonance Raman and UV-VIS spectroscopy, as well as coupling efficiency and substrate-binding experiments, have been performed. The computational findings, supported by the experimental results, enable structural rationalizations of the mutants. The substrates used in this study are known to be metabolized by human cytochrome P450 2D6. Interestingly, the major metabolites formed by the P450 BM3 mutants differ from those formed by human cytochrome P450 2D6. The computational findings, supported by resonance Raman data, suggest a conformational change of one of the heme propionate groups. The modeling results furthermore suggest that this conformational change allows for an interaction between the negatively charged carboxylate of the heme substituent and the positively charged nitrogen of the substrates. This allows for an orientation of the substrates favorable for formation of the major metabolite by P450 BM3.