Resolving the energy paradox of chaperone/usher-mediated fibre assembly

Periplasmic chaperone/usher machineries are used for assembly of filamentous adhesion organelles of Gram-negative pathogens in a process that has been suggested to be driven by folding energy. Structures of mutant chaperone-subunit complexes revealed a final folding transition (condensation of the subunit hydrophobic core) on the release of organelle subunit from the chaperone-subunit pre-assembly complex and incorporation into the final fibre structure. However, in view of the large interface between chaperone and subunit in the pre-assembly complex and the reported stability of this complex, it is difficult to understand how final folding could release sufficient energy to drive assembly. In the present paper, we show the X-ray structure for a native chaperone-fibre complex that, together with thermodynamic data, shows that the final folding step is indeed an essential component of the assembly process. We show that completion of the hydrophobic core and incorporation into the fibre results in an exceptionally stable module, whereas the chaperone-subunit preassembly complex is greatly destabilized by the high-energy conformation of the bound subunit. This difference in stabilities creates a free energy potential that drives fibre formation.

Treatment of PVC using an alternative low energy ion bombardment procedure

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Article Accurate Experimental Values for the Free Energies of Hydration of H+, OH-, and H3O+

Accurate experimental values for the free energies of hydration, or the free energies of solvation, of the H+, OH-, and H3O+ ions are of fundamental importance. By use of the most accurate value for the free energy of solvation of H+, the known value for the free energy of solvation of water, and the known values for the gas phase and aqueous phase deprotonation of water, the corresponding experimental free energy of solvation for OH- is −106.4 ± 0.5 kcal/mol. Similarly, by use of the known values for ΔGf 0 for H3O, H2O+, and OH-, the known values for ΔGs for H+ and OH-, and the known value for the aqueous phase autoionization of water, we obtain an experimental free energy of solvation value for H3O+ of −103.4 ± 0.5 kcal/mol. These values are in excellent agreement with the commonly accepted values and with the value for ΔGs(OH-) obtained from embedding clusters of OH-(H2O)n in a dielectric continuum.

Global Search for Minimum Energy (H2O)n Clusters, n = 3−5

The Gaussian-3 (G3) model chemistry method has been used to calculate the relative ΔG° values for all possible conformers of neutral clusters of water, (H2O)n, where n = 3−5. A complete 12-fold conformational search around each hydrogen bond produced 144, 1728, and 20 736 initial starting structures of the water trimer, tetramer, and pentamer. These structures were optimized with PM3, followed by HF/6-31G* optimization, and then with the G3 model chemistry. Only two trimers are present on the G3 potential energy hypersurface. We identified 5 tetramers and 10 pentamers on the potential energy and free-energy hypersurfaces at 298 K. None of these 17 structures were linear; all linear starting models folded into cyclic or three-dimensional structures. The cyclic pentamer is the most stable isomer at 298 K. On the basis of this and previous studies, we expect the cyclic tetramers and pentamers to be the most significant cyclic water clusters in the atmosphere.

Hydration of the Bisulfate Ion: Atmospheric Implications

Using molecular dynamics configurational sampling combined with ab initio energy calculations, we determined the low energy isomers of the bisulfate hydrates. We calculated the CCSD(T) complete basis set (CBS) binding electronic and Gibbs free energies for 53 low energy isomers of HSO4–(H2O)n=1–6 and derived the thermodynamics of adding waters sequentially to the bisulfate ion and its hydrates. Comparing the HSO4–/H2O system to the neutral H2SO4/H2O cluster, water binds more strongly to the anion than it does to the neutral molecules. The difference in the binding thermodynamics of HSO4–/H2O and H2SO4/H2O systems decreases with increasing number of waters. The thermodynamics for the formation of HSO4–(H2O)n=1–5 is favorable at 298.15 K, and that of HSO4–(H2O)n=1–6 is favorable for T < 273.15 K. The HSO4– ion is almost always hydrated at temperatures and relative humidity values encountered in the troposphere. Because the bisulfate ion binds more strongly to sulfuric acid than it does to water, it is expected to play a role in ion-induced nucleation by forming a strong complex with sulfuric acid and water, thus facilitating the formation of a critical nucleus.

Use of intrinsic binding energy for catalysis by an RNA enzyme

The contribution of several individual ribozyme⋅substrate base pairs to binding and catalysis has been investigated using hammerhead ribozyme substrates that were truncated at their 3′ or 5′ ends. The base pairs at positions 1.1–2.1 and 15.2–16.2, which flank the conserved core, each contribute 104-fold in the chemical step, without affecting substrate binding. In contrast, base pairs distal to the core contribute to substrate binding but have no effect on the chemical step. These results suggest a “fraying model” in which each ribozyme⋅substrate helix can exist in either an unpaired (“open”) state or a helical (“closed”) state, with the closed state required for catalysis. The base pairs directly adjacent to the conserved core contribute to catalysis by allowing the closed state to form. Once the number of base pairs is sufficient to ensure that the closed helical state predominates, additional residues provide stabilization of the helix, and therefore increase binding, but have no further effect on the chemical step. Remarkably, the >5 kcal/mol free energy contribution to catalysis from each of the internal base pairs is considerably greater than the free energy expected for formation of a base pair. It is suggested that this unusually large energetic contribution arises because free energy that is typically lost in constraining residues within a base pair is expressed in the transition state, where it is used for positioning. This extends the concept of “intrinsic binding energy” from protein to RNA enzymes, suggesting that intrinsic binding energy is a fundamental feature of biological catalysis.

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.

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.

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.

Equilibrium fluid-solid coexistence of hard spheres

We present a tethered Monte Carlo simulation of the crystallization of hard spheres. Our method boosts the traditional umbrella sampling to the point of making practical the study of constrained Gibbs’ free energies depending on several crystalline order parameters. We obtain high-accuracy estimates of the fluid-crystal coexistence pressure for up to 2916 particles (enough to accommodate fluid-solid interfaces). We are able to extrapolate to infinite volume the coexistence pressure [p_(co) = 11.5727(10)k_(B)T/σ^(3)] and the interfacial free energy [γ_({100}) = 0.636(11)k_(B)T/σ^(2)].

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.

Investigating routes toward atomic layer deposition of silicon carbide: Ab initio screening of potential silicon and carbon precursors

Silicon carbide (SiC) is a promising material for electronics due to its hardness, and ability to carry high currents and high operating temperature. SiC films are currently deposited using chemical vapor deposition (CVD) at high temperatures 1500–1600 °C. However, there is a need to deposit SiC-based films on the surface of high aspect ratio features at low temperatures. One of the most precise thin film deposition techniques on high-aspect-ratio surfaces that operates at low temperatures is atomic layer deposition (ALD). However, there are currently no known methods for ALD of SiC. Herein, the authors present a first-principles thermodynamic analysis so as to screen different precursor combinations for SiC thin films. The authors do this by calculating the Gibbs energy ΔGΔG of the reaction using density functional theory and including the effects of pressure and temperature. This theoretical model was validated for existing chemical reactions in CVD of SiC at 1000 °C. The precursors disilane (Si2H6), silane (SiH4), or monochlorosilane (SiH3Cl) with ethyne (C2H2), carbontetrachloride (CCl4), or trichloromethane (CHCl3) were predicted to be the most promising for ALD of SiC at 400 °C.

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.

In situ imaging of catalytic etching on silver during methanol oxidation conditions by environmental scanning electron microscopy

Polycrystalline silver is used to catalytically oxidise methanol to formaldehyde. This paper reports the results of extensive investigations involving the use of environmental scanning electron microscopy (ESEM) to monitor structural changes in silver during simulated industrial reaction conditions. The interaction of oxygen, nitrogen, and water, either singly or in combination, with a silver catalyst at temperatures up to 973 K resulted in the appearance of a reconstructed silver surface. More spectacular was the effect an oxygen/methanol mixture had on the silver morphology. At a temperature of ca. 713 K pinholes were created in the vicinity of defects as a consequence of subsurface explosions. These holes gradually increased in size and large platelet features were created. Elevation of the catalyst temperature to 843 K facilitated the wholescale oxygen induced restructuring of the entire silver surface. Methanol reacted with subsurface oxygen to produce subsurface hydroxyl species which ultimately formed water in the subsurface layers of silver. The resultant hydrostatic pressure forced the silver surface to adopt a "hill and valley" conformation in order to minimise the surface free energy. Upon approaching typical industrial operating conditions widespread explosions occurred on the catalyst and it was also apparent that the silver surface was extremely mobile under the applied conditions. The interaction of methanol alone with silver resulted in the initial formation of pinholes primarily in the vicinity of defects, due to reaction with oxygen species incorporated in the catalyst during electrochemical synthesis. However, dramatic reduction in the hole concentration with time occurred as all the available oxygen became consumed. A remarkable correlation between formaldehyde production and hole concentration was found.