Contrasting the denaturing effect of guanidinium chloride with the stabilizing effect of guanidinium sulfate

Guanidinium chloride, GdmCl, is a strong denaturing agent of globular proteins, whereas guanidinium sulfate, Gdm(2)SO(4), is a stabilizing agent of globular proteins. The stabilizing activity of Gdm(2)SO(4) is unexpected because the denaturant capability of GdmCl is due to direct interactions of Gdm(+) ions with protein surface groups. It is shown that the statistical thermodynamic approach devised to explain the molecular origin of cold denaturation [G. Graziano, Phys. Chem. Chem. Phys., 2010, 12, 14245-14252] can provide a rationalization of the different behaviour of GdmCl and Gdm(2)SO(4) towards globular proteins. The fundamental quantity is the reversible work to create in the aqueous solution a cavity suitable to host the D-state and a cavity suitable to host the N-state. In aqueous GdmCl solutions, this contribution is not large enough to overwhelm the conformational entropy gain upon unfolding and the direct attractions between Gdm(+) ions and protein surface groups; in aqueous Gdm(2)SO(4) solutions, it is so large that it overwhelms the two destabilizing contributions. Sulfate ions, due to their high charge density, interact strongly with water molecules producing a number density increase, that, in turn, renders the cavity creation process very costly, reversing the denaturing power of Gdm(+) ions and stabilizing the N-state of globular proteins.     

Nanometer-scale ion aggregates in aqueous electrolyte solutions: guanidinium carbonate

Neutron diffraction with isotopic substitution (NDIS) experiments and molecular dynamics (MD) simulations have been used to characterize the structure of aqueous guanidinium carbonate (Gdm2CO3) solutions. The MD simulations found very strong hetero-ion pairing in Gdm2CO3 solution and were used to determine the best structural experiment to demonstrate this ion pairing. The NDIS experiments confirm the most significant feature of the MD simulation, which is the existence of strong hetero-ion pairing between the Gdm+ and CO3(2-) ions. The neutron structural data also support the most interesting feature of the MD simulation, that the hetero-ion pairing is sufficiently strong as to lead to nanometer-scale aggregation of the ions. The presence of such clustering on the nanometer length scale was then confirmed using small-angle neutron scattering experiments. Taken together, the experiment and simulation suggest a molecular-level explanation for the contrasting denaturant properties of guanidinium salts in solution.     

Nanometer-scale ion aggregates in aqueous electrolyte solutions: guanidinium sulfate and guanidinium thiocyanate

Neutron diffraction experiments and molecular dynamics simulations are used to study the structure of aqueous solutions of two electrolytes: guanidinium sulfate (a mild protein conformation stabilizer) and guanidinium thiocyanate (a powerful denaturant). The MD simulations find the unexpected result that in the Gdm2SO4 solution the ions aggregated into mesoscopic (nanometer-scale) clusters, while no such aggregation is found in the GdmSCN solution. The neutron diffraction studies, the most direct experimental probe of solution structure, provide corroborating evidence that the predicted very strong ion pairing does occur in solutions of 1.5 m Gdm2SO4 but not in 3 m solutions of GdmSCN. A mechanism is proposed as to how this mesoscopic solution structure affects solution denaturant properties and suggests an explanation for the Hofmeister ordering of these solutions in terms of this ion pairing and the ability of sulfate to reverse the denaturant power of guanidinium.     

Complex ion effects on polypeptide conformational stability: chloride and sulfate salts of guanidinium and tetrapropylammonium

The effects of chloride and sulfate salts of tetrapropylammonium (TPA(+)) and guanidinium (Gdm(+)) on the conformational stabilities of tryptophan zipper (trpzip) and α-helical (alahel) peptides were measured by circular dichroism spectroscopy. Like Gdm(+), TPA(+) interacts with the planar tryptophan indole group, perturbing the conformational stability of trpzip peptides. TPA(+) effects are largely unaffected by sulfate, indicating an absence of the heteroion pairing that is observed in concentrated Gdm(2)SO(4) solutions. TPA(+) stabilizes helical conformations in alahel peptides, indicating exclusion from the peptide bond. The observations are broadly consistent with predictions of molecular dynamics simulations [Mason, P. E.; et al. J. Phys. Chem. B2009, 113, 3227-3234], indicating that the effects of complex ions on proteins are increasingly predictable in terms of ion hydration, complementary interactions with specific protein groups, and ion-pairing contributions.     

The reversal by sulfate of the denaturant activity of guanidinium

Guanidinium (Gdm+) chloride is a powerful protein denaturant, whereas the sulfate dianion (SO42-) is a strong stabilizer of folded protein states; Gdm2SO4 is effectively neutral in its effects on protein stability. While the "neutralizing" effects of protein-stabilizing solutes on the activity of denaturants can be broadly interpreted in terms of additive effects of the solutes, recent experimental and simulation studies support a role for hetero-ion interactions in the effect of sulfate on Gdm+ denaturation [Mason, P. E.; et al. J. Phys. Chem. B 2005, 109, 24185-24196]. Here we describe an experimental strategy for testing this mechanism that involves spectroscopic analysis of the separate effects of alkali metal sulfates (Na2SO4, Rb2SO4), GdmCl, and Gdm2SO4 on the folded populations of several peptides chosen to dissect specific noncovalent contributions to the conformational stability of proteins [alanine-based helical peptides stabilized by hydrogen bonds, tryptophan zipper (trpzip) peptides stabilized largely by cross-strand indole-indole interactions]. While the trpzip peptides are highly sensitive to GdmCl denaturation, they are unaffected by NaCl, Na2SO4, or Gdm2SO4, indicating that the reversal of the denaturant activity of Gdm+ by sulfate in this case is not due to competing stabilizing (sulfate) and destabilizing (Gdm+) interactions. Gdm2SO4 was found to retain considerable denaturant activity against alanine-based alpha-helical peptides. The differences in the effects of Gdm2SO4 on the two peptide types can be understood in terms of the different mechanisms of Gdm+ denaturation of trpzip peptides and helical peptides, respectively, and the specific nature of Gdm+ and SO42- ionic "clustering" that differentially affects the ability of Gdm+ to make the molecular interactions with the peptides that underlie its denaturant activity.     

Insight into the effects of graphene oxide sheets on the conformation and activity of glucose oxidase: towards developing a nanomaterial-based protein conformation assay

Protein conformation associates with particular properties of proteins and relates to protein-mediated diseases. Detailed elucidation of secondary and tertiary formation, stability, and the structural and dynamic properties of proteins has been one of the main topics studied in chemistry and biology. In this work, the conformation changes in glucose oxidase (GOx) induced by the graphene oxide (GO) sheets were studied in detail by various spectroscopic techniques including ultraviolet-visible (UV-vis) absorption, fluorescence, and circular dichroism (CD) spectroscopy. The results indicated that GOx underwent substantial conformation changes after assembling on the surface of GO. The interaction of GOx with GO could induce the exposure of the FAD (flavin adenine dinucleotide) moiety to solvent and transfer tryptophan (Trp) residues to a more hydrophobic environment. The calculation from CD spectra showed that GO could induce the conversion of α-helix to β-sheet structures, even unfolding of the protein. These alterations in the conformation of GOx resulted in a significant decrease in the catalytic activity of the enzyme in glucose oxidation. The possible reasons for these conformation changes in GOx are also discussed. This work not only provides insight into the interaction between atomically flat graphitic structures and various biological systems but also creates a framework for analyzing the biosafety of nanomaterials in terms of the biological behavior of biomacromolecules. The results obtained here can direct the further applications of the nanomaterials.     

Guanidinium chloride molecular diffusion in aqueous and mixed water-ethanol solutions

Solutions containing guanidinium chloride (GdmCl), or equivalently guanidine hydrochloride (GdnHCl), are commonly used to denature macromolecules such as proteins and DNA in, for example, microfluidics studies of protein unfolding. To design and study such applications, it is necessary to know the diffusion coefficients for GdmCl in the solution. To this end, we use molecular dynamics simulations to calculate the diffusion coefficients of GdmCl in water and in water-ethanol solutions, for which no direct experimental measurements exist. The fully atomistic simulations show that the guandinium cation Gdm (+) diffusion decreases as the concentration of both Gdm (+) and ethanol in the solution increases. The simulations are validated against available literature data, both transformed measured viscosity values and computed diffusion coefficients, and we show that a prudent choice of water model, namely TIP4P-Ew, gives calculated diffusion coefficients in good agreement with the transformed measured viscosity values. The calculated Gdm (+) diffusion behavior is explained as a dynamic mixture of free cation, stacked cation, and ion-paired species in solution, with weighted contributions to Gdm (+) diffusion from the stacked and paired states helping explain measured viscosity data in terms of atom-scale dynamics.     

Nonisomerizable non-retinal chromophores initiate light-induced conformational alterations in bacterioopsin.

The photoactivation of retinal proteins is usually interpreted in terms of C=C photoisomerization of the retinal moiety, which triggers appropriate conformational changes in the protein. In this work several dye molecules, characterized by a completely rigid structure in which no double-bond isomerization is possible, were incorporated into the binding site of bacteriorhodopsin (bR). Using a light-induced chemical reaction of a labeled EPR probe, it was observed that specific conformational alterations in the protein are induced following light absorption by the dye molecules occupying the binding site. The exact nature of these changes and their relationship to those occurring in the bR photocycle are still unclear. Nevertheless, their occurrence proves that C=C or C=NH(+) isomerization is not a prerequisite for protein conformational changes in a retinal protein. More generally, we show that conformational changes, leading to changes in reactivity, may be induced in proteins by optical excitation of simple nonisomerizable dyes located in the macromolecular matrix.     

Fluorescence probe of Trp-cage protein conformation in solution and in gas phase

Measurements of protein unfolding in the absence of solvent, when combined with unfolding studies in solution, offer a unique opportunity to measure the effects of solvent on protein structure and dynamics. The experiments presented here rely on the fluorescence of an attached dye to probe the local conformational dynamics through interactions with a Trp residue and fields originating on charge sites. We present fluorescence measurements of thermal fluctuations accompanying conformational change of a miniprotein, Trp-cage, in solution and in gas phase. Molecular dynamics (MD) simulations are performed as a function of temperature, charge state, and charge location to elucidate the dye-protein conformational dynamics leading to the changes in measured fluorescence. The results indicate that the stability of the unsolvated protein is dominated by hydrogen bonds. Substituting asparagine for aspartic acid at position 9 results in a dramatic alteration of the solution unfolding curve, indicating that the salt bridge involving Lys8, Asp9, and Arg16 (+ - +) is essential for Trp-cage stability in solution. In contrast, this substitution results in minor changes in the unfolding curve of the unsolvated protein, showing that hydrogen bonds are the major contributor to the stability of Trp-cage in gas phase. Consistent with this hypothesis, the decrease in the number of hydrogen bonds with increasing temperature indicated by MD simulations agrees reasonably well with the experimentally derived enthalpies of conformational change. The simulation results display relatively compact conformations compared with NMR structures that are generally consistent with experimental results. The measured unfolding curves of unsolvated Trp-cage ions are invariant with the acetonitrile content of the solution from which they are formed, possibly as a result of conformational relaxation during or after desolvation. This work demonstrates the power of combined solution and gas-phase studies and of single-point mutations to identify specific noncovalent interactions which contribute to protein-fold stability. The combination of experiment and simulation is particularly useful because these approaches yield complementary information which can be used to deduce the details of structural changes of proteins in the gas phase.     

Engineering the interface between glucose oxidase and nanoparticles

The behavior of glucose oxidase (GOx) on gold nanoparticles (NPs) was investigated as a function of (1) NP surface chemistry, (2) stabilizing protein additives, and (3) protein microenvironment. GOx secondary structure and unfolding was probed by circular dichroism (CD) spectroscopy and fluorescence, and GOx enzymatic activity was measured by a colorimetric assay. We also examined the activity and structure of GOx after displacement from the NP surface. Generally, GOx behavior was negatively impacted by conjugation to the NP, and conjugation conditions could vary the influence of the NP. Surface chemistry and protein microenvironment could improve behavior, but addition of stabilizing proteins negatively influenced activity. After displacement from the NPs, GOx tended to remain unfolded, indicating that the interactions with the NP were irreversible.     

Complex interactions between molecular ions in solution and their effect on protein stability

Protein stability in ionic solutions depends on the delicate balance between protein-ion and ion-ion interactions. For molecular ions containing multiple charged groups, the role of ion-ion interactions is particularly important. In this study, we show how the interplay between homo- and heteroion pairing influences protein stability using polyarginine salts as a model system. For the chloride salts, protein thermostability decreases as the size of the peptide increases, indicating enhanced binding to the protein. Moreover, it indicates reduced homoion pairing between Gdm(+) and carboxylate groups that is largely responsible for aggregation suppression, rather than denaturation, in monomeric arginine solutions. However, for the sulfate salts, strong heteroion pairing between the Gdm(+) groups and the sulfate counterions compensates for the loss of homoion pairing and, in return, leads to enhanced thermostability and a dramatically reduced (up to 10-30 times) rate of protein aggregation. Molecular dynamics simulations reveal how this ion pairing enhances conformational stability and, at the same time, reduces protein association. This study provides insight into complex ion effects on protein stability and serves as an example of how these intrasolvent interactions can be leveraged to enhance protein stability.     

Dissection of the difference between the group I metal ions in inhibiting GSK3β: a computational study

Glycogen synthase kinase 3β (GSK3β) is a serine/threonine kinase that requires two cofactor Mg(2+) ions for catalysis in regulating many important cellular signals. Experimentally, Li(+) is a competitive inhibitor of GSK3β relative to Mg(2+), while this mechanism is not experienced with other group I metal ions. Herein, we use native Mg(2)(2+)-Mg(1)(2+) GSK3β and its Mg(2)(2+)-M(1)(+) (M = Li, Na, K, and Rb) derivatives to investigate the effect of metal ion substitution on the mechanism of inhibition through two-layer ONIOM-based quantum mechanics/molecular mechanics (QM/MM) calculations and molecular dynamics (MD) simulations. The results of ONIOM calculations elucidate that the interaction of Na(+), K(+), and Rb(+) with ATP is weaker compared to that of Mg(2+) and Li(+) with ATP, and the critical triphosphate moiety of ATP undergoes a large conformational change in the Na(+), K(+), and Rb(+) substituted systems. As a result, the three metal ions (Na(+), K(+), and Rb(+)) are not stable and depart from the active site, while Mg(2+) and Li(+) can stabilize in the active site, evident in MD simulations. Comparisons of Mg(2)(2+)-Mg(1)(2+) and Mg(2)(2+)-Li(1)(+) systems reveal that the inline phosphor-transfer of ATP and the two conserved hydrogen bonds between Lys85 and ATP, together with the electrostatic potential at the Li(1)(+) site, are disrupted in the Mg(2)(2+)-Li(1)(+) system. These computational results highlight the possible mechanism why Li(+) inhibits GSK3β.     

Infrared fingerprint spectroscopy and theoretical studies of potassium ion tagged amino acids and peptides in the gas phase

Infrared multiple-photon dissociation spectroscopy is effected on the K(+) tagged aromatic amino acids tyrosine and phenylalanine, as well as the K(+) tagged peptides bradykinin fragment 1-5 and [Leu]-enkephalin. The fingerprint (800-1800 cm(-1)) infrared spectra of these species are compared to density-functional theory (DFT) calculated spectra to determine whether the complex is in the charge solvation (CS) or salt bridge (SB) (i.e. zwitterionic) configuration. For the aromatic amino acids the CS structure is favored and the tridentate N/O/ring structure is found to be the preferred binding geometry for K(+). The experimental and theoretical evidence for bradykinin fragment 1-5 tagged with K(+) suggests that the SB structure is favored; the calculations indicate a head-to-tail looped structure stabilized by a salt bridge between the protonated guanidine group and the deprotonated C-terminus, which allows K(+) to sit in a binding pocket with five C=O electrostatic interactions. For K(+) tagged [Leu]-enkephalin the spectroscopic evidence is not as clear. While the calculations clearly favor a CS structure and the observation of a weak carboxylic acid C=O stretching band in the infrared spectrum matches this finding, the prominence of a band at 1600 cm(-1) renders the analysis more ambiguous, and hence the presence of some salt bridge ions cannot be excluded. Another striking feature in the [Leu]-enkephalin spectrum is the high infrared activity of the tyrosine side-chain modes, which can be clearly identified from comparison to the [Tyr + K](+) experimental spectrum, but which is not reproduced by the DFT calculations.     

Parametrization, molecular dynamics simulation, and calculation of electron spin resonance spectra of a nitroxide spin label on a polyalanine alpha-helix

The nitroxide spin label 1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl-methanethiosulfonate (MTSSL), commonly used in site-directed spin labeling of proteins, is studied with molecular dynamics (MD) simulations. After developing force field parameters for the nitroxide moiety and the spin label linker, we simulate MTSSL attached to a polyalanine alpha-helix in explicit solvent to elucidate the factors affecting its conformational dynamics. Electron spin resonance spectra at 9 and 250 GHz are simulated in the time domain using the MD trajectories and including global rotational diffusion appropriate for the tumbling of T4 Lysozyme in solution. Analysis of the MD simulations reveals the presence of significant hydrophobic interactions of the spin label with the alanine side chains.     

Salt-specific stability and denaturation of a short salt-bridge-forming alpha-helix

The structure of a single alanine-based Ace-AEAAAKEAAAKA-Nme peptide in explicit aqueous electrolyte solutions (NaCl, KCl, NaI, and KF) at large salt concentrations (3-4 M) is investigated using approximately 1 mus molecular dynamics (MD) computer simulations. The peptide displays 71% alpha-helical structure without salt and destabilizes with the addition of NaCl in agreement with experiments of a somewhat longer version. It is mainly stabilized by direct and indirect (" i + 4")EK salt bridges between the Lys and Glu side chains and a concomitant backbone shielding mechanism. NaI is found to be a stronger denaturant than NaCl, while the potassium salts hardly show influence. Investigation of the molecular structures reveals that consistent with recent experiments Na (+) has a much stronger affinity to side chain carboxylates and backbone carbonyls than K (+), thereby weakening salt bridges and secondary structure hydrogen bonds. At the same time, the large I (-) has a considerable affinity to the nonpolar alanine in line with recent observations of a large propensity of I (-) to adsorb to simple hydrophobes, and thereby "assists" Na (+) in its destabilizing action. In the denatured states of the peptide, novel long-lived (10-20 ns) "loop" configurations are observed in which single Na (+) ions and water molecules are hydrogen-bonded to multiple backbone carbonyls. In an attempt to analyze the denaturation behavior within the preferential interaction formalism, we find indeed that for the strongest denaturant, NaI, the protein is least hydrated. Additionally, a possible indication for protein denaturation might be a preferential solvation of the peptide backbone by the destabilizing cosolute (sodium). The mechanisms found in this work may be of general importance to understand salt effects on protein secondary structure stability.     

Synthesis,structure, spectral and coordination properties of a crown ether derivative of 1,3-indandione. A new structural evidence for the versatile reactivity of 2-acetyl-1,3-indandione

The current study describes the synthesis of a new 1,3-indandione derivative with conjugated N-phenylaza-15-crown-5 moiety (4). The crystal structure of compound 4 was solved and its optical properties were studied in various solvents and in presence of alkaline and alkaline-earth metal ions. Quantum chemical (DFT) methods were employed to describe the structure and the optical properties of the studied compound and its complexes. The obtained results indicated that the synthesis of compound 4 using acid-catalyzed aldol reaction between 2-acetyl-1,3-indandione and the corresponding aldehyde is accompanied with an unexpected deacetylation step. In this way, the N-phenylaza-15-crown-5 moiety is directly conjugated with the 1,3-indandione fragment, known as a very strong electron acceptor. Therefore, the absorption spectra of 4 are only slightly influenced by complexation with Ba²⁺ and Sr²⁺ ions.     

Identification of active sites of biomolecules. 1. Methyl-alpha-mannopyranoside and Fe(III)

Car-Parrinello molecular dynamics (CPMD) simulations, DFT chemical reactivity index calculations, and mass spectrometric measurements are combined in an integrated effort to elucidate the details of the coordination of a transition-metal ion to a carbohydrate. The impact of the interaction with the FeIII ion on the glycosidic linkage conformation of methyl-alpha-d-mannopyranoside is studied by classical molecular dynamics (MD) and CPMD simulations. This study shows that FeIII interacts with specific hydroxyl oxygen atoms of the carbohydrate, affecting the ground state carbohydrate conformation. These conformational details are discussed in terms of a set of supporting experiments involving electrospray ionization mass spectrometry, and CPMD simulations clearly indicate that the specific conformational preference is due to intramolecular hydrogen bonding. Classical MD simulations proved insensitive to these important chemical properties. Thus, we demonstrate the importance of chemical reactivity calculations and CPMD simulations in predicting the active sites of biological molecules toward metal cations.     

Solvent-induced conformational changes of polypeptides probed by electrospray-ionization mass spectrometry.

Electrospray-ionization (ESI) mass spectrometry is used to monitor higher order structural changes of polypeptides induced by alteration of the pH or organic solvent composition in the protein solution environment. A bimodal charge-state distribution is observed in the ESI mass spectrum of ubiquitin (relative molecular mass 8565) in solutions containing small amounts (less than 20%) of organic solvents. The distribution of peaks at high m/z (low-charge state) is found to represent the protein in its native, globular state; the higher-charge-state distribution is characteristic for a more extended conformation. Addition of methanol denaturant in excess of 40% v/v is needed to eliminate the low-charge-state distribution completely. Lesser amounts of acetonitrile, acetone, or isopropanol (approximately 20%) are required to denature the ubiquitin protein. Other proteins showing conformational effects in their ESI mass spectra are also illustrated. While the ESI spectra are related to solution phase structure, ESI-tandem mass spectrometry of multiply charged molecular ions of different conformation is suggested as a probe of gas-phase protein three-dimensional structure.     

Structural changes of ultrasonicated bovine serum albumin revealed by hydrogen–deuterium exchange and mass spectrometry

The structural changes of bovine serum albumin (BSA) under high-intensity ultrasonication were investigated by fluorescence spectroscopy and mass spectrometry. Evidence for the ultrasonication-induced conformational changes of BSA was provided by the intensity changes and maximum-wavelength shift in fluorescence spectrometry. Matrix-assisted laser desorption–ionization time-of-flight mass spectroscopy (MALDI-TOF MS) revealed the increased intensity of the peak at the charge state +5 and a newly emerged peak at charge state +6, indicating that the protein became unfolded after ultrasonication. Prevalent unfolding of BSA after ultrasonication was revealed by hydrogen–deuterium exchange coupled with mass spectrometry (HDX-MS). Increased intensity and duration of ultrasonication further promoted the unfolding of the protein. The unfolding induced by ultrasonication goes through an intermediate state similar to that induced by a low concentration of denaturant.     

Role of ions on structure and stability of a synthetic gramicidin ion channel in solution. A molecular dynamics study

We performed a molecular dynamics (MD) simulation to the investigate structure and stability of a synthetic gramicidin-like peptide in solution with and without ions. The starting structures of the MD simulations were taken from two recently solved NMR structures of this peptide in isotropic solution, which forms stable monomers or dimers in the presence or absence of ions, respectively. The monomeric structure is channel-like and is assumed to be stabilized by the presence of two Cs(+) ions bound in the channel, each one close to one channel entrance. In our MD simulations, we observed how the Cs(+) ions bind in the channel formed by the monomeric gramicidin-like peptide using implicit solvent and explicit ions with a concentration of 2 M. MD simulations were performed with and without explicit ions but with an implicit solvent model defined by the generalized Born approximation, which was used to mimic the dielectric properties of the solvent and to speed up the computations.