Computational chemistry on Fujitsu vector-parallel processors: Development and performance of applications software

In this and a preceding paper, we provide an introduction to the Fujitsu VPP range of vector-parallel supercomputers and to some of the computational chemistry software available for the VPP. Here, we consider the implementation and performance of seven popular chemistry application packages. The codes discussed range from classical molecular dynamics to semiempirical and ab initio quantum chemistry. All have evolved from sequential codes, and have typically been parallelised using a replicated data approach. As such they are well suited to the large-memory/fast-processor architecture of the VPP. For one code, CASTEP, a distributed-memory data-driven parallelisation scheme is presented. (C) 2000 Published by Elsevier Science B.V. All rights reserved.

Novel silica gel supported TiO2 photocatalyst synthesized by CVD method

TiO2 in anatase crystal phase is a very effective catalyst in the photocatalytic oxidation of organic compounds in water. To improve the recovery rate of TiO2 photocatalysts, which in most cases are in fine powder form, the chemical vapor deposition (CVD) method was used to load TiO2 onto a bigger particle support, silica gel. The amount of titania coating was found to depend strongly on the synthesis parameters of carrier gas flow rate and coating time. XPS and nitrogen ads/desorption results showed that most of the TiO2 particles generated from CVD were distributed on the external surface of the support and the coating was stable. The photocatalytic activities of TiO2/silica gel with different amounts of titania were evaluated for the oxidation of phenol aqueous solution and compared with that of Degussa P25. The optimum titania loading rate was found around 6 wt % of the TiO2 bulk concentration. Although the activity of the best TiO2/silica gel sample was still lower than that of P25, the synthesized TiO2/silica gel catalyst can be easily separated from the treated water and was found to maintain its TiO2 content and catalytic activity.

Conotoxin TVIIA, a novel peptide from the venom of Conus tulipa - 2. Three-dimensional solution structure

The three-dimensional solution structure of conotoxin TVIIA, a 30-residue polypeptide from the venom of the piscivorous cone snail Conus tulipa, has been determined using 2D H-1 NMR spectroscopy. TVIIA contains six cysteine residues which form a 'four-loop' structural framework common to many peptides from Conus venoms including the omega-, delta-, kappa-, and mu O-conotoxins. However, TVIIA does not belong to these well-characterized pharmacological classes of conotoxins, but displays high sequence identity with conotoxin GS, a muscle sodium channel blocker from Conus geographus. Structure calculations were based on 562 interproton distance restraints inferred from NOE data, together with 18 backbone and nine side-chain torsion angle restraints derived from spin-spin coupling constants. The final family of 20 structures had mean pairwise rms differences over residues 2-27 of 0.18 +/- 0.05 Angstrom for the backbone atoms and 1.39 +/- 0.33 Angstrom for all heavy atoms. The structure consists of a triple-stranded, antiparallel beta sheet with +2x, -1 topology (residues 7-9, 16-20 and 23-27) and several beta turns. The core of the molecule is formed by three disulfide bonds which form a cystine knot motif common to many toxic and inhibitory polypeptides. The global fold, molecular shape and distribution of amino-acid sidechains in TVIIA is similar to that previously reported for conotoxin GS, and comparison with other four-loop conotoxin structures provides further indication that TVIIA and GS represent a new and distinct subgroup of this structural family. The structure of TVIIA determined in this study provides the basis for determining a structure-activity relationship for these molecules and their interaction with target receptors.

Role of phosphorylation in the conformation of tau peptides implicated in Alzheimer's disease

A series of peptides corresponding to isolated regions of Tau (tau) protein have been synthesized and their conformations determined by H-1 NMR spectroscopy. Immunodominant peptides corresponding to tau(224-240) and a bisphosphorylated derivative in which a single Thr and a single Ser are phosphorylated at positions 231 and 235 respectively, and which are recognized by an Alzheimer's disease-specific monoclonal antibody, were the main focus of the study. The nonphosphorylated peptide adopts essentially a random coil conformation in aqueous solution, but becomes slightly more ordered into P-type structure as the hydrophobicity of the solvent is increased by adding up to 50% trifluoroethanol (TFE). Similar trends are observed for the bisphosphorylated peptide, with a somewhat stronger tendency to form an extended structure, There is tentative NMR evidence for a small population of species containing a turn at residues 229-231 in the phosphorylated peptide, and this is strongly supported by CD spectroscopy. A proposal that the selection of a bioactive conformation from a disordered solution ensemble may be an important step (in either tubulin binding or in the formation of PHF) is supported by kinetic data on Pro isomerization. A recent study showed that Thr231 phosphorylation affected the rate of prolyl isomerization and abolished tubulin binding. This binding was restored by the action of the prolyl isomerase Pin1. In the current study, we find evidence for the existence of both trans and cis forms of tau peptides in solution but no difference in the equilibrium distribution of cis-trans isomers upon phosphorylation. Increasing hydrophobicity decreases the prevalence of cis forms and increases the major trans conformation of each of the prolines present in these molecules. We also synthesized mutant peptides containing Tyr substitutions preceding the Pro residues and found that phosphorylation of Tyr appears to have an effect on the equilibrium ratio of cis-trans isomerization and decreases the cis content.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

'Tetol': a stereo-rigid four-strand motif for alkali and alkaline earth metal ion coordination

The tetraalcohol 2,3,5,6-endo,endo,endo,endo-tetrakis(hydroxymethyl]bicyclo[2.2.1]heptane (tetol, 1) has been prepared and crystallises readily as the lithium(I) complex [Li(1)(2)]Cl, forming an oligomeric multi-chain structure in which pairs of alcohols from two crystallographically independent tetol molecules bind lithium ions tetrahedrally. However, formation of monomeric structures in solution is inferred from electrospray mass spectroscopy, which has also shown evidence of exchange of lithium ion in the complexed species by added alkaline earth ions. (C) 2000 Elsevier Science S.A. All rights reserved.

Water wave-driven seepage in marine sediments

The action of water waves moving over a porous seabed drives a seepage flux into and out of the marine sediments. The volume of fluid exchange per wave cycle may affect the rate of contaminant transport in the sediments. In this paper, the dynamic response of the seabed to ocean waves is treated analytically on the basis of pore-elastic theory applied to a porous seabed. The seabed is modelled as a semi-infinite, isotropic, homogeneous material. Most previous investigations on the wave-seabed interaction problem have assumed quasi-static conditions within the seabed, although dynamic behaviour often occurs in natural environments. Furthermore, wave pressures used in the previous approaches were obtained from conventional ocean wave theories: which are based on the assumption of an impermeable rigid seabed. By introducing a complex wave number, we derive a new wave dispersion equation, which includes the seabed characteristics (such as soil permeability, shear modulus, etc.). Based on the new closed-form analytical solution, the relative differences of the wave-induced seabed response under dynamic and quasi-static conditions are examined. The effects of wave and soil parameters on the seepage flux per wave cycle are also discussed in detail. (C) 2000 Elsevier Science Ltd. All rights reserved.

Mechanisms of substitution reactions on cyclometallated platinum(IV) complexes: "Quasi-labile" systems

The substitution reactions of SMe2 by phosphines (PMePh2, PEtPh2, PPh3, P(4-MeC6H4)(3), P(3-MeC6H4)(3), PCy3) on Pt-IV complexes having a cyclometalated imine ligand, two methyl groups in a cis-geometrical arrangement, a halogen, and a dimethyl sulfide as ligands, [Pt(CN)(CH3)(2)(X)(SMe2)], have been studied as a function of temperature, solvent, and electronic and steric characteristics of the phosphines and the X and CN ligands. In all cases, a limiting dissociative mechanism has been found, where the dissociation of the SMe2 ligand corresponds to the rate-determining step. The pentacoordinated species formed behaves as a true pentacoordinated Pt-IV compound in a steady-state concentration, given the solvent independence of the rate constant. The X-ray crystal structures of two of the dimethyl sulfide complexes and a derivative of the pentacoordinate intermediate have been determined. Differences in the individual rate constants for the entrance of the phosphine ligand can only be estimated as reactivity ratios. In all cases an effect of the phosphine size is detected, indicating that an associative step takes place from the pentacoordinated intermediate. The nature of the (CN) imine and X ligands produces differences in the dimethyl sulfide dissociation reactions rates, which can be quantified by the corresponding DeltaS double dagger values (72, 64, 48, 31, and 78 J K-1 mol(-1) for CN/X being C6H4CHNCH2C6H5/Br, C6H4CHNCH2-(2,4,6-(CH3)(3))C6H2/Br, C6H4CHNCH2C6H5/Cl, C6Cl4CHNCH2C6H5/Cl, and C6W4CH2NCHC6H5/ Pr, respectively). As a whole, the donor character of the coordinated C-aromatic and X atoms have the greatest influence on the dissociativeness of the rate-determining step.

Characterization of the structural and surface properties of chemically modified MCM-41 material

Mesoporous Mobil catalytic materials of number 41 (MCM-41) silica was chemically modified using both inorganic and organic precursors and characterized using the techniques, XRD, XPS, MAS NMR, FTIR, W-Vis, and physical adsorption of nitrogen, hydrocarbons (hexane, benzene, acetone, and methanol) and water vapor. Modification using organic reagents was found to result in a significant loss in porosity and a shape change of surface properties (increased hydrophobicity and decreased acidity). With inorganic modifying reagents, the decrease in porosity was also observed while the surface properties were not significantly altered as reflected by the adsorption isotherms of organics and water vapors. Chemical modifications can greatly improve the hydrothermal stability of MCM-41 material because of the enhanced surface hydrophobicity (with organic modifiers) or increased pore wall thickness (with inorganic modifiers). (C) 2000 Elsevier Science B.V. All rights reserved.

Encoding combinatorial libraries: A novel application of fluorescent silica colloids

A major challenge associated with using large chemical libraries synthesized on microscopic solid support beads is the rapid discrimination of individual compounds in these libraries. This challenge can be overcome by encoding the beads with 1 mum silica colloidal particles (reporters) that contain specific and identifiable combinations of fluorescent byes. The colored bar code generated on support beads during combinatorial library synthesis can be easily, rapidly, and inexpensively decoded through the use of fluorescence microscopy. All reporters are precoated with polyelectrolytes [poly(acrylic acid), PAA, poly(sodium 4-styrenesulfonate PSSS, polyethylenimine, PEI, and/or poly(diallyldimethylammonium chloride), PDADMAC] with the aim of enhancing surface charge, promoting electrostatic attraction to the bead, and facilitating polymer bridging between the bead and reporter for permanent adhesion. As shown in this article, reporters coated with polyelectrolytes clearly outperform uncoated reporters with regard to quantity of attached reporters per bead (54 +/- 23 in 2500 mum(2) area for PEI/PAA coated and 11 +/- 6 for uncoated reporters) and minimization of cross-contamination (1 red reporter in 2500 mum(2) area of green-labeled bead for PEI/PAA coated and 26 +/- 15 red reporters on green-labeled beads for uncoated reporters after 10 days). Examination of various polyelectrolyte systems shows that the magnitude of the xi -potential of polyelectrolyte-coated reporters (-64 mV for PDADMAC/PSSS and -42 mV for PEI/PAA-coated reporters) has no correlation with the number of reporters that adhere to the solid support beads (21 +/- 16 in 2500 mum(2) area for PDADMAC/PSSS and 54 +/- 23 for PEI/PAA-coated reporters). The contribution of polymer bridging to the adhesion has a far greater influence than electrostatic attraction and is demonstrated by modification of the polyelectrolyte multilayers using gamma irradiation of precoated reporters either in aqueous solution or in polyelectrolyte solution.

Commentary: Using the convection-dispersion model and transit time density functions in the analysis of organ distribution kinetics

The convection-dispersion model and its extended form have been used to describe solute disposition in organs and to predict hepatic availabilities. A range of empirical transit-time density functions has also been used for a similar purpose. The use of the dispersion model with mixed boundary conditions and transit-time density functions has been queried recently by Hisaka and Sugiyanaa in this journal. We suggest that, consistent with soil science and chemical engineering literature, the mixed boundary conditions are appropriate providing concentrations are defined in terms of flux to ensure continuity at the boundaries and mass balance. It is suggested that the use of the inverse Gaussian or other functions as empirical transit-time densities is independent of any boundary condition consideration. The mixed boundary condition solutions of the convection-dispersion model are the easiest to use when linear kinetics applies. In contrast, the closed conditions are easier to apply in a numerical analysis of nonlinear disposition of solutes in organs. We therefore argue that the use of hepatic elimination models should be based on pragmatic considerations, giving emphasis to using the simplest or easiest solution that will give a sufficiently accurate prediction of hepatic pharmacokinetics for a particular application. (C) 2000 Wiley-Liss Inc. and the American Pharmaceutical Association J Pharm Sci 89:1579-1586, 2000.

Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

Adsorption hysteresis and criticality in regular mesoporous materials

In the present work, various theories predicting the critical diameter for the absence of capillary condensation and hysteresis are applied to experimental adsorption isotherms of vapors on regular mesoporous materials. Among the various theories studied, the tensile strength approximation proposed by the authors was found to be the most successful. Reversibility of nitrogen adsorption at 77.4 K was studied on pure MCM-41 of various pore sizes, as well as mixtures of pure MCM-41 samples in a 1:1 ratio. The results of PSD and hysteresis on MCM-41 mixtures are close to that expected from studies of the pure materials. The estimates of hysteresis critical temperature and diameter of MCM-41, HMS, FSM and KIT materials are also provided.