Geometry of nonbonded interactions involving planar groups in proteins

On the binding affinity of macromolecular interactions: daring to ask why proteins interact

Interactions between proteins are orchestrated in a precise and time-dependent manner, underlying cellular function. The binding affinity, defined as the strength of these interactions, is translated into physico-chemical terms in the dissociation constant (K_d), the latter being an experimental measure that determines whether an interaction will be formed in solution or not. Predicting binding affinity from structural models has been a matter of active research for more than 40 years because of its fundamental role in drug development. However, all available approaches are incapable of predicting the binding affinity of protein–protein complexes from coordinates alone. Here, we examine both theoretical and experimental limitations that complicate the derivation of structure–affinity relationships. Most work so far has concentrated on binary interactions. Systems of increased complexity are far from being understood. The main physico-chemical measure that relates to binding affinity is the buried surface area, but it does not hold for flexible complexes. For the latter, there must be a significant entropic contribution that will have to be approximated in the future. We foresee that any theoretical modelling of these interactions will have to follow an integrative approach considering the biology, chemistry and physics that underlie protein–protein recognition.

Biophysical and computational fragment-based approaches to targeting protein-protein interactions: applications in structure-guided drug discovery

Drug discovery has classically targeted the active sites of enzymes or ligand-binding sites of receptors and ion channels. In an attempt to improve selectivity of drug candidates, modulation of protein-protein interfaces (PPIs) of multiprotein complexes that mediate conformation or colocation of components of cell-regulatory pathways has become a focus of interest. However, PPIs in multiprotein systems continue to pose significant challenges, as they are generally large, flat and poor in distinguishing features, making the design of small molecule antagonists a difficult task. Nevertheless, encouragement has come from the recognition that a few amino acids - so-called hotspots - may contribute the majority of interaction-free energy. The challenges posed by protein-protein interactions have led to a wellspring of creative approaches, including proteomimetics, stapled α-helical peptides and a plethora of antibody inspired molecular designs. Here, we review a more generic approach: fragment-based drug discovery. Fragments allow novel areas of chemical space to be explored more efficiently, but the initial hits have low affinity. This means that they will not normally disrupt PPIs, unless they are tethered, an approach that has been pioneered by Wells and co-workers. An alternative fragment-based approach is to stabilise the uncomplexed components of the multiprotein system in solution and employ conventional fragment-based screening. Here, we describe the current knowledge of the structures and properties of protein-protein interactions and the small molecules that can modulate them. We then describe the use of sensitive biophysical methods - nuclear magnetic resonance, X-ray crystallography, surface plasmon resonance, differential scanning fluorimetry or isothermal calorimetry - to screen and validate fragment binding. Fragment hits can subsequently be evolved into larger molecules with higher affinity and potency. These may provide new leads for drug candidates that target protein-protein interactions and have therapeutic value.

Isoniazid and thioacetazone may exhibit antitubercular activity by binding directly with the active site of mycolic acid cyclopropane synthase: Hypothesis based on computational analysis

Isoniazid and thioacetazone are the two important antitubercular drugs. In case of thioacetazone it is established that it inhibits mycolic acid cyclopropane synthase but the exact binding site accounting for such inhibition is presently unknown. In case of isoniazid its action on the said enzyme is unexplored. In this work we have analyzed the binding of isoniazid and thioacetazone with mycolic acid cyclopropane synthase (CmaA1 and CmaA2) using tools of computational biology. We have observed that thioacetazone fits well at the active site of CmaA1 and CmaA2 while isoniazid binds at the active site of CmaA1 only. We have recommended experimental validation of such results. If such results are proved to be fact it will explore the exact binding site of thioacetazone and discover a new mechanism of anti-tubercular action of isoniazid.

Hypervalent nonbonded interactions of a divalent sulfur atom. Implications in protein architecture and the functions

In organic molecules a divalent sulfur atom sometimes adopts weak coordination to a proximate heteroatom (X). Such hypervalent nonbonded S···X interactions can control the molecular structure and chemical reactivity of organic molecules, as well as their assembly and packing in the solid state. In the last decade, similar hypervalent interactions have been demonstrated by statistical database analysis to be present in protein structures. In this review, weak interactions between a divalent sulfur atom and an oxygen or nitrogen atom in proteins are highlighted with several examples. S···O interactions in proteins showed obviously different structural features from those in organic molecules (i.e., π(o) → σ(s)* versus n(o) → σ(s)* directionality). The difference was ascribed to the HOMO of the amide group, which expands in the vertical direction (π(o)) rather than in the plane (n(o)). S···X interactions in four model proteins, phospholipase A₂ (PLA₂), ribonuclease A (RNase A), insulin, and lysozyme, have also been analyzed. The results suggested that S···X interactions would be important factors that control not only the three-dimensional structure of proteins but also their functions to some extent. Thus, S···X interactions will be useful tools for protein engineering and the ligand design.

Structural specializations of α(4)β(7), an integrin that mediates rolling adhesion

Electron microscopy and crystallography studies of α4β7 integrin reveal the mechanism by which this atypical integrin enables rolling adhesion prior to integrin activation.

The lymphocyte homing receptor integrin α₄β₇ is unusual for its ability to mediate both rolling and firm adhesion. α₄β₁ and α₄β₇ are targeted by therapeutics approved for multiple sclerosis and Crohn’s disease. Here, we show by electron microscopy and crystallography how two therapeutic Fabs, a small molecule (RO0505376), and mucosal adhesion molecule-1 (MAdCAM-1) bind α₄β₇. A long binding groove at the α₄–β₇ interface for immunoglobulin superfamily domains differs in shape from integrin pockets that bind Arg-Gly-Asp motifs. RO0505376 mimics an Ile/Leu-Asp motif in α₄ ligands, and orients differently from Arg-Gly-Asp mimics. A novel auxiliary residue at the metal ion–dependent adhesion site in α₄β₇ is essential for binding to MAdCAM-1 in Mg²⁺ yet swings away when RO0505376 binds. A novel intermediate conformation of the α₄β₇ headpiece binds MAdCAM-1 and supports rolling adhesion. Lack of induction of the open headpiece conformation by ligand binding enables rolling adhesion to persist until integrin activation is signaled.

Structural insights into antibody sequestering and neutralizing of Na+ channel α-type modulator from old world scorpion venom

The Old World scorpion Androctonus australis hector (Aah) produces one of the most lethal venoms for humans. Peptidic α-toxins AahI to AahIV are responsible for its potency, with AahII accounting for half of it. All four toxins are high affinity blockers of the fast inactivation phase of mammalian voltage-activated Na(+) channels. However, the high antigenic polymorphism of α-toxins prevents production of a polyvalent neutralizing antiserum, whereas the determinants dictating their trapping by neutralizing antibodies remain elusive. From an anti-AahII mAb, we generated an antigen binding fragment (Fab) with high affinity and selectivity for AahII and solved a 2.3 Å-resolution crystal structure of the complex. Sequestering of the C-terminal region of the bound toxin within a groove formed by the Fab combining loops is associated with a toxin orientation and main and side chain conformations that dictate the AahII antigenic specificity and efficient neutralization. From an anti-AahI mAb, we also preformed and crystallized a high affinity AahI-Fab complex. The 1.6 Å-resolution structure solved revealed a Fab molecule devoid of a bound AahI and with combining loops involved in packing interactions, denoting expulsion of the bound antigen upon crystal formation. Comparative analysis of the groove-like combining site of the toxin-bound anti-AahII Fab and planar combining surface of the unbound anti-AahI Fab along with complementary data from a flexible docking approach suggests occurrence of distinctive trapping orientations for the two toxins relative to their respective Fab. This study provides complementary templates for designing new molecules aimed at capturing Aah α-toxins and suitable for immunotherapy.

Residual structure in unfolded proteins

The denatured state ensemble (DSE) of unfolded proteins, once considered to be well-modeled by an energetically featureless random coil, is now well-known to contain flickering elements of residual structure. The position and nature of DSE residual structure may provide clues toward deciphering the protein folding code. This review focuses on recent advances in our understanding of the nature of DSE collapse under folding conditions, the quantification of the stability of residual structure in the DSE, the determination of the location and types of residues involved in thermodynamically significant residual structure and advances in detection of long-range interactions in the DSE.

Aromatic interactions at atom-to-atom contact and just beyond: a case study of protein interactions of NAD⁺/NADP⁺

We probed aromatic-protein interactions based on specificity of enrichment of protein residues across a contact-based cutoff. Thus, 155 protein-NAD(+)/NADP(+) complexes were analyzed for enrichments within 10Å of centroids of aromatic groups of the ligand when the residues were contacted and not contacted with the aromatic ligand. Specifically, neutral-adenine and cationic-nicotinamide groups of the oxidized coenzymes evoked interest to know whether the contrast of charge or the shared aromaticity will manifest in the enrichments across the cutoff. We found that when in contact, the enrichments are highly specific for nicotinamide and adenine-aromatic structures, and thus possibly complex in the basis, but when not in contact, they are generic for charge and aromaticity of the structures, and thus possibly specific in the basis. The order of enrichments over the contacted residues is Tyr>Cys>Thr>His>Asn>Ser>Met>Ile>Phe against nicotinamide-π(+) structure and Asp>Ile>Thr>His>Arg>Tyr>Gly>Val against adenine-π structure, while the order over the non-contacted residues is Trp>Gly>His>Asn>Cys>Met>Tyr>Ser>Thr>Phe against nicotinamide-π(+) structure and Asn>Thr>Ser>Gly>Cys>His>Val against adenine-π structure. Neutral Trp, His, Tyr, and Phe, but not cationic Arg, are thus the non-contacted residues enriched specifically against nicotinamide-π(+) structure, while Asn, Gly, Thr, Ser, and Cys are the non-contacted residues enriched generically against both the nicotinamide-π(+) and adenine-π aromatic structures. By analyzing the enriched groups in their geometric specificities, we found that, the enrichments against nicotinamide cation manifest the specificity expected of cation-π interaction and against nicotinamide- and adenine-aromatic groups manifest the specificity expected of dipole-π interaction. The cutoff-based method is proven valuable in probing protein-ligand interactions in the physics involved.

Modeling ADAMTS13-von Willebrand factor interaction: Implications for oxidative stress-related cardiovascular diseases and type 2A von Willebrand disease

The haemostatic potential of von Willebrand factor, a glycoprotein expressed by endothelial cells as ultra-large polymers (UL-vWF)(1), increases with its length, which in turn is regulated proteolytically by ADAMTS13, a zinc-metalloprotease selectively cleaving vWF at the Tyr1605-Met1606 bond. We have recently shown that in vitro oxidation of Met1606, under conditions mimicking those found in diseases characterized by high oxidative stress, severely impairs proteolysis by ADAMTS13, with a resulting pro-thrombotic effect caused by the accumulation of UL-vWF species. Conversely, Val1607Asp mutation, found in vWF from patients with type 2A von Willebrand disease, accelerates proteolysis of vWF, with a final hemorrhagic effect. Considering the physio-pathological importance of ADAMTS13-vWF interaction and the absence of experimental structural data, here we produced by homology modeling techniques a three-dimensional model of ADAMTS13 metalloprotease domain (M13). Thereafter, the vWF(1604-1607) peptide, containing the cleavable Tyr1605-Met1606 bond, was manually docked into the protease active site and the resulting model complex provided us key information for interpreting on structural grounds the variable effects that chemical modifications/mutations in vWF have on proteolysis by ADAMTS13.

Synthesis, activity, and structural analysis of novel α-hydroxytropolone inhibitors of human immunodeficiency virus reverse transcriptase-associated ribonuclease H

The α-hydroxytroplone, manicol (5,7-dihydroxy-2-isopropenyl-9-methyl-1,2,3,4-tetrahydro-benzocyclohepten-6-one), potently and specifically inhibits ribonuclease H (RNase H) activity of human immunodeficiency virus reverse transcriptase (HIV RT) in vitro. However, manicol was ineffective in reducing virus replication in culture. Ongoing efforts to improve the potency and specificity over the lead compound led us to synthesize 14 manicol derivatives that retain the divalent metal-chelating α-hydroxytropolone pharmacophore. These efforts were augmented by a high resolution structure of p66/p51 HIV-1 RT containing the nonnucleoside reverse transcriptase inhibitor (NNRTI), TMC278 and manicol in the DNA polymerase and RNase H active sites, respectively. We demonstrate here that several modified α-hydroxytropolones exhibit antiviral activity at noncytotoxic concentrations. Inclusion of RNase H active site mutants indicated that manicol analogues can occupy an additional site in or around the DNA polymerase catalytic center. Collectively, our studies will promote future structure-based design of improved α-hydroxytropolones to complement the NRTI and NNRTI currently in clinical use.

A compact review on the comparison of conventional and non-conventional interactions on the structural stability of therapeutic proteins

Therapeutic proteins carry out the most difficult tasks in living cells. They do so by interacting specifically with other molecules. This requires that they fold to a unique and more stable conformation. A prerequisite for comprehending the folding processes in their immense complexity entails a thorough understanding of many weak interactions. The purpose of this review is to systematically study the role of weak interactions such as cation-π, C-H......π, N-H......π and O-H......π, in the set of 49 therapeutic proteins. The importance of many of these interactions (for example, cationic residues interacting with π system) is revealed by the higher degree of conservation observed for them in protein structures. These interactions are mainly formed by long-range contacts and significant percentage of cation-π, C-H......π, N-H......π and O-H......π interacting residues had one or more stabilization centers. Further, a comparison of conventional and nonconventional interactions in the present data set unambiguously highlights the significance of these weak interactions in the structural stability of therapeutic proteins. We propose that the incorporation of the entirety of these interactions leads to a more complete description of the problem, and that this could provide new perspectives and new possible answers.

Aromatic rings in chemical and biological recognition: energetics and structures

This review describes a multidimensional treatment of molecular recognition phenomena involving aromatic rings in chemical and biological systems. It summarizes new results reported since the appearance of an earlier review in 2003 in host-guest chemistry, biological affinity assays and biostructural analysis, data base mining in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB), and advanced computational studies. Topics addressed are arene-arene, perfluoroarene-arene, S⋅⋅⋅aromatic, cation-π, and anion-π interactions, as well as hydrogen bonding to π systems. The generated knowledge benefits, in particular, structure-based hit-to-lead development and lead optimization both in the pharmaceutical and in the crop protection industry. It equally facilitates the development of new advanced materials and supramolecular systems, and should inspire further utilization of interactions with aromatic rings to control the stereochemical outcome of synthetic transformations.

Tryptophan residues: scarce in proteins but strong stabilizers of β-hairpin peptides

Tryptophan plays important roles in protein stability and recognition despite its scarcity in proteins. Except as fluorescent groups, they have been used rarely in peptide design. Nevertheless, Trp residues were crucial for the stability of some designed minimal proteins. In 2000, Trp-Trp pairs were shown to contribute more than any other hydrophobic interaction to the stability of β-hairpin peptides. Since then, Trp-Trp pairs have emerged as a paradigm for the design of stable β-hairpins, such as the Trpzip peptides. Here, we analyze the nature of the stabilizing capacity of Trp-Trp pairs by reviewing the β-hairpin peptides containing Trp-Trp pairs described up to now, the spectroscopic features and geometry of the Trp-Trp pairs, and their use as binding sites in β-hairpin peptides. To complete the overview, we briefly go through the other relevant β-hairpin stabilizing Trp-non-Trp interactions and illustrate the use of Trp in the design of short peptides adopting α-helical and mixed α/β motifs. This review is of interest in the field of rational design of proteins, peptides, peptidomimetics, and biomaterials.

Amino acids-nucleotides biomolecular recognition: from biological occurrence to affinity chromatography

In this review, the protein-DNA interactions are discussed considering different perspectives, and the biological occurrence of this interaction is explained at atomic level. The evaluation of the amino acid-nucleotide recognition has been investigated analysing datasets for predicting the association preferences and the geometry that favours the interaction. Based on this knowledge, an affinity chromatographic method was developed also exploiting this biological favoured contact. In fact, the implementation of this technique brings the possibility to apply the concept of molecular interactions to the development of new purification methodologies. In addition, the integration of the information recovered by all the different perspectives can bring new insights about some biological mechanisms, though not totally clarified.

Systematic mutation and thermodynamic analysis of central tyrosine pairs in polyspecific NKG2D receptor interactions

The homodimeric, activating natural killer cell receptor NKG2D interacts with multiple monomeric ligands polyspecifically, yet without central conformational flexibility. Crystal structures of multiple NKG2D-ligand interactions have identified the NKG2D tyrosine pair Tyr 152 and Tyr 199 as forming multiple specific but diverse interactions with MICA and related proteins. Here we systematically altered each tyrosine to tryptophan, phenylalanine, isoleucine, leucine, valine, serine, and alanine to measure the effect of mutation on affinity and thermodynamics for binding a range of similar ligands: MICA, the higher-affinity ligand MICB, and MICdesign, a high-affinity version of MICA that shares all NKG2D contact residues with MICA. Affinity and residue size were related: tryptophan could often substitute for tyrosine without loss of affinity; loss of the tyrosine hydroxyl through mutation to phenylalanine was tolerated more at position 152 than 199; and the smallest residues coincide with lowest affinities in general. NKG2D mutant van't Hoff binding thermodynamics generally show that substitution of other residues for tyrosine causes a moderate positive or flat van't Hoff slope consistent with moderate loss of binding enthalpy. One set of NKG2D mutations caused MICA to adopt a positive van't Hoff slope corresponding to absorption of heat, and another set caused MICB to adopt a negative slope of greater heat release than wild-type. MICdesign shared one example of the first set with MICA and one of the second set with MICB. When the NKG2D mutation affinities were arranged according to change in nonpolar surface area and compared to results from specific antibody-antigen and protein-peptide interactions, it was found that hydrophobic surface loss in NKG2D reduced binding affinity less than reported in the other contexts. The hydrophobic effect at the center of the NKG2D binding appears more similar to that at the periphery of an antibody-antigen binding site than at its center. Therefore the polyspecific NKG2D binding site is more tolerant of structural alteration in general than either an antibody-antigen or protein-peptide binding site, and this tolerance may adapt NKG2D to a broad range of protein surfaces with micromolar affinity.

[Molecular docking: role of intermolecular contacts in formation of complexes of proteins with nucleotides and peptides]

Knowledge of 3D-structure of protein-ligand complex is a major prerequisite for understanding the functioning mechanism of cellular proteins and membrane receptors. This is also of a great help in rational drug design projects. In the present paper we briefly review the molecular docking approaches used to predict possible orientation of a ligand in the protein binding site. The recent trends to improve the accuracy and efficiency of docking algorithms are demonstrated with the results obtained in Laboratory of Biomolecular Modeling. Particular attention is paid to protein-ligand hydrophobic and stacking interactions responsible for molecular recognition of ligand fragments. Such type of interactions are not always adequately represented in scoring criteria of docking applications that leads to mismatch in 3D-structure complexes predictions. That is why further inquiry of methods to account for these interactions is now the area of active research.

Propensities of aromatic amino acids versus leucine and proline to induce residual structure in the denatured-state ensemble of iso-1-cytochrome c

Histidine-heme loop formation in the denatured state of a protein is a sensitive means for probing residual structure under unfolding conditions. In this study, we use a host-guest approach to investigate the relative tendencies of different amino acids to promote residual structure under denaturing conditions. The host for this work is a 6-amino-acid insert of five alanines, followed by a lysine engineered immediately following a unique histidine near the N-terminus of yeast iso-1-cytochrome c. We substitute the fourth alanine in this sequence HAAAXAK (with X=Trp, Phe, Tyr, and Leu). The effects of proline are tested with substitutions at positions 1 and 5 in the insert (HPAAAAK and HAAAAPK, respectively). Thermodynamic studies on His-heme loop formation in 3 M guanidine hydrochloride reveal significant stabilization of residual structure by aromatic amino acids, particularly Trp and Phe, and minimal stabilization of residual structure by Leu. Prolines slightly disfavor His-heme loop formation, presumably due to enhanced chain stiffness. Kinetic studies reveal that much of the change in His-heme loop stability for the aromatic amino acids is caused by a slowdown in the rate of His-heme loop breakage, indicating that residual structure is preferentially stabilized in the closed-loop form of the denatured state.

Structure-activity relationships of 11 new congeners of the SMTP plasminogen modulator

The fungal metabolite Stachybotrys microspora triprenyl phenols (SMTPs) are small-molecule plasminogen modulators that enhance plasminogen activation. The SMTP molecule consists of a tricyclic γ-lactam moiety, an isoprene side-chain and an N-linked side-chain. Previous investigations have demonstrated that the N-linked side-chain is crucial for its activity. In this study, we have isolated 11 new SMTP congeners with a variety of N-linked side-chain structures, to investigate structure-activity relationships. Active compounds included congeners with a carboxyl or a sulfonic acid group in the N-linked side-chain, whereas not all the congeners with a carboxyl group were active. Of these congeners, that with methionine or tyrosine as the N-linked side-chain moiety was more active than that with an aliphatic amino acid. Congeners without ionizable group in the N-linked side-chain were essentially inactive.

Discriminating the native structure from decoys using scoring functions based on the residue packing in globular proteins

BACKGROUND

Setting the rules for the identification of a stable conformation of a protein is of utmost importance for the efficient generation of structures in computer simulation. For structure prediction, a considerable number of possible models are generated from which the best model has to be selected.

RESULTS

Two scoring functions, Rs and Rp, based on the consideration of packing of residues, which indicate if the conformation of an amino acid sequence is native-like, are presented. These are defined using the solvent accessible surface area (ASA) and the partner number (PN) (other residues that are within 4.5 A) of a particular residue. The two functions evaluate the deviation from the average packing properties (ASA or PN) of all residues in a polypeptide chain corresponding to a model of its three-dimensional structure. While simple in concept and computationally less intensive, both the functions are at least as efficient as any other energy functions in discriminating the native structure from decoys in a large number of standard decoy sets, as well as on models submitted for the targets of CASP7. Rs appears to be slightly more effective than Rp, as determined by the number of times the native structure possesses the minimum value for the function and its separation from the average value for the decoys.

CONCLUSION

Two parameters, Rs and Rp, are discussed that can very efficiently recognize the native fold for a sequence from an ensemble of decoy structures. Unlike many other algorithms that rely on the use of composite scoring function, these are based on a single parameter, viz., the accessible surface area (or the number of residues in contact), but still able to capture the essential attribute of the native fold.

An allosteric transition trapped in an intermediate state of a new kinesin-inhibitor complex

Human kinesin Eg5 plays an essential role in mitosis by separating duplicated centrosomes and establishing the bipolar spindle. Eg5 is an interesting drug target for the development of cancer chemotherapy, with seven inhibitors already in clinical trials. In the present paper, we report the crystal structure of the Eg5 motor domain complexed with a potent antimitotic inhibitor STLC (S-trityl-L-cysteine) to 2.0 A (1 A=0.1 nm) resolution. The Eg5-STLC complex crystallizes in space group P3(2) with three molecules per asymmetric unit. Two of the molecules reveal the final inhibitor-bound state of Eg5, whereby loop L5 has swung downwards to close the inhibitor-binding pocket, helix alpha4 has rotated by approx. 15 degrees and the neck-linker has adopted a docked conformation. The third molecule, however, revealed an unprecedented intermediate state, whereby local changes at the inhibitor-binding pocket have not propagated to structural changes at the switch II cluster and neck-linker. This provides structural evidence for the sequence of drug-induced conformational changes.

Solvent control of the soft angular potential in hydroxyl-pi hydrogen bonds: inertial orientational dynamics

Ultrafast polarization and wavelength selective IR pump-probe spectroscopy is used to measure the inertial and long time orientational dynamics of pi-hydrogen bonding complexes. Inertial orientational relaxation is sensitive to the angular potential associated with the hydrogen bond. The complexes studied are composed of phenol-OD (hydroxyl hydrogen replaced by deuterium) and various pi-base solvents with different electron donating or withdrawing substituents (chlorobenzene, bromobenzene, benzene, toluene, p-xylene, mesitylene, 1-pentyne). The different substituents provide experimental control of the hydrogen bond strength. The inertial orientational relaxation of the complexes, measured at the center frequency of each line, is independent of the hydrogen bond strength, demonstrating the insensitivity of the OD inertial dynamics, and therefore the H-bond angular potential, to the hydrogen bond strength. OD stretch absorption bands are inhomogeneously broadened through interactions with the solvent. The hydrogen bonding complexes all have similar wavelength dependent inertial orientational relaxation across their inhomogeneously broadened OD stretch absorption lines. The wavelength dependence of the inertial reorientation across each line arises because of a correlation between local solvent structure and the angular potential. These two results imply that local solvent structure acts as the controlling influence in determining the extent of inertial orientational relaxation, and therefore the angular potential, and that variation in the pi-hydrogen bond strength is of secondary importance.

Protein ionizable groups: pK values and their contribution to protein stability and solubility

The structure, stability, solubility, and function of proteins depend on their net charge and on the ionization state of the individual residues. Consequently, biochemists are interested in the pK values of the ionizable groups in proteins and how these pK values depend on their environment. We review what has been learned about pK values of ionizable groups in proteins from experimental studies and discuss the important contributions they make to protein stability and solubility.

Molecular basis of antimony treatment in leishmaniasis

Leishmaniasis is a disease that affects 2 million people and kills 70000 persons every year. It is caused by Leishmania species, which are human protozoan parasites of the trypanosomatidae family. Trypanosomatidae differ from the other eukaryotes in their specific redox metabolism because the glutathione/glutathione reductase system is replaced by the unique trypanothione/trypanothione reductase system. The current treatment of leishmaniasis relies mainly on antimonial drugs. The crystal structures of oxidized trypanothione reductase (TR) from Leishmania infantum and of the complex of reduced TR with NADPH and Sb(III), reported in this paper, disclose for the first time the molecular mechanism of action of antimonial drugs against the parasite. Sb(III), which is coordinated by the two redox-active catalytic cysteine residues (Cys52 and Cys57), one threonine residue (Thr335), and His461' of the 2-fold symmetry related subunit in the dimer, strongly inhibits TR activity. Because TR is essential for the parasite survival and virulence and it is absent in mammalian cells, these findings provide insights toward the design of new more affordable and less toxic drugs against Leishmaniasis.

Very short peptides with stable folds: building on the interrelationship of Trp/Trp, Trp/cation, and Trp/backbone-amide interaction geometries

By combining a favorable turn sequence with a turn flanking Trp/Trp interaction and a C-terminal H-bonding interaction between a backbone amide and an i-2 Trp ring, a particularly stable (DeltaG(U) > 7 kJ/mol) truncated hairpin, Ac-WI-(D-Pro-D-Asn)-KWTG-NH(2), results. In this construct and others with a W-(4-residue turn)-W motif in severely truncated hairpins, the C-terminal Trp is the edge residue in a well-defined face-to-edge (FtE) aryl/aryl interaction. Longer hairpins and those with six-residue turns retain the reversed "edge-to-face" (EtF) Trp/Trp geometry first observed for the trpzip peptides. Mutational studies suggest that the W-(4-residue turn)-W interaction provides at least 3 kJ/mol of stabilization in excess of that due to the greater beta-propensity of Trp. The pi-cation, and Trp/Gly-H(N) interactions have been defined. The latter can give rise to >3 ppm upfield shifts for the Gly-H(N) in -WX(n)G- units both in turns (n = 2) and at the C-termini (n = 1) of hairpins. Terminal YTG units result in somewhat smaller shifts (extrapolated to 2 ppm for 100% folding). In peptides with both the EtF and FtE W/W interaction geometries, Trp to Tyr mutations indicate that Trp is the preferred "face" residue in aryl/aryl pairings, presumably because of its greater pi basicity.

Accessibility and partner number of protein residues, their relationship and a webserver, ContPlot for their display

Background

Depending on chemical features residues have preferred locations – interior or exterior – in protein structures, which also determine how many other residues are found around them. The close packing of residues is the hallmark of protein interior and protein-protein interaction sites.

Results

The average values of accessible surface area (ASA) and partner number (PN, the number of other residues within a distance of 4.5 Å from any atom of a given residue) of different residues have been determined and a webserver, ContPlot has been designed to display these values (relative to the average values) along the protein sequence. This would be useful to visually identify residues that are densely packed, or those involved in protein-protein interactions. The skewness observed in the distribution of PNs is indicative of the hydrophobic or hydrophilic nature of the residue. The variation of ASA with PN can be analytically expressed in terms of a cubic equation. These equations (one for each residue) can be used to estimate the ASA of a polypeptide chain using the PNs of the individual residues in the structure.

Conclusion

The atom-based PNs (obtained by counting surrounding atoms) are highly correlated to the residue-based PN, indicating that the latter can adequately capture the atomic details of packing. The average values of ASA and PN associated with each residue should be useful in protein structure prediction or fold-recognition algorithm. ContPlot would provide a handy tool to assess the importance of a residue in the protein structure or interaction site.

Urea, but not guanidinium, destabilizes proteins by forming hydrogen bonds to the peptide group

The mechanism by which urea and guanidinium destabilize protein structure is controversial. We tested the possibility that these denaturants form hydrogen bonds with peptide groups by measuring their ability to block acid- and base-catalyzed peptide hydrogen exchange. The peptide hydrogen bonding found appears sufficient to explain the thermodynamic denaturing effect of urea. Results for guanidinium, however, are contrary to the expectation that it might H-bond. Evidently, urea and guanidinium, although structurally similar, denature proteins by different mechanisms.

Probing the lower size limit for protein-like fold stability: ten-residue microproteins with specific, rigid structures in water

Mutational optimization of two long-range interactions first observed in Ac-WINGKWT-NH2, (a) bifurcated H-bonding involving the threonine amide H(N) and side chain OH and the N-terminal acetyl carbonyl and (b) an H-bond between the entgegen-H(N) of the C-terminal amide and the indole ring of Trp6 that stabilizes a face-to-edge indole/indole interaction between Trp1 and Trp6, has afforded < or = 10 residue systems that yield a remarkably stable fold in water. Optimization was achieved by designing a hydrophobic cluster that sequesters these H-bonds from solvent exposure. The structures and extent of amide H/D exchange protection for CH3CH2CO-WI pGXWTGPS (p = D-Pro, X = Leu or Ile) were determined. These two systems are greater than 94% folded at 298 K (97.5% at 280 K) with melting temperatures > 75 degrees C. The fold appears to display minimal fluxionality; a well-converged NMR structure rationalizes all of the large structuring shifts observed, and we suggest that these designed constructs can be viewed as microproteins.

Efficiency of histidine-associating compounds for blocking the alzheimer's Abeta channel activity and cytotoxicity

The opening of the Alzheimer's Abeta channel permits the flux of calcium into the cell, thus critically disturbing intracellular ion homeostasis. Peptide segments that include the characteristic histidine (His) diad, His(13) and His(14), efficiently block the Abeta channel activity, blocking Abeta cytotoxicity. We hypothesize that the vicinal His-His peptides coordinate with the rings of His in the mouth of the pore, thus blocking the flow of calcium ions through the channel, with consequent blocking of Abeta cytotoxicity. To test this hypothesis, we studied Abeta ion channel activity and cytotoxicity after the addition of compounds that are known to have His association capacity, such as Ni(2+), imidazole, His, and a series of His-related compounds. All compounds were effective at blocking both Abeta channel and preventing Abeta cytotoxicity. The efficiency of protection of His-related compounds was correlated with the number of imidazole side chains in the blocker compounds. These data reinforce the premise that His residues within the Abeta channel sequence are in the pathway of ion flow. Additionally, the data confirm the contribution of the Abeta channel to the cytotoxicity of exogenous Abeta.