A new ultrasensitive scanning calorimeter

Effects of Osmolytes on Ligand Binding to Dihydropteroate Synthase from Bacillus anthracis

Osmolyte interactions with ligands can affect their affinity for proteins and are dependent upon the cosolute and the functional groups of the ligand. Here, we explored ligand binding to Bacillus anthracis dihydropteroate synthase (BaDHPS) under osmotic stress conditions. Osmolyte effects were specific to the cosolute and ligand, suggesting interaction of the osmolytes with the free ligands in solution. The association rates of pterin pyrophosphate were mostly unaffected by the osmolytes, except for a 2-fold decrease in the presence of 1 M trehalose, while the dissociation rates decreased in most osmolyte solutions. The viscosity and dielectric constant of the solution did not correlate with the effects of the osmolytes. Experimental results were compared with predicted preferential interaction coefficients (Δμ₂₃/RT) between the osmolytes and ligands. The Δμ₂₃/RT were able to predict the experimental data for most of the osmolytes. Trehalose and proline effects did not correlate with the predicted values, indicating that these two osmolytes may affect binding in more complex ways than simple preferential interactions. Additionally, osmolytes weakly interacted with the sulfa drug sulfathiazole, which altered its affinity for BaDHPS, suggesting that these types of weak interactions can also impact drug binding. As osmolytes affect ligands binding to two different folate cycle enzymes (DHFRs and DHPS), we predicted how ligand binding to other folate cycle enzymes will be altered by the presence of osmolytes.

Drosophila methionine sulfoxide reductase A (MSRA) lacks methionine oxidase activity

Mouse, human, and E. coli methionine sulfoxide reductase A (MSRA) stereospecifically catalyze both the reduction of S-methionine sulfoxide to methionine and the oxidation of methionine to S-methionine sulfoxide. Calmodulin has 9 methionine residues, but only Met77 is oxidized by MSRA, and this is completely reversed when MSRA operates in the reductase direction. Given the powerful genetic tools available for Drosophila, we selected this model organism to identify the in vivo calmodulin targets regulated by redox modulation of Met77. The active site sequences of mammalian and Drosophila MSRA are identical, and both contain two cysteine residues in their carboxy terminal domains. We produced recombinant Drosophila MSRA and studied its biochemical and biophysical properties. The enzyme is active as a methionine sulfoxide reductase, but it cannot function as a methionine oxidase. The first step in the mammalian oxidase reaction is formation of a sulfenic acid at the active site, and the second step is the reaction of the sulfenic acid with a carboxy terminal domain cysteine to form a disulfide bond. The third step regenerates the active site through a disulfide exchange reaction with a second carboxy terminal domain cysteine. Drosophila MSRA carries out the first and second steps, but it cannot regenerate the active site in the third step. Thus, unlike the E. coli and mammalian enzymes, Drosophila MSRA catalyzes only the reduction of methionine sulfoxide and not the oxidation of methionine.

Differential Scanning Calorimetry - A Method for Assessing the Thermal Stability and Conformation of Protein Antigen

Differential scanning calorimetry (DSC) is an analytical technique that measures the molar heat capacity of samples as a function of temperature. In the case of protein samples, DSC profiles provide information about thermal stability, and to some extent serves as a structural "fingerprint" that can be used to assess structural conformation. It is performed using a differential scanning calorimeter that measures the thermal transition temperature (melting temperature; Tm) and the energy required to disrupt the interactions stabilizing the tertiary structure (enthalpy; ∆H) of proteins. Comparisons are made between formulations as well as production lots, and differences in derived values indicate differences in thermal stability and structural conformation. Data illustrating the use of DSC in an industrial setting for stability studies as well as monitoring key manufacturing steps are provided as proof of the effectiveness of this protocol. In comparison to other methods for assessing the thermal stability of protein conformations, DSC is cost-effective, requires few sample preparation steps, and also provides a complete thermodynamic profile of the protein unfolding process.

Exploring the Balance between DNA Pressure and Capsid Stability in Herpesviruses and Phages

We have recently shown in both herpesviruses and phages that packaged viral DNA creates a pressure of tens of atmospheres pushing against the interior capsid wall. For the first time, using differential scanning microcalorimetry, we directly measured the energy powering the release of pressurized DNA from the capsid. Furthermore, using a new calorimetric assay to accurately determine the temperature inducing DNA release, we found a direct influence of internal DNA pressure on the stability of the viral particle. We show that the balance of forces between the DNA pressure and capsid strength, required for DNA retention between rounds of infection, is conserved between evolutionarily diverse bacterial viruses (phages λ and P22), as well as a eukaryotic virus, human herpes simplex 1 (HSV-1). Our data also suggest that the portal vertex in these viruses is the weakest point in the overall capsid structure and presents the Achilles heel of the virus's stability. Comparison between these viral systems shows that viruses with higher DNA packing density (resulting in higher capsid pressure) have inherently stronger capsid structures, preventing spontaneous genome release prior to infection. This force balance is of key importance for viral survival and replication. Investigating the ways to disrupt this balance can lead to development of new mutation-resistant antivirals.

IMPORTANCE

A virus can generally be described as a nucleic acid genome contained within a protective protein shell, called the capsid. For many double-stranded DNA viruses, confinement of the large DNA molecule within the small protein capsid results in an energetically stressed DNA state exerting tens of atmospheres of pressures on the inner capsid wall. We show that stability of viral particles (which directly relates to infectivity) is strongly influenced by the state of the packaged genome. Using scanning calorimetry on a bacterial virus (phage λ) as an experimental model system, we investigated the thermodynamics of genome release associated with destabilizing the viral particle. Furthermore, we compare the influence of tight genome confinement on the relative stability for diverse bacterial and eukaryotic viruses. These comparisons reveal an evolutionarily conserved force balance between the capsid stability and the density of the packaged genome.

Dendritic domains with hexagonal symmetry formed by x-shaped bolapolyphiles in lipid membranes

A novel class of bolapolyphile (BP) molecules are shown to integrate into phospholipid bilayers and self-assemble into unique sixfold symmetric domains of snowflake-like dendritic shapes. The BPs comprise three philicities: a lipophilic, rigid, π–π stacking core; two flexible lipophilic side chains; and two hydrophilic, hydrogen-bonding head groups. Confocal microscopy, differential scanning calorimetry, XRD, and solid-state NMR spectroscopy confirm BP-rich domains with transmembrane-oriented BPs and three to four lipid molecules per BP. Both species remain well organized even above the main 1,2-dipalmitoyl-sn-glycero-3-phosphocholine transition. The BP molecules only dissolve in the fluid membrane above 70 °C. Structural variations of the BP demonstrate that head-group hydrogen bonding is a prerequisite for domain formation. Independent of the head group, the BPs reduce membrane corrugation. In conclusion, the BPs form nanofilaments by π stacking of aromatic cores, which reduce membrane corrugation and possibly fuse into a hexagonal network in the dendritic domains.

Blind prediction of charged ligand binding affinities in a model binding site

Predicting absolute protein-ligand binding affinities remains a frontier challenge in ligand discovery and design. This becomes more difficult when ionic interactions are involved because of the large opposing solvation and electrostatic attraction energies. In a blind test, we examined whether alchemical free-energy calculations could predict binding affinities of 14 charged and 5 neutral compounds previously untested as ligands for a cavity binding site in cytochrome c peroxidase. In this simplified site, polar and cationic ligands compete with solvent to interact with a buried aspartate. Predictions were tested by calorimetry, spectroscopy, and crystallography. Of the 15 compounds predicted to bind, 13 were experimentally confirmed, while 4 compounds were false negative predictions. Predictions had a root-mean-square error of 1.95 kcal/mol to the experimental affinities, and predicted poses had an average RMSD of 1.7Å to the crystallographic poses. This test serves as a benchmark for these thermodynamically rigorous calculations at predicting binding affinities for charged compounds and gives insights into the existing sources of error, which are primarily electrostatic interactions inside proteins. Our experiments also provide a useful set of ionic binding affinities in a simplified system for testing new affinity prediction methods.

Chemical and mechanical impact of silica nanoparticles on the phase transition behavior of phospholipid membranes in theory and experiment

The interaction of nanoparticles (NPs) with lipid membranes is an integral step in the interaction of NPs and living cells. During particle uptake, the membrane has to bend. Due to the nature of their phase diagram, the modulus of compression of these membranes can vary by more than one order of magnitude, and thus both the thermodynamic and mechanical aspects of the membrane have to be considered simultaneously. We demonstrate that silica NPs have at least two independent effects on the phase transition of phospholipid membranes: 1), a chemical effect resulting from the finite instability of the NPs in water; and 2), a mechanical effect that originates from a bending of the lipid membrane around the NPs. Here, we report on recent experiments that allowed us to clearly distinguish both effects, and present a thermodynamic model that includes the elastic energy of the membranes and correctly predicts our findings both quantitatively and qualitatively.

Thermal denaturation assays in chemical biology

Thermal denaturation-based methods are becoming increasingly used to characterize protein stability and interactions. Recent technical advances have made these methods more suitable for high throughput screening. Reasonable throughput and the ability to perform these screens using commonly used instruments, such as RT-PCR machines or simple plate readers equipped with heating devices, facilitate these experiments in almost any laboratory. Introducing an aggregation-based monitoring approach as well as alternative fluorophores has allowed the screening of a wider range of proteins, including membrane proteins, against large chemical libraries. Thermal denaturation-based methods are independent of protein function, which is especially useful for the identification of orphan protein function. Here, we review applications of thermal denaturation-based methods in characterizing protein stability and ligand binding, and also provide information on protocol modifications that may further increase throughput.

The discovery of compounds that stimulate the activity of kallikrein-related peptidase 3 (KLK3)

Kallikrein-related peptidase 3 (KLK3), also known as prostate-specific antigen (PSA), is the most useful biomarker for prostate cancer (PCa). KLK3 is suggested to play a role in regulating cancer growth through anti-angiogenic activity in vivo and in vitro. This feature, together with its specificity for prostate tissue, makes KLK3 an intriguing target for the design of new therapies for PCa. 3D pharmacophores for KLK3-stimulating compounds were generated based on peptides that bind specifically to KLK3 and increase its enzymatic activity. As a result of pharmacophore-based virtual screening, four small, drug-like compounds with affinity for KLK3 were discovered and validated by capillary differential scanning calorimetry. One of the compounds also stimulated the activity of KLK3, and is therefore the first published small molecule with such an activity.

The second coordination sphere of FIH controls hydroxylation

The factor inhibiting HIF (FIH) is a proximate oxygen sensor for human cells, hydroxylating Asn(803) within the α-subunit of the hypoxia inducible factor (HIF). FIH is an α-ketoglutatrate (αKG)-dependent, non-heme Fe(II) dioxygenase, in which Fe(II) is coordinated by a (His(2)Asp) facial triad, αKG, and H(2)O. Hydrogen bonding among the facial triad, the HIF-Asn(803) side chain, and various second-sphere residues suggests a functional role for the second coordination sphere in tuning the chemistry of the Fe(II) center. Point mutants of FIH were prepared to test the functional role of the αKG-centered (Asn(205) and Asn(294)) or HIF-Asn(803)-centered (Arg(238) and Gln(239)) second-sphere residues. The second sphere was tested for local effects on priming Fe(II) to react with O(2), oxidative decarboxylation, and substrate positioning. Steady-sate kinetics were used to test for overall catalytic effects; autohydroxylation rates were used to test for priming and positioning, and electronic spectroscopy was used to assess the primary coordination sphere and the electrophilicity of αKG. Asn(205) → Ala and Asn(294) → Ala mutants exhibited diminished rates of steady-state turnover, while minimally affecting autohydroxylation, consistent with impaired oxidative decarboxylation. Blue-shifted metal to ligand charge transfer transitions for (Fe+αKG)FIH indicated that these point mutations destabilized the π* orbitals of αKG, further supporting a slowed rate of oxidative decarboxylation. The Arg(238) → Met mutant exhibited steady-state rates too low to measure and diminished product yields, suggesting impaired substrate positioning or priming; the Arg(238) → Met mutant was capable of O(2) activation for the autohydroxylation reaction. The Gln(239) → Asn mutant exhibited significantly slowed steady-state kinetics and diminished product yields, suggesting impaired substrate positioning or priming. As HIF binding to the Gln(239) → Asn mutant stimulated autohydroxylation, it is more likely that this point mutant simply mispositions the HIF-Asn(803) side chain. This work combines kinetics and spectroscopy to show that these second-sphere hydrogen bonds play roles in promoting oxidative decarboxylation, priming Fe(II) to bind O(2), and positioning HIF-Asn(803).

The thermodynamics of simple biomembrane mimetic systems

Insight into the forces governing a system is essential for understanding its behavior and function. Thermodynamic investigations provide a wealth of information that is not, or is hardly, available from other methods. This article reviews thermodynamic approaches and assays to measure collective properties such as heat adsorption / emission and volume variations. These methods can be successfully applied to the study of lipid vesicles (liposomes) and biological membranes. With respect to instrumentation, differential scanning calorimetry, pressure perturbation calorimetry, isothermal titration calorimetry, dilatometry, and acoustic techniques aimed at measuring the isothermal and adiabatic processes, two- and three-dimensional compressibilities are considered. Applications of these techniques to lipid systems include the measurement of different thermodynamic parameters and a detailed characterization of thermotropic, barotropic, and lyotropic phase behavior. The membrane binding and / or partitioning of solutes (proteins, peptides, drugs, surfactants, ions, etc.) can also be quantified and modeled. Many thermodynamic assays are available for studying the effect of proteins and other additives on membranes, characterizing non-ideal mixing, domain formation, bilayer stability, curvature strain, permeability, solubilization, and fusion. Studies of membrane proteins in lipid environments elucidate lipid-protein interactions in membranes. Finally, a plethora of relaxation phenomena toward equilibrium thermodynamic structures can be also investigated. The systems are described in terms of enthalpic and entropic forces, equilibrium constants, heat capacities, partial volume changes, volume and area compressibility, and so on, also shedding light on the stability of the structures and the molecular origin and mechanism of the structural changes.

Thermomechanic-electrical coupling in phospholipid monolayers near the critical point

Lipid monolayers have been shown to represent a powerful tool in studying mechanical and thermodynamic properties of lipid membranes as well as their interaction with proteins. Using Einstein's theory of fluctuations we here demonstrate that an experimentally derived linear relationship both between transition entropy S and area A as well as between transition entropy and charge q implies a linear relationships between compressibility κT, heat capacity cπ, thermal expansion coefficient αT, and electric capacity CT. We demonstrate that these couplings have strong predictive power as they allow calculating electrical and thermal properties from mechanical measurements. The precision of the prediction increases as the critical point TC is approached.

Thermodynamic study of the BRCT domain of BARD1 and its interaction with the -pSER-X-X-Phe- motif-containing BRIP1 peptide

The BRCA1-associated RING domain protein 1 (BARD1) is the heterodimeric partner of BRCA1. The BRCA1/BARD1 complex demonstrates ubiquitin ligase activity and has been implicated in genomic stability and tumor suppression. Both proteins possess a structurally conserved C-terminal domain (BRCT). While BRCA1-BRCT has been shown to mediate BRCA1 interactions with phosphoproteins such as BRIP1 by recognizing the pSer-X-X-Phe motif, attempts to demonstrate analogous interactions of its dimeric counterpart BARD1-BRCT, have so far been unsuccessful. In this study, chemical-denaturation experiments of BARD1-BRCT domain suggest that its low thermodynamic stability (DeltaG=2.5 kcal/mol) at room temperature, may affect some of its biochemical properties, such as its interaction with phosphopeptides. The stability of BARD1-BRCT domain at 10 degrees C, increases to 7.5 kcal/mol and isothermal titration calorimetry (ITC) experiments at this lower temperature showed binding to the BRIP1 phosphopeptide via an enthalpy-driven interaction, which appears to be specific to the pSer-X-X-Phe peptide-binding motif. Substitution of either pSer at position 0 with Ser (non-phosphorylated peptide) or Phe with Val at position +3, leads to no-binding ITC results. While these findings are indicative that BRIP1 is a potential BARD1 binding partner, it becomes evident that in vitro binding assays involving the entire BARD1 protein and in vivo experiments are also needed to establish its binding partners and its potential role in tumor suppression pathways.

Thiol-disulfide exchange between glutaredoxin and glutathione

Glutaredoxins are ubiquitous thiol-disulfide oxidoreductases which catalyze the reduction of glutathione-protein mixed disulfides. Belonging to the thioredoxin family, they contain a conserved active site CXXC motif. The N-proximal active site cysteine can form a mixed disulfide with glutathione or an intramolecular disulfide with the C-proximal cysteine. The C-proximal cysteine is not known to be involved in the catalytic mechanism. The stability of the mixed disulfide with glutathione has been investigated in detail using a mutant variant of yeast glutaredoxin 1, in which the C-proximal active site cysteine has been replaced with serine. The exchange reaction between the reduced protein and oxidized glutathione leading to formation of the mixed disulfide could readily be monitored by isothermal titration calorimetry (ITC) due to the enthalpic contributions from the noncovalent interactions and the protonation of glutathione thiolate. An algorithm for the analysis of this type of reaction by ITC was developed and showed that the interaction is enthalpy driven with a large entropy penalty. The applicability of the method was verified by a mass spectrometry-based approach, which gave a standard reduction potential of -295 mV for the mixed disulfide. In another set of experiments, the pK(a) value of the active site cysteine was determined. In line with what has been observed for other glutaredoxins, this cysteine was found to have a very low pK(a) value. The glutathionylation of glutaredoxin was shown to have a substantial effect on the thermal stability of the protein as revealed by differential scanning calorimetry.

Analysis of lipoplex structure and lipid phase changes

Efficient delivery of genetic material to cells is needed for tasks of utmost importance in the laboratory and clinic, such as gene transfection and gene silencing. Synthetic cationic lipids can be used as delivery vehicles for nucleic acids and are now considered the most promising nonviral gene carriers. They form complexes (lipoplexes) with the polyanionic nucleic acids. A critical obstacle for clinical application of the lipid-mediated DNA delivery (lipofection) is its unsatisfactory efficiency for many cell types. Understanding the mechanism of lipid-mediated DNA delivery is essential for their successful application, as well as for a rational design and synthesis of novel cationic lipoid compounds for enhanced gene delivery. A viewpoint now emerging is that the critical factor in lipid-mediated transfection is the structural evolution of lipoplexes within the cell, upon interacting and mixing with cellular lipids. In particular, recent studies showed that the phase evolution of lipoplex lipids upon interaction and mixing with membrane lipids appears to be decisive for transfection success: specifically, lamellar lipoplex formulations, which were readily susceptible to undergoing lamellar-nonlamellar phase transition upon mixing with cellular lipids and were found rather consistently associated with superior transfection potency, presumably as a result of facilitated DNA release. Thus, understanding the lipoplex structure and the phase changes upon interacting with membrane lipids is important for the successful application of the cationic lipids as gene carriers.

Quantitative analysis of small molecule-nucleic acid interactions with a biosensor surface and surface plasmon resonance detection

Surface plasmon resonance (SPR) technology with biosensor surfaces has become a widely-used tool for the study of nucleic acid interactions without any labeling requirements. The method provides simultaneous kinetic and equilibrium characterization of the interactions of biomolecules as well as small molecule-biopolymer binding. SPR monitors molecular interactions in real time and provides significant advantages over optical or calorimetic methods for systems with strong binding coupled to small spectroscopic signals and/or reaction heats. A detailed and practical guide for nucleic acid interaction analysis using SPR-biosensor methods is presented. Details of the SPR technology and basic fundamentals are described with recommendations on the preparation of the SPR instrument, sensor chips, and samples, as well as extensive information on experimental design, quantitative and qualitative data analysis and presentation. A specific example of the interaction of a minor-groove-binding agent with DNA is evaluated by both kinetic and steady-state SPR methods to illustrate the technique. Since the molecules that bind cooperatively to specific DNA sequences are attractive for many applications, a cooperative small molecule-DNA interaction is also presented.

Predicting ligand binding affinity with alchemical free energy methods in a polar model binding site

We present a combined experimental and modeling study of organic ligand molecules binding to a slightly polar engineered cavity site in T4 lysozyme (L99A/M102Q). For modeling, we computed alchemical absolute binding free energies. These were blind tests performed prospectively on 13 diverse, previously untested candidate ligand molecules. We predicted that eight compounds would bind to the cavity and five would not; 11 of 13 predictions were correct at this level. The RMS error to the measurable absolute binding energies was 1.8 kcal/mol. In addition, we computed "relative" binding free energies for six phenol derivatives starting from two known ligands: phenol and catechol. The average RMS error in the relative free energy prediction was 2.5 kcal/mol (phenol) and 1.1 kcal/mol (catechol). To understand these results at atomic resolution, we obtained x-ray co-complex structures for nine of the diverse ligands and for all six phenol analogs. The average RMSD of the predicted pose to the experiment was 2.0 A (diverse set), 1.8 A (phenol-derived predictions), and 1.2 A (catechol-derived predictions). We found that predicting accurate affinities and rank-orderings required near-native starting orientations of the ligand in the binding site. Unanticipated binding modes, multiple ligand binding, and protein conformational change all proved challenging for the free energy methods. We believe that these results can help guide future improvements in physics-based absolute binding free energy methods.

Substrate-induced conformational changes in the essential peripheral membrane-associated mannosyltransferase PimA from mycobacteria: implications for catalysis

Phosphatidyl-myo-inositol mannosyltransferase A (PimA) is an essential glycosyltransferase (GT) involved in the biosynthesis of phosphatidyl-myo-inositol mannosides (PIMs), which are key components of the mycobacterial cell envelope. PimA is the paradigm of a large family of peripheral membrane-binding GTs for which the molecular mechanism of substrate/membrane recognition and catalysis is still unknown. Strong evidence is provided showing that PimA undergoes significant conformational changes upon substrate binding. Specifically, the binding of the donor GDP-Man triggered an important interdomain rearrangement that stabilized the enzyme and generated the binding site for the acceptor substrate, phosphatidyl-myo-inositol (PI). The interaction of PimA with the beta-phosphate of GDP-Man was essential for this conformational change to occur. In contrast, binding of PI had the opposite effect, inducing the formation of a more relaxed complex with PimA. Interestingly, GDP-Man stabilized and PI destabilized PimA by a similar enthalpic amount, suggesting that they formed or disrupted an equivalent number of interactions within the PimA complexes. Furthermore, molecular docking and site-directed mutagenesis experiments provided novel insights into the architecture of the myo-inositol 1-phosphate binding site and the involvement of an essential amphiphatic alpha-helix in membrane binding. Altogether, our experimental data support a model wherein the flexibility and conformational transitions confer the adaptability of PimA to the donor and acceptor substrates, which seems to be of importance during catalysis. The proposed mechanism has implications for the comprehension of the peripheral membrane-binding GTs at the molecular level.

Differential scanning calorimetry

Differential scanning calorimetry (DSC) has emerged as a powerful experimental technique for determining thermodynamic properties of biomacromolecules. The ability to monitor unfolding or phase transitions in proteins, polynucleotides, and lipid assemblies has not only provided data on thermodynamic stability for these important molecules, but also made it possible to examine the details of unfolding processes and to analyze the characteristics of intermediate states involved in the melting of biopolymers. The recent improvements in DSC instrumentation and software have generated new opportunities for the study of the effects of structure and changes in environment on the behavior of proteins, nucleic acids, and lipids. This review presents some of the details of application of DSC to the examination of the unfolding of biomolecules. After a brief introduction to DSC instrumentation used for the study of thermal transitions, the methods for obtaining basic thermodynamic information from the DSC curve are presented. Then, using DNA unfolding as an example, methods for the analysis of the melting transition are presented that allow deconvolution of the DSC curves to determine more subtle characteristics of the intermediate states involved in unfolding. Two types of transitions are presented for analysis, the first example being the unfolding of two large synthetic polynucleotides, which display high cooperativity in the melting process. The second example shows the application of DSC for the study of the unfolding of a simple hairpin oligonucleotide. Details of the data analysis are presented in a simple spreadsheet format.

Adhesion of giant vesicles mediated by weak binding of sialyl-LewisX to E-selectin in the presence of repelling poly(ethylene glycol) molecules

Prior to establishing tight contact with the endothelium, cells such as leukocytes or cancer cells use the recognition between sialyl-LewisX ligands and E-selectin receptors to establish weak, reversible adhesion and to roll along the vessel wall. We study the physical aspects of this process by constructing a mimetic system that consists of a giant fluid vesicle with incorporated lipid-anchored sialyl-LewisX molecules that bind to E-selectin that is immobilized on the flat substrate. The vesicles also carry a certain fraction of repelling PEG2000 molecules. We analyze the equilibrium state of adhesion in detail by means of reflection interference contrast microscopy and find that the adhesion process relies purely on the formation of one or more adhesion domains within the vesicle-substrate contact zone. We find that the content of ligands in the vesicle must be above 5 mol % to establish specific contacts. All concentrations of sialyl-LewisX above 8 mol % provide a very similar final state of adhesion. However, the size and shape of the adhesion domains strongly depend on both the concentrations of E-selectin (0-3500 molecules/microm2) and PEG2000 (0-5 mol %). At 3500 E-selectin molecules/microm2 and small concentrations of PEG2000, the vesicle-substrate contact is maximized and fully occupied by a single adhesion domain. At concentrations of 5 mol %, PEG2000 completely impedes the specific binding to any substrate. Lastly, an increase in the adhesion strength is observed in systems with identical compositions if the reduced volume of the vesicles is larger.

Predicting absolute ligand binding free energies to a simple model site

A central challenge in structure-based ligand design is the accurate prediction of binding free energies. Here we apply alchemical free energy calculations in explicit solvent to predict ligand binding in a model cavity in T4 lysozyme. Even in this simple site, there are challenges. We made systematic improvements, beginning with single poses from docking, then including multiple poses, additional protein conformational changes, and using an improved charge model. Computed absolute binding free energies had an RMS error of 1.9 kcal/mol relative to previously determined experimental values. In blind prospective tests, the methods correctly discriminated between several true ligands and decoys in a set of putative binders identified by docking. In these prospective tests, the RMS error in predicted binding free energies relative to those subsequently determined experimentally was only 0.6 kcal/mol. X-ray crystal structures of the new ligands bound in the cavity corresponded closely to predictions from the free energy calculations, but sometimes differed from those predicted by docking. Finally, we examined the impact of holding the protein rigid, as in docking, with a view to learning how approximations made in docking affect accuracy and how they may be improved.

Thermal unfolding of human BRCA1 BRCT-domain variants

Missense mutations at the BRCT domain of human BRCA1 protein have been associated with an elevated risk for hereditary breast/ovarian cancer. They have been shown to affect the binding site and they have also been proposed to affect domain stability, severely hampering the protein's tumor suppressor function. In order to assess the impact of various such mutations upon the stability and the function of the BRCT domain, heat-induced denaturation has been employed to study the thermal unfolding of variants M1775R and R1699W, which have been linked with the disease, as well as of V1833M, which has been reported for patients with a family history. Calorimetric and circular dichroism results reveal that in pH 9.0, 5 mM borate buffer, 200 mM NaCl, analogously to wild type BRCT, all three variants undergo partial thermal unfolding to a denatured state, which retains most of the native's structural characteristics. With respect to wild-type BRCT, the mutation M1775R induces the most severe effects especially upon the thermostability, while R1699W also has a strong impact. On the other hand, the thermal unfolding of variant V1833M is only moderately affected relative to wild-type BRCT. Moreover, isothermal titration calorimetric measurements reveal that contrary to M1775R and R1699W variants, V1833M binds to BACH1 and CtIP phosphopeptides.

Defect formation of lytic peptides in lipid membranes and their influence on the thermodynamic properties of the pore environment

We present an experimental study of the pore formation processes of small amphipathic peptides in model phosphocholine lipid membranes. We used atomic force microscopy to characterize the spatial organization and structure of alamethicin- and melittin-induced defects in lipid bilayer membranes and the influence of the peptide on local membrane properties. Alamethicin induced holes in gel DPPC membranes were directly visualized at different peptide concentrations. We found that the thermodynamic state of lipids in gel membranes can be influenced by the presence of alamethicin such that nanoscopic domains of fluid lipids form close to the peptide pores, and that the elastic constants of the membrane are altered in their vicinity. Melittin-induced holes were visualized in DPPC and DLPC membranes at room temperature in order to study the influence of the membrane state on the peptide induced hole formation. Also differential scanning calorimetry was used to investigate the effect of alamethicin on the lipid membrane phase behaviour.

Protein phase diagrams II: nonideal behavior of biochemical reactions in the presence of osmolytes

In the age of biochemical systems biology, proteomics, and high throughput methods, the thermodynamic quantification of cytoplasmatic reaction networks comes into reach of the current generation of scientists. What is needed to efficiently extract the relevant information from the raw data is a robust tool for evaluating the number and stoichiometry of all observed reactions while providing a good estimate of the thermodynamic parameters that determine the molecular behavior. The recently developed phase-diagram method, strictly speaking a graphical representation of linkage or Maxwell Relations, offers such capabilities. Here, we extend the phase diagram method to nonideal conditions. For the sake of simplicity, we choose as an example a reaction system involving the protein RNase A, its inhibitor CMP, the osmolyte urea, and water. We investigate this system as a function of the concentrations of inhibitor and osmolyte at different temperatures ranging from 280 K to 340 K. The most interesting finding is that the protein-inhibitor binding equilibrium depends strongly on the urea concentration--by orders-of-magnitude more than expected from urea-protein interaction alone. Moreover, the m-value of ligand binding is strongly concentration-dependent, which is highly unusual. It is concluded that the interaction between small molecules like urea and CMP can significantly contribute to cytoplasmic nonideality. Such a finding is highly significant because of its impact on renal tissue where high concentrations of cosolutes occur regularly.

Chemical and Thermal Stability of Insulin: Effects of Zinc and Ligand Binding to the Insulin Zinc-Hexamer

PURPOSE

To study the correlation between the thermal and chemical stability of insulin formulations with various insulin hexamer ligands.

MATERIALS AND METHODS

The thermal stability was investigated using differential scanning calorimetry (DSC) and near-UV circular dichroism (NUV-CD). The formation of chemical degradation products was studied with reversed-phase and size-exclusion chromatography and mass spectrometry.

RESULTS

An excellent correlation between the thermal stabilization by ligand binding and the deamidation of Asn(B3) was observed. The correlation between thermal stability and the formation of covalent dimer and other insulin related products was less clear. Zinc was found to specifically increase the deamidation and covalent dimer formation rate when the insulin hexamer was not further stabilized by phenolic ligand. Thiocyanate alone had no effect on the thermal stability of the insulin zinc-hexamer but significantly improved the chemical stability at 37 degrees C. At low temperatures thiocyanate induced a conformational change in the insulin hexamer. NUV-CD thermal scans revealed that this effect decreased with temperature; when the thermal denaturation temperature was reached, the effect was eliminated.

CONCLUSIONS

Thermal stability can be used to predict the rate of Asn(B3) deamidation in human insulin. Chemical degradation processes that do not rely on the structural stability of the protein do not necessarily correlate to the thermal stability.

Specific DNA binding by the homeodomain Nkx2.5(C56S): detection of impaired DNA or unfolded protein by isothermal titration calorimetry

Titrations of specific 18-bp duplex DNA with the cardiac-specific homeodomain Nkx2.5(C56S) have utilized an ultrasensitive isothermal titration calorimeter (ITC). As the free DNA nears depletion, we observe large apparent decreases in the binding enthalpy when the DNA is impaired or when the temperature is sufficiently high to produce some unfolding of the free protein. Either effect can be attributed to refolding of the biopolymer that occurs as a result of stabilization due to the large favorable change in free energy on the homeodomain binding to DNA (-49.4 kJ/mol at 298 K). In either case, thermodynamic parameters obtained in such ITC experiments are unreliable. By using a lower temperature (85 vs. 95 degrees C) during the annealing of complementary DNA strands, damage of the 18-bp duplex DNA (T(m) = 72 degrees C) is avoided, and titrations with the homeodomain are normal at temperatures from 10 to 40 degrees C when >95% of the protein is folded. Under the latter conditions, the heat capacity plot is linear with a DeltaC(p) value of -0.80 +/- 0.03 kJ K(-1) mol(-1), which is more negative than that calculated from the burial of solvent accessible surface areas (-0.64 +/- 0.05 kJ K(-1) mol(-1)), consistent with water structures being at the protein-DNA interfaces.

Nonideal mixing in multicomponent lipid/detergent systems

A detailed understanding of the mixing properties of membranes to which detergents are added is mandatory for improving the application and interpretation of detergent based protein or lipid extraction assays. For Triton X-100 (TX-100), a nonionic detergent frequently used in the process of solubilizing and purifying membrane proteins and lipids, we present here a detailed study of the mixing properties of binary and ternary lipid mixtures by means of high-sensitivity isothermal titration calorimetry (ITC). To this end the partitioning thermodynamics of TX-100 molecules from the aqueous phase to lipid bilayers composed of various mixtures of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), egg-sphingomyelin (SM), and cholesterol (cho) are characterized. Composition-dependent partition coefficients K are analysed within the frame of a thermodynamic model developed to describe nonideal mixing in multicomponent lipid/detergent systems. The results imply that POPC, fluid SM, and TX-100 mix almost ideally (nonideality parameters |ρ(α/β)|<RT). However, favourable SM/cho (ρ(SM/cho)≤-6RT) and unfavourable PC/cho interactions (ρ(PC/cho) = 2RT) may under certain conditions cause POPC/TX-100-enriched domains to segregate from SM/cho-enriched ones. TX-100/cho contacts are unfavourable (ρ(cho/TX) = 4RT), so the system tends to avoid them. That means, addition of TX-100 promotes the separation of SM/cho-rich from PC/TX-100-rich domains. It appears that cho/detergent interactions are crucial governing the abundance and composition of detergent-resistant membrane patches.

Direct determination of the thermodynamics of polyelectrolyte complexation and implications thereof for electrostatic layer-by-layer assembly of multilayer films

Interpolyelectrolyte complex (IPEC) formation between poly(styrene sulfonate) (PSS) and poly(allylamine hydrochloride) (PAH) has been studied over a range of ionic strengths by isothermal titration calorimetry (ITC), turbidity titration, and electrostatic layer-by-layer assembly (ELBL). The results indicate that IPEC formation of PSS/PAH in aqueous solution is predominantly entropy-driven. The thermodynamic parameters suggest the formation of different types of complexes and aggregates due to salt-induced conformational changes in the polyelectrolyte conformation. Differences in polyelectrolyte behavior in the different salt-concentration regimes are described in terms of changes in the Debye screening length of the polyelectrolytes. The relationship of the results to the effect of salt concentration on the assembly of polyelectrolyte multilayer films (PEMs) is discussed.

Equilibrium heat-induced denaturation of chitinase 40 from Streptomyces thermoviolaceus

High-precision differential scanning calorimetry (DSC) and circular dichroism (CD) have been employed to study the thermal unfolding of chitinase 40 (Chi40) from Streptomyces thermoviolaceus. Chi40 belongs to family 18 of glycosyl hydrolase superfamily bearing a catalytic domain with a "TIM barrel"-like fold, which exhibits deviations from the (beta/alpha)8 fold. The thermal unfolding is reversible at pH = 8.0 and 9.0. The denatured state is characterized by extensive structural changes with respect to the native. The process is characterized by slow relaxation kinetics. Even slower refolding rates are recorded upon cooling. It is shown that the denaturation calorimetric data obtained at slow heating rate (0.17 K/min) are in excellent agreement with equilibrium data obtained by extrapolation of the experimental results to zero scanning rate. Analysis of the DSC results reveals that the experimental data can be successfully fitted using either a non-two-state sequential model involving one equilibrium intermediate, or an independent transitions model involving the unfolding of two Chi40 energetic domains to intermediate states. The stability of the native state with respect to the final denatured state is estimated, deltaG = 24.0 kcal/mol at 25 degrees C. The thermal results are in agreement with previous findings from chemical denaturation studies of a wide variety of (beta/alpha)8 barrel proteins, that their unfolding is a non-two-state process, always involving at least one unfolding intermediate.

Effect of pH on the structure and stability of Bacillus circulans ssp. alkalophilus phosphoserine aminotransferase: Thermodynamic and crystallographic studies

pH is one of the key parameters that affect the stability and function of proteins. We have studied the effect of pH on the pyridoxal-5'-phosphate-dependent enzyme phosphoserine aminotransferase produced by the facultative alkaliphile Bacillus circulans ssp. alkalophilus using thermodynamic and crystallographic analysis. Enzymatic activity assay showed that the enzyme has maximum activity at pH 9.0 and relative activity less than 10% at pH 7.0. Differential scanning calorimetry and circular dichroism experiments revealed variations in the stability and denaturation profiles of the enzyme at different pHs. Most importantly, release of pyridoxal-5'-phosphate and protein thermal denaturation were found to occur simultaneously at pH 6.0 in contrast to pH 8.5 where denaturation preceded cofactor's release by approximately 3 degrees C. To correlate the observed differences in thermal denaturation with structural features, the crystal structure of phosphoserine aminotransferase was determined at 1.2 and 1.5 A resolution at two different pHs (8.5 and 4.6, respectively). Analysis of the two structures revealed changes in the vicinity of the active site and in surface residues. A conformational change in a loop involved in substrate binding at the entrance of the active site has been identified upon pH change. Moreover, the number of intramolecular ion pairs was found reduced in the pH 4.6 structure. Taken together, the presented kinetics, thermal denaturation, and crystallographic data demonstrate a potential role of the active site in unfolding and suggest that subtle but structurally significant conformational rearrangements are involved in the stability and integrity of phosphoserine aminotransferase in response to pH changes.

Pressure-Modulated Differential Scanning Calorimetry: Theoretical Background

We demonstrate in this work and in the accompanying paper that it is possible to measure simultaneously heat capacity and expansibility of biomolecules in a single DSC experiment. In this study, we provide the theoretical basis for this new method based on rigorous statistical thermodynamics. The theoretical treatment presented here demonstrates that there are two additive contributions to the heat capacity at variable pressure, viz. (1) the isobaric heat capacity and (2) an expansibility term, and that these contributions can be experimentally separated to obtain simultaneously both heat capacity and expansibility in continuous DSC temperature scans preformed under pressure modulation. Equations that describe the mixed heat capacity and expansibility signal are derived, and experimental strategies as well as data extraction procedures are discussed.

Calorimetric and crystallographic analysis of the oligomeric structure of Escherichia coli GMP kinase

Guanosine monophosphate kinases (GMPKs), which catalyze the phosphorylation of GMP and dGMP to their diphosphate form, have been characterized as monomeric enzymes in eukaryotes and prokaryotes. Here, we report that GMPK from Escherichia coli (ecGMPK) assembles in solution and in the crystal as several different oligomers. Thermodynamic analysis of ecGMPK using differential scanning calorimetry shows that the enzyme is in equilibrium between a dimer and higher order oligomers, whose relative amounts depend on protein concentration, ionic strength, and the presence of ATP. Crystallographic structures of ecGMPK in the apo, GMP and GDP-bound forms were solved at 3.2A, 2.9A and 2.4A resolution, respectively. ecGMPK forms a hexamer with D3 symmetry in all crystal forms, in which the two nucleotide-binding domains are able to undergo closure comparable to that of monomeric GMPKs. The 2-fold and 3-fold interfaces involve a 20-residue C-terminal extension and a sequence signature, respectively, that are missing from monomeric eukaryotic GMPKs, explaining why ecGMPK forms oligomers. These signatures are found in GMPKs from proteobacteria, some of which are human pathogens. GMPKs from these bacteria are thus likely to form the same quaternary structures. The shift of the thermodynamic equilibrium towards the dimer at low ecGMPK concentration together with the observation that inter-subunit interactions partially occlude the ATP-binding site in the hexameric structure suggest that the dimer may be the active species at physiological enzyme concentration.

Oscillatory instabilities in phase separation of binary mixtures: Fixing the thermodynamic driving

Binary liquid mixtures can show pronounced oscillations in the differential scanning calorimeter signal for the specific heat and in the turbidity when phase separation is induced by continuously ramping the temperature. For a fixed ramp rate, i.e., a linear temporal drift of temperature, only a small number of oscillations have been observed. In the present manuscript we describe an experimental setup where simultaneous video-microscopy and shadow-graph measurements can be performed on mixtures subjected to an arbitrary temporal temperature evolution. In particular, it can be adjusted to fix the thermodynamic driving force, which characterizes the rate of change of the composition of the coexisting phases. With this novel technique both the number of oscillations and the temperature interval where oscillations are observed increase significantly. This technique can easily be applied to a great variety of binary mixtures, permitting detailed investigations of their phase-separation kinetics under slowly ramping temperature.

Formation of self-supporting reversible cellular networks in suspensions of colloids and liquid crystals

In mixtures of thermotropic liquid crystals with spherical poly (methyl methacrylate) particles, self-supporting networklike structures are formed during slow cooling past the isotropic-to-nematic phase transformation. To characterize the process of network formation in terms of morphology, phase transformation kinetics, and mechanical properties, we have combined data from polarization and laser scanning confocal microscopy with calorimetric, NMR, and rheological results. Our data suggest that the mechanism of network formation is dominated by a broadened temperature and time interval of phase transformation rather than by particle size or concentration. The observation that the width of the transformation interval strongly depends on sample preparation supports the hypothesis that a third component, most likely alkane remnants slowly liberated from the particles, plays a crucial role. In addition, calorimetric findings for liquid crystal/colloid mixtures, heated and cooled up to 13 times, point to separation of the liquid crystal into two compartments with different phase transformation kinetics. This could be explained by redistribution and enrichment of alkane in the particle-composed network walls. A further increase of the storage modulus, G', and incomplete dissolution of the networks in the isotropic state indicate that the formation of two compartments during repeated temperature cycles stabilizes the network and confers strong memory effects.

Folding of a designed simple ankyrin repeat protein

Ankyrin repeats (AR) are 33-residue motifs containing a beta-turn, followed by two alpha-helices connected by a loop. AR occur in tandem arrangements and stack side-by-side to form elongated domains involved in very different cellular tasks. Recently, consensus libraries of AR repeats were constructed. Protein E1_5 represents a member of the shortest library, and consists of only a single consensus repeat flanked by designed N- and C-terminal capping repeats. Here we present a biophysical characterization of this AR domain. The protein is compactly folded, as judged from the heat capacity of the native state and from the specific unfolding enthalpy and entropy. From spectroscopic data, thermal and urea-induced unfolding can be modeled by a two-state transition. However, scanning calorimetry experiments reveal a deviation from the two-state behavior at elevated temperatures. Folding and unfolding at 5 degrees C both follow monoexponential kinetics with k(folding) = 28 sec(-1) and k(unfolding) = 0.9 sec(-1). Kinetic and equilibrium unfolding parameters at 5 degrees C agree very well. We conclude that E1_5 folds in a simple two-state manner at low temperatures while equilibrium intermediates become populated at higher temperatures. A chevron-plot analysis indicates that the protein traverses a very compact transition state along the folding/unfolding pathway. This work demonstrates that a designed minimal ankyrin repeat protein has the thermodynamic and kinetic properties of a compactly folded protein, and explains the favorable properties of the consensus framework.

Thermal denaturation of the BRCT tandem repeat region of human tumour suppressor gene product BRCA1

Reduced stability of the tandem BRCT domains of human BReast CAncer 1 (BRCA1) due to missense mutations may be critical for loss of function in DNA repair and damage-induced checkpoint control. In the present thermal denaturation study of the BRCA1 BRCT region, high-precision differential scanning calorimetry (DSC) and circular dichroism (CD) spectroscopy provide evidence for the existence of a denatured state that is structurally very similar to the native. Consistency between theoretical structure-based estimates of the enthalpy (DeltaH) and heat capacity change (DeltaCp) and the calorimetric results is obtained when considering partial thermal unfolding contained in the region of the conserved hydrophobic pocket formed at the interface of the two BRCT repeats. The structural integrity of this region has been shown to be crucial for the interaction of BRCA1 with phosphorylated peptides. In addition, cancer-causing missense mutations located at the inter-BRCT-repeat interface have been linked to the destabilization of the tandem BRCT structure.

Evaluation of riboflavin binding protein domain interaction using differential scanning calorimetry

Riboflavin binding (or carrier) protein (RfBP) is a monomeric, two-domain protein, originally purified from hens' egg white. RfBP contains nine disulfide bridges; as a result, the protein forms a compact structure and undergoes reversible three-state thermal denaturation. This was demonstrated using a differential scanning calorimetry (DSC) method [Wasylewski M. (2000) J. Prot. Chem. 19(6), 523-528]. It has been shown that the RfBP complex with riboflavin denaturates in a three-state process which may be attributed to sequential unfolding of the RfBP domains. In case of apo RfBP, the ligand binding domain denaturates at a lower temperature than the C-terminal domain. Ligand binding greatly enhances the thermostability of the N-terminal domain, whereas the C-terminal domain thermostability is only slightly affected and, in case of the examined holo RfBPs, the denaturation peaks of both domains merge or cross over. The magnitude of the changes depends on ligand structure. A detailed study of protein concentration effects carried out in this work allowed to estimate not only the thermostability of both domains but also the strength of domain interactions. The DeltaCp, of denaturation was found for C-terminus and N-terminus of RfBP-riboflavin complex to amount to 2.5 and -1.9 kcal mol(-1), respectively. The calculated domain interaction free energy, DeltaGCN, was estimated to be approximately -1580 cal mol(-1) at 67.0 degrees C. This value indicates that the interdomain interaction is of medium strength.