Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below

In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton-Rogers-Lapwood problem in a fluid-saturated porous medium. The Horton-Rogers-Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton-Rogers-Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton-Rogers-Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton-Rogers-Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. (C) 1997 by John Wiley & Sons, Ltd.

Integrable eight-state supersymmetric U model with boundary terms and its Bethe ansatz solution

A class of integrable boundary terms for the eight-state supersymmetric U model are presented by solving the graded reflection equations. The boundary model is solved by using the coordinate Bethe ansatz method and the Bethe ansatz equations are obtained. (C) 1998 Elsevier Science B.V.

Analytical solution for the Mollow and Autler-Townes probe absorption spectra of a three-level atom in a squeezed vacuum

The Mellow and Autler-Townes probe absorption spectra of a three-level atom in a cascade configuration with the lower transition coherently driven and also coupled to a narrow bandwidth squeezed-vacuum field are studied. Analytical studies of the modifications caused by the finite squeezed-vacuum bandwidth to the spectra are made for the case when the Rabi frequency of the driving field is much larger than the natural linewidth. The squeezed vacuum center frequency and the driving laser frequency are assumed equal. We show that the spectral features depend on the bandwidth of a squeezed vacuum field and whether the sources of the squeezing field are degenerate (DPA) or nondegenerate (NDPA) parametric amplifiers. In a broadband or narrow bandwidth squeezed vacuum generated by a NDPA, the central component of the Mellow spectrum can be significantly narrower than that in the normal vacuum. When the source of the squeezed vacuum is a DPA, the central feature is insensitive to squeezing. The Rabi sidebands, however, can be significantly narrowed only in the squeezed vacuum produced by the DPA. The two lines of the Autler-Townes absorption spectrum can be narrowed only in a narrow bandwidth squeezed vacuum, whereas they are independent of the phase and are always broadened in a broadband squeezed vacuum.

Quantum integrability and exact solution of the supersymmetric U model with boundary terms

Quantum integrability is established for the one-dimensional supersymmetric U model with boundary terms by means of the quantum inverse-scattering method. The boundary supersymmetric U chain is solved by using the coordinate-space Bethe-ansatz technique and Bethe-ansatz equations are derived. This provides us with a basis for computing the finite-size corrections to the low-lying energies in the system. [S0163-1829(98)00425-1].

Solution of adsorption problems involving steep moving profiles

The moving finite element collocation method proposed by Kill et al. (1995) Chem. Engng Sci. 51 (4), 2793-2799 for solution of problems with steep gradients is further developed to solve transient problems arising in the field of adsorption. The technique is applied to a model of adsorption in solids with bidisperse pore structures. Numerical solutions were found to match the analytical solution when it exists (i.e. when the adsorption isotherm is linear). The method is simple yet sufficiently accurate for use in adsorption problems, where global collocation methods fail. (C) 1998 Elsevier Science Ltd. All rights reserved.

Solution structure of amyloid beta-peptide(1-40) in a water-micelle environment. Is the membrane-spanning domain where we think it is?

The three-dimensional solution structure of the 40 residue amyloid beta-peptide, A beta(1-40), has been determined using NMR spectroscopy at pH 5.1, in aqueous sodium dodecyl sulfate (SDS) micelles, In this environment, which simulates to some extent a water-membrane medium, the peptide is unstructured between residues 1 and 14 which are mainly polar and likely solvated by water. However, the rest of the protein adopts an alpha-helical conformation between residues 15 and 36 with a kink or hinge at 25-27. This largely hydrophobic region is likely solvated by SDS. Based on the derived structures, evidence is provided in support of a possible new location for the transmembrane domain of A beta within the amyloid precursor protein (APP). Studies between pH 4.2 and 7.9 reveal a pH-dependent helix-coil conformational switch. At the lower pH values, where the carboxylate residues are protonated, the helix is uncharged, intact, and lipid-soluble. As the pH increases above 6.0, part of the helical region (15-24) becomes less structured, particularly near residues E22 and D23 where deprotonation appears to facilitate unwinding of the helix. This pH-dependent unfolding to a random coil conformation precedes any tendency of this peptide to aggregate to a beta-sheet as the pH increases. The structural biology described herein for A beta(1-40) suggests that (i) the C-terminal two-thirds of the peptide is an alpha-helix in membrane-like environments, (ii) deprotonation of two acidic amino acids in the helix promotes a helix-coil conformational transition that precedes aggregation, (iii) a mobile hinge exists in the helical region of A beta(1-40) and this may be relevant to its membrane-inserting properties and conformational rearrangements, and (iv) the location of the transmembrane domain of amyloid precursor proteins may be different from that accepted in the Literature. These results may provide new insight to the structural properties of amyloid beta-peptides of relevance to Alzheimer's disease.

Small molecular probes for G-protein-coupled C5a receptors: Conformationally constrained antagonists derived from the C terminus of the human plasma protein C5a

Activation of the human complement system of plasma proteins in response to infection or injury produces a 4-helix bundle glycoprotein (74 amino acids) known as C5a. C5a binds to G-protein-coupled receptors on cell surfaces triggering receptor-ligand internalization, signal transduction, and powerful inflammatory responses. Since excessive levels of C5a are associated with autoimmune and chronic inflammatory disorders, inhibitors of receptor activation may have therapeutic potential. We now report solution structures and receptor-binding and antagonist activities for some of the first small molecule antagonists of C5a derived from its hexapeptide C terminus. The antagonist NMe-Phe-Lys-Pro-D-Cha-Trp-D-Arg-CO2H (1) surprisingly shows an unusually well-defined solution structure as determined by H-1 NMR spectroscopy. This is one of the smallest acyclic peptides found to possess a defined solution conformation, which can be explained by the constraining role of intramolecular hydrogen bonding. NOE and coupling constant data, slow deuterium exchange, and a low dependence on temperature for the chemical shift of the D-Cha-NH strongly indicate an inverse gamma turn stabilized by a D-Cha-NH ... OC-Lys hydrogen bond. Smaller conformational populations are associated with a hydrogen bond between Trp-NH ... OC-Lys, defining a type II beta turn distorted by the inverse gamma turn incorporated within it. An excellent correlation between receptor-affinity and antagonist activity is indicated for a limited set of synthetic peptides. Conversion of the C-terminal carboxylate of 1 to an amide decreases antagonist potency 5-fold, but potency is increased up to 10-fold over 1 if the amide bond is made between the C-terminal carboxylate and a Lys/Orn side chain to form a cyclic analogue. The solution structure of cycle 6 also shows gamma and beta turns; however, the latter occurs in a different position, and there are clear conformational changes in 6 vs 1 that result in enhanced activity. These results indicate that potent C5a antagonists can be developed by targeting site 2 alone of the C5a receptor and define a novel pharmacophore for developing powerful receptor probes or drug candidates.

Expokit: A software package for computing matrix exponentials

Expokit provides a set of routines aimed at computing matrix exponentials. More precisely, it computes either a small matrix exponential in full, the action of a large sparse matrix exponential on an operand vector, or the solution of a system of linear ODEs with constant inhomogeneity. The backbone of the sparse routines consists of matrix-free Krylov subspace projection methods (Arnoldi and Lanczos processes), and that is why the toolkit is capable of coping with sparse matrices of large dimension. The software handles real and complex matrices and provides specific routines for symmetric and Hermitian matrices. The computation of matrix exponentials is a numerical issue of critical importance in the area of Markov chains and furthermore, the computed solution is subject to probabilistic constraints. In addition to addressing general matrix exponentials, a distinct attention is assigned to the computation of transient states of Markov chains.

Solution structure of methionine-oxidized amyloid beta-peptide (1-40). Does oxidation affect conformational switching?

The solution structure of A beta(1-40)Met(O), the methionine-oxidized form of amyloid beta-peptide A beta(1-40), has been investigated by CD and NMR spectroscopy. Oxidation of Met35 may have implications in the aetiology of Alzheimer's disease. Circular dichroism experiments showed that whereas A beta(1-40) and A beta(1-40)Met(O) both adopt essentially random coil structures in water (pH 4) at micromolar concentrations, the former aggregates within several days while the latter is stable for at least 7 days under these conditions. This remarkable difference led us to determine the solution structure of A beta(1-40)Met(O) using H-1 NMR spectroscopy. In a water-SDS micelle medium needed to solubilize both peptides at the millimolar concentrations required to measure NMR spectra, chemical shift and NOE data for A beta(1-40)Met(O) strongly suggest the presence of a helical region between residues 16 and 24. This is supported by slow H-D exchange of amide protons in this region and by structure calculations using simulated annealing with the program XPLOR. The remainder of the structure is relatively disordered. Our previously reported NMR data for A beta(1-40) in the same solvent shows that helices are present over residues 15-24 (helix 1) and 28-36 (helix 2), Oxidation of Met35 thus causes a local and selective disruption of helix 2. In addition to this helix-coil rearrangement in aqueous micelles, the CD data show that oxidation inhibits a coil-to-beta-sheet transition in water. These significant structural rearrangements in the C-terminal region of A beta may be important clues to the chemistry and biology of A beta(1-40) and A beta(1-42).

The solution structure of a copper(II) compound of a new cyclic octapeptide by EPR spectroscopy and force field calculations

A new cyclic octapeptide, cyclo(Ile-Ser-(Gly)Thz-Ile-Thr-(Gly)Thz) (PatN), related to patellamide A, has been synthesized and reacted with copper(II) and base to form mono- and dinuclear complexes. The coordination environments around copper(TI) have been characterized by EPR spectroscopy. The solution structure of the thermodynamically most stable product, a purple dicopper(TI) compound, has been examined by simulating weakly dipole-dipole coupled EPR spectra based upon structural parameters obtained from force field (MM and MD) calculations. The MM-EPR method produces a saddle-shaped structure for [Cu-2(PatN)(OH2)(6)] that is similar to the known solution structure of patellamide A and the known solid-state structure of [Cu-2(AscidH(2))CO3(OH2)(2)]. Compared with the latter, [Cu-2(PatN)] has no carbonate bridge and a significantly flatter topology. The MM-EPR approach to solution-structure determination for paramagnetic metallopeptides may find wide applications to other metallopeptides and metalloproteins.

The solution structure of the N-terminal zinc finger of GATA-1 reveals a specific binding face for the transcriptional co-factor FOG

Zinc fingers (ZnFs) are generally regarded as DNA-binding motifs. However, a number of recent reports have implicated particular ZnFs in the mediation of protein-protein interactions. The N-terminal ZnF of GATA-1 (NF) is one such finger, having been shown to interact with a number of other proteins, including the recently discovered transcriptional co-factor FOG. Here we solve the three-dimensional structure of the NF in solution using multidimensional H-1/N-15 NMR spectroscopy, and we use H-1/N-15 spin relation measurements to investigate its backbone dynamics. The structure consists of two distorted beta-hairpins and a single alpha-helix, and is similar to that of the C-terminal ZnF of chicken GATA-1. Comparisons of the NF structure with those of other C-4-type zinc binding motifs, including hormone receptor and LIM domains, also reveal substantial structural homology. Finally, we use the structure to map the spatial locations of NF residues shown by mutagenesis to be essential for FOG binding, and demonstrate that these residues all lie on a single face of the NE Notably, this face is well removed from the putative DNA-binding face of the NE an observation which is suggestive of simultaneous roles for the NF; that is, stabilisation of GATA-1 DNA complexes and recruitment of FOG to GATA-1-controlled promoter regions.

Detection of lead(II) and silver(I) in aqueous solution using a piezoelectric device

The gold surface of a quartz crystal microbalance was modified by the attachment of silica particles derivatised with N-[(3-trimethoxysilyl)propyl] ethylenediaminetriacetic acid. The device was employed to study the kinetics of the interaction of aqueous solutions of lead(II) nitrate and silver(I) nitrate with the surface and for the selective separation of the metal ions.

Structure-activity relationships of omega-conotoxins MVIIA, MVIIC and 14 loop splice hybrids at N and P/Q-type calcium channels

The omega-conotoxins are a set of structurally related, four-loop, six cysteine containing peptides, that have a range of selectivities for different subtypes of the voltage-sensitive calcium channel (VSCC). To investigate the basis of the selectivity displayed by these peptides, we have studied the binding affinities of two naturally occurring omega-conotoxins, MVIIA and MVIIC and a series of 14 MVIIA/MVIIC loop hybrids using radioligand binding assays for N and P/Q-type Ca2+ channels in rat brain tissue. A selectivity profile was developed from the ratio of relative potencies at N-type VSCCs (using [I-125]GVIA radioligand binding assays) and P/Q-type VSCCs (using [I-125]MVIIC radioligand binding assays). in these peptides, loops 2 and 4 make the greatest contribution to VSCC subtype selectivity, while the effects of loops 1 and 3 are negligible. Peptides with homogenous combinations of loop 2 and 4 display clear selectivity preferences, while those with heterogeneous combinations of loops 2 and 4 are less discriminatory. H-1 NMR spectroscopy revealed that the global folds of MVIIA, MVIIC and the 14 loop hybrid peptides were similar; however, several differences in local structure were identified. Based on the binding data and the 3D structures of MVIIA, GVIA and MVIIC, we have developed a preliminary pharmacophore based on the omega-conotoxin residues most Likely to interact with the N-type VSCC. (C) 1999 Academic Press.

Solution structure of alpha-conotoxin ImI by H-1 nuclear magnetic resonance

alpha-Conotoxin ImI derives from the venom of Conus imperialis and is the first and only small-peptide ligand that selectively binds to the neuronal alpha(7) homopentameric subtype of the nicotinic acetylcholine receptor (nAChR). This receptor subtype is a possible drug target for several neurological disorders. The cysteines are connected in the pairs Cys2-Cys8 and Cys3-Cys12, To date it is the only alpha-conotoxin with a 4/3 residue spacing between the cysteines, The structure of ImI has been determined by H-1 NMR spectroscopy in aqueous solution, The NMR structure is of high quality, with a backbone pairwise rmsd of 0.34 Angstrom for a family of 19 structures, and comprises primarily a series of nested beta turns. Addition of organic solvent does not perturb the solution structure. The first eight residues of ImI are identical to the larger, but related, conotoxin EpI and adopt a similar structure, despite a truncated second loop. Residues important for binding of ImI to the alpha 7 nAChR are all clustered on one face of the molecule. Once further binding data for EPI and ImI are available, the ImI structure will allow for design of novel alpha(7) nAChR-specific agonists and antagonists with a wide range of potential pharmaceutical applications.