Finite element modelling of dissipative structures for nonequilibrium chemical reactions in fluid-saturated porous media

We use the finite element method to model and predict the dissipative structures of chemical species for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. In particular, we explore the conditions under which dissipative structures of the species may exist in the Brusselator type of nonequilibrium chemical reaction. Since this is the first time the finite element method and related strategies have been used to study the chemical instability problems in a fluid-saturated porous medium, it is essential to validate the method and strategies before they are put into application. For this purpose, we have rigorously derived the analytical solutions for dissipative structures of chemical species in a benchmark problem, which geometrically is a square. Comparison of the numerical solutions with the analytical ones demonstrates that the proposed numerical method and strategy are robust enough to solve chemical instability problems in a fluid-saturated porous medium. Finally, the related numerical results from two application examples indicate that both the regime and the magnitude of pore-fluid flow have significant effects on the nature of the dissipative structures that developed for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. The motivation for this study is that self-organization under conditions of pore-fluid flow in a porous medium is a potential mechanism of the orebody formation and mineralization in the upper crust of the Earth. (C) 2000 Elsevier Science S.A. All rights reserved.

The use of NMR for determination of new structures in irradiated TFE/PMVE fluoropolymers

The radiation chemistry of two TFE/PMVE copolymers with TFE mole fractions of 0.66 and 0.81 has been studied under vacuum using Co-60 gamma -radiation over absorbed dose ranges up to 4.2 MGy. The radiolysis temperature was 313 K for both TFE/PMVE copolymers. New structure formation in the copolymers was identified by solid-state F-19 NMR and the G-values for new chain ends of 2.1 and 0.5 and for branching sites of 0.9 and 0.2 have been obtained for the TFE/PMVE with TFE mole fractions of 0.66 and 0.81, respectively. The relative yields of-O-CF3 and -CF2-CF3 chain ends were found to be proportional to the copolymer composition, but the yields of the -CF2-CF3 chain ends and -CF- branch points mere not linearly related ia the composition. rather they wets correlated with the radical yields measured at 77 K. (C) 2001 Elsevier Science Ltd. All rights reserved.

MR image-based measurement of rates of change in volumes of brain structures. Part 1: Method and validation

A detailed analysis procedure is described for evaluating rates of volumetric change in brain structures based on structural magnetic resonance (MR) images. In this procedure, a series of image processing tools have been employed to address the problems encountered in measuring rates of change based on structural MR images. These tools include an algorithm for intensity non-uniforniity correction, a robust algorithm for three-dimensional image registration with sub-voxel precision and an algorithm for brain tissue segmentation. However, a unique feature in the procedure is the use of a fractional volume model that has been developed to provide a quantitative measure for the partial volume effect. With this model, the fractional constituent tissue volumes are evaluated for voxels at the tissue boundary that manifest partial volume effect, thus allowing tissue boundaries be defined at a sub-voxel level and in an automated fashion. Validation studies are presented on key algorithms including segmentation and registration. An overall assessment of the method is provided through the evaluation of the rates of brain atrophy in a group of normal elderly subjects for which the rate of brain atrophy due to normal aging is predictably small. An application of the method is given in Part 11 where the rates of brain atrophy in various brain regions are studied in relation to normal aging and Alzheimer's disease. (C) 2002 Elsevier Science Inc. All rights reserved.

MR image-based measurement of rates of change in volumes of brain structures. Part II: Application to a study of Alzheimer's disease and normal aging

We present global and regional rates of brain atrophy measured on serially acquired T1-weighted brain MR images for a group of Alzheimer's disease (AD) patients and age-matched normal control (NC) subjects using the analysis procedure described in Part I. Three rates of brain atrophy: the rate of atrophy in the cerebrum, the rate of lateral ventricular enlargement and the rate of atrophy in the region of temporal lobes, were evaluated for 14 AD patients and 14 age-matched NC subjects. All three rates showed significant differences between the two groups, However, the greatest separation of the two groups was obtained when the regional rates were combined. This application has demonstrated that rates of brain atrophy, especially in specific regions of the brain, based on MR images can provide sensitive measures for evaluating the progression of AD. These measures will be useful for the evaluation of therapeutic effects of novel therapies for AD. (C) 2002 Elsevier Science Inc. All rights reserved.

Mimicking the action of GroEL in molecular dynamics simulations: Application to the refinement of protein structures

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results Sur port the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.

Solution structures by H-1 NMR of the novel cyclic trypsin inhibitor SFTI-1 from sunflower seeds and an acyclic permutant

SFTI-1 is a recently discovered cyclic peptide trypsin inhibitor from sunflower seeds comprising 14 amino acid residues. It is the most potent known Bowman-Birk inhibitor and the only naturally occurring cyclic one. The solution structure of SFTI-1 has been determined by H-1-NMR spectroscopy and compared with a synthetic acyclic permutant. The solution structures of both are remarkably similar. The lowest energy structures from each family of 20 structures of cyclic and acyclic SFTI-1 have an rmsd over the backbone and heavy atoms of 0.29 Angstrom and 0.66 Angstrom, respectively. The structures consist of two short antiparallel beta -strands joined by an extended loop containing the active site at one end. Cyclic SFTI-1 also has a hairpin turn completing the cycle. Both molecules contain particularly stable arrangements of cross-linking hydrogen bonds between the beta -strands and a single disulfide bridge, making them rigid and well defined in solution. These stable arrangements allow both the cyclic and acyclic variants of SFTI-1 to inhibit trypsin with very high potencies (0.5 nM and 12.1 nM, respectively). The cyclic nature of SFTI-1 appears to have evolved to provide higher trypsin inhibition as well as higher stability. The solution structures are similar to the crystal structure of the cyclic inhibitor in complex with trypsin. The lack of a major conformational change upon binding suggests that the structure of SFTI-1 is rigid and already pre-organized for maximal binding due to minimization of entropic losses compared to a more flexible ligand. These properties make SFTI-1 an ideal platform for the design of small peptidic pharmaceuticals or pesticides. (C) 2001 Academic Press.

Solution Structures of the cis- and trans-Pro30 Isomers of a Novel 38-Residue Toxin from the Venom of Hadronyche Infensa sp. that Contains a Cystine-Knot Motif within Its Four Disulfide Bonds

The primary sequence and three-dimensional structure of a novel peptide toxin isolated from the Australian funnel-web spider Hadronyche infensa sp. is reported. ACTX-HI:OB4219 contains 38 amino acids, including eight-cysteine residues that form four disulfide bonds. The connectivities of these disulfide bonds were previously unknown but have been unambiguously determined in this study. Three of these disulfide bonds are arranged in an inhibitor cystine-knot (ICK) motif, which is observed in a range of other disulfide-rich peptide toxins. The motif incorporates an embedded ring in the structure formed by two of the disulfides and their connecting backbone segments penetrated by a third disulfide bond. Using NMR spectroscopy, we determined that despite the isolation of a single native homologous product by RP-HPLC, ACTX-HI:OB4219 possesses two equally populated conformers in solution. These two conformers were determined to arise from cis/trans isomerization of the bond preceding Pro30. Full assignment of the NMR spectra for both conformers allowed for the calculation of their structures, revealing, the presence of a triple-stranded antiparallel sheet consistent with the inhibitor cystine-knot (ICK) motif.

Towards identifying the new structures formed on the γ-radiolysis of Ultem

Ultem irradiated up to 10.0 MGy has been analysed using C-13, H-1 and D-2 proton-carbon and proton-proton correlation NMR spectroscopy to shed light on the formation of new structures. Chemical shifts and correlation data were used to determine the structure or partial structures of several new components. The spectra indicated the presence of new groups and structures involving the isopropylidene group, the imide ring, and hydrogen-abstraction reactions. Possible pathways for formation of the new structures are proposed and the G-values for their formation have been estimated. (C) 2003 Elsevier Science Ltd. All rights reserved.

Potent cyclic antagonists of the complement C5a receptor on human polymorphonulcear leukocytes. Relationships between structures and activity

Human C5a is a plasma protein with potent chemoattractant and pro-inflammatory properties, and its overexpression correlates with severity of inflammatory diseases. C5a binds to its G protein-coupled receptor (C5aR) on polymorphonuclear leukocytes (PMNLs) through a high-affinity helical bundle and a low-affinity C terminus, the latter being solely responsible for receptor activation. Potent and selective C5a antagonists are predicted to be effective anti-inflammatory drugs, but no pharmacophore for small molecule antagonists has yet been developed, and it would significantly aid drug design. We have hypothesized that a turn conformation is important for activity of the C terminus of C5a and herein report small cyclic peptides that are stable turn mimics with potent antagonism at C5aR on human PMNLs. A comparison of solution structures for the C terminus of C5a, small acyclic peptide ligands, and cyclic antagonists supports the importance of a turn for receptor binding. Competition between a cyclic antagonist and either C5a or an acyclic agonist for C5aR on PMNLs supports a common or overlapping binding site on the C5aR. Structure-activity relationships for 60 cyclic analogs were evaluated by competitive radioligand binding with C5a (affinity) and myeloperoxidase release (antagonist potency) from human PMNLs, with 20 compounds having high antagonist potencies (IC50, 20 nM(-1) muM). Computer modeling comparisons reveal that potent antagonists share a common cyclic backbone shape, with affinity-determining side chains of defined volume projecting from the cyclic scaffold. These results define a new pharmacophore for C5a antagonist development and advance our understanding of ligand recognition and receptor activation of this G protein-coupled receptor.