Cucurbitacins from Ecballium elaterium juice increase the binding of bilirubin and ibuprofen to albumin in human plasma

Ecballium elaterium, a medicinal plant, whose fruit juice is used for the treatment of jaundice in folk medicine, has been reported as being capable of decreasing bilirubinemia in animals with jaundice [H.H. Elayan, M.N. Garaibeh, S.M. Zmeili, S.A. Salhab, Effects of Ecballium elaterium juice on serum bilirubin concentration in male rats, Int. J. Crude Drug Res. 27 (1989) 227-234]. The aim of this study is to identify the Ecballium elaterium components, which are able to modify the binding of bilirubin to albumin. The juice is fiber-free but contains proteins, lipids, sugars, and minerals. The extract of the juice, analyzed by liquid chromatography-electrospray ionization-mass spectrometry (LC-ESI-MS), contains cucurbitacins (Cuc) B, D, E, and I as well as several glycosylated compounds. Human plasma containing no or serial concentrations of Ecballium elaterium components were prepared and the direct bilirubin (DB) and total bilirubin (TB) were determined by the Jendrassik and Grof method. Our results showed that Cuc D, E, and B decreased the levels of DB and TB in plasma, while Cuc I, glycosyl derivatives, and proteins of the juice did not modify the bilirubin levels. The binding of domain specific ligands to HSA, bilirubin (domain IIA), and ibuprofen (domain IIIA), were studied in the absence and presence of Cuc D, E, and I, by fluorescence spectroscopy. The values of binding constant K(a) and binding site number n, determined by Scatchard method, increased for the both ligands only in the presence of Cuc E and D. Cuc I decreased slightly the K(a) of ibuprofen, suggesting an interaction with the domain IIIA of the protein. As a conclusion, Cuc E, D, and B produce rearrangement in the structure of albumin leading to increase the binding of domain specific ligands, ibuprofen and bilirubin.     

Quantitative determination of bilirubin by inhibition of chemiluminescence from lucigenin.

Bilirubin is a metabolic breakdown product of blood haem, of great biological and diagnostic importance. A new chemiluminescence (CL) method has been developed for the quantification of bilirubin. The method is combined with the flow injection analysis (FIA) technique and based on the inhibition effect of bilirubin on the CL from the lucigenin-hydrogen peroxide system in an alkaline medium. Under the optimum conditions, the decreased CL intensity was proportional to the concentration of bilirubin, in the range 0.0585-58.47 microg/mL. The detection limit estimated from the calibration graph was about 7.8826 ng/mL. The relative standard deviation (RSD) of 10 parallel measurements (1 x 10(-4) mol/L bilirubin) was 2.5%. Recoveries of bilirubin were found to fall in the range 94-97.5% using control sera. The method is interference-free, fast and easy to carry out.     

ATR-FTIR study of the protonation states of the Glu residue in the multicopper oxidases, CueO and bilirubin oxidase

Redox-induced protonation state changes of the Glu residue in the multicopper oxidases, CueO and bilirubin oxidase (BO), were studied by attenuated total reflectance-Fourier transform infrared spectroscopy. By monitoring IR bands of the carboxylic acid CO stretch in the wild-type and Glu-to-Gln mutant enzymes the Glu506 of CueO (Glu463 of BO) was found to be unprotonated in the oxidised and protonated in the reduced forms. The results provided direct evidence for proton uptake by the Glu, suggesting it plays a key role in the proton donation to the activated oxygen species in the catalytic cycle.     

The water-dependence of the catalytic activity of bilirubin oxidase suspensions in low-water systems.

We investigated the enzymic activity of bilirubin oxidase when it is suspended as a lyophilized powder in a low-water system. The enzyme required buffer salts and a source of water to show activity. This study investigated the complete range of water thermodynamic activity (a(w)) by combining the use of salt hydrates and two-phase systems with concentrated solutes in the aqueous phase. When free water was added, activity reached a maximum at a defined water content, but this maximum increased with buffer content, suggesting that there was competition for water with the buffer salts from which the enzyme was lyophilized. Alternatively, a range of salt hydrates was used, each able to fix the water activity (a(w)) at a different value. By providing water to the organic solvent phase in this way, the dependency of enzyme activity upon a(w) was investigated and shown to be independent of buffer concentration. However, the optimum a(w) was uncertain because the available a(w) range for salt hydrates is < or = 0.90. Investigation of the remaining water activity range was made possible by using an a(w) depressor (sorbitol) to lower the a(w) of a two-phase system. The optimum a(w) for the bilirubin oxidase activity in this two-phase system was a(w) = 0.936, independent of buffer concentration. The study therefore confirmed the need to control the water 'available' to low-water systems and the dependence of enzyme activity on water thermodynamic activity (a(w)) not water content.     

The use of nonhuman primates in the study of bilirubin metabolism and bile secretion

Close similarities exist in the chemical structure of fetal and adult bile pigments and mechanisms for hepatic bilirubin conjugation and transport in human and nonhuman primates. Newborn monkeys, unlike other animals, develop physiologic hyperbilirubinemia and are ideal animal models in which to investigate this developmental human disorder. Aberrations in bile pigment metabolism and compartmentalization which occur in hemolytic, hepatocellular, and obstructive diseases, closely resemble those reported in counterpart syndromes in man. Similarities also exist in the mechanisms responsible for the secretion of canalicular and ductular/ductal bile and the composition of bile. Nonhuman primates offer distinct advantages for comparative studies concerning bile pigment metabolism and the secretion of bile.     

Retardation of the unfolding process by single N-glycosylation of ribonuclease A based on molecular dynamics simulations

The conformational characteristics of glycosylated- and unglycosylated bovine pancreatic ribonuclease A (RNaseA) were traced with unfolding molecular dynamics simulations using CHARMM program at 470 K. The glycosylated RNase (Glc_RNase) possesses nearly identical protein structure with RNaseA, differing only by presence of a single acetylglucosamine residue N-linked to Asn34 in the RNaseA. Attaching of monomeric N-acetylglucosamine residue to the Asn34 in RNaseA resulted in a change of denaturing process of Glc_RNase. Simulations showed that the unfolding of RNaseA involved significant weakening of nonlocal interactions whereas the glycosylation led Glc_RNase to preserve the nonlocal interactions even in its denatured form. Even in simulations over 8 ns at 470 K, Glc_RNase remained relatively stable as a less denatured conformation. However, conformation of RNaseA was changed to a fully unfolded state before 3 ns of the simulations at 470 K. This difference was due to fact that formation of hydrogen bond bridges and nonlocal contacts induced by the attached N-acetylglucosamine of Glc_RNase showing in the unfolding simulations. These high-temperature unfolding MD simulations provided a theoretical basis for the previous experimental work in which Glc_RNase showed slower unfolding kinetics compared with unglycosylated RNaseA, suggesting that single N-glycosylation induced retardation of unfolding process of the ribonuclease protein.     

Theoretical evidence of a salt bridge disruption as the initiating process for the alpha1d-adrenergic receptor activation: a molecular dynamics and docking study

This study reports the building of the three-dimensional structure of the rat alpha1d-adrenergic receptor through a topology approach based on the structure of the rhodopsin receptor from cryoelectron microscopy. The validity and reliability of the receptor model were assessed through exhaustive molecular dynamics and docking studies. Some interesting ligand-receptor interactions were identified along with significant differences between the binding mode of agonists and antagonists. The importance of the disruption of a salt bridge as a possible initial event leading to receptor activation is discussed on the basis of data from mutagenesis and molecular dynamics studies.     

Osmolyte‐induced conformational changes in the Hsp90 molecular chaperone

Osmolytes are small molecules that play a central role in cellular homeostasis and the stress response by maintaining protein thermodynamic stability at controlled levels. The underlying physical chemistry that describes how different osmolytes impact folding free energy is well understood, however little is known about their influence on other crucial aspects of protein behavior, such as native‐state conformational changes. Here we investigate this issue with the Hsp90 molecular chaperone, a large dimeric protein that populates a complex conformational equilibrium. Using small angle X‐ray scattering we observe dramatic osmolyte‐dependent structural changes within the native ensemble. The degree to which different osmolytes affect the Hsp90 conformation strongly correlates with thermodynamic metrics of their influence on stability. This observation suggests that the well‐established osmolyte principles that govern stability also apply to large‐scale conformational changes, a proposition that is corroborated by structure‐based fitting of the scattering data, surface area comparisons and m‐value analysis. This approach shows how osmolytes affect a highly cooperative open/closed structural transition between two conformations that differ by a domain‐domain interaction. Hsp90 adopts an additional ligand‐specific conformation in the presence of ATP and we find that osmolytes do not significantly affect this conformational change. Together, these results extend the scope of osmolytes by suggesting that they can maintain protein conformational heterogeneity at controlled levels using similar underlying principles that allow them to maintain protein stability; however the relative impact of osmolytes on different structural states can vary significantly.     

Polyelectrolytic aspects of the thermodynamics of conformational transition: κ‐carrageenan in formamide

A recently published method for the determination of the enthalpy and entropy changes of nonionic origin upon conformational transition of linear biopolyelectrolytes in solution [J. C. Benegas, A. Cesàro, R. Rizzo and S. Paoletti (1998) Biopolymers, Vol. 45, pp. 203–216] has been extended from the case of aqueous salt solutions to that of the organic solvent formamide (FA). The calculation have been applied to the case of the intramolecular transition of the K⁺ salt form of the sulfated polysaccharide κ‐carrageenan. The method proved to be effective in providing the desired data in FA, as it has been previously been successful for the water cases. The comparison between the predicted enthalpy change of transition and the microcalorimetric experimental one turned out to be excellent, thereby ensuring on the validity of the approach. © 1999 John Wiley & Sons, Inc. Biopoly 49: 127–130, 1999     

Redox potentials of the blue copper sites of bilirubin oxidases

The redox potentials of the multicopper redox enzyme bilirubin oxidase (BOD) from two organisms were determined by mediated and direct spectroelectrochemistry. The potential of the T1 site of BOD from the fungus Myrothecium verrucaria was close to 670 mV, whereas that from Trachyderma tsunodae was > 650 mV vs. NHE. For the first time, direct electron transfer was observed between gold electrodes and BODs. The redox potentials of the T2 sites of both BODs were near 390 mV vs. NHE, consistent with previous finding for laccase and suggesting that the redox potentials of the T2 copper sites of most blue multicopper oxidases are similar, about 400 mV.     

Upregulation in the expression of multidrug resistance protein Mrp1 mRNA and protein by increased bilirubin production in rat

Earlier studies suggest that Mrp1 may mediate ATP-dependent cellular extrusion of unconjugated bilirubin (UCB). We studied the serial responses of expression of Mrp1 mRNA and protein in rats with increased bilirubin production due to hemolysis induced by phenylhydrazine (PHZ) treatment. Mrp1 mRNA was analyzed by quantitative PCR and protein by Western blot. Hepatic expression of Mrp1 mRNA and protein peaked at day 3 of PHZ treatment. Splenic expression of Mrp1 mRNA peaked within 24h and returned to baseline at day 5 whereas Mrp1 protein expression peaked at day 3. Pretreatment with heme-oxygenase inhibitor, tin mesoporphyrin, blunted the increase in serum UCB and erased the overexpression of Mrp1 both in liver and spleen. Thus, the upregulation of Mrp1 in hemolysis is mediated by UCB and/or other products of heme oxygenase, further supporting a role of Mrp1 in UCB transport and protection from its cellular toxicity.     

Conformational states of the glucocorticoid receptor DNA-binding domain from molecular dynamics simulations

Molecular dynamics simulations (MD) have been performed on variant crystal and NMR-derived structures of the glucocorticoid receptor DNA-binding domain (GR DBD). A loop region five residues long, the so-called D-box, exhibits significant flexibility, and transient perturbations of the tetrahedral geometry of two structurally important Cys4 zinc finger are seen, coupled to conformational changes in the D-box. In some cases, one of the Cys ligands to zinc exchanges with water, although no global distortion of the protein structure is observed. Thus, from MD simulation, dynamics of the D-box could partly be explained by solvent effects in conjunction with structural reformation of the zinc finger.     

The fluorescence quenching of Ru(bipy)32+: an application for the determination of bilirubin in biological samples

A simple, sensitive and efficient fluorescence method has been established for the quantitative analysis of bilirubin. The fluorometric determination method was based on the kinetic quenching of ruthenium(II) fluorescence. The quenching effect may be due to the complexation reaction of bilirubin with ruthenium(II). Therefore, the effects of ruthenium concentrations and different surfactants have been studied. Under the optimized experimental parameters, the fluorescence intensity decreased proportionally with the bilirubin concentration and linearity was established in the range of 3.3 × 10⁻⁷ to 3.0 × 10⁻⁴ M bilirubin. The detection limit calculated from the calibration graph was found to be 5.2 × 10⁻⁸ M. The relative standard deviation (RSD) of 10 consecutive measurements of 8.0 × 10⁻⁶ M bilirubin was 3.0%, while the recoveries of bilirubin in both human serum and urine samples were obtained in the range 94.0–99.5%. The interference study shows that the developed fluorescence based technique is fast, easy to carry out and shows negligible interference. The developed technique was successfully applied for the analysis of bilirubin in human urine and serum samples. All the experimental results and quality parameters confirmed the sensitivity and reproducibility of the proposed technique for bilirubin determination in human urine and serum samples .     

Insight into the intrinsic flexibility of the SL1 stem-loop from genomic RNA of HIV-1 as probed by molecular dynamics simulation

The SL1 stem-loop is the dimerization initiation site for linking the two copies of the RNA forming the HIV-1 genome. The 26 nucleotides stem contains a defect consisting on a highly conserved G-rich 1-3 asymmetrical internal loop, which is a major site for nucleocapsid protein binding. Several NMR attempts were undertaken to determine the internal loop structure in the SL1 monomer. However, the RNA constructs used in the different studies were largely mutated, in particular with replacement of the nine nucleotides apical loop by a tetraloop, and divergent results were obtained ranging from a rigid structure to a particularly large flexibility. To investigate the reasons for such discrepancies, we used molecular dynamics simulation of the SL1 monomer to probe the effect of mutations on the conformational stability of the internal loop and of the whole stem. It is found that in the wild-type sequence, the internal loop displays conformational variability originating mainly from the nine nucleotides apical loop flexibility that causes large conformational fluctuations (without changing the average structure) in the 7 bp duplex linking the apical and internal loops. The large amplitude atomic motions in the duplex are transmitted to the internal loop in which they induce conformational changes characterized by a labile hydrogen bond network such as G5 successively H-bonded to A29 and G30. On the contrary, with a four nucleotides apical loop, conformational fluctuations in the duplex are reduced by a factor of 2 and are not sufficiently energizing for promoting changes in the internal loop structure at the time scale of the simulations.     

Why is Leu55-->Pro55 transthyretin variant the most amyloidogenic: insights from molecular dynamics simulations of transthyretin monomers

Transthyretin (TTR) is one of the known human amyloidogenic proteins. Its native state is a homotetramer with each monomer having a beta-sandwich structure. Strong experimental evidence suggests that TTR dissociates into monomeric intermediates and that the monomers subsequently self-assemble to form amyloid deposits and insoluble fibrils. However, details on the early steps along the pathway of TTR amyloid formation are unclear, although various experimental approaches with resolutions at the molecular or residue level have provided some clues. It is highly likely that the stability and flexibility of monomeric TTR play crucial roles in the early steps of amyloid formation; thereby, it is essential to characterize initial conformational changes of TTR monomers. In this article we probe the possibility that the differences in the monomeric forms of wild-type (WT) TTR and its variants are responsible for differential amyloidogenesis. We begin with the simulations of WT, Val30-->Met (V30M), and Leu55-->Pro (L55P) TTR monomers. Nanosecond time scale molecular dynamics simulations at 300 K were performed using AMBER. The results indicate that the L55P-TTR monomer undergoes substantial structural changes relative to fluctuations observed in the WT and V30M TTR monomers. The observation supports earlier speculation that the L55P mutation may lead to disruption of the beta-sheet structure through the disorder of the "edge strands" that might facilitate amyloidogenesis.     

Elucidating the interaction of ticlopidine with serum albumin and its role in bilirubin displacement in vitro

Ticlopidine is an anti-platelet drug that functions as a P2Y₁₂ receptor antagonist. The present study provides a detailed characterization of interaction of ticlopidine with a model transport protein, bovine serum albumin (BSA) as well as an assessment of its bilirubin displacing ability using a multi-spectroscopic approach in combination with isothermal titration calorimetry. The value of binding constant determined using ITC studies was found to be 3.03 × 10³ M^?1 with a binding stoichiometry of approximately 1:1. Competitive site marker experiments indicate that ticlopidine binds to Sudlow site I, located in subdomain IIA of BSA. In addition, Circular dichroism and 3D fluorescence spectroscopy indicated structural and conformational changes in BSA on interaction with ticlopidine. Thermodynamic parameters suggested that the reaction was spontaneous, exothermic, entropically driven, and involved hydrophobic interactions. These results were well supported by those obtained through molecular docking studies. Additionally, the effect of ticlopidine on bilirubin and albumin interaction was evaluated using the peroxidase method as well as through fluorescence spectroscopy. Ticlopidine was found to displace bilirubin from serum albumin. Moreover, the binding constant of bilirubin–serum albumin interaction also decreased in presence of ticlopidine. The results indicated that ticlopidine is a competitive displacer of bilirubin in vitro and may contribute to the incidences hyperbilirubinemia associated with the usage of this drug.     

Conformational analysis of endomorphin-2 by molecular dynamics methods

Endomorphin-2 (EM2, H-Tyr-Pro-Phe-Phe-NH(2)) is a highly potent and selective mu-opioid receptor agonist. A conformational analysis of EM2 was carried out by simulated annealing (SA) and molecular dynamics (MD) methods. Molecular modeling was conducted on both neutral (N-terminal NH(2)) and charged (N-terminal NH(3) (+)) molecules. Based on the results of NMR investigations showing an equilibrium mixture of cis and trans Tyr(1)-Pro(2) peptide bonds for EM2 in solution, simulations were performed with restrained cis-Pro and trans-Pro peptide bonds, too. A separate SA study with unrestrained Pro peptide bonds was also conducted. Preferred conformational states are presented in Ramachandran plots. The g(+), g(-), and trans populations of the aromatic amino acid residue side chains were determined in chi(1) space. The distances between the N-terminal N atom and the other backbone N and O atoms, and the distances between the centers of the aromatic rings and the Pro(2) ring, were determined. The energy distribution of the structures obtained by SA was calculated. The preferred secondary structural elements were different kinds of beta-turns, an inverse gamma-turn located in the N-terminal region, and regular and inverse gamma-turns located in the C-terminal region. These turns were stabilized by intramolecular H-bonds and bifurcated H-bonds.     

The role of conformational selection in the molecular recognition of the wild type and mutants XPA67‐80 peptides by ERCC1

Molecular recognition is a fundamental step in the coordination of biomolecular pathways. Understanding how recognition and binding occur between highly flexible protein domains is a complex task. The conformational selection theory provides an elegant rationalization of the recognition mechanism, especially valid in cases when unstructured protein regions are involved. The recognition of a poorly structured peptide, namely XPA_67‐80, by its target receptor ERCC1, falls in this challenging study category. The microsecond molecular dynamics (MD) simulations, discussed in this work, show that the conformational propensity of the wild type XPA_67‐80 peptide in solution supports conformational selection as the key mechanism driving its molecular recognition by ERCC1. Moreover, all the mutations of the XPA_67‐80 peptide studied here cause a significant increase of its conformational disorder, relative to the wild type. Comparison to experimental data suggests that the loss of the recognized structural motifs at the microscopic time scale can contribute to the critical decrease in binding observed for one of the mutants, further substantiating the key role of conformational selection in recognition. Ultimately, because of the high sequence identity and analogy in binding, it is conceivable that the conclusions of this study on the XPA_67‐80 peptide also apply to the ERCC1‐binding domain of the XPA protein. Proteins 2015; 83:1341–1351. © 2015 Wiley Periodicals, Inc.