Thermodynamic consequences of disrupting a water-mediated hydrogen bond network in a protein:pheromone complex

Exploring Glycan Binding Specificity of Odorranalectin by Alanine Scanning Library

Fluorescently labelled alanine scan analogues of odorranalectin (OL), a cyclic peptide that exhibits lectin like properties, were screened for binding BSA-conjugated monosaccharides using an enzyme-linked lectin assay (ELLA). Results revealed that Lys⁵, Phe⁷, Tyr⁹, Gly¹², Leu¹⁴, and Thr¹⁷ were crucial for binding BSA-L-fucose, BSA-D-galactose and BSA-N-acetyl-D-galactosamine. Notably, Ala substitution of Ser³, Pro⁴, and Val¹³ resulted in higher binding affinities compared to the native OL. The obtained data also indicated that Arg⁸ plays an important role in differentiation of binding for BSA-L-fucose/D-galactose from BSA-N-acetyl-D-galactosamine. The thermodynamics of binding of the selected alanine analogues was evaluated by isothermal titration calorimetry. Low to moderate binding affinities were determined for the tetravalent MUC1 glycopeptide and asialofetuin, respectively, and high for the fucose rich polysaccharide, fucoidan. The thermodynamic profile of interactions with asialofetuin exhibits shift to an entropy-driven mechanism compared to the fucoidan, which displayed an enthalpyentropy compensation, typically associated with the carbohydratelectin recognition process.

Computational and Structure-Based Development of High Potent Cell-Active Covalent Inhibitor Targeting the Peptidyl-Prolyl Isomerase NIMA-Interacting-1 (Pin1)

The unique proline isomerase peptidyl-prolyl isomerase NIMA-interacting-1 (Pin1) is reported to activate numerous cancer-driving pathways simultaneously, and aberrant Pin1 activation is present in many human cancers. Here, we identified a novel hit compound, ZL-Pin01, that covalently modified Pin1 at Cys113 with an half-maximal inhibitory concentration (IC₅₀) of 1.33 ± 0.07 μM through screening an in-house library. Crystallographic study drove the process of structure-guided optimization and led to the potent inhibitor ZL-Pin13 with an IC₅₀ of 0.067 ± 0.03 μM. We obtained four co-crystal structures of Pin1 complexed with inhibitors that elucidated the detailed binding mode of the derivatives with Pin1. Interestingly, the co-crystal of Pin1 with ZL-Pin13 obtained by co-crystallization revealed the conformational change of Gln129 induced by the inhibitor. Furthermore, ZL-Pin13 effectively inhibited the proliferation and downregulated the Pin1 substrates in MDA-MB-231 cells. Collectively, we developed a potent covalent inhibitor of Pin1, ZL-Pin13, which could be an effective probe for studying the functional roles of Pin1.

Cationic polythiophene-anionic fullerene pair in water and water-dioxane: studies on hydrogen bonding capabilities, kinetic and thermodynamic properties

Despite the vast array of solution- and solid-state bio-analytical, bioelectronic and optoelectronic applications of cationic polythiophenes (CPTs), the number of studies focused on the role of hydrogen bonding (H-bonding) between these and other molecules is scarce, regardless of whether H-bonding is expected to play an important role in several such applications. Also, despite the advantages of using cosolvents to systematically examine the molecular interactions, there are no such studies for CPTs to our knowledge. This work presents a steady-state UV-vis/fluorescence spectroscopic, kinetic and thermodynamic study on the H-bonding interactions between a water-soluble, cationic-anionic (isothiouronium-tetraphosphonate), polythiophene-fullerene donor-acceptor pair with two-point, charge-assisted H-bonding (CAHB) capabilities, tuned using water or a 1,4-dioxane-water mixture (W-DI). Both solvents generate photoinduced electron transfer (PET), fluorescence resonance energy transfer (FRET), spontaneous binding, H-bonding, ground-state complexing via multiple site binding, formation of micelle-like aggregates and equivalence points at a similar concentration of the quencher. However, in comparison with water, W-DI promotes less-ordered, less packed micellar aggregates, due to hydrophobic desolvation of the H-bond and larger solvent displacement during the PT1-4Fo complexation. This would decrease the extent of charge-transfer and the size of the sphere-of-quenching, mainly by displacements or rotations of the H-bonds, instead of elongations, together with a possible larger extent of diffusion-controlled static quenching. At [4Fo] larger than the equivalence point the micelles formed in water do not have available binding sites due to a tighter aggregation, causing a decrease in the quenching efficiency, while the micelles formed in W-DI start showing larger quenching efficiencies, possibly due to an increase in entropy that overcomes the desolvation of the H-bonding. These results could be useful when analyzing outputs from systems including CPTs with H-bonding capabilities, operating in (or casted from) solvents with clear differences in polarity and/or H-bonding capacity.

Nonadditivity in public and inhouse data: implications for drug design

Numerous ligand-based drug discovery projects are based on structure-activity relationship (SAR) analysis, such as Free-Wilson (FW) or matched molecular pair (MMP) analysis. Intrinsically they assume linearity and additivity of substituent contributions. These techniques are challenged by nonadditivity (NA) in protein-ligand binding where the change of two functional groups in one molecule results in much higher or lower activity than expected from the respective single changes. Identifying nonlinear cases and possible underlying explanations is crucial for a drug design project since it might influence which lead to follow. By systematically analyzing all AstraZeneca (AZ) inhouse compound data and publicly available ChEMBL25 bioactivity data, we show significant NA events in almost every second assay among the inhouse and once in every third assay in public data sets. Furthermore, 9.4% of all compounds of the AZ database and 5.1% from public sources display significant additivity shifts indicating important SAR features or fundamental measurement errors. Using NA data in combination with machine learning showed that nonadditive data is challenging to predict and even the addition of nonadditive data into training did not result in an increase in predictivity. Overall, NA analysis should be applied on a regular basis in many areas of computational chemistry and can further improve rational drug design.

The allergen Mus m 1.0102: Cysteine residues and molecular allergology

Mus m 1.0102 is a member of the mouse Major Urinary Protein family, belonging to the Lipocalins superfamily. Major Urinary Proteins (MUPs) are characterized by highly conserved structural motifs. These include a disulphide bond, involved in protein oxidative folding and protein structure stabilization, and a free cysteine residue, substituted by serine only in the pheromonal protein Darcin (MUP20). The free cysteine is recognized as responsible for the onset of inter- or intramolecular thiol/disulphide exchange, an event that favours protein aggregation. Here we show that the substitution of selected cysteine residues modulates Mus m 1.0102 protein folding, fold stability and unfolding reversibility, while maintaining its allergenic potency. Recombinant allergens used for immunotherapy or employed in allergy diagnostic kits require, as essential features, conformational stability, sample homogeneity and proper immunogenicity. In this perspective, recombinant Mus m 1.0102 might appear reasonably adequate as lead molecule because of its allergenic potential and thermal stability. However, its modest resistance to aggregation renders the protein unsuitable for pharmacological preparations. Point mutation is considered a winning strategy. We report that, among the tested mutants, C138A mutant acquires a structure more resistant to thermal stress and less prone to aggregation, two events that act positively on the protein shelf life. Those features make that MUP variant an attractive lead molecule for the development of a diagnostic kit and/or a vaccine.

Effects of point mutations in the binding pocket of the mouse major urinary protein MUP20 on ligand affinity and specificity

The mouse Major Urinary Proteins (MUPs) contain a conserved β-barrel structure with a characteristic central hydrophobic pocket that binds a variety of volatile compounds. After release of urine, these molecules are slowly emitted in the environment where they play an important role in chemical communication. MUPs are highly polymorphic and conformationally stable. They may be of interest in the construction of biosensor arrays capable of detection of a broad range of analytes. In this work, 14 critical amino acids in the binding pocket involved in ligand interactions were identified in MUP20 using in silico techniques and 7 MUP20 mutants were synthesised and characterised to produce a set of proteins with diverse ligand binding profiles to structurally different ligands. A single amino acid substitution in the binding pocket can dramatically change the MUPs binding affinity and ligand specificity. These results have great potential for the design of new biosensor and gas-sensor recognition elements.

Triazole resistance mediated by mutations of a conserved active site tyrosine in fungal lanosterol 14α-demethylase

Emergence of fungal strains showing resistance to triazole drugs can make treatment of fungal disease problematic. Triazole resistance can arise due to single mutations in the drug target lanosterol 14α-demethylase (Erg11p/CYP51). We have determined how commonly occurring single site mutations in pathogenic fungi affect triazole binding using Saccharomyces cerevisiae Erg11p (ScErg11p) as a target surrogate. The mutations Y140F/H were introduced into full-length hexahistidine-tagged ScErg11p. Phenotypes and high-resolution X-ray crystal structures were determined for the mutant enzymes complexed with short-tailed (fluconazole and voriconazole) or long-tailed (itraconazole and posaconazole) triazoles and wild type enzyme complexed with voriconazole. The mutations disrupted a water-mediated hydrogen bond network involved in binding of short-tailed triazoles, which contain a tertiary hydroxyl not present in long-tailed triazoles. This appears to be the mechanism by which resistance to these short chain azoles occurs. Understanding how these mutations affect drug affinity will aid the design of azoles that overcome resistance.

Rational Design of Benzylidenehydrazinyl-Substituted Thiazole Derivatives as Potent Inhibitors of Human Dihydroorotate Dehydrogenase with in Vivo Anti-arthritic Activity

Human dihydroorotate dehydrogenase (hDHODH) is an attractive therapeutic target for the treatment of rheumatoid arthritis, transplant rejection and other autoimmune diseases. Based on the X-ray structure of hDHODH in complex with lead compound 7, a series of benzylidenehydrazinyl-substituted thiazole derivatives as potent inhibitors of hDHODH were designed and synthesized, of which 19 and 30 were the most potent with IC₅₀ values in the double-digit nanomolar range. Moreover, compound 19 displayed significant anti-arthritic effects and favorable pharmacokinetic profiles in vivo. Further X-ray structure and SAR analyses revealed that the potencies of the designed inhibitors were partly attributable to additional water-mediated hydrogen bond networks formed by an unexpected buried water between hDHODH and the 2-(2-methylenehydrazinyl)thiazole scaffold. This work not only elucidates promising scaffolds targeting hDHODH for the treatment of rheumatoid arthritis, but also demonstrates that the water-mediated hydrogen bond interaction is an important factor in molecular design and optimization.

Molecular dynamics study of the conformations of glycosidic linkages in sialic acid modified ganglioside GM3 analogues

The objective of the present study is to model the analogues of monosialoganglioside (GM3) by making modifications in its sialic acid residue with different substitutions in aqueous environment and to determine their structural stability based upon computational molecular dynamics. Molecular mechanics and molecular dynamics investigation was carried out to study the conformational preferences of the analogues of GM3. Dynamic simulations were carried out on the analogues of GM3 varying in the substituents at C-1, C-4, C-5, C-8 and C-9 positions of their sialic acid or Neuraminic acid (NeuAc) residue. The analogues are soaked in a periodic box of TIP3P water as solvent and subjected to a 10 ns molecular dynamics (MD) simulation using AMBER ff03 and gaff force fields with 30 ps equilibration. The analogue of GM3 with 9-N-succNeuAc (analogue5, C9 substitution) was observed to have the lowest energy of -6112.5 kcal/mol. Graphical analysis made on the MD trajectory reveals the direct and water mediated hydrogen bonds existing in these sialic acid analogues. The preferable conformations for glycosidic linkages of GM3 analogues found in different minimum energy regions in the conformational maps were identified. This study sheds light on the conformational preferences of GM3 analogues which may be essential for the design of GM3 analogues as inhibitors for different ganglioside specific pathogenic proteins such as bacterial toxins, influenza toxins and neuraminidases.

Quantitative analysis of pheromone-binding protein specificity

Many pheromones have very low water solubility, posing experimental difficulties for quantitative binding measurements. A new method is presented for determining thermodynamically valid dissociation constants for ligands binding to pheromone-binding proteins, using β-cyclodextrin as a solubilizer and transfer agent. The method is applied to LUSH, a Drosophila odorant-binding protein that binds the pheromone 11-cis vaccenyl acetate (cVA). Refolding of LUSH expressed in Escherichia coli was assessed by measuring N-phenyl-1-naphthylamine (NPN) binding and Förster resonance energy transfer between LUSH tryptophan 123 (W123) and NPN. Binding of cVA was measured from quenching of W123 fluorescence as a function of cVA concentration. The equilibrium constant for transfer of cVA between β-cyclodextrin and LUSH was determined from a linked equilibria model. This constant, multiplied by the β-cyclodextrin-cVA dissociation constant, gives the LUSH-cVA dissociation constant: ∼100 nM. It was also found that other ligands quench W123 fluorescence. The LUSH-ligand dissociation constants were determined to be ∼200 nM for the silk moth pheromone bombykol and ∼90 nM for methyl oleate. The results indicate that the ligand-binding cavity of LUSH can accommodate a variety ligands with strong binding interactions. Implications of this for the Laughlin, Ha, Jones and Smith model of pheromone reception are discussed.

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.

In search of a vaccine for mouse allergy: significant reduction of Mus m 1 allergenicity by structure-guided single-point mutations

BACKGROUND

Mouse urinary proteins are relevant allergens from mice urine. We used the recombinant protein Mus m 1 as an allergen model to identify if, by altering Mus m 1 architecture via single-point mutations, we could effectively modify its allergenicity.

METHODS

Based on structural considerations, we synthesized two single-point mutants, Mus m 1-Y120A and Mus m 1-Y120F, which were expected to harbor large structural alterations. Circular dichroism and fluorescence analysis showed significant conformational rearrangements of the aromatic side chains in the internal cavity of Mus m 1-Y120A when compared to Mus m 1-Y120F and Mus m 1. Evaluation of the allergenic potential of the recombinant molecules was performed in vitro with both immunochemical approaches and assays based on the measurement of basophil degranulation. Moreover, to assess the integrity of the T cell epitopes and as an in vitro measure of immunogenicity, we tested the reactivity of T lymphocytes from subjects allergic to mouse urine against proteins and synthetic peptides encompassing the immunodominant linear epitope containing the mutation.

RESULTS

We found that the selected point mutation was able to modulate the protein allergenicity, and to severely impair the recognition of Mus m 1 by IgE, while T cell reactivity was fully maintained.

CONCLUSIONS

In silico predicted, minimum selected structural modifications allowed to design one protein with reduced allergenicity and preserved immunogenicity. Structurally guided mutations can direct the design of proteins with reduced allergenicity which can be used as vaccines for a safer and more effective immunotherapy of allergic disorders.

Insights into scFv:drug binding using the molecular dynamics simulation and free energy calculation

Molecular dynamics simulations and free energy calculation have been performed to study how the single-chain variable fragment (scFv) binds methamphetamine (METH) and amphetamine (AMP). The structures of the scFv:METH and the scFv:AMP complexes are analyzed by examining the time-dependence of their RMSDs, by analyzing the distance between some key atoms of the selected residues, and by comparing the averaged structures with their corresponding crystallographic structures. It is observed that binding an AMP to the scFv does not cause significant changes to the binding pocket of the scFv:ligand complex. The binding free energy of scFv:AMP without introducing an extra water into the binding pocket is much stronger than scFv:METH. This is against the first of the two scenarios postulated in the experimental work of Celikel et al. (Protein Science 18, 2336 (2009)). However, adding a water to the AMP (at the position of the methyl group of METH), the binding free energy of the scFv:AMP-H2O complex, is found to be significantly weaker than scFv:METH. This is consistent with the second of the two scenarios given by Celikel et al. Decomposition of the binding energy into ligand-residue pair interactions shows that two residues (Tyr175 and Tyr177) have nearly-zero interactions with AMP in the scFv:AMP-H2O complex, whereas their interactions with METH in the scFv:METH complex are as large as -0.8 and -0.74 kcal mol(-1). The insights gained from this study may be helpful in designing more potent antibodies in treating METH abuse.

High resolution X-ray structures of mouse major urinary protein nasal isoform in complex with pheromones

In mice, the major urinary proteins (MUP) play a key role in pheromonal communication by binding and transporting semiochemicals. MUP-IV is the only isoform known to be expressed in the vomeronasal mucosa. In comparison with the MUP isoforms that are abundantly excreted in the urine, MUP-IV is highly specific for the male mouse pheromone 2-sec-butyl-4,5-dihydrothiazole (SBT). To examine the structural basis of this ligand preference, we determined the X-ray crystal structure of MUP-IV bound to three mouse pheromones: SBT, 2,5-dimethylpyrazine, and 2-heptanone. We also obtained the structure of MUP-IV with 2-ethylhexanol bound in the cavity. These four structures show that relative to the major excreted MUP isoforms, three amino acid substitutions within the binding calyx impact ligand coordination. The F103 for A along with F54 for L result in a smaller cavity, potentially creating a more closely packed environment for the ligand. The E118 for G substitution introduces a charged group into a hydrophobic environment. The sidechain of E118 is observed to hydrogen bond to polar groups on all four ligands with nearly the same geometry as seen for the water-mediated hydrogen bond network in the MUP-I and MUP-II crystal structures. These differences in cavity size and interactions between the protein and ligand are likely to contribute to the observed specificity of MUP-IV.

Long-timescale molecular-dynamics simulations of the major urinary protein provide atomistic interpretations of the unusual thermodynamics of ligand binding

The mouse major urinary protein (MUP) has proved to be an intriguing test bed for detailed studies on protein-ligand recognition. NMR, calorimetric, and modeling investigations have revealed that the thermodynamics of ligand binding involve a complex interplay between competing enthalpic and entropic terms. We performed six independent, 1.2 micros molecular-dynamics simulations on MUP--three replicates on the apo-protein, and three on the complex with the pheromone isobutylmethoxypyrazine. Our findings provide the most comprehensive picture to date of the structure and dynamics of MUP, and how they are modulated by ligand binding. The mechanical pathways by which amino acid side chains can transmit information regarding ligand binding to surface loops and either increase or decrease their flexibility (entropy-entropy compensation) are identified. Dewetting of the highly hydrophobic binding cavity is confirmed, and the results reveal an aspect of ligand binding that was not observed in earlier, shorter simulations: bound ligand retains extensive rotational freedom. Both of these features have significant implications for interpretations of the entropic component of binding. More generally, these simulations test the ability of current molecular simulation methods to produce a reliable and reproducible picture of protein dynamics on the microsecond timescale.

Role of interfacial water molecules in proline-rich ligand recognition by the Src homology 3 domain of Abl

The interaction of Abl-Src homology 3 domain (SH3) with the high affinity peptide p41 is the most notable example of the inconsistency existing between the currently accepted description of SH3 complexes and their binding thermodynamic signature. We had previously hypothesized that the presence of interfacial water molecules is partially responsible for this thermodynamic behavior. We present here a thermodynamic, structural, and molecular dynamics simulation study of the interaction of p41 with Abl-SH3 and a set of mutants designed to alter the water-mediated interaction network. Our results provide a detailed description of the dynamic properties of the interfacial water molecules and a molecular interpretation of the thermodynamic effects elicited by the mutations in terms of the modulation of the water-mediated hydrogen bond network. In the light of these results, a new dual binding mechanism is proposed that provides a better description of proline-rich ligand recognition by Abl-SH3 and that has important implications for rational design.

MM-PBSA Captures Key Role of Intercalating Water Molecules at a Protein-Protein Interface

The calculation of protein interaction energetics is of fundamental interest, yet accurate quantities are difficult to obtain due to the complex and dynamic nature of protein interfaces. This is further complicated by the presence of water molecules, which can exhibit transient interactions of variable duration and strength with the protein surface. The T-cell receptor (TCR) and its staphylococcal enterotoxin 3 (SEC3) binding partner are well-characterized examples of a protein−protein interaction system exhibiting interfacial plasticity, cooperativity, and additivity among mutants. Specifically engineered mutants induce intercalating interfacial water molecules, which subsequently enhance protein−protein binding affinity. In this work, we perform a set of molecular mechanics (MM) Poisson−Boltzmann (PB) surface area (SA) calculations on the wild type and two mutant TCR-SEC3 systems and show that the method is able to discriminate between weak and strong binders only when key explicit water molecules are included in the analysis. The results presented here point to the promise of MM-PBSA toward rationalizing molecular recognition at protein−protein interfaces, while establishing a general approach to handle explicit interfacial water molecules in such calculations.

Involvement of water in carbohydrate-protein binding: concanavalin A revisited

Ordered water molecules bound to protein surfaces, or in protein-ligand interfaces, are frequently observed by crystallography. The investigation of the impact of such conserved water molecules on protein stability and ligand affinity requires detailed structural, dynamic, and thermodynamic analyses. Several crystal structures of the legume lectin concanavalin A (Con A) bound to closely related carbohydrate ligands show the presence of a conserved water molecule that mediates ligand binding. Experimental thermodynamic and theoretical studies have examined the role of this conserved water in the complexation of Con A with a synthetic analog of the natural trisaccharide, in which a hydroxyethyl side chain replaces the hydroxyl group at the C-2 position in the central mannosyl residue. Molecular modeling earlier indicated (Clarke, C.; Woods, R. J.; Glushka, J.; Cooper, A.; Nutley, M. A.; Boons, G.-J. J. Am. Chem. Soc. 2001, 123, 12238-12247) that the hydroxyl group in this synthetic side chain could occupy a position equivalent to that of the conserved water, and thus might displace it. An interpretation of the experimental thermodynamic data, which was consistent with the displacement of the conserved water, was also presented. The current work reports the crystal structure of Con A with this synthetic ligand and shows that even though the position and interactions of the conserved water are distorted, this key water is not displaced by the hydroxyethyl moiety. This new structural data provides a firm basis for molecular dynamics simulations and thermodynamic integration calculations whose results indicate that differences in van der Waals contacts (insertion energy), rather than electrostatic interactions (charging energy) are fundamentally responsible for the lower affinity of the synthetic ligand. When combined with the new crystallographic data, this study provides a straightforward interpretation for the lower affinity of the synthetic analog; specifically, that it arises primarily from weaker interactions with the protein via the positionally perturbed conserved water. This interpretation is fully consistent with the experimental observations that the free energy of binding is enthalpy driven, that there is both less enthalpic gain and less entropic penalty for binding the synthetic ligand, relative to the natural trisaccharide, and that the entropic component does not arise from releasing an ordered water molecule from the protein surface to the bulk solvent.

Conserved water molecules stabilize the Omega-loop in class A beta-lactamases

A set of 49 high-resolution (<or=2.2 A) structures of the TEM, SHV, and CTX-M class A beta-lactamase families was systematically analyzed to investigate the role of conserved water molecules in the stabilization of the Omega-loop. Overall, 13 water molecules were found to be conserved in at least 45 structures, including two water positions which were found to be conserved in all structures. Of the 13 conserved water molecules, 6 are located at the Omega-loop, forming a dense cluster with hydrogen bonds to residues at the Omega-loop as well as to the rest of the protein. This layer of conserved water molecules is packed between the Omega-loop and the rest of the protein and acts as structural glue, which could reduce the flexibility of the Omega-loop. A correlation between conserved water molecules and conserved protein residues could in general not be detected, with the exception of the conserved water molecules at the Omega-loop. Furthermore, the evolutionary relationship between the three families, derived from the number of conserved water molecules, is similar to the relationship derived from phylogenetic analysis.

Modulation of androgen receptor activation function 2 by testosterone and dihydrotestosterone

The androgen receptor (AR) is transcriptionally activated by high affinity binding of testosterone (T) or its 5alpha-reduced metabolite, dihydrotestosterone (DHT), a more potent androgen required for male reproductive tract development. The molecular basis for the weaker activity of T was investigated by determining T-bound ligand binding domain crystal structures of wild-type AR and a prostate cancer somatic mutant complexed with the AR FXXLF or coactivator LXXLL peptide. Nearly identical interactions of T and DHT in the AR ligand binding pocket correlate with similar rates of dissociation from an AR fragment containing the ligand binding domain. However, T induces weaker AR FXXLF and coactivator LXXLL motif interactions at activation function 2 (AF2). Less effective FXXLF motif binding to AF2 accounts for faster T dissociation from full-length AR. T can nevertheless acquire DHT-like activity through an AR helix-10 H874Y prostate cancer mutation. The Tyr-874 mutant side chain mediates a new hydrogen bonding scheme from exterior helix-10 to backbone protein core helix-4 residue Tyr-739 to rescue T-induced AR activity by improving AF2 binding of FXXLF and LXXLL motifs. Greater AR AF2 activity by improved core helix interactions is supported by the effects of melanoma antigen gene protein-11, an AR coregulator that binds the AR FXXLF motif and targets AF2 for activation. We conclude that T is a weaker androgen than DHT because of less favorable T-dependent AR FXXLF and coactivator LXXLL motif interactions at AF2.

Survey of the year 2005: literature on applications of isothermal titration calorimetry

Isothermal titration calorimetry (ITC) can provide a full thermodynamic characterization of an interaction. Its usage does not suffer from constraints of molecular size, shape or chemical constitution. Neither is there any need for chemical modification or attachment to solid support. This ease of use has made it an invaluable instrumental resource and led to its appearance in many laboratories. Despite this, the value of the thermodynamic parameterization has, only quite recently, become widely appreciated. Although our understanding of the correlation between thermodynamic data and structural details continues to be somewhat naïve, a large number of publications have begun to improve the situation. In this overview of the literature for 2005, we have attempted to highlight works of interest and novelty. Furthermore, we draw attention to those works which we feel have provided a route to better analysis and increased our ability to understand the meaning of thermodynamic change on binding.

The mutations Lys 114 → Gln and Asp 126 → Asn disrupt an intersubunit salt bridge and convert Listeria innocua Dps into its natural mutant Listeria monocytogenes Dps. Effects on protein stability at Low pH

The stability of the dodecameric Listeria monocytogenes Dps has been compared with that of the Listeria innocua protein. The two proteins differ only in two amino acid residues that form an intersubunit salt-bridge in L. innocua Dps. This salt-bridge is replaced by a hydrogen bonding network in L. monocytogenes Dps as revealed by the X-ray crystal structure. The resistance to low pH and high temperature was assayed for both Dps proteins under equilibrium conditions and kinetically. Despite the identical equilibrium behavior, significant differences in the kinetic stability and activation energy of the unfolding process are apparent at pH 1.5. The higher stability of L. monocytogenes Dps has been accounted for in terms of the persistence of the hydrogen bonding network at this low pH value. In contrast, the salt-bridge between Lys 114 and Asp 126 characteristic of L. innocua Dps is most likely abolished due to protonation of Asp 126.

Water at biomolecular binding interfaces

Water molecules are often found at the binding interface of biomolecular complexes mediating the interaction between polar groups via hydrogen bonds, or simply filling space providing van der Waals interactions. Recent studies have demonstrated the importance of taking such water molecules into account in docking and binding affinity prediction. Here, we review the recent experimental and theoretical work aimed at quantifying the influence of interfacial water on the thermodynamic properties of binding. We highlight especially our recent results obtained by inhomogeneous fluid solvation theory in several systems and the prediction of the thermodynamic consequences of displacement of the bound water molecule by ligand modification. Finally, we discuss possible directions for further progress in this field.

CIRSE: a solvation energy estimator compatible with flexible protein docking and design applications

We present the Coordinate Internal Representation of Solvation Energy (CIRSE) for computing the solvation energy of protein configurations in terms of pairwise interactions between their atoms with analytic derivatives. Currently, CIRSE is trained to a Poisson/surface-area benchmark, but CIRSE is not meant to fit this benchmark exclusively. CIRSE predicts the overall solvation energy of protein structures from 331 NMR ensembles with 0.951+/-0.047 correlation and predicts relative solvation energy changes between members of individual ensembles with an accuracy of 15.8+/-9.6 kcal/mol. The energy of individual atoms in any of CIRSE's 17 types is predicted with at least 0.98 correlation. We apply the model in energy minimization, rotamer optimization, protein design, and protein docking applications. The CIRSE model shows some propensity to accumulate errors in energy minimization as well as rotamer optimization, but these errors are consistent enough that CIRSE correctly identifies the relative solvation energies of designed sequences as well as putative docked complexes. We analyze the errors accumulated by the CIRSE model during each type of simulation and suggest means of improving the model to be generally useful for all-atom simulations.

Anesthetic interaction with ketosteroid isomerase: insights from molecular dynamics simulations

The nature and the sites of interactions between anesthetic halothane and homodimeric Delta5-3-ketosteroid isomerase (KSI) are characterized by flexible ligand docking and confirmed by 1H-15N NMR. The dynamics consequence of halothane interaction and the implication of the dynamic changes to KSI function are studied by multiple 5-ns molecular dynamics simulations in the presence and absence of halothane. Both docking and MD simulations show that halothane prefer the amphiphilic dimeric interface to the hydrophobic active site of KSI. Halothane occupancy at the dimer interface disrupted the intersubunit hydrogen bonding formed either directly through side chains of polar residues or indirectly through the mediation of the interfacial water molecules. Moreover, in the presence of halothane, the exchange rate of the bound waters with bulk water was increased. Halothane perturbation to the dimer interface affected the overall flexibility of the active site. This action is likely to contribute to the halothane-induced reduction of the KSI activity. The allosteric halothane modulation of the dynamics-function relationship of KSI without direct competition at the enzymatic active sites may be generalized to offer a unifying explanation of anesthetic action on a diverse range of multidomain neuronal proteins that are potentially relevant to clinical general anesthesia.

The design, synthesis and evaluation of high affinity macrocyclic carbohydrate inhibitors

Carbohydrate-protein interactions have been investigated for a model system of a monoclonal antibody, SYA/J6, which binds a trisaccharide epitope of the O-polysaccharide of the Shigella flexneri variant Y lipopolysaccharide. The thermodynamics of binding for the methyl glycoside of the native trisaccharide epitope, Rha-Rha-GlcNAc () to SYA/J6 over a range of temperatures exhibits strong, linear enthalpy-entropy compensation and a negative heat capacity change (DeltaC(p)=-152 cal mol(-1) degree(-1)). At 293 K the free energy of association is the sum of favourable enthalpy and entropy contributions (DeltaH=-3.9 kcal mol(-1) and -TDeltaS=-2.9 kcal mol(-1)). Crystal structures for SYA/J6 Fab detailed the position of the native trisaccharide epitope, Rha-Rha-GlcNAc, and facilitated a strategy to design a tighter binding, low molecular weight ligand. This involved pre-organization of the native trisaccharide in its bound conformation by addition of intramolecular constraints (a beta-alanyl or glycinyl tether). ELISA measurements indicated that the glycinyl tethered trisaccharide was not an optimal candidate for further analysis, while microcalorimetry provided data showing that the beta-alanyl tethered trisaccharide displayed a 15-fold increase in affinity for SYA/J6. Tethering resulted in a favourable entropic contribution to binding, relative to the native trisaccharide (-TDeltaDeltaS=-1.2 kcal mol(-1)). Potential energy and dynamics calculations using the AMBER Plus force fields indicated that trisaccharide adopted a rigid conformation similar to that of the bound conformation of the native trisaccharide epitope. While this strategy resulted in modest free energy gains by minimizing losses due to conformational entropy, thermodynamic data are consistent with significant contributions from solvent reorganization.