Protein fold recognition score functions: Unusual construction strategies

We describe two ways of optimizing score functions for protein sequence to structure threading. The first method adjusts parameters to improve sequence to structure alignment. The second adjusts parameters so as to improve a score function's ability to rank alignments calculated in the first score function. Unlike those functions known as knowledge-based force fields, the resulting parameter sets do not rely on Boltzmann statistics, have no claim to representing free energies and are purely constructions for recognizing protein folds. The methods give a small improvement, but suggest that functions can be profitably optimized for very specific aspects of protein fold recognition, Proteins 1999;36:454-461. (C) 1999 Wiley-Liss, Inc.

Structure optimization combining soft-core interaction functions, the diffusion equation method, and molecular dynamics

Smoothing the potential energy surface for structure optimization is a general and commonly applied strategy. We propose a combination of soft-core potential energy functions and a variation of the diffusion equation method to smooth potential energy surfaces, which is applicable to complex systems such as protein structures; The performance of the method was demonstrated by comparison with simulated annealing using the refinement of the undecapeptide Cyclosporin A as a test case. Simulations were repeated many times using different initial conditions and structures since the methods are heuristic and results are only meaningful in a statistical sense.

Protein fold recognition without Boltzmann statistics or explicit physical basis

We present a fast method for finding optimal parameters for a low-resolution (threading) force field intended to distinguish correct from incorrect folds for a given protein sequence. In contrast to other methods, the parameterization uses information from >10(7) misfolded structures as well as a set of native sequence-structure pairs. In addition to testing the resulting force field's performance on the protein sequence threading problem, results are shown that characterize the number of parameters necessary for effective structure recognition.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Load-unload response ratio and accelerating moment/energy release critical region scaling and earthquake prediction

The main idea of the Load-Unload Response Ratio (LURR) is that when a system is stable, its response to loading corresponds to its response to unloading, whereas when the system is approaching an unstable state, the response to loading and unloading becomes quite different. High LURR values and observations of Accelerating Moment/Energy Release (AMR/AER) prior to large earthquakes have led different research groups to suggest intermediate-term earthquake prediction is possible and imply that the LURR and AMR/AER observations may have a similar physical origin. To study this possibility, we conducted a retrospective examination of several Australian and Chinese earthquakes with magnitudes ranging from 5.0 to 7.9, including Australia's deadly Newcastle earthquake and the devastating Tangshan earthquake. Both LURR values and best-fit power-law time-to-failure functions were computed using data within a range of distances from the epicenter. Like the best-fit power-law fits in AMR/AER, the LURR value was optimal using data within a certain epicentral distance implying a critical region for LURR. Furthermore, LURR critical region size scales with mainshock magnitude and is similar to the AMR/AER critical region size. These results suggest a common physical origin for both the AMR/AER and LURR observations. Further research may provide clues that yield an understanding of this mechanism and help lead to a solid foundation for intermediate-term earthquake prediction.

Theory of strongly correlated electron systems. I. Intersite coulomb interaction and the approximation of renormalized fermions in total energy calculations

The diagrammatic strong-coupling perturbation theory (SCPT) for correlated electron systems is developed for intersite Coulomb interaction and for a nonorthogonal basis set. The construction is based on iterations of exact closed equations for many - electron Green functions (GFs) for Hubbard operators in terms of functional derivatives with respect to external sources. The graphs, which do not contain the contributions from the fluctuations of the local population numbers of the ion states, play a special role: a one-to-one correspondence is found between the subset of such graphs for the many - electron GFs and the complete set of Feynman graphs of weak-coupling perturbation theory (WCPT) for single-electron GFs. This fact is used for formulation of the approximation of renormalized Fermions (ARF) in which the many-electron quasi-particles behave analogously to normal Fermions. Then, by analyzing: (a) Sham's equation, which connects the self-energy and the exchange- correlation potential in density functional theory (DFT); and (b) the Galitskii and Migdal expressions for the total energy, written within WCPT and within ARF SCPT, a way we suggest a method to improve the description of the systems with correlated electrons within the local density approximation (LDA) to DFT. The formulation, in terms of renormalized Fermions LIDA (RF LDA), is obtained by introducing the spectral weights of the many electron GFs into the definitions of the charge density, the overlap matrices, effective mixing and hopping matrix elements, into existing electronic structure codes, whereas the weights themselves have to be found from an additional set of equations. Compared with LDA+U and self-interaction correction (SIC) methods, RF LDA has the advantage of taking into account the transfer of spectral weights, and, when formulated in terms of GFs, also allows for consideration of excitations and nonzero temperature. Going beyond the ARF SCPT, as well as RF LIDA, and taking into account the fluctuations of ion population numbers would require writing completely new codes for ab initio calculations. The application of RF LDA for ab initio band structure calculations for rare earth metals is presented in part 11 of this study (this issue). (c) 2005 Wiley Periodicals, Inc.

C-terminal charged cluster of MscL, RKKEE, functions as a pH sensor

The RKKEE cluster of charged residues located within the cytoplasmic helix of the bacterial mechanosensitive channel, MscL, is essential for the channel function. The structure of MscL determined by x-ray crystallography and electron paramagnetic resonance spectroscopy has revealed discrepancies toward the C-terminus suggesting that the structure of the C-terminal helical bundle differs depending on the pH of the cytoplasm. In this study we examined the effect of pH as well as charge reversal and residue substitution within the RKKEE cluster on the mechanosensitivity of Escherichia coli MscL reconstituted into liposomes using the patch-clamp technique. Protonation of either positively or negatively charged residues within the cluster, achieved by changing the experimental pH or residue substitution within the RKKEE cluster, significantly increased the free energy of activation for the MscL channel due to an increase in activation pressure. Our data suggest that the orientation of the C-terminal helices relative to the aqueous medium is pH dependent, indicating that the RKKEE cluster functions as a proton sensor by adjusting the channel sensitivity to membrane tension in a pH-dependent fashion. A possible implication of our results for the physiology of bacterial cells is briefly discussed.

Titanium phosphate for Fuel Cell Proton Conduction Membranes

Environmental issues due to increases in emissions of air pollutants and greenhouse gases are driving the development of clean energy delivery technologies such as fuel cells. Low temperature Proton Exchange Membrane Fuel Cells (PEMFC) use hydrogen as a fuel and their only emission is water. While significant advances have been made in recent years, a major limitation of the current technology is the cost and materials limitations of the proton conduction membrane. The proton exchange membrane performs three critical functions in the PEMFC membrane electrode assembly (MEA): (i) conduction of protons with minimal resistance from the anode (where they are generated from hydrogen) to the cathode (where they combine with oxygen and electrons, from the external circuit or load), (ii) providing electrical insulation between the anode and cathode to prevent shorting, and (iii) providing a gas impermeable barrier to prevent mixing of the fuel (hydrogen) and oxidant. The PFSA (perfluorosulphonic acid) family of membranes is currently the best developed proton conduction membrane commercially available, but these materials are limited to operation below 100oC (typically 80oC, or lower) due to the thermochemical limitations of this polymer. For both mobile and stationary applications, fuel cell companies require more durable, cost effective membrane technologies capable of delivering enhanced performance at higher temperatures (typically 120oC, or higher. This is driving research into a wide range of novel organic and inorganic materials with the potential to be good proton conductors and form coherent membranes. There are several research efforts recently reported in the literature employing inorganic nanomaterials. These include functionalised silica phosphates [1,2], fullerene [3] titania phosphates [4], zirconium pyrophosphate [5]. This work addresses the functionalisation of titania particles with phosphoric acid. Proton conductivity measurements are given together with structural properties.

Modal functions of future tenses in French

The paper disputes two influential claims in the Romance Linguistics literature. The first is that the synthetic future tenses in spoken Western Romance are now rivalled, if not supplanted, as temporal functors by the more recently developed GO futures. The second is that these synthetic futures now have modal rather than temporal meanings in spoken Romance. These claims are seen as reflecting a universal cycle of diachronic change, in which verb forms originally expressing modal (or aspectual) values take on future temporal reference, becoming tenses. The new modal meanings supplant the temporal, which are then taken up by new forms. Challenges to this theory for French are raised on the basis of empirical evidence of two sorts. Positively, future tenses in spoken Romance continue to be used with temporal meaning. Negatively, evidence of modal meaning for these forms is lacking. The evidence comes froma corpora of spoken French, native speaker judgements and verb data from a daily broadsheet. Cumulatively, it points to the reverse of the claims noted above: the synthetic future in spoken French has temporal but little modal meaning.

Correcting physical activity energy expenditure for body size in young children

Being able to compare the energy cost of physical activity across and between populations is important. However, energy expenditure is related to body size, so it is necessary to appropriately adjust for differences in body size when comparisons are made. This study examined the relationship between the daily energy cost of activity and body weight in 47 children aged 6-10 years. Log-log regression showed weight(1.0) to be an inappropriate adjustment for activity energy expenditure in children, with a more valid adjustment being weight(0.3). Clearly, both weight dependent and non-weight dependent activities are part of everyday living in children. This balance influences how energy expenditure is correctly adjusted for body size. Investigators interpreting data of energy expenditure in children from children of different body sizes need to take this into consideration.

A first-principles density-functional calculation of the electronic and vibrational structure of the key melanin monomers

We report first-principles density-functional calculations for hydroquinone (HQ), indolequinone (IQ), and semiquinone (SQ). These molecules are believed to be the basic building blocks of the eumelanins, a class of biomacromolecules with important biological functions (including photoprotection) and with the potential for certain bioengineering applications. We have used the difference of self-consistent fields method to study the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital, HL. We show that HL is similar in IQ and SQ, but approximately twice as large in HQ. This may have important implications for our understanding of the observed broadband optical absorption of the eumelanins. The possibility of using this difference in HL to molecularly engineer the electronic properties of eumelanins is discussed. We calculate the infrared and Raman spectra of the three redox forms from first principles. Each of the molecules have significantly different infrared and Raman signatures, and so these spectra could be used in situ to nondestructively identify the monomeric content of macromolecules. It is hoped that this may be a helpful analytical tool in determining the structure of eumelanin macromolecules and hence in helping to determine the structure-property-function relationships that control the behavior of the eumelanins.

Hydraulic Design of Stepped Spillways and Downstream Energy Dissipators for Embankment Dams

In recent years, the design flows of many dams were re-evaluated, often resulting in discharges larger than the original design. In many cases, the occurrence of the revised flows could result in dam overtopping because of insufficient storage and spillway capacity. An experimental study was conducted herein to gain a better understanding of the flow properties in stepped chutes with slopes typical of embankment dams. The work was based upon a Froude similitude in large-size experimental facilities. A total of 10 configurations were tested including smooth steps, steps equipped with devices to enhance energy dissipation and rough steps. The present results yield a new design procedure. The design method includes some key issues not foreseen in prior studies : e.g., gradually varied flow, type of flow regime, flow resistance. It is believed that the outcomes are valid for a wide range of chute geometry and flow conditions typical of embankment chutes.

Least energy solutions of a critical Neumann problem with a weight

We investigate the effect of the coefficient of the critical nonlinearity for the Neumann problem on the existence of least energy solutions. As a by-product we establish a Sobolev inequality with interior norm.

Algebraic Bethe ansatz method for the exact calculation of energy spectra and form factors: applications to models of Bose-Einstein condensates and metallic nanograins

In this review we demonstrate how the algebraic Bethe ansatz is used for the calculation of the-energy spectra and form factors (operator matrix elements in the basis of Hamiltonian eigenstates) in exactly solvable quantum systems. As examples we apply the theory to several models of current interest in the study of Bose-Einstein condensates, which have been successfully created using ultracold dilute atomic gases. The first model we introduce describes Josephson tunnelling between two coupled Bose-Einstein condensates. It can be used not only for the study of tunnelling between condensates of atomic gases, but for solid state Josephson junctions and coupled Cooper pair boxes. The theory is also applicable to models of atomic-molecular Bose-Einstein condensates, with two examples given and analysed. Additionally, these same two models are relevant to studies in quantum optics; Finally, we discuss the model of Bardeen, Cooper and Schrieffer in this framework, which is appropriate for systems of ultracold fermionic atomic gases, as well as being applicable for the description of superconducting correlations in metallic grains with nanoscale dimensions.; In applying all the above models to. physical situations, the need for an exact analysis of small-scale systems is established due to large quantum fluctuations which render mean-field approaches inaccurate.