Efficient accommodation of local minima in watershed model calibration

The Gauss-Marquardt-Levenberg (GML) method of computer-based parameter estimation, in common with other gradient-based approaches, suffers from the drawback that it may become trapped in local objective function minima, and thus report optimized parameter values that are not, in fact, optimized at all. This can seriously degrade its utility in the calibration of watershed models where local optima abound. Nevertheless, the method also has advantages, chief among these being its model-run efficiency, and its ability to report useful information on parameter sensitivities and covariances as a by-product of its use. It is also easily adapted to maintain this efficiency in the face of potential numerical problems (that adversely affect all parameter estimation methodologies) caused by parameter insensitivity and/or parameter correlation. The present paper presents two algorithmic enhancements to the GML method that retain its strengths, but which overcome its weaknesses in the face of local optima. Using the first of these methods an intelligent search for better parameter sets is conducted in parameter subspaces of decreasing dimensionality when progress of the parameter estimation process is slowed either by numerical instability incurred through problem ill-posedness, or when a local objective function minimum is encountered. The second methodology minimizes the chance of successive GML parameter estimation runs finding the same objective function minimum by starting successive runs at points that are maximally removed from previous parameter trajectories. As well as enhancing the ability of a GML-based method to find the global objective function minimum, the latter technique can also be used to find the locations of many non-global optima (should they exist) in parameter space. This can provide a useful means of inquiring into the well-posedness of a parameter estimation problem, and for detecting the presence of bimodal parameter and predictive probability distributions. The new methodologies are demonstrated by calibrating a Hydrological Simulation Program-FORTRAN (HSPF) model against a time series of daily flows. Comparison with the SCE-UA method in this calibration context demonstrates a high level of comparative model run efficiency for the new method. (c) 2006 Elsevier B.V. All rights reserved.

The cost of uniqueness in groundwater model calibration

Calibration of a groundwater model requires that hydraulic properties be estimated throughout a model domain. This generally constitutes an underdetermined inverse problem, for which a Solution can only be found when some kind of regularization device is included in the inversion process. Inclusion of regularization in the calibration process can be implicit, for example through the use of zones of constant parameter value, or explicit, for example through solution of a constrained minimization problem in which parameters are made to respect preferred values, or preferred relationships, to the degree necessary for a unique solution to be obtained. The cost of uniqueness is this: no matter which regularization methodology is employed, the inevitable consequence of its use is a loss of detail in the calibrated field. This, ill turn, can lead to erroneous predictions made by a model that is ostensibly well calibrated. Information made available as a by-product of the regularized inversion process allows the reasons for this loss of detail to be better understood. In particular, it is easily demonstrated that the estimated value for an hydraulic property at any point within a model domain is, in fact, a weighted average of the true hydraulic property over a much larger area. This averaging process causes loss of resolution in the estimated field. Where hydraulic conductivity is the hydraulic property being estimated, high averaging weights exist in areas that are strategically disposed with respect to measurement wells, while other areas may contribute very little to the estimated hydraulic conductivity at any point within the model domain, this possibly making the detection of hydraulic conductivity anomalies in these latter areas almost impossible. A study of the post-calibration parameter field covariance matrix allows further insights into the loss of system detail incurred through the calibration process to be gained. A comparison of pre- and post-calibration parameter covariance matrices shows that the latter often possess a much smaller spectral bandwidth than the former. It is also demonstrated that, as all inevitable consequence of the fact that a calibrated model cannot replicate every detail of the true system, model-to-measurement residuals can show a high degree of spatial correlation, a fact which must be taken into account when assessing these residuals either qualitatively, or quantitatively in the exploration of model predictive uncertainty. These principles are demonstrated using a synthetic case in which spatial parameter definition is based oil pilot points, and calibration is Implemented using both zones of piecewise constancy and constrained minimization regularization. (C) 2005 Elsevier Ltd. All rights reserved.

Optimum coal properties and engineering processes for CO2 geosequestration in Coal

CO2 Geosequestration is seen by many worldwide scientists and engineers as a leading prospective solution to the global warming problem arising from excessive CO2 and other greenhouse gas emissions. CO2 geosequestration in coal seams has two important strategic benefits: the process has an extremely low risk of leakage, due to the adsorbed state of the CO2 and the known reservoir context of essentially-zero leakage into which it is be injected; the second benefit arises from the valuable by-product, clean burning coalbed methane gas. This paper presents the authors’ experience, knowledge and perspective on what coal properties and engineering processes would favour implementing a demonstration or commercial CO2 storage-in-coal project, in Queensland, Australia. As such, it may be considered a template for screening studies to select the optimum coal seam reservoir, and for preliminary studies in designing the injection system and predicting production response to the technology. The paper concludes by examining the current knowledge gaps of CO2 geosequestration in coal, identifying further basic and applied research topics.

Quantum symmetries and the Weyl–Wigner product of group representations

In the usual formulation of quantum mechanics, groups of automorphisms of quantum states have ray representations by unitary and antiunitary operators on complex Hilbert space, in accordance with Wigner's theorem. In the phase-space formulation, they have real, true unitary representations in the space of square-integrable functions on phase space. Each such phase-space representation is a Weyl–Wigner product of the corresponding Hilbert space representation with its contragredient, and these can be recovered by 'factorizing' the Weyl–Wigner product. However, not every real, unitary representation on phase space corresponds to a group of automorphisms, so not every such representation is in the form of a Weyl–Wigner product and can be factorized. The conditions under which this is possible are examined. Examples are presented.

Pilot-scale extraction of PHB from recombinant E-coli by homogenization and centrifugation

A new method of poly-beta-hydroxybutyrate (PHB) extraction from recombinant E. coli is proposed, using homogenization and centrifugation coupled with sodium hypochlorite treatment. The size of PHB granules and cell debris in homogenates was characterised as a function of the number of homogenization passes. Simulation was used to develop the PHB and cell debris fractionation system, enabling numerical examination of the effects of repeated homogenization and centrifuge-feedrate variation. The simulation provided a good prediction of experimental performance. Sodium hypochlorite treatment was necessary to optimise PHB fractionation. A PHB recovery of 80% at a purity of 96.5% was obtained with the final optimised process. Protein and DNA contained in the resultant product were negligible. The developed process holds promise for significantly reducing the recovery cost associated with PHB manufacture.

Murine DEP-1, a receptor protein tyrosine phosphatase, is expressed in macrophages and is regulated by CSF-1 and LPS

The spectrum of protein tyrosine phosphatases (PTPs) expressed in bone marrow-derived murine macrophages (BMMs) was examined using reverse transcriptase-polymerase chain reaction. Ten different PTP cDNAs were isolated and in this study we focus on mDEP-1, a type III receptor PTP. Three mDEP-1 transcripts were expressed in primary macrophages and macrophage cell lines and were induced during macrophage differentiation of M1 myeloid leukemia cells. A valiant mRNA Tvas identified that encodes an alternate carboxyl-terminus and 3' UTR. The expression of mDEP-1 was down-regulated by CSF-1 (macrophage colony-stimulating factor) and up-regulated by bacterial lipopolysaccharide, an important physiological regulator of macrophage function that opposes CSF-1 action. Whole mount irt situ hybridization, and immunolocalization of the protein, confirmed that mDEP-1 is expressed by a subset of embryonic macrophages in the liver and mesenchyme. mDEP-1 was also detected in the eye and peripheral nervous system of the developing embryo. Attempts to express mDEP-1 constitutively in the macrophage cell line RAW264 were unsuccessful, with results suggesting that the gene product inhibits cell proliferation.

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.

Molecular recognition of macrocyclic peptidomimetic inhibitors by HIV-1 protease

High-resolution crystal structures are described for seven macrocycles complexed with HIV-1 protease (HIVPR). The macrocycles possess two amides and an aromatic group within 15-17 membered rings designed to replace N- or C-terminal tripeptides from peptidic inhibitors of HIVPR. Appended to each macrocycle is a transition state isostere and either an acyclic peptide, nonpeptide, or another macrocycle. These cyclic analogues are potent inhibitors of HIVPR, and the crystal structures show them to be structural mimics of acyclic peptides, binding in the active site of HIVPR via the same interactions. Each macrocycle is restrained to adopt a P-strand conformation which is preorganized for protease binding. An unusual feature of the binding of C-terminal macrocyclic inhibitors is the interaction between a positively charged secondary amine and a catalytic aspartate of HIVPR. A bicyclic inhibitor binds similarly through its secondary amine that lies between its component N-terminal and C-terminal macrocycles. In contrast, the corresponding tertiary amine of the N-terminal macrocycles does not interact with the catalytic aspartates. The amine-aspartate interaction induces a 1.5 Angstrom N-terminal translation of the inhibitors in the active site and is accompanied by weakened interactions with a water molecule that bridges the ligand to the enzyme, as well as static disorder in enzyme flap residues. This flexibility may facilitate peptide cleavage and product dissociation during catalysis. Proteases [Aba(67,95)]HIVPR and [Lys(7),Ile(33),Aba(67,95)]- HIVPR used in this work were shown to have very similar crystal structures.

Oxidation of oxa and thia fatty acids and related compounds catalysed by 5-and 15-lipoxygenase

The modified fatty acids, (Z,Z,Z)-(octadeca-6,9,12-trienyloxy)acetic acid, (Z,Z,Z)-(octadeca-9,12,15-trienyloxy)acetic acid, (all-Z)-(eicosa-5,8,11,14-tetraenyloxy)acetic acid, (all-Z)-(eicosa-5,8,11,14-tetraenylthio)acetic acid, 3-[(all-Z)-(eicosa-5,8,11,14-tetraenylthio)]propionic acid, (all-Z)-(eicosa-5,8,11,14-tetraenylthio)succinic acid, N-[(all-Z)-(eicosa-5,8,11,14-tetraenoyl)]glycine and N-[(all-Z)-(eicosa-5,8,11,14-tetraenoyl)]aspartic acid, all react with soybean 15-lipoxygenase. The products were treated with triphenylphosphine to give alcohols, which were isolated using HPLC. Analysis of the alcohols using negative ion tandem electrospray mass spectrometry, and by comparison with compounds obtained by autoxidation of arachidonic acid, shows that each enzyme catalysed oxidation occurs at the omega -6 position of the substrate. In a similar fashion, it has been found that (Z,Z,Z)-(octadeca-6,9,12-trienyloxy)acetic acid, (Z,Z,Z)-(octadeca-9,12,15-trienyloxy)acetic acid, (all-Z)-(eicosa-5,8,11,14-tetraenylthio)acetic acid and N-[(all-Z)-(eicosa-5,8, 11.14-tetraenylthio)]propionic acid each undergoes regioselective oxidation at the carboxyl end of the polyene moiety on treatment with potato 5-lipoxygenase. Neither (all-Z)-(eicosa-5,8,11,14-tetraenylthio)succinic acid nor N-[(all-Z)-(eicosa-5,8,11,14-tetraenoyl)]aspartic acid reacts in the presence of this enzyme, while N-[(all-Z)-(eicosa-5,8,11,14-tetraenoyl)]glycine affords the C11' oxidation product. The alcohol derived from (Z,Z,Z)-(octadeca-6,9, 12-trienyloxy)acetic acid using the 15-lipoxygenase reacts at the C6' position with the 5-lipoxygenase. (C) 2001 Elsevier Science Ltd. All rights reserved.

In vivo protein cyclization promoted by a circularly permuted Synechocystis sp PCC6803 DnaB mini-intein

A synthetic Synechocystis sp. PCC6803 DnaB split mini-intein gene was constructed for the in vivo cyclization of recombinant proteins expressed in Escherichia coli. The system was used to cyclize the NH2-terminal domain of E. coli DnaB, the structure of which had been determined previously by NMR spectroscopy. Cyclization was found to proceed efficiently, with little accumulation of precursor, and the product was purified in high yield. The solution structure of cyclic DnaB-N is not significantly different from that of linear DnaB-N and it unfolds reversibly at temperatures similar to14 degreesC higher. Improved hydrogen bonding was observed in the first and last helices, and the length of the last helix was increased, while the 9-amino acid linker used to join the NH2 and COOH termini was found to be highly mobile. The measured thermodynamic stabilization of the structure (DeltaDeltaG approximate to 2 kcal/mol) agrees well with the value estimated from the reduced conformational entropy in the unfolded form. Simple polymer theory can be used to predict likely free energy changes resulting from protein cyclization and how the stabilization depends on the size of the protein and the length of the linker used to connect the termini.

A coiled-coil region of an insect immune suppressor protein is involved in binding and uptake by hemocytes

Polydnaviruses are associated with certain parasitoid wasps and are introduced into the body cavity of the host caterpillar during oviposition. Some of the viral genes are expressed in host tissues and corresponding proteins are secreted into the hemocoel causing suppression of the host immune system. The Cotesia rubecula polydnavirus gene product, CrV1, effectively inactivates hemocytes by mediating cytoskeleton break-down. A precondition for the CrV1 function is the incorporation of the extracellular protein by hemocytes. Here, we show that a coiled-coil domain containing a putative leucine zipper is required for CrV1 function, since removal of this domain abolishes binding and uptake of the CrV1 protein by hemocytes. (C) 2002 Elsevier Science Ltd. All rights reserved.

p53 Suppresses c-Myb-induced trans-Activation and Transformation by Recruiting the Corepressor mSin3A

p53 is known to repress transcription of a number of genes, but the mechanism of p53 recruitment to these target genes is unknown. The c-myb proto-oncogene product (c-Myb) positively regulates proliferation of immature hematopoietic cells, whereas p53 blocks cell cycle progression. Here, we demonstrate that p53 inhibits c-Myb-induced transcription and transformation by directly binding to c-Myb. The ability of c-Myb to maintain the undifferentiated state of M1 cells was also suppressed by p53. p53 did not affect the ability of c-Myb to bind to DNA but formed a ternary complex with the corepressor mSin3A and c-Myb. Thus, p53 antagonizes c-Myb by recruiting mSin3A to down-regulate specific Myb target genes.

H-1 NMR spectroscopic studies of the stability of an oxazolidine condensation product

Condensation of (-)-norephedrine with excess formaldehyde under mild conditions leads to formation of the 2:1 condensation product N,N'-methylenebis(4-methyl-5-phenyl)oxazolidine compared with the reaction with 1 mol of formaldehyde, which leads to 4-methyl-5-phenyloxazolidine. H-1 and C-13 NMR spectroscopy was used to monitor the stability of this compound and its decomposition products. The 2:1 condensation product is found to be stable in CDC1(3) but breaks down rapidly in CD3OD to yield a 50:50 mixture of 4-methyl-5-phenyloxazolidine and 3-hydroxymethyl-4-methyl-5-phenyloxazolidine. Upon addition of D2O to this equimolar mixture, the latter compound decomposes to norephedrine and formaldehyde, whereas the former compound is stable. (C) 1997 by John Wiley & Sons, Ltd.