Entropy in biological binding processes: estimation of translational entropy loss

Tetra- and Hexavalent Siglec-8 Ligands Modulate Immune Cell Activation

Carbohydrate-binding proteins are generally characterized by poor affinities for their natural glycan ligands, predominantly due to the shallow and solvent-exposed binding sites. To overcome this drawback, nature has exploited multivalency to strengthen the binding by establishing multiple interactions simultaneously. The development of oligovalent structures frequently proved to be successful, not only for proteins with multiple binding sites, but also for proteins that possess a single recognition domain. Herein we present the syntheses of a number of oligovalent ligands for Siglec-8, a monomeric I-type lectin found on eosinophils and mast cells, alongside the thermodynamic characterization of their binding. While the enthalpic contribution of each binding epitope was within a narrow range to that of the monomeric ligand, the entropy penalty increased steadily with growing valency. Additionally, we observed a successful agonistic binding of the tetra- and hexavalent and, to an even larger extent, multivalent ligands to Siglec-8 on immune cells and modulation of immune cell activation. Thus, triggering a biological effect is not restricted to multivalent ligands but could be induced by low oligovalent ligands as well, whereas a monovalent ligand, despite binding with similar affinity, showed an antagonistic effect.

Nuclear localisation sequences of chloride intracellular channels 1 and 4 facilitate nuclear import via interactions with import mediator importin-α: An empirical and theoretical perspective

Chloride intracellular channel proteins (CLICs) display ubiquitous expression, with each member exhibiting specific subcellular localisation. While all CLICs, except CLIC3, exhibit a highly conserved putative nuclear localisation sequence (NLS), only CLIC1, CLIC3 and CLIC4 exist within the nucleus. The CLIC4 NLS, 199-KVVAKKYR-206, appears crucial for nuclear entry and interacts with mouse nuclear import mediator Impα isoform 1, omitting the IBB domain (mImpα1ΔIBB). The essential nature of the basic residues in the CLIC4 NLS has been established by the fact that mutating out these residues inhibits nuclear import, which in turn is linked to cutaneous squamous cell cancer. Given the conservation of the CLIC NLS, CLIC1 likely follows a similar import pathway to CLIC4. Peptides of the CLIC1 (Pep1; Pep1_S C/S mutant) and CLIC4 (Pep4) NLSs were designed to examine binding to human Impα isoform 1, omitting the IBB domain (hImpα1ΔIBB). Molecular docking indicated that the core CLIC NLS region (KKYR) forms a similar binding pattern to both mImpα1ΔIBB and hImpα1ΔIBB. Fluorescence quenching demonstrated that Pep1_S (K_d  ≈ 237 μM) and Pep4 (K_d  ≈ 317 μM) bind hImpα1ΔIBB weakly. Isothermal titration calorimetry confirmed the weak binding interaction between Pep4 and hImpα1ΔIBB (K_d  ≈ 130 μM) and the presence of a proton-linked effect. This weak interaction may be due to regions distal from the CLIC NLS needed to stabilise and strengthen hImpα1ΔIBB binding. Additionally, this NLS may preferentially bind another hImpα isoform with different flexibility properties.

Handling the Hurdles on the Way to Anti-tuberculosis Drug Development

The global epidemic of tuberculosis (TB) imposes a sustained epidemiologic vigilance and investments in research by governments. Mycobacterium tuberculosis, the main causative agent of TB in human beings, is a very successful pathogen, being the main cause of death in the population among infectious agents. In 2018, ~10 million individuals were contaminated with this bacillus and became ill with TB, and about 1.2 million succumbed to the disease. Most of the success of the M. tuberculosis to linger in the population comes from its ability to persist in an asymptomatic latent state into the host and, in fact, the majority of the individuals are unaware of being contaminated. Even though TB is a treatable disease and is curable in most cases, the treatment is lengthy and laborious. In addition, the rise of resistance to first-line anti-TB drugs elicits a response from TB research groups to discover new chemical entities, preferably with novel mechanisms of action. The pathway to find a new TB drug, however, is arduous and has many barriers that are difficult to overcome. Fortunately, several approaches are available today to be pursued by scientists interested in anti-TB drug development, which goes from massively testing chemical compounds against mycobacteria, to discovering new molecular targets by genetic manipulation. This review presents some difficulties found along the TB drug development process and illustrates different approaches that might be used to try to identify new molecules or targets that are able to impair M. tuberculosis survival.

Ligand Binding Rate Constants in Heme Proteins Using Markov State Models and Molecular Dynamics Simulations

Computer simulation studies of the molecular basis for ligand migration in proteins allow the description of key events such as the transition between docking sites, displacement of existing ligands and solvent molecules, and open/closure of specific "gates", among others. In heme proteins, ligand migration from the solvent to the active site preludes the binding to the heme iron and triggers different functions. In this work, molecular dynamics simulations, a Markov State Model of migration and empirical kinetic equations are combined to study the migration of O₂ and NO in two truncated hemoglobins of Mycobacterium tuberculosis (Mt-TrHbN and Mt-TrHbO). For Mt-TrHbN, we show that the difference in the association constant in the oxy and deoxy states relies mainly in the displacement of water molecules anchored in the distal cavity in the deoxy form. The results here provide a valuable approach to study ligand migration in globins.

Functional Role of Solvent Entropy and Conformational Entropy of Metal Binding in a Dynamically Driven Allosteric System

Allostery is a regulatory phenomenon whereby ligand binding to one site influences the binding of the same or a different ligand to another site on a macromolecule. The physical origins of allosteric regulation remain under intense investigation. In general terms, ligand-induced structural changes, perturbations of residue-specific dynamics, and surrounding solvent molecules all potentially contribute to the global energetics of allostery. While the role of solvent is generally well understood in regulatory events associated with major protein structural rearrangements, the degree to which protein dynamics impact solvent degrees of freedom is unclear, particularly in cases of dynamically driven allostery. With the aid of new crystal structures, extensive calorimetric and residue-specific dynamics studies over a range of time scales and temperatures, we dissect for the first time the relative degree to which changes in solvent entropy and residue-specific dynamics impact dynamically driven, allosteric inhibition of DNA binding by Zn in the zinc efflux repressor, CzrA (chromosomal zinc-regulated repressor). We show that non-native residue-specific dynamics in allosterically impaired CzrA mutants are accompanied by significant perturbations in solvent entropy that cannot be predicted from crystal structures. We conclude that functional dynamics are not necessarily restricted to protein residues but involve surface water molecules that may be responding to ligand (Zn)-mediated perturbations in protein internal motions that define the conformational ensemble, rather than major structural rearrangements.

Conformational switch of the bacterial adhesin FimH in the absence of the regulatory domain: Engineering a minimalistic allosteric system

For many biological processes such as ligand binding, enzymatic catalysis, or protein folding, allosteric regulation of protein conformation and dynamics is fundamentally important. One example is the bacterial adhesin FimH, where the C-terminal pilin domain exerts negative allosteric control over binding of the N-terminal lectin domain to mannosylated ligands on host cells. When the lectin and pilin domains are separated under shear stress, the FimH-ligand interaction switches in a so-called catch-bond mechanism from the low- to high-affinity state. So far, it has been assumed that the pilin domain is essential for the allosteric propagation within the lectin domain that would otherwise be conformationally rigid. To test this hypothesis, we generated mutants of the isolated FimH lectin domain and characterized their thermodynamic, kinetic, and structural properties using isothermal titration calorimetry, surface plasmon resonance, nuclear magnetic resonance, and X-ray techniques. Intriguingly, some of the mutants mimicked the conformational and kinetic behaviors of the full-length protein and, even in absence of the pilin domain, conducted the cross-talk between allosteric sites and the mannoside-binding pocket. Thus, these mutants represent a minimalistic allosteric system of FimH, useful for further mechanistic studies and antagonist design.

The price of flexibility - a case study on septanoses as pyranose mimetics

Seven-membered ring mimetics of mannose were studied as ligands for the mannose-specific bacterial lectin FimH, which plays an essential role in the first step of urinary tract infections (UTI). A competitive binding assay and isothermal titration calorimetry (ITC) experiments indicated an approximately ten-fold lower affinity for the seven-membered ring mannose mimetic 2-O-n-heptyl-1,6-anhydro-d-glycero-d-galactitol (7) compared to n-heptyl α-d-mannopyranoside (2), resulting exclusively from a loss of conformational entropy. Investigations by solution NMR, X-ray crystallography, and molecular modeling revealed that 7 establishes a superimposable H-bond network compared to mannoside 2, but at the price of a high entropic penalty due to the loss of its pronounced conformational flexibility. These results underscore the importance of having access to the complete thermodynamic profile of a molecular interaction to "rescue" ligands from entropic penalties with an otherwise perfect fit to the protein binding site.

Structural Stability of the Coiled-Coil Domain of Tumor Susceptibility Gene (TSG)-101

The tumor susceptibility gene-101 coiled coil domain (TSG101cc) is an integral component of the endosomal maturation machinery and cytokinesis, and also interacts with several transcription factors. The TSG101cc has been crystallized as a homotetramer but is known to interact with two of its binding partners as a heterotrimer. To investigate this apparent discrepancy, we examined the solution thermodynamics of the TSG101cc. Here, we use circular dichroism, differential scanning calorimetry, analytical ultracentrifugation, fluorescence, and structural thermodynamic analysis to investigate the structural stability and the unfolding of the TSG101cc. We demonstrate that TSG101cc exists in solution primarily as a tetramer, which unfolds in a two-state manner. Surprisingly, no homodimeric or homotrimeric species were detected. Structural thermodynamic analysis of the homotetrameric structure and comparison with known oligomeric coiled-coils suggests that the TSG101cc homotetramer is comparatively unstable on a per residue basis. Furthermore, the homotrimeric coiled-coil is predicted to be much less stable than the functional heterotrimeric coiled-coil in the endosomal sorting complex required for transport 1 (ESCRT1). These results support a model whereby the tetramer–monomer equilibrium of TSG101 serves as the cellular reservoir of TSG101, which is effectively outcompeted when its binding partners are present and the heteroternary complex can form.

Bacterial protease uses distinct thermodynamic signatures for substrate recognition

Porphyromonas gingivalis and Porphyromonas endodontalis are important bacteria related to periodontitis, the most common chronic inflammatory disease in humans worldwide. Its comorbidity with systemic diseases, such as type 2 diabetes, oral cancers and cardiovascular diseases, continues to generate considerable interest. Surprisingly, these two microorganisms do not ferment carbohydrates; rather they use proteinaceous substrates as carbon and energy sources. However, the underlying biochemical mechanisms of their energy metabolism remain unknown. Here, we show that dipeptidyl peptidase 11 (DPP11), a central metabolic enzyme in these bacteria, undergoes a conformational change upon peptide binding to distinguish substrates from end products. It binds substrates through an entropy-driven process and end products in an enthalpy-driven fashion. We show that increase in protein conformational entropy is the main-driving force for substrate binding via the unfolding of specific regions of the enzyme (“entropy reservoirs”). The relationship between our structural and thermodynamics data yields a distinct model for protein-protein interactions where protein conformational entropy modulates the binding free-energy. Further, our findings provide a framework for the structure-based design of specific DPP11 inhibitors.

Application of a BOSS-Gaussian interface for QM/MM simulations of Henry and methyl transfer reactions

Hybrid quantum mechanics and molecular mechanics (QM/MM) computer simulations have become an indispensable tool for studying chemical and biological phenomena for systems too large to treat with QM alone. For several decades, semiempirical QM methods have been used in QM/MM simulations. However, with increased computational resources, the introduction of ab initio and density function methods into on-the-fly QM/MM simulations is being increasingly preferred. This adaptation can be accomplished with a program interface that tethers independent QM and MM software packages. This report introduces such an interface for the BOSS and Gaussian programs, featuring modification of BOSS to request QM energies and partial atomic charges from Gaussian. A customizable C-shell linker script facilitates the interprogram communication. The BOSS-Gaussian interface also provides convenient access to Charge Model 5 (CM5) partial atomic charges for multiple purposes including QM/MM studies of reactions. In this report, the BOSS-Gaussian interface is applied to a nitroaldol (Henry) reaction and two methyl transfer reactions in aqueous solution. Improved agreement with experiment is found by determining free-energy surfaces with MP2/CM5 QM/MM simulations than previously reported investigations using semiempirical methods.

Free Energy Diagram for the Heterogeneous Enzymatic Hydrolysis of Glycosidic Bonds in Cellulose

Kinetic and thermodynamic data have been analyzed according to transition state theory and a simplified reaction scheme for the enzymatic hydrolysis of insoluble cellulose. For the cellobiohydrolase Cel7A from Hypocrea jecorina (Trichoderma reesei), we were able to measure or collect relevant values for all stable and activated complexes defined by the reaction scheme and hence propose a free energy diagram for the full heterogeneous process. For other Cel7A enzymes, including variants with and without carbohydrate binding module (CBM), we obtained activation parameters for the association and dissociation of the enzyme-substrate complex. The results showed that the kinetics of enzyme-substrate association (i.e. formation of the Michaelis complex) was almost entirely entropy-controlled and that the activation entropy corresponded approximately to the loss of translational and rotational degrees of freedom of the dissolved enzyme. This implied that the transition state occurred early in the path where the enzyme has lost these degrees of freedom but not yet established extensive contact interactions in the binding tunnel. For dissociation, a similar analysis suggested that the transition state was late in the path where most enzyme-substrate contacts were broken. Activation enthalpies revealed that the rate of dissociation was far more temperature-sensitive than the rates of both association and the inner catalytic cycle. Comparisons of one- and two-domain variants showed that the CBM had no influence on the transition state for association but increased the free energy barrier for dissociation. Hence, the CBM appeared to promote the stability of the complex by delaying dissociation rather than accelerating association.

A universal entropy-driven mechanism for thioredoxin-target recognition

Cysteine residues in cytosolic proteins are maintained in their reduced state, but can undergo oxidation owing to posttranslational modification during redox signaling or under conditions of oxidative stress. In large part, the reduction of oxidized protein cysteines is mediated by a small 12-kDa thiol oxidoreductase, thioredoxin (Trx). Trx provides reducing equivalents for central metabolic enzymes and is implicated in redox regulation of a wide number of target proteins, including transcription factors. Despite its importance in cellular redox homeostasis, the precise mechanism by which Trx recognizes target proteins, especially in the absence of any apparent signature binding sequence or motif, remains unknown. Knowledge of the forces associated with the molecular recognition that governs Trx-protein interactions is fundamental to our understanding of target specificity. To gain insight into Trx-target recognition, we have thermodynamically characterized the noncovalent interactions between Trx and target proteins before S-S reduction using isothermal titration calorimetry (ITC). Our findings indicate that Trx recognizes the oxidized form of its target proteins with exquisite selectivity, compared with their reduced counterparts. Furthermore, we show that recognition is dependent on the conformational restriction inherent to oxidized targets. Significantly, the thermodynamic signatures for multiple Trx targets reveal favorable entropic contributions as the major recognition force dictating these protein-protein interactions. Taken together, our data afford significant new insight into the molecular forces responsible for Trx-target recognition and should aid the design of new strategies for thiol oxidoreductase inhibition.

Single water entropy: hydrophobic crossover and application to drug binding

Entropy of water plays an important role in both chemical and biological processes e.g. hydrophobic effect, molecular recognition etc. Here we use a new approach to calculate translational and rotational entropy of the individual water molecules around different hydrophobic and charged solutes. We show that for small hydrophobic solutes, the translational and rotational entropies of each water molecule increase as a function of its distance from the solute reaching finally to a constant bulk value. As the size of the solute increases (0.746 nm), the behavior of the translational entropy is opposite; water molecules closest to the solute have higher entropy that reduces with distance from the solute. This indicates that there is a crossover in translational entropy of water molecules around hydrophobic solutes from negative to positive values as the size of the solute is increased. Rotational entropy of water molecules around hydrophobic solutes for all sizes increases with distance from the solute, indicating the absence of crossover in rotational entropy. This makes the crossover in total entropy (translation + rotation) of water molecule happen at much larger size (>1.5 nm) for hydrophobic solutes. Translational entropy of single water molecule scales logarithmically (Str(QH) = C + kB ln V), with the volume V obtained from the ellipsoid of inertia. We further discuss the origin of higher entropy of water around water and show the possibility of recovering the entropy loss of some hypothetical solutes. The results obtained are helpful to understand water entropy behavior around various hydrophobic and charged environments within biomolecules. Finally, we show how our approach can be used to calculate the entropy of the individual water molecules in a protein cavity that may be replaced during ligand binding.

Thermodynamics of the GTP-GDP-operated conformational switch of selenocysteine-specific translation factor SelB

SelB is a specialized translation factor that binds GTP and GDP and delivers selenocysteyl-tRNA (Sec-tRNA(Sec)) to the ribosome. By analogy to elongation factor Tu (EF-Tu), SelB is expected to control the delivery and release of Sec-tRNA(Sec) to the ribosome by the structural switch between GTP- and GDP-bound conformations. However, crystal structures of SelB suggested a similar domain arrangement in the apo form and GDP- and GTP-bound forms of the factor, raising the question of how SelB can fulfill its delivery function. Here, we studied the thermodynamics of guanine nucleotide binding to SelB by isothermal titration calorimetry in the temperature range between 10 and 25 °C using GTP, GDP, and two nonhydrolyzable GTP analogs, guanosine 5'-O-(γ-thio)triphosphate (GTPγS) and guanosine 5'-(β,γ-imido)-triphosphate (GDPNP). The binding of SelB to either guanine nucleotide is characterized by a large heat capacity change (-621, -467, -235, and -275 cal × mol(-1) × K(-1), with GTP, GTPγS, GDPNP, and GDP, respectively), associated with compensatory changes in binding entropy and enthalpy. Changes in heat capacity indicate a large decrease of the solvent-accessible surface area in SelB, amounting to 43 or 32 amino acids buried upon binding of GTP or GTPγS, respectively, and 15-19 amino acids upon binding GDP or GDPNP. The similarity of the GTP and GDP forms in the crystal structures can be attributed to the use of GDPNP, which appears to induce a structure of SelB that is more similar to the GDP than to the GTP-bound form.

FimH antagonists: structure-activity and structure-property relationships for biphenyl α-D-mannopyranosides

Urinary tract infections (UTIs) are caused primarily by uropathogenic Escherichia coli (UPEC), which encode filamentous surface-adhesive organelles called type 1 pili. FimH is located at the tips of these pili. The initial attachment of UPEC to host cells is mediated by the interaction of the carbohydrate recognition domain (CRD) of FimH with oligomannosides on urothelial cells. Blocking these lectins with carbohydrates or analogues thereof prevents bacterial adhesion to host cells and therefore offers a potential therapeutic approach for prevention and/or treatment of UTIs. Although numerous FimH antagonists have been developed so far, few of them meet the requirement for clinical application due to poor pharmacokinetics. Additionally, the binding mode of an antagonist to the CRD of FimH can switch from an in-docking mode to an out-docking mode, depending on the structure of the antagonist. In this communication, biphenyl α-D-mannosides were modified to improve their binding affinity, to explore their binding mode, and to optimize their pharmacokinetic properties. The inhibitory potential of the FimH antagonists was measured in a cell-free competitive binding assay, a cell-based flow cytometry assay, and by isothermal titration calorimetry. Furthermore, pharmacokinetic properties such as log D, solubility, and membrane permeation were analyzed. As a result, a structure-activity and structure-property relationships were established for a series of biphenyl α-D-mannosides.

On the role of flexibility in protein-ligand interactions: the example of p53 tetramerization domain

The recognition of protein surfaces by designed ligands has become an attractive approach in drug discovery. However, the variable nature and irregular behavior of protein surfaces defy this new area of research. The easy to understand "lock-and-key" model is far from being the ideal paradigm in biomolecular interactions and, hence, any new finding on how proteins and ligands behave in recognition events paves a step of the way. Herein, we illustrate a clear example on how an increase in flexibility of both protein and ligand can result in an increase in the stability of the macromolecular complex. The biophysical study of the interaction between a designed flexible tetraguanidinium-calix[4]arene and the tetramerization domain of protein p53 (p53TD) and its natural mutant p53TD-R337H shows how the floppy mutant domain interacts more tightly with the ligand than the well-packed wild-type protein. Moreover, the flexible calixarene ligand interacts with higher affinity to both wild-type and mutated protein domains than a conformationally rigid calixarene analog previously reported. These findings underscore the crucial role of flexibility in molecular recognition processes, for both small ligands and large biomolecular surfaces.

Calculation of configurational entropy with a Boltzmann-quasiharmonic model: the origin of high-affinity protein-ligand binding

Accurate assessment of configurational entropy remains a large challenge in biology. While many methods exist to calculate configurational entropy, there is a balance between accuracy and computational demands. Here we calculate ligand and protein conformational entropies using the Boltzmann-quasiharmonic (BQH) method, which treats the first-order entropy term by the Boltzmann expression for entropy while determining correlations using the quasiharmonic model. This method is tested by comparison with the exact Clausius expression for entropy on a range of test molecules ranging from small ligands to a protein. Using the BQH method, we then analyze the rotational and translational (R/T) entropy change upon ligand binding for five protein complexes to explore the origins of extremely tight affinity. The results suggest that in these systems such affinity is achieved by a combination of simultaneously maintaining good protein-ligand contacts while allowing significant residual R/T motion of the ligand through suitable protein motions.

High affinity antigen recognition of the dual specific variants of herceptin is entropy-driven in spite of structural plasticity

The antigen-binding site of Herceptin, an anti-human Epidermal Growth Factor Receptor 2 (HER2) antibody, was engineered to add a second specificity toward Vascular Endothelial Growth Factor (VEGF) to create a high affinity two-in-one antibody bH1. Crystal structures of bH1 in complex with either antigen showed that, in comparison to Herceptin, this antibody exhibited greater conformational variability, also called "structural plasticity". Here, we analyzed the biophysical and thermodynamic properties of the dual specific variants of Herceptin to understand how a single antibody binds two unrelated protein antigens. We showed that while bH1 and the affinity-improved bH1-44, in particular, maintained many properties of Herceptin including binding affinity, kinetics and the use of residues for antigen recognition, they differed in the binding thermodynamics. The interactions of bH1 and its variants with both antigens were characterized by large favorable entropy changes whereas the Herceptin/HER2 interaction involved a large favorable enthalpy change. By dissecting the total entropy change and the energy barrier for dual interaction, we determined that the significant structural plasticity of the bH1 antibodies demanded by the dual specificity did not translate into the expected increase of entropic penalty relative to Herceptin. Clearly, dual antigen recognition of the Herceptin variants involves divergent antibody conformations of nearly equivalent energetic states. Hence, increasing the structural plasticity of an antigen-binding site without increasing the entropic cost may play a role for antibodies to evolve multi-specificity. Our report represents the first comprehensive biophysical analysis of a high affinity dual specific antibody binding two unrelated protein antigens, furthering our understanding of the thermodynamics that drive the vast antigen recognition capacity of the antibody repertoire.

SH2 domains recognize contextual peptide sequence information to determine selectivity

Selective ligand recognition by modular protein interaction domains is a primary determinant of specificity in signaling pathways. Src homology 2 (SH2) domains fulfill this capacity immediately downstream of tyrosine kinases, acting to recruit their host polypeptides to ligand proteins harboring phosphorylated tyrosine residues. The degree to which SH2 domains are selective and the mechanisms underlying selectivity are fundamental to understanding phosphotyrosine signaling networks. An examination of interactions between 50 SH2 domains and a set of 192 phosphotyrosine peptides corresponding to physiological motifs within FGF, insulin, and IGF-1 receptor pathways indicates that individual SH2 domains have distinct recognition properties and exhibit a remarkable degree of selectivity beyond that predicted by previously described binding motifs. The underlying basis for such selectivity is the ability of SH2 domains to recognize both permissive amino acid residues that enhance binding and non-permissive amino acid residues that oppose binding in the vicinity of the essential phosphotyrosine. Neighboring positions affect one another so local sequence context matters to SH2 domains. This complex linguistics allows SH2 domains to distinguish subtle differences in peptide ligands. This newly appreciated contextual dependence substantially increases the accessible information content embedded in the peptide ligands that can be effectively integrated to determine binding. This concept may serve more broadly as a paradigm for subtle recognition of physiological ligands by protein interaction domains.

In vivo distribution of polymeric nanoparticles at the whole-body, tumor, and cellular levels

PURPOSE

Block copolymer micelles (BCMs) were functionalized with indium-111 and/or epidermal growth factor (EGF), which enabled investigation of the in vivo transport of passively and actively targeted BCMs. The integration of conventional and image-based techniques afforded novel quantitative means to achieve an in-depth insight into the fate of polymeric nanoparticles in vivo.

METHODS

Pharmacokinetics and biodistribution studies were performed in athymic mice bearing human breast xenografts to evaluate the whole-body transport of NT-BCMs (non-targeted, EGF-) and T-BCMs (targeted, EGF+). The intratumoral distribution of BCMs was investigated using MicroSPECT/CT and autoradiographic imaging, complemented with quantitative MATLAB® analyses. Tumors were fractionated for quantifying intracellular uptake of BCMs via γ-counting.

RESULTS

The intratumoral distribution of NT-BCMs and T-BCMs were found to be heterogeneous, and positively correlated with tumor vascularization (r>0.68 ± 0.04). The enhanced in vivo cell uptake and cell membrane binding of T-BCMs were found to delay their clearance from tumors overexpressing EGFR, and therefore resulted in enhanced tumor accumulation for the T-BCMs in comparison to the NT-BCMs.

CONCLUSIONS

Adequate passive targeting is required in order to achieve effective active targeting. Tumor physiology has a significant impact on the transvascular and intratumoral transport of passively and actively targeted BCMs.

Using multi-objective computational design to extend protein promiscuity

Many enzymes possess, besides their native function, additional promiscuous activities. Proteins with several activities (multipurpose catalysts) may have a wide range of biotechnological and biomedical applications. Natural promiscuity, however, appears to be of limited scope in this context, because the latent (promiscuous) function is often related to the evolved one (sharing the active site and even the chemical mechanism) and its enhancement upon suitable mutations usually brings about a decrease in the native activity. Here we explore the use of computational protein design to overcome these limitations. The high-plasticity positions close to the original ("native") active-site are the most promising candidates for mutations that create a second active-site associated to a new function. To avoid compromising protein folding and native activity, we propose a minimal-perturbation approach based on the combinatorial optimization of, both the de novo catalytic activity and the folding free-energy: essentially, we construct the Pareto Set of optimal stability/promiscuous-function solutions. We validate our approach by introducing a promiscuous esterase activity in E. coli thioredoxin on the basis of mutations at positions close to the native-active-site disulfide-bridge. Native oxidoreductase activity is not compromised and it is, in fact, found to be 1.5-fold enhanced, as determined by an insulin-reduction assay. This work provides general guidelines as to how computational design can be used to expand the scope and applications of protein promiscuity. From a more general viewpoint, it illustrates the potential of multi-objective optimization as the computational analogue of multi-feature natural selection.

Thermodynamic analysis of ligand binding and ligand binding-induced tertiary structure formation by the thiamine pyrophosphate riboswitch

The thi-box riboswitch regulates gene expression in response to the intracellular concentration of thiamine pyrophosphate (TPP) in archaea, bacteria, and eukarya. To complement previous biochemical, genetic, and structural studies of this phylogenetically widespread RNA domain, we have characterized its interaction with TPP by isothermal titration calorimetry. This shows that TPP binding is highly dependent on Mg(2+) concentration. The dissociation constant decreases from approximately 200 nM at 0.5 mM Mg(2+) concentration to approximately 9 nM at 2.5 mM Mg(2+) concentration. Binding is enthalpically driven, but the unfavorable entropy of binding decreases as Mg(2+) concentration rises, suggesting that divalent cations serve to pre-organize the RNA. Mutagenesis, biochemical analysis, and a new crystal structure of the riboswitch suggest that a critical element that participates in organizing the riboswitch structure is the tertiary interaction formed between the P3 and L5 regions. This tertiary contact is distant from the TPP binding site, but calorimetric analysis reveals that even subtle mutations in L5 can have readily detectable effects on TPP binding. The thermodynamic signatures of these mutations, namely decreased favorable enthalpy of binding and small effects on entropy of binding, are consistent with the P3-L5 association contributing allosterically to TPP-induced compaction of the RNA.

Thermodynamics of Protein–Ligand Interactions: History, Presence, and Future Aspects

The understanding of molecular recognition processes of small ligands and biological macromolecules requires a complete characterization of the binding energetics and correlation of thermodynamic data with interacting structures involved. A quantitative description of the forces that govern molecular associations requires determination of changes of all thermodynamic parameters, including free energy of binding (deltaG), enthalpy (deltaH), and entropy (deltaS) of binding and the heat capacity change (deltaCp). A close insight into the binding process is of significant and practical interest, since it provides the fundamental know-how for development of structure-based molecular design-strategies. The only direct method to measure the heat change during complex formation at constant temperature is provided by isothermal titration calorimetry (ITC). With this method one binding partner is titrated into a solution containing the interaction partner, thereby generating or absorbing heat. This heat is the direct observable that can be quantified by the calorimeter. The use of ITC has been limited due to the lack of sensitivity, but recent developments in instrument design permit to measure heat effects generated by nanomol (typically 10-100) amounts of reactants. ITC has emerged as the primary tool for characterizing interactions in terms of thermodynamic parameters. Because heat changes occur in almost all chemical and biochemical processes, ITC can be used for numerous applications, e.g., binding studies of antibody-antigen, protein-peptide, protein-protein, enzyme-inhibitor or enzyme-substrate, carbohydrate-protein, DNA-protein (and many more) interactions as well as enzyme kinetics. Under appropriate conditions data analysis from a single experiment yields deltaH, K(B), the stoichiometry (n), deltaG and deltaS of binding. Moreover, ITC experiments performed at different temperatures yield the heat capacity change (deltaCp). The informational content of thermodynamic data is large, and it has been shown that it plays an important role in the elucidation of binding mechanisms and, through the link to structural data, also in rational drug design. In this review we will present a comprehensive overview to ITC by giving some historical background to calorimetry, outline some critical experimental and data analysis aspects, discuss the latest developments, and give three recent examples of studies published with respect to macromolecule-ligand interactions that have utilized ITC technology.

Entropic cost of protein-ligand binding and its dependence on the entropy in solution

Two theoretical formulations are proposed and compared for the loss of translational and rotational entropy upon protein-ligand binding in water. The two theories share the same approach to evaluate the translational and rotational entropy of the ligand when bound. The potential of the bound ligand is modeled by six harmonic oscillators that are parametrized from the force and torque magnitudes measured in a molecular dynamics simulation, yielding vibrational and librational entropies. In the aqueous phase, the theories differ because there is no unique way to assign the total entropy to molecules in solution. In one approach, the ligand is allowed unrestricted access to the full solution volume at the standard concentration and is assigned the same translational and rotational entropy as if it were an ideal gas. We term this a "molecule-frame" (MF) theory because it considers configurational space in the reference frame of the molecule of interest. The entropy of the solvent is penalized because it is excluded from the molecule's volume. In the second theory, all molecules including the solvent are confined by their neighbors in mean-field configurational volumes. This we term a "system-frame" (SF) theory because the configurational space available to all molecules is considered in the reference frame of the whole system. Molecules have vibrational and librational entropy in the same way as they do when bound. In addition, the discrete size of the solvent molecules quantizes the configurational space into an effective number of minima according to the solute molecule's standard concentration and the mean volume of a solvent molecule. This leads to the cratic entropy expressed in terms of the solute molecule's mole fraction. The equivalent number of minima in rotational space depends on both the solute molecule's volume and the solvent molecule's volume. This leads to an equation for the orientational entropy based on the proposed concept of "angle fraction". The MF and SF theories are applied to calculate the translational and rotational entropy losses involved in the formation of six different protein-ligand complexes, in two of which the ligand is water. The MF entropy losses range from -80 to -142 J K(-1) mol(-1) for ligands at the 1 M standard-state concentration and from -52 to -63 J K(-1) mol(-1) for water at the 55.6 M standard-state concentration. They depend logarithmically on both the number and strength of interactions between the ligand and protein through the forces and torques. This is observed to lead to moderate dependencies on the ligands' moments of inertia and masses. The SF entropy losses are smaller and range from -50 to -75 J K(-1) mol(-1) for ligands at the 1 M standard-state concentration and from 0 to -12 J K(-1) mol(-1) for water. They depend logarithmically on the ligand solvent's molecular volume and weakly on the relative strengths of the ligand's interactions with the protein and water. The cratic entropy loss in water at the standard concentration is constant and is also demonstrated to be implicit in MF theories. Entropy losses from the two approaches are also compared with those from other computational approaches and with experiment. The use of the force and torque magnitudes leads to smaller bound volumes than are obtained from ligand-displacement approaches. The general agreement of the SF entropy losses with those from experiment suggests that the SF theory is more consistent with the assumptions made in experimental measurements than the MF solvation theories, which would require a compensating entropy gain in the solvent in order to agree.

Characterizing the role of ensemble modulation in mutation-induced changes in binding affinity

Protein conformational fluctuations are key contributors to biological function, mediating important processes such as enzyme catalysis, molecular recognition, and allosteric signaling. To better understand the role of conformational fluctuations in substrate/ligand recognition, we analyzed, experimentally and computationally, the binding reaction between an SH3 domain and the recognition peptide of its partner protein. The fluctuations in this SH3 domain were enumerated by using an algorithm based on the hard sphere collision model, and the binding energetics resulting from these fluctuations were calculated using a structure-based energy function parametrized to solvent accessible surface areas. Surprisingly, this simple model reproduced the effects of mutations on the experimentally determined SH3 binding energetics, within the uncertainties of the measurements, indicating that conformational fluctuations in SH3, and in particular the RT loop region, are structurally diverse and are well-approximated by the randomly configured states. The mutated positions in SH3 were distant to the binding site and involved Ala and Gly substitutions of solvent exposed positions in the RT loop. To characterize these fluctuations, we applied principal coordinate analysis to the computed ensembles, uncovering the principal modes of conformational variation. It is shown that the observed differences in binding affinity between each mutant, and thus the apparent coupling between the mutated sites, can be described in terms of the changes in these principal modes. These results indicate that dynamic loops in proteins can populate a broad conformational ensemble and that a quantitative understanding of molecular recognition requires consideration of the entire distribution of states.

The structure of the anti-c-myc antibody 9E10 Fab fragment/epitope peptide complex reveals a novel binding mode dominated by the heavy chain hypervariable loops

The X-ray structure of the Fab fragment from the anti-c-myc antibody 9E10 was determined both as complex with its epitope peptide and for the free Fab. In the complex, two Fab molecules adopt an unusual head to head orientation with the epitope peptide arranged between them. In contrast, the free Fab forms a dimer with different orientation. In the Fab/peptide complex the peptide is bound to one of the two Fabs at the "back" of its extended CDR H3, in a cleft with CDR H1, thus forming a short, three-stranded antiparallel beta-sheet. The N- and C-terminal parts of the peptide are also in contact with the neighboring Fab fragment. Comparison between the CDR H3s of the two Fab molecules in complex with the peptide and those from the free Fab reveals high flexibility of this loop. This structural feature is in line with thermodynamic data from isothermic titration calorimetry.

The solution structure of BMPR-IA reveals a local disorder-to-order transition upon BMP-2 binding

The structure of the extracellular domain of BMP receptor IA was determined in solution by NMR spectroscopy and compared to its structure when bound to its ligand BMP-2. While most parts of the secondary structure are highly conserved between the bound and unbound forms, large conformational rearrangements can be observed in the beta4beta5 loop of BMPR-IA, which is in contact with BMP-2 and harbors the main binding determinants for the BMPR-IA-BMP-2 interaction. In its unbound form, helix alpha1 in BMPR-IA, which is in the center of the binding epitope for BMP-2, is missing. Since BMP-2 also shows conformational changes in the type I receptor epitope upon binding to BMPR-IA, both binding partners pass through an induced fit mechanism to adapt their binding interfaces to a given interaction surface. The inherent flexibility of both partners possibly explains the promiscuous ligand-receptor interaction observed in the BMP protein superfamily.

Structural basis of the differential binding of the SH3 domains of Grb2 adaptor to the guanine nucleotide exchange factor Sos1

Grb2-Sos1 interaction, mediated by the canonical binding of N-terminal SH3 (nSH3) and C-terminal SH3 (cSH3) domains of Grb2 to a proline-rich sequence in Sos1, provides a key regulatory switch that relays signaling from activated receptor tyrosine kinases to downstream effector molecules such as Ras. Here, using isothermal titration calorimetry in combination with site-directed mutagenesis, we show that the nSH3 domain binds to a Sos1-derived peptide containing the proline-rich consensus motif PPVPPR with an affinity that is nearly threefold greater than that observed for the binding of cSH3 domain. We further demonstrate that such differential binding of nSH3 domain relative to the cSH3 domain is largely due to the requirement of a specific acidic residue in the RT loop of the beta-barrel fold to engage in the formation of a salt bridge with the arginine residue in the consensus motif PPVPPR. While this role is fulfilled by an optimally positioned D15 in the nSH3 domain, the chemically distinct and structurally non-equivalent E171 substitutes in the case of the cSH3 domain. Additionally, our data suggest that salt tightly modulates the binding of both SH3 domains to Sos1 in a thermodynamically distinct manner. Our data further reveal that, while binding of both SH3 domains to Sos1 is under enthalpic control, the nSH3 binding suffers from entropic penalty in contrast to entropic gain accompanying the binding of cSH3, implying that the two domains employ differential thermodynamic mechanisms for Sos1 recognition. Our new findings are rationalized in the context of 3D structural models of SH3 domains in complex with the Sos1 peptide. Taken together, our study provides structural basis of the differential binding of SH3 domains of Grb2 to Sos1 and a detailed thermodynamic profile of this key protein-protein interaction pertinent to cellular signaling and cancer.

Exploring the energy landscape of antibody-antigen complexes: protein dynamics, flexibility, and molecular recognition

The production of antibodies that selectively bind virtually any foreign compound is the hallmark of the immune system. While much is understood about how sequence diversity contributes to this remarkable feat of molecular recognition, little is known about how sequence diversity impacts antibody dynamics, which is also expected to contribute to molecular recognition. Toward this goal, we examined a panel of antibodies elicited to the chromophoric antigen fluorescein. On the basis of isothermal titration calorimetry, we selected six antibodies that bind fluorescein with diverse binding entropies, suggestive of varying contributions of dynamics to molecular recognition. Sequencing revealed that two pairs of antibodies employ homologous heavy chains that were derived from common germline genes, while the other two heavy chains and all six of the light chains were derived from different germline genes and are not homologous. Interestingly, more than half of all the somatic mutations acquired during affinity maturation among the six antibodies are located in positions unlikely to contact fluorescein directly. To quantify and compare the dynamics of the antibody-fluorescein complexes, three-pulse photon echo peak shift and transient grating spectroscopy were employed. All of the antibodies exhibited motions on three distinct time scales, ultrafast motions on the <100 fs time scale, diffusive motions on the picosecond time scale, and motions that occur on time scales longer than nanoseconds and thus appear static. However, the exact frequency of the picosecond time scale motion and the relative contribution of the different motions vary significantly among the antibody-chromophore complexes, revealing a high level of dynamic diversity. Using a hierarchical model, we relate the data to features of the antibodies' energy landscapes as well as their flexibility in terms of elasticity and plasticity. In all, the data provide a consistent picture of antibody flexibility, which interestingly appears to be correlated with binding entropy as well as with germline gene use and the mutations introduced during affinity maturation. The data also provide a gauge of the dynamic diversity of the antibody repertoire and suggest that this diversity might contribute to molecular recognition by facilitating the recognition of the broadest range of foreign molecules.

Coupling of folding and DNA-binding in the bZIP domains of Jun-Fos heterodimeric transcription factor

In response to mitogenic stimuli, the heterodimeric transcription factor Jun-Fos binds to the promoters of a diverse array of genes involved in critical cellular responses such as cell growth and proliferation, cell cycle regulation, embryogenic development and cancer. In so doing, Jun-Fos heterodimer regulates gene expression central to physiology and pathology of the cell in a specific and timely manner. Here, using the technique of isothermal titration calorimetry (ITC), we report detailed thermodynamics of the bZIP domains of Jun-Fos heterodimer to synthetic dsDNA oligos containing the TRE and CRE consensus promoter elements. Our data suggest that binding of the bZIP domains to both TRE and CRE is under enthalpic control and accompanied by entropic penalty at physiological temperatures. Although the bZIP domains bind to both TRE and CRE with very similar affinities, the enthalpic contributions to the free energy of binding to CRE are more favorable than TRE, while the entropic penalty to the free energy of binding to TRE is smaller than CRE. Despite such differences in their thermodynamic signatures, enthalpy and entropy of binding of the bZIP domains to both TRE and CRE are highly temperature-dependent and largely compensate each other resulting in negligible effect of temperature on the free energy of binding. From the plot of enthalpy change versus temperature, the magnitude of heat capacity change determined is much larger than that expected from the direct association of bZIP domains with DNA. This observation is interpreted to suggest that the basic regions in the bZIP domains are largely unstructured in the absence of DNA and only become structured upon interaction with DNA in a coupled folding and binding manner. Our new findings are rationalized in the context of 3D structural models of bZIP domains of Jun-Fos heterodimer in complex with dsDNA oligos containing the TRE and CRE consensus sequences. Taken together, our study demonstrates that enthalpy is the major driving force for a key protein-DNA interaction pertinent to cellular signaling and that protein-DNA interactions with similar binding affinities may be accompanied by differential thermodynamic signatures. Our data corroborate the notion that the DNA-induced protein structural changes are a general feature of the bZIP family of transcription factors.

Modeling fatty acid delivery from intestinal fatty acid binding protein to a membrane

Intestinal fatty acid binding protein (IFABP) interacts with biological membranes and delivers fatty acid (FA) into them via a collisional mechanism. However, the membrane-bound structure of the protein and the pathway of FA transfer are not precisely known. We used molecular dynamics (MD) simulations with an implicit membrane model to determine the optimal orientation of apo- and holo-IFABP (bound with palmitate) on an anionic membrane. In this orientation, the helical portal region, delimited by the alphaII helix and the betaC-betaD and betaE-betaF turns, is oriented toward the membrane whereas the putative beta-strand portal, delimited by the betaB-betaC, betaF-betaG, betaH-betaI turns and the N terminus, is exposed to solvent. Starting from the MD structure of holo-IFABP in the optimal orientation relative to the membrane, we examined the release of palmitate via both pathways. Although the domains can widen enough to allow the passage of palmitate, fatty acid release through the helical portal region incurs smaller conformational changes and a lower energetic cost.

The contribution of conformational adjustments and long-range electrostatic forces to the CD2/CD58 interaction

CD2 is a T cell surface molecule that enhances T and natural killer cell function by binding its ligands CD58 (humans) and CD48 (rodents) on antigen-presenting or target cells. Here we show that the CD2/CD58 interaction is enthalpically driven and accompanied by unfavorable entropic changes. Taken together with structural studies, this indicates that binding is accompanied by energetically significant conformational adjustments. Despite having a highly charged binding interface, neither the affinity nor the rate constants of the CD2/CD58 interaction were affected by changes in ionic strength, indicating that long-range electrostatic forces make no net contribution to binding.

Global changes in local protein dynamics reduce the entropic cost of carbohydrate binding in the arabinose-binding protein

Protein dynamics make important but poorly understood contributions to molecular recognition phenomena. To address this, we measure changes in fast protein dynamics that accompany the interaction of the arabinose-binding protein (ABP) with its ligand, d-galactose, using NMR relaxation and molecular dynamics simulation. These two approaches present an entirely consistent view of the dynamic changes that occur in the protein backbone upon ligand binding. Increases in the amplitude of motions are observed throughout the protein, with the exception of a few residues in the binding site, which show restriction of dynamics. These counter-intuitive results imply that a localised binding event causes a global increase in the extent of protein dynamics on the pico- to nanosecond timescale. This global dynamic change constitutes a substantial favourable entropic contribution to the free energy of ligand binding. These results suggest that the structure and dynamics of ABP may be adapted to exploit dynamic changes to reduce the entropic costs of binding.

Hydrophobic interactions are the driving force for the binding of peptide mimotopes and Staphylococcal protein A to recombinant human IgG1

We studied the interaction of several nona-peptide mimotopes of different sequence and Staphylococcal protein A (SpA) with a recombinant human IgG1 antibody using isothermal titration calorimetry (ITC). The amino acid primary structure of the peptides was varied in order to identify the specific antibody-peptide binding sites. Additionally, the influence of temperature and salt concentration was investigated. An attempt was made to elucidate the structural changes upon complex formation using the determined thermodynamic parameters. The amino acid composition of the mimotopes determined their binding affinity. The binding constant K (a) of the mimotopes was in the range 1 x 10(4) to 1 x 10(6) M(-1). The binding constant of SpA was on the average about three orders of magnitude higher than that of the peptides. The binding constant of the peptides and of SpA decreased with temperature and the binding process was connected with negative changes in enthalpy, entropy, and heat capacity. The binding of the mimotopes to the Fab part of the IgG1 antibody and binding of SpA to the Fc part of the IgG1 antibody were mainly driven by hydrophobic effects and associated with a relatively large change in water-accessible surface area. Determinants for a strong/reduced antibody-peptide binding were identified.

Energetics in correlation with structural features: the case of micellization

Understanding micellization processes at the molecular level has direct relevance for biological self-assembly, folding, and association processes. As such, it requires complete characterization of the micellization thermodynamics, including its correlation with the corresponding structural features. In this context, micellization of a series of model non-ionic surfactants (poly(ethylene glycol) monooctyl ethers, C(8)E(gamma)) was studied by isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC). The corresponding structural properties of C(8)E(gamma) micelles were investigated by small-angle X-ray scattering (SAXS). The C(8)E(gamma) micellization, characterized independently from ITC, DSC, and structural data, reveals that deltaH(M)(o) > 0, deltaS(M)(o) > 0, and deltaC(P)(M)(o) < 0, while the dissection of its energetics shows that it is primarily governed by the transfer of 20-30 C(8) alkyl chains from aqueous solution into the nonpolar core (r approximately 1.3 nm) of the spherical micelle. Moreover, thermodynamic parameters of micellization, estimated from the structural features related to the changes in solvent-accessible surface areas upon micellization, are in a good agreement with the corresponding parameters obtained from the analysis of ITC and DSC data. We have shown that the contributions to deltaS(M)(o) other than from hydration (deltaS(M)(other)(o)), estimated from experimental data, appear to be small (deltaS(M)(other)(o) < 0.1 deltaS(M)(other)(o)) and agree well with the theoretical estimates expressed as a sum of the corresponding translational, conformational, and size contributions. These deltaS(M)(other)(o) contributions are much less unfavorable than those estimated for a rigid-body association, which indicates the dynamic nature of the C(8)E(gamma) micellar aggregates. the dynamic nature of the C8EY micellar aggregates.

Temperature-dependent variations of ligand-receptor contact points in hAT1

Photoaffinity labelling is regularly used to investigate proteins, including peptidergic G protein-coupled receptors (GPCR). To this purpose benzophenone photolabels have been widely used to identify many contact residues in ligand-binding pockets. The three-dimensional binding environment of the human angiotensin II type 1 receptor hAT(1) has been determined using an iterative methionine mutagenesis strategy based on the photochemical properties and preferential incorporation of benzophenone onto methionine. This has led to the construction of a ligand-bound receptor structure. The present study investigated the effect of temperature on the accessibility of some of these contact points. The hAT(1) receptor and two representative Met mutants (H256M-hAT(1) and F293M-hAT(1)) from the iterative mutagenesis study were photolabelled with the benzophenone-ligand (125)I-[Sar(1), Bpa(8)]AngII at temperatures ranging from - 15 degrees C to 37 degrees C. Labelled receptors were partially purified and digested with cyanogen bromide to identify the contact points or segments. There were no changes in receptor contacts or labelling in the 7th transmembrane domains (TMD) of hAT(1) and F293M-hAT(1) across the temperature range. However, a temperature-dependent change in the ligand-receptor contact of H256M-hAT(1) was observed. At - 15 degrees C, H256M labelling was identical to that of hAT(1), indicating that the interaction was specific to the 7th TMD. Significant labelling changes were observed at higher temperatures and at 37 degrees C labelling occurred almost exclusively at mutated residue H256M-hAT(1) in the 6th TMD. Simultaneous competitive labelling of different areas of this target protein indicated that the ligand-receptor structure became increasingly fluctual at physiological temperatures, while a more compact, low mobility, and low energy conformation prevailed at low temperatures.

Rigidification of a flexible protease inhibitor variant upon binding to trypsin

The Tyr35-->Gly replacement in bovine pancreatic trypsin inhibitor (BPTI) has previously been shown to dramatically enhance the flexibility of the trypsin-binding region of the free inhibitor and to destabilize the interaction with the protease by about 3 kcal/mol. The effects of this replacement on the enzyme-inhibitor interaction were further studied here by X-ray crystallography and isothermal titration calorimetry (ITC). The co-crystal structure of Y35G BPTI bound to trypsin was determined using 1.65 A resolution X-ray diffraction data collected from cryopreserved crystals, and a new structure of the complex with wild-type BPTI under the same conditions was determined using 1.62 A data. These structures reveal that, in contrast to the free protein, Y35G BPTI adopts a conformation nearly identical with that of the wild-type protein, with a water-filled cavity in place of the missing Tyr side-chain. The crystallographic temperature factors for the two complexes indicate that the mutant inhibitor is nearly as rigid as the wild-type protein when bound to trypsin. Calorimetric measurements show that the change in enthalpy upon dissociation of the complex is 2.5 kcal/mol less favorable for the complex containing Y35G BPTI than for the complex with the wild-type inhibitor. Thus, the destabilization of the complex resulting from the Y35G replacement is due to a more favorable change in entropy upon dissociation. The heat capacity changes for dissociation of the mutant and wild-type complexes were very similar, suggesting that the entropic effects probably do not arise from solvation effects, but are more likely due to an increase in protein conformational entropy upon dissociation of the mutant inhibitor. These results define the biophysical role of a highly conserved core residue located outside of a protein-binding interface, demonstrating that Tyr35 has little impact on the trypsin-bound BPTI structure and acts primarily to define the structure of the free protein so as to maximize binding affinity.

Calculations of solute and solvent entropies from molecular dynamics simulations

The translational, rotational and conformational (vibrational) entropy contributions to ligand-receptor binding free energies are analyzed within the standard formulation of statistical thermodynamics. It is shown that the partitioning of the binding entropy into different components is to some extent arbitrary, but an appropriate method to calculate both translational and rotational entropy contributions to noncovalent association is by estimating the configurational volumes of the ligand in the bound and free states. Different approaches to calculating solute entropies using free energy perturbation calculations, configurational volumes based on root-mean-square fluctuations and covariance matrix based quasiharmonic analysis are illustrated for some simple molecular systems. Numerical examples for the different contributions demonstrate that theoretically derived results are well reproduced by the approximations. Calculation of solvent entropies, either using total potential energy averages or van't Hoff plots, are carried out for the case of ion solvation in water. Although convergence problems will persist for large and complex simulation systems, good agreement with experiment is obtained here for relative and absolute ion hydration entropies. We also outline how solvent and solute entropic contributions are taken into account in empirical binding free energy calculations using the linear interaction energy method. In particular it is shown that empirical scaling of the nonpolar intermolecular ligand interaction energy effectively takes into account size dependent contributions to the binding free energy.