On uniquely determining cell-wall material properties with the compression experiment

A finite element model (FEM) of the cell-compression experiment has been developed in dimensionless form to extract the fundamental cell-wall-material properties (i.e. the constitutive equation and its parameters) from experiment force-displacement data. The FEM simulates the compression of a thin-walled, liquid-filled sphere between two flat surfaces. The cell-wall was taken to be permeable and the FEM therefore accounts for volume loss during compression. Previous models assume an impermeable wall and hence a conserved cell volume during compression. A parametric study was conducted for structural parameters representative of yeast. It was shown that the common approach of assuming reasonable values for unmeasured parameters (e.g. cell-wall thickness, initial radial stretch) can give rise to nonunique solutions for both the form and constants in the cell-wall constitutive relationship. Similarly, measurement errors can also lead to an incorrectly defined cell-wall constitutive relationship. Unique determination of the fundamental wall properties by cell compression requires accurate and precise measurement of a minimum set of parameters (initial cell radius, initial cell-wall thickness, and the volume loss during compression). In the absence of such measurements the derived constitutive relationship may be in considerable error, and should be evaluated against its ability to predict the outcome of other mechanical experiments. (C) 1998 Elsevier Science Ltd. All rights reserved.

Activation of beta(2)-adrenergic receptors hastens relaxation and mediates phosphorylation of phospholamban, troponin I, and C-protein in ventricular myocardium from patients with terminal heart failure

Background-Catecholamines hasten cardiac relaxation through beta-adrenergic receptors, presumably by phosphorylation of several proteins, but it is unknown which receptor subtypes are involved in human ventricle. We assessed the role of beta(1)- and beta(2)-adrenergic receptors in phosphorylating proteins implicated in ventricular relaxation. Methods and Results-Right ventricular trabeculae, obtained from freshly explanted hearts of patients with dilated cardiomyopathy (n=5) or ischemic cardiomyopathy (n=5), were paced at 60 bpm. After measurement of the contractile and relaxant effects of epinephrine (10 mu mol/L) or zinterol (10 mu mol/L), mediated through beta(2)-adrenergic receptors, and of norepinephrine (10 mu mol/L), mediated through beta(1)-adrenergic receptors, tissues were freeze clamped. We assessed phosphorylation of phospholamban, troponin I, and C-protein, as well as specific phosphorylation of phospholamban at serine 16 and threonine 17, Data did not differ between the 2 disease groups and were therefore pooled. Epinephrine, zinterol, and norepinephrine increased contractile force to approximately the same extent, hastened the onset of relaxation by 15+/-3%, 5+/-2%, and 20+/-3%, respectively, and reduced the time to half-relaxation by 26+/-3%, 21+/-3%, and 37+/-3%. These effects of epinephrine, zinterol, and norepinephrine were associated with phosphorylation (pmol phosphate/mg protein) of phospholamban 14+/-3, 12+/-4, and 12+/-3, troponin I 40+/-7, 33+/-7, and 31+/-6; and C-protein 7.2+/-1.9, 9.3 +/- 1.4, and 7.5 +/- 2.0. Phosphorylation of phospholamban occurred at both Ser16 and Thr17 residues through both beta(1)- and beta(2)-adrenergic receptors. Conclusions-Norepinephrine and epinephrine hasten human ventricular relaxation and promote phosphorylation of implicated proteins through both beta(1)- and beta(2)-adrenergic receptors, thereby potentially improving diastolic function.

Computational binding assays of antigenic peptides

Computer models can be combined with laboratory experiments for the efficient determination of (i) peptides that bind MHC molecules and (ii) T-cell epitopes. For maximum benefit, the use of computer models must be treated as experiments analogous to standard laboratory procedures. This requires the definition of standards and experimental protocols for model application. We describe the requirements for validation and assessment of computer models. The utility of combining accurate predictions with a limited number of laboratory experiments is illustrated by practical examples. These include the identification of T-cell epitopes from IDDM-, melanoma- and malaria-related antigens by combining computational and conventional laboratory assays. The success rate in determining antigenic peptides, each in the context of a specific HLA molecule, ranged from 27 to 71%, while the natural prevalence of MHC-binding peptides is 0.1-5%.

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.

The mechanical properties of Saccharomyces cerevisiae

Cell-wall mechanical properties play an integral part in the growth and form of Saccharomyces cerevisiae, In contrast to the tremendous knowledge on the genetics of S. cerevisiae, almost nothing is known about its mechanical properties. We have developed a micromanipulation technique to measure the force required to burst single cells and have recently established a mathematical model to extract the mechanical properties of the cell wall from such data, Here we determine the average surface modulus of the S, cerevisiae cell wall to be 11.1 +/- 0.6 N/m and 12.9 +/- 0.7 N/m in exponential and stationary phases, respectively, giving corresponding Young's moduli of 112 +/- 6 MPa and 107 +/- 6 MPa, This result demonstrates that yeast cell populations strengthen as they enter stationary phase by increasing wall thickness and hence the surface modulus, without altering the average elastic properties of the cell-wall material. We also determined the average breaking strain of the cell wall to be 82% +/- 3% in exponential phase and 80% +/- 3% in stationary phase, This finding provides a failure criterion that can be used to predict when applied stresses (e,g,, because of fluid flow) will lead to wall rupture, This work analyzes yeast compression experiments in different growth phases by using engineering methodology.

How thick is thick? Multicenter study of the rheological and material property characteristics of mealtime fluids and videofluoroscopy fluids

Objective rheological assessment of fluids given to dysphagic patients at mealtime and during videofluoroscopy was carried out using a multicenter format. Thin, quarter-thick, half-thick and full-thick fluids were examined for the degree of correlation between mealtime fluids and their allegedly matched videofluoroscopy counterparts. The study was carried out to determine whether perceived subjective differences between mealtime fluids and videofluoroscopy fluids could be quantified using the rheological parameters of viscosity, density, and yield stress. The results showed poor correlation between mealtime fluids and videofluoroscopy fluids over all parameters. In general, the videofluoroscopy fluids were more viscous, more dense, and showed higher yield stress values than their mealtime counterparts. Given these results, it is reasonable to assume that the fluids used during videofluoroscopy do not provide an accurate indication of swallowing ability at mealtime. Therefore, it is suggested that clinicians use objective methods to theologically match videofluoroscopy fluids to mealtime fluids.

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.

A computational analysis of substrate binding strength by phosphorylase kinase and protein kinase A

Protein kinases exhibit various degrees of substrate specificity. The large number of different protein kinases in the eukaryotic proteomes makes it impractical to determine the specificity of each enzyme experimentally. To test if it were possible to discriminate potential substrates from non-substrates by simple computational techniques, we analysed the binding enthalpies of modelled enzyme-substrate complexes and attempted to correlate it with experimental enzyme kinetics measurements. The crystal structures of phosphorylase kinase and cAMP-dependent protein kinase were used to generate models of the enzyme with a series of known peptide substrates and non-substrates, and the approximate enthalpy of binding assessed following energy minimization. We show that the computed enthalpies do not correlate closely with kinetic measurements, but the method can distinguish good substrates from weak substrates and non-substrates. Copyright (C) 2002 John Wiley Sons, Ltd.

Analytical solution to the zero-inertia problem for surge flow phenomena in nonprismatic channels

Surge flow phenomena. e.g.. as a consequence of a dam failure or a flash flood, represent free boundary problems. ne extending computational domain together with the discontinuities involved renders their numerical solution a cumbersome procedure. This contribution proposes an analytical solution to the problem, It is based on the slightly modified zero-inertia (ZI) differential equations for nonprismatic channels and uses exclusively physical parameters. Employing the concept of a momentum-representative cross section of the moving water body together with a specific relationship for describing the cross sectional geometry leads, after considerable mathematical calculus. to the analytical solution. The hydrodynamic analytical model is free of numerical troubles, easy to run, computationally efficient. and fully satisfies the law of volume conservation. In a first test series, the hydrodynamic analytical ZI model compares very favorably with a full hydrodynamic numerical model in respect to published results of surge flow simulations in different types of prismatic channels. In order to extend these considerations to natural rivers, the accuracy of the analytical model in describing an irregular cross section is investigated and tested successfully. A sensitivity and error analysis reveals the important impact of the hydraulic radius on the velocity of the surge, and this underlines the importance of an adequate description of the topography, The new approach is finally applied to simulate a surge propagating down the irregularly shaped Isar Valley in the Bavarian Alps after a hypothetical dam failure. The straightforward and fully stable computation of the flood hydrograph along the Isar Valley clearly reflects the impact of the strongly varying topographic characteristics on the How phenomenon. Apart from treating surge flow phenomena as a whole, the analytical solution also offers a rigorous alternative to both (a) the approximate Whitham solution, for generating initial values, and (b) the rough volume balance techniques used to model the wave tip in numerical surge flow computations.

Computational immunology: The coming of age

The explosive growth in biotechnology combined with major advancesin information technology has the potential to radically transformimmunology in the postgenomics era. Not only do we now have readyaccess to vast quantities of existing data, but new data with relevanceto immunology are being accumulated at an exponential rate. Resourcesfor computational immunology include biological databases and methodsfor data extraction, comparison, analysis and interpretation. Publiclyaccessible biological databases of relevance to immunologists numberin the hundreds and are growing daily. The ability to efficientlyextract and analyse information from these databases is vital forefficient immunology research. Most importantly, a new generationof computational immunology tools enables modelling of peptide transportby the transporter associated with antigen processing (TAP), modellingof antibody binding sites, identification of allergenic motifs andmodelling of T-cell receptor serial triggering.

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.

Effect of material anisotropy on the onset of convective flow in three-dimensional fluid-saturated faults

In order to investigate the effect of material anisotropy on convective instability of three-dimensional fluid-saturated faults, an exact analytical solution for the critical Rayleigh number of three-dimensional convective flow has been obtained. Using this critical Rayleigh number, effects of different permeability ratios and thermal conductivity ratios on convective instability of a vertically oriented three-dimensional fault have been examined in detail. It has been recognized that (1) if the fault material is isotropic in the horizontal direction, the horizontal to vertical permeability ratio has a significant effect on the critical Rayleigh number of the three-dimensional fault system, but the horizontal to vertical thermal conductivity ratio has little influence on the convective instability of the system, and (2) if the fault material is isotropic in the fault plane, the thermal conductivity ratio of the fault normal to plane has a considerable effect on the critical Rayleigh number of the three-dimensional fault system, but the effect of the permeability ratio of the fault normal to plane on the critical Rayleigh number of three-dimensional convective flow is negligible.

Conductivity, NMR and crystallographic study of N,N,N,N-tetramethylammonium dicyanamide plastic crystal phases: an archetypal ambient temperature plastic electrolyte material

N,N,N,N-Tetramethylammonium dicyanamide (Me(4)NDCA) has been examined via differential scanning calorimetry (DSC), thermogravimetric analysis, conductivity, single crystal X-ray diffraction and H-1 nuclear magnetic resonance (NMR) analyses, and was found to be highly conductive in the solid state (sigma = 10(-3) S cm(-2) at 420 K) and to also exhibit unusual plastic crystal behaviour. To investigate the correlation between such behaviour and the occurrence of molecular rotations in the crystal, H-1 NMR second moment measurements are compared with calculated values predicted from the crystal structure. While DSC analysis indicates a number of solid-solid transitions at ambient temperatures, subsequent H-1 NMR analysis of the Me4N+ cation shows that a variety of rotational motions become active at low (