Similarities and Differences between Crystal and Enzyme Environmental Effects on the Electron Density of Drug Molecules.

The crystal interaction density is generally assumed to be a suitable measure of the polarization of a low-molecular weight ligand inside an enzyme, but this approximation has seldomly been tested and has never been quantified before. In this study, we compare the crystal interaction density and the interaction electrostatic potential for a model compound of loxistatin acid (E64c) with those inside cathepsin B, in solution, and in vacuum. We apply QM/MM calculations and experimental quantum crystallography to show that the crystal interaction density is indeed very similar to the enzyme interaction density. Less than 0.1 e are shifted between these two environments in total. However, this difference has non-negligible consequences for derived properties.

ReFlex3D: Refined Flexible Alignment of Molecules Using Shape and Electrostatics.

We present an algorithm, ReFlex3D, for the refinement of flexible molecular alignments based on their three-dimensional shape and electrostatic properties. The algorithm is designed to be used with fast conformer generators to refine an initial overlay between two molecules and thus to obtain improved overlaps as judged by an increase in calculated similarity values. ReFlex3D is open-source and built as a python package working in combination with the OEChem Toolkit. As such it can readily be implemented in existing workflows ranging from the selection of compounds from a virtual screening campaign to the construction of similarity based prediction models to estimate binding affinities. We evaluate ReFlex3D against the AstraZeneca Validation Test Set and illustrate its potential within a predictive model compared to an established method (Posit).

Discovery and Evaluation of Anti-Fibrinolytic Plasmin Inhibitors Derived from 5-(4-Piperidyl)isoxazol-3-ol (4-PIOL).

Inhibition of plasmin has been found to effectively reduce fibrinolysis and to avoid hemorrhage. This can be achieved by addressing its kringle 1 domain with the known drug and lysine analogue tranexamic acid. Guided by shape similarities toward a previously discovered lead compound, 5-(4-piperidyl)isoxazol-3-ol, a set of 16 structurally similar compounds was assembled and investigated. Successfully, in vitro measurements revealed one compound, 5-(4-piperidyl)isothiazol-3-ol, superior in potency compared to the initial lead. Furthermore, a strikingly high correlation (R2 = 0.93) between anti-fibrinolytic activity and kringle 1 binding affinity provided strong support for the hypothesized inhibition mechanism, as well as revealing opportunities to fine-tune biological effects through minor structural modifications. Several different ligand-based (Freeform, shape, and electrostatic-based similarities) and structure-based methods (e.g., Posit, MM/GBSA, FEP+) were used to retrospectively predict the binding affinities. A combined method, molecular alignment using Posit and scoring with Tcombo, lead to the highest coefficient of determination (R2 = 0.6).

QM/MM investigations of organic chemistry oriented questions.

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

Protocol for rational design of covalently interacting inhibitors.

The inhibition potencies of covalent inhibitors mainly result from the formation of a covalent bond to the enzyme during the inhibition mechanism. This class of inhibitors has essentially been ignored in previous target-directed drug discovery projects because of concerns about possible side effects. However, their advantages, such as higher binding energies and longer drug-target residence times moved them into the focus of recent investigations. While the rational design of non-covalent inhibitors became standard the corresponding design of covalent inhibitors is still in its early stages. Potent covalent inhibitors can be retrieved from large compound libraries by covalent docking approaches but protocols are missing that can reliably predict the influence of variations in the substitution pattern on the affinity and/or reactivity of a given covalent inhibitor. Hence, the wanted property profile can only be obtained from trial-and-error proceedings. This paper presents an appropriate protocol which is able to predict improved covalent inhibitors. It uses hybrid approaches, which mix quantum mechanical (QM) and molecular mechanical (MM) methods to predict variations in the reactivity of the inhibitor. They are also used to compute the required information about the non-covalent enzyme-inhibitor complex. Docking tools are employed to improve the inhibitor with respect to the non-covalent interactions formed in the binding site.