The problem of the spreading of a liquid film along a solid surface: A new mathematical formulation

A new mathematical model is proposed for the spreading of a liquid film on a solid surface. The model is based on the standard lubrication approximation for gently sloping films (with the no-slip condition for the fluid at the solid surface) in the major part of the film where it is not too thin. In the remaining and relatively small regions near the contact lines it is assumed that the so-called autonomy principle holds—i.e., given the material components, the external conditions, and the velocity of the contact lines along the surface, the behavior of the fluid is identical for all films. The resulting mathematical model is formulated as a free boundary problem for the classical fourth-order equation for the film thickness. A class of self-similar solutions to this free boundary problem is considered.

Geometry of circular vectors and pattern recognition of shape of a boundary

This paper deals with pattern recognition of the shape of the boundary of closed figures on the basis of a circular sequence of measurements taken on the boundary at equal intervals of a suitably chosen argument with an arbitrary starting point. A distance measure between two boundaries is defined in such a way that it has zero value when the associated sequences of measurements coincide by shifting the starting point of one of the sequences. Such a distance measure, which is invariant to the starting point of the sequence of measurements, is used in identification or discrimination by the shape of the boundary of a closed figure. The mean shape of a given set of closed figures is defined, and tests of significance of differences in mean shape between populations are proposed.

The thermal explosion revisited

The classical problem of the thermal explosion in a long cylindrical vessel is modified so that only a fraction α of its wall is ideally thermally conducting while the remaining fraction 1−α is thermally isolated. Partial isolation of the wall naturally reduces the critical radius of the vessel. Most interesting is the case when the structure of the boundary is a periodic one, so that the alternating conductive α and isolated 1−α parts of the boundary occupy together the segments 2π/N (N is the number of segments) of the boundary. A numerical investigation is performed. It is shown that at small α and large N, the critical radius obeys a scaling law with the coefficients depending on N. For large N, the result is obtained that in the central core of the vessel the temperature distribution is axisymmetric. In the boundary layer near the wall having the thickness ≈2πr0/N (r0 is the radius of the vessel), the temperature distribution varies sharply in the peripheral direction. The temperature distribution in the axisymmetric core at the critical value of the vessel radius is subcritical.

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.

A gain-of-function of an amyotrophic lateral sclerosis-associated Cu,Zn-superoxide dismutase mutant: An enhancement of free radical formation due to a decrease in Km for hydrogen peroxide.

Cu,Zn-superoxide dismutase (SOD) is known to be a locus of mutation in familial amyotrophic lateral sclerosis (FALS). Transgenic mice that express a mutant Cu,Zn-SOD, Gly-93--> Ala (G93A), have been shown to develop amyotrophic lateral sclerosis (ALS) symptoms. We cloned the FALS mutant, G93A, and wild-type cDNA of human Cu,Zn-SOD, overexpressed them in Sf9 insect cells, purified the proteins, and studied their enzymic activities for catalyzing the dismutation of superoxide anions and the generation of free radicals with H2O2 as substrate. Our results showed that both enzymes contain one copper ion per subunit and have identical dismutation activity. However, the free radical-generating function of the G93A mutant, as measured by the spin trapping method, is enhanced relative to that of the wild-type enzyme, particularly at lower H2O2 concentrations. This is due to a small, but reproducible, decrease in the value of Km for H2O2 for the G93A mutant, while the kcat is identical for both enzymes. Thus, the ALS symptoms observed in G93A transgenic mice are not caused by the reduction of Cu,Zn-SOD activity with the mutant enzyme; rather, it is induced by a gain-of-function, an enhancement of the free radical-generating function. This is consistent with the x-ray crystallographic studies showing the active channel of the FALS mutant is slightly larger than that of the wild-type enzyme; thus, it is more accessible to H2O2. This gain-of-function, in part, may provide an explanation for the association between ALS and Cu,Zn-SOD mutants.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.