STABILITY IN THE PROPORTION OF HARBOR SEALS HAULED OUT UNDER LOCALLY IDEAL CONDITIONS

We monitored the haul-out behavior of 68 radio-tagged harbor seals ( Phoca vitulina ) during the molt season at two Alaskan haul-out sites (Grand Island, August-September 1994; Nanvak Bay, August-September 2000). For each site, we created a statistical model of the proportion of seals hauled out as a function of date, time of day, tide, and weather covariates. Using these models, we identified the conditions that would result in the greatest proportion of seals hauled out. Although those "ideal conditions" differed between sites, the proportion of seals predicted to be hauled out under those conditions was very similar (81.3% for Grand Island and 85.7% for Nanvak Bay). The similar estimates for both sites suggest that haul-out proportions under locally ideal conditions may be constant between years and geographic regions, at least during the molt season.     

Using Molecular Dynamics to Predict Factors Affecting Binding Strength and Magnetic Relaxivity of MRI Contrast Agents

We demonstrate the use of molecular dynamics and molecular mechanics methods to calculate properties and behavior of metal-chelate complexes that can be used as MRI contrast agents. Static and dynamic properties of several known agents were calculated and compared with experiment. We calculated the static properties such as the q-values (number of inner shell waters) and binding distances of chelate atoms to the metal ion for a set of chelates with known X-ray structure. The dynamic flexibility of the chelate arms was also calculated. These computations were extended to a series of exploratory chelate structures in order to estimate their potential as MRI contrast agents. We have also calculated for the first time the NMR relaxivity of an MRI contrast agent using a long (5 nsec) molecular dynamics simulation. Our predictions are promising enough that the method should prove useful for evaluating novel candidate compounds before they are synthesized. One novel static property, the projected area of chelate atoms onto a virtual surface centered on the metal ion (gnomonic projection), was found to give an effective measure of how well the chelate atoms use the free space around the metal ion.     

Conformational analysis and clustering of short and medium size loops connecting regular secondary structures: a database for modeling and prediction.

Loops are regions of nonrepetitive conformation connecting regular secondary structures. We identified 2,024 loops of one to eight residues in length, with acceptable main-chain bond lengths and peptide bond angles, from a database of 223 protein and protein-domain structures. Each loop is characterized by its sequence, main-chain conformation, and relative disposition of its bounding secondary structures as described by the separation between the tips of their axes and the angle between them. Loops, grouped according to their length and type of their bounding secondary structures, were superposed and clustered into 161 conformational classes, corresponding to 63% of all loops. Of these, 109 (51% of the loops) were populated by at least four nonhomologous loops or four loops sharing a low sequence identity. Another 52 classes, including 12% of the loops, were populated by at least three loops of low sequence similarity from three or fewer nonhomologous groups. Loop class suprafamilies resulting from variations in the termini of secondary structures are discussed in this article. Most previously described loop conformations were found among the classes. New classes included a 2:4 type IV hairpin, a helix-capping loop, and a loop that mediates dinucleotide-binding. The relative disposition of bounding secondary structures varies among loop classes, with some classes such as beta-hairpins being very restrictive. For each class, sequence preferences as key residues were identified; those most frequently at these conserved positions than in proteins were Gly, Asp, Pro, Phe, and Cys. Most of these residues are involved in stabilizing loop conformation, often through a positive phi conformation or secondary structure capping. Identification of helix-capping residues and beta-breakers among the highly conserved positions supported our decision to group loops according to their bounding secondary structures. Several of the identified loop classes were associated with specific functions, and all of the member loops had the same function; key residues were conserved for this purpose, as is the case for the parvalbumin-like calcium-binding loops. A significant number, but not all, of the member loops of other loop classes had the same function, as is the case for the helix-turn-helix DNA-binding loops. This article provides a systematic and coherent conformational classification of loops, covering a broad range of lengths and all four combinations of bounding secondary structure types, and supplies a useful basis for modelling of loop conformations where the bounding secondary structures are known or reliably predicted.     

Ligand recognition in μ opioid receptor: experimentally based modeling of μ opioid receptor binding sites and their testing by ligand docking

For three-dimensional understanding of the mechanisms that control potency and selectivity of the ligand binding at the atomic level, we have analysed opioid receptor-ligand interaction based on the receptor's 3D model. As a first step, we have constructed molecular models for the multiple opioid receptor subtypes using bacteriorhodopsin as a template. The S-activated dihydromorphine derivatives should serve as powerful tools in mapping the three-dimensional structure of the μ opioid receptor, including the nature of the agonist-mediated conformational change that permits G protein-coupling to ‘second messenger’ effector molecules, and in identifying specific ligand-binding contacts with the μ opioid receptor. The analyses of the interactions of some opioid ligands with the predicted ligand binding sites are consistent with the results of the affinity labeling experiments.     

Coding of odor intensity in a steady-state deterministic model of an olfactory receptor neuron

The coding of odor intensity by an olfactory receptor neuron model was studied under steady-state stimulation. Our model neuron is an elongated cylinder consisting of the following three components: a sensory dendritic region bearing odorant receptors, a passive region consisting of proximal dendrite and cell body, and an axon. First, analytical solutions are given for the three main physiological responses: (1) odorant-dependent conductance change at the sensory dendrite based on the Michaelis-Menten model, (2) generation and spreading of the receptor potential based on a new solution of the cable equation, and (3) firing frequency based on a Lapicque model. Second, the magnitudes of these responses are analyzed as a function of odorant concentration. Their dependence on chemical, electrical, and geometrical parameters is examined. The only evident gain in magnitude results from the activation-to-conductance conversion. An optimal encoder neuron is presented that suggests that increasing the length of the sensory dendrite beyond about 0.3 space constant does not increase the magnitude of the receptor potential. Third, the sensivities of the responses are examined as functions of (1) the concentration at half-maximum response, (2) the lower and upper concentrations actually discriminated, and (3) the width of the dynamic range. The overall gain in sensitivity results entirely from the conductance-to-voltage conversion. The maximum conductance at the sensory dendrite appears to be the main tuning constant of the neuron because it determines the shift toward low concentrations and the increase in dynamic range. The dynamic range of the model cannot exceed 5.7 log units, for a sensitivity increase at low odor concentration is compensated by a sensitivity decrease at high odor concentration.     

基于结构比较的蛋白质同源模建系统及其评估Ⅱ侧链的安装

至今,有关蛋白质侧链的同源模建,除了在本体模板上安装侧链和少数限制条件下在同源模板上安装侧链的报道外,系统的研究和实施似乎还未见报道。本软件系统ＰＭＯＤＥＬＩＮＧ采用同源移植和“死端排除“相结合的侧链安装策略,对与模板蛋白相应践基具有相似大小和形状的目标残基采用直接移植的方法。其余铡链则用广义“死端排除定则”安装。经众多蛋白的测试,达到了较好的模建品质。     

Thermal isoelectric precipitation of alpha-lactalbumin from a whey protein concentrate: Influence of protein-calcium complexation

The selective precipitation of alpha-lactalbumin (alpha-LA) at a pH around its isoelectric point (4.2) under heat treatment is the basis for a fractionation process of whey proteins. As precipitation is a phenomenon dependent on the protein hydrophobicity, and as the release of the tightly bound calcium occurring at pH around 4 modifies the alpha-LA hydrophobicity, the specific role of calcium on isoelectric precipitation is investigated. A study of the extent of alpha-LA precipitation in a whey protein concentrate under various operating conditions of pH, temperature, protein concentration, and calcium content is presented. We propose a mechanism for this phenomenon as a combination of a complexation equilibrium and of an irreversible precipitation, to account for the influence of temperature, alpha-LA concentration total ionic content, and calcium concentration, and also to estimate the complexation equilibrium constant. (c) 1995 John Wiley & Sons, Inc.     

A modeling study of anaerobic biofilm systems: II. Reactor modeling

A rigorous steady-state model of anaerobic biofilm reactors taking into account acid-base and gas-phase equilibria in the reactor in conjunction with detailed chemical equilibria and mass transfer in acetate-utilizing methanogenic biofilms is presented. The performances of ideal completely stirred tank reactors (CSTRs) and plug-flow reactors, as well as reactors with nonideal hydraulic conditions, are simulated. Decreasing the surface loading rate increases the acetate removal efficiency, while decreasing the influent pH and increasing the buffering capacity improves the removal efficiency only if the bulk pH of the reactor shifts toward more optimal values between 6.8 to 7.0. The reactor can have negative or positive removal efficiencies depending on the start-up conditions. The respiration coefficient plays a critical role in determining the minimum influent pH required for reactor recovery after failure. Having multiple CSTRs-in-series generally increases the overall removal efficiency for the influent conditions investigated. Monitoring of the influent feed quality is critical for plug-flow reactors, becasue failure of the initial sections of the reactor may cause a cascading effect that may lead to a rapid reactor failure. (c) 1995 John Wiley & Sons, Inc.     

Changes in proteoglycan synthesis of chondrocytes in articular cartilage are associated with the time-dependent changes in their mechanical environment

Explant loading experiments were conducted to investigate the effect of load duration on proteoglycan synthesis. A compressive load of 0.1 MPa applied for 10 min was found to stimulate proteoglycan synthesis, while the same load applied for 20 h suppressed synthesis. This bimodal response suggests that the cells are responding to different mechanical stimuli as time progresses. A theoretical model has therefore been developed to describe the mechanical environment perceived by cells within soft hydrated tissues (e.g. articular cartilage) while the tissue is being loaded. The cells are modeled, using the biphasic theory, as fluid-solid inclusions embedded in and attached to a biphasic extracellular matrix of distinct material properties. A method of solution is developed which is valid for any axisymmetric loading configuration, provided that the cell radius, a, is small relative to the tissue height, h (i.e. h/a 1). A closed-form analytical solution for this inclusion problem is then presented for the confined compression configuration. Results from this model show that the mechanical environment in and around the cells is time dependent and inhomogeneous, and can be significantly influenced by differences in properties between the cell and the extracellular matrix.     

Simulation and dynamic optimisation of animal cell culture

In animal cell culture, there are some 25 substrates that both have a significant effect on the culture performance and which can be measured with relative ease. A detailed dynamic simulation for such a culture has been produced and an optimisation policy that use this model to identify ideal media conditions has been developed. This paper describes an extension of that work to include the dynamic optimisation of cultures under fed-batch operation. Two different types of feeding policy were considered – in the first, discrete shots of feed were supplied, while in the second, feed was added continuously. Both policies offered significant improvements in the predicted productivity of the culture - up to 30% that of an experimentally optimisedbatch culture.