On the fully-nonlinear shallow-water generalized Serre equations

A fully-nonlinear weakly dispersive system for the shallow water wave regime is presented. In the simplest case the model was first derived by Serre in 1953 and rederived various times since then. Two additions to this system are considered: the effect of surface tension, and that of using a different reference fluid level to describe the velocity field. It is shown how the system can be further expanded by consistent exchanges of spatial and time derivatives. Properties of the solitary waves of the resulting system as well as a symmetric splitting of the equations based on the Riemann invariants of the hyperbolic shallow water system are presented. The latter leads to a fully-nonlinear one-way model and, upon further approximations, existing weakly nonlinear models. Our study also helps clarify the differences or similarities between existing models.     

An adaptive multi-grid chemistry (AMC) model for efficient simulation of HCCI and DI engine combustion

There is a need to reduce the computational expense of practical multidimensional combustion simulations. Simulation of Homogeneous Charge Compression Ignition (HCCI) engine processes requires consideration of detailed chemistry in order to capture the ignition and combustion characteristics. Even with relatively coarse numerical meshes and reduced chemistry mechanisms, calculation times are still unacceptably long. For the simulation of Direct Injection (DI) engines, fine meshes are needed to achieve the resolution required by the spray and mixing models, and they are computationally expensive even with reduced chemistry. In addition, the increasing application of CFD for engine design optimization is pushing the demand to reduce computational time. In current design optimizations, depending on the size of the parametric space, hundreds of individual simulations are needed.

This work presents an efficient Adaptive Multi-grid Chemistry (AMC) model that can be used in engine CFD codes for simulations of HCCI and DI engines with detailed chemistry. It was found that the number of cells computed with the chemistry solver can be reduced by two orders of magnitude for HCCI engines. The results predicted by the present KIVA AMC code are also consistent with those calculated by the original code using every cell.

In the method, progressively coarser grids are used for cells with similar gas properties in the chemistry calculation (up to four neighbour levels) or in the global method, cells are grouped without regard for their locations in the cylinder. Averaged and gradient-preserving remapping techniques used in multi-zone engine simulations were also explored. A parametric study was conducted for determining the model variables, such as the degree of local homogeneity for the multi-grid solvers.

The simulation results were compared with experimental data obtained from a Honda engine operated with n-heptane under HCCI conditions for which directly measured in-cylinder temperature and H₂O mole fraction data are available. In addition, simulation results were found to agree well with experimental data from a DI diesel engine operated under PCCI conditions with ultra-high EGR rates. It was found that computer time was reduced by a factor of ten for HCCI cases and two to three for DI cases without losing prediction accuracy.     

Visualizing abnormal climate changes in Central America from 1995 to 2000

This article describes statistical approaches that led to our main results regarding the abnormal climate changes of Central America from 1995 to 2000. These results are depicted in the poster () submitted at the contest of 2006 ASA Data Exposition.     

Interaction of nucleic acids with the glycocalyx

Mammalian cells resist the uptake of nucleic acids. The lipid bilayer of the plasma membrane presents one barrier. Here, we report on a second physicochemical barrier for uptake. To create a sensitive probe for nucleic acid-cell interactions, we synthesized fluorescent conjugates in which lipids are linked to DNA oligonucleotides. We found that these conjugates incorporate readily into the plasma membrane but are not retained there. Expulsion of lipid-oligonucleotide conjugates from the plasma membrane increases with oligonucleotide length. Conversely, the incorporation of conjugates increases markedly in cells that lack the major anionic components of the glycocalyx, sialic acid and glycosaminoglycans, and in cells that had incorporated highly cationic lipids into their plasma membrane. We conclude that anionic oligosaccharides provide a formidable barrier to the uptake of nucleic acids by mammalian cells. This conclusion has implications for genomic stability, as well as the delivery of genes and siRNAs into mammalian cells.     

Parallel β-sheet secondary structure is stabilized and terminated by interstrand disulfide cross-linking

Disulfide bonds between Cys residues in adjacent strands of parallel β-sheets are rare among proteins, which suggests that parallel β-sheet structure is not stabilized by such disulfide cross-links. We report experimental results that show, surprisingly, that an interstrand disulfide bond can stabilize parallel β-sheets formed by an autonomously folding peptide in aqueous solution. NMR analysis reveals that parallel β-sheet structure is terminated beyond the disulfide bond, which causes deviation from the extended backbone conformation at one of the Cys residues.     

Enantiomeric interactions between liquid crystals and organized monolayers of tyrosine-containing dipeptides

We have examined the orientational ordering of nematic liquid crystals (LCs) supported on organized monolayers of dipeptides with the goal of understanding how peptide-based interfaces encode intermolecular interactions that are amplified into supramolecular ordering. By characterizing the orientations of nematic LCs (4-cyano-4'-pentylbiphenyl and TL205 (a mixture of mesogens containing cyclohexane-fluorinated biphenyls and fluorinated terphenyls)) on monolayers of l-cysteine-l-tyrosine, l-cysteine-l-phenylalanine, or l-cysteine-l-phosphotyrosine formed on crystallographically textured films of gold, we conclude that patterns of hydrogen bonds generated by the organized monolayers of dipeptides are transduced via macroscopic orientational ordering of the LCs. This conclusion is supported by the observation that the ordering exhibited by the achiral LCs is specific to the enantiomers used to form the dipeptide-based monolayers. The dominant role of the -OH group of tyrosine in dictating the patterns of hydrogen bonds that orient the LCs was also evidenced by the effects of phosphorylation of the tyrosine on the ordering of the LCs. Overall, these results reveal that crystallographic texturing of gold films can direct the formation of monolayers of dipeptides with long-range order, thus unmasking the influence of hydrogen bonding, chirality, and phosphorylation on the macroscopic orientational ordering of LCs supported on these surfaces. These results suggest new approaches based on supramolecular assembly for reporting the chemical functionality and stereochemistry of synthetic and biological peptide-based molecules displayed at surfaces.     

A potent, versatile disulfide-reducing agent from aspartic acid

Dithiothreitol (DTT) is the standard reagent for reducing disulfide bonds between and within biological molecules. At neutral pH, however, >99% of DTT thiol groups are protonated and thus unreactive. Herein, we report on (2S)-2-amino-1,4-dimercaptobutane (dithiobutylamine or DTBA), a dithiol that can be synthesized from l-aspartic acid in a few high-yielding steps that are amenable to a large-scale process. DTBA has thiol pK(a) values that are ~1 unit lower than those of DTT and forms a disulfide with a similar E°' value. DTBA reduces disulfide bonds in both small molecules and proteins faster than does DTT. The amino group of DTBA enables its isolation by cation-exchange and facilitates its conjugation. These attributes indicate that DTBA is a superior reagent for reducing disulfide bonds in aqueous solution.     

Spectroscopic characterization of active-site variants of the PduO-type ATP:corrinoid adenosyltransferase from Lactobacillus reuteri: insights into the mechanism of four-coordinate Co(II)corrinoid formation

The PduO-type adenosine 5'-triphosphate (ATP):corrinoid adenosyltransferase from Lactobacillus reuteri (LrPduO) catalyzes the transfer of the adenosyl-group of ATP to Co(1+)cobalamin (Cbl) and Co(1+)cobinamide (Cbi) substrates to synthesize adenosylcobalamin (AdoCbl) and adenosylcobinamide (AdoCbi(+)), respectively. Previous studies revealed that to overcome the thermodynamically challenging Co(2+) → Co(1+) reduction, the enzyme drastically weakens the axial ligand-Co(2+) bond so as to generate effectively four-coordinate (4c) Co(2+)corrinoid species. To explore how LrPduO generates these unusual 4c species, we have used magnetic circular dichroism (MCD) and electron paramagnetic resonance (EPR) spectroscopic techniques. The effects of active-site amino acid substitutions on the relative yield of formation of 4c Co(2+)corrinoid species were examined by performing eight single-amino acid substitutions at seven residues that are involved in ATP-binding, an intersubunit salt bridge, and the hydrophobic region surrounding the bound corrin ring. A quantitative analysis of our MCD and EPR spectra indicates that the entire hydrophobic pocket below the corrin ring, and not just residue F112, is critical for the removal of the axial ligand from the cobalt center of the Co(2+)corrinoids. Our data also show that a higher level of coordination among several LrPduO amino acid residues is required to exclude the dimethylbenzimidazole moiety of Co(II)Cbl from the active site than to remove the water molecule from Co(II)Cbi(+). Thus, the hydrophilic interactions around and above the corrin ring are more critical to form 4c Co(II)Cbl than 4c Co(II)Cbi(+). Finally, when ATP analogues were used as cosubstrate, only "unactivated" five-coordinate (5c) Co(II)Cbl was observed, disclosing an unexpectedly large role of the ATP-induced active-site conformational changes with respect to the formation of 4c Co(II)Cbl. Collectively, our results indicate that the level of control exerted by LrPduO over the timing for the formation of the 4c Co(2+)corrinoid intermediates is even more exquisite than previously anticipated.     

The μ3 model of acids and bases: extending the Lewis theory to intermetallics

A central challenge in the design of new metallic materials is the elucidation of the chemical factors underlying the structures of intermetallic compounds. Analogies to molecular bonding phenomena, such as the Zintl concept, have proven very productive in approaching this goal. In this Article, we extend a foundational concept of molecular chemistry to intermetallics: the Lewis theory of acids and bases. The connection is developed through the method of moments, as applied to DFT-calibrated Hückel calculations. We begin by illustrating that the third and fourth moments (μ(3) and μ(4)) of the electronic density of states (DOS) distribution tune the properties of a pseudogap. μ(3) controls the balance of states above and below the DOS minimum, with μ(4) then determining the minimum's depth. In this way, μ(3) predicts an ideal occupancy for the DOS distribution. The μ(3)-ideal electron count is used to forge a link between the reactivity of transition metals toward intermetallic phase formation, and that of Lewis acids and bases toward adduct formation. This is accomplished through a moments-based definition of acidity which classifies systems that are electron-poor relative to the μ(3)-ideal as μ(3)-acidic, and those that are electron-rich as μ(3)-basic. The reaction of μ(3) acids and bases, whether in the formation of a Lewis acid/base adduct or an intermetallic phase, tends to neutralize the μ(3) acidity or basicity of the reactants. This μ(3)-neutralization is traced to the influence of electronegativity differences at heteroatomic contacts on the projected DOS curves of the atoms involved. The role of μ(3)-acid/base interactions in intermetallic phases is demonstrated through the examination of 23 binary phases forming between 3d metals, the stability range of the CsCl type, and structural trends within the Ti-Ni system.