Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Solution structure of CnErg1 (Ergtoxin), a HERG specific scorpion toxin

The three-dimensional structure of chemically synthesized CnErg1 (Ergtoxin), which specifically blocks HERG (human ether-a-go-go-related gene) K+ channels, was determined by nuclear magnetic resonance spectroscopy. CnErg1 consists of a triple-stranded beta-sheet and an a-helix, as is typical of K+ channel scorpion toxins. The peptide structure differs from the canonical structures in that the first beta-strand is shorter and is nearer to the second beta-strand rather than to the third beta-strand on the C-terminus. There is also a large hydrophobic patch on the surface of the toxin, surrounding a central lysine residue, Lys13. We postulate that this hydrophobic patch is likely to form part of the binding surface of the toxin. (C) 2003 Published by Elsevier Science B.V. on behalf of the Federation of European Biochemical Societies.

A review of drainage and spontaneous rupture in free standing thin films with tangentially immobile interfaces

A review of spontaneous rupture in thin films with tangentially immobile interfaces is presented that emphasizes the theoretical developments of film drainage and corrugation growth through the linearization of lubrication theory in a cylindrical geometry. Spontaneous rupture occurs when corrugations from adjacent interfaces become unstable and grow to a critical thickness. A corrugated interface is composed of a number of waveforms and each waveform becomes unstable at a unique transition thickness. The onset of instability occurs at the maximum transition thickness, and it is shown that only upper and lower bounds of this thickness can be predicted from linear stability analysis. The upper bound is equivalent to the Freakel criterion and is obtained from the zeroth order approximation of the H-3 term in the evolution equation. This criterion is determined solely by the film radius, interfacial tension and Hamaker constant. The lower bound is obtained from the first order approximation of the H-3 term in the evolution equation and is dependent on the film thinning velocity A semi-empirical equation, referred to as the MTR equation, is obtained by combining the drainage theory of Manev et al. [J. Dispersion Sci. Technol., 18 (1997) 769] and the experimental measurements of Radoev et al. [J. Colloid Interface Sci. 95 (1983) 254] and is shown to provide accurate predictions of film thinning velocity near the critical thickness of rupture. The MTR equation permits the prediction of the lower bound of the maximum transition thickness based entirely on film radius, Plateau border radius, interfacial tension, temperature and Hamaker constant. The MTR equation extrapolates to Reynolds equation under conditions when the Plateau border pressure is small, which provides a lower bound for the maximum transition thickness that is equivalent to the criterion of Gumerman and Homsy [Chem. Eng. Commun. 2 (1975) 27]. The relative accuracy of either bound is thought to be dependent on the amplitude of the hydrodynamic corrugations, and a semiempirical correlation is also obtained that permits the amplitude to be calculated as a function of the upper and lower bound of the maximum transition thickness. The relationship between the evolving theoretical developments is demonstrated by three film thickness master curves, which reduce to simple analytical expressions under limiting conditions when the drainage pressure drop is controlled by either the Plateau border capillary pressure or the van der Waals disjoining pressure. The master curves simplify solution of the various theoretical predictions enormously over the entire range of the linear approximation. Finally, it is shown that when the Frenkel criterion is used to assess film stability, recent studies reach conclusions that are contrary to the relevance of spontaneous rupture as a cell-opening mechanism in foams. (C) 2003 Elsevier Science B.V. All rights reserved.

Cytotoxic iron chelators: Characterization of the structure, solution chemistry and redox activity of ligands and iron complexes of the di-2-pyridyl ketone isonicotinoyl hydrazone (HPKIH) analogues

Di-2-pyridyl ketone isonicotinoyl hydrazone (HPKIH) and a range of its analogues comprise a series of monobasic acids that are capable of binding iron (Fe) as tridentate (N,N,O) ligands. Recently, we have shown that these chelators are highly cytotoxic, but show selective activity against cancer cells. Particularly interesting was the fact that cytotoxicity of the HPKIH analogues is maintained even after complexation with Fe. To understand the potent anti-tumor activity of these compounds, we have fully characterized their chemical properties. This included examination of the solution chemistry and X-ray crystal structures of both the ligands and Fe complexes from this class and the ability of these complexes to mediate redox reactions. Potentiometric titrations demonstrated that all chelators are present predominantly in their charge-neutral form at physiological pH (7.4), allowing access across biological membranes. Keto-enol tautomerism of the ligands was identified, with the tautomers exhibiting distinctly different protonation constants. Interestingly, the chelators form low-spin (diamagnetic) divalent Fe complexes in solution. The chelators form distorted octahedral complexes with Fe-II, with two tridentate ligands arranged in a meridional fashion. Electrochemistry of the Fe complexes in both aqueous and non-aqueous solutions revealed that the complexes are oxidized to their ferric form at relatively high potentials, but this oxidation is coupled to a rapid reaction with water to form a hydrated (carbinolamine) derivative, leading to irreversible electrochemistry. The Fe complexes of the HPKIH analogues caused marked DNA degradation in the presence of hydrogen peroxide. This observation confirms that Fe complexes from the HPKIH series mediate Fenton chemistry and do not repel DNA. Collectively, studies on the solution chemistry and structure of these HPKIH analogues indicate that they can bind cellular Fe and enhance its redox activity, resulting in oxidative damage to vital biomolecules.