Thermodynamics of the binding of L-arabinose and of D-galactose to the L-arabinose-binding protein of Escherichia coli

From X-ray crystallographic structure to intrinsic thermodynamics of protein-ligand binding using carbonic anhydrase isozymes as a model system

Carbonic anhydrase (CA) was among the first proteins whose X-ray crystal structure was solved to atomic resolution. CA proteins have essentially the same fold and similar active centers that differ in only several amino acids. Primary sulfonamides are well defined, strong and specific binders of CA. However, minor variations in chemical structure can significantly alter their binding properties. Over 1000 sulfonamides have been designed, synthesized and evaluated to understand the correlations between the structure and thermodynamics of their binding to the human CA isozyme family. Compound binding was determined by several binding assays: fluorescence-based thermal shift assay, stopped-flow enzyme activity inhibition assay, isothermal titration calorimetry and competition assay for enzyme expressed on cancer cell surfaces. All assays have advantages and limitations but are necessary for deeper characterization of these protein-ligand interactions. Here, the concept and importance of intrinsic binding thermodynamics is emphasized and the role of structure-thermodynamics correlations for the novel inhibitors of CA IX is discussed - an isozyme that is overexpressed in solid hypoxic tumors, and thus these inhibitors may serve as anticancer drugs. The abundant structural and thermodynamic data are assembled into the Protein-Ligand Binding Database to understand general protein-ligand recognition principles that could be used in drug discovery.

Investigating Protein Binding with the Isothermal Ligand-induced Resolubilization Assay

Target engagement assays typically detect and quantify the direct physical interaction of a protein of interest and its ligand through stability changes upon ligand binding. Commonly used target engagement methods detect ligand-induced stability by subjecting samples to thermal or proteolytic stress. Here we describe a new variation to these approaches called Isothermal Ligand-induced Resolubilization Assay (ILIRA), which utilizes lyotropic solubility stress to measure ligand binding through changes in target protein solubility. We identified distinct buffer systems and salt concentrations that compromised protein solubility for four diverse proteins: dihydrofolate reductase (DHFR), nucleoside diphosphate-linked moiety X motif 5 (NUDT5), poly [ADP-ribose] polymerase 1 (PARP1), and protein arginine N-methyltransferase 1 (PRMT1). Ligand-induced solubility rescue was demonstrated for these proteins, suggesting that ILIRA can be used as an additional target engagement technique. Differences in ligand-induced protein solubility were assessed by Coomassie blue staining for SDS-PAGE and dot blot, as well as by NanoOrange, Thioflavin T, and Proteostat fluorescence, thus offering flexibility for readout and assay throughput.

Calculation of Protein Folding Thermodynamics Using Molecular Dynamics Simulations

Despite advances in artificial intelligence methods, protein folding remains in many ways an enigma to be solved. Accurate computation of protein folding energetics could help drive fields such as protein and drug design and genetic interpretation. However, the challenge of calculating the state functions governing protein folding from first-principles remains unaddressed. We present here a simple approach that allows us to accurately calculate the energetics of protein folding. It is based on computing the energy of the folded and unfolded states at different temperatures using molecular dynamics simulations. From this, two essential quantities (ΔH and ΔCp) are obtained and used to calculate the conformational stability of the protein (ΔG). With this approach, we have successfully calculated the energetics of two- and three-state proteins, representatives of the major structural classes, as well as small stability differences (ΔΔG) due to changes in solution conditions or variations in an amino acid residue.

Isothermal chemical denaturation assay for monitoring protein stability and inhibitor interactions

Thermal shift assay (TSA) with altered temperature has been the most widely used method for monitoring protein stability for drug research. However, there is a pressing need for isothermal techniques as alternatives. This urgent demand arises from the limitations of TSA, which can sometimes provide misleading ranking of protein stability and fail to accurately reflect protein stability under physiological conditions. Although differential scanning fluorimetry has significantly improved throughput in comparison to differential scanning calorimetry and differential static light scattering throughput, all these methods exhibit moderate sensitivity. In contrast, current isothermal chemical denaturation (ICD) techniques may not offer the same throughput capabilities as TSA, but it provides more precise information about protein stability and interactions. Unfortunately, ICD also suffers from limited sensitivity, typically in micromolar range. We have developed a novel method to overcome these challenges, namely throughput and sensitivity. The novel Förster Resonance Energy Transfer (FRET)-Probe as an external probe is highly applicable to isothermal protein stability monitoring but also to conventional TSA. We have investigated ICD for multiple proteins with focus on KRASG12C with covalent inhibitors and three chemical denaturants performed at nanomolar protein concentration. Data showed corresponding inhibitor-induced stabilization of KRASG12C to those reported by nucleotide exchange assay.

Retracing the evolution of a modern periplasmic binding protein

Investigating the evolution of structural features in modern multidomain proteins helps to understand their immense diversity and functional versatility. The class of periplasmic binding proteins (PBPs) offers an opportunity to interrogate one of the main processes driving diversification: the duplication and fusion of protein sequences to generate new architectures. The symmetry of their two-lobed topology, their mechanism of binding, and the organization of their operon structure led to the hypothesis that PBPs arose through a duplication and fusion event of a single common ancestor. To investigate this claim, we set out to reverse the evolutionary process and recreate the structural equivalent of a single-lobed progenitor using ribose-binding protein (RBP) as our model. We found that this modern PBP can be deconstructed into its lobes, producing two proteins that represent possible progenitor halves. The isolated halves of RBP are well folded and monomeric proteins, albeit with a lower thermostability, and do not retain the original binding function. However, the two entities readily form a heterodimer in vitro and in-cell. The x-ray structure of the heterodimer closely resembles the parental protein. Moreover, the binding function is fully regained upon formation of the heterodimer with a ligand affinity similar to that observed in the modern RBP. This highlights how a duplication event could have given rise to a stable and functional PBP-like fold and provides insights into how more complex functional structures can evolve from simpler molecular components.

Thermodynamic and kinetic analysis of the LAO binding protein and its isolated domains reveal non-additivity in stability, folding and function

Substrate-binding proteins (SBPs) are used by organisms from the three domains of life for transport and signalling. SBPs are composed of two domains that collectively trap ligands with high affinity and selectivity. To explore the role of the domains and the integrity of the hinge region between them in the function and conformation of SBPs, here, we describe the ligand binding, conformational stability and folding kinetics of the Lysine Arginine Ornithine (LAO) binding protein from Salmonella thiphimurium and constructs corresponding to its two independent domains. LAO is a class II SBP formed by a continuous and a discontinuous domain. Contrary to the expected behaviour based on their connectivity, the discontinuous domain shows a stable native-like structure that binds l-arginine with moderate affinity, whereas the continuous domain is barely stable and shows no detectable ligand binding. Regarding folding kinetics, studies of the entire protein revealed the presence of at least two intermediates. While the unfolding and refolding of the continuous domain exhibited only a single intermediate and simpler and faster kinetics than LAO, the folding mechanism of the discontinuous domain was complex and involved multiple intermediates. These findings suggest that in the complete protein the continuous domain nucleates folding and that its presence funnels the folding of the discontinuous domain avoiding nonproductive interactions. The strong dependence of the function, stability and folding pathway of the lobes on their covalent association is most likely the result of the coevolution of both domains as a single unit.

Ion Mobility Mass Spectrometry (IM-MS) for Structural Biology: Insights Gained by Measuring Mass, Charge, and Collision Cross Section

The investigation of macromolecular biomolecules with ion mobility mass spectrometry (IM-MS) techniques has provided substantial insights into the field of structural biology over the past two decades. An IM-MS workflow applied to a given target analyte provides mass, charge, and conformation, and all three of these can be used to discern structural information. While mass and charge are determined in mass spectrometry (MS), it is the addition of ion mobility that enables the separation of isomeric and isobaric ions and the direct elucidation of conformation, which has reaped huge benefits for structural biology. In this review, where we focus on the analysis of proteins and their complexes, we outline the typical features of an IM-MS experiment from the preparation of samples, the creation of ions, and their separation in different mobility and mass spectrometers. We describe the interpretation of ion mobility data in terms of protein conformation and how the data can be compared with data from other sources with the use of computational tools. The benefit of coupling mobility analysis to activation via collisions with gas or surfaces or photons photoactivation is detailed with reference to recent examples. And finally, we focus on insights afforded by IM-MS experiments when applied to the study of conformationally dynamic and intrinsically disordered proteins.

Discovery of first-in-class nanomolar inhibitors of heptosyltransferase I reveals a new aminoglycoside target and potential alternative mechanism of action

A clinically relevant inhibitor for Heptosyltransferase I (HepI) has been sought after for many years because of its critical role in the biosynthesis of lipopolysaccharides on bacterial cell surfaces. While many labs have discovered or designed novel small molecule inhibitors, these compounds lacked the bioavailability and potency necessary for therapeutic use. Extensive characterization of the HepI protein has provided valuable insight into the dynamic motions necessary for catalysis that could be targeted for inhibition. Structural inspection of Kdo2-lipid A suggested aminoglycoside antibiotics as potential inhibitors for HepI. Multiple aminoglycosides have been experimentally validated to be first-in-class nanomolar inhibitors of HepI, with the best inhibitor demonstrating a Ki of 600 ± 90 nM. Detailed kinetic analyses were performed to determine the mechanism of inhibition while circular dichroism spectroscopy, intrinsic tryptophan fluorescence, docking, and molecular dynamics simulations were used to corroborate kinetic experimental findings. While aminoglycosides have long been described as potent antibiotics targeting bacterial ribosomes' protein synthesis leading to disruption of the stability of bacterial cell membranes, more recently researchers have shown that they only modestly impact protein production. Our research suggests an alternative and novel mechanism of action of aminoglycosides in the inhibition of HepI, which directly leads to modification of LPS production in vivo. This finding could change our understanding of how aminoglycoside antibiotics function, with interruption of LPS biosynthesis being an additional and important mechanism of aminoglycoside action. Further research to discern the microbiological impact of aminoglycosides on cells is warranted, as inhibition of the ribosome may not be the sole and primary mechanism of action. The inhibition of HepI by aminoglycosides may dramatically alter strategies to modify the structure of aminoglycosides to improve the efficacy in fighting bacterial infections.

Ligand Binding Introduces Significant Allosteric Shifts in the Locations of Protein Fluctuations

Allostery is usually considered to be a mechanism for transmission of signals associated with physical or dynamic changes in some part of a protein. Here, we investigate the changes in fluctuations across the protein upon ligand binding based on the fluctuations computed with elastic network models. These results suggest that binding reduces the fluctuations at the binding site but increases fluctuations at remote sites, but not to fully compensating extents. If there were complete conservation of entropy, then only the enthalpies of binding would matter and not the entropies; however this does not appear to be the case. Experimental evidence also suggests that energies and entropies of binding can compensate but that the extent of compensation varies widely from case to case. Our results do however always show transmission of an allosteric signal to distant locations where the fluctuations are increased. These fluctuations could be used to compute entropies to improve evaluations of the thermodynamics of binding. We also show the allosteric relationship between peptide binding in the GroEL trans-ring that leads directly to the release of GroES from the GroEL-GroES cis-ring. This finding provides an example of how calculating these changes to protein dynamics induced by the binding of an allosteric ligand can regulate protein function and mechanism.

Engineering Cofactor Specificity of a Thermostable Phosphite Dehydrogenase for a Highly and Robust NADPH Regeneration System

Nicotinamide adenine dinucleotide phosphate (NADP)-dependent dehydrogenases catalyze a range of chemical reactions useful for practical applications. However, their dependence on the costly cofactor, NAD(P)H remains a challenge which must be addressed. Here, we engineered a thermotolerant phosphite dehydrogenase from Ralstonia sp. 4506 (RsPtxD) by relaxing the cofactor specificity for a highly efficient and robust NADPH regeneration system. The five amino acid residues, Cys174-Pro178, located at the C-terminus of β7-strand region in the Rossmann-fold domain of RsPtxD, were changed by site-directed mutagenesis, resulting in four mutants with a significantly increased preference for NADP. The catalytic efficiency of mutant RsPtxD_HARRA for NADP (K _cat/K _M)^NADP was 44.1 μM^-1 min^-1, which was the highest among the previously reported phosphite dehydrogenases. Moreover, the RsPtxD_HARRA mutant exhibited high thermostability at 45°C for up to 6 h and high tolerance to organic solvents, when bound with NADP. We also demonstrated the applicability of RsPtxD_HARRA as an NADPH regeneration system in the coupled reaction of chiral conversion of 3-dehydroshikimate to shikimic acid by the thermophilic shikimate dehydrogenase of Thermus thermophilus HB8 at 45°C, which could not be supported by the parent RsPtxD enzyme. Therefore, the RsPtxD_HARRA mutant might be a promising alternative NADPH regeneration system for practical applications.

Acrylamide-hemoglobin adduct: A spectroscopic study

Acrylamide is a neurotoxic and carcinogenic organic compound that is able to bind to several biomolecules and form adducts, through nucleophilic addition and in vivo by the Maillard Reaction, interfering with the biological functions of these molecules. Hemoglobin is one of the most abundant intracellular blood proteins, and thus it is of high interest to understand whether the binding of acrylamide can alter its properties. The interaction of acrylamide with hemoglobin was assessed in a 20:1 ratio, and after a 72 h-incubation period, a decrease of ca. 50% in the absorbance of the hemoglobin's Soret band was observed at 37 °C. This together with the analysis of circular dichroism spectra indicate that acrylamide binds in close proximity to the heme group. These perturbations were confirmed to not correspond to the loss of the heme group and were mostly reverted after passing the protein through a size-exclusion chromatographic matrix, suggesting a dominant non-covalent interaction for the observed effect. The thermodynamic parameters of unfolding in the absence and presence of acrylamide, suggest an interaction based on H-bonds and van der Waals forces that slightly stabilizes hemoglobin. The oxygen binding capacity of hemoglobin does not seem to be hindered, as no differences in the Q bands were observed in the adduct.

Searching for Antipneumococcal Targets: Choline-Binding Modules as Phagocytosis Enhancers

Choline-binding proteins (CBPs) from Streptococcus pneumoniae comprise a family of modular polypeptides involved in essential events of this pathogen. They recognize the choline residues present in the teichoic and lipoteichoic acids of the cell wall using the so-called choline-binding modules (CBMs). The importance of CBPs in pneumococcal physiology points to them as novel targets to combat antimicrobial resistances shown by this organism. In this work we have tested the ability of exogenously added CBMs to act as CBP inhibitors by competing with the latter for the binding to the choline molecules in the bacterial surface. First, we carried out a thorough physicochemical characterization of three native CBMs, namely C-LytA, C-Cpl1, and C-CbpD, and assessed their affinity for choline and macromolecular, pneumococcal cell-wall mimics. The interaction with these substrates was evaluated by molecular modeling, analytical ultracentrifugation, surface plasmon resonance, and fluorescence and circular dichroism spectroscopies. Van't Hoff thermal analyses unveiled the existence of one noncanonical choline binding site in each of the C-Cpl1 and C-CbpD proteins, leading in total to 5 ligand-binding sites per dimer and 4 sites per monomer, respectively. Remarkably, the binding affinities of the CBMs do not directly correlate with their native oligomeric state or with the number of choline-binding sites, suggesting that choline recognition by these modules is a complex phenomenon. On the other hand, the exogenous addition of CBMs to pneumococcal planktonic cultures caused extensive cell-chaining probably as a consequence of the inhibition of CBP attachment to the cell wall. This was accompanied by bacterial aggregation and sedimentation, causing an enhancement of bacterial phagocytosis by peritoneal macrophages. In addition, the rational design of an oligomeric variant of a native CBM led to a substantial increase in its antibacterial activity by multivalency effects. These results suggest that CBMs might constitute promising nonlytic antimicrobial candidates based on the natural induction of the host defense system.

Structural analysis of β-L-arabinobiose-binding protein in the metabolic pathway of hydroxyproline-rich glycoproteins in Bifidobacterium longum

Bifidobacterium longum is a symbiotic human gut bacterium that has a degradation system for β-arabinooligosaccharides, which are present in the hydroxyproline-rich glycoproteins of edible plants. Whereas microbial degradation systems for α-linked arabinofuranosyl carbohydrates have been extensively studied, little is understood about the degradation systems targeting β-linked arabinofuranosyl carbohydrates. We functionally and structurally analyzed a substrate-binding protein (SBP) of a putative ABC transporter (BLLJ_0208) in the β-arabinooligosaccharide degradation system. Thermal shift assays and isothermal titration calorimetry revealed that the SBP specifically bound Araf-β1,2-Araf (β-Ara₂ ) with a K_d of 0.150 μm, but did not bind L-arabinose or methyl-β-Ara₂ . Therefore, the SBP was termed β-arabinobiose-binding protein (BABP). Crystal structures of BABP complexed with β-Ara₂ were determined at resolutions of up to 1.78 Å. The findings showed that β-Ara₂ was bound to BABP within a short tunnel between two lobes as an α-anomeric form at its reducing end. BABP forms extensive interactions with β-Ara₂ , and its binding mode was unique among SBPs. A molecular dynamics simulation revealed that the closed conformation of substrate-bound BABP is stable, whereas substrate-free form can adopt a fully open and two distinct semi-open states. The importer system specific for β-Ara₂ may contribute to microbial survival in biological niches with limited amounts of digestible carbohydrates. DATABASE: Atomic coordinates and structure factors (codes 6LCE and 6LCF) have been deposited in the Protein Data Bank (http://wwpdb.org/).

Exploring the combination characteristics of lumefantrine, an antimalarial drug and human serum albumin through spectroscopic and molecular docking studies

Binding of lumefantrine (LUM), an antimalarial drug to human serum albumin (HSA), the main carrier protein in human blood circulation was investigated using fluorescence quenching titration, UV-vis absorption and circular dichroism (CD) spectroscopy as well as molecular docking. LUM-induced quenching of the protein (HSA) fluorescence was characterized as static quenching, as revealed by the decrease in the value of the Stern-Volmer quenching constant, K_sv with increasing temperature, thus suggesting LUM-HSA complex formation. This was also confirmed from the UV-vis absorption spectral results. Values of the association constant, K_a for LUM-HSA interaction were found to be within the range, 7.27-5.01 × 10⁴ M^-1 at three different temperatures, i.e. 288 K, 298 K and 308 K, which indicated moderate binding affinity between LUM and HSA. The LUM-HSA complex was stabilized by hydrophobic interactions, H-bonds, as well as van der Waals forces, as predicted from the thermodynamic data (ΔS = +50.34 J mol^-1 K^-1 and ΔH = -12.3 kJ mol^-1) of the binding reaction. Far-UV and near-UV CD spectral results demonstrated smaller changes in both secondary and tertiary structures of HSA upon LUM binding, while three-dimensional fluorescence spectra suggested alterations in the microenvironment around protein fluorophores (Trp and Tyr). LUM binding to HSA offered stability to the protein against thermal stress. Competitive drug displacement results designated Sudlow's Site I, located in subdomain IIA of HSA as the preferred binding site of LUM on HSA, which was well supported by molecular docking analysis.Communicated by Ramaswamy H. Sarma.

Unfolding of CBP21, a lytic polysaccharide monooxygenase, without dissociation of its copper ion cofactor

Chitin-binding protein 21 (CBP21) from Serratia marcescens is a lytic polysaccharide monooxygenase that contains a copper ion as a cofactor. We aimed to elucidate the unfolding mechanism of CBP21 and the effects of Cu²⁺ on its structural stability at pH 5.0. Thermal unfolding of both apo- and holoCBP21 was reversible. ApoCBP21 unfolded in a simple two-state transition manner. The peak temperature of the DSC curve, t_p , for holoCBP21 (74.4°C) was about nine degrees higher than that for apoCBP21 (65.6°C). The value of t_p in the presence of excess Cu²⁺ was around 75°C, indicating that Cu²⁺ does not dissociate from the protein molecule during unfolding. The unfolding mechanism of holoCBP21 was considered to be as follows: N∙Cu²⁺ ⇌ U∙Cu²⁺ , where N and U represent the native and unfolded states, respectively. Urea-induced equilibrium unfolding analysis showed that holoCBP21 was stabilized by 35 kJ mol^-1 in terms of the Gibbs energy change for unfolding (pH 5.0, 25°C), compared with apoCBP21. The increased stability of holoCBP21 was considered to result from the structural stabilization of the protein-Cu²⁺ complex itself.

Technical Capabilities and Limitations of Optical Spectroscopy and Calorimetry Using Water-Miscible Solvents: The Case of Dimethyl Sulfoxide, Acetonitrile, and 1,4-Dioxane

In drug development, water-miscible solvents are commonly used to dissolve drug substances. Typical routine procedures in drug development include dilution of the stock drug solution into an aqueous solution containing target macromolecules for drug binding assays. However, water-miscible solvents impose some technical limitations on the assays on account of their light absorption and heat capacity. Here, we examined the effects of the dilution of 3 water-miscible solvents, that is, dimethyl sulfoxide, acetonitrile, and 1,4-dioxane, on the baseline stability and signal/noise ratio in circular dichroism spectroscopy, isothermal titration calorimetry, and differential scanning calorimetry. Dimethyl sulfoxide and 1,4-dioxane affect the signal/noise ratio of circular dichroism spectra at typically used concentrations due to their light absorbance. The water-miscible solvents generate interfering signals in the isothermal titration calorimetry due to their mixing heat. They show negative or positive slope in the differential scanning calorimetry. Such interfering effects of the solvents are reduced by appropriate dilution according to the analytical techniques. Because the water-miscible solvents have solubilization capacity for alkyl chain moieties and aromatic moieties of chemicals, drug substances containing these moieties can be dissolved into the solvents and then subjected to the analyses to examine their interactions with target proteins after appropriate dilution of the drug solutions.

Small-Molecule Inhibitors of the SOX18 Transcription Factor

Pharmacological modulation of transcription factors (TFs) has only met little success over the past four decades. This is mostly due to standard drug discovery approaches centered on blocking protein/DNA binding or interfering with post-translational modifications. Recent advances in the field of TF biology have revealed a central role of protein-protein interaction in their mode of action. In an attempt to modulate the activity of SOX18 TF, a known regulator of vascular growth in development and disease, we screened a marine extract library for potential small-molecule inhibitors. We identified two compounds, which inspired a series of synthetic SOX18 inhibitors, able to interfere with the SOX18 HMG DNA-binding domain, and to disrupt HMG-dependent protein-protein interaction with RBPJ. These compounds also perturbed SOX18 transcriptional activity in a cell-based reporter gene system. This approach may prove useful in developing a new class of anti-angiogenic compounds based on the inhibition of TF activity.

Ligand binding specificity of the Escherichia coli periplasmic histidine binding protein, HisJ

The HisJ protein from Escherichia coli and related Gram negative bacteria is the periplasmic component of a bacterial ATP-cassette (ABC) transporter system. Together these proteins form a transmembrane complex that can take up L-histidine from the environment and translocate it into the cytosol. We have studied the specificity of HisJ for binding L-His and many related naturally occurring compounds. Our data confirm that L-His is the preferred ligand, but that 1-methyl-L-His and 3-methyl-L-His can also bind, while the dipeptide carnosine binds weakly and D-histidine and the histidine degradation products, histamine, urocanic acid and imidazole do not bind. L-Arg, homo-L-Arg, and post-translationally modified methylated Arg-analogs also bind with reasonable avidity, with the exception of symmetric dimethylated-L-Arg. In contrast, L-Lys and L-Orn have considerably weaker interactions with HisJ and methylated and acetylated Lys variants show relatively poor binding. It was also observed that the carboxylate group of these amino acids and their variants was very important for proper recognition of the ligand. Taken together our results are a key step towards designing HisJ as a specific protein-based reagentless biosensor.

Molecular Determinants of Mutant Phenotypes, Inferred from Saturation Mutagenesis Data

Understanding how mutations affect protein activity and organismal fitness is a major challenge. We used saturation mutagenesis combined with deep sequencing to determine mutational sensitivity scores for 1,664 single-site mutants of the 101 residue Escherichia coli cytotoxin, CcdB at seven different expression levels. Active-site residues could be distinguished from buried ones, based on their differential tolerance to aliphatic and charged amino acid substitutions. At nonactive-site positions, the average mutational tolerance correlated better with depth from the protein surface than with accessibility. Remarkably, similar results were observed for two other small proteins, PDZ domain (PSD95^pdz3) and IgG-binding domain of protein G (GB1). Mutational sensitivity data obtained with CcdB were used to derive a procedure for predicting functional effects of mutations. Results compared favorably with those of two widely used computational predictors. In vitro characterization of 80 single, nonactive-site mutants of CcdB showed that activity in vivo correlates moderately with thermal stability and solubility. The inability to refold reversibly, as well as a decreased folding rate in vitro, is associated with decreased activity in vivo. Upon probing the effect of modulating expression of various proteases and chaperones on mutant phenotypes, most deleterious mutants showed an increased in vivo activity and solubility only upon over-expression of either Trigger factor or SecB ATP-independent chaperones. Collectively, these data suggest that folding kinetics rather than protein stability is the primary determinant of activity in vivo. This study enhances our understanding of how mutations affect phenotype, as well as the ability to predict fitness effects of point mutations.

Thermal stability and binding energetics of thymidylate synthase ThyX

The bacterial thymidylate synthase ThyX is a multisubstrate flavoenzyme that takes part in the de novo synthesis of thymidylate in a variety of microorganisms. Herein we study the effect of FAD and dUMP binding on the thermal stability of wild type (WT) ThyX from the mesophilic Paramecium bursaria chlorella virus-1 (PBCV-1) and from the thermophilic bacterium Thermotoga maritima (TmThyX), and from two variants of TmThyX, Y91F and S88W, using differential scanning calorimetry. The energetics underlying these processes was characterized by isothermal titration calorimetry. The PBCV-1 protein is significantly less stable against the thermal challenge than the TmThyX WT. FAD exerted stabilizing effect greater for PBCV-1 than for TmThyX and for both mutants, whereas binding of dUMP to FAD-loaded proteins stabilized further only TmThyX. Different thermodynamic signatures describe the FAD binding to the WT ThyX proteins. While TmThyX binds FAD with a low μM binding affinity in a process characterized by a favorable entropy change, the assembly of PBCV-1 with FAD is governed by a large enthalpy change opposed by an unfavorable entropy change resulting in a relatively strong nM binding. An enthalpy-driven formation of a high affinity ternary ThyX/FAD/dUMP complex was observed only for TmThyX.

Assessment of the Dissociation Energetics of Some Selected Ligand Drugs Bound on Human Serum Albumin by Differential Scanning Calorimetry

Drug-protein binding may play a role in the thermal energetics of protein denaturation and could lead to the determination of its equilibrium dissociation parameter. The aim of this study was to assess the energetics of a drug that was bound to human serum albumin (HSA) during thermal denaturation. Drugs that were bound at a single high-affinity primary binding site on HSA, including diazepam and ibuprofen, were employed. Commercial HSA was treated with charcoal to remove stabilizers and adjusted to 20% w/v in a pH 7.4 buffered solution. Serial concentrations of individual drugs up to 0.16 mmole/g-protein were added to the cleaned HSA solutions whereas diazepam was added to a commercial HSA solution. Samples were subjected to differential scanning calorimetry (DSC) set to run from 37 to 90°C at 3.0°C/min. Binding of the drug slightly increased the denaturing temperature of the cleaned HSA due to a shift in the equilibrium toward the native protein bound with the drug. Diazepam depressed the denaturing temperature of the commercial HSA by competing with the stabilizers already bound to the primary site of the HSA. This yielded not only the HSA-stabilizer but also the HSA-diazepam type complexes that exhibited a different denaturation process. A rational approximation of the Lumry-Eyring protein denaturation model was used to treat the DSC endotherms. The approximated scheme: [Formula: see text] was successfully fitted to the data. It was used to determine the dissociation parameters for diazepam and ibuprofen bound to the HSA. These results were comparable to those obtained from other methods.

The Protein Acetyltransferase PatZ from Escherichia coli Is Regulated by Autoacetylation-induced Oligomerization

Lysine acetylation is an important post-translational modification in the metabolic regulation of both prokaryotes and eukaryotes. In Escherichia coli, PatZ (formerly YfiQ) is the only known acetyltransferase protein and is responsible for acetyl-CoA synthetase acetylation. In this study, we demonstrated PatZ-positive cooperativity in response to acetyl-CoA and the regulation of acetyl-CoA synthetase activity by the acetylation level. Furthermore, functional analysis of an E809A mutant showed that the conserved glutamate residue is not relevant for the PatZ catalytic mechanism. Biophysical studies demonstrated that PatZ is a stable tetramer in solution and is transformed to its octameric form by autoacetylation. Moreover, this modification is reversed by the sirtuin CobB. Finally, an in silico PatZ tetramerization model based on hydrophobic and electrostatic interactions is proposed and validated by three-dimensional hydrodynamic analysis. These data reveal, for the first time, the structural regulation of an acetyltransferase by autoacetylation in a prokaryotic organism.

Thermal, chemical and pH induced unfolding of turmeric root lectin: modes of denaturation

Curcuma longa rhizome lectin, of non-seed origin having antifungal, antibacterial and α-glucosidase inhibitory activities, forms a homodimer with high thermal stability as well as acid tolerance. Size exclusion chromatography and dynamic light scattering show it to be a dimer at pH 7, but it converts to a monomer near pH 2. Circular dichroism spectra and fluorescence emission maxima are virtually indistinguishable from pH 7 to 2, indicating secondary and tertiary structures remain the same in dimer and monomer within experimental error. The tryptophan environment as probed by acrylamide quenching data yielded very similar data at pH 2 and pH 7, implying very similar folding for monomer and dimer. Differential scanning calorimetry shows a transition at 350.3 K for dimer and at 327.0 K for monomer. Thermal unfolding and chemical unfolding induced by guanidinium chloride for dimer are both reversible and can be described by two-state models. The temperatures and the denaturant concentrations at which one-half of the protein molecules are unfolded, are protein concentration-dependent for dimer but protein concentration-independent for monomer. The free energy of unfolding at 298 K was found to be 5.23 Kcal mol⁻¹ and 14.90 Kcal mol⁻¹ for the monomer and dimer respectively. The value of change in excess heat capacity upon protein denaturation (ΔC_p) is 3.42 Kcal mol⁻¹ K⁻¹ for dimer. The small ΔC_p for unfolding of CLA reflects a buried hydrophobic core in the folded dimeric protein. These unfolding experiments, temperature dependent circular dichroism and dynamic light scattering for the dimer at pH 7 indicate its higher stability than for the monomer at pH 2. This difference in stability of dimeric and monomeric forms highlights the contribution of inter-subunit interactions in the former.

Force spectroscopy studies on protein-ligand interactions: a single protein mechanics perspective

Protein-ligand interactions are ubiquitous and play important roles in almost every biological process. The direct elucidation of the thermodynamic, structural and functional consequences of protein-ligand interactions is thus of critical importance to decipher the mechanism underlying these biological processes. A toolbox containing a variety of powerful techniques has been developed to quantitatively study protein-ligand interactions in vitro as well as in living systems. The development of atomic force microscopy-based single molecule force spectroscopy techniques has expanded this toolbox and made it possible to directly probe the mechanical consequence of ligand binding on proteins. Many recent experiments have revealed how ligand binding affects the mechanical stability and mechanical unfolding dynamics of proteins, and provided mechanistic understanding on these effects. The enhancement effect of mechanical stability by ligand binding has been used to help tune the mechanical stability of proteins in a rational manner and develop novel functional binding assays for protein-ligand interactions. Single molecule force spectroscopy studies have started to shed new lights on the structural and functional consequence of ligand binding on proteins that bear force under their biological settings.

XIAP antagonist embelin inhibited proliferation of cholangiocarcinoma cells

Cholangiocarcinoma cells are dependent on antiapoptotic signaling for survival and resistance to death stimuli. Recent mechanistic studies have revealed that increased cellular expression of the E3 ubiquitin-protein ligase X-linked inhibitor of apoptosis (XIAP) impairs TRAIL- and chemotherapy-induced cytotoxicity, promoting survival of cholangiocarcinoma cells. This study was undertaken to determine if pharmacologic antagonism of XIAP protein was sufficient to sensitize cholangiocarcinoma cells to cell death. We employed malignant cholangiocarcinoma cell lines and used embelin to antagonize XIAP protein. Embelin treatment resulted in decreased XIAP protein levels by 8 hours of treatment with maximal effect at 16 hours in KMCH and Mz-ChA-1 cells. Assessment of nuclear morphology demonstrated a concentration-dependent increase in nuclear staining. Interestingly, embelin induced nuclear morphology changes as a single agent, independent of the addition of TNF-related apoptosis inducing ligand (TRAIL). However, caspase activity assays revealed that increasing embelin concentrations resulted in slight inhibition of caspase activity, not activation. In addition, the use of a pan-caspase inhibitor did not prevent nuclear morphology changes. Finally, embelin treatment of cholangiocarcinoma cells did not induce DNA fragmentation or PARP cleavage. Apoptosis does not appear to contribute to the effects of embelin on cholangiocarcinoma cells. Instead, embelin caused inhibition of cell proliferation and cell cycle analysis indicated that embelin increased the number of cells in S and G2/M phase. Our results demonstrate that embelin decreased proliferation in cholangiocarcinoma cell lines. Embelin treatment resulted in decreased XIAP protein expression, but did not induce or enhance apoptosis. Thus, in cholangiocarcinoma cells the mechanism of action of embelin may not be dependent on apoptosis.

Non-allosteric enzyme switches possess larger effector-induced changes in thermodynamic stability than their non-switch analogs

The ability to regulate cellular protein activity offers a broad range of biotechnological and biomedical applications. Such protein regulation can be achieved by modulating the specific protein activity or through processes that regulate the amount of protein in the cell. We have previously demonstrated that the nonhomologous recombination of the genes encoding maltose binding protein (MBP) and TEM1 β-lactamase (BLA) can result in genes that confer maltose-dependent resistance to β-lactam antibiotics even though the encoded proteins are not allosteric enzymes. We showed that these phenotypic switches-named based on their conferral of a switching phenotype to cells-resulted from a specific interaction with maltose in the cell that increased the switches cellular accumulation. Since phenotypic switches represent an important class of engineered proteins for basic science and biotechnological applications in vivo, we sought to elucidate the phenomena behind the increased accumulation and switching properties. Here, we demonstrate the key role for the linker region between the two proteins. Experimental evidence supports the hypothesis that in the absence of their effector, some phenotypic switches possess an increased rate of unfolding, decreased conformational stability, and increased protease susceptibility. These factors alone or in combination serve to decrease cellular accumulation. The effector functions to increase cellular accumulation by alleviating one or more of these defects. This perspective on the mechanism for phenotypic switching will aid the development of design rules for switch construction for applications and inform the study of the regulatory mechanisms of natural cellular proteins.

Systematic interaction analysis of human lipocalin-type prostaglandin D synthase with small lipophilic ligands

L-PGDS [lipocalin-type PG (prostaglandin) D synthase] is a multi-functional protein, acting as a PGD2-producing enzyme and a lipid-transporter. In the present study, we focus on the function of L-PGDS as an extracellular transporter for small lipophilic molecules. We characterize the binding mechanism of human L-PGDS for the molecules, especially binding affinity stoichiometry and driving force, using tryptophan fluorescence quenching, ICD (induced circular dichroism) and ITC (isothermal titration calorimetry). The tryptophan fluorescence quenching measurements revealed that haem metabolites such as haemin, biliverdin and bilirubin bind to L-PGDS with significantly higher affinities than the other small lipophilic ligands examined, showing dissociation constant (Kd) values from 17.0 to 20.9 nM. We focused particularly on the extra-specificities of haem metabolites and L-PGDS. The ITC and ICD data revealed that two molecules of the haem metabolites bind to L-PGDS with high and low affinities, showing Kd values from 2.8 to 18.1 nM and from 0.209 to 1.63 μM respectively. The thermodynamic parameters for the interactions revealed that the contributions of enthalpy and entropy change were considerably different for each haem metabolite even when the Gibbs energy change was the same. Thus we believe that the binding energy of haem metabolites to L-PGDS is optimized by balancing enthalpy and entropy change.

Linking the structure and thermal stability of beta-galactoside-binding protein galectin-1 to ligand binding and dimerization equilibria

The stability of proteins involves a critical balance of interactions of different orders of magnitude. In this work, we present experimental evidence of an increased thermal stability of galectin-1, a multifunctional beta-galactoside-binding protein, upon binding to the disaccharide lactose. Analysis of structural changes occurring upon binding of lectin to its specific glycans and thermal denaturation of the protein and the complex were analyzed by circular dichroism. On the other hand, we studied dimerization as another factor that may induce structural and thermal stability changes. The results were then complemented with molecular dynamics simulations followed by a detailed computation of thermodynamic properties, including the internal energy, solvation free energy, and conformational entropy. In addition, an energetic profile of the binding and dimerization processes is also presented. Whereas binding and cross-linking of lactose do not alter galectin-1 structure, this interaction leads to substantial changes in the flexibility and internal energy of the protein which confers increased thermal stability to this endogenous lectin. Given that an improved understanding of the physicochemical properties of galectin-glycan lattices may contribute to the dissection of their biological functions and prediction of their therapeutic applications, our study suggests that galectin binding to specific disaccharide ligands may increase the thermal stability of this glycan-binding protein, an effect that could influence its critical biological functions.

Analysis of the Escherichia coli glucosamine-6-phosphate synthase activity by isothermal titration calorimetry and differential scanning calorimetry

Glucosamine-6-phosphate synthase (GlmS) is responsible for the first and rate-limiting step in the hexosamine biosynthetic pathway. It catalyzes the conversion of D-fructose-6P (F6P) into D-glucosamine-6P (GlcN6P) using L-glutamine (Gln) as nitrogen donor (synthase activity) according to an ordered bi-bi process where F6P binds first. In the absence of F6P, the enzyme exhibits a weak hydrolyzing activity of Gln into Glu and ammonia (glutaminase activity), whereas the presence of F6P strongly stimulates it (hemi-synthase activity). Until now, these different activities were indirectly measured using either coupled enzyme or colorimetric methods. In this work, we have developed a direct assay monitoring the heat released by the reaction. Isothermal titration calorimetry and differential scanning calorimetry were used to determine kinetic and thermodynamic parameters of GlmS. The direct determination at 37 degrees C of kinetic parameters and affinity constants for both F6P and Gln demonstrated that part of the ammonia produced by Gln hydrolysis in the presence of both substrates is not used for the formation of the GlcN6P. The full characterization of this phenomenon allowed to identify experimental conditions where this leak of ammonia is negligible. Enthalpy measurements at 25 degrees C in buffers of various heats of protonation demonstrated that no proton exchange with the medium occurred during the enzyme-catalyzed glutaminase or synthase reaction suggesting for the first time that both products are released as a globally neutral pair composed by the Glu carboxylic side chain and the GlcN6P amine function. Finally we showed that the oligomerization state of GlmS is concentration-dependent.

Influence of the H-site residue 108 on human glutathione transferase P1-1 ligand binding: structure-thermodynamic relationships and thermal stability

The effect of the Y108V mutation of human glutathione S-transferase P1-1 (hGST P1-1) on the binding of the diuretic drug ethacrynic acid (EA) and its glutathione conjugate (EASG) was investigated by calorimetric, spectrofluorimetric, and crystallographic studies. The mutation Tyr 108 --> Val resulted in a 3D-structure very similar to the wild type (wt) enzyme, where both the hydrophobic ligand binding site (H-site) and glutathione binding site (G-site) are unchanged except for the mutation itself. However, due to a slight increase in the hydrophobicity of the H-site, as a consequence of the mutation, an increase in the entropy was observed. The Y108V mutation does not affect the affinity of EASG for the enzyme, which has a higher affinity (K(d) approximately 0.5 microM) when compared with those of the parent compounds, K(d) (EA) approximately 13 microM, K(d) (GSH) approximately 25 microM. The EA moiety of the conjugate binds in the H-site of Y108V mutant in a fashion completely different to those observed in the crystal structures of the EA or EASG wt complex structures. We further demonstrate that the Delta C(p) values of binding can also be correlated with the potential stacking interactions between ligand and residues located in the binding sites as predicted from crystal structures. Moreover, the mutation does not significantly affect the global stability of the enzyme. Our results demonstrate that calorimetric measurements maybe useful in determining the preference of binding (the binding mode) for a drug to a specific site of the enzyme, even in the absence of structural information.

Ligand binding strategies of human serum albumin: how can the cargo be utilized?

Human serum albumin (HSA), being the most abundant carrier protein in blood and a modern day clinical tool for drug delivery, attracts high attention among biologists. Hence, its unfolding/refolding strategies and exogenous/endogenous ligand binding preference are of immense use in therapeutics and clinical biochemistry. Among its fellow proteins albumin is known to carry almost every small molecule. Thus, it is a potential contender for being a molecular cargo/or nanovehicle for clinical, biophysical and industrial purposes. Nonetheless, its structure and function are largely regulated by various chemical and physical factors to accommodate HSA to its functional purpose. This multifunctional protein also possesses enzymatic properties which may be used to convert prodrugs to active therapeutics. This review aims to highlight current overview on the binding strategies of protein to various ligands that may be expected to lead to significant clinical applications.

Effects of pH and iminosugar pharmacological chaperones on lysosomal glycosidase structure and stability

Human lysosomal enzymes acid-beta-glucosidase (GCase) and acid-alpha-galactosidase (alpha-Gal A) hydrolyze the sphingolipids glucosyl- and globotriaosylceramide, respectively, and mutations in these enzymes lead to the lipid metabolism disorders Gaucher and Fabry disease, respectively. We have investigated the structure and stability of GCase and alpha-Gal A in a neutral-pH environment reflective of the endoplasmic reticulum and an acidic-pH environment reflective of the lysosome. These details are important for the development of pharmacological chaperone therapy for Gaucher and Fabry disease, in which small molecules bind mutant enzymes in the ER to enable the mutant enzyme to meet quality control requirements for lysosomal trafficking. We report crystal structures of apo GCase at pH 4.5, at pH 5.5, and in complex with the pharmacological chaperone isofagomine (IFG) at pH 7.5. We also present thermostability analysis of GCase at pH 7.4 and 5.2 using differential scanning calorimetry. We compare our results with analogous experiments using alpha-Gal A and the chaperone 1-deoxygalactonijirimycin (DGJ), including the first structure of alpha-Gal A with DGJ. Both GCase and alpha-Gal A are more stable at lysosomal pH with and without their respective iminosugars bound, and notably, the stability of the GCase-IFG complex is pH sensitive. We show that the conformations of the active site loops in GCase are sensitive to ligand binding but not pH, whereas analogous galactose- or DGJ-dependent conformational changes in alpha-Gal A are not seen. Thermodynamic parameters obtained from alpha-Gal A unfolding indicate two-state, van't Hoff unfolding in the absence of the iminosugar at neutral and lysosomal pH, and non-two-state unfolding in the presence of DGJ. Taken together, these results provide insight into how GCase and alpha-Gal A are thermodynamically stabilized by iminosugars and suggest strategies for the development of new pharmacological chaperones for lysosomal storage disorders.

Thermodynamics of single strand DNA base stacking

The thermodynamics of the stacking to unstacking transitions of 24 single-stranded DNA sequences (ssDNA), 10-12 bases in length, in sodium phosphate buffer were determined from 10 to 95 degrees C, using differential scanning calorimetry (DSC). An additional 22 ssDNA sequences did not exhibit an S<=>U transition in this temperature range. The transition properties of the ssDNA sequences with <or=60% self-complementarily in the reverse direction were independent of concentration with transition temperatures ranging from 15 to 70 degrees C, van't Hoff transition enthalpies from 92 to 201 kJ mol(-1) and transition enthalpies from 5 to 75% of the corresponding van't Hoff transition enthalpies. Since all the 16 doublets and 60 of the 64 triplets are present in both the transition and the non-transition ssDNA sequences, it is unlikely that the nucleation subset initiating stacking of the sequence is a specific doublet or triplet subset. Of the 141 quadruplet subsets of the 46 sequences, each transition ssDNA sequence contained at least one or more quadruplets not found in the non-transition ssDNA sequences. It could be concluded that the thermal stability of the stacked conformation was dependent on the presence of a possible nucleation quadruplet and the length of the ssDNA sequence and not on the G or C content of the ssDNA sequence, nor on the number of purine bases in the sequence.

Stability and structural recovery of the tetramerization domain of p53-R337H mutant induced by a designed templating ligand

Protein p53 is a transcription factor crucial for cell cycle and genome integrity. It is able to induce both cell arrest when DNA is damaged and the expression of DNA repair machinery. When the damage is irreversible, it triggers apoptosis. Indeed, the protein, which is a homotetramer, is mutated in most human cancers. For instance, the inherited mutation p53-R337H results in destabilization of the tetramer and, consequently, leads to an organism prone to tumor setup. We describe herein a rational designed molecule capable of holding together the four monomers of the mutated p53-R337H protein, recovering the tetramer integrity as in the wild-type structure. Two ligand molecules, based on a conical calix[4]arene with four cationic guanidiniomethyl groups at the wider edge (upper rim) and hydrophobic loops at the narrower edge (lower rim), fit nicely and cooperatively into the hydrophobic clefts between two of the monomers at each side of the protein and keep the tetrameric structure, like molecular templates, by both ion-pair and hydrophobic interactions. We found a good agreement between the structure of the complex and the nature of the interactions involved by a combination of theory (molecular dynamics) and experiments (circular dichroism, differential scanning calorimetry and (1)H saturation transfer difference NMR).

The relative binding affinities of PDZ partners for CFTR: a biochemical basis for efficient endocytic recycling

The cystic fibrosis transmembrane conductance regulator (CFTR) is an epithelial chloride channel mutated in patients with cystic fibrosis. Its expression and functional interactions in the apical membrane are regulated by several PDZ (PSD-95, discs large, zonula occludens-1) proteins, which mediate protein-protein interactions, typically by binding C-terminal recognition motifs. In particular, the CFTR-associated ligand (CAL) limits cell-surface levels of the most common disease-associated mutant DeltaF508-CFTR. CAL also mediates degradation of wild-type CFTR, targeting it to lysosomes following endocytosis. Nevertheless, wild-type CFTR survives numerous cycles of uptake and recycling. In doing so, how does it repeatedly avoid CAL-mediated degradation? One mechanism may involve competition between CAL and other PDZ proteins including Na (+)/H (+) exchanger-3 regulatory factors 1 and 2 (NHERF1 and NHERF2), which functionally stabilize cell-surface CFTR. Thus, to understand the biochemical basis of WT-CFTR persistence, we need to know the relative affinities of these partners. However, no quantitative binding data are available for CAL or the individual NHERF2 PDZ domains, and published estimates for the NHERF1 PDZ domains conflict. Here we demonstrate that the affinity of the CAL PDZ domain for the CFTR C-terminus is much weaker than those of NHERF1 and NHERF2 domains, enabling wild-type CFTR to avoid premature entrapment in the lysosomal pathway. At the same time, CAL's affinity is evidently sufficient to capture and degrade more rapidly cycling mutants, such as DeltaF508-CFTR. The relatively weak affinity of the CAL:CFTR interaction may provide a pharmacological window for stabilizing rescued DeltaF508-CFTR in patients with cystic fibrosis.

A Gibbs free energy correlation for automated docking of carbohydrates

Thermodynamic information can be inferred from static atomic configurations. To model the thermodynamics of carbohydrate binding to proteins accurately, a large binding data set has been assembled from the literature. The data set contains information from 262 unique protein-carbohydrate crystal structures for which experimental binding information is known. Hydrogen atoms were added to the structures and training conformations were generated with the automated docking program AutoDock 3.06, resulting in a training set of 225,920 all-atom conformations. In all, 288 formulations of the AutoDock 3.0 free energy model were trained against the data set, testing each of four alternate methods of computing the van der Waals, solvation, and hydrogen-bonding energetic components. The van der Waals parameters from AutoDock 1 produced the lowest errors, and an entropic model derived from statistical mechanics produced the only models with five physically and statistically significant coefficients. Eight models predict the Gibbs free energy of binding with an error of less than 40% of the error of any similar models previously published.

Calorimetric and Spectroscopic Studies on the Interaction of Methimazole with Bovine Serum Albumin

The potential binding interaction(s) of the anti-thyroid drug methimazole (MMZ) with the protein bovine serum albumin (BSA) has been studied using isothermal titration calorimetry (ITC) and UV-Visible, fluorescence and circular dichroism (CD) spectroscopic techniques. The binding of MMZ to BSA has been studied in both the presence and absence of added surfactants. Since, the ITC thermograms show the molar enthalpy of binding of MMZ and BSA to be zero within experimental error, either the enthalpy change of the binding interaction is zero or there is no binding occurring. The CD and the intrinsic fluorescence and life time spectra of BSA were unchanged by the addition of MMZ. This is also indicative of the absence of any significant interaction of MMZ with BSA.