Equilibrium constants and free energies in unfolding of proteins in urea solutions

A novel thermodynamic approach to the reversible unfolding of proteins in aqueous urea solutions has been developed based on the premise that urea ligands are bound cooperatively to the macromolecule. When successive stoichiometric binding constants have values larger than expected from statistical effects, an equation for moles of bound urea can be derived that contains imaginary terms. For a very steep unfolding curve, one can then show that the fraction of protein unfolded, B̄, depends on the square of the urea concentration, U, and is given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}\bar {B}=\frac{{\mathit{A}}^{{\mathit{2}}}_{{\mathit{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}n\bar {B}}{\mathit{U}}^{{\mathit{2}}}}{{\mathrm{1\hspace{.167em}+\hspace{.167em}}}{\mathit{A}}^{{\mathrm{2}}}_{{\mathrm{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}\bar {B}}{\mathit{U}}^{{\mathrm{2}}}}{\mathrm{.}}\end{equation*}\end{document} A12 is the binding constant as B̄→ 0, and λ is a parameter that reflects the augmentation in affinities of protein for urea as the moles bound increases to the saturation number, n. This equation provides an analytic expression that reproduces the unfolding curve with good precision, suggests a simple linear graphical procedure for evaluating A12 and λ, and leads to the appropriate standard free energy changes. The calculated ΔG° values reflect the coupling of urea binding with unfolding of the protein. Some possible implications of this analysis to protein folding in vivo are described.

Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities

The alanine helix provides a model system for studying the energetics of interaction between water and the helical peptide group, a possible major factor in the energetics of protein folding. Helix formation is enthalpy-driven (−1.0 kcal/mol per residue). Experimental transfer data (vapor phase to aqueous) for amides give the enthalpy of interaction with water of the amide group as ≈−11.5 kcal/mol. The enthalpy of the helical peptide hydrogen bond, computed for the gas phase by quantum mechanics, is −4.9 kcal/mol. These numbers give an enthalpy deficit for helix formation of −7.6 kcal/mol. To study this problem, we calculate the electrostatic solvation free energy (ESF) of the peptide groups in the helical and β-strand conformations, by using the delphi program and parse parameter set. Experimental data show that the ESF values of amides are almost entirely enthalpic. Two key results are: in the β-strand conformation, the ESF value of an interior alanine peptide group is −7.9 kcal/mol, substantially less than that of N-methylacetamide (−12.2 kcal/mol), and the helical peptide group is solvated with an ESF of −2.5 kcal/mol. These results reduce the enthalpy deficit to −1.5 kcal/mol, and desolvation of peptide groups through partial burial in the random coil may account for the remainder. Mutant peptides in the helical conformation show ESF differences among nonpolar amino acids that are comparable to observed helix propensity differences, but the ESF differences in the random coil conformation still must be subtracted.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Deuterium isotope effect on the intramolecular electron transfer in Pseudomonas aeruginosa azurin

Intramolecular electron transfer in azurin in water and deuterium oxide has been studied over a broad temperature range. The kinetic deuterium isotope effect, kH/kD, is smaller than unity (0.7 at 298 K), primarily caused by the different activation entropies in water (−56.5 J K−1 mol−1) and in deuterium oxide (−35.7 J K−1 mol−1). This difference suggests a role for distinct protein solvation in the two media, which is supported by the results of voltammetric measurements: the reduction potential (E0′) of Cu2+/+ at 298 K is 10 mV more positive in D2O than in H2O. The temperature dependence of E0′ is also different, yielding entropy changes of −57 J K−1 mol−1 in water and −84 J K−1 mol−1 in deuterium oxide. The driving force difference of 10 mV is in keeping with the kinetic isotope effect, but the contribution to ΔS‡ from the temperature dependence of E0′ is positive rather than negative. Isotope effects are, however, also inherent in the nuclear reorganization Gibbs free energy and in the tunneling factor for the electron transfer process. A slightly larger thermal protein expansion in H2O than in D2O (0.001 nm K−1) is sufficient both to account for the activation entropy difference and to compensate for the different temperature dependencies of E0′. Thus, differences in driving force and thermal expansion appear as the most straightforward rationale for the observed isotope effect.

DNA-strand exchange promoted by RecA protein in the absence of ATP: implications for the mechanism of energy transduction in protein-promoted nucleic acid transactions.

DNA-strand exchange promoted by Escherichia coli RecA protein normally requires the presence of ATP and is accompanied by ATP hydrolysis, thereby implying a need for ATP hydrolysis. Previously, ATP hydrolysis was shown not to be required; here we demonstrate furthermore that a nucleoside triphosphate cofactor is not required for DNA-strand exchange. A gratuitous allosteric effector consisting of the noncovalent complex of ADP and aluminum fluoride, ADP.AIF4-, can both induce the high-affinity DNA-binding state of RecA protein and support the homologous pairing and exchange of up to 800-900 bp of DNA. These results demonstrate that induction of the functionally active, high-affinity DNA-binding state of RecA protein is needed for RecA protein-promoted DNA-strand exchange and that there is no requirement for a high-energy nucleotide cofactor for the exchange of DNA strands. Consequently, the free energy needed to activate the DNA substrates for DNA-strand exchange is not derived from ATP hydrolysis. Instead, the needed free energy is derived from ligand binding and is transduced to the DNA via the associated ligand-induced structural transitions of the RecA protein-DNA complex; ATP hydrolysis simply destroys the effector ligand. This concept has general applicability to the mechanism of energy transduction by proteins.

Tuning the metal-to-metal charge transfer energy of cyano-bridged dinuclear complexes

The metal-to-metal charge transfer (MMCT) transitions of a series of Class II mixed valence dinuclear complexes bearing cyano bridging ligands may be varied systematically by variations to either the hexacyanometallate(II) donor or Co-III acceptor moieties. Specifically, the new dinuclear species trans-[(LCoNCFe)-Co-14S(CN)(5)](-) (L-14S = 6-methyl-1,11-diaza-4,8-dithia- cyclotetradecane-6-amine) and trans-[(LCoNCRu)-Co-14(CN)(5)]-(L-14 = 6-methyl-1,4,8,11-tetraazacyclotetradecane-6-amine) have been prepared and their spectroscopic and electrochemical properties are compared with the relative trans-[(LCoNCFe)-Co-14(CN)(5)](-). The crystal structures of Na{trans-[(LCoNCFe)-Co-14S(CN)(5)]}.51/2H(2)O.1/2EtOH, Na{trans-[(LCoNCRu)-Co-14(CN)(5)]}.3H(2)O and Na{trans-[(LCoNCRu)-Co-14(CN)(5)]}.8H(2)O are also reported. The ensuing changes to the MMCT energy have been examined within the framework of Hush theory, and it was found that the free energy change between the redox isomers was the dominant effect in altering the energy of the MMCT transition.

Charge-balanced copolymer hydrogels:potential as biomedical materials

Polyzwitterionic-containing hydrogel materials been proposed for use in biomaterial applications. Polyzwitterions contain anions and cations in the same monomeric unit, unlike polyampholytes which contain them in different monomeric units. The use of cationic and anionic monomers in stoichiometrically equivalent proportions produces charge-balanced polyampholytes (PA) copolymers. Membranes prepared using either betaine-containing (BT) polyzwitterionic copolymers or PA copolymers can share similar properties, but the range of EWCs offered by membranes incorporating BT and PA monomers is greater than that for conventional neutral hydrogels and methacrylic acid-based systems. Here we compare properties of BT-containing and PA-containing copolymer membranes, relevant to their potential as biomedical materials. Membranes of the copolymers were prepared as previously described. Surface energy was determined using a GBX Digidrop (GBX Scientific Instruments), with diidomethane and water as probes. The absorption of proteins was determined by soaking the membranes in 1mg/ml protein solutions for a predetermined time, and measuring UV absorption of the membranes at certain wavelengths. The BT and PA copolymer membranes displayed similar values for the polar components and dispersive components of total surface free energy. This was perhaps not surprising when the structures of the monomers were considered. The BT and PA copolymer membranes displayed differences in their protein absorption over time, with the PA demonstrating higher uptake of protein than the BT. In addition to the aforementioned greater EWC range, the use of BT and PA copolymer membranes also avoids some of the problems associated with net anionicity. Comparison of the BT copolymer with the “pseudo” zwitterionic PA copolymers shows that controlled molecular architecture is required to gain the benefits of balancing the charges present in the copolymers in a way that will make them beneficial to hydrogel design.

Applicability of contact angle techniques used in the analysis of contact lenses, part 1:comparative methodologies

Objectives and Methods: Contact angle, as a representative measure of surface wettability, is often employed to interpret contact lens surface properties. The literature is often contradictory and can lead to confusion. This literature review is part of a series regarding the analysis of hydrogel contact lenses using contact angle techniques. Here we present an overview of contact angle terminology, methodology, and analysis. Having discussed this background material, subsequent parts of the series will discuss the analysis of contact lens contact angles and evaluate differences in published laboratory results. Results: The concepts of contact angle, wettability and wetting are presented as an introduction. Contact angle hysteresis is outlined and highlights the advantages in using dynamic analytical techniques over static methods. The surface free energy of a material illustrates how contact angle analysis is capable of providing supplementary surface characterization. Although single values are able to distinguish individual material differences, surface free energy and dynamic methods provide an improved understanding of material behavior. The frequently used sessile drop, captive bubble, and Wilhelmy plate techniques are discussed. Their use as both dynamic and static methods, along with the advantages and disadvantages of each technique, is explained. Conclusions: No single contact angle technique fully characterizes the wettability of a material surface, and the application of complimenting methods allows increased characterization. At present, there is not an ISO standard method designed for soft materials. It is important that each contact angle technique has a standard protocol, as small protocol differences between laboratories often contribute to a variety of published data that are not easily comparable. © 2013 Contact Lens Association of Ophthalmologists.

Can surface energy measurements predict the impact of catalyst hydrophobicity upon fatty acid esterification over sulfonic acid functionalised periodic mesoporous organosilicas?

Sulfonic acid functionalised periodic mesoporous organosilicas (PrSO3 H-PMOs) with tunable hydrophobicity were synthesised via a surfactant-templating route, and characterised by porosimetry, TEM, XRD, XPS, inverse gas chromatography (IGC) and ammonia pulse chemisorption. IGC reveals that incorporation of ethyl or benzyl moieties into a mesoporous SBA-15 silica framework significantly increases the non-specific dispersive surface energy of adsorption for alkane adsorption, while decreasing the free energy of adsorption of methanol, reflecting increased surface hydrophobicity. The non-specific dispersive surface energy of adsorption of PMO-SO3H materials is strongly correlated with their activity towards palmitic acid esterification with methanol, demonstrating the power of IGC as an analytical tool for identifying promising solid acid catalysts for the esterification of free fatty acids. A new parameter [-ΔGCNP-P], defined as the per carbon difference in Gibbs free energy of adsorption between alkane and polar probe molecules, provides a simple predictor of surface hydrophobicity and corresponding catalyst activity in fatty acid esterification. © 2014 Elsevier B.V.

Anisotropic surface chemistry of aspirin crystals

The wettability of the (001), (100), and (011) crystallographic facets of macroscopic aspirin crystals has been experimentally investigated using a sessile drop contact angle (θ) method. θ for a nonpolar liquid was very similar for all three facets, though significant θ differences were observed for three polar probe liquids. The observed hydrophobicity of the (001) and (100) facets is ascribed to a reduced hydrogen bonding potential at these surfaces, whilst the observed hydrophilicity of facet (011) may be attributed to presence of surface carboxylic functionalities as confirmed by X-ray photoelectron spectroscopy (XPS). The dispersive component of the surface free energy (γ) was similar for all three facets (35 ± 2 mJ/m). The total surface energy, γs varied between 46 and 60 mJ/m due to significant variations in the polar/acid-base components of γ for all facets. Surface polarity as determined by γ measurements and XPS data were in good agreement, linking the variations in wettability to the concentration of oxygen containing surface functional groups. In conclusion, the wettability and the surface energy of a crystalline organic solid, such as aspirin, was found to be anisotropic and facet dependant, and in this case, related to the presence of surface carboxylic functionalities. © 2007 Wiley-Liss, Inc. and the American Pharmacists Association.

Isolation and characterization of the U1 and U2 snRNAs in the silk gland of the 5th instar larval stage silk moth Bombyx mori

In isolation and characterization studies, expression level U1 and U2 snRNA isoforms were obtained from the 5th instar larval stage silk gland (SG). The DNA content of the SG cells is approximately 200,000-fold higher compared to the usual (2N) somatic cells of B. mori due to endoreduplication. In this study, the existence of U1 and U2 snRNA isoforms in the SG of the organism is investigated. Bombyx mori U1 and U2-specific RT-PCR libraries from the silk gland were generated. Five U1 and eight U2 isoforms were isolated and characterized. Nucleotide differences, structural alterations, as well as protein and RNA interaction sites were analyzed in these variants. For the U1 snRNA variants, they were compared to the previously reported BmN isoforms. In all these U-snRNA variants, polymorphic sites do not predominate at the core of known functional sequences, which were interspecifically conserved. Variant sites and inter-species differences are located in moderately conserved regions. Free energy (ΔG) values for the entire U1 and U2 snRNA secondary structures and for the individual stem/loops domains of the isoforms were generated and compared to determine their structural stability. This will be the first time that U1 and U2 variants are shown specific for a development stage (larval) other than embryonic or adult. ^ Using phylogenetic analysis, evolutionary trees were generated for the U1 and U2 snRNAs using animal, plant, protista and fungal species. The resulting trees were boostrapped for robustness and rooted with the self-splicing RNA group II intron sequence from the cyanobacterium Calothrix. Using phylogenetic analyses, possible structural and functional evolutionary interdependence between the U1 and U2 snRNAs was investigated. ^

Electrostatic Contributions to the Thermodynamics of Ribonuclease P Protein Folding

Electrostatic interactions are of fundamental importance in determining the structure and stability of macromolecules. For example, charge-charge interactions modulate the folding and binding of proteins and influence protein solubility. Electrostatic interactions are highly variable and can be both favorable and unfavorable. The ability to quantify these interactions is challenging but vital to understanding the detailed balance and major roles that they have in different proteins and biological processes. Measuring pKa values of ionizable groups provides a sensitive method for experimentally probing the electrostatic properties of a protein.

pKa values report the free energy of site-specific proton binding and provide a direct means of studying protein folding and pH-dependent stability. Using a combination of NMR, circular dichroism, and fluorescence spectroscopy along with singular value decomposition, we investigated the contributions of electrostatic interactions to the thermodynamic stability and folding of the protein subunit of Bacillus subtilis ribonuclease P, P protein. Taken together, the results suggest that unfavorable electrostatics alone do not account for the fact that P protein is intrinsically unfolded in the absence of ligand because the pKa differences observed between the folded and unfolded state are small. Presumably, multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.

Topics in equilibrium and nonequilibrium thermodynamics: computing crystalline free energies and engineering Maxwell’s demon.

This dissertation covers two separate topics in statistical physics. The first part of the dissertation focuses on computational methods of obtaining the free energies (or partition functions) of crystalline solids. We describe a method to compute the Helmholtz free energy of a crystalline solid by direct evaluation of the partition function. In the many-dimensional conformation space of all possible arrangements of N particles inside a periodic box, the energy landscape consists of localized islands corresponding to different solid phases. Calculating the partition function for a specific phase involves integrating over the corresponding island. Introducing a natural order parameter that quantifies the net displacement of particles from lattices sites, we write the partition function in terms of a one-dimensional integral along the order parameter, and evaluate this integral using umbrella sampling. We validate the method by computing free energies of both face-centered cubic (FCC) and hexagonal close-packed (HCP) hard sphere crystals with a precision of $10^{-5}k_BT$ per particle. In developing the numerical method, we find several scaling properties of crystalline solids in the thermodynamic limit. Using these scaling properties, we derive an explicit asymptotic formula for the free energy per particle in the thermodynamic limit. In addition, we describe several changes of coordinates that can be used to separate internal degrees of freedom from external, translational degrees of freedom. The second part of the dissertation focuses on engineering idealized physical devices that work as Maxwell's demon. We describe two autonomous mechanical devices that extract energy from a single heat bath and convert it into work, while writing information onto memory registers. Additionally, both devices can operate as Landauer's eraser, namely they can erase information from a memory register, while energy is dissipated into the heat bath. The phase diagrams and the efficiencies of the two models are solved and analyzed. These two models provide concrete physical illustrations of the thermodynamic consequences of information processing.

Interactions of phenyldithioesters with gold nanoparticles (AuNPs): Implications for AuNP functionalization and molecular barcoding of AuNP assemblies

The interactions of phenyldithioesters with gold nanoparticles (AuNPs) have been studied by monitoring changes in the surface plasmon resonance (SPR), depolarised light scattering, and surface enhanced Raman spectroscopy (SERS). Changes in the SPR indicated that an AuNP-phenyldithioester charge transfer complex forms in equilibrium with free AuNPs and phenyldithioester. Analysis of the Langmuir binding isotherms indicated that the equilibrium adsorption constant, Kads, was 2.3 ± 0.1 × 106 M−1, which corresponded to a free energy of adsorption of 36 ± 1 kJ mol−1. These values are comparable to those reported for interactions of aryl thiols with gold and are of a similar order of magnitude to moderate hydrogen bonding interactions. This has significant implications in the application of phenyldithioesters for the functionalization of AuNPs. The SERS results indicated that the phenyldithioesters interact with AuNPs through the C═S bond, and the molecules do not disassociate upon adsorption to the AuNPs. The SERS spectra are dominated by the portions of the molecule that dominate the charge transfer complex with the AuNPs. The significance of this in relation to the use of phenyldithioesters for molecular barcoding of nanoparticle assemblies is discussed.