Comparison between computational alanine scanning and per-residue binding free energy decomposition for protein-protein association using MM-GBSA: application to the TCR-p-MHC complex.

Recognition by the T-cell receptor (TCR) of immunogenic peptides (p) presented by Class I major histocompatibility complexes (MHC) is the key event in the immune response against virus-infected cells or tumor cells. A study of the 2C TCR/SIYR/H-2K(b) system using a computational alanine scanning and a much faster binding free energy decomposition based on the Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) method is presented. The results show that the TCR-p-MHC binding free energy decomposition using this approach and including entropic terms provides a detailed and reliable description of the interactions between the molecules at an atomistic level. Comparison of the decomposition results with experimentally determined activity differences for alanine mutants yields a correlation of 0.67 when the entropy is neglected and 0.72 when the entropy is taken into account. Similarly, comparison of experimental activities with variations in binding free energies determined by computational alanine scanning yields correlations of 0.72 and 0.74 when the entropy is neglected or taken into account, respectively. Some key interactions for the TCR-p-MHC binding are analyzed and some possible side chains replacements are proposed in the context of TCR protein engineering. In addition, a comparison of the two theoretical approaches for estimating the role of each side chain in the complexation is given, and a new ad hoc approach to decompose the vibrational entropy term into atomic contributions, the linear decomposition of the vibrational entropy (LDVE), is introduced. The latter allows the rapid calculation of the entropic contribution of interesting side chains to the binding. This new method is based on the idea that the most important contributions to the vibrational entropy of a molecule originate from residues that contribute most to the vibrational amplitude of the normal modes. The LDVE approach is shown to provide results very similar to those of the exact but highly computationally demanding method.

Functional dynamics of PDZ binding domains: a normal-mode analysis.

Postsynaptic density-95/disks large/zonula occludens-1 (PDZ) domains are relatively small (80-120 residues) protein binding modules central in the organization of receptor clusters and in the association of cellular proteins. Their main function is to bind C-terminals of selected proteins that are recognized through specific amino acids in their carboxyl end. Binding is associated with a deformation of the PDZ native structure and is responsible for dynamical changes in regions not in direct contact with the target. We investigate how this deformation is related to the harmonic dynamics of the PDZ structure and show that one low-frequency collective normal mode, characterized by the concerted movements of different secondary structures, is involved in the binding process. Our results suggest that even minimal structural changes are responsible for communication between distant regions of the protein, in agreement with recent NMR experiments. Thus, PDZ domains are a very clear example of how collective normal modes are able to characterize the relation between function and dynamics of proteins, and to provide indications on the precursors of binding/unbinding events.

Versatile binding properties of a di-2-pyridyl ketone nicotinoylhydrazone ligand: Crystal structure of a Cu(II) complex

Six new copper complexes of di-2-pyridyl ketone nicotinoylhydrazone (HDKN) have been synthesized. The complexes have been characterized by a variety of spectroscopic techniques and the structure of [Cu(DKN)2]·H2O has been determined by single crystal X-ray diffraction. The compound [Cu(DKN)2]·H2O crystallized in the monoclinic space group P21 and has a distorted octahedral geometry. The IR spectra revealed the presence of variable modes of chelation for the investigated ligand. The EPR spectra of compounds [Cu2(DKN)2( -N3)2] and [Cu2(DKN)2( -NCS)2] in polycrystalline state suggest a dimeric structure as they exhibited a half field signal, which indicate the presence of a weak interaction between two Cu(II) ions in these complexes

Consequences of N,C,N '- and C,N,N '-coordination modes on electronic and photophysical properties of cyclometalated aryl ruthenium(II) complexes

The effects of isoelectronic replacement of a neutral nitrogen donor atom by an anionic carbon atom in terpyridine ruthenium(II) complexes on the electronic and photophysical properties of the resulting N,C,N'- and C,N,N'-cyclometalated aryl ruthenium(II) complexes were investigated. To this end, a series of complexes was prepared either with ligands containing exclusively nitrogen donor atoms, that is, [Ru(R-1-tpy)(R-2-tpy)](2+) (R-1, R-2 = H, CO2Et), or bearing either one N,C,N'- or C,N,N'-cyclometalated ligand and one tpy ligand, that is, [Ru(R-1-(NCN)-C-Lambda-N-Lambda)(R-2-tpy)](+) and [Ru(R-1-(CNN)-N-Lambda-N-Lambda)(R-2-tpy)](+), respectively. Single-crystal X-ray structure determinations showed that cyclometalation does not significantly alter the overall geometry of the complexes but does change the bond lengths around the ruthenium(II) center, especially the nitrogen-to-ruthenium bond length trans to the carbanion. Substitution of either of the ligands with electron-withdrawing ester functionalities fine-tuned the electronic properties and resulted in the presence of an IR probe. Using trends obtained from redox potentials, emission energies, IR spectroelectrochemical responses, and the character of the lowest unoccupied molecular orbitals from DFT studies, it is shown that the first reduction process and luminescence are associated with the ester-substituted C,N,N'-cyclometalated ligand in [Ru(EtO2C-(CNN)-N-Lambda-N-Lambda)(tpy)](+). Cyclometalation in an N,C,N'-bonding motif changed the energetic order of the ruthenium d(zx), d(yz), and d(xy) orbitals. The red-shifted absorption in the N,C,N'-cyclometalated complexes is assigned to MLCT transitions to the tpy ligand. The red shift observed upon introduction of the ester moiety is associated with an increase in intensity of low-energy transitions, rather than a red shift of the main transition. Cyclometalation in the C,N,N'-binding motif also red-shifts the absorption, but the corresponding transition is associated with both ligand types. Luminescence of the cyclometalated complexes is relatively independent of the mode of cyclometalation, obeying the energy gap law within each individual series.

Crystal structures of Λ-[Ru(phen)2dppz]2+ with oligonucleotides containing TA/TA and AT/AT steps show two intercalation modes

The ruthenium complex [Ru(phen)2(dppz)] (where phen is a phenanthroline and dppz a dipyridyl–phenazine ligand) is known as a ‘light switch’ complex because its luminescence in solution is significantly enhanced in the presence of DNA. This property is poised to serve in diagnostic and therapeutic applications, but its binding mode with DNA needs to be elucidated further. Here, we describe the crystal structures of the L enantiomer bound to two oligonucleotide duplexes. The dppz ligand intercalates symmetrically and perpendicularly from the minor groove of the d(CCGGTACCGG)2 duplex at the central TA/TA step, but not at the central AT/AT step of d(CCGGATCCGG)2. In both structures, however, a second ruthenium complex links the duplexes through the combination of a shallower angled intercalation into the C1C2/G9G10 step at the end of the duplex, and semi-intercalation into the G3G4 step of an adjacent duplex. The TA/TA specificity of the perpendicular intercalation arises from the packing of phenanthroline ligands against the adenosine residue.

Identification of a novel ligand binding motif in the transthyretin channel

The design of therapeutic compounds targeting transthyretin (TTR) is challenging due to the low specificity of interaction in the hormone binding site. Such feature is highlighted by the interactions of TTR with diclofenac, a compound with high affinity for TTR, in two dissimilar modes, as evidenced by crystal structure of the complex. We report here structural analysis of the interactions of TTR with two small molecules, 1-amino-5-naphthalene sulfonate (1,5-AmNS) and 1-anilino-8-naphthalene sulfonate (1,8-ANS). Crystal structure of TTR: 1,8-ANS complex reveals a peculiar interaction, through the stacking of the naphthalene ring between the side-chain of Lys15 and Leu17. The sulfonate moiety provides additional interaction with Lys15` and a water-mediated hydrogen bond with Thr119`. The uniqueness of this mode of ligand recognition is corroborated by the crystal structure of TTR in complex with the weak analogue 1,5-AmNS, the binding of which is driven mainly by hydrophobic partition and one electrostatic interaction between the sulfonate group and the Lys15. The ligand binding motif unraveled by 1,8-ANS may open new possibilities to treat TTR amyloid diseases by the elucidation of novel candidates for a more specific pharmacophoric pattern. (C) 2009 Published by Elsevier Ltd.

Flexibility and inhibitor binding in Cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti-cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo- or sulfate-bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer-generated conformational ensembles shown here indicate that the last 30-40 residues in the C-terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C-terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B-inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins.

The Role of Nanometer-Scaled Ligand Patterns in Polyvalent Binding by Large Mannan-Binding Lectin Oligomers

Mannan-binding lectin (MBL) is an important protein of the innate immune system and protects the body against infection through opsonization and activation of the complement system on surfaces with an appropriate presentation of carbohydrate ligands. The quaternary structure of human MBL is built from oligomerization of structural units into polydisperse complexes typically with three to eight structural units, each containing three lectin domains. Insight into the connection between the structure and ligand-binding properties of these oligomers has been lacking. In this article, we present an analysis of the binding to neoglycoprotein-coated surfaces by size-fractionated human MBL oligomers studied with small-angle x-ray scattering and surface plasmon resonance spectroscopy. The MBL oligomers bound to these surfaces mainly in two modes, with dissociation constants in the micro to nanomolar order. The binding kinetics were markedly influenced by both the density of ligands and the number of ligand-binding domains in the oligomers. These findings demonstrated that the MBL-binding kinetics are critically dependent on structural characteristics on the nanometer scale, both with regard to the dimensions of the oligomer, as well as the ligand presentation on surfaces. Therefore, our work suggested that the surface binding of MBL involves recognition of patterns with dimensions on the order of 10-20 nm. The recent understanding that the surfaces of many microbes are organized with structural features on the nanometer scale suggests that these properties of MBL ligand recognition potentially constitute an important part of the pattern-recognition ability of these polyvalent oligomers. The Journal of Immunology, 2012, 188: 1292-1306.

Molecular architecture of fur binding to iron-induced and - repressed genes in Helicobacter pylori

The ferric uptake regulator protein Fur regulates iron-dependent gene expression in bacteria. In the human pathogen Helicobacter pylori, Fur has been shown to regulate iron-induced and iron-repressed genes. Herein we investigate the molecular mechanisms that control this differential iron-responsive Fur regulation. Hydroxyl radical footprinting showed that Fur has different binding architectures, which characterize distinct operator typologies. On operators recognized with higher affinity by holo-Fur, the protein binds to a continuous AT-rich stretch of about 20 bp, displaying an extended protection pattern. This is indicative of protein wrapping around the DNA helix. DNA binding interference assays with the minor groove binding drug distamycin A, point out that the recognition of the holo-operators occurs through the minor groove of the DNA. By contrast, on the apo-operators, Fur binds primarily to thymine dimers within a newly identified TCATTn10TT consensus element, indicative of Fur binding to one side of the DNA, in the major groove of the double helix. Reconstitution of the TCATTn10TT motif within a holo-operator results in a feature binding swap from an holo-Fur- to an apo-Fur-recognized operator, affecting both affinity and binding architecture of Fur, and conferring apo-Fur repression features in vivo. Size exclusion chromatography indicated that Fur is a dimer in solution. However, in the presence of divalent metal ions the protein is able to multimerize. Accordingly, apo-Fur binds DNA as a dimer in gel shift assays, while in presence of iron, higher order complexes are formed. Stoichiometric Ferguson analysis indicates that these complexes correspond to one or two Fur tetramers, each bound to an operator element. Together these data suggest that the apo- and holo-Fur repression mechanisms apparently rely on two distinctive modes of operator-recognition, involving respectively the readout of a specific nucleotide consensus motif in the major groove for apo-operators, and the recognition of AT-rich stretches in the minor groove for holo-operators, whereas the iron-responsive binding affinity is controlled through metal-dependent shaping of the protein structure in order to match preferentially the major or the minor groove.

Computational Study of the Fe(CN)(2)CO Cofactor and Its Binding to HypC Protein.

In the intricate maturation process of [NiFe]-hydrogenases, the Fe(CN)2CO cofactor is first assembled in a HypCD complex with iron coordinated by cysteines from both proteins and CO is added after ligation of cyanides. The small accessory protein HypC is known to play a role in delivering the cofactor needed for assembling the hydrogenase active site. However, the chemical nature of the Fe(CN)2CO moiety and the stability of the cofactor–HypC complex are open questions. In this work, we address geometries, properties, and the nature of bonding of all chemical species involved in formation and binding of the cofactor by means of quantum calculations. We also study the influence of environmental effects and binding to cysteines on vibrational frequencies of stretching modes of CO and CN used to detect the presence of Fe(CN)2CO. Carbon monoxide is found to be much more sensitive to sulfur binding and the polarity of the medium than cyanides. The stability of the HypC–cofactor complex is analyzed by means of molecular dynamics simulation of cofactor-free and cofactor-bound forms of HypC. The results show that HypC is stable enough to carry the cofactor, but since its binding cysteine is located at the N-terminal unstructured tail, it presents large motions in solution, which suggests the need for a guiding interaction to achieve delivery of the cofactor.

Structural factors controlling ligand binding to myoglobin: A kinetic hole-burning study

Using temperature-derivative spectroscopy in the temperature range below 100 K, we have studied the dependence of the Soret band on the recombination barrier in sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation at 12 K. The spectra were separated into contributions from the photodissociated species, Mb*CO, and CO-bound myoglobin. The line shapes of the Soret bands of both photolyzed and liganded myoglobin were analyzed with a model that takes into account the homogeneous bandwidth, coupling of the electronic transition to vibrational modes, and static conformational heterogeneity. The analysis yields correlations between the activation enthalpy for rebinding and the model parameters that characterize the homogeneous subensembles within the conformationally heterogeneous ensemble. Such couplings between spectral and functional parameters arise when they both originate from a common structural coordinate. This effect is frequently denoted as “kinetic hole burning.” The study of these correlations gives direct insights into the structure–function relationship in proteins. On the basis of earlier work that assigned spectral parameters to geometric properties of the heme, the connections with the heme geometry are discussed. We show that two separate structural coordinates influence the Soret line shape, but only one of the two is coupled to the enthalpy barrier for rebinding. We give evidence that this coordinate, contrary to widespread belief, is not the iron displacement from the mean heme plane.

Identification of a family of low-affinity insulin-like growth factor binding proteins (IGFBPs): Characterization of connective tissue growth factor as a member of the IGFBP superfamily

The insulin-like growth factor (IGF) binding proteins (IGFBPs) modulate the actions of the insulin-like growth factors in endocrine, paracrine, and autocrine settings. Additionally, some IGFBPs appear to exhibit biological effects that are IGF independent. The six high-affinity IGFBPs that have been characterized to date exhibit 40–60% amino acid sequence identity overall, with the most conserved sequences in their NH2 and COOH termini. We have recently demonstrated that the product of the mac25/IGFBP-7 gene, which shows significant conservation in the NH2 terminus, including an “IGFBP motif” (GCGCCXXC), exhibits low-affinity IGF binding. The closely related mammalian genes connective tissue growth factor (CTGF) gene, nov, and cyr61 encode secreted proteins that also contain the conserved sequences and IGFBP motifs in their NH2 termini. To ascertain if these genes, along with mac25/IGFBP-7, encode a family of low-affinity IGFBPs, we assessed the IGF binding characteristics of recombinant human CTGF (rhCTGF). The ability of baculovirus-synthesized rhCTGF to bind IGFs was demonstrated by Western ligand blotting, affinity cross-linking, and competitive affinity binding assays using 125I-labeled IGF-I or IGF-II and unlabeled IGFs. CTGF, like mac25/IGFBP-7, specifically binds IGFs, although with relatively low affinity. On the basis of these data, we propose that CTGF represents another member of the IGFBP family (IGFBP-8) and that the CTGF gene, mac25/IGFBP-7, nov, and cyr61 are members of a family of low-affinity IGFBP genes. These genes, along with those encoding the high-affinity IGFBPs 1–6, together constitute an IGFBP superfamily whose products function in IGF-dependent or IGF-independent modes to regulate normal and neoplastic cell growth.

A transcriptional activating region with two contrasting modes of protein interaction

A C-terminal segment of the yeast activator Gal4 manifests two functions: When tethered to DNA, it elicits gene activation, and it binds the inhibitor Gal80. Here we examine the effects on these two functions of cysteine and proline substitutions. We find that, although certain cysteine substitutions diminish interaction with Gal80, those substitutions have little effect on the activating function in vivo and interaction with TATA box-binding protein (TBP) in vitro. Proline substitutions introduced near residues critical for Gal80 binding abolish that interaction but once again have no effect on the activating function. Crosslinking experiments show that a defined position in the activating peptide is in close proximity to TBP and Gal80 in the two separate reactions and show that binding of the inhibitor blocks binding to TBP. Thus, the same stretch of amino acids are involved in two quite different protein–protein interactions: binding to Gal80, which depends on a precise sequence and the formation of a defined secondary structure, or interactions with the transcriptional machinery in vivo, which are not impaired by perturbations of either sequence or structure.

Identification of a second aryl phosphate-binding site in protein-tyrosine phosphatase 1B: A paradigm for inhibitor design

The structure of the catalytically inactive mutant (C215S) of the human protein-tyrosine phosphatase 1B (PTP1B) has been solved to high resolution in two complexes. In the first, crystals were grown in the presence of bis-(para-phosphophenyl) methane (BPPM), a synthetic high-affinity low-molecular weight nonpeptidic substrate (Km = 16 μM), and the structure was refined to an R-factor of 18.2% at 1.9 Å resolution. In the second, crystals were grown in a saturating concentration of phosphotyrosine (pTyr), and the structure was refined to an R-factor of 18.1% at 1.85 Å. Difference Fourier maps showed that BPPM binds PTP1B in two mutually exclusive modes, one in which it occupies the canonical pTyr-binding site (the active site), and another in which a phosphophenyl moiety interacts with a set of residues not previously observed to bind aryl phosphates. The identification of a second pTyr molecule at the same site in the PTP1B/C215S–pTyr complex confirms that these residues constitute a low-affinity noncatalytic aryl phosphate-binding site. Identification of a second aryl phosphate binding site adjacent to the active site provides a paradigm for the design of tight-binding, highly specific PTP1B inhibitors that can span both the active site and the adjacent noncatalytic site. This design can be achieved by tethering together two small ligands that are individually targeted to the active site and the proximal noncatalytic site.

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)2 and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.