Stretching single-domain proteins: Phase diagram and kinetics of force-induced unfolding

Single-molecule force spectroscopy reveals unfolding of domains in titin on stretching. We provide a theoretical framework for these experiments by computing the phase diagrams for force-induced unfolding of single-domain proteins using lattice models. The results show that two-state folders (at zero force) unravel cooperatively, whereas stretching of non-two-state folders occurs through intermediates. The stretching rates of individual molecules show great variations reflecting the heterogeneity of force-induced unfolding pathways. The approach to the stretched state occurs in a stepwise “quantized” manner. Unfolding dynamics and forces required to stretch proteins depend sensitively on topology. The unfolding rates increase exponentially with force f till an optimum value, which is determined by the barrier to unfolding when f = 0. A mapping of these results to proteins shows qualitative agreement with force-induced unfolding of Ig-like domains in titin. We show that single-molecule force spectroscopy can be used to map the folding free energy landscape of proteins in the absence of denaturants.

Bond orientational order, molecular motion, and free energy of high-density DNA mesophases.

By equilibrating condensed DNA arrays against reservoirs of known osmotic stress and examining them with several structural probes, it has been possible to achieve a detailed thermodynamic and structural characterization of the change between two distinct regions on the liquid-crystalline phase diagram: (i) a higher density hexagonally packed region with long-range bond orientational order in the plane perpendicular to the average molecular direction and (ii) a lower density cholesteric region with fluid-like positional order. X-ray scattering on highly ordered DNA arrays at high density and with the helical axis oriented parallel to the incoming beam showed a sixfold azimuthal modulation of the first-order diffraction peak that reflects the macroscopic bond-orientational order. Transition to the less-dense cholesteric phase through osmotically controlled swelling shows the loss of this bond orientational order, which had been expected from the change in optical birefringence patterns and which is consistent with a rapid onset of molecular positional disorder. This change in order was previously inferred from intermolecular force measurements and is now confirmed by 31P NMR. Controlled reversible swelling and compaction under osmotic stress, spanning a range of densities between approximately 120 mg/ml to approximately 600 mg/ml, allow measurement of the free-energy changes throughout each phase and at the phase transition, essential information for theories of liquid-crystalline states.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Free, long-chain, polyunsaturated fatty acids reduce membrane electrical excitability in neonatal rat cardiac myocytes.

Because previous studies showed that polyunsaturated fatty acids can reduce the contraction rate of spontaneously beating heart cells and have antiarrhythmic effects, we examined the effects of the fatty acids on the electrophysiology of the cardiac cycle in isolated neonatal rat cardiac myocytes. Exposure of cardiomyocytes to 10 microM eicosapentaenoic acid for 2-5 min markedly increased the strength of the depolarizing current required to elicit an action potential (from 18.0 +/- 2.4 pA to 26.8 +/- 2.7 pA, P < 0.01) and the cycle length of excitability (from 525 ms to 1225 ms, delta = 700 +/- 212, P < 0.05). These changes were due to an increase in the threshold for action potential (from -52 mV to -43 mV, delta = 9 +/- 3, P < 0.05) and a more negative resting membrane potential (from -52 mV to -57 mV, delta = 5 +/- 1, P < 0.05). There was a progressive prolongation of intervals between spontaneous action potentials and a slowed rate of phase 4 depolarization. Other polyunsaturated fatty acids--including docosahexaenoic acid, linolenic acid, linoleic acid, arachidonic acid, and its nonmetabolizable analog eicosatetraynoic acid, but neither the monounsaturated oleic acid nor the saturated stearic acid--had similar effects. The effects of the fatty acids could be reversed by washing with fatty acid-free bovine serum albumin. These results show that free polyunsaturated fatty acids can reduce membrane electrical excitability of heart cells and provide an electrophysiological basis for the antiarrhythmic effects of these fatty acids.