Calculation of Relative Binding Free Energy in the Water-Filled Active Site of Oligopeptide-Binding Protein A.

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.

Free energy calculations give insight into the stereoselective hydroxylation of α-ionones by engineered cytochrome P450 BM3 mutants.

Previously, stereoselective hydroxylation of α-ionone by Cytochrome P450 BM3 mutants M01 A82W and M11 L437N was observed. While both mutants hydroxylate α-ionone in a regioselective manner at the C3 position, M01 A82W catalyzes formation of trans-3-OH-α-ionone products whereas M11 L437N exhibits opposite stereoselectivity, producing trans-(3S,6S)-OH-α-ionone and cis-(3S,6R)-OH-α-ionone. Here, we explore the stereoselective C3 hydroxylation of α-ionone by Cytochrome P450 BM3 mutants M01 A82W and M11 L437N using molecular dynamics-based free energy calculations to study the interaction between the enzyme and both the substrates and the products. The one-step perturbation approach is applied using an optimized reference state for substrates and products. While the free energy differences between the substrates free in solution amount to ~0 kJ mol(-1), the differences in mutant M01 A82W agree with the experimentally obtained dissociation constants K(d). Moreover, a correlation with experimentally observed trends in product formation is found in both mutants. The trans isomers show the most favorable relative binding free energy in the range of all four possible hydroxylated diastereomers for mutant M01 A82W, while the trans product from (6S)-α-ionone and the cis product from (6R)-α-ionone show highest affinity for mutant M11 L437N. Marcus theory is subsequently used to relate the thermodynamic stability to transition state energies and rates of formation.

The role of protein plasticity in computational rationalization studies on regioselectivity in testosterone hydroxylation by cytochrome P450 BM3 mutants.

Recently, it was found that mutations in the binding cavity of drug-metabolizing Cytochrome P450 BM3 mutants can result in major changes in regioselectivity in testosterone (TES) hydroxylation. In the current work, we report the intrinsic reactivity of TES' C-H bonds and our attempts to rationalize experimentally observed changes in TES hydroxylation using a protein structure-based in silico approach, by setting up and employing a combined Molecular Dynamics (MD) and ligand docking approach to account for the flexibility and plasticity of BM3 mutants. Using this approach, about 100,000 TES binding poses were obtained per mutant. The predicted regioselectivity in TES hydroxylation by the mutants was found to be in disagreement with experiment. As revealed in a detailed structural analysis of the obtained docking poses, this disagreement is due to limitations in correctly scoring hydrogen-bonding and steric interactions with specific active-site residues, which could explain the experimentally observed trends in regioselectivity in TES hydroxylation.

Biophysical and physicochemical methods differentiate highly ligand-efficient human D-amino acid oxidase inhibitors.

Many early drug research efforts are too reductionist thereby not delivering key parameters such as kinetics and thermodynamics of target-ligand binding. A set of human D-Amino Acid Oxidase (DAAO) inhibitors 1-6 was applied to demonstrate the impact of key biophysical techniques and physicochemical methods in the differentiation of chemical entities that cannot be adequately distinguished on the basis of their normalized potency (ligand efficiency) values. The resulting biophysical and physicochemical data were related to relevant pharmacodynamic and pharmacokinetic properties. Surface Plasmon Resonance data indicated prolonged target-ligand residence times for 5 and 6 as compared to 1-4, based on the observed k(off) values. The Isothermal Titration Calorimetry-derived thermodynamic binding profiles of 1-6 to the DAAO enzyme revealed favorable contributions of both ΔH and ΔS to their ΔG values. Surprisingly, the thermodynamic binding profile of 3 elicited a substantially higher favorable contribution of ΔH to ΔG in comparison with the structurally closely related fused bicyclic acid 4. Molecular dynamics simulations and free energy calculations of 1, 3, and 4 led to novel insights into the thermodynamic properties of the binding process at an atomic level and in the different thermodynamic signatures of 3 and 4. The presented holistic approach is anticipated to facilitate the identification of compounds with best-in-class properties at an early research stage.

Molecular dynamics simulations and free energy calculations on the enzyme 4-hydroxyphenylpyruvate dioxygenase.

4-Hydroxyphenylpyruvate dioxygenase is a relevant target in both pharmaceutical and agricultural research. We report on molecular dynamics simulations and free energy calculations on this enzyme, in complex with 12 inhibitors for which experimental affinities were determined. We applied the thermodynamic integration approach and the more efficient one-step perturbation. Even though simulations seem well converged and both methods show excellent agreement between them, the correlation with the experimental values remains poor. We investigate the effect of slight modifications on the charge distribution of these highly conjugated systems and find that accurate models can be obtained when using improved force field parameters. This study gives insight into the applicability of free energy methods and current limitations in force field parameterization.

The role of water molecules in computational drug design.

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.

Structural rationalization of novel drug metabolizing mutants of cytochrome P450 BM3.

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.