Connectivity and binding‐site recognition: Applications relevant to drug design

Here, we describe a family of methods based on residue–residue connectivity for characterizing binding sites and apply variants of the method to various types of protein–ligand complexes including proteases, allosteric‐binding sites, correctly and incorrectly docked poses, and inhibitors of protein–protein interactions. Residues within ligand‐binding sites have about 25% more contact neighbors than surface residues in general; high‐connectivity residues are found in contact with the ligand in 84% of all complexes studied. In addition, a k‐means algorithm was developed that may be useful for identifying potential binding sites with no obvious geometric or connectivity features. The analysis was primarily carried out on 61 protein–ligand structures from the MEROPS protease database, 250 protein–ligand structures from the PDBSelect (25%), and 30 protein–protein complexes. Analysis of four proteases with crystal structures for multiple bound ligands has shown that residues with high connectivity tend to have less variable side‐chain conformation. The relevance to drug design is discussed in terms of identifying allosteric‐binding sites, distinguishing between alternative docked poses and designing protein interface inhibitors. Taken together, this data indicate that residue–residue connectivity is highly relevant to medicinal chemistry. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Multiple ligand simultaneous docking: Orchestrated dancing of ligands in binding sites of protein

Present docking methodologies simulate only one single ligand at a time during docking process. In reality, the molecular recognition process always involves multiple molecular species. Typical protein–ligand interactions are, for example, substrate and cofactor in catalytic cycle; metal ion coordination together with ligand(s); and ligand binding with water molecules. To simulate the real molecular binding processes, we propose a novel multiple ligand simultaneous docking (MLSD) strategy, which can deal with all the above processes, vastly improving docking sampling and binding free energy scoring. The work also compares two search strategies: Lamarckian genetic algorithm and particle swarm optimization, which have respective advantages depending on the specific systems. The methodology proves robust through systematic testing against several diverse model systems: E. coli purine nucleoside phosphorylase (PNP) complex with two substrates, SHP2NSH2 complex with two peptides and Bcl‐xL complex with ABT‐737 fragments. In all cases, the final correct docking poses and relative binding free energies were obtained. In PNP case, the simulations also capture the binding intermediates and reveal the binding dynamics during the recognition processes, which are consistent with the proposed enzymatic mechanism. In the other two cases, conventional single‐ligand docking fails due to energetic and dynamic coupling among ligands, whereas MLSD results in the correct binding modes. These three cases also represent potential applications in the areas of exploring enzymatic mechanism, interpreting noisy X‐ray crystallographic maps, and aiding fragment‐based drug design, respectively. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Two‐dimensional replica‐exchange method for predicting protein–ligand binding structures

We have developed a two‐dimensional replica‐exchange method for the prediction of protein–ligand binding structures. The first dimension is the umbrella sampling along the reaction coordinate, which is the distance between a protein binding pocket and a ligand. The second dimension is the solute tempering, in which the interaction between a ligand and a protein and water is weakened. The second dimension is introduced to make a ligand follow the umbrella potential more easily and enhance the binding events, which should improve the sampling efficiency. As test cases, we applied our method to two protein‐ligand complex systems (MDM2 and HSP 90‐alpha). Starting from the configuration in which the protein and the ligand are far away from each other in each system, our method predicted the ligand binding structures in excellent agreement with the experimental data from Protein Data Bank much faster with the improved sampling efficiency than the replica‐exchange umbrella sampling method that we have previously developed. © 2013 Wiley Periodicals, Inc.     

Affinity Screening Using Competitive Binding with Fluorine‐19 Hyperpolarized Ligands

Fluorine‐19 NMR and hyperpolarization form a powerful combination for drug screening. Under a competitive equilibrium with a selected fluorinated reporter ligand, the dissociation constant (K_D) of other ligands of interest is measurable using a single‐scan Carr–Purcell–Meiboom–Gill (CPMG) experiment, without the need for a titration. This method is demonstrated by characterizing the binding of three ligands with different affinities for the serine protease trypsin. Monte Carlo simulations show that the highest accuracy is obtained when about one‐half of the bound reporter ligand is displaced in the binding competition. Such conditions can be achieved over a wide range of affinities, allowing for rapid screening of non‐fluorinated compounds when a single fluorinated ligand for the binding pocket of interest is known.     

Binding constant determination of high-affinity protein-ligand complexes by electrospray ionization mass spectrometry and ligand competition

We describe an approach for the determination of binding constants for protein-ligand complexes with electrospray ionization mass spectrometry, based on the observation of unbound ligands competing for binding to a protein target. For the first time, dissociation constants lower than picomolar could be determined with good accuracy by electrospray ionization mass spectrometry. The presented methodology relies only on the determination of signal intensity ratios for free ligands in the low mass region. Therefore, all the advantages of measuring low masses with mass spectrometry, such as high resolution are preserved. By using a reference ligand with known binding affinity, the affinity of a second ligand can be determined. Since no noncovalently bound species are observed, assumptions about response factors are not necessary. The method is validated with ligands binding to avidin and applied to ligands binding to p38 mitogen-activated protein kinase.     

Electronic Olfactory Sensor Based on A. mellifera Odorant‐Binding Protein 14 on a Reduced Graphene Oxide Field‐Effect Transistor

An olfactory biosensor based on a reduced graphene oxide (rGO) field‐effect transistor (FET), functionalized by the odorant‐binding protein 14 (OBP14) from the honey bee (Apis mellifera) has been designed for the in situ and real‐time monitoring of a broad spectrum of odorants in aqueous solutions known to be attractants for bees. The electrical measurements of the binding of all tested odorants are shown to follow the Langmuir model for ligand–receptor interactions. The results demonstrate that OBP14 is able to bind odorants even after immobilization on rGO and can discriminate between ligands binding within a range of dissociation constants from K_d=4 μM to K_d=3.3 mM . The strongest ligands, such as homovanillic acid, eugenol, and methyl vanillate all contain a hydroxy group which is apparently important for the strong interaction with the protein.     

Evaluation of the searching abilities of HBOP and HBSITE for binding pocket detection

We evaluated the pocket‐searching abilities of the computer programs HBOP and HBSITE, in which hydrophobic potentials calculated on grid points generated around a protein function as an indicator of the pocket‐like property‐using a test set of 458 experimentally observed structures of protein–ligand complexes. The results were compared with those obtained using the alternative algorithms PASS and SiteID, which only consider geometric properties, and Q‐SiteFinder, which considers not only geometry but also physicochemical properties. The comparison revealed that HBOP and HBSITE could detect experimentally observed ligand‐binding pockets for more test complexes than PASS and SiteID. In addition, the success rate of HBOP for detecting binding pockets was higher than that of Q‐SiteFinder, and that of HBSITE was comparable with that of Q‐SiteFinder. Results of tests for ligand‐unbound state proteins indicated that HBSITE was more appropriate than Q‐SiteFinder for pocket searches of ligand‐unbound proteins, and HBSITE was more robust than Q‐SiteFinder against structural differences between ligand‐bound and ‐unbound state proteins. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

Synthesis,spectroscopy, thermal and biological aspect of novel six‐coordinated dimeric iron(III) mixed‐ligand complexes

The mixed‐ligand complexes of iron(III) with 1‐cyclopropyl‐6‐fluoro‐4‐oxo‐7‐piperazin‐1‐yl‐1,4‐dihydroquinoline‐3‐carboxylic acid and various neutral bidentate Schiff base ligands were prepared. The structure of mixed‐ligand complexes was investigated using spectral, physicochemical and elemental analyses. Biocidal activity was determined using agar plate technique against Staphylococcus aureus, Bacillus subtilis, Bacillus cereus, Salmonella typhi, Escherichia coli and Serratia marcescens . The result showed a significant increase in a biocidal activity compared with parent ligands, metal salts and standard drugs (ofloxacin, levofloxacin). DNA binding and cleavage studies were carried out using absorption titration and gel electrophoresis techniques, respectively. The binding constant of Fe(III) complexes was obtained in the range 2.5–4.0 × 10⁴ M ^?1. The DNA binding and cleavage efficacy were raised in mixed‐ligand complexes as compared with parental ligands and metal salts. Copyright © 2008 John Wiley & Sons, Ltd.     

Converging ligand‐binding free energies obtained with free‐energy perturbations at the quantum mechanical level

In this article, the convergence of quantum mechanical (QM) free‐energy simulations based on molecular dynamics simulations at the molecular mechanics (MM) level has been investigated. We have estimated relative free energies for the binding of nine cyclic carboxylate ligands to the octa‐acid deep‐cavity host, including the host, the ligand, and all water molecules within 4.5 Å of the ligand in the QM calculations (158–224 atoms). We use single‐step exponential averaging (ssEA) and the non‐Boltzmann Bennett acceptance ratio (NBB) methods to estimate QM/MM free energy with the semi‐empirical PM6‐DH2X method, both based on interaction energies. We show that ssEA with cumulant expansion gives a better convergence and uses half as many QM calculations as NBB, although the two methods give consistent results. With 720,000 QM calculations per transformation, QM/MM free‐energy estimates with a precision of 1 kJ/mol can be obtained for all eight relative energies with ssEA, showing that this approach can be used to calculate converged QM/MM binding free energies for realistic systems and large QM partitions. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

STOCK: Structure mapper and online coarse‐graining kit for molecular simulations

We present a web toolkit STructure mapper and Online Coarse‐graining Kit for setting up coarse‐grained molecular simulations. The kit consists of two tools: structure mapping and Boltzmann inversion tools. The aim of the first tool is to define a molecular mapping from high, for example, all‐atom, to low, that is, coarse‐grained, resolution. Using a graphical user interface it generates input files, which are compatible with standard coarse‐graining packages, for example, Versatile Object‐oriented Toolkit for Coarse‐graining Applications and DL_CGMAP. Our second tool generates effective potentials for coarse‐grained simulations preserving the structural properties, for example, radial distribution functions, of the underlying higher resolution model. The required distribution functions can be provided by any simulation package. Simulations are performed on a local machine and only the distributions are uploaded to the server. The applicability of the toolkit is validated by mapping atomistic pentane and polyalanine molecules to a coarse‐grained representation. Effective potentials are derived for systems of TIP3P (transferable intermolecular potential 3 point) water molecules and salt solution. The presented coarse‐graining web toolkit is available at http://stock.cmm.ki.si . © 2014 Wiley Periodicals, Inc.     

Semiempirical calculations of binding enthalpy for protein–ligand complexes

A systematic semiempirical quantum mechanical study of the interactions between proteins and ligands has been performed to determine the ability of this approach for the accurate estimation of the enthalpic contribution to the binding free energy of the protein–ligand systems. This approach has been applied for eight test protein–ligand complexes with experimentally known binding enthalpies. The calculations were performed using the semiempirical PM3 approach incorporated in the MOPAC 97, ZAVA originally elaborated in Algodign, and MOPAC 2002 with MOZYME facility packages. Special attention was paid to take into account structural water molecules, which were located in the protein–ligand binding site. It was shown that the results of binding enthalpy calculations fit experimental data within ～2 kcal/mol in the presented approach. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Fast,metadynamics‐based method for prediction of the stereochemistry‐dependent relative free energies of ligand–receptor interactions

The computational approach applicable for the molecular dynamics (MD)‐based techniques is proposed to predict the ligand–protein binding affinities dependent on the ligand stereochemistry. All possible stereoconfigurations are expressed in terms of one set of force‐field parameters [stereoconfiguration‐independent potential (SIP)], which allows for calculating all relative free energies by only single simulation. SIP can be used for studying diverse, stereoconfiguration‐dependent phenomena by means of various computational techniques of enhanced sampling. The method has been successfully tested on the β₂‐adrenergic receptor (β₂‐AR) binding the four fenoterol stereoisomers by both metadynamics simulations and replica‐exchange MD. Both the methods gave very similar results, fully confirming the presence of stereoselective effects in the fenoterol‐β₂‐AR interactions. However, the metadynamics‐based approach offered much better efficiency of sampling which allows for significant reduction of the unphysical region in SIP. © 2014 Wiley Periodicals, Inc.     

Theoretical studies on the structural rearrangement of ligand binding pocket in human vitamin D receptor

The structural rearrangement of the ligand binding domain (LBD) of human Vitamin D receptor (hVDR) complexed with 1α, 25‐dihydroxyvitamin D₃ (natural ligand) and its analogues (denoted as b and c ) was studied by molecular dynamics (MD) simulations. MD simulations revealed that these ligands could induce different structural changes of LBD, in which 1α, 25‐dihydroxyvitamin D₃ only led to a minute change, suggesting that LBD adopted its canonical active conformation upon binding the natural ligand, while b and c could provoke a clear structural rearrangement of the LBD. In complex of hVDR‐LBD/ b , it is found that helix 6 (H6) and subsequent loop 6‐7 shift outward and the last turn of H11 shifts away from H12, which generate a new cavity at the bottom of binding pocket to accommodate the extra butyl group on the side chain of ligand b . As for hVDR‐LBD/ c , the steric exclusion of the second side chain of ligand c makes the N‐terminal of H7 move outsides and C‐terminal of H11 close to H12, expanding the bottom of the pocket. These calculation results agree well the experimental observations. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Efficient estimation of binding free energies between peptides and an MHC class II molecule using coarse‐grained molecular dynamics simulations with a weighted histogram analysis method

We estimate the binding free energy between peptides and an MHC class II molecule using molecular dynamics (MD) simulations with the weighted histogram analysis method (WHAM). We show that, owing to its more thorough sampling in the available computational time, the binding free energy obtained by pulling the whole peptide using a coarse‐grained (CG) force field (MARTINI) is less prone to significant error induced by inadequate‐sampling than using an atomistic force field (AMBER). We further demonstrate that using CG MD to pull 3–4 residue peptide segments while leaving the remaining peptide segments in the binding groove and adding up the binding free energies of all peptide segments gives robust binding free energy estimations, which are in good agreement with the experimentally measured binding affinities for the peptide sequences studied. Our approach thus provides a promising and computationally efficient way to rapidly and reliably estimate the binding free energy between an arbitrary peptide and an MHC class II molecule. © 2017 Wiley Periodicals, Inc.     

PLIMSTEX: a novel mass spectrometric method for the quantification of protein–ligand interactions in solution

Protein–ligand interactions by mass spectrometry, titration, and H/D exchange (PLIMSTEX) is a new mass spectrometric method for determining association constants and binding stoichiometry for interactions of proteins with various ligands, as well as for quantifying the conformational changes associated with ligand binding to proteins. The association constants determined with PLIMSTEX agree with literature values within a factor of six, establishing its validity for protein interactions involving metal ions, small organic molecules, peptides, and proteins. PLIMSTEX provides solution, not gas-phase, properties by taking advantage of ESI and MALDI mass spectrometry to measure accurately the mass of a protein as it undergoes amide H/D exchange. The approach sidesteps the problem of relating gas-phase abundances of the protein or protein–ligand complex ions to their solution concentrations. With on-column concentration and desalting, high picomole quantities of proteins are sufficient for reproducible mass detection, and the concentration of the protein can be as low as 10⁻⁸ M. It is amenable to different protein/ligand systems in physiologically relevant media. No specially labeled protein or ligand is needed. PLIMSTEX offers minimal perturbation of the binding equilibrium because it uses no denaturants, no additional spectroscopy or reaction probes, and no physical separation of ligand and protein during binding.     

How to obtain statistically converged MM/GBSA results

The molecular mechanics/generalized Born surface area (MM/GBSA) method has been investigated with the aim of achieving a statistical precision of 1 kJ/mol for the results. We studied the binding of seven biotin analogues to avidin, taking advantage of the fact that the protein is a tetramer with four independent binding sites, which should give the same estimated binding affinities. We show that it is not enough to use a single long simulation (10 ns), because the standard error of such a calculation underestimates the difference between the four binding sites. Instead, it is better to run several independent simulations and average the results. With such an approach, we obtain the same results for the four binding sites, and any desired precision can be obtained by running a proper number of simulations. We discuss how the simulations should be performed to optimize the use of computer time. The correlation time between the MM/GBSA energies is ～5 ps and an equilibration time of 100 ps is needed. For MM/GBSA, we recommend a sampling time of 20–200 ps for each separate simulation, depending on the protein. With 200 ps production time, 5–50 separate simulations are required to reach a statistical precision of 1 kJ/mol (800–8000 energy calculations or 1.5–15 ns total simulation time per ligand) for the seven avidin ligands. This is an order of magnitude more than what is normally used, but such a number of simulations is needed to obtain statistically valid results for the MM/GBSA method. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Water molecules inside protein structure affect binding of monosaccharides with HIV‐1 antibody 2G12

Water molecules inside biomolecules constitute integral parts of their structure and participate in the functions of the proteins. Some of the X‐ray crystallographic data are insufficient for analyzing a series of ligand–protein complexes in the same condition. We theoretically investigated antibody binding abilities of saccharide ligands and the effects of the inner water molecules of ligand–antibody complexes. Classical molecular dynamics and quantum chemical simulations using a model with possible water molecules inside the protein were performed with saccharide ligands and Human Immunodeficiency Virus 1 neutralizing antibody 2G12 complexes to estimate how inner water molecules of the protein affect the dynamics of the complexes as well as the ligand–antibody interaction. Our results indicate the fact that d ‐fructose's strong affinity to the antibody was partly due to the good retentiveness of solvent water molecules of the ligand and its stability of the ligand's conformation and relative position in the active site. © 2016 Wiley Periodicals, Inc.     

Covalent Lipid Pocket Ligands Targeting p38α MAPK Mutants

A chemical genetic approach is presented to covalently target a unique lipid binding pocket in the protein kinase p38α, whose function is not yet known. Based on a series of cocrystal structures, a library of 2‐arylquinazolines that were decorated with electrophiles were designed and synthesized to covalently target tailored p38α mutants containing artificially introduced cysteine residues. Matching protein–ligand pairs were identified by MS analysis and further validated by MS/MS studies and protein crystallography. The covalent ligands that emerged from this approach showed excellent selectivity towards a single p38α mutant and will be applicable as suitable probes in future studies of biological systems to dissect the function of the lipid pocket by means of pharmacological perturbations.     

Unraveling Complex Small‐Molecule Binding Mechanisms by Using Simple NMR Spectroscopy

Heteronuclear NMR spectroscopy provides a unique way to obtain site‐specific information about protein–ligand interactions. Usually, such studies rely on the availability of isotopically labeled proteins, thereby allowing both editing of the spectra and ligand signals to be filtered out. Herein, we report that the use of the methyl SOFAST correlation experiment enables the determination of site‐specific equilibrium binding constants by using unlabeled proteins. By using the binding of L ‐ and D ‐tryptophan to serum albumin as a test case, we determined very accurate dissociation constants for both the high‐ and low‐affinity sites present at the protein surface. The values of site‐specific dissociation constants were closer to those obtained by isothermal titration calorimetry than those obtained from ligand‐observed methods, such as saturation transfer difference. The possibility of measuring ligand binding to serum albumin at physiological concentrations with unlabeled proteins may open up new perspectives in the field of drug discovery.     

Exploring Weak Ligand–Protein Interactions by Long‐Lived NMR States: Improved Contrast in Fragment‐Based Drug Screening

Ligands that have an affinity for protein targets can be screened very effectively by exploiting favorable properties of long‐lived states (LLS) in NMR spectroscopy. In this work, we describe the use of LLS for competitive binding experiments to measure accurate dissociation constants of fragments that bind weakly to the ATP binding site of the N‐terminal ATPase domain of heat shock protein 90 (Hsp90), a therapeutic target for cancer treatment. The LLS approach allows one to characterize ligands with an exceptionally wide range of affinities, since it can be used for ligand concentrations [L] that are several orders of magnitude smaller than the dissociation constants K_D. This property makes the LLS method particularly attractive for the initial steps of fragment‐based drug screening, where small molecular fragments that bind weakly to a target protein must be identified, which is a difficult task for many other biophysical methods.