Ligand-biased ensemble receptor docking (LigBEnD): a hybrid ligand/receptor structure-based approach

Ligand docking to flexible protein molecules can be efficiently carried out through ensemble docking to multiple protein conformations, either from experimental X-ray structures or from in silico simulations. The success of ensemble docking often requires the careful selection of complementary protein conformations, through docking and scoring of known co-crystallized ligands. False positives, in which a ligand in a wrong pose achieves a better docking score than that of native pose, arise as additional protein conformations are added. In the current study, we developed a new ligand-biased ensemble receptor docking method and composite scoring function which combine the use of ligand-based atomic property field (APF) method with receptor structure-based docking. This method helps us to correctly dock 30 out of 36 ligands presented by the D3R docking challenge. For the six mis-docked ligands, the cognate receptor structures prove to be too different from the 40 available experimental Pocketome conformations used for docking and could be identified only by receptor sampling beyond experimentally explored conformational subspace.     

Heterocyclic Chemistry Applied to the Design of N-Acyl Homoserine Lactone Analogues as Bacterial Quorum Sensing Signals Mimics

N-acyl homoserine lactones (AHLs) are small signaling molecules used by many Gram-negative bacteria for coordinating their behavior as a function of their population density. This process, based on the biosynthesis and the sensing of such molecular signals, and referred to as Quorum Sensing (QS), regulates various gene expressions, including growth, virulence, biofilms formation, and toxin production. Considering the role of QS in bacterial pathogenicity, its modulation appears as a possible complementary approach in antibacterial strategies. Analogues and mimics of AHLs are therefore biologically relevant targets, including several families in which heterocyclic chemistry provides a strategic contribution in the molecular design and the synthetic approach. AHLs consist of three main sections, the homoserine lactone ring, the central amide group, and the side chain, which can vary in length and level of oxygenation. The purpose of this review is to summarize the contribution of heterocyclic chemistry in the design of AHLs analogues, insisting on the way heterocyclic building blocks can serve as replacements of the lactone moiety, as a bioisostere for the amide group, or as an additional pattern appended to the side chain. A few non-AHL-related heterocyclic compounds with AHL-like QS activity are also mentioned.     

Predicting bioactive conformations and binding modes of macrocycles

Macrocyclic compounds experience increasing interest in drug discovery. It is often thought that these large and chemically complex molecules provide promising candidates to address difficult targets and interfere with protein–protein interactions. From a computational viewpoint, these molecules are difficult to treat. For example, flexible docking of macrocyclic compounds is hindered by the limited ability of current docking approaches to optimize conformations of extended ring systems for pose prediction. Herein, we report predictions of bioactive conformations of macrocycles using conformational search and binding modes using docking. Conformational ensembles generated using specialized search technique of about 70 % of the tested macrocycles contained accurate bioactive conformations. However, these conformations were difficult to identify on the basis of conformational energies. Moreover, docking calculations with limited ligand flexibility starting from individual low energy conformations rarely yielded highly accurate binding modes. In about 40 % of the test cases, binding modes were approximated with reasonable accuracy. However, when conformational ensembles were subjected to rigid body docking, an increase in meaningful binding mode predictions to more than 50 % of the test cases was observed. Electrostatic effects did not contribute to these predictions in a positive or negative manner. Rather, achieving shape complementarity at macrocycle-target interfaces was a decisive factor. In summary, a combined computational protocol using pre-computed conformational ensembles of macrocycles as a starting point for docking shows promise in modeling binding modes of macrocyclic compounds.     

Relationship among ligand conformations in solution, in the solid state, and at the Hsp90 binding site: geldanamycin and radicicol

The unknown effects of a receptor's environment on a ligand's conformation presents a difficult challenge in predicting feasible bioactive conformations, particularly if the receptor is ill-defined. The primary hypothesis of this work is that a structure's conformational ensemble in solution presents viable candidates for protein binding. The experimental solution profile can be achieved with the NAMFIS (NMR analysis of molecular flexibility in solution) method, which deconvolutes the average NMR spectrum of small flexible molecules into individual contributing conformations with varying populations. Geldanamycin and radicicol are structurally different macrocycles determined by X-ray crystallography to bind to a common site on the cellular chaperone heat shock protein 90 (Hsp90). Without benefit of a receptor structure, NAMFIS has identified the bioactive conformers of geldanamycin and radicicol in CDCl3 solution with populations of 4% and 21%, respectively. Conversely, docking the set of NAMFIS conformers into the unliganded proteins with GLIDE followed by MM-GBSA scoring reproduces the experimental crystallographic binding poses.     

Response surface methodology in drug design: A case study on docking analysis of a potent antifungal fluconazole

Molecular docking is a valuable in silico technique for discovery/design of bioactive compounds. A current challenge within docking simulations is the incorporation of receptor flexibility. A useful strategy toward solving such problem would be the docking of a typical ligand into the multiple conformations of the target. In this study, a multifactor response surface model was constructed to estimate the AutoDock based binding free energy of fluconazole within multiple conformations of 14α-demethylase (CYP51) (cross docking) as a validated antifungal target. On the basis of developed models, individual and interactive effects of important experimental parameters on cross docking of fluconazole were elucidated. For this purpose, a set of high-resolution holo crystallographic structures from CYP51 of human pathogen Trypanosoma cruzi were retrieved to statistically model the binding mode and affinity of fluconazole. The changes of AutoDock binding free energy for the complexes of CYP51-fluconazole were elucidated with the simultaneous variations of six independent variables including grid size, grid spacing, number of genetic algorithm (GA) runs, maximum number of energy evaluations, torsion degrees for ligand and quaternion degrees for ligand. It was revealed that grid spacing (distance between adjacent grid points) and maximum number of energy evaluations were two significant model terms. It was also revealed that grid size, torsion degrees for ligand and quaternion degrees for ligand had insignificant effects on estimated binding energy while the effect of GA runs was non-significant. The interactive effect of “torsion degrees for ligand” with number of GA runs was found to be the significant factor. Other important interactive effects were the interaction of “number of GA runs” with “grid spacing” and “number of energy evaluations” with “grid size”. Furthermore; results of modeling studies within several CYP51 conformations exhibited that “number of GA runs” and “number of energy evaluations” were less sensitive to varied target conformations.     

GalaxyDock2: Protein–ligand docking using beta‐complex and global optimization

In this article, an enhanced version of GalaxyDock protein–ligand docking program is introduced. GalaxyDock performs conformational space annealing (CSA) global optimization to find the optimal binding pose of a ligand both in the rigid‐receptor mode and the flexible‐receptor mode. Binding pose prediction has been improved compared to the earlier version by the efficient generation of high‐quality initial conformations for CSA using a predocking method based on a beta‐complex derived from the Voronoi diagram of receptor atoms. Binding affinity prediction has also been enhanced by using the optimal combination of energy components, while taking into consideration the energy of the unbound ligand state. The new version has been tested in terms of binding mode prediction, binding affinity prediction, and virtual screening on several benchmark sets, showing improved performance over the previous version and AutoDock, on which the GalaxyDock energy function is based. GalaxyDock2 also performs better than or comparable to other state‐of‐the‐art docking programs. GalaxyDock2 is freely available at http://galaxy.seoklab.org/softwares/galaxydock.html . © 2013 Wiley Periodicals, Inc.     

Expanding dialogues: from natural autoinducers to non-natural analogues that modulate quorum sensing in Gram-negative bacteria

Bacteria are capable of "communicating" their local population densities via a process termed quorum sensing (QS). Gram-negative bacteria use N-acylated l-homoserine lactones (AHLs), in conjunction with their cognate LuxR-type receptors, as their primary signalling circuit for QS. In this critical review, we examine AHL signalling in Gram-negative bacteria with a primary focus on the design of non-natural AHLs, their structure-activity relationships, and their application in chemical biological approaches to study QS (72 references).     

Antibody interference with N-acyl homoserine lactone-mediated bacterial quorum sensing

Many bacterial pathogens coordinate their virulence factor expression in a cell density-dependent manner. This population-dependent coordination of gene expression in bacteria has been termed "quorum sensing" (QS). N-Acyl homoserine lactones (AHLs) are used by over 70 Gram-negative bacterial species as autoinducers. Inhibition of QS signaling might represent a new target for antimicrobial therapy. Here we report the hapten design, synthesis, generation of monoclonal antibodies (mAbs) against AHLs, and the evaluation of these mAbs for their ability to blunt QS signaling and inhibit virulence factor expression in P. aeruginosa. The mAbs can be envisioned as a tool for future investigations into AHL-based QS, which may aid in gaining new insights into the pathogenesis of P. aeruginosa and may ultimately lead to the development of new strategies to combat bacterial diseases.     

Use of solid-phase extraction to enable enhanced detection of acyl homoserine lactones (AHLs) in environmental samples

A challenge for understanding the role of bacterial cell–cell signalling in the environment is the detection of those signals, which are often present in low (nmol L⁻¹) concentrations. We describe here a simple purification method, solid-phase extraction (SPE), for increasing the sensitivity of detection for one such group of signals, acyl homoserine lactones (AHLs), in environmental samples. Spiking of dried marine sponge tissue (Stylinos sp.) with AHLs resulted in detection down to 0.01 ppm for 3-oxo-hexanoyl homoserine lactone (3-oxo C₆-HSL) and 1 ppm for hexanoyl homoserine lactone (C₆-HSL). Compared with liquid extraction methods use of SPE resulted in twofold and tenfold improvements in sensitivity, respectively.     

Efficient inclusion of receptor flexibility in grid-based protein-ligand docking

Accounting for receptor flexibility is an essential component of successful protein-ligand docking but still marks a major computational challenge. For many target molecules of pharmaceutical relevance, global backbone conformational changes are relevant during the ligand binding process. However, popular methods that represent the protein receptor molecule as a potential grid typically assume a rigid receptor structure during ligand-receptor docking. A new approach has been developed that combines inclusion of global receptor flexibility with the efficient potential grid representation of the receptor molecule. This is achieved using interpolation between grid representations of the receptor protein deformed in selected collective degrees of freedom. The method was tested on the docking of three ligands to apo protein kinase A (PKA), an enzyme that undergoes global structural changes upon inhibitor binding. Structural variants of PKA were generated along the softest normal mode of an elastic network representation of apo PKA. Inclusion of receptor deformability during docking resulted in a significantly improved docking performance compared with rigid PKA docking, thus allowing for systematic virtual screening applications at small additional computational cost.     

A robust force field based method for calculating conformational energies of charged drug-like molecules

The binding affinity of a drug-like molecule depends among other things on the availability of the bioactive conformation. If the bioactive conformation has a significantly higher energy than the global minimum energy conformation, then the molecule is unlikely to bind to its target. Determination of the global minimum energy conformation and calculation of conformational penalties of binding is a prerequisite for prediction of reliable binding affinities. Here, we present a simple and computationally efficient procedure to estimate the global energy minimum for a wide variety of structurally diverse molecules, including polar and charged compounds. Identifying global energy minimum conformations of such compounds with force field methods is problematic due to the exaggeration of intramolecular electrostatic interactions. We demonstrate that the global energy minimum conformations of zwitterionic compounds generated by conformational analysis with modified electrostatics are good approximations of the conformational distributions predicted by experimental data and with molecular dynamics performed in explicit solvent. Finally the method is used to calculate conformational penalties for zwitterionic GluA2 agonists and to filter false positives from a docking study.     

蛋白质分子与配体结合过程构象变化的研究

蛋白质分子与配体的作用模式主要有直接的环区结合及铰链式结合两种方式。针对这两种不同的作用方式,我们提出采用不同的策略进行结合过程的构象研究。对于直接的环区结合模式,通过建立环区主链构象库,来实现蛋白质环区与配体的准柔性对接,并以链霉抗生物素蛋白体系为例对构象库建立的可行性进行了验证计算。对铰链结合方式,采用分步对接的方法进行计算,并具体应用于HIV蛋白酶与其小分子配体的结合过程。计算结果表明,这两种处理方法分别能较好地模拟不同类型的蛋白质与配体结合的的构象变化。     

Single drop microextraction of homoserine lactones based quorum sensing signal molecules, and the separation of their enantiomers using gas chromatography mass spectrometry in the presence of biological matrices

N-acylated homoserine lactones (AHLs) are produced by Gram-negative bacteria as communication signals and are frequently studied as mediators of the “quorum sensing” response of bacterial communities. AHLs are optically active components and the L-form is known to have activity in the quorum sensing. However, the knowledge regarding the stereospecific production of the AHLs in bacterial cultures is limited; therefore, there is a need for a fast and easy method for their chiral analysis. A method was developed for the preconcentration of the AHLs using single drop microextraction or liquid-liquid microphase extraction in toluene for their analysis using GC-MS. The performance of the method were determined and discussed for the chiral separation of these autoinducers using a capillary column coated with heptakis-(2,3-di-O-acetyl-6-O-t-butyldimethyl-silyl)-β-cyclodextrin. The salient feature of this study is the demonstration, that Burkholderia cepacia LA3 produced D-decanoyl-homoserine lactone beside L-decanoyl- and L-octaonyl- enantiomers.     

Fully automated molecular mechanics based induced fit protein-ligand docking method

We describe a method for docking a ligand into a protein receptor while allowing flexibility of the protein binding site. The method employs a multistep procedure that begins with the generation of protein and ligand conformations. An initial placement of the ligand is then performed by computing binding site hotspots. This initial placement is followed by a protein side-chain refinement stage that models protein flexibility. The final step of the process is an energy minimization of the ligand pose in the presence of the rigid receptor. Thus the algorithm models flexibility of the protein at two stages, before and after ligand placement. We validated this method by performing docking and cross docking studies of eight protein systems for which crystal structures were available for at least two bound ligands. The resulting rmsd values of the 21 docked protein-ligand complexes showed values of 2 A or less for all but one of the systems examined. The method has two critical benefits for high throughput virtual screening studies. First, no user intervention is required in the docking once the initial binding site selection has been made in the protein. Second, the initial protein conformation generation needs to be performed only once for a given binding region. Also, the method may be customized in various ways depending on the particular scenario in which dockings are being performed. Each of the individual steps of the method is fully independent making it straightforward to explore different variants of the high level workflow to further improve accuracy and performance.     

A fast and efficient method to generate biologically relevant conformations

Summary Mutual binding between a ligand of low molecular weight and its macromolecular receptor demands structural complementarity of both species at the recognition site. To predict binding properties of new molecules before synthesis, information about possible conformations of drug molecules at the active site is required, especially if the 3D structure of the receptor is not known. The statistical analysis of small-molecule crystal data allows one to elucidate conformational preferences of molecular fragments and accordingly to compile libraries of putative ligand conformations. A comparison of geometries adopted by corresponding fragments in ligands bound to proteins shows similar distributions in conformation space. We have developed an automatic procedure that generates different conformers of a given ligand. The entire molecule is decomposed into its individual ring and open-chain torsional fragments, each used in a variety of favorable conformations. The latter ones are produced according to the library information about conformational preferences. During this building process, an extensive energy ranking is applied. Conformers ranked as energetically favorable are subjected to an optimization in torsion angle space. During minimization, unfavorable van der Waals interactions are removed while keeping the open-chain torsion angles as close as possible to the experimentally most frequently observed values. In order to assess how well the generated conformers map conformation space, a comparison with experimental data has been performed. This comparison gives some confidence in the efficiency and completeness of this approach. For some ligands that had been structurally characterized by protein crystallography, the program was used to generate sets of some 10 to 100 conformers. Among these, geometries are found that fall convincingly close to the conformations actually adopted by these ligands at the binding site.     

Identification of bacterial <Emphasis Type="Italic">N</Emphasis>-acylhomoserine lactones (AHLs) with a combination of ultra-performance liquid chromatography (UPLC), ultra-high-resolution mass spectrometry,and in-situ biosensors

N-Acylated homoserine lactones (AHLs) are produced by Gram-negative bacteria as communication signals and are frequently studied as mediators of the “quorum sensing” response of bacterial communities. Several reports have recently been published on the identification of AHLs from different species and attempts have been made to study their role in natural habitats, for example the surface of plant roots in the rhizosphere. In this article, different analytical methods, including bacterial biosensors and chromatographic techniques, are reviewed. A concept for assignment of the structures of AHLs is also presented. The retention behaviour of derivatives of AHLs containing β-keto or hydroxyl groups and/or double bonds has been evaluated in relation to the separation behaviour of AHLs with saturated and unsubstituted alkanoyl chains. Samples have also been analysed by high resolution mass spectrometry (Fourier-transform ion-cyclotron-resonance mass spectrometry, FTICR-MS), nano liquid chromatography–electrospray ionization ion trap mass spectrometry (nano-LC–MS) and by the aid of a biosensor. The results obtained from ultra performance liquid chromatography (UPLC), FTICR-MS, nano-LC–MS, and bioassays have been compared to attempt structural characterisation of AHL without chemical synthesis of analytical standards. The method was used to identify the major AHL compound produced by the rhizosphere bacterium Acidovorax sp. N35 as N-(3-hydroxydecanoyl)homoserine lactone.     

Simultaneous quantitative profiling of N-acyl-l-homoserine lactone and 2-alkyl-4(1H)-quinolone families of quorum-sensing signaling molecules using LC-MS/MS

An LC-MS/MS method, using positive mode electrospray ionization, for the simultaneous, quantitative and targeted profiling of the N-acyl-L-homoserine lactone (AHL) and 2-alkyl 4-(1H)-quinolone (AQ) families of bacterial quorum-sensing signaling molecules (QSSMs) is presented. This LC-MS/MS technique was applied to determine the relative molar ratios of AHLs and AQs produced by Pseudomonas aeruginosa and the consequences of mutating individual or multiple QSSM synthase genes (lasI, rhlI, pqsA) on AHL and AQ profiles and concentrations. The AHL profile of P. aeruginosa was dominated by N-butanoyl-L-homoserine lactone (C4-HSL) with lesser concentrations of N-hexanoyl-L-homoserine lactone (C6-HSL) and 3-oxo-substituted longer chain AHLs including N-(3-oxodecanoyl)-L-homoserine lactone (3-oxo-C10-HSL) and N-(3-oxododecanoyl)-L-homoserine lactone (3-oxo-C12-HSL). The AQ profile of P. aeruginosa comprised the C7 and C9 long alkyl chain AQs including 2-heptyl-4-hydroxyquinoline (HHQ), 2-nonyl-4-hydroxyquinoline, the "pseudomonas quinolone signal" (2-heptyl-3-hydroxy-4-quinolone) and the N-oxides, 2-heptyl-4-hydroxyquinoline N-oxide and 2-nonyl-4-hydroxyquinoline N-oxide. Application of the method showed significant effects of growth medium type on the ratio and the nature of the QSSMs synthesized and the dramatic effect of single, double and triple mutations in the P. aeruginosa QS synthase genes. The LC-MS/MS methodology is applicable in organisms where either or both AHL and AQ QSSMs are produced and can provide comprehensive profiles and concentrations from a single sample.     

Predicting binding energies of CDK6 inhibitors in the hit-to-lead process

The main challenge for the ??hit-to-lead?? stage in the drug discovery process relies on the accuracy of existing docking methods. In fact, accuracy of docking methods depends not only on the scoring function used to rank the poses but also on the ability of the docking method to reproduce the experimental binding mode. At this purpose, the performance of different approximations to properly dock and score compounds with known activity in a narrow range of IC₅₀ values was analyzed. A set of five ATP-competitive CDK6 inhibitors and three receptor conformations for CDK6 were considered for analysis, and three methodologies were used and analyzed in order to include different degrees of receptor flexibility. Thus, a completely rigid receptor is considered when using Glide, while the so-called Induced Fit Docking Protocol accounts for receptor sidechain rearrangements. Finally, force field calculations were also performed in order to consider a completely flexible receptor.