Structure of a protein–detergent complex: the balance between detergent cohesion and binding

Despite the major interest in membrane proteins at functional, genomic, and therapeutic levels, their biochemical and structural study remains challenging, as they require, among other things, solubilization in detergent micelles. The complexity of this task derives from the dependence of membrane protein structure on their anisotropic environment, influenced by a delicate balance between many different physicochemical properties. To study such properties in a small protein–detergent complex, we used fluorescence measurements and molecular dynamics (MD) simulations on the transmembrane part of glycophorin A (GpAtm) solubilized in micelles of dihexanoylphosphatidylcholine (DHPC) detergent. Fluorescence measurements show that DHPC has limited ability to solubilize the peptide, while MD provides a possible molecular explanation for this. We observe that the detergent molecules are balanced between two different types of interactions: cohesive interactions between detergent molecules that hold the micelle together, and adhesive interactions with the peptide. While the cohesive interactions are detergent mediated, the adhesion to the peptide depends on the specific interactions between the hydrophobic parts of the detergent and the topography of the peptide dictated by the amino acids. The balance between these two parameters results in a certain frustration of the system and rather slow equilibration. These observations suggest how molecular properties of detergents could influence membrane protein stabilization and solubilization.     

Variation of protein binding cavity volume and ligand volume in protein–ligand complexes

We have systematically analyzed the variation of protein binding cavity volume of 200 protein–ligand complexes belonging to eight protein families. Wide variation in protein binding cavity volume for the same protein is observed on binding different ligands. Analysis of individual protein families shows high correlation between atom–atom interactions in binding site and ligand volume. This study implies the significance of protein flexibility in docking small molecule inhibitors on the basis of protein binding cavity volume with respect to ligand volume.     

The protein–protein interaction network of the human Sirtuin family

Protein–protein interaction networks are useful for studying human diseases and to look for possible health care through a holistic approach. Networks are playing an increasing and important role in the understanding of physiological processes such as homeostasis, signaling, spatial and temporal organizations, and pathological conditions. In this article we show the complex system of interactions determined by human Sirtuins (Sirt) largely involved in many metabolic processes as well as in different diseases. The Sirtuin family consists of seven homologous Sirt-s having structurally similar cores but different terminal segments, being rather variable in length and/or intrinsically disordered. Many studies have determined their cellular location as well as biological functions although molecular mechanisms through which they act are actually little known therefore, the aim of this work was to define, explore and understand the Sirtuin-related human interactome. As a first step, we have integrated the experimentally determined protein–protein interactions of the Sirtuin-family as well as their first and second neighbors to a Sirtuin-related sub-interactome. Our data showed that the second-neighbor network of Sirtuins encompasses 25% of the entire human interactome, and exhibits a scale-free degree distribution and interconnectedness among top degree nodes. Moreover, the Sirtuin sub interactome showed a modular structure around the core comprising mixed functions. Finally, we extracted from the Sirtuin sub-interactome subnets related to cancer, aging and post-translational modifications for information on key nodes and topological space of the subnets in the Sirt family network.     

GPCRs through the keyhole: the role of protein flexibility in ligand binding to β-adrenoceptors

G protein-coupled receptors (GPCRs) are proteins of pharmaceutical importance, with over 30% of all drugs in clinical use targeting them. Increasing numbers of X-ray crystal (XRC) structures of GPCRs offer a wealth of data relating to ligand binding. For the β-adrenoceptors (β-ARs), XRC structures are available for human β₂- and turkey β₁-subtypes, in complexes with a range of ligands. While these structures provide insight into the origins of ligand structure-activity relationships (SARs), questions remain. The ligands in all published complexed XRC structures lack extensive substitution, with no obvious way the ligand-binding site can accommodate β₁-AR-selective antagonists with extended side-chains para- to the common aryloxypropanolamine pharmacophore. Using standard computational docking tools with such ligands generally returns poses that fail to explain known SARs. Application of our Active Site Pressurisation modelling method to β-AR XRC structures and homology models, however, reveals a dynamic area in the ligand-binding pocket that, through minor changes in amino acid side chain orientations, opens a fissure between transmembrane helices H4 and H5, exposing intra-membrane space. This fissure, which we term the “keyhole”, is ideally located to accommodate extended moieties present in many high-affinity β₁-AR-selective ligands, allowing the rest of the ligand structure to adopt a canonical pose in the orthosteric binding site. We propose the keyhole may be a feature of both β₁- and β₂-ARs, but that subtle structural differences exist between the two, contributing to subtype-selectivity. This has consequences for the rational design of future generations of subtype-selective ligands for these therapeutically important targets.     

Impact of protein and ligand impurities on ITC-derived protein–ligand thermodynamics

Background

The thermodynamic characterization of protein–ligand interactions by isothermal titration calorimetry (ITC) is a powerful tool in drug design, giving valuable insight into the interaction driving forces. ITC is thought to require protein and ligand solutions of high quality, meaning both the absence of contaminants as well as accurately determined concentrations.

Methods

Ligands synthesized to deviating purity and protein of different pureness were titrated by ITC. Data curation was attempted also considering information from analytical techniques to correct stoichiometry.

Results and conclusions

We used trypsin and tRNA-guanine transglycosylase (TGT), together with high affinity ligands to investigate the effect of errors in protein concentration as well as the impact of ligand impurities on the apparent thermodynamics. We found that errors in protein concentration did not change the thermodynamic properties obtained significantly. However, most ligand impurities led to pronounced changes in binding enthalpy. If protein binding of the respective impurity is not expected, the actual ligand concentration was corrected for and the thus revised data compared to thermodynamic properties obtained with the respective pure ligand. Even in these cases, we observed differences in binding enthalpy of about 4 kJ ⋅ mol^− 1, which is considered significant.

General significance

Our results indicate that ligand purity is the critical parameter to monitor if accurate thermodynamic data of a protein–ligand complex are to be recorded. Furthermore, artificially changing fitting parameters to obtain a sound interaction stoichiometry in the presence of uncharacterized ligand impurities may lead to thermodynamic parameters significantly deviating from the accurate thermodynamic signature.     

Probing the structure of Leishmania donovani chagasi DHFR-TS: comparative protein modeling and protein–ligand interaction studies

Dihydrofolate reductase (DHFR) has been used successfully as a drug target in the area of anti-bacterial, anti-cancer and anti-malarial therapy. It also acts as a drug target for Leishmaniasis. Inhibition of DHFR leads to cell death through lack of thymine (nucleotide metabolism). Although the crystal structures of Leishmania major and Trypanosoma cruzi DHFR-thymidylate synthase (TS) have been resolved, to date there is no three-dimensional (3D)-structural information on DHFR-TS of Leishmania donovani chagasi, which causes visceral leishmaniasis. Our aim in this study was to model the 3D structure of L. donovani chagasi DHFR-TS, and to investigate the structural requirements for its inhibition. In this paper we describe a highly refined homology model of L. donovani chagasi DHFR-TS based on available crystallographic structures by using the Homology module of Insight II. Structural refinement and minimization of the generated L. donovani chagasi DHFR-TS model employed the Discover 3 module of Insight II and molecular dynamic simulations. The model was further validated through use of the PROCHECK, Verify_3D, PROSA, PSQS and ERRAT programs, which confirm that the model is reliable. Superimposition of the model structure with the templates L. major A chain, L. major B chain And T. cruzi A chain showed root mean square deviations of 0.69 Å, 0.71 Å and 1.11 Å, respectively. Docking analysis of the L. donovani chagasi DHFR-TS model with methotrexate enabled us to identify specific residues, viz. Val156, Val30, Lys95, Lys75 and Arg97, within the L. donovani chagasi DHFR-TS binding pocket, that play an important role in ligand or substrate binding. Docking studies clearly indicated that these five residues are important determinants for binding as they have strong hydrogen bonding interactions with the ligand.     

Identification of a protein–protein interaction between KCNE1 and the activation gate machinery of KCNQ1

KCNQ1 channels assemble with KCNE1 transmembrane (TM) peptides to form voltage-gated K⁺ channel complexes with slow activation gate opening. The cytoplasmic C-terminal domain that abuts the KCNE1 TM segment has been implicated in regulating KCNQ1 gating, yet its interaction with KCNQ1 has not been described. Here, we identified a protein–protein interaction between the KCNE1 C-terminal domain and the KCNQ1 S6 activation gate and S4–S5 linker. Using cysteine cross-linking, we biochemically screened over 300 cysteine pairs in the KCNQ1–KCNE1 complex and identified three residues in KCNQ1 (H363C, P369C, and I257C) that formed disulfide bonds with cysteine residues in the KCNE1 C-terminal domain. Statistical analysis of cross-link efficiency showed that H363C preferentially reacted with KCNE1 residues H73C, S74C, and D76C, whereas P369C showed preference for only D76C. Electrophysiological investigation of the mutant K⁺ channel complexes revealed that the KCNQ1 residue, H363C, formed cross-links not only with KCNE1 subunits, but also with neighboring KCNQ1 subunits in the complex. Cross-link formation involving the H363C residue was state dependent, primarily occurring when the KCNQ1–KCNE1 complex was closed. Based on these biochemical and electrophysiological data, we generated a closed-state model of the KCNQ1–KCNE1 cytoplasmic region where these protein–protein interactions are poised to slow activation gate opening.     

Estimation of binding parameters for the protein–protein interaction using a site-directed spin labeling and EPR spectroscopy

Sensitivity of the electron paramagnetic resonance (CW EPR) to molecular tumbling provides potential means for studying processes of molecular association. It uses spin-labeled macromolecules, whose CW EPR spectra may change upon binding to other macromolecules. When a spin-labeled molecule is mixed with its liganding partner, the EPR spectrum constitutes a linear combination of spectra of the bound and unbound ligand (as seen in our example of spin-labeled cytochrome c ₂ interacting with cytochrome bc ₁ complex). In principle, the fraction of each state can be extracted by the numerical decomposition of the spectrum; however, the accuracy of such decomposition may often be compromised by the lack of the spectrum of the fully bound ligand, imposed by the equilibrium nature of molecular association. To understand how this may affect the final estimation of the binding parameters, such as stoichiometry and affinity of the binding, a series of virtual titration experiments was conducted. Our non-linear regression analysis considered a case in which only a single class of binding sites exists, and a case in which classes of both specific and non-specific binding sites co-exist. The results indicate that in both models, the error due to the unknown admixture of the unbound ligand component in the EPR spectrum causes an overestimation of the bound fraction leading to the bias in the dissociation constant. At the same time, the stoichiometry of the binding remains relatively unaffected, which overall makes the decomposition of the EPR spectrum an attractive method for studying protein–protein interactions in equilibrium. Our theoretical treatment appears to be valid for any spectroscopic techniques dealing with overlapping spectra of free and bound component.     

Analysis of the aplyronine A-induced protein–protein interaction between actin and tubulin by surface plasmon resonance

The antitumor macrolide aplyronine A induces protein–protein interaction (PPI) between actin and tubulin to exert highly potent biological activities. The interactions and binding kinetics of these molecules were analyzed by the surface plasmon resonance with biotinylated aplyronines or tubulin as ligands. Strong binding was observed for tubulin and actin with immobilized aplyronine A. These PPIs were almost completely inhibited by one equivalent of either aplyronine A or C, or mycalolide B. In contrast, a non-competitive actin-depolymerizing agent, latrunculin A, highly accelerated their association. Significant binding was also observed for immobilized tubulin with an actin–aplyronine A complex, and the dissociation constant K_D was 1.84 μM. Our method could be used for the quantitative analysis of the PPIs between two polymerizing proteins stabilized with small agents.     

The role of the anticodon in the interaction between methionyl-tRNA synthetase and bacterial initiator tRNA?

Complementary and antiparallel oligonucleotides bind to exposed regions of the tRNA molecule. Aminoacylation in the presence of triplets has been used to determine the role of the anticodon in the interaction between methionyl-tRNA synthetase and initiator tRNA. ApUpG has no effect on the charging even when 70% of the tRNA is bound to the triplet, whereas in the presence of GpGpU which binds to the A-C-C sequence adjacent to the 3' terminal adenosine that fraction of the tRNA which is bound to the triplet is completely unavailable for charging. Hence the anticodon is probably not involved in a primary interaction while the A-C-C-A-OH clearly is. This conclusion is supported by the failure of the isolated anticodon loop and stem oligonucleotides to inhibit the aminoacylation reaction.     

Evaluation of small ligand–protein interaction by ligation reaction with DNA-modified ligand

A method for the evaluation of interactions between protein and ligand using DNA-modified ligands, including signal enhancement of the DNA ligation reactions, is described. For proof of principle, a DNA probe modified by biotin was used. Two DNA probes were prepared with complementary sticky-ends. While one DNA probe was modified at the 5′-end of the sticky-end, the other was not modified. The probes could be ligated together by T4 DNA ligase along the strand without biotin modification. However, in the presence of streptavidin or anti-biotin Fab, the ligation reaction joining the two probes could not occur on either strand.     

The fine-tuning of TRAF2–GSTP1-1 interaction: effect of ligand binding and in situ detection of the complex

We provide the first biochemical evidence of a direct interaction between the glutathione transferase P1-1 (GSTP1-1) and the TRAF domain of TNF receptor-associated factor 2 (TRAF2), and describe how ligand binding modulates such an equilibrium. The dissociation constant of the heterocomplex is K_d=0.3 μM; however the binding affinity strongly decreases when the active site of GSTP1-1 is occupied by the substrate GSH (K_d≥2.6 μM) or is inactivated by oxidation (K_d=1.7 μM). This indicates that GSTP1-1''s TRAF2-binding region involves the GSH-binding site. The GSTP1-1 inhibitor NBDHEX further decreases the complex''s binding affinity, as compared with when GSH is the only ligand; this suggests that the hydrophobic portion of the GSTP1-1 active site also contributes to the interaction. We therefore hypothesize that TRAF2 binding inactivates GSTP1-1; however, analysis of the data, using a model taking into account the dimeric nature of GSTP1-1, suggests that GSTP1-1 engages only one subunit in the complex, whereas the second subunit maintains the catalytic activity or binds to other proteins. We also analyzed GSTP1-1''s association with TRAF2 at the cellular level. The TRAF2–GSTP1-1 complex was constitutively present in U-2OS cells, but strongly decreased in S, G2 and M phases. Thus the interaction appears regulated in a cell cycle-dependent manner. The variations in the levels of individual proteins seem too limited to explain the complex''s drastic decline observed in cells progressing from the G0/G1 to the S–G2–M phases. Moreover, GSH''s intracellular content was so high that it always saturated GSTP1-1. Interestingly, the addition of NBDHEX maintains the TRAF2–GSTP1-1 complex at low levels, thus causing a prolonged cell cycle arrest in the G2/M phase. Overall, these findings suggest that a reversible sequestration of TRAF2 into the complex may be crucial for cell cycle progression and that multiple factors are involved in the fine-tuning of this interaction.     

The role of temperament in the relationship between morningness–eveningness and mood

Apart from differences in circadian phase position, individuals with different morningness–eveningness levels vary in many more characteristics. Particularly consistent relationships have been observed between morningness–eveningness and mood. Eveningness has been associated with disadvantageous mood, e.g. depressiveness in healthy individuals, and mood disorders. A concept of social jetlag suggests that evening subjects function in less advantageous environments due to discrepancies between internal and social time (societies promote morning-oriented functioning), which results in their lowered mood. Individual temperament, as defined by the Regulative Theory of Temperament (RTT), refers to the capacity of the human organism to meet environmental requirements – the greater the capacity, the less negative impact of external conditions. Thus, the aim of this study is to determine which RTT traits are linked to both morningness–eveningness and mood dimensions and to test whether they account for the relationship between morningness–eveningness and mood. A sample of 386 university students (267 female) aged between 19 and 47 (M?=?21.15, SD?=?4.23) years completed the University of Wales Institute of Science and Technology (UWIST) Mood Adjective Check List, Morningness–Eveningness Questionnaire and Formal Characteristics of Behaviour – Temperament Inventory. Analyses revealed lower endurance (EN) and higher emotional reactivity (ER) related to eveningness as well as to lower hedonic tone (HT), energetic arousal (EA) and to higher tense arousal (TA). Moreover, eveningness was associated with lower HT, EA and higher TA. Among RTT traits, EN was most strongly related to eveningness, and mediation analyses revealed that this temperamental trait fully mediated the relationship between eveningness and the three mood dimensions. The remaining RTT traits did not provide more explanation of the association between morningness–eveningness and mood than EN itself. If subjects did not differ in EN, the association between morningness–eveningness and mood was absent. EN is discussed as a protective factor against negative consequences of social jetlag and particularly lowered mood in evening individuals.     

Interactome of signaling networks in wheat: the protein–protein interaction between TaRAR1 and TaSGT1

RAR1 and SGT1 are required for development and disease resistance in plants. In many cases, RAR1 and SGT1 regulate the resistance (R)-gene-mediated defense signaling pathways. Lr21 is the first identified NBS-LRR-type R protein in wheat and is required for resistance to the leaf rust pathogen. The Lr21-mediated signaling pathways require the wheat homologs of RAR1, SGT1, and HSP90. However, the molecular mechanisms of the Lr21-mediated signaling networks remain unknown. Here I present the DNA and protein sequences of TaRAR1 and TaSGT1, and demonstrate for the first time a direct protein-protein interaction between them.     

Accuracy and precision of protein–ligand interaction kinetics determined from chemical shift titrations

NMR-monitored chemical shift titrations for the study of weak protein?Cligand interactions represent a rich source of information regarding thermodynamic parameters such as dissociation constants (K _D ) in the micro- to millimolar range, populations for the free and ligand-bound states, and the kinetics of interconversion between states, which are typically within the fast exchange regime on the NMR timescale. We recently developed two chemical shift titration methods wherein co-variation of the total protein and ligand concentrations gives increased precision for the K _D value of a 1:1 protein?Cligand interaction (Markin and Spyracopoulos in J Biomol NMR 53: 125?C138, 2012). In this study, we demonstrate that classical line shape analysis applied to a single set of ¹H?C¹⁵N 2D HSQC NMR spectra acquired using precise protein?Cligand chemical shift titration methods we developed, produces accurate and precise kinetic parameters such as the off-rate (k _off ). For experimentally determined kinetics in the fast exchange regime on the NMR timescale, k _off ?~?3,000?s^?1 in this work, the accuracy of classical line shape analysis was determined to be better than 5?% by conducting quantum mechanical NMR simulations of the chemical shift titration methods with the magnetic resonance toolkit GAMMA. Using Monte Carlo simulations, the experimental precision for k _off from line shape analysis of NMR spectra was determined to be 13?%, in agreement with the theoretical precision of 12?% from line shape analysis of the GAMMA simulations in the presence of noise and protein concentration errors. In addition, GAMMA simulations were employed to demonstrate that line shape analysis has the potential to provide reasonably accurate and precise k _off values over a wide range, from 100 to 15,000?s^?1. The validity of line shape analysis for k _off values approaching intermediate exchange (~100?s^?1), may be facilitated by more accurate K _D measurements from NMR-monitored chemical shift titrations, for which the dependence of K _D on the chemical shift difference (????) between free and bound states is extrapolated to ?????=?0. The demonstrated accuracy and precision for k _off will be valuable for the interpretation of biological kinetics in weakly interacting protein?Cprotein networks, where a small change in the magnitude of the underlying kinetics of a given pathway may lead to large changes in the associated downstream signaling cascade.     

Benzimidazole inhibitors of the protein kinase CHK2: Clarification of the binding mode by flexible side chain docking and protein–ligand crystallography

Two closely related binding modes have previously been proposed for the ATP-competitive benzimidazole class of checkpoint kinase 2 (CHK2) inhibitors; however, neither binding mode is entirely consistent with the reported SAR. Unconstrained rigid docking of benzimidazole ligands into representative CHK2 protein crystal structures reveals an alternative binding mode involving a water-mediated interaction with the hinge region; docking which incorporates protein side chain flexibility for selected residues in the ATP binding site resulted in a refinement of the water-mediated hinge binding mode that is consistent with observed SAR. The flexible docking results are in good agreement with the crystal structures of four exemplar benzimidazole ligands bound to CHK2 which unambiguously confirmed the binding mode of these inhibitors, including the water-mediated interaction with the hinge region, and which is significantly different from binding modes previously postulated in the literature.     

Influence of the protein conformation on the interaction between α-lactalbumin and dimyristoylphosphatidylcholine vesicles

α-Lactalbumin is a globular protein containing helical regions with highly amphiphathic character. In this work, the interaction between bovine α-lactalbumin and sonicated dimyristoylphosphatidylcholine vesicles has been compared in different circumstances which influence the protein conformation i.e., pH, ionic strength, decalcification, guanidine hydrochloride denaturation. Above the isoelectric point the interaction is mainly electrostatic; improved electrostatic interaction results in better contact with the apolar lipid phase. Below the isoelectric point, hydrophobic forces dominate the interaction and the vesicles are solubilized. The mode of interaction is not determined to a great extent by the demetallization of the protein. However, by a more explicit unfolding of the globular structure with guanidine hydrochloride, micellar complexes can be formed with the lipid, even at neutral pH. From this study it is obvious that the presence or capability for formation of helices with high amphipathic character is not a sufficient condition for lipid solubilization by a globular protein. Also, the capability of a globular protein to unfold its tertiary structure seems to be a prerequisite for its capability to lipid solubilization.     

Energetics of ligand recognition and self-association of bovine β-lactoglobulin: differences between variants A and B

An understanding of the interplay between structure and energetics is crucial for the optimization of modern protein engineering techniques. In this context, the study of natural isoforms is a subject of major interest, as it provides the scenario for analyzing mutations that have endured during biological evolution. In this study, we performed a comparative analysis of the ligand-recognition and homodimerization energetics of bovine β-lactoglobulin variants A (βlgA) and B (βlgB). These variants differ by only two amino-acid substitutions: 64th (Asp(A) → Gly(B)), which is fully exposed to the solvent, and 118th (Val(A) → Ala(B)), immersed in the hydrophobic core of the protein. Calorimetric measurements revealed significant enthalpic and entropic differences between the isoforms in both binding processes. A structural comparison suggests that a variation in the conformation of the loop C-D, induced by mutation Asp/Gly, could be responsible for the differences in ligand-binding energetics. While recognition of lauric acid was entropically driven, recognition of sodium dodecyl sulfate was both entropically and enthalpically driven, confirming the key role of the ligand polar moiety. Because of a more favorable enthalpy, the dimerization equilibrium constant of βlgB was larger than that of βlgA at room temperature, while the two dimers became similarly stable at 35 °C. The isoforms exchanged the same number of structural water molecules and protons and shared similar stereochemistry at the dimer interface. MD simulations revealed that the subunits of both variants become more flexible upon dimer formation. It is hypothesized that a larger increase of βlgA mobility could account for the dimerization energetic differences observed.