Relative stability of the open and closed conformations of the active site loop of streptavidin

The eight-residue surface loop, 45-52 (Ser, Ala, Val, Gly, Asn, Ala, Glu, Ser), of the homotetrameric protein streptavidin has a "closed" conformation in the streptavidin-biotin complex, where the corresponding binding affinity is one of the strongest found in nature (ΔG ～ -18 kcal∕mol). However, in most of the crystal structures of apo (unbound) streptavidin, the loop conformation is "open" and typically exhibits partial disorder and high B-factors. Thus, it is plausible to assume that the loop structure is changed from open to closed upon binding of biotin, and the corresponding difference in free energy, ΔF = F(open) - F(closed) in the unbound protein, should therefore be considered in the total absolute free energy of binding. ΔF (which has generally been neglected) is calculated here using our "hypothetical scanning molecular-dynamics" (HSMD) method. We use a protein model in which only the atoms closest to the loop are considered (the "template") and they are fixed in the x-ray coordinates of the free protein; the x-ray conformation of the closed loop is attached to the same (unbound) template and both systems are capped with the same sphere of TIP3P water. Using the force field of the assisted model building with energy refinement (AMBER), we carry out two separate MD simulations (at temperature T = 300 K), starting from the open and closed conformations, where only the atoms of the loop and water are allowed to move (the template-water and template-loop interactions are considered). The absolute F(open) and F(closed) (of loop + water) are calculated from these trajectories, where the loop and water contributions are obtained by HSMD and a thermodynamic integration (TI) process, respectively. The combined HSMD-TI procedure leads to total (loop + water) ΔF = -27.1 ± 2.0 kcal∕mol, where the entropy TΔS constitutes 34% of ΔF, meaning that the effect of S is significant and should not be ignored. Also, ΔS is positive, in accord with the high flexibility of the open loop observed in crystal structures, while the energy ΔE is unexpectedly negative, thus also adding to the stability of the open loop. The loop and the 250 capped water molecules are the largest system studied thus far, which constitutes a test for the efficiency of HSMD-TI; this efficiency and technical issues related to the implementation of the method are also discussed. Finally, the result for ΔF is a prediction that will be considered in the calculation of the absolute free energy of binding of biotin to streptavidin, which constitutes our next project.     

Calculation of the entropy and free energy of peptides by molecular dynamics simulations using the hypothetical scanning molecular dynamics method

Hypothetical scanning (HS) is a method for calculating the absolute entropy S and free energy F from a sample generated by any simulation technique. With this approach each sample configuration is reconstructed with the help of transition probabilities (TPs) and their product leads to the configuration's probability, hence to the entropy. Recently a new way for calculating the TPs by Monte Carlo (MC) simulations has been suggested, where all system interactions are taken into account. Therefore, this method--called HSMC--is in principle exact where the only approximation is due to insufficient sampling. HSMC has been applied very successfully to liquid argon, TIP3P water, self-avoiding walks on a lattice, and peptides. Because molecular dynamics (MD) is considered to be significantly more efficient than MC for a compact polymer chain, in this paper HSMC is extended to MD simulations as applied to peptides. Like before, we study decaglycine in vacuum but for the first time also a peptide with side chains, (Val)(2)(Gly)(6)(Val)(2). The transition from MC to MD requires implementing essential changes in the reconstruction process of HSMD. Results are calculated for three microstates, helix, extended, and hairpin. HSMD leads to very stable differences in entropy TDeltaS between these microstates with small errors of 0.1-0.2 kcal/mol (T=100 K) for a wide range of calculation parameters with extremely high efficiency. Various aspects of HSMD and plans for future work are discussed.     

pKa analysis for the zinc-bound water in human carbonic anhydrase II: Benchmark for "multiscale" QM/MM simulations and mechanistic implications

To quantitatively explore the applicability of the generalized solvent boundary potential (GSBP) based QM/MM approach as a "multiscale" framework for studying chemical reactions in biomolecules, the structural and energetic properties of the Human Carbonic Anhydrase II (CAII) are analyzed and compared to those from periodic boundary condition (PBC) simulations and available experimental data. Although the atomic fluctuations from GSBP based simulations are consistently lower compared to those from PBC simulations or crystallographic data, the fluctuations and internal coordinate distributions for residues in the proximity of the active site as well as diffusion constants of active-site water molecules are fairly well described by GSBP simulations. The pKa of the zinc-bound water, calculated with a SCC-DFTB/MM-GSBP based thermodynamic integration approach, agrees well with experiments for the wild type CAII. For the E106Q mutant, however, a 9 pKa unit downward shift relative to the wild type is found in contrast with previous experiments that found little change. This dramatic discrepancy signals a possible change in the mechanism for the interconversion between CO2/HCO3- in the E106Q mutant, which may be similar to the bicarbonate mediated mechanism proposed for the Co2+ substituted CAII (J. Am. Chem. Soc. 2001, 123, 5861).1 The study highlights pKa analyses as a valuable approach for quantitatively validating the computational model for complex biomolecules as well as for revealing energetic properties intimately related to the chemical process of interest.     

Conformational sensitivity analysis of FKBP–FK506/rapamycin complexes

Sensitivity analysis techniques are applied to the FKBP–FK506 and FKBP–rapamycin complexes to quantify the conformational relationships between FKBP and its ligands. Crystal structures of the two FKBP complexes are energy minimized in the Amber force field using a continuum solvent model, and derived Green's function sensitivity coefficients are developed to describe the relationship between the ?, ψ, and χ₁ torsional angles of the FKBP residues and the bound ligand macrocycle torsional angles. Sensitivity analysis is applied to the entire FKBP structure and reveals that the local conformation of the residues of the 80s and 50s loops and of the active site are sensitive to the ligand conformation. The analysis also reveals that the torsional angles controlling the orientation of the amide and keto carbonyls of FK506 are sensitive to the aromatic side chains in the FKBP carbonyl binding pocket. © 1994 by John Wiley & Sons, Inc.     

WATsite: Hydration site prediction program with PyMOL interface

Water molecules that mediate protein–ligand interactions or are released from the binding site on ligand binding can contribute both enthalpically and entropically to the free energy of ligand binding. To elucidate the thermodynamic profile of individual water molecules and their potential contribution to ligand binding, a hydration site analysis program WATsite was developed together with an easy‐to‐use graphical user interface based on PyMOL. WATsite identifies hydration sites from a molecular dynamics simulation trajectory with explicit water molecules. The free energy profile of each hydration site is estimated by computing the enthalpy and entropy of the water molecule occupying a hydration site throughout the simulation. The results of the hydration site analysis can be displayed in PyMOL. A key feature of WATsite is that it is able to estimate the protein desolvation free energy for any user specified ligand. The WATsite program and its PyMOL plugin are available free of charge from http://people.pnhs.purdue.edu/~mlill/software . © 2014 Wiley Periodicals, Inc.     

Spectroscopic and density functional studies of the red copper site in nitrosocyanin: role of the protein in determining active site geometric and electronic structure

The electronic structure of the red copper site in nitrosocyanin is defined relative to that of the well understood blue copper site of plastocyanin by using low-temperature absorption, circular dichroism, magnetic circular dichroism, resonance Raman, EPR and X-ray absorption spectroscopies, combined with DFT calculations. These studies indicate that the principal electronic structure change in the red copper site is the sigma rather than the pi donor interaction of the cysteine sulfur with the Cu 3d(x2-y2) redox active molecular orbital (RAMO). Further, MCD data show that there is an increase in ligand field strength due to an increase in coordination number, whereas resonance Raman spectra indicate a weaker Cu-S bond. The latter is supported by the S K-edge data, which demonstrate a less covalent thiolate interaction with the RAMO of nitrosocyanin at 20% relative to plastocyanin at 38%. EXAFS results give a longer Cu-S(Cys) bond distance in nitrosocyanin (2.28 A) compared to plastocyanin (2.08 A) and also show a large change in structure with reduction of the red copper site. The red copper site is the only presently known blue copper-related site with an exogenous water coordinated to the copper. Density functional calculations reproduce the experimental properties and are used to determine the specific protein structure contributions to exogenous ligand binding in red copper. The relative orientation of the CuNNS and the CuSC(beta) planes (determined by the protein sequence) is found to be key in generating an exchangeable coordination position at the red copper active site. The exogenous water ligation at the red copper active site greatly increases the reorganization energy (by approximately 1.0 eV) relative to that of the blue copper protein site, making the red site unfavorable for fast outer-sphere electron transfer, while providing an exchangeable coordination position for inner-sphere electron transfer.     

Spectrometric study of the interaction between Alpinetin and bovine serum albumin using chemometrics approaches

The binding interaction of Alpinetin (APT) with bovine serum albumin (BSA) was studied by fluorescence, UV-visible and synchronous fluorescence spectroscopy (SFS) under simulated physiological conditions. The measured complex spectra were resolved by multivariate curve resolution-alternating least squares (MCR-ALS), yielding a host of data and information, which otherwise would have been impossible to obtain. The extracted profiles corresponded to the spectra of the single species in the APT/BSA mixture. In addition, the presence of the APT-BSA complex was demonstrated, and it was shown that the associated quenching of the fluorescence from the BSA protein resulted from the formation of APT-BSA complex via a static mechanism. The binding constant (K_a(ave) = 2.34 × 10⁶ L mol⁻¹) and the number of sites (n = 1) were obtained by fluorescence methods as were the thermodynamic parameters (ΔH⁰, ΔS⁰ and ΔG⁰). This work suggested that the principal binding between APT to BSA was facilitated by hydrophobic interactions. The thermodynamic parameters for APT were compared to those from the structurally similar Chrysin and Wogonin molecules. It appeared that the entropy parameters were relatively more affected by the small structural changes. SFS from the interaction of BSA and APT showed that the ligand affected the conformation of BSA. The competitive interaction of APT and site makers with BSA indicated site I as the binding area of APT in BSA.     

Improved mapping of protein binding sites

Computational mapping methods place molecular probes – small molecules or functional groups – on a protein surface in order to identify the most favorable binding positions by calculating an interaction potential. Mapping is an important step in a number of flexible docking and drug design algorithms. We have developed improved algorithms for mapping protein surfaces using small organic molecules as molecular probes. The calculations reproduce the binding of eight organic solvents to lysozyme as observed by NMR, as well as the binding of four solvents to thermolysin, in good agreement with x-ray data. Application to protein tyrosine phosphatase 1B shows that the information provided by the mapping can be very useful for drug design. We also studied why the organic solvents bind in the active site of proteins, in spite of the availability of alternative pockets that can very tightly accommodate some of the probes. A possible explanation is that the binding in the relatively large active site retains a number of rotational states, and hence leads to smaller entropy loss than the binding elsewhere else. Indeed, the mapping reveals that the clusters of the ligand molecules in the protein's active site contain different rotational-translational conformers, which represent different local minima of the free energy surface. In order to study the transitions between different conformers, reaction path and molecular dynamics calculations were performed. Results show that most of the rotational states are separated by low free energy barriers at the experimental temperature, and hence the entropy of binding in the active site is expected to be high.     

Ligand coordinate analysis of SC-558 from the active site to the surface of COX-2: a molecular dynamics study

We have performed a ligand coordinate analysis to monitor the movement of the inhibitor SC-558 from the active site of the COX-2 protein to the exterior using molecular dynamics techniques. This study provides an insight into the intermolecular interactions formed by the ligand during this journey. The published crystal structure of COX-2 with SC-558 in the active site (1cx2) was taken, and the ligand was moved incrementally in 13 steps. At each of these points on the path, exhaustive minimization and dynamics calculations were performed. The role of water was found to be important in these computations. An average structure was obtained from 250 conformations at each point and minimized. At each point on the path, the 10 lowest-energy conformations were also selected; a consideration of the average and lowest conformations provides fine details on the consistency of specific and strong interactions, and also on the geometry of the ligand. The movement of the ligand through the protein may be divided into three stages that are distinguished from each other because of energy and geometry discontinuities in both the ligand and the protein. The first of these covers the region between the active site and the point at 5.8 A from it. The second, which covers the distance between 8.2 and 10.0 A and is associated with maximum energetic and structural instability, is of critical importance. The third stage covers the distance between 10.5 A and the exterior and represents a stage of increasing hydration and expulsion of the ligand from the protein. Our results provide a confirmation for the existence of a shallow cavity near the protein surface in which the ligand is bound reversibly. By examining the residues that show maximum mobility, one obtains an idea of the gating mechanism that governs the entry and exit of the protein into or from the deep pocket that contains the active site. We note, however, that the variation of the root-mean-square deviation of all residues begins to increase almost as soon as the ligand leaves the active site, and even before there are any changes in the gate inter-residue distances. This loosening of the protein even before the gate opens might be a part of the enthalpy-entropy balance that accompanies the ligand's passage through the protein. Our results provide an energy profile of the ligand during its entry/exit into/from the protein and can, in principle, enable one to assess the residence time, which in turn may be associated or indirectly correlated with adverse cardiovascular side effects of nonsteroidal anti-inflammatory drugs. We believe that similar analyses for other selected COX-2-specific inhibitors will provide a measure (or prediction) of possible toxicity effects.     

Native FKBP12 engineering by ligand-directed tosyl chemistry: labeling properties and application to photo-cross-linking of protein complexes in vitro and in living cells

The ability to modify target "native" (endogenous) proteins selectively in living cells with synthetic molecules should provide powerful tools for chemical biology. To this end, we recently developed a novel protein labeling technique termed ligand-directed tosyl (LDT) chemistry. This method uses labeling reagents in which a protein ligand and a synthetic probe are connected by a tosylate ester group. We previously demonstrated its applicability to the selective chemical labeling of several native proteins in living cells and mice. However, many fundamental features of this chemistry remain to be studied. In this work, we investigated the relationship between the LDT reagent structure and labeling properties by using native FK506-binding protein 12 (FKBP12) as a target protein. In vitro experiments revealed that the length and rigidity of the spacer structure linking the protein ligand and the tosylate group have significant effects on the overall labeling yield and labeling site. In addition to histidine, which we reported previously, tyrosine and glutamate residues were identified as amino acids that are modified by LDT-mediated labeling. Through the screening of various spacer structures, piperazine was found to be optimal for FKBP12 labeling in terms of labeling efficiency and site specificity. Using a piperazine-based LDT reagent containing a photoreactive probe, we successfully demonstrated the labeling and UV-induced covalent cross-linking of FKBP12 and its interacting proteins in vitro and in living cells. This study not only furthers our understanding of the basic reaction properties of LDT chemistry but also extends the applicability of this method to the investigation of biological processes in mammalian cells.     

The effect of water displacement on binding thermodynamics: concanavalin A

Interactions at the binding interface of biomolecular complexes are often mediated by ordered water molecules. In this work, we considered two concanavalin A-carbohydrate complexes. In the first, a water molecule is buried at the binding interface. In the second, this water molecule is displaced by a modification of the ligand (Clarke, C.; Woods, R. J.; Gluska, J.; Cooper, A.; Nutley, M. A.; Boons, G. J. J. Am. Chem. Soc. 2001, 123, 12238-12247). We computed the contribution of this water molecule to the thermodynamic properties using statistical mechanical formulas for the energy and entropy and molecular dynamics simulations. Other contributions to the binding affinity, including desolvation, entropy of conformational restriction, and interaction between the ligand and protein, were also computed. The thermodynamic consequences of displacement of the ordered water molecule by ligand modification were in qualitative agreement with experimental data. The free energy contribution of the water molecule (-17.2 kcal/mol; -19.2 enthalpic and +2 entropic) was nearly equivalent to the additional protein-ligand interactions in trimannoside 2 (-18.9 kcal/mol). The two structural ions interact more strongly with the water than with the hydroxyl of trimannoside 2, thus favoring trimannoside 1. The contributions from desolvation and conformational entropy are much smaller but significant, compared to the binding free energy difference. The picture that emerges is that the final outcome of water displacement is sensitive to the details of the binding site and cannot be predicted by simple empirical rules.     

Spectroscopic definition of the ferroxidase site in M ferritin: comparison of binuclear substrate vs cofactor active sites

Maxi ferritins, 24 subunit protein nanocages, are essential in humans, plants, bacteria, and other animals for the concentration and storage of iron as hydrated ferric oxide, while minimizing free radical generation or use by pathogens. Formation of the precursors to these ferric oxides is catalyzed at a nonheme biferrous substrate site, which has some parallels with the cofactor sites in other biferrous enzymes. A combination of circular dichroism (CD), magnetic circular dichroism (MCD), and variable-temperature, variable-field MCD (VTVH MCD) has been used to probe Fe(II) binding to the substrate active site in frog M ferritin. These data determined that the active site within each subunit consists of two inequivalent five-coordinate (5C) ferrous centers that are weakly antiferromagnetically coupled, consistent with a mu-1,3 carboxylate bridge. The active site ligand set is unusual and likely includes a terminal water bound to each Fe(II) center. The Fe(II) ions bind to the active sites in a concerted manner, and cooperativity among the sites in each subunit is observed, potentially providing a mechanism for the control of ferritin iron loading. Differences in geometric and electronic structure--including a weak ligand field, availability of two water ligands at the biferrous substrate site, and the single carboxylate bridge in ferritin--coincide with the divergent reaction pathways observed between this substrate site and the previously studied cofactor active sites.     

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.     

A water sluice is generated in the active site of bovine lens leucine aminopeptidase

Nature has provided the binuclear zinc based active site of bovine lens leucine aminopeptidase (blLAP) with two water channels: one for substrate docking and a much smaller one (function unknown) above Zn1. In addition, Zn1 possesses an unusual pentacoordinate geometry with a loosely bound carbonyl ligand (Ala333). Extensive DFT calculations on a model of the active site provide first mechanistic implications for these structural features. The weakly bound carbonyl ligand is capable of functioning as a "traffic cop" to direct water molecules coming from the small channel into the heart of the active site. A water sluice is thus generated that is capable of repeatedly providing a series of nucleophilic active "Zn-OH" functionalities.     

SOS-NMR: a saturation transfer NMR-based method for determining the structures of protein-ligand complexes

An NMR-based alternative to traditional X-ray crystallography and NMR methods for structure-based drug design is described that enables the structure determination of ligands complexed to virtually any biomolecular target regardless of size, composition, or oligomeric state. The method utilizes saturation transfer difference (STD) NMR spectroscopy performed on a ligand complexed to a series of target samples that have been deuterated everywhere except for specific amino acid types. In this way, the amino acid composition of the ligand-binding site can be defined, and, given the three-dimensional structure of the protein target, the three-dimensional structure of the protein-ligand complex can be determined. Unlike earlier NMR methods for solving the structures of protein-ligand complexes, no protein resonance assignments are necessary. Thus, the approach has broad potential applications--especially in cases where X-ray crystallography and traditional NMR methods have failed to produce structural data. The method is called SOS-NMR for structural information using Overhauser effects and selective labeling and is validated on two protein-ligand complexes: FKBP complexed to 2-(3'-pyridyl)-benzimidazole and MurA complexed to uridine diphosphate N-acetylglucosamine.     

Computation of binding free energy with molecular dynamics and grand canonical Monte Carlo simulations

The binding of a ligand to a receptor is often associated with the displacement of a number of bound water molecules. When the binding site is exposed to the bulk region, this process may be sampled adequately by standard unbiased molecular dynamics trajectories. However, when the binding site is deeply buried and the exchange of water molecules with the bulk region may be difficult to sample, the convergence and accuracy in free energy perturbation (FEP) calculations can be severely compromised. These problems are further compounded when a reduced system including only the region surrounding the binding site is simulated. To address these issues, we couple molecular dynamics (MD) with grand canonical Monte Carlo (GCMC) simulations to allow the number of water to fluctuate during an alchemical FEP calculation. The atoms in a spherical inner region around the binding pocket are treated explicitly while the influence of the outer region is approximated using the generalized solvent boundary potential (GSBP). At each step during thermodynamic integration, the number of water in the inner region is equilibrated with GCMC and energy data generated with MD is collected. Free energy calculations on camphor binding to a deeply buried pocket in cytochrome P450cam, which causes about seven water molecules to be expelled, are used to test the method. It concluded that solvation free energy calculations with the GCMC/MD method can greatly improve the accuracy of the computed binding free energy compared to simulations with fixed number of water.     

Accurate predictions of nonpolar solvation free energies require explicit consideration of binding-site hydration

Continuum solvation methods are frequently used to increase the efficiency of computational methods to estimate free energies. In this paper, we have evaluated how well such methods estimate the nonpolar solvation free-energy change when a ligand binds to a protein. Three different continuum methods at various levels of approximation were considered, viz., the polarized continuum model (PCM), a method based on cavity and dispersion terms (CD), and a method based on a linear relation to the solvent-accessible surface area (SASA). Formally rigorous double-decoupling thermodynamic integration was used as a benchmark for the continuum methods. We have studied four protein-ligand complexes with binding sites of varying solvent exposure, namely the binding of phenol to ferritin, a biotin analogue to avidin, 2-aminobenzimidazole to trypsin, and a substituted galactoside to galectin-3. For ferritin and avidin, which have relatively hidden binding sites, rather accurate nonpolar solvation free energies could be obtained with the continuum methods if the binding site is prohibited to be filled by continuum water in the unbound state, even though the simulations and experiments show that the ligand replaces several water molecules upon binding. For the more solvent exposed binding sites of trypsin and galectin-3, no accurate continuum estimates could be obtained, even if the binding site was allowed or prohibited to be filled by continuum water. This shows that continuum methods fail to give accurate free energies on a wide range of systems with varying solvent exposure because they lack a microscopic picture of binding-site hydration as well as information about the entropy of water molecules that are in the binding site before the ligand binds. Consequently, binding affinity estimates based upon continuum solvation methods will give absolute binding energies that may differ by up to 200 kJ/mol depending on the method used. Moreover, even relative energies between ligands with the same scaffold may differ by up to 75 kJ/mol. We have tried to improve the continuum solvation methods by adding information about the solvent exposure of the binding site or the hydration of the binding site, and the results are promising at least for this small set of complexes.