Site-Specific Thermodynamics: Understanding Cooperativity in Molecular Recognition

Biomolecular dynamics with machine-learned quantum-mechanical force fields trained on diverse chemical fragments

The GEMS method enables molecular dynamics simulations of large heterogeneous systems at ab initio quality.

Epitopes, paratopes, and other topes 30 years on: Understanding what we are talking about

The question of which protein antigens, such as HLA class I or class II molecules, will bind, and how well, to a given antibody is often assumed to depend exclusively on the details of protein surface structure. These structures are usually based on static models resulting from X-ray crystallography. While these notions are useful, the ultimate causal factors determining how well a given antigen binds a given antibody are based in thermodynamics and can include atomic mobility and the time-varying conformations of proteins. In this article, fundamental biophysical principles of antibody-antigen interaction are discussed, concepts critical for a deeper understanding of the pertinent molecular phenomena are highlighted, and common misunderstandings are identified and debunked.

Phase Transitions of Associative Biomacromolecules

Multivalent proteins and nucleic acids, collectively referred to as multivalent associative biomacromolecules, provide the driving forces for the formation and compositional regulation of biomolecular condensates. Here, we review the key concepts of phase transitions of aqueous solutions of associative biomacromolecules, specifically proteins that include folded domains and intrinsically disordered regions. The phase transitions of these systems come under the rubric of coupled associative and segregative transitions. The concepts underlying these processes are presented, and their relevance to biomolecular condensates is discussed.

Cation-π Interactions in Graphene-Containing Systems for Water Treatment and Beyond

Cation-π interactions are common in nature, especially in organisms. Their profound influences in chemistry, physics, and biology have been continuously investigated since they were discovered in 1981. However, the importance of cation-π interactions in materials science, regarding carbonaceous nanomaterials, has just been realized. The interplay between cations and delocalized polarizable π electrons of graphene would bring about significant changes to the intrinsic characteristics of graphene and greatly affect the device performance based on graphene and its derivatives. Here, the cation-π interactions in graphene containing systems for water treatment applications (e.g., separation membranes, adsorbents) are highlighted. The cross-linking effects caused by cation-π interactions contribute to membrane stability and selectivity and enhanced adsorption. Their roles in dominating the performance of graphene-based structures for other specific applications are also discussed. Relevant theoretical modeling and calculations are summarized to offer an in-depth understanding of the underlying mechanisms which can help in designing more functional materials and structures. Perspectives on the potential directions that deserve effort are also presented.

Cooperativity effect of the ππ interaction between drug and DNA on intercalative binding induced by H-bonds: a QM/QTAIM investigation of the curcuminadenineH2O model system

In order to reveal the nature of intercalative binding of drug to DNA, the cooperativity effect of the ππ interaction was investigated in the curcuminadenineH2O model system by applying a combined QM and QTAIM computational approach. The H-bonds between the electron-donating group of curcumin and adenine induce the formation of the ππ stacking. The introduction of H2O weakens the H-bonding and ππ interactions, leading to an anti-cooperativity effect, as is confirmed by the AIM (atoms in molecules) and RDG (reduced density gradient) analysis. Thus, it can be inferred that the anti-cooperative effect is the main driving force for the intercalative binding of drug to DNA bases, which is in agreement with many experimental phenomena. Therefore, the designed DNA-targeted intercalating drugs should possess not only hydrophobic moieties, but also strong electron-donating groups bound to the DNA bases with H-bonds, which can slow the variation rates of the strengths of the H-bonding and ππ interactions between drug and DNA bases in the anti-cooperative process, leading to the intercalation formation. The enthalpy change is the major factor driving the positive thermodynamic cooperativity.

Complexity in pH-Dependent Ribozyme Kinetics: Dark pKa Shifts and Wavy Rate-pH Profiles

Charged bases occur in RNA enzymes, or ribozymes, where they play key roles in catalysis. Cationic bases donate protons and perform electrostatic catalysis, while anionic bases accept protons. We previously published simulations of rate-pH profiles for ribozymes in terms of species plots for the general acid and general base that have been useful for understanding how ribozymes respond to pH. In that study, we did not consider interaction between the general acid and general base or interaction with other species on the RNA. Since that report, diverse small ribozyme classes have been discovered, many of which have charged nucleobases or metal ions in the active site that can either directly interact and participate in catalysis or indirectly interact as "influencers". Herein, we simulate experimental rate-pH profiles in terms of species plots in which reverse protonated charged nucleobases interact. These analyses uncover two surprising features of pH-dependent enzyme kinetics. (1) Cooperativity between the general acid and general base enhances population of the functional forms of a ribozyme and manifests itself as hidden or "dark" pK_a shifts, real pK_a shifts that accelerate the reaction but are not readily observed by standard experimental approaches, and (2) influencers favorably shift the pK_as of proton-transferring nucleobases and manifest themselves as "wavy" rate-pH profiles. We identify parallels with the protein enzyme literature, including reverse protonation and wavelike behavior, while pointing out that RNA is more prone to reverse protonation. The complexities uncovered, which arise from simple pairwise interactions, should aid deconvolution of complex rate-pH profiles for RNA and protein enzymes and suggest veiled catalytic devices for promoting catalysis that can be tested by experiment and calculation.

Using competition assays to quantitatively model cooperative binding by transcription factors and other ligands

BACKGROUND

The affinities of DNA binding proteins for target sites can be used to model the regulation of gene expression. These proteins can bind to DNA cooperatively, strongly impacting their affinity and specificity. However, current methods for measuring cooperativity do not provide the means to accurately predict binding behavior over a wide range of concentrations.

METHODS

We use standard computational and mathematical methods, and develop novel methods as described in Results.

RESULTS

We explore some complexities of cooperative binding, and develop an improved method for relating in vitro measurements to in vivo function, based on ternary complex formation. We derive expressions for the equilibria among the various complexes, and explore the limitations of binding experiments that model the system using a single parameter. We describe how to use single-ligand binding and ternary complex formation in tandem to determine parameters that have thermodynamic relevance. We develop an improved method for finding both single-ligand dissociation constants and concentrations simultaneously. We show how the cooperativity factor can be found when only one of the single-ligand dissociation constants can be measured.

CONCLUSIONS

The methods that we develop constitute an optimized approach to accurately model cooperative binding.

GENERAL SIGNIFICANCE

The expressions and methods we develop for modeling and analyzing DNA binding and cooperativity are applicable to most cases where multiple ligands bind to distinct sites on a common substrate. The parameters determined using these methods can be fed into models of higher-order cooperativity to increase their predictive power.

Deciphering Selectivity in Organic Reactions: A Multifaceted Problem

Computational chemistry has made a sustained contribution to the understanding of chemical reactions. In earlier times, half a century ago, the goal was to distinguish allowed from forbidden reactions (e.g., Woodward-Hoffmann rules), that is, reactions with low or high to very high activation barriers. A great achievement of computational chemistry was also to contribute to the determination of structures with the bonus of proposing a rationalization (e.g., anomeric effect, isolobal analogy, Gillespie valence shell pair electron repulsion rules and counter examples, Wade-Mingos rules for molecular clusters). With the development of new methods and the constant increase in computing power, computational chemists move to more challenging problems, close to the daily concerns of the experimental chemists, in determining the factors that make a reaction both efficient and selective: a key issue in organic synthesis. For this purpose, experimental chemists use advanced synthetic and analytical techniques to which computational chemists added other ways of determining reaction pathways. The transition states and intermediates contributing to the transformation of reactants into the desired and undesired products can now be determined, including their geometries, energies, charges, spin densities, spectroscopy properties, etc. Such studies remain challenging due to the large number of chemical species commonly present in the reactive media whose role may have to be determined. Calculating chemical systems as they are in the experiment is not always possible, bringing its own share of complexity through the large number of atoms and the associated large number of conformers to consider. Modeling the chemical species with smaller systems is an alternative that historically led to artifacts. Another important topic is the choice of the computational method. While DFT is widely used, the vast diversity of functionals available is both an opportunity and a challenge. Though chemical knowledge helps, the relevant computational method is best chosen in conjunction with the nature of the experimental systems and many studies have been concerned with this topic. We will not address this aspect but give references in the text. Usually, a computational study starts with the validation of the method by means of benchmark calculations vs accurate experimental data or state-of-the-art calculations. Finally, computational chemists can bring more than the sole determination of the reaction pathways through the analysis of the electronic structure. In our case, we have privileged the NBO analysis, which has the advantage of describing interactions on the basis of terms and concepts that are shared within the chemical community. In this Account, we have chosen to select representative reactions from our own work to highlight the diversity of situations than can be addressed nowadays. These include selective activation of C(sp(3))-H bonds, selective reactions with low energy barriers, involving closed shell or radical species, the role of noncovalent interactions, and the importance of considering side reactions.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

A self-consistent approach for determining pairwise interactions that underlie channel activation

Net free-energy measurements can be combined with mutant cycle analysis to determine interaction energies between specific amino acid pairs during channel activation.

Signaling proteins such as ion channels largely exist in two functional forms, corresponding to the active and resting states, connected by multiple intermediates. Multiparametric kinetic models based on sophisticated electrophysiological experiments have been devised to identify molecular interactions of these conformational transitions. However, this approach is arduous and is not suitable for large-scale perturbation analysis of interaction pathways. Recently, we described a model-free method to obtain the net free energy of activation in voltage- and ligand-activated ion channels. Here we extend this approach to estimate pairwise interaction energies of side chains that contribute to gating transitions. Our approach, which we call generalized interaction-energy analysis (GIA), combines median voltage estimates obtained from charge-voltage curves with mutant cycle analysis to ascertain the strengths of pairwise interactions. We show that, for a system with an arbitrary gating scheme, the nonadditive contributions of amino acid pairs to the net free energy of activation can be computed in a self-consistent manner. Numerical analyses of sequential and allosteric models of channel activation also show that this approach can measure energetic nonadditivities even when perturbations affect multiple transitions. To demonstrate the experimental application of this method, we reevaluated the interaction energies of six previously described long-range interactors in the Shaker potassium channel. Our approach offers the ability to generate detailed interaction energy maps in voltage- and ligand-activated ion channels and can be extended to any force-driven system as long as associated “displacement” can be measured.

A B3LYP and MP2(full) theoretical investigation into cooperativity effects, aromaticity and thermodynamic properties in the Na(+)⋯benzonitrile⋯H2O ternary complex

The cooperativity effects between H-bonding and Na(+)⋯π or Na(+)⋯σ interactions in Na(+)⋯benzonitrile⋯H2O complexes were investigated using the B3LYP and MP2(full) methods with 6-311++G(2d,p) and aug-cc-pVTZ basis sets. The thermodynamic cooperativity and the influence of this cooperativity on aromaticity was evaluated by nucleus-independent chemical shifts (NICS). The results showed that the influence of the Na(+)⋯σ or Na(+)⋯π interaction on the hydrogen bond is more pronounced than that of the latter on the former. The cooperativity effect appeared in the Na(+)⋯σ interaction complex while the anti-cooperativity effect tended to be in the Na(+)⋯π system. The change in enthalpy is the major factor driving cooperativity. Thermodynamic cooperativity is not in accordance with the cooperativity effect evaluated by the change of interaction energy. The ring aromaticity of is weakened while the bond dissociation energy (BDE) of the C-CN bond increases upon ternary complex formation. The cooperativity effect (E coop) correlates with R c (NICS(1)ternary/NICS(1)binary) and ΔΔδ (Δδ ternary - Δδ binary) involving the ring and C ≡ N bond, as well as R BDE(C-CN) [BDE(C-CN)ternary/BDE(C-CN)binary], respectively. AIM (atoms in molecules) analysis confirms the existence of cooperativity.

Identification of serpin determinants of specificity and selectivity for furin inhibition through studies of α1PDX (α1-protease inhibitor Portland)-serpin B8 and furin active-site loop chimeras

α1-Protease inhibitor Portland (α1PDX) is an engineered serpin family inhibitor of the proprotein convertase (PC), furin, that exhibits high specificity but limited selectivity for inhibiting furin over other PC family proteases. Here, we characterize serpin B8, a natural inhibitor of furin, together with α1PDX-serpin B8 and furin-PC chimeras to identify determinants of serpin specificity and selectivity for furin inhibition. Replacing reactive center loop (RCL) sequences of α1PDX with those of serpin B8 demonstrated that both the P4-P1 RXXR recognition sequence as well as the P1'-P5' sequence are critical determinants of serpin specificity for furin. Alignments of PC catalytic domains revealed four variable active-site loops whose role in furin reactivity with serpin B8 was tested by engineering furin-PC loop chimeras. The furin(298-300) loop but not the other loops differentially affected furin reactivity with serpin B8 and α1PDX in a manner that depended on the serpin RCL-primed sequence. Modeling of the serpin B8-furin Michaelis complex identified serpin exosites in strand 3C close to the 298-300 loop whose substitution in α1PDX differentially affected furin reactivity depending on the furin loop and serpin RCL-primed sequences. These studies demonstrate that RCL-primed residues, strand 3C exosites, and the furin(298-300) loop are critical determinants of serpin reactivity with furin, which may be exploited in the design of specific and selective α1PDX inhibitors of PCs.

Water mediated ligand functional group cooperativity: the contribution of a methyl group to binding affinity is enhanced by a COO(-) group through changes in the structure and thermodynamics of the hydration waters of ligand-thermolysin complexes

Ligand functional groups can modulate the contributions of one another to the ligand-protein binding thermodynamics, producing either positive or negative cooperativity. Data presented for four thermolysin phosphonamidate inhibitors demonstrate that the differential binding free energy and enthalpy caused by replacement of a H with a Me group, which binds in the well-hydrated S2' pocket, are more favorable in presence of a ligand carboxylate. The differential entropy is however less favorable. Dissection of these differential thermodynamic parameters, X-ray crystallography, and density-functional theory calculations suggest that these cooperativities are caused by variations in the thermodynamics of the complex hydration shell changes accompanying the H→Me replacement. Specifically, the COO(-) reduces both the enthalpic penalty and the entropic advantage of displacing water molecules from the S2' pocket and causes a subsequent acquisition of a more enthalpically, less entropically, favorable water network. This study contributes to understanding the important role water plays in ligand-protein binding.

A novel type of allosteric regulation: functional cooperativity in monomeric proteins

Cooperative functional properties and allosteric regulation in cytochromes P450 play an important role in xenobiotic metabolism and define one of the main mechanisms of drug-drug interactions. Recent experimental results suggest that ability to bind simultaneously two or more small organic molecules can be the essential feature of cytochrome P450 fold, and often results in rich and complex pattern of allosteric behavior. Manifestations of non-Michaelis kinetics include homotropic and heterotropic activation and inhibition effects depending on the stoichiometric ratios of substrate and effector, changes in the regio- and stereospecificity of catalytic transformations, and often give rise to the clinically important drug-drug interactions. In addition, functional response of P450 systems is modulated by the presence of specific and non-specific effector molecules, metal ions, membrane incorporation, formation of homo- and hetero-oligomers, and interactions with the protein redox partners. In this article we briefly overview the main factors contributing to the allosteric effects in cytochromes P450 with the main focus on the sources of cooperative behavior in xenobiotic metabolizing monomeric heme enzymes with their conformational flexibility and extremely broad substrate specificity. The novel mechanism of functional cooperativity in P450 enzymes does not require substantial binding cooperativity, rather it implies the presence of one or more binding sites with higher affinity than the single catalytically active site in the vicinity of the heme iron.

Membrane anchoring and interaction between transmembrane domains are crucial for K+ channel function

The small viral channel Kcv is a Kir-like K(+) channel of only 94 amino acids. With this simple structure, the tetramer of Kcv represents the pore module of all complex K(+) channels. To examine the structural contribution of the transmembrane domains (TMDs) to channel function, we performed Ala scanning mutagenesis of the two domains and tested the functionality of the mutants in a yeast complementation assay. The data reveal, in combination with computational models, that the upper halves of both TMDs, which face toward the external medium, are rather rigid, whereas the inner parts are more flexible. The rigidity of the outer TMD is conferred by a number of essential aromatic amino acids that face the membrane and probably anchor this domain in the bilayer. The inner TMD is intimately connected with the rigid part of the outer TMD via π···π interactions between a pair of aromatic amino acids. This structural principle is conserved within the viral K(+) channels and also present in Kir2.2, implying a general importance of this architecture for K(+) channel function.

Deconstructing thermodynamic parameters of a coupled system from site-specific observables

Cooperative interactions mediate information transfer between structural domains of a protein molecule and are major determinants of protein function and modulation. The prevalent theories to understand the thermodynamic origins of cooperativity have been developed to reproduce the complex behavior of a global thermodynamic observable such as ligand binding or enzyme activity. However, in most cases the measurement of a single global observable cannot uniquely define all the terms that fully describe the energetics of the system. Here we establish a theoretical groundwork for analyzing protein thermodynamics using site-specific information. Our treatment involves extracting a site-specific parameter (defined as χ value) associated with a structural unit. We demonstrate that, under limiting conditions, the χ value is related to the direct interaction terms associated with the structural unit under observation and its intrinsic activation energy. We also introduce a site-specific interaction energy term (χ(diff)) that is a function of the direct interaction energy of that site with every other site in the system. When combined with site-directed mutagenesis and other molecular level perturbations, analyses of χ values of site-specific observables may provide valuable insights into protein thermodynamics and structure.

Role of unsatisfied hydrogen bond acceptors in RNA energetics and specificity

RNA plays essential roles in much of biology. These functions are dictated by structures mediated by hydrogen bonding, stacking, electrostatics, and steric interactions. Roles of unsatisfied hydrogen bond functionalities in these structures are less well understood. Herein, we evaluated the energetic contributions of unsatisfied hydrogen bonding groups by placing chemically modified substituents in select internal positions in RNA helices and conducting thermodynamic studies. We find that unsatisfied carbonyl groups make exceptional contributions to structure formation (approximately 3 kcal/mol in free energy), most likely due to a combination of strain and dehydration effects. Thus, unsatisfied hydrogen bonding groups are likely key determinants in the folding energetics and specificity of many RNA and DNA molecules and may be especially important in tertiary structure interactions.

Cooperative properties of cytochromes P450

Cytochromes P450 form a large and important class of heme monooxygenases with a broad spectrum of substrates and corresponding functions, from steroid hormone biosynthesis to the metabolism of xenobiotics. Despite decades of study, the molecular mechanisms responsible for the complex non-Michaelis behavior observed with many members of this superfamily during metabolism, often termed 'cooperativity', remain to be fully elucidated. Although there is evidence that oligomerization may play an important role in defining the observed cooperativity, some monomeric cytochromes P450, particularly those involved in xenobiotic metabolism, also display this behavior due to their ability to simultaneously bind several substrate molecules. As a result, formation of distinct enzyme-substrate complexes with different stoichiometry and functional properties can give rise to homotropic and heterotropic cooperative behavior. This review aims to summarize the current understanding of cooperativity in cytochromes P450, with a focus on the nature of cooperative effects in monomeric enzymes.

Thinking inside the box: designing, implementing, and interpreting thermodynamic cycles to dissect cooperativity in RNA and DNA folding

Double and triple mutant thermodynamic cycles provide a means to dissect the cooperativity of RNA and DNA folding at both the secondary and tertiary structural levels through use of the thermodynamic box or cube. In this article, we describe three steps for applying thermodynamic cycles to nucleic acid folding, with considerations of both conceptual and experimental features. The first step is design of an appropriate system and development of hypotheses regarding which residues might interact. Next is implementing this design in terms of a tractable experimental strategy, with an emphasis on UV melting. The final step, and the one we emphasize the most, is interpreting mutant cycles in terms of coupling between specific residues in the RNA or DNA. Coupling free energy in the absence and presence of changes elsewhere in the molecule is discussed in terms of specific folding models, including stepwise folding and concerted changes. Last, we provide a practical section on the use of commercially available software (KaleidaGraph) to fit melting data, along with a consideration of error propagation. Along the way, specific examples are chosen from the literature to illustrate the methods. This article is intended to be accessible to the biochemist or biologist without extensive thermodynamics background.

Structures, kinetics, thermodynamics, and biological functions of RNA hairpins

Most RNA comprises one strand and therefore can fold back on itself to form complex structures. At the heart of these structures is the hairpin, which is composed of a stem having Watson-Crick base pairing and a loop wherein the backbone changes directionality. First, we review the structure of hairpins including diversity in the stem, loop, and closing base pair. The function of RNA hairpins in biology is discussed next, including roles for isolated hairpins, as well as hairpins in the context of complex tertiary structures. We describe the kinetics and thermodynamics of hairpin folding including models for hairpin folding, folding transition states, and the cooperativity of folding. Lastly, we discuss some ways in which hairpins can influence the folding and function of tertiary structures, both directly and indirectly. RNA hairpins provide a simple means of controlling gene expression that can be understood in the language of physical chemistry.

Coarse-grained lattice model for investigating the role of cooperativity in molecular recognition

Equilibrium aspects of the molecular recognition of rigid biomolecules are investigated using coarse-grained lattice models. The analysis is carried out in two stages. First, an ensemble of probe molecules is designed with respect to the target biomolecule. The recognition ability of the probe ensemble is then investigated by calculating the free energy of association. The influence of cooperative and anticooperative effects accompanying the association of the target and probe molecules is studied. Numerical findings are presented and compared to analytical results which can be obtained in the limit of dominating cooperativity and in the mean-field formulation of the models.

A thermodynamic ligand binding study of the third PDZ domain (PDZ3) from the mammalian neuronal protein PSD-95

The thermodynamic parameters associated with the binding of several series of linear peptides to the third PDZ domain (PDZ3) of the postsynaptic density 95 protein (PSD-95) have been measured using isothermal titration calorimetry (ITC). Two strategies were pursued in developing these binding ligands: (1) systematic N-terminal truncation of sequences derived from the C-terminal regions of identified PDZ3-binding proteins (CRIPT, neuroligin-1, and citron) and (2) selective mutation of specific positions within a consensus hexapeptide (KKETEV) known to bind PDZ3. Each synthetically prepared peptide was used to titrate PDZ3, which yielded the changes in Gibbs free energy (DeltaG), enthalpy (DeltaH), and entropy (TDeltaS) for the binding event. Selected peptides were subjected to additional analysis, which entailed (1) measuring the change in heat capacity (DeltaCp) upon association, to assess the character of the binding interface, and (2) constructing thermodynamic double mutant cycles, to determine the presence of cooperative effects. From the first series, the CRIPT protein proved to be the better source for higher affinity sequences. From the second series, enhanced binding was associated with peptides that closely adhered to the established motif for class I PDZ domain C-termini, X-(T/S)-X-(V/I/L), and more specifically to a narrower motif of X-T-X-V. Further, in both series a length of six residues was necessary and sufficient to capture maximal affinity. In addition, there were significant influences upon binding by modifying the abutting "X" positions. The cumulative results provide greater detail into the specific nature of ligand binding to PDZ3 and will assist in the development of selective molecular probes for the study of this and structurally homologous PDZ domains.

Folding cooperativity in RNA and DNA is dependent on position in the helix

Secondary structural motifs play essential roles in the folding and function of RNA and DNA molecules. Previous work from our lab compared the folding of small DNA and RNA hairpin loops containing a sheared GA pair [Moody, E. M., Feerar, J. C., and Bevilacqua, P. C. (2004) Biochemistry 43, 7992-7998]. We found that the small DNA hairpins fold in a highly cooperative manner with indirect coupling, while their RNA counterparts fold in a much less cooperative fashion and display direct coupling. Herein, we extend this study to the double-stranded helix. We carried out double mutant cycles on base pairs having identical nearest-neighbor contexts but located in either external or internal helical registers. In the external register, both RNA and DNA exhibit extensive folding cooperativity between the penultimate and terminal base pair, which is independent of mismatch identity. In contrast, DNA exhibits virtually no folding cooperativity in the center of the helix, while RNA maintains substantial coupling, which is dependent on mismatch identity. Two models account for these non-nearest-neighbor effects: one involves the unfavorable entropy of helix initiation common to DNA and RNA, and the other involves steric and electrostatic strain peculiar to RNA. These data show that RNA can display cooperativity less than, greater than, or equal to that of DNA depending on context and position.

Chemical double-mutant cycles: dissecting non-covalent interactions

Thermodynamic double-mutant cycles and triple-mutant boxes are widely employed for the experimental quantification of non-covalent interactions and cooperative effects in proteins. This review describes the application of these powerful methodologies to the study of non-covalent interactions in synthetic systems.

A study of collective atomic fluctuations and cooperativity in the U1A-RNA complex based on molecular dynamics simulations

Cooperative interactions play an important role in recognition and binding in macromolecular systems. In this study, we find that cross-correlated atomic fluctuations can be used to identify cooperative networks in a protein-RNA system. The dynamics of the RRM-containing protein U1A-stem loop 2 RNA complex have been calculated theoretically from a 10 ns molecular dynamics (MD) simulation. The simulation was analyzed by calculating the covariance matrix of all atomic fluctuations. These matrix elements are then presented in the form of a two-dimensional grid, which displays fluctuations on a per residue basis. The results indicate the presence of strong, selective cross-correlated fluctuations throughout the RRM in U1A-RNA. The atomic fluctuations correspond well with previous biophysical studies in which a multiplicity of cooperative networks have been reported and indicate that the various networks identified in separate individual experiments are fluctuationally correlated into a hyper-network encompassing most of the RRM. The calculated results also correspond well with independent results from a statistical covariance analysis of 330 aligned RRM sequences. This method has significant implications as a predictive tool regarding cooperativity in the protein-nucleic acid recognition process.

Coarse-grained lattice model for molecular recognition

We present a simple model which allows us to investigate the equilibrium aspects of molecular recognition between rigid biomolecules on a generic level. Using a two-stage approach, which consists of a design and a testing step, the role of cooperativity and of varying bond strength in molecular recognition is investigated. Cooperativity is found to enhance selectivity. In complexes which require a high binding flexibility, a small number of strong bonds seems to be favored compared to a situation with many but weak bonds.

Probing the Ligand-Binding Specificity and Analyzing the Folding State of SPOT-Synthesized FBP28 WW Domain Variants

The WW domains are known as the smallest naturally occurring, monomeric, triple-stranded, antiparallel beta-sheet domains. Hence, we chose the FBP28 WW domain as a model to investigate the stability of the beta-sheet structure at the amino acid level in the context of its function (ligand binding). The structure-function relationship was investigated through a complete substitution analysis of the FBP28 WW domain, with variants synthesized as a cellulose-bound peptide array. The functionality of the FBP28 WW domain variants was examined by probing the peptide array for ligand binding. In addition, selected FBP28 WW domain variants were investigated by CD measurements to determine the stability of the antiparallel beta-sheet structure. We discuss the correlation between structure stability and functionality for the FBP28 WW domain, as well as the effect of ligand-induced structure stabilization.

Alanine scanning of MS2 coat protein reveals protein-phosphate contacts involved in thermodynamic hot spots

The co-crystal structure of the MS2 coat protein dimer with its RNA operator reveals eight amino acid side-chains contacting seven of the RNA phosphates. These eight amino acids and five nearby control positions were individually changed to an alanine residue and the binding affinities of the mutant proteins to the RNA were determined. In general, the data agreed well with the crystal structure and previous RNA modification data. Interestingly, amino acid residues that are energetically most important for complex formation cluster in the middle of the RNA binding interface, forming thermodynamic hot spots, and are surrounded by energetically less relevant amino acids. In order to evaluate whether or not a given alanine mutation causes a global change in the RNA-protein interface, the affinities of the mutant proteins to RNAs containing one of 14 backbone modifications spanning the entire interface were determined. In three of six protein mutations tested, thermodynamic coupling between the site of the mutation and RNA groups that can be even more than 16 A away was detected. This suggests that, in some cases, the mutation may subtly alter the entire protein-RNA interface.

Historical review: Negative efficacy and the constitutive activity of G-protein-coupled receptors

The idea that a receptor can produce signalling without agonist intervention and that several antagonists can be 'active' in repressing such spontaneous activity is contained in the concept of ligand-induced conformational changes. Yet, this idea was neglected by pharmacologists for many years. In this article, we review the events that brought inverse agonism and constitutive activity to general attention and made this phenomenon a topic of current research. We also suggest a classification of antagonists based on the cooperativity that links their primary site of interaction with other functional domains of the receptor.

Role of Anions on the Crystal Structures of Copper(II) and Zinc(II) Complexes of a Tunable Butterfly Cyclophane Macrocycle

Three crystal structures of a ditopic cyclophane ligand (L) in which two 1,5,8,12-tetraamine molecules have been attached through methylene spacers to the ortho positions of a benzene ring are reported. The first one (1) corresponds to the tetraprotonated free macrocycle (H4L4+) having two tetrachlorozincate(II) counteranions (C24H54O2N8Cl8Zn2, a = 9.1890(2) A, b = 14.0120(3) A, c = 15.3180(3) A, alpha = 89.2320(7) degrees , beta = 82.0740(6) degrees , gamma = 83.017(1) degrees , Z = 2.00, triclinic, P); the second one (2) is of a binuclear Cu2+ complex having coordinated chloride anions and perchlorate counteranions (C24H58O14N8Cl4Cu2 a = 9.9380(2) A, b = 30.2470(6) A, c = 53.143(1) A, orthorhombic, F2dd, Z = 18), and the third one (3) corresponds to an analogous Zn2+ complex that has been crystallized using triflate as counteranion (C26H(51.2)O(6.6)N8Cl2F6S2Zn2 a = 8.472(5) A, b = 9.310(5), c = 13.745(5) A, alpha = 84.262(5) degrees , beta = 77.490(5) degrees , gamma = 73.557(5) degrees , triclinic, P, Z = 2). The analysis of the crystallographic data clearly shows that the conformation of the macrocycle and, in consequence, the overall architecture of the crystals are controlled by the anions present in the moiety, pi-pi-stacking associations, and hydrogen bonding interactions. The protonation and stability constants for the formation of the Cu2+ and Zn2+ complexes in aqueous solution have been determined potentiometrically in 0.15 mol dm(-3) NaClO4 at 298.1 K. Intramolecular hydrogen bonding defines the protonation behavior of the compound. Positive cooperativity is observed in the formation of the Cu2+ complexes.

Hydrogen Bonding and Cooperativity in Isolated and Hydrated Sugars: Mannose, Galactose, Glucose, and Lactose

The conformation of phenyl-substituted monosaccharides (mannose, galactose, and glucose) and their singly hydrated complexes has been investigated in the gas phase by means of a combination of mass selected, conformer specific ultraviolet and infrared double resonance hole burning spectroscopy experiments, and ab initio quantum chemistry calculations. In each case, the water molecule inserts into the carbohydrate at a position where it can replace a weak intramolecular interaction by two stronger intermolecular hydrogen bonds. The insertion can produce significant changes in the conformational preferences of the carbohydrates, and there is a clear preference for structures where cooperative effects enhance the stability of the monosaccharide conformers to which the water molecule chooses to bind. The conclusions drawn from the study of monosaccharide-water complexes are extended to the disaccharide lactose and discussed in the light of the underlying mechanisms that may be involved in the binding of carbohydrate assemblies to proteins and the involvement, or not, of key structural water molecules.

Structural and energetic consequences of expanding a highly cooperative stable DNA hairpin loop

Many hairpin loops are expanded versions of smaller, stable ones. Herein we investigate the extent to which the energetics and structure of d(cGNAg) hairpin loops will tolerate sequence variation. Changing the closing base pair from CG to GC was found to completely eliminate loop-loop interactions; in contrast, expanding the loop at the 3'-end resulted in similar energetics and nonadditivity parameters as the parent loop, suggesting that loop-loop interactions remain intact and highly coupled upon expansion. Together, these data suggest that the CG closing base pair forms an essential platform upon which a stable d(GNA) hairpin loop can fold and that this loop can undergo 3'-expansion with little effect to its structure or energetics.

Coupling interactions between voltage sensors of the sodium channel as revealed by site-specific measurements

The voltage-sensing S4 segments in the sodium channel undergo conformational rearrangements in response to changes in the electric field. However, it remains unclear whether these structures move independently or in a coordinated manner. Previously, site-directed fluorescence measurements were shown to track S4 transitions in each of the four domains (Chanda, B., and F. Bezanilla. 2002. J. Gen. Physiol. 120:629–645). Here, using a similar technique, we provide direct evidence of coupling interactions between voltage sensors in the sodium channel. Pairwise interactions between S4s were evaluated by comparing site-specific conformational changes in the presence and absence of a gating perturbation in a distal domain. Reciprocity of effect, a fundamental property of thermodynamically coupled systems, was measured by generating converse mutants. The magnitude of a local gating perturbation induced by a remote S4 mutation depends on the coupling strength and the relative equilibrium positions of the two voltage sensors. In general, our data indicates that the movement of all four voltage sensors in the sodium channel are coupled to a varying extent. Moreover, a gating perturbation in S4-DI has the largest effect on the activation of S4-DIV and vice versa, demonstrating an energetic linkage between S4-DI and S4-DIV. This result suggests a physical mechanism by which the activation and inactivation process may be coupled in voltage-gated sodium channels. In addition, we propose that cooperative interactions between voltage sensors may be the mechanistic basis for the fast activation kinetics of the sodium channel.