Hydration in drug design. 2. Influence of local site surface shape on water binding

Summary If water molecules are strongly bound at a protein-ligand interface, they are unlikely to be displaced during ligand binding. Such water molecules can change the shape of the ligand binding site and thus affect strategies for drug design. To understand the nature of water binding, and factors influencing it, water molecules at the ligand binding sites of 26 high-resolution protein-ligand complexes have been examined here. Water molecules bound in deep grooves and cavities between the protein and the ligand are located in the indentations on the protein-site surface, but not in the indentations on the ligand surface. The majority of the water molecules bound in deep indentations on the protein-site surface make multiple polar contacts with the protein surface. This may indicate a strong binding of water molecules in deep indentations on protein-site surfaces. The local shape of the site surface may influence the binding of water molecules that mediate protein-ligand interactions.     

MC3T3-E1 cell adhesion to hydroxyapatite with adsorbed bone sialoprotein, bone osteopontin, and bovine serum albumin

Native bone tissue is composed of a complex matrix of collagen, non-collagenous proteins, and hydroxyapatite (HAP). Bone sialoprotein (BSP) and bone osteopontin (OPN) are members of the non-collagenous protein family termed the SIBLING (small integrin-binding ligand, N-linked glycoproteins) proteins, which are primarily found in mineralized tissues. Previously, OPN was shown to exhibit a preferential orientation for MC3T3-E1 cell adhesion when it was specifically bound to collagen, while the MC3T3-E1 cell adhesion was shown to be dependant on the conformational flexibility of BSP specifically bound to collagen. Additionally, OPN was shown to play a greater role than BSP for cell binding to collagen. In this work, the orientations and conformations of BSP and OPN specifically bound to HAP are probed under similar conditions. Radiolabeled adsorption isotherms were obtained for BSP and OPN on HAP formed from a simulated body fluid, and the results show that HAP has the capacity to bind significantly more BSP than OPN. An in vitro MC3T3-E1 cell adhesion assay was then performed to compare the cell binding ability of adsorbed BSP and OPN specifically bound to HAP. It was found that there is a preference for cell binding to HAP with adsorbed BSP as compared to OPN, but not to a statistically significant level. However, the maximum cell binding was observed on HAP substrates with adsorbed heat denatured bovine serum albumin (BSA). The influence of BSA on cell binding was shown to be concentration dependant and it is believed that the adsorbed BSA modulates the proliferation state of the bound cells.     

Refinement of the conformation of UDP-galactose bound to galactosyltransferase using the STD NMR intensity-restrained CORCEMA optimization

The STD NMR technique has originally been described as a tool for screening large compound libraries to identify the lead compounds that are specific to target proteins of interest. The application of this technique in the qualitative epitope mapping of ligands weakly binding to proteins, virus capsid shells, and nucleic acids has also been described. Here we describe the application of the STD NMR intensity-restrained CORCEMA optimization (SICO) procedure for refining the bound conformation of UDP-galactose in galactosyltransferase complex using STD NMR intensities recorded at 500 MHz as the experimental constraints. A comparison of the SICO structure for the bound UDP-galactose in solution with that in the crystal structure for this complex shows some differences in ligand torsion angles and V253 side-chain orientation in the protein. This work describes the first application of an STD NMR intensity-restrained CORCEMA optimization procedure for refining the torsion angles of a bound ligand structure. This method is likely to be useful in structure-based drug design programs since most initial lead compounds generally exhibit weak affinity (millimolar to micromolar) to target proteins of pharmaceutical interest, and the bound conformation of these lead compounds in the protein binding pocket can be determined by the CORCEMA-ST refinement.     

Method for computing protein binding affinity

A Monte Carlo method is given to compute the binding affinity of a ligand to a protein. The method involves extending configuration space by a discrete variable indicating whether the ligand is bound to the protein and a special Monte Carlo move, which allows transitions between the unbound and bound states. Provided that an accurate protein structure is given, that the protein-ligand binding site is known, and that an accurate chemical force field together with a continuum solvation model is used, this method provides a quantitative estimate of the free energy of binding.     

Quantifying labile protein—Ligand interactions using electrospray ionization mass spectrometry

A new electrospray ionization mass spectrometry (ES-MS) approach for quantifying protein—ligand complexes that are prone to in-source (gas-phase) dissociation is described. The method, referred to here as the reference ligand ES-MS method, is based on the direct ES-MS assay and competitive ligand binding. A reference ligand (L_ref), which binds specifically to the protein (P), at the same binding site as the ligand (L) of interest, with known affinity and forms a stable protein—ligand complex in the gas phase, is added to the solution. The fraction of P bound to L_ref, which is determined directly from the ES mass spectrum, is sensitive to the fraction of P bound to L in solution and enables the affinity of P for L to be determined. A mathematical framework for the implementation of the method in cases where P has one or two specific ligand binding sites is given. Affinities of two carbohydrate-binding proteins, a single chain fragment of a monoclonal antibody and the lectin concanavalin A, for monosaccharide ligands are reported and the results are shown to agree with values obtained using isothermal titration calorimetry.     

Structure determination of protein-ligand complexes by transferred paramagnetic shifts

Rational drug design depends on the knowledge of the three-dimensional (3D) structure of complexes between proteins and lead compounds of low molecular weight. A novel nuclear magnetic resonance (NMR) spectroscopy strategy based on the paramagnetic effects from lanthanide ions allows the rapid determination of the 3D structure of a small ligand molecule bound to its protein target in solution and, simultaneously, its location and orientation with respect to the protein. The method relies on the presence of a lanthanide ion in the protein target and on fast exchange between bound and free ligand. The binding affinity of the ligand and the paramagnetic effects experienced in the bound state are derived from concentration-dependent (1)H and (13)C spectra of the ligand at natural isotopic abundance. Combined with prior knowledge of the crystal or solution structure of the protein and of the magnetic susceptibility tensor of the lanthanide ion, the paramagnetic data define the location and orientation of the bound ligand molecule with respect to the protein from simple 1D NMR spectra. The method was verified with the ternary 30 kDa complex between the lanthanide-labeled N-terminal domain of the epsilon exonuclease subunit from the Escherichia coli DNA polymerase III, the subunit theta, and thymidine. The binding mode of thymidine was found to be very similar to that of thymidine monophosphate present in the crystal structure.     

Structural chemistry of a green fluorescent protein Zn biosensor

We designed a green fluorescent protein mutant (BFPms1) that preferentially binds Zn(II) (enhancing fluorescence intensity) and Cu(II) (quenching fluorescence) directly to a chromophore ligand that resembles a dipyrrole unit of a porphyrin. Crystallographic structure determination of apo, Zn(II)-bound, and Cu(II)-bound BFPms1 to better than 1.5 A resolution allowed us to refine metal centers without geometric restraints, to calculate experimental standard uncertainty errors for bond lengths and angles, and to model thermal displacement parameters anisotropically. The BFPms1 Zn(II) site (KD = 50 muM) displays distorted trigonal bipyrimidal geometry, with Zn(II) binding to Glu222, to a water molecule, and tridentate to the chromophore ligand. In contrast, the BFPms1 Cu(II) site (KD = 24 muM) exhibits square planar geometry similar to metalated porphyrins, with Cu(II) binding to the chromophore chelate and Glu222. The apo structure reveals a large electropositive region near the designed metal insertion channel, suggesting a basis for the measured metal cation binding kinetics. The preorganized tridentate ligand is accommodated in both coordination geometries by a 0.4 A difference between the Zn and Cu positions and by distinct rearrangements of Glu222. The highly accurate metal ligand bond lengths reveal different protonation states for the same oxygen bound to Zn vs Cu, with implications for the observed metal ion specificity. Crystallographic anisotropic thermal factor analysis validates metal ion rigidification of the chromophore in enhancement of fluorescence intensity upon Zn(II) binding. Thus, our high-resolution structures reveal how structure-based design has effectively linked selective metal binding to changes in fluorescent properties. Furthermore, this protein Zn(II) biosensor provides a prototype suitable for further optimization by directed evolution to generate metalloprotein variants with desirable physical or biochemical properties.     

AFM imaging of ligand binding to platelet integrin alphaIIbbeta3 receptors reconstituted into planar lipid bilayers

The platelet integrin alphaIIbbeta3 plays a key role in platelet adhesion, activation, and aggregation at the subendothelium and at protein-coated synthetic biomaterials. In this study, interactions between alphaIIbbeta3 and both protein and peptide ligands for the receptor were imaged under physiological conditions by high-resolution atomic force microscopy (AFM). To directly image the ligand-receptor interactions, alphaIIbbeta3 receptors were reconstituted into a supported lipid bilayer formed on a mica surface in the AFM fluid cell assembly and subsequently activated with Mn2+. Fibrinogen, the natural protein ligand for the integrin, as well as a nanogold-labeled peptide ligand (an RGD-containing heptamer) were infused into the AFM fluid cell, incubated with the reconstituted and activated receptors, and imaged under buffer. Height images illustrating topographical features showed the integrin reconstituted in the bilayer. Fibrinogen molecules binding to the receptors were easily observed in the height images, with fibrinogen showing its characteristic trinodular structure and occasionally bridging integrin receptors. Fibrinogen was observed to bind to integrins at the D-domain consistent with the location of the gamma-chain dodecapeptide, while fibrinogen bridging integrins bound to receptors on opposite sides of the protein consistent with a 2-fold axis of symmetry. Peptide ligands were not visible in height images; however, phase images that map the mechanical properties detected the nanogold labels and demonstrated the presence of peptide ligands bound to the receptors. The results demonstrate the ability of this high-resolution microscopy technique to directly visualize single ligand/receptor interactions in a dynamic and physiologically relevant environment, and establish a framework for future fundamental studies of single protein/receptor interactions during normal pathological processes as well as biomaterial surface-induced thrombosis.     

Site‐Specific Investigation of the Steady‐State Kinetics and Dynamics of the Multistep Binding of Bile Acid Molecules to a Lipid Carrier Protein

The investigation of multi‐site ligand–protein binding and multi‐step mechanisms is highly demanding. In this work, advanced NMR methodologies such as 2D ¹H–¹⁵N line‐shape analysis, which allows a reliable investigation of ligand binding occurring on micro‐ to millisecond timescales, have been extended to model a two‐step binding mechanism. The molecular recognition and complex uptake mechanism of two bile salt molecules by lipid carriers is an interesting example that shows that protein dynamics has the potential to modulate the macromolecule–ligand encounter. Kinetic analysis supports a conformational selection model as the initial recognition process in which the dynamics observed in the apo form is essential for ligand uptake, leading to conformations with improved access to the binding cavity. Subsequent multi‐step events could be modelled, for several residues, with a two‐step binding mechanism. The protein in the ligand‐bound state still exhibits a conformational rearrangement that occurs on a very slow timescale, as observed for other proteins of the family. A global mechanism suggesting how bile acids access the macromolecular cavity is thus proposed.     

De novo ligand design with explicit water molecules: an application to bacterial neuraminidase

Most computer-aided drug design methods ignore the presence of crystallographically-determined water molecules in the binding site of a target protein. In this paper, our de novo ligand design methods are applied to the X-ray crystal structure of bacterial neuraminidase in the presence of some selected water molecules. We have found that, for this particular protein, the complete removal of all bound water molecules leads to difficulties in generating any potential ligands if the unsatisfied hydrogen-bonding sitepoints left by removing these water molecules are to be satisfied by a ligand. As more of the crystallographically determined water molecules are allowed in the binding site, it becomes much easier to generate ligands in larger numbers and with wider chemical diversity. This example shows that, in some cases, bound water molecules can be more accessible for hydrogen bonding to an incoming ligand than the actual protein binding sitepoints associated with them. From the point of view of de novo ligand design, water molecules can thus act as versatile amphiprotic hydrogen-bonding sitepoints and reduce the conformational constraints of a particular binding site.     

A carrier protein strategy yields the structure of dalbavancin

Many large natural product antibiotics act by specifically binding and sequestering target molecules found on bacterial cells. We have developed a new strategy to expedite the structural analysis of such antibiotic-target complexes, in which we covalently link the target molecules to carrier proteins, and then crystallize the entire carrier-target-antibiotic complex. Using native chemical ligation, we have linked the Lys-D-Ala-D-Ala binding epitope for glycopeptide antibiotics to three different carrier proteins. We show that recognition of this peptide by multiple antibiotics is not compromised by the presence of the carrier protein partner, and use this approach to determine the first-ever crystal structure for the new therapeutic dalbavancin. We also report the first crystal structure of an asymmetric ristocetin antibiotic dimer, as well as the structure of vancomycin bound to a carrier-target fusion. The dalbavancin structure reveals an antibiotic molecule that has closed around its binding partner; it also suggests mechanisms by which the drug can enhance its half-life by binding to serum proteins, and be targeted to bacterial membranes. Notably, the carrier protein approach is not limited to peptide ligands such as Lys-D-Ala-D-Ala, but is applicable to a diverse range of targets. This strategy is likely to yield structural insights that accelerate new therapeutic development.     

A water-swap reaction coordinate for the calculation of absolute protein-ligand binding free energies

The accurate prediction of absolute protein-ligand binding free energies is one of the grand challenge problems of computational science. Binding free energy measures the strength of binding between a ligand and a protein, and an algorithm that would allow its accurate prediction would be a powerful tool for rational drug design. Here we present the development of a new method that allows for the absolute binding free energy of a protein-ligand complex to be calculated from first principles, using a single simulation. Our method involves the use of a novel reaction coordinate that swaps a ligand bound to a protein with an equivalent volume of bulk water. This water-swap reaction coordinate is built using an identity constraint, which identifies a cluster of water molecules from bulk water that occupies the same volume as the ligand in the protein active site. A dual topology algorithm is then used to swap the ligand from the active site with the identified water cluster from bulk water. The free energy is then calculated using replica exchange thermodynamic integration. This returns the free energy change of simultaneously transferring the ligand to bulk water, as an equivalent volume of bulk water is transferred back to the protein active site. This, directly, is the absolute binding free energy. It should be noted that while this reaction coordinate models the binding process directly, an accurate force field and sufficient sampling are still required to allow for the binding free energy to be predicted correctly. In this paper we present the details and development of this method, and demonstrate how the potential of mean force along the water-swap coordinate can be improved by calibrating the soft-core Coulomb and Lennard-Jones parameters used for the dual topology calculation. The optimal parameters were applied to calculations of protein-ligand binding free energies of a neuraminidase inhibitor (oseltamivir), with these results compared to experiment. These results demonstrate that the water-swap coordinate provides a viable and potentially powerful new route for the prediction of protein-ligand binding free energies.     

Supramolecular displacement-mediated activation of a silent fluorescence probe for label-free ligand screening

We report a new approach for the rapid screening of analyte binding affinities for a target protein. We demonstrate that a molecular probe, with a pro-fluorophore substrate and ligand moieties, can be hindered from enzymatic access when bound to the target protein. When analytes displace the probe from the protein's binding pocket, a fluorescence profile is generated. This profile is used to discriminate analytes based on their relative binding affinities.     

Intercalation of a Heterocyclic Ligand between Quartets in a G-Rich Tetrahelical Structure

YES G-rich oligonucleotide VK2 folds into an AGCGA-quadruplex tetrahelical structure distinct and significantly different from G-quadruplexes, even though it contains four G₃ tracts. Herein, a bis-quinolinium ligand 360A with high affinity for G-quadruplex structures and selective telomerase inhibition is shown to strongly bind to VK2. Upon binding, 360A does not induce a conformational switch from VK2 to an expected G-quadruplex. In contrast, NMR structural study revealed formation of a well-defined VK2–360A complex with a 1:1 binding stoichiometry, in which 360A intercalates between GAGA- and GCGC-quartets in the central cavity of VK2. This is the first high-resolution structure of a G-quadruplex ligand intercalating into a G-rich tetrahelical fold. This unique mode of ligand binding into tetrahelical DNA architecture offers insights into the stabilization of an AGCGA-quadruplex by a heterocyclic ligand and provides guidelines for rational design of novel VK2 binding molecules with selectivity for different DNA secondary structures.     

Investigating bidentate and tridentate carbamoylmethylphosphine oxide ligand interactions with rare-Earth elements using electrospray ionization quadrupole ion trap mass spectrometry

Electrospray ionization (ESI) quadrupole ion trap mass spectrometry (QIT-MS) and collisionally activated dissociation (CAD) were used to evaluate the rare-earth binding properties of two hydrophobic carbamoylmethylphosphine oxide (CMPO) ligands, the normal bidentate variety, (t-BuC6H4)2P(O)CH2C(O)N(i-Bu)2 (A), a new potentially tridentate extractant, (t-BuC6H4)2P(O)CH[CH2C(O)N(i-Bu)2]C(O)N(i-Bu)2 (B), and tributyl phosphate. The mass spectral results obtained from analysis of 1% HNO3/methanol solution containing the ligands and dissolved lanthanide salts reveal that the favorable stoichiometries of the ligand/metal/nitrate complexes are 2:1:2 for the bidentate ligand A, 1:1:2 for the tridentate ligand B, and 3:1:2 for the monodentate tributyl phosphate. These observed stoichiometries correlate with the number of available binding sites on each ligand as well as with potential steric effects. Energy-variable collisionally activated dissociation experiments showed that for the 2:1:2 complexes involving ligand A or B, as the ionic radius of the bound metal decreased, the removal of nitric acid required less energy and resulted in less extensive spontaneous solvent coordination. This experimental trend suggests that, as the ionic radius of the lanthanide ion decreases, a pair of the carbamoylmethylphosphine ligands is able to more completely solvate the bound metal ion thereby weakening the nitrate-metal interaction.     

Conformational selection of protein kinase A revealed by flexible-ligand flexible-protein docking

Protein kinases have high structural plasticity: their structure can change significantly, depending on what ligands are bound to them. Rigid-protein docking methods are not capable of describing such effects. Here, we present a new flexible-ligand flexible-protein docking model in which the protein can adopt conformations between two extremes observed experimentally. The model utilized a molecular dynamics-based simulated annealing cycling protocol and a distance-dependent dielectric model to perform docking. By testing this model on docking four diverse ligands to protein kinase A, we found that the ligands were able to dock successfully to the protein with the proper conformations of the protein induced. By imposing relatively soft conformational restraints to the protein during docking, this model reduced computational costs yet permitted essential conformational changes that were essential for these inhibitors to dock properly to the protein. For example, without adequate movement of the glycine-rich loop, it was difficult for the ligands to move from the surface of the protein to the binding site. In addition, these simulations called for better ways to compare simulation results with experiment other than using the popular root-mean-square deviation between the structure of a ligand in a docking pose and that in experiment because the structure of the protein also changed. In this work, we also calculated the correlation coefficient between protein-ligand/protein-protein distances in the docking structure and those in the crystal structure to check how well a ligand docked into the binding site of the protein and whether the proper conformation of the protein was induced.