ADME prediction is an extremely challenging area as many of the properties we try to predict are a result of multiple physiological processes. In this review we consider how in-silico predictions of ADME processes can be used to help bias medicinal chemistry into more ideal areas of property space, minimizing the number of compounds needed to be synthesized to obtain the required biochemical/physico-chemical profile. While such models are not sufficiently accurate to act as a replacement for in-vivo or in-vitro methods, in-silico methods nevertheless can help us to understand the underlying physico-chemical dependencies of the different ADME properties, and thus can give us inspiration on how to optimize them. Many global in-silico ADME models (i.e generated on large, diverse datasets) have been reported in the literature. In this paper we selectively review representatives from each distinct class and discuss their relative utility in drug discovery. For each ADME parameter, we limit our discussion to the most recent, most predictive or most insightful examples in the literature to highlight the current state of the art. In each case we briefly summarize the different types of models available for each parameter (i.e simple rules, physico-chemical and 3D based QSAR predictions), their overall accuracy and the underlying SAR. We also discuss the utility of the models as related to lead generation and optimization phases of discovery research.

Kalata B1 has been demonstrated to have bioactivity relating to membrane disruption. In this study, we conducted coarse-grained molecular dynamics simulations to gain further insight into kB1 bioactivity. The simulations were performed at various concentrations of kB1 to capture the overall progression of its activity. Two configurations of kB1 oligomers, termed tower-like and wall-like clusters, were detected. The conjugation between the wall-like oligomers resulted in the formation of a ring-like hollow in the kB1 cluster on the membrane surface. Our results indicated that the molecules of kB1 were trapped at the membrane-water interface. The interfacial membrane binding of kB1 induced a positive membrane curvature, and the lipids were eventually extracted from the membrane through the kB1 ring-like hollow into the space inside the kB1 cluster. These findings provide an alternative view of the mechanism of kB1 bioactivity that corresponds with the concept of an interfacial bioactivity model.

In 2004, we used NMR to solve the structure of the minor groove binder thiazotropsin A bound in a 2:1 complex to the DNA duplex, d(CGACTAGTCG)2. In this current work, we have combined theory and experiment to confirm the binding thermodynamics of this system. Molecular dynamics simulations that use polarizable or non-polarizable force fields with single and separate trajectory approaches have been used to explore complexation at the molecular level. We have shown that the binding process invokes large conformational changes in both the receptor and ligand, which is reflected by large adaptation energies. This is compensated for by the net binding free energy, which is enthalpy driven and entropically opposed. Such a conformational change upon binding directly impacts on how the process must be simulated in order to yield accurate results. Our MM-PBSA binding calculations from snapshots obtained from MD simulations of the polarizable force field using separate trajectories yield an absolute binding free energy (-15.4 kcal mol(-1)) very close to that determined by isothermal titration calorimetry (-10.2 kcal mol(-1)). Analysis of the major energy components reveals that favorable non-bonded van der Waals and electrostatic interactions contribute predominantly to the enthalpy term, whilst the unfavorable entropy appears to be driven by stabilization of the complex and the associated loss of conformational freedom. Our results have led to a deeper understanding of the nature of side-by-side minor groove ligand binding, which has significant implications for structure-based ligand development.

No abstract

No abstract

Two-layered and three-layered ONIOM calculations were performed to compare the binding energies of 8-Cl TIBO inhibitor when bound into the human immunodeficiency virus reverse transcriptase binding pocket and a Y181C variant. Both consisted of 20 residues within a radius of 15 A. A combination of different methods [MP2/6-31G(d), B3LYP/6-31G(d,p), and PM3] were performed to take advantage of ONIOM's layering strategy analysis. The obtained results clearly indicate that the Y181C mutation reduces the binding affinity and stability of the inhibitor by approximately 8-9 kcal/mol as obtained from different combined MO:MO methods. Analyses regarding the energetic components of the interaction and deformation energies for 8-Cl TIBO inhibitor upon binding were also examined extensively. Additional calculations involving the interaction energies between 8-Cl TIBO with individual residues surrounding the binding pocket were performed at MP2/6-31G(d,p) and B3LYP/6-31G(d,p) levels of theory to gain more insight into the energetic differences of wild-type and Y181C mutant type at the atomistic level.

Comparative molecular field analysis (CoMFA) has been applied to a large set of 1-[(2-hydroxyethoxy)methyl]-6-(phenylthio)thymine (HEPT) analogues. The starting geometry of HEPT was obtained from crystallographic data of HEPT/HIV-1 reverse transcriptase (RT) complexes. The structures of 101 HEPT derivatives were considered and fully optimized by ab initio molecular orbital calculations at the HF/3-21G level. The best CoMFA model is satisfactory in both statistical significance and predictive ability. It shows excellent, high predictive ability as r2cv = 0.858. The derived model indicates the importance of steric contributions (64.4%) as well as electrostatic interactions for the HIV-1 RT inhibition. In addition, steric and electrostatic contour maps from this analysis agree well with the experimentally observed trend that there are steric interactions between the side chain of HEPT and an aromatic ring of Tyr181. It is concluded that a moderately sized group at C5 enhances contact with Tyr181 enough to push it into a position which renders the protein nonfunctional, but a smaller group has insufficient steric requirements to do this and a larger group renders the ligand too large for the cavity. The mutation-induced resistance of reverse transcriptase is explained by this analysis. The obtained results not only lead to a better understanding of structural requirements of this set of compounds for the inhibition but also enable the suggestions for new and more potent drugs.

No abstract