Validation of tautomeric and protomeric binding modes by free energy calculations. A case study for the structure based optimization of D-amino acid oxidase inhibitors

Optimization of fragment size D-amino acid oxidase (DAAO) inhibitors was investigated using a combination of computational and experimental methods. Retrospective free energy perturbation (FEP) calculations were performed for benzo[d]isoxazole derivatives, a series of known inhibitors with two potential binding modes derived from X-ray structures of other DAAO inhibitors. The good agreement between experimental and computed binding free energies in only one of the hypothesized binding modes strongly support this bioactive conformation. Then, a series of 1-H-indazol-3-ol derivatives formerly not described as DAAO inhibitors was investigated. Binding geometries could be reliably identified by structural similarity to benzo[d]isoxazole and other well characterized series and FEP calculations were performed for several tautomers of the deprotonated and protonated compounds since all these forms are potentially present owing to the experimental pKa values of representative compounds in the series. Deprotonated compounds are proposed to be the most important bound species owing to the significantly better agreement between their calculated and measured affinities compared to the protonated forms. FEP calculations were also used for the prediction of the affinities of compounds not previously tested as DAAO inhibitors and for a comparative structure-activity relationship study of the benzo[d]isoxazole and indazole series. Selected indazole derivatives were synthesized and their measured binding affinity towards DAAO was in good agreement with FEP predictions.

Force generation by titin folding

Titin is a giant protein that provides elasticity to muscle. As the sarcomere is stretched, titin extends hierarchically according to the mechanics of its segments. Whether titin's globular domains unfold during this process and how such unfolded domains might contribute to muscle contractility are strongly debated. To explore the force-dependent folding mechanisms, here we manipulated skeletal-muscle titin molecules with high-resolution optical tweezers. In force-clamp mode, after quenching the force (<10 pN), extension fluctuated without resolvable discrete events. In position-clamp experiments, the time-dependent force trace contained rapid fluctuations and a gradual increase of average force, indicating that titin can develop force via dynamic transitions between its structural states en route to the native conformation. In 4 M urea, which destabilizes H-bonds hence the consolidated native domain structure, the net force increase disappeared but the fluctuations persisted. Thus, whereas net force generation is caused by the ensemble folding of the elastically-coupled domains, force fluctuations arise due to a dynamic equilibrium between unfolded and molten-globule states. Monte-Carlo simulations incorporating a compact molten-globule intermediate in the folding landscape recovered all features of our nanomechanics results. The ensemble molten-globule dynamics delivers significant added contractility that may assist sarcomere mechanics, and it may reduce the dissipative energy loss associated with titin unfolding/refolding during muscle contraction/relaxation cycles.

Binding thermodynamics discriminates fragments from druglike compounds: a thermodynamic description of fragment-based drug discovery

Small is beautiful - reducing the size and complexity of chemical starting points for drug design allows better sampling of chemical space, reveals the most energetically important interactions within protein-binding sites and can lead to improvements in the physicochemical properties of the final drug. The impact of fragment-based drug discovery (FBDD) on recent drug discovery projects and our improved knowledge of the structural and thermodynamic details of ligand binding has prompted us to explore the relationships between ligand-binding thermodynamics and FBDD. Information on binding thermodynamics can give insights into the contributions to protein-ligand interactions and could therefore be used to prioritise compounds with a high degree of specificity in forming key interactions.

Ensemble docking-based virtual screening yields novel spirocyclic JAK1 inhibitors

Small molecule inhibition of Janus kinases (JAKs) has been demonstrated as a viable strategy for the treatment of various inflammatory conditions and continues to emerge in cancer-related indications. In this study, a large supplier database was screened to identify novel chemistry starting points for JAK1. The docking-based screening was followed up by testing ten hit compounds experimentally, out of which five have displayed single-digit micromolar and submicromolar IC₅₀ values on JAK1. Thus, the study was concluded with the discovery of five novel JAK inhibitors from a tiny screening deck with a remarkable hitrate of 50%. The results have highlighted spirocyclic pyrrolopyrimidines with submicromolar JAK1 IC₅₀ values and a preference for JAK1 over JAK2 as potential starting points in developing a novel class of JAK1 inhibitors.

Glossary of terms used in computational drug design, part II (IUPAC Recommendations 2015)

Computational drug design is a rapidly changing field that plays an increasingly important role in medicinal chemistry. Since the publication of the first glossary in 1997, substantial changes have occurred in both medicinal chemistry and computational drug design. This has resulted in the use of many new terms and the consequent necessity to update the previous glossary. For this purpose a Working Party of eight experts was assembled. They produced explanatory definitions of more than 150 new and revised terms.

Structure-Based Consensus Scoring Scheme for Selecting Class A Aminergic GPCR Fragments

Aminergic G-protein coupled receptors (GPRCs) represent well-known targets of central nervous-system related diseases. In this study a structure-based consensus virtual screening scheme was developed for designing targeted fragment libraries against class A aminergic GPCRs. Nine representative aminergic GPCR structures were selected by first clustering available X-ray structures and then choosing the one in each cluster that performs best in self-docking calculations. A consensus scoring protocol was developed using known promiscuous aminergic ligands and decoys as a training set. The consensus score (FrACS-fragment aminergic consensus score) calculated for the optimized protein ensemble showed improved enrichments in most cases as compared to stand-alone structures. Retrospective validation was carried out on public screening data for aminergic targets (5-HT1 serotonin receptor, TA1 trace-amine receptor) showing 8-17-fold enrichments using an ensemble of aminergic receptor structures. The performance of the structure based FrACS in combination with our ligand-based prefilter (FrAGS) was investigated both in a retrospective validation on the ChEMBL database and in a prospective validation on an in-house fragment library. In prospective validation virtual fragment hits were tested on 5-HT6 serotonin receptors not involved in the development of FrACS. Six out of the 36 experimentally tested fragments exhibited remarkable antagonist efficacies, and 4 showed IC50 values in the low micromolar or submicromolar range in a cell-based assay. Both retrospective and prospective validations revealed that the methodology is suitable for designing focused class A GPCR fragment libraries from large screening decks, commercial compound collections, or virtual databases.

On the enthalpic preference of fragment binding

Detailed thermodynamic analysis of fragment binding revealed that unlike drug-like compounds, fragments bind with significant enthalpic preference.

Property-based characterization of kinase-like ligand space for library design and virtual screening

A property-based desirability scoring scheme has been developed for kinase-focused library design and ligand-based pre-screening of large compound sets.

Calculation of wave-functions with frozen orbitals in mixed quantum mechanics/molecular mechanics methods. Part I. Application of the Huzinaga equation

Mixed quantum mechanics/quantum mechanics (QM/QM) and quantum mechanics/molecular mechanics (QM/MM) methods make computations feasible for extended chemical systems by separating them into subsystems that are treated at different level of sophistication. In many applications, the subsystems are covalently bound and the use of frozen localized orbitals at the boundary is a possible way to separate the subsystems and to ensure a sensible description of the electronic structure near to the boundary. A complication in these methods is that orthogonality between optimized and frozen orbitals has to be warranted and this is usually achieved by an explicit orthogonalization of the basis set to the frozen orbitals. An alternative to this approach is proposed by calculating the wave-function from the Huzinaga equation that guaranties orthogonality to the frozen orbitals without basis set orthogonalization. The theoretical background and the practical aspects of the application of the Huzinaga equation in mixed methods are discussed. Forces have been derived to perform geometry optimization with wave-functions from the Huzinaga equation. Various properties have been calculated by applying the Huzinaga equation for the central QM subsystem, representing the environment by point charges and using frozen strictly localized orbitals to connect the subsystems. It is shown that a two to three bond separation of the chemical or physical event from the frozen bonds allows a very good reproduction (typically around 1 kcal/mol) of standard Hartree-Fock-Roothaan results. The proposed scheme provides an appropriate framework for mixed QM/QM and QM/MM methods.

Calculation of wave-functions with frozen orbitals in mixed quantum mechanics/molecular mechanics methods. II. Application of the local basis equation

The application of the local basis equation (Ferenczy and Adams, J. Chem. Phys. 2009, 130, 134108) in mixed quantum mechanics/molecular mechanics (QM/MM) and quantum mechanics/quantum mechanics (QM/QM) methods is investigated. This equation is suitable to derive local basis nonorthogonal orbitals that minimize the energy of the system and it exhibits good convergence properties in a self-consistent field solution. These features make the equation appropriate to be used in mixed QM/MM and QM/QM methods to optimize orbitals in the field of frozen localized orbitals connecting the subsystems. Calculations performed for several properties in divers systems show that the method is robust with various choices of the frozen orbitals and frontier atom properties. With appropriate basis set assignment, it gives results equivalent with those of a related approach [G. G. Ferenczy previous paper in this issue] using the Huzinaga equation. Thus, the local basis equation can be used in mixed QM/MM methods with small size quantum subsystems to calculate properties in good agreement with reference Hartree-Fock-Roothaan results. It is shown that bond charges are not necessary when the local basis equation is applied, although they are required for the self-consistent field solution of the Huzinaga equation based method. Conversely, the deformation of the wave-function near to the boundary is observed without bond charges and this has a significant effect on deprotonation energies but a less pronounced effect when the total charge of the system is conserved. The local basis equation can also be used to define a two layer quantum system with nonorthogonal localized orbitals surrounding the central delocalized quantum subsystem.

Thermodynamics of fragment binding

The ligand binding pockets of proteins have preponderance of hydrophobic amino acids and are typically within the apolar interior of the protein; nevertheless, they are able to bind low complexity, polar, water-soluble fragments. In order to understand this phenomenon, we analyzed high resolution X-ray data of protein-ligand complexes from the Protein Data Bank and found that fragments bind to proteins with two near optimal geometry H-bonds on average. The linear extent of the fragment binding site was found not to be larger than 10 Å, and the H-bonding region was found to be restricted to about 5 Å on average. The number of conserved H-bonds in proteins cocrystallized with multiple different fragments is also near to 2. These fragment binding sites that are able to form limited number of strong H-bonds in a hydrophobic environment are identified as hot spots. An estimate of the free-energy gain of H-bond formation versus apolar desolvation supports that fragment sized compounds need H-bonds to achieve detectable binding. This suggests that fragment binding is mostly enthalpic that is in line with their observed binding thermodynamics documented in Isothermal Titration Calorimetry (ITC) data sets and gives a thermodynamic rationale for fragment based approaches. The binding of larger compounds tends to more rely on apolar desolvation with a corresponding increase of the entropy content of their binding free-energy. These findings explain the reported size-dependence of maximal available affinity and ligand efficiency both behaving differently in the small molecule region featured by strong H-bond formation and in the larger molecule region featured by apolar desolvation.

Classical Polarization in Hybrid QM/MM Methods

We have presented a method for modeling polarization in hybrid QM/MM calculations. The method, which expresses the induced dipoles as a set of "induced" charges, is based on the induced dipole approach and methodology for calculating potential-derived point charges from distributed multipole series. The method has the advantage that the same methodology can be used to determine the induced charges and the potential derived charges and so both sets of charges are rigorously defined within the same framework. This underlying link with the wave function makes the method particularly suitable for use in hybrid QM/MM calculations. Here we assess the importance of explicit polarization in the classical part of a QM/MM system with regard to improving the classical description and the consequent effects on the quantum description. The main advantages of the induced charge approach are that the method is readily interfaced with quantum mechanical methods and that induced charges are more readily interpreted than induced dipoles. The ease of interpretation is illustrated by analysis of the charges involved in dimeric and trimeric hydrogen bonded systems. The method for treating the MM polarization has been validated by a regression analysis of the charges induced in both the QM and MM systems against those derived from full quantum mechanical calculations. The method has also been validated using two energy decomposition approaches, which show that MM polarization makes a significant and reliable contribution to the QM - MM interaction energy in a hybrid system. The distance dependency of the induced charges is investigated in calculations on methylsuccinyl-Ala-Ala-Pro-Ala chlormethyl ketone interacting with human neutrophil elastase and propranolol interacting with asparagine residues in a model of the beta(2)-adrenergic receptor.