Absolute and relative binding free energy calculations of the interaction of biotin and its analogs with streptavidin using molecular dynamics/free energy perturbation approaches

Structure and dynamics of thioguanine-modified duplex DNA

Mercaptopurine and thioguanine, two of the most widely used antileukemic agents, exert their cytotoxic, therapeutic effects by being incorporated into DNA as deoxy-6-thioguanosine. However, the molecular mechanism(s) by which incorporation of these thiopurines into DNA translates into cytotoxicity is unknown. The solution structure of thioguanine-modified duplex DNA presented here shows that the effects of the modification on DNA structure were subtle and localized to the modified base pair. Specifically, thioguanine existed in the keto form, formed weakened Watson-Crick hydrogen bonds with cytosine and caused a modest approximately 10 degrees opening of the modified base pair toward the major groove. In contrast, thioguanine significantly altered base pair dynamics, causing an approximately 80-fold decrease in the base pair lifetime with cytosine compared with normal guanine. This perturbation was consistent with the approximately 6 degrees C decrease in DNA melting temperature of the modified oligonucleotide, the 1.13 ppm upfield shift of the thioguanine imino proton resonance, and the large increase in the exchange rate of the thioguanine imino proton with water. Our studies provide new mechanistic insight into the effects of thioguanine incorporation into DNA at the level of DNA structure and dynamics, provide explanations for the effects of thioguanine incorporation on the activity of DNA-processing enzymes, and provide a molecular basis for the specific recognition of thioguanine-substituted sites by proteins. These combined effects likely cooperate to produce the cellular responses that underlie the therapeutic effects of thiopurines.

Dissecting the streptavidin-biotin interaction by phage-displayed shotgun scanning

Shotgun scanning the streptavidin-biotin interaction identifies long-range hydrophobic interactions that contribute to one of the strongest naturally occurring noncovalent protein-ligand interactions. The femtomolar dissociation constant for this interaction makes it a useful model system to dissect the forces that govern high-affinity molecular recognition between proteins and small molecules. Shotgun scanning combines the diversity and in vitro binding selection of phage-displayed libraries with a binomial mutagenesis strategy. Libraries consist of proteins with the residues in multiple positions mutated to give a 1:1 ratio of alanine:wild type. Here, we use shotgun scanning to determine the functional contribution of the 38 C-terminal residues of streptavidin to the high-affinity interaction with biotin. The library pools were subjected to three rounds of selection for functional streptavidin variants that bind biotin and statistical analysis was used to assess side-chain contributions to biotin binding. The results demonstrate the utility of shotgun scanning for the dissection of receptor-small-molecule interactions. While shotgun scanning results were largely consistent with previous single-point, site-directed mutagenesis studies for residues in direct contact with biotin, residues distant from the biotin binding site have not previously been explored. Key streptavidin residues identified by shotgun scanning as contributors to the interaction with biotin include those with side chains that fill the beta barrel, residues at the tetramer interface, and second-sphere residues, which are reinforced by long-distance propagation of hydrophobic interactions.

Binding free energy differences in a TCR-peptide-MHC complex induced by a peptide mutation: a simulation analysis

Recognition by the T-cell receptor (TCR) of immunogenic peptides presented by class I major histocompatibility complexes (MHCs) is the determining event in the specific cellular immune response against virus-infected cells or tumor cells. It is of great interest, therefore, to elucidate the molecular principles upon which the selectivity of a TCR is based. These principles can in turn be used to design therapeutic approaches, such as peptide-based immunotherapies of cancer. In this study, free energy simulation methods are used to analyze the binding free energy difference of a particular TCR (A6) for a wild-type peptide (Tax) and a mutant peptide (Tax P6A), both presented in HLA A2. The computed free energy difference is 2.9 kcal/mol, in good agreement with the experimental value. This makes possible the use of the simulation results for obtaining an understanding of the origin of the free energy difference which was not available from the experimental results. A free energy component analysis makes possible the decomposition of the free energy difference between the binding of the wild-type and mutant peptide into its components. Of particular interest is the fact that better solvation of the mutant peptide when bound to the MHC molecule is an important contribution to the greater affinity of the TCR for the latter. The results make possible identification of the residues of the TCR which are important for the selectivity. This provides an understanding of the molecular principles that govern the recognition. The possibility of using free energy simulations in designing peptide derivatives for cancer immunotherapy is briefly discussed.

On the calculation of absolute macromolecular binding free energies

The standard framework for calculating the absolute binding free energy of a macromolecular association reaction A + B --> AB with an association constant K(AB) is to equate chemical potentials of the species on the left- and right-hand sides of this reaction and evaluate the chemical potentials from theory. This theory involves (usually hidden) assumptions about what constitutes the bound species, AB, and where the contribution of the solvent appears. We present here an alternative derivation that can be traced back to Bjerrum, in which the expectation value of K(AB) is obtained directly through the statistical mechanical method of evaluating its ensemble (Boltzmann-weighted) average. The generalized Bjerrum approach more clearly delineates: (i) the different contributions to binding; (ii) the origin of the much-discussed and somewhat controversial association entropy term; and (iii) where the solvent contribution appears. This approach also allows approximations required for practical evaluation of the binding constant in complex macromolecular systems, to be introduced in a well defined way. We provide an example, with application to test cases that illustrate a range of binding behavior.

A fast empirical approach to binding free energy calculations based on protein interface information

Three useful variables from the interfaces of 20 protein-protein complexes were investigated. These variables are the side-chain accessible number (N(b)), the number of hydrophilic pairs (N(pair)) and buried a polar solvent accessible surface areas (DeltaDeltaASA(apol)). An empirical model based on the three variables was developed to describe the free energy of protein associations. As the results show, the side-chain accessible numbers characterize the loss of side-chain conformational entropy of protein interactions and the effective empirical function presented here has great capability for estimating the binding free energy. It was found that the variables of interface information capture most of the significant features of protein-protein association. Also, we applied the model based on the variables as a rescoring function to docking simulations and found that it has the potential to distinguish the 'true' binding mode. It is clear that the simple and empirical scale developed here is an attractive target function for calculating binding free energy for various biological processes to rational protein design.

Free energy simulations come of age: protein-ligand recognition

In recent years, molecular dynamics simulations of biomolecular free energy differences have benefited from significant methodological advances and increased computer power. Applications to molecular recognition provide an understanding of the interactions involved that goes beyond, and is an important complement to, experimental studies. Poisson-Boltzmann electrostatic models provide a faster and simpler free energy method in cases where electrostatic interactions are important. We illustrate both molecular dynamics and Poisson-Boltzmann methods with a detailed study of amino acid recognition by aspartyl-tRNA synthetase, whose specificity is important for maintaining the integrity of the genetic code.

Contributions to the binding free energy of ligands to avidin and streptavidin

The free energy of binding of a ligand to a macromolecule is here formally decomposed into the (effective) energy of interaction, reorganization energy of the ligand and the macromolecule, conformational entropy change of the ligand and the macromolecule, and translational and rotational entropy loss of the ligand. Molecular dynamics simulations with implicit solvation are used to evaluate these contributions in the binding of biotin, biotin analogs, and two peptides to avidin and streptavidin. We find that the largest contribution opposing binding is the protein reorganization energy, which is calculated to be from 10 to 30 kcal/mol for the ligands considered here. The ligand reorganization energy is also significant for flexible ligands. The translational/rotational entropy is 4.5-6 kcal/mol at 1 M standard state and room temperature. The calculated binding free energies are in the correct range, but the large statistical uncertainty in the protein reorganization energy precludes precise predictions. For some complexes, the simulations show multiple binding modes, different from the one observed in the crystal structure. This finding is probably due to deficiencies in the force field but may also reflect considerable ligand flexibility.

Molecular dynamics simulations

Molecular dynamics simulations have become a standard tool for the investigation of biomolecules. Simulations are performed of ever bigger systems using more realistic boundary conditions and better sampling due to longer sampling times. Recently, realistic simulations of systems as complex as transmembrane channels have become feasible. Simulations aid our understanding of biochemical processes and give a dynamic dimension to structural data; for example, the transformation of harmless prion protein into the disease-causing agent has been modeled.

Theoretical and experimental studies of biotin analogues that bind almost as tightly to streptavidin as biotin

We have used a newly developed qualitative computational approach, PROFEC (Pictorial Representation of Free Energy Changes), to visualize the areas of the ligand biotin where modifications of its structure might lead to tighter binding to the protein streptavidin. The PROFEC analysis, which includes protein flexibility and ligand solvation/desolvation, led to the suggestion that the pro-9R hydrogen atom of biotin, which is in alpha-position to the CO(2)(-) group, might be changed to a larger group and lead to better binding with streptavidin and avidin. Free energy calculations supported this suggestion and predicted that the methyl analogue should bind approximately 3 kcal/mol more tightly to streptavidin, with this difference coming exclusively from the relative desolvation free energy of the ligand. The PROFEC analysis further suggested little or no improvement for changing the pro-9S hydrogen atom to a methyl group, and great reduction in changing the ureido N-H groups to N-CH(3). Stimulated by these results, we synthesized 9R-methylbiotin and 9S-methylbiotin, and their binding free energies and enthalpies were measured for interaction with streptavidin and avidin, respectively. In contrast to the calculated results, experiments found both 9-methylbiotin isomers to bind more weakly to streptavidin than biotin. The calculated preference for the binding of the 9R- over the 9S-stereoisomer was observed. In addition, 9-methylbiotin is considerably less soluble in water than biotin, as predicted by the calculation, and the 9R isomer is, to our knowledge, thus far the tightest binding analogue of biotin to streptavidin. Subsequently, X-ray structures of the complexes between streptavidin and both 9R- and 9S-methylbiotin were determined, and the structures were consistent with those used in the free energy calculations. Thus, the reason for the discrepancy between the calculated and experimental binding free energy does not lie in unusual binding modes for the 9-methylbiotins.

Contribution of translational and rotational motions to molecular association in aqueous solution

Much uncertainty and controversy exist regarding the estimation of the enthalpy, entropy, and free energy of overall translational and rotational motions of solute molecules in aqueous solutions, quantities that are crucial to the understanding of molecular association/recognition processes and structure-based drug design. A critique of the literature on this topic is given that leads to a classification of the various views. The major stumbling block to experimentally determining the translational/rotational enthalpy and entropy is the elimination of vibrational perturbations from the measured effects. A solution to this problem, based on a combination of energy equi-partition and enthalpy-entropy compensation, is proposed and subjected to verification. This method is then applied to analyze experimental data on the dissociation/unfolding of dimeric proteins. For one translational/rotational unit at 1 M standard state in aqueous solution, the results for enthalpy (H degrees (tr)), entropy (S degrees (tr)), and free energy (G degrees (tr)) are H (degrees) (tr) = 4.5 +/- 1.5RT, S (degrees) (tr) = 5 +/- 4R, and G (degrees) (tr) = 0 +/- 5RT. Therefore, the overall translational and rotational motions make negligible contribution to binding affinity (free energy) in aqueous solutions at 1 M standard state.

Free energy decomposition of protein-protein interactions

A free energy decomposition scheme has been developed and tested on antibody-antigen and protease-inhibitor binding for which accurate experimental structures were available for both free and bound proteins. Using the x-ray coordinates of the free and bound proteins, the absolute binding free energy was computed assuming additivity of three well-defined, physical processes: desolvation of the x-ray structures, isomerization of the x-ray conformation to a nearby local minimum in the gas-phase, and subsequent noncovalent complex formation in the gas phase. This free energy scheme, together with the Generalized Born model for computing the electrostatic solvation free energy, yielded binding free energies in remarkable agreement with experimental data. Two assumptions commonly used in theoretical treatments; viz., the rigid-binding approximation (which assumes no conformational change upon complexation) and the neglect of vdW interactions, were found to yield large errors in the binding free energy. Protein-protein vdW and electrostatic interactions between complementary surfaces over a relatively large area (1400--1700 A(2)) were found to drive antibody-antigen and protease-inhibitor binding.

Evaluation of protein-protein association energies by free energy perturbation calculations

The association energy upon binding of different amino acids in the specificity pocket of trypsin was evaluated by free energy perturbation calculations on complexes between bovine trypsin (BT) and bovine pancreatic trypsin inhibitor (BPTI). Three simulations of mutations of the primary binding residue (P(1)) were performed (P(1)-Ala to Gly, P(1)-Met to Gly and P(1)-Met to Ala) and the resulting differences in association energy (DeltaDeltaG(a)) are 2. 28, 5.08 and 2.93 kcal/mol for P(1)-Ala to Gly, P(1)-Met to Gly and to Ala with experimental values of 1.71, 4.62 and 2.91 kcal/mol, respectively. The calculated binding free energy differences are hence in excellent agreement with the experimental binding free energies. The binding free energies, however, were shown to be highly dependent on water molecules at the protein-protein interface and could only be quantitatively estimated if the correct number of such water molecules was included. Furthermore, the cavities that were formed when a large amino acid side-chain is perturbed to a smaller one seem to create instabilities in the systems and had to be refilled with water molecules in order to obtain reliable results. In addition, if the protein atoms that were perturbed away were not replaced by water molecules, the simulations dramatically overestimated the initial state of the free energy perturbations.

Streptavidin-biotin binding energetics

The high affinity energetics in the streptavidin-biotin system provide an excellent model system for studying how proteins balance enthalpic and entropic components to generate an impressive overall free energy for ligand binding. We review here concerted site-directed mutagenesis, biophysical, and computational studies of aromatic and hydrogen bonding interaction energetics between streptavidin and biotin. These results also have provided insight into how streptavidin builds a large activation barrier to dissociation by managing the enthalpic and entropic activation components. Finally, we review recent studies of the biotin dissociation pathway that address the fundamental question of how ligands exit protein binding pockets.

Streptavidin-binding and -dimerizing ligands discovered by phage display, topochemistry, and structure-based design

Structural and mechanistic determinants of affinity of streptavidin-binding peptide ligands discovered by phage display are reviewed along with the use of streptavidin as a paradigm for structure-based design. A novel way of producing protein-dimerizing ligands in the streptavidin model system is discussed, in which crystal packing topochemically mediates or even catalyzes dimerization of adjacent bound ligands whose reactive ligating groups are presented toward one another in productive orientations in the crystal lattice. Finally, through crystallography on a set of streptavidin complexes with small molecule and peptide ligands at multiple pHs in two space groups, the mechanism by which ligands enhance intersubunit stabilization of the streptavidin tetramer is probed.

Direct force measurements of the streptavidin-biotin interaction

The interaction between streptavidin and its ligand, biotin, were studied by direct force measurements. The complimentary approaches of surface force apparatus (SFA) and atomic force microscopy (AFM) were used to elucidate both long-range and short-range adhesive interactions of the streptavidin biotin interaction. The high spatial resolution of the SFA provided a detailed profile of the intersurface forces of apposing surfaces functionalized with streptavidin and biotin. Measurements obtained by the SFA corresponded to long and intermediate-range forces that are important in determining ligand receptor association. AFM was used to measure the unbinding force of individual streptavidin biotin complexes. These measurements revealed the short-range interactions (i.e. hydrophobic and hydrogen bonding forces) that stabilize the intermolecular bond.

Predicting binding affinities of protein ligands from three-dimensional models: application to peptide binding to class I major histocompatibility proteins

A simple and fast free energy scoring function (Fresno) has been developed to predict the binding free energy of peptides to class I major histocompatibility (MHC) proteins. It differs from existing scoring functions mainly by the explicit treatment of ligand desolvation and of unfavorable protein-ligand contacts. Thus, it may be particularly useful in predicting binding affinities from three-dimensional models of protein-ligand complexes. The Fresno function was independently calibrated for two different training sets: (a) five HLA-A0201-peptide structures, which had been determined by X-ray crystallography, and (b) three-dimensional models of 37 H-2K(k)-peptide structures, which had been obtained by knowledge-based homology modeling. For both training sets, a good cross-validated fit to experimental binding free energies was obtained with predictive errors of 3-3.5 kJ/mol. As expected, lipophilic interactions were found to contribute the most to HLA-A0201-peptide interactions, whereas H-bonding predominates in H-2K(k) recognition. Both cross-validated models were afterward used to predict the binding affinity of a test set of 26 peptides to HLA-A0204 (an HLA allele closely related to HLA-A0201) and of a series of 16 peptides to H-2K(k). Predictions were more accurate for HLA-A2-binding peptides as the training set had been built from experimentally determined structures. The average error in predicting the binding free energy of the test peptides was 3.1 kJ/mol. For the homology model-derived equation, the average error in predicting the binding free energy of peptides to K(k) was significantly higher (5.4 kJ/mol) but still very acceptable. The present scoring function is thus able to predict with a good accuracy binding free energies from three-dimensional models, at the condition that the backbone coordinates of the MHC-bound peptide have first been determined with an accuracy of about 1-1.5 A. Furthermore, it may be easily recalibrated for any protein-ligand complex.

Free energy calculations and molecular dynamics simulations of wild-type and variants of the DNA-EcoRI complex

Molecular dynamics simulations and free energy calculations of the wild-type EcoRI-DNA complex and several variants have been performed in aqueous solvent. In general, he theoretical estimations of the free energy differences (DeltaDeltaA) qualitatively agree well with the corresponding experimental data. The modifications which were experimentally found unfavorable compared to the wild-type complex were also found to be so in theoretical estimates. The mutant where the amino group of the base Ade(6) was replaced by a hydrogen atom eliminating one H-bond between the DNA and the protein, was experimentally found to be more stable than the wild-type complex. It was speculated that the modification also caused a structural relaxation in the DNA making DeltaDeltaA favorable. Our theoretical estimate yields a positive DeltaDeltaA in this case, but the difference is small, and no significant local structural relaxation was observed. The major H-bonds between the DNA and the protein in the wild-type complex are found to be maintained in the different mutants although the specific and non-specific interaction energies between the interacting the DNA bases and the protein residues are different in different mutants. The interaction pattern of the other nearby nucleotides are significantly influenced by each modification. Thus, the alteration of the non-specific interactions may also play an indirect role in determining the specificity of the complex. The interaction of the Gua(4) of the DNA with the protein is found to be most sensitive to any alteration in the recognition site. Because Gua(4) is the nucleotide closest to the scissile bond, this extra sensitivity seems to play an important role in altering the functional activity of the complex.

The maximal affinity of ligands

We explore the question of what are the best ligands for macromolecular targets. A survey of experimental data on a large number of the strongest-binding ligands indicates that the free energy of binding increases with the number of nonhydrogen atoms with an initial slope of approximately -1.5 kcal/mol (1 cal = 4.18 J) per atom. For ligands that contain more than 15 nonhydrogen atoms, the free energy of binding increases very little with relative molecular mass. This nonlinearity is largely ascribed to nonthermodynamic factors. An analysis of the dominant interactions suggests that van der Waals interactions and hydrophobic effects provide a reasonable basis for understanding binding affinities across the entire set of ligands. Interesting outliers that bind unusually strongly on a per atom basis include metal ions, covalently attached ligands, and a few well known complexes such as biotin-avidin.

A structural snapshot of an intermediate on the streptavidin-biotin dissociation pathway

It is currently unclear whether small molecules dissociate from a protein binding site along a defined pathway or through a collection of dissociation pathways. We report herein a joint crystallographic, computational, and biophysical study that suggests the Asp-128 --> Ala (D128A) streptavidin mutant closely mimics an intermediate on a well-defined dissociation pathway. Asp-128 is hydrogen bonded to a ureido nitrogen of biotin and also networks with the important aromatic binding contacts Trp-92 and Trp-108. The Asn-23 hydrogen bond to the ureido oxygen of biotin is lengthened to 3.8 A in the D128A structure, and a water molecule has moved into the pocket to replace the missing carboxylate interaction. These alterations are accompanied by the coupled movement of biotin, the flexible binding loop containing Ser-45, and the loop containing the Ser-27 hydrogen bonding contact. This structure closely parallels a key intermediate observed in a potential of mean force-simulated dissociation pathway of native streptavidin, where the Asn-23 hydrogen bond breaks first, accompanied by the replacement of the Asp-128 hydrogen bond by an entering water molecule. Furthermore, both biotin and the flexible loop move in a concerted conformational change that closely approximates the D128A structural changes. The activation and thermodynamic parameters for the D128A mutant were measured and are consistent with an intermediate that has traversed the early portion of the dissociation reaction coordinate through endothermic bond breaking and concomitant gain in configurational entropy. These composite results suggest that the D128A mutant provides a structural "snapshot" of an early intermediate on a relatively well-defined dissociation pathway for biotin.

Scoring functions: a view from the bench

Computational approaches to drug design are presently hindered by the complexity of the physical chemistry which underlies weak, non-covalent interactions between protein targets and small molecule ligands. Although a number of programs are now available for the design of novel potential ligands, it remains a key problem to rank these rapidly and reliably by estimated binding affinity. Such a step is necessary to select only the most promising candidates for synthesis and experimental characterisation. To calculate ligand affinity quickly and reliably is an extremely difficult problem, but it may well prove possible to estimate sufficiently accurately given an appropriate set of parameters to 'score' individual protein-ligand interactions. Improvements in the situation will require a wider set of thermodynamically characterised systems than is currently available.

Protein-protein recognition: an experimental and computational study of the R89K mutation in Raf and its effect on Ras binding

Binding of the protein Raf to the active form of Ras promotes activation of the MAP kinase signaling pathway, triggering cell growth and differentiation. Raf/Arg89 in the center of the binding interface plays an important role determining Ras-Raf binding affinity. We have investigated experimentally and computationally the Raf-R89K mutation, which abolishes signaling in vivo. The binding to [gamma-35S]GTP-Ras of a fusion protein between the Raf-binding domain (RBD) of Raf and GST was reduced at least 175-fold by the mutation, corresponding to a standard binding free energy decrease of at least 3.0 kcal/mol. To compute this free energy and obtain insights into the microscopic interactions favoring binding, we performed alchemical simulations of the RBD, both complexed to Ras and free in solution, in which residue 89 is gradually mutated from Arg into Lys. The simulations give a standard binding free energy decrease of 2.9+/-1.9 kcal/mol, in agreement with experiment. The use of numerous runs with three different force fields allows insights into the sources of uncertainty in the free energy and its components. The binding decreases partly because of a 7 kcal/mol higher cost to desolvate Lys upon binding, compared to Arg, due to better solvent interactions with the more concentrated Lys charge in the unbound state. This effect is expected to be general, contributing to the lower propensity of Lys to participate in protein-protein interfaces. Large contributions to the free energy change also arise from electrostatic interactions with groups up to 8 A away, namely residues 37-41 in the conserved effector domain of Ras (including 4 kcal/mol from Ser39 which loses a bifurcated hydrogen bond to Arg89), the conserved Lys84 and Lys87 of Raf, and 2-3 specific water molecules. This analysis will provide insights into the large experimental database of Ras-Raf mutations.

A streptavidin mutant with altered ligand-binding specificity

The biotin-binding site of streptavidin was modified to alter its ligand-binding specificity. In natural streptavidin, the side chains of N23 and S27 make two of the three hydrogen bonds with the ureido oxygen of biotin. These two residues were mutated to severely weaken biotin binding while attempting to maintain the affinity for two biotin analogs, 2-iminobiotin and diaminobiotin. Redesigning of the biotin-binding site used the difference in local electrostatic charge distribution between biotin and these biotin analogs. Free energy calculations predicted that the introduction of a negative charge at the position of S27 plus the mutation N23A should disrupt two of the three hydrogen bonds between natural streptavidin and the ureido oxygen of biotin. In contrast, the imino hydrogen of 2-iminobiotin should form a hydrogen bond with the side chain of an acidic amino acid at position 27. This should reduce the biotin-binding affinity by approximately eight orders of magnitude, while leaving the affinities for these biotin analogs virtually unaffected. In good agreement with these predictions, a streptavidin mutant with the N23A and S27D substitutions binds 2-iminobiotin with an affinity (Ka) of 1 x 10(6) M-1, two orders of magnitude higher than that for biotin (1 x 10(4) M-1). In contrast, the binding affinity of this streptavidin mutant for diaminobiotin (2.7 x 10(4) M-1) was lower than predicted (2.9 x 10(5) M-1), suggesting the position of the diaminobiotin in the biotin-binding site was not accurately determined by modeling.

Structural studies of binding site tryptophan mutants in the high-affinity streptavidin-biotin complex

Previous thermodynamic and computational studies have pointed to the important energetic role of aromatic contacts in generating the exceptional binding free energy of streptavidin-biotin association. We report here the crystallographic characterization of single site tryptophan mutants in investigating structural consequences of alterations in these aromatic contacts. Four tryptophan residues, Trp79, Trp92, Trp108 and Trp120, play an important role in the hydrophobic binding contributions, which along with a hydrogen bonding network and a flexible binding loop give rise to tight ligand binding (Ka approximately 10(13) M-1). The crystal structures of ligand-free and biotin-bound mutants, W79F, W108F, W120F and W120A, in the resolution range from 1.9 to 2.3 A were determined. Nine data sets for these four different mutants were collected, and structural models were refined to R-values ranging from 0.15 to 0.20. The major question addressed here is how these mutations influence the streptavidin binding site and in particular how they affect the binding mode of biotin in the complex. The overall folding of streptavidin was not significantly altered in any of the tryptophan mutants. With one exception, only minor deviations in the unbound structures were observed. In one crystal form of unbound W79F, there is a coupled shift in the side-chains of Phe29 and Tyr43 toward the mutation site, although in a different crystal form these shifts are not observed. In the bound structures, the orientation of biotin in the binding pocket was not significantly altered in the mutant complex. Compared with the wild-type streptavidin-biotin complex, there were no additional crystallographic water molecules observed for any of the mutants in the binding pocket. These structural studies thus suggest that the thermodynamic alterations can be attributed to the local alterations in binding residue composition, rather than a rearrangement of binding site architectures.

Computational carbohydrate chemistry: what theoretical methods can tell us

Computational methods have had a long history of application to carbohydrate systems and their development in this regard is discussed. The conformational analysis of carbohydrates differs in several ways from that of other biomolecules. Many glycans appear to exhibit numerous conformations coexisting in solution at room temperature and a conformational analysis of a carbohydrate must address both spatial and temporal properties. When solution nuclear magnetic resonance data are used for comparison, the simulation must give rise to ensemble-averaged properties. In contrast, when comparing to experimental data obtained from crystal structures a simulation of a crystal lattice, rather than of an isolated molecule, is appropriate. Molecular dynamics simulations are well suited for such condensed phase modeling. Interactions between carbohydrates and other biological macromolecules are also amenable to computational approaches. Having obtained a three-dimensional structure of the receptor protein, it is possible to model with accuracy the conformation of the carbohydrate in the complex. An example of the application of free energy perturbation simulations to the prediction of carbohydrate-protein binding energies is presented.

Towards structural models of molecular recognition in olfactory receptors

The G protein coupled receptors (GPCR) are an important class of proteins that act as signal transducers through the cytoplasmic membrane. Understanding the structure and activation mechanism of these proteins is crucial for understanding many different aspects of cellular signalling. The olfactory receptors correspond to the largest family of GPCRs. Very little is known about how the structures of the receptors govern the specificity of interaction which enables identification of particular odorant molecules. In this paper, we review recent developments in two areas of molecular modelling: methods for modelling the configuration of trans-membrane helices and methods for automatic docking of ligands into receptor structures. We then show how a subset of these methods can be combined to construct a model of a rat odorant receptor interacting with lyral for which experimental data are available. This modelling can help us make progress towards elucidating the specificity of interactions between receptors and odorant molecules.

Binding of biotin to streptavidin stabilizes intersubunit salt bridges between Asp61 and His87 at low pH

The remarkable stability of the streptavidin tetramer towards subunit dissociation becomes even greater upon binding of biotin. At two equivalent extensive monomer-monomer interfaces, monomers tightly associate into dimers that in turn associate into the tetramer at a less extensive dimer-dimer interface. To probe the structural basis for the enhancement of the stability of streptavidin by biotin, the crystal structures of apostreptavidin and its complexes with biotin and other small molecule and cyclic peptide ligands were determined and compared at resolutions as high as 1.36 A over a range of pH values from as low as 1.39. At low pH dramatic changes occur in the conformation and intersubunit hydrogen bonds involving the loop comprising Asp61 to Ser69. The hydrogen-bonded salt bridge between Asp61 Odelta2 and His87 Ndelta1, observed at higher pH, is replaced with a strong hydrogen bond between Asp61 Odelta1 and Asn85 Odelta1. Through crystallography at multiple pH values, the pH where this conformational change occurs, and thus the pKa of Asp61, was determined in crystals of space group I222 and/or I4122 of apostreptavidin and complexes. A range in pKa values for Asp61 was observed in these structures, the lowest being 1.78+/-0.19 for I222 streptavidin-biotin in 2.9 M (NH4)2SO4. At low pH the decrease in pKa of Asp61 and preservation of the intersubunit Asp61 Odelta2-Ndelta1 His87 hydrogen-bonded salt bridge in streptavidin-biotin versus apostreptavidin or streptavidin-peptide complexes is associated with an ordering of the flexible flap comprising residues Ala46 to Glu51, that in turn orders the Arg84 side-chain of a neighboring loop through resulting hydrogen bonds. Ordering of Arg84 in close proximity to the strong intersubunit interface appears to stabilize the conformation associated with the Asp61 Odelta2-Ndelta1 His87 hydrogen-bonded salt bridge. Thus, in addition to the established role of biotin in tetramer stabilization by direct mediation of intersubunit interactions at the weak interface through contact with Trp120, biotin may enhance tetramer stability at the strong interface more indirectly by ordering loop residues.

Constructing protein models for ligand-receptor binding thermodynamic simulations: an application to a set of peptidometic renin inhibitors

Structure-based design is the application of ligand-receptor modeling to predict the activity of a series of molecules that bind to a common receptor for which the molecular geometry is available. Successful structure-based design requires an accurate receptor model which can be economically employed in the design calculations. One goal of the work reported here has been to reduce the size of a model structure of a macromolecular receptor to allow multiple ligand-receptor molecular dynamic (MD) simulations to be computationally economical yet still provide meaningful binding thermodynamic data. A scaled-down 10 A receptor model of the enzyme renin, when subjected to an alternate atomic mass constraint, maintains the structural integrity of the composite parent crystal structure. A second goal of the work has been to develop schemes to explore and characterize the protonation states of receptors and ligand-receptor systems. Application of the charge state characterization schemes to the hydroxyethylene and statine transition state inhibitors of renin in the training set suggests a monoprotonation state of the two active-site aspartate residues, where the lone proton resides on the outer carboxylate oxygen of Asp226 is most likely. For the reduced amide transition state inhibitors an active site consisting of both aspartates in the totally ionized state, and the ligand carrying a net +1.0 charge, is most stable and consistent with experimental data.

The kinetics of protein-protein recognition

We examine a simple kinetic model for association that incorporates the basic features of protein-protein recognition within the rigid body approximation, that is, when no large conformation change occurs. Association starts with random collision at the rate k(coll) predicted by the Einstein-Smoluchowski equation. This creates an encounter pair that can evolve into a stable complex if and only if the two molecules are correctly oriented and positioned, which has a probability p(r). In the absence of long-range interactions, the bimolecular rate of association is p(r) k(coll). Long-range electrostatic interactions affect both k(coll) and p(r). The collision rate is multiplied by q(t), a factor larger than 1 when the molecules carry net charges of opposite sign as coulombic attraction makes collisions more frequent, and less than 1 in the opposite case. The probability p(r) is multiplied by a factor q(r) that represents the steering effect of electric dipoles, which preorient the molecules before they collide. The model is applied to experimental data obtained by Schreiber and Fersht (Nat. Struct. Biol. 3:427-431, 1996) on the kinetics of barnase-barstar association. When long-range electrostatic interactions are fully screened or mutated away, q(t)q(r) approximately 1, and the observed rate of productive collision is p(r) k(coll) approximately 10(5) M(-1) x s(-1). Under these conditions, p(r) approximately 1.5 x 10(-5) is determined by geometric constraints corresponding to a loss of rotational freedom. Its value is compatible with computer docking simulations and implies a rotational entropy loss deltaS(rot) approximately 22 e.u. in the transition state. At low ionic strength, long-range electrostatic interactions accelerate barnase-barstar association by a factor q(t)q(r) of up to 10(5) as favorable charge-charge and charge-dipole interactions work together to make it much faster than free diffusion would allow.

On the calculation of binding free energies using continuum methods: application to MHC class I protein-peptide interactions

This paper describes a methodology to calculate the binding free energy (delta G) of a protein-ligand complex using a continuum model of the solvent. A formal thermodynamic cycle is used to decompose the binding free energy into electrostatic and non-electrostatic contributions. In this cycle, the reactants are discharged in water, associated as purely nonpolar entities, and the final complex is then recharged. The total electrostatic free energies of the protein, the ligand, and the complex in water are calculated with the finite difference Poisson-Boltzmann (FDPB) method. The nonpolar (hydrophobic) binding free energy is calculated using a free energy-surface area relationship, with a single alkane/water surface tension coefficient (gamma aw). The loss in backbone and side-chain configurational entropy upon binding is estimated and added to the electrostatic and the nonpolar components of delta G. The methodology is applied to the binding of the murine MHC class I protein H-2Kb with three distinct peptides, and to the human MHC class I protein HLA-A2 in complex with five different peptides. Despite significant differences in the amino acid sequences of the different peptides, the experimental binding free energy differences (delta delta Gexp) are quite small (< 0.3 and < 2.7 kcal/mol for the H-2Kb and HLA-A2 complexes, respectively). For each protein, the calculations are successful in reproducing a fairly small range of values for delta delta Gcalc (< 4.4 and < 5.2 kcal/mol, respectively) although the relative peptide binding affinities of H-2Kb and HLA-A2 are not reproduced. For all protein-peptide complexes that were treated, it was found that electrostatic interactions oppose binding whereas nonpolar interactions drive complex formation. The two types of interactions appear to be correlated in that larger nonpolar contributions to binding are generally opposed by increased electrostatic contributions favoring dissociation. The factors that drive the binding of peptides to MHC proteins are discussed in light of our results.

Synthesis of conjugates for a barbiturate screening assay

Novel derivatives of barbiturates functionalized with free carboxylic acids were designed and synthesized. Coupling of 5-cyclopentyl-5-carboxycrotylbarbituric acid via its active ester to an aminofluorescein derivative produced a fluorescent tracer. Conjugation of the 5-cyclopentenyl-5-carboxyethylbarbituric acid via its mixed anhydride to thyroglobulin allowed for subsequent development of a polyclonal antibody which was evaluated for binding in a fluorescence polarization immunoassay format with various barbiturates. The binding studies showed good cross-reactivity of a variety of barbiturates containing both aromatic and aliphatic 5-substituents with the tested antisera. The relationship between the immunogen architecture, the chemical structure of the binding analytes, and the characteristics of the antisera is also presented.

Molecular dynamics study of unbinding of the avidin-biotin complex

We report molecular dynamics simulations that induce, over periods of 40-500 ps, the unbinding of biotin from avidin by means of external harmonic forces with force constants close to those of AFM cantilevers. The applied forces are sufficiently large to reduce the overall binding energy enough to yield unbinding within the measurement time. Our study complements earlier work on biotin-streptavidin that employed a much larger harmonic force constant. The simulations reveal a variety of unbinding pathways, the role of key residues contributing to adhesion as well as the spatial range over which avidin binds biotin. In contrast to the previous studies, the calculated rupture forces exceed by far those observed. We demonstrate, in the framework of models expressed in terms of one-dimensional Langevin equations with a schematic binding potential, the associated Smoluchowski equations, and the theory of first passage times, that picosecond to nanosecond simulation of ligand unbinding requires such strong forces that the resulting protein-ligand motion proceeds far from the thermally activated regime of millisecond AFM experiments, and that simulated unbinding cannot be readily extrapolated to the experimentally observed rupture.