Membrane-derived fluorinated radicals detected by electron spin resonance in UV-irradiated Nafion and Dow ionomers: effect of counterions and H2O2

Electron spin resonance (ESR) spectroscopy was used to detect and identify radicals formed by UV irradiation of Nafion and Dow perfluorinated membranes partially or fully neutralized by Cu(II), Fe(II), and Fe(III). This method allowed the monitoring of ESR signals from the paramagnetic counterions together with the appearance of membrane-derived radical species. The most surprising aspect of this study was the formation of membrane-derived radical species only in the neutralized membranes, and even in the absence of H2O2 in the case of Nafion/Cu(II) and Nafion/Fe(III). In Nafion/Cu(II), ESR spectra from radicals exhibiting hyperfine interactions with three equivalent 19F nuclei (the "quartet") and with four equivalent 19F nuclei (the "quintet") were detected. In Nafion/Fe(II) exposed to H2O2 solutions, the formation of Fe(III) was detected. Upon UV irradiation, strong signals from the chain-end radical ROCF2CF2* were detected first, followed by the appearance, upon annealing above 200 K, of the quartet signal observed in Nafion/Cu(II). In subsequent experiments with Nafion and Dow membranes neutralized by Fe(III), the ROCF2CF2* radicals were formed even in the absence of H2O2, indicating that the role of H2O2 is oxidation of Fe(II) to Fe(III); moreover, in these systems small amounts of the chain-end radicals were detected even without UV irradiation. This result validates the method used to form the radicals: the role of UV irradiation is to accelerate the formation of a signal that is produced, albeit slowly, even in the dark, and possibly during fuel cell operation. The major conclusion is that cations are involved in degradation processes; the point of attack appears to be at or near the pendant chain of the ionomer. Therefore when studying membrane stability, it is important to consider not only the formation of oxygen radicals, such as HO*, HOO*, and O2*-, that can attack the membrane but also the specific reactivity of counterions.     

A recipe for evaluating and differentiating cos ϕ expressions

We present a simple recipe for calculating and differentiating cosine of bond angle and dihedral angle expressions. The resulting formulas can be incorporated in a straightforward manner into the bond angle and dihedral angle components of potential energy functions. These formulas rely only on expressions and derivatives of dot products, and, in particular, they avoid cross products as well as excessive Fortran function references. Consequently, the expressions derived in this article can be written compactly and evaluated rapidly.     

A powerful truncated Newton method for potential energy minimization

With advances in computer architecture and software, Newton methods are becoming not only feasible for large-scale nonlinear optimization problems, but also reliable, fast and efficient. Truncated Newton methods, in particular, are emerging as a versatile subclass. In this article we present a truncated Newton algorithm specifically developed for potential energy minimization. The method is globally convergent with local quadratic convergence. Its key ingredients are: (1) approximation of the Newton direction far away from local minima, (2) solution of the Newton equation iteratively by the linear Conjugate Gradient method, and (3) preconditioning of the Newton equation by the analytic second-derivative components of the “local” chemical interactions: bond length, bond angle and torsional potentials. Relaxation of the required accuracy of the Newton search direction diverts the minimization search away from regions where the function is nonconvex and towards physically interesting regions. The preconditioning strategy significantly accelerates the iterative solution for the Newton search direction, and therefore reduces the computation time for each iteration. With algorithmic variations, the truncated Newton method can be formulated so that storage and computational requirements are comparable to those of the nonlinear Conjugate Gradient method. As the convergence rate of nonlinear Conjugate Gradient methods is linear and performance less predictable, the application of the truncated Newton code to potential energy functions is promising.     

Inherent speedup limitations in multiple time step/particle mesh Ewald algorithms

Multiple time step (MTS) algorithms present an effective integration approach to reduce the computational cost of dynamics simulations. By using force splitting to allow larger time steps for the more slowly varying force components, computational savings can be realized. The Particle-Mesh-Ewald (PME) method has been independently devised to provide an effective and efficient treatment of the long-range electrostatics interactions. Here we examine the performance of a combined MTS/PME algorithm previously developed for AMBER on a large polymerase beta/DNA complex containing 40,673 atoms. Our goal is to carefully combine the robust features of the Langevin/MTS (LN) methodology implemented in CHARMM-which uses position rather than velocity Verlet with stochasticity to make possible outer time steps of 150 fs-with the PME formulation. The developed MTS/PME integrator removes fast terms from the reciprocal-space Ewald component by using switch functions. We analyze the advantages and limitations of the resulting scheme by comparing performance to the single time step leapfrog Verlet integrator currently used in AMBER by evaluating different time-step protocols using three assessors for accuracy, speedup, and stability, all applied to long (i.e., nanosecond) simulations to ensure proper energy conservation. We also examine the performance of the algorithm on a parallel, distributed shared-memory computer (SGI Origin 2000 with 8 300-MHz R12000 processors). Good energy conservation and stability behavior can be demonstrated, for Newtonian protocols with outer time steps of up to 8 fs and Langevin protocols with outer time steps of up to 16 fs. Still, we emphasize the inherent limitations imposed by the incorporation of MTS methods into the PME formulation that may not be widely appreciated. Namely, the limiting factor on the largest outer time-step size, and hence speedup, is an intramolecular cancellation error inherent to PME. This error stems from the excluded-nonbonded correction term contained in the reciprocal-space component. This cancellation error varies in time and introduces artificial frequencies to the governing dynamics motion. Unfortunately, we find that this numerical PME error cannot be easily eliminated by refining the PME parameters (grid resolution and/or order of interpolating polynomial). We suggest that methods other than PME for fast electrostatics may allow users to reap the full advantages from MTS algorithms.     

Constructing irregular surfaces to enclose macromolecular complexes for mesoscale modeling using the discrete surface charge optimization (DISCO) algorithm

Salt-mediated electrostatics interactions play an essential role in biomolecular structures and dynamics. Because macromolecular systems modeled at atomic resolution contain thousands of solute atoms, the electrostatic computations constitute an expensive part of the force and energy calculations. Implicit solvent models are one way to simplify the model and associated calculations, but they are generally used in combination with standard atomic models for the solute. To approximate electrostatics interactions in models on the polymer level (e.g., supercoiled DNA) that are simulated over long times (e.g., milliseconds) using Brownian dynamics, Beard and Schlick have developed the DiSCO (Discrete Surface Charge Optimization) algorithm. DiSCO represents a macromolecular complex by a few hundred discrete charges on a surface enclosing the system modeled by the Debye-Hückel (screened Coulombic) approximation to the Poisson-Boltzmann equation, and treats the salt solution as continuum solvation. DiSCO can represent the nucleosome core particle (>12,000 atoms), for example, by 353 discrete surface charges distributed on the surfaces of a large disk for the nucleosome core particle and a slender cylinder for the histone tail; the charges are optimized with respect to the Poisson-Boltzmann solution for the electric field, yielding a approximately 5.5% residual. Because regular surfaces enclosing macromolecules are not sufficiently general and may be suboptimal for certain systems, we develop a general method to construct irregular models tailored to the geometry of macromolecules. We also compare charge optimization based on both the electric field and electrostatic potential refinement. Results indicate that irregular surfaces can lead to a more accurate approximation (lower residuals), and the refinement in terms of the electric field is more robust. We also show that surface smoothing for irregular models is important, that the charge optimization (by the TNPACK minimizer) is efficient and does not depend on the initial assigned values, and that the residual is acceptable when the distance to the model surface is close to, or larger than, the Debye length. We illustrate applications of DiSCO's model-building procedure to chromatin folding and supercoiled DNA bound to Hin and Fis proteins. DiSCO is generally applicable to other interesting macromolecular systems for which mesoscale models are appropriate, to yield a resolution between the all-atom representative and the polymer level.     

Principal component analysis combined with truncated-Newton minimization for dimensionality reduction of chemical databases

 The similarity and diversity sampling problems are two challenging optimization tasks that arise in the analysis of chemical databases. As a first step to their solution, we propose an efficient projection/ refinement protocol based on the principal component analysis (PCA) and the truncated-Newton minimization method implemented by our package TNPACK (PCA/TNPACK). We show that PCA can provide the same initial guess as the singular value decomposition (SVD) for the optimization task of solving the distance-geometry optimization problem if each column of a database matrix has a mean of zero. Hence, our PCA/TNPACK approach is analogous to the SVD/TNPACK projection/refinement protocol that we developed recently for visualizing large chemical databases. Using PCA/TNPACK and the Merck MDDR database (MDL Drug Data Report), we further investigate the projection/refinement procedure with regards to the preservation of the original clusters of chemical compounds, the accuracy of similarity and diversity sampling of chemical compounds, and the potential application in the study of structure activity relationships. We also explore by simple experiments accuracy and efficiency aspects of the PCA/TNPACK procedure compared to those of a global optimization algorithm (simulated annealing, as implemented by the program package SIMANN) in terms of producing the projection mapping of a database. Numerical results show that the 2D PCA/TNPACK mapping can preserve the distance relationships of the original database and is thus valuable as a first step in similarity and diversity applications. Of course, the generation of a global rather than local minimizer and its interpretation in terms of pharameceutical applications remains a challenge. Since all numerical tests are performed on the Merck MDDR database, results are representative of realistic cases encountered in the field of drug design, and may help analyze properties of medicinal compounds. Received: December 21, 2000 / Accepted: August 19, 2002 Published online: September 27, 2002 Mathematics Subject Classification (2000): 65K10, 62H25, 92C50     

Molecular dynamics by the Backward-Euler method

This paper introduces a new computational method for molecular dynamics. The method combines the Backward-Euler scheme for the solution of stiff differential equations with a Langevin-equation approach to the establishment of thermal equilibrium. The method allows the user to choose a cutoff frequency ω_c. Vibrational modes with frequencies below ω_c will be fully excited (receive a mean energy of kT per mode), while modes with frequencies greater than ω_c will be effectively frozen by the method. By setting ω_c = kT/h, one can obtain reasonable agreement with the quantum-mechanical energy distribution among the various modes, despite the classical character of the computation.