Thermodynamic databases for proteins and protein-nucleic acid interactions

Thermodynamic data regarding proteins and their interactions are important for understanding the mechanisms of protein folding, protein stability, and molecular recognition. Although there are several structural databases available for proteins and their complexes with other molecules, databases for experimental thermodynamic data on protein stability and interactions are rather scarce. Thus, we have developed two electronically accessible thermodynamic databases. ProTherm, Thermodynamic Database for Proteins and Mutants, contains numerical data of several thermodynamic parameters of protein stability, experimental methods and conditions, along with structural, functional, and literature information. ProNIT, Thermodynamic Database for Protein-Nucleic Acid Interactions, contains thermodynamic data for protein-nucleic acid binding, experimental conditions, structural information of proteins, nucleic acids and the complex, and literature information. These data have been incorporated into 3DinSight, an integrated database for structure, function, and properties of biomolecules. A WWW interface allows users to search for data based on various conditions, with different display and sorting options, and to visualize molecular structures and their interactions. These thermodynamic databases, together with structural databases, help researchers gain insight into the relationship among structure, function, and thermodynamics of proteins and their interactions, and will become useful resources for studying proteins in the postgenomic era.     

Molecular recognition and information gain

Molecular recognition is considered as an information transduction process, and information theory is applied to analyze the efficiency of the information transduction during molecular recognition. It is shown that a certain amount of information is gained in the recognition process, and it can be expressed in terms of thermodynamic quantities in the binding process of molecules. The information gain, which can be interpreted as the amount of information extracted from the sequence or structure by a molecule, represents the intrinsic ability of a molecule to recognize specific sequence or structure out of large ambient ensemble based on physical interactions. In contrast with thermodynamic quantities themselves, the information gain is a normalized quantity, and thus serves as a good measure of specificity in the molecular recognition. The information gain, which can be evaluated experimentally, enables us to identify the specific interactions and compare the specificity among different molecules quantitatively. In contrast with macroscopic recognition, specificity of molecular recognition may be severely restricted by the thermodynamic nature of molecular interactions. The thermodynamic meaning of the information gain and its biological implications to molecular recognition are discussed.     

Thermodynamic database supports deciphering protein–nucleic acid interactions

The thermodynamics of protein–nucleic acid interactions (PNIs) is crucial for elucidating the mechanisms of molecular recognition and pathological consequences. The Protein–Nucleic Acid Thermodynamics Database (PNATDB) is a database containing experimentally determined thermodynamic parameters along with sequence, structural, and function data, which is available free online.     

Mass spectrometry‐based methods to study protein architecture and dynamics

Mass spectrometry is now an indispensable tool in the armamentarium of molecular biophysics, where it is used for tasks ranging from protein sequencing and mapping of post‐translational modifications to studies of higher order structure, conformational dynamics, and interactions of proteins with small molecule ligands and other biopolymers. This mini‐review highlights several popular mass spectrometry‐based tools that are now commonly used for structural studies of proteins beyond their covalent structure with a particular emphasis on hydrogen exchange and direct electrospray ionization mass spectrometry.     

The thermodynamics of protein interactions with essential first row transition metals

BackgroundThe binding of metal ions to proteins is a crucial process required for their catalytic activity, structural stability and/or functional regulation. Isothermal titration calorimetry provides a wealth of fundamental information which when combined with structural data allow for a much deeper understanding of the underlying molecular mechanism.Scope of reviewA rigorous understanding of any molecular interaction requires in part an in-depth quantification of its thermodynamic properties. Here, we provide an overview of recent studies that have used ITC to quantify the interaction of essential first row transition metals with relevant proteins and highlight major findings from these thermodynamic studies.General significanceThe thermodynamic characterization of metal ion–protein interactions is one important step to understanding the role that metal ions play in living systems. Such characterization has important implications not only to elucidating proteins' structure-function relationships and biological properties but also in the biotechnology sector, medicine and drug design particularly since a number of metal ions are involved in several neurodegenerative diseases.Major conclusionsIsothermal titration calorimetry measurements can provide complete thermodynamic profiles of any molecular interaction through the simultaneous determination of the reaction binding stoichiometry, binding affinity as well as the enthalpic and entropic contributions to the free energy change thus enabling a more in-depth understanding of the nature of these interactions. This article is part of a Special Issue entitled Microcalorimetry in the BioSciences — Principles and Applications, edited by Fadi Bou-Abdallah.     

A molecular thermodynamic view of DNA-drug interactions: a case study of 25 minor-groove binders

Developing a molecular view of the thermodynamics of DNA recognition is essential to the design of ligands for regulating gene expression. In a first comprehensive attempt at sketching an atlas of DNA-drug energetics, we present here a detailed thermodynamic view of minor-groove recognition by small molecules via a computational study on 25 DNA-drug complexes. The studies are configured in the MMGBSA (Molecular Mechanics-Generalized Born-Solvent Accessibility) framework at the current state of the art and facilitate a structure-energy component correlation. Analyses were conducted on both energy minimized structures of DNA-drug complexes and molecular dynamics trajectories developed for the purpose of this study. While highlighting the favorable role of packing, shape complementarity, and van der Waals and hydrophobic interactions of the drugs in the minor groove in conformity with experiment, the studies reveal an interesting annihilation of favorable electrostatics by desolvation. Structural modifications attempted on the ligands point to the requisite physico-chemical factors for obtaining improved binding energies. Hydrogen bonds predicted to be important for specificity based on structural considerations do not always turn out to be significant to binding in post facto analyses of molecular dynamics trajectories, which treat thermal averaging, solvent, and counterion effects rigorously. The strength of the hydrogen bonds retained between the DNA and drug during the molecular dynamics simulations is approximately 1kcal/mol. Overall, the study reveals the compensatory nature of the diverse binding free energy components, possible threshold limits for some of these properties, and the availability of a computationally viable free energy methodology which could be of value in drug-design endeavors.     

Modeling protein recognition of carbohydrates

We have a limited understanding of the details of molecular recognition of carbohydrates by proteins, which is critical to a multitude of biological processes. Furthermore, carbohydrate-modifying proteins such as glycosyl hydrolases and phosphorylases are of growing importance as potential drug targets. Interactions between proteins and carbohydrates have complex thermodynamics, and in general the specific positioning of only a few hydroxyl groups determines their binding affinities. A thorough understanding of both carbohydrate and protein structures is thus essential to predict these interactions. An atomic-level view of carbohydrate recognition through structures of carbohydrate-active enzymes complexed with transition-state inhibitors reveals some of the distinctive molecular features unique to protein-carbohydrate complexes. However, the inherent flexibility of carbohydrates and their often water-mediated hydrogen bonding to proteins makes simulation of their complexes difficult. Nonetheless, recent developments such as the parameterization of specific force fields and docking scoring functions have greatly improved our ability to predict protein-carbohydrate interactions. We review protein-carbohydrate complexes having defined molecular requirements for specific carbohydrate recognition by proteins, providing an overview of the different computational techniques available to model them.     

Conformational Dynamics and Cooperativity Drive the Specificity of a Protein-Ligand Interaction

Molecular recognition is critical for the fidelity of signal transduction in biology. Conversely, the disruption of protein-protein interactions can lead to disease. Thus, comprehension of the molecular determinants of specificity is essential for understanding normal biological signaling processes and for the development of precise therapeutics. Although high-resolution structures have provided atomic details of molecular interactions, much less is known about the influence of cooperativity and conformational dynamics. Here, we used the Tiam2 PSD-95/Dlg/ZO-1 (PDZ) domain and a quadruple mutant (QM), engineered by swapping the identity of four residues important for specificity in the Tiam1 PDZ into the Tiam2 PDZ domain, as a model system to investigate the role of cooperativity and dynamics in PDZ ligand specificity. Surprisingly, equilibrium binding experiments found that the ligand specificity of the Tiam2 QM was switched to that of the Tiam1 PDZ. NMR-based studies indicated that Tiam2 QM PDZ, but not other mutants, had extensive microsecond to millisecond motions distributed throughout the entire domain suggesting structural cooperativity between the mutated residues. Thermodynamic analyses revealed energetic cooperativity between residues in distinct specificity subpockets that was dependent upon the identity of the ligand, indicating a context-dependent binding mechanism. Finally, isothermal titration calorimetry experiments showed distinct entropic signatures along the mutational trajectory from the Tiam2 wild-type to the QM PDZ domain. Collectively, our studies provide unique insights into how structure, conformational dynamics, and thermodynamics combine to modulate ligand-binding specificity and have implications for the evolution, regulation, and design of protein-ligand interactions.     

Ligand-mediated dimerization of a carbohydrate-binding molecule reveals a novel mechanism for protein-carbohydrate recognition

The structural and thermodynamic basis for carbohydrate-protein recognition is of considerable importance. NCP-1, which is a component of the Piromyces equi cellulase/hemicellulase complex, presents a provocative model for analyzing how structural and mutational changes can influence the ligand specificity of carbohydrate-binding proteins. NCP-1 contains two "family 29" carbohydrate-binding modules designated CBM29-1 and CBM29-2, respectively, that display unusually broad specificity; the proteins interact weakly with xylan, exhibit moderate affinity for cellulose and mannan, and bind tightly to the beta-1,4-linked glucose-mannose heteropolymer glucomannan. The crystal structure of CBM29-2 in complex with cellohexaose and mannohexaose identified key residues involved in ligand recognition. By exploiting this structural information and the broad specificity of CBM29-2, we have used this protein as a template to explore the evolutionary mechanisms that can lead to significant changes in ligand specificity. Here, we report the properties of the E78R mutant of CBM29-2, which displays ligand specificity that is different from that of wild-type CBM29-2; the protein retains significant affinity for cellulose but does not bind to mannan or glucomannan. Significantly, E78R exhibits a stoichiometry of 0.5 when binding to cellohexaose, and both calorimetry and ultracentrifugation show that the mutant protein displays ligand-mediated dimerization in solution. The three-dimensional structure of E78R in complex with cellohexaose reveals the intriguing molecular basis for this "dimeric" binding mode that involves the lamination of the oligosaccharide between two CBM molecules. The 2-fold screw axis of the ligand is mirrored in the orientation of the two protein domains with adjacent sugar rings stacking against the equivalent aromatic residues in the binding site of each protein molecule of the molecular sandwich. The sandwiching of an oligosaccharide chain between two protein modules, leading to ligand-induced formation of the binding site, represents a completely novel mechanism for protein-carbohydrate recognition that may mimic that displayed by naturally dimeric protein-carbohydrate interactions.     

Multiple molecular recognition mechanisms. Cytochrome P450--a case study

Biomolecular recognition is complex. The balance between the different molecular properties that contribute to molecular recognition, such as shape, electrostatics, dynamics and entropy, varies from case to case. This, along with the extent of experimental characterization, influences the choice of appropriate computational approaches to study biomolecular interactions. Here, we present computational studies of cytochrome P450 enzymes and their interactions with small molecules and with other proteins. These interactions exemplify some of the diversity of molecular determinants of binding affinity and specificity observed for proteins and we discuss some of the challenges that they pose for molecular modelling and simulation.     

Thermodynamics of T-cell receptor-peptide/MHC interactions: progress and opportunities

alphabeta T-cell receptors (TCRs) recognize peptide antigens presented by class I or class II major histocompatibility complex molecules (pMHC). Here we review the use of thermodynamic measurements in the study of TCR-pMHC interactions, with attention to the diversity in binding thermodynamics and how this is related to the variation in TCR-pMHC interfaces. We show that there is no enthalpic or entropic signature for TCR binding; rather, enthalpy and entropy changes vary in a compensatory manner that reflects a narrow free energy window for the interactions that have been characterized. Binding enthalpy and entropy changes do not correlate with structural features such as buried surface area or the number of hydrogen bonds within TCR-pMHC interfaces, possibly reflecting the myriad of contributors to binding thermodynamics, but likely also reflecting a reliance on van't Hoff over calorimetric measurements and the unaccounted influence of equilibria linked to binding. TCR-pMHC binding heat capacity changes likewise vary considerably. In some cases, the heat capacity changes are consistent with conformational differences between bound and free receptors, but there is little data indicating these conformational differences represent the need to organize disordered CDR loops. In this regard, we discuss how thermodynamics may provide additional insight into conformational changes occurring upon TCR binding. Finally, we highlight opportunities for the further use of thermodynamic measurements in the study of TCR-pMHC interactions, not only for understanding TCR binding in general, but also for understanding specifics of individual interactions and the engineering of TCRs with desired molecular recognition properties.     

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.     

Protein homeostasis disorders of key enzymes of amino acids metabolism: mutation-induced protein kinetic destabilization and new therapeutic strategies

Many inborn errors of amino acids metabolism are caused by single point mutations affecting the ability of proteins to fold properly (i.e., protein homeostasis), thus leading to enzyme loss-of-function. Mutations may affect protein homeostasis by altering intrinsic physical properties of the polypeptide (folding thermodynamics, and rates of folding/unfolding/misfolding) as well as the interaction of partially folded states with elements of the protein homeostasis network (such as molecular chaperones and proteolytic machineries). Understanding these mutational effects on protein homeostasis is required to develop new therapeutic strategies aimed to target specific features of the mutant polypeptide. Here, I review recent work in three different diseases of protein homeostasis associated to inborn errors of amino acids metabolism: phenylketonuria, inherited homocystinuria and primary hyperoxaluria type I. These three different genetic disorders involve proteins operating in different cell organelles and displaying different structural complexities. Mutations often decrease protein kinetic stability of the native state (i.e., its half-life for irreversible denaturation), which can be studied using simple kinetic models amenable to biophysical and biochemical characterization. Natural ligands and pharmacological chaperones are shown to stabilize mutant enzymes, thus supporting their therapeutic application to overcome protein kinetic destabilization. The role of molecular chaperones in protein folding and misfolding is also discussed as well as their potential pharmacological modulation as promising new therapeutic approaches. Since current available treatments for these diseases are either burdening or only successful in a fraction of patients, alternative treatments must be considered covering studies from protein structure and biophysics to studies in animal models and patients.     

Specific interactions for ab initio folding of protein terminal regions with secondary structures

Proteins fold into unique three-dimensional structures by specific, orientation-dependent interactions between amino acid residues. Here, we extract orientation-dependent interactions from protein structures by treating each polar atom as a dipole with a direction. The resulting statistical energy function successfully refolds 13 out of 16 fully unfolded secondary-structure terminal regions of 10-23 amino acid residues in 15 small proteins. Dissecting the orientation-dependent energy function reveals that the orientation preference between hydrogen-bonded atoms is not enough to account for the structural specificity of proteins. The result has significant implications on the theoretical and experimental searches for specific interactions involved in protein folding and molecular recognition between proteins and other biologically active molecules.     

Structural and functional analysis of the PDZ domains of human HtrA1 and HtrA3

High-temperature requirement A (HtrA) and its homologs contain a serine protease domain followed by one or two PDZ domains. Bacterial HtrA proteins and the mitochondrial protein HtrA2/Omi maintain cell function by acting as both molecular chaperones and proteases to manage misfolded proteins. The biological roles of the mammalian family members HtrA1 and HtrA3 are less clear. We report a detailed structural and functional analysis of the PDZ domains of human HtrA1 and HtrA3 using peptide libraries and affinity assays to define specificity, structural studies to view the molecular details of ligand recognition, and alanine scanning mutagenesis to investigate the energetic contributions of individual residues to ligand binding. In common with HtrA2/Omi, we show that the PDZ domains of HtrA1 and HtrA3 recognize hydrophobic polypeptides, and while C-terminal sequences are preferred, internal sequences are also recognized. However, the details of the interactions differ, as different domains rely on interactions with different residues within the ligand to achieve high affinity binding. The results suggest that mammalian HtrA PDZ domains interact with a broad range of hydrophobic binding partners. This promiscuous specificity resembles that of bacterial HtrA family members and suggests a similar function for recognizing misfolded polypeptides with exposed hydrophobic sequences. Our results support a common activation mechanism for the HtrA family, whereby hydrophobic peptides bind to the PDZ domain and induce conformational changes that activate the protease. Such a mechanism is well suited to proteases evolved for the recognition and degradation of misfolded proteins.     

Mutational analysis of a sequence-specific ssDNA binding lupus autoantibody

11F8 is a murine anti-ssDNA monoclonal autoantibody isolated from a lupus prone autoimmune mouse. This mAb binds sequence specifically, and prior studies have defined the thermodynamic and kinetic basis for sequence-specific recognition of ssDNA (Ackroyd, P. C., et al. (2001) Biochemistry 40, 2911-2922; Beckingham, J. A. and Glick, G. D. (2001) Bioorg. Med. Chem. 9, 2243-2252). Here we present experiments designed to identify the residues on 11F8 that mediate sequence-specific, noncognate, and nonspecific recognition of ssDNA and their contribution to the overall binding thermodynamics. Site-directed mutagenesis of an 11F8 single-chain construct reveals that six residues within the complementarity determining regions of 11F8 account for ca. 80% of the binding free energy and that there is little cooperativity between these residues. Germline-encoded aromatic and hydrophobic side chains provides the basis for nonspecific recognition of single-stranded thymine nucleobases. Sequence-specific recognition is controlled by a tyrosine in the heavy chain along with a somatically mutated arginine residue. Our data show that the manner in which 11F8 achieves sequence-specific recognition more closely resembles RNA-binding proteins such as U1A than other types of nucleic acid binding proteins. In addition, comparing the primary sequence of 11F8 with clonally related antibodies that differ by less than five amino acids suggests that somatic mutations which confer sequence specificity may be a feature that distinguishes glomerulotrophic pathogenic anti-DNA from those that are benign.     

On the Importance of Polar Interactions for Complexes Containing Intrinsically Disordered Proteins

There is a growing recognition for the importance of proteins with large intrinsically disordered (ID) segments in cell signaling and regulation. ID segments in these proteins often harbor regions that mediate molecular recognition. Coupled folding and binding of the recognition regions has been proposed to confer high specificity to interactions involving ID segments. However, researchers recently questioned the origin of the interaction specificity of ID proteins because of the overrepresentation of hydrophobic residues in their interaction interfaces. Here, we focused on the role of polar and charged residues in interactions mediated by ID segments. Making use of the extended nature of most ID segments when in complex with globular proteins, we first identified large numbers of complexes between globular proteins and ID segments by using radius-of-gyration-based selection criteria. Consistent with previous studies, we found the interfaces of these complexes to be enriched in hydrophobic residues, and that these residues contribute significantly to the stability of the interaction interface. However, our analyses also show that polar interactions play a larger role in these complexes than in structured protein complexes. Computational alanine scanning and salt-bridge analysis indicate that interfaces in ID complexes are highly complementary with respect to electrostatics, more so than interfaces of globular proteins. Follow-up calculations of the electrostatic contributions to the free energy of binding uncovered significantly stronger Coulombic interactions in complexes harbouring ID segments than in structured protein complexes. However, they are counter-balanced by even higher polar-desolvation penalties. We propose that polar interactions are a key contributing factor to the observed high specificity of ID segment-mediated interactions.     

Thermodynamics of drug-DNA interactions

Many anticancer, antibiotic, and antiviral drugs exert their primary biological effects by reversibly interacting with nucleic acids. Therefore, these biomolecules represent a major target in drug development strategies designed to produce next generation therapeutics for diseases such as cancer. In order to improve the clinical efficacy of existing drugs and also to design new ones it is necessary to understand the molecular basis of drug-DNA interactions in structural, thermodynamic, and kinetic detail. The past decade has witnessed an increase in the number of rigorous biophysical studies of drug-DNA systems and considerable knowledge has been gained in the energetics of these binding reactions. This is, in part, due to the increased availability of high-sensitivity calorimetric techniques, which have allowed the thermodynamics of drug-DNA interactions to be probed directly and accurately. The focus of this article is to review thermodynamic approaches to examining drug-DNA recognition. Specifically, an overview of a recently developed method of analysis that dissects the binding free energy of these reactions into five component terms is presented. The results of applying this analysis to the DNA binding interactions of both minor groove drugs and intercalators are discussed. The solvent water plays a key role in nucleic acid structure and consequently in the binding of ligands to these biomolecules. Any rational approach to DNA-targeted drug design requires an understanding of how water participates in recognition and binding events. Recent studies examining hydration changes that accompany DNA binding by intercalators will be reviewed. Finally some aspects of cooperativity in drug-DNA interactions are described and the importance of considering cooperative effects when examining these reactions is highlighted.