定量分析调控人类免疫缺陷病毒潜伏与激活的动力学机制

目的 治疗艾滋病最大的障碍在于无法根除人类免疫缺陷病毒（HIV）潜伏于人体细胞所形成的病毒存储库。构建描述病毒存储库建立分子机制的动力学模型需考虑生物体内的噪声环境和多重影响因素,本文通过一种全新的动力学结构分解方法将随机微分方程的确定性部分与随机性噪声分开,从而在仅需分析常微分方程不动点的情况下即可判断不同药物靶点的作用效果。方法 使用连续的随机微分方程构建了HIV转录过程的动力学模型,简化了描述系统所需方程的维度,增大了模型的可探索空间,在此基础上,通过计算得到的势能函数和概率分布函数直观表示病毒潜伏与激活的不同表达状态以及它们之间的关系。结果 定量分析了不同动力学参数对系统稳态和势函数的影响程度,分别得到了系统处于双稳态和单稳态时的参数范围,并将不同因素对动力系统分岔的影响程度与生物学实验结果对比,验证了本工作的理论基础。结论 本文突破了以往离散、随机的方法,可以通过常微分方程定量分析HIV转录调控的动力学机制,有利于推广到处理高维情况,进一步研究艾滋病在生物体内的发生发展,从而指导设计实验寻找临床上的治疗方案。     

Evaluation of Fourier-transform infrared spectroscopy for data capture in predictive microbiology

Chemical changes in the medium, induced by the fermentative species Lactobacillus plantarum and Lactobacillus brevis and by the enzymatic action of a proteolytic, spoilage species, Yarrowia lipolytica, were analysed using Fourier-transform i.r. spectroscopy (FTIR). Changes in the absorbance data over time could be modelled using one of the more current predictive, mathematical models of microbial growth, such as the Gompertz equation. Moreover, a linear correlation between FTIR data (expressed as absorbance of some selected peaks) and viability data (expressed as log₁₀ c.f.u./g or ml) was observed during the fermentation process, both for L. plantarum and L. brevis.     

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.     

生物种群动态微分方程模型参数估计方法

本文以数值分析和最优化技术的有机结合为基础,提出了一种新的对动态微分方程模型直接进行数据拟合和参数估计方法,并以Logistic微分方程、生物种间竞争关系微分方程以及一种复合形态的Logistic微分方程为例进行了数据拟合试验．结果表明,该方法对各种动态微分方程模型均能进行最优拟合分析并求解其参数．同时发现,以前有的作者〔1,2,3,4,5〕提出的方法所得到的参数估计值存在系统误差且误差较大．     

Prediction of titration properties of structures of a protein derived from molecular dynamics trajectories.

This paper explores the dependence of the molecular dynamics (MD) trajectory of a protein molecule on the titration state assigned to the molecule. Four 100-ps MD trajectories of bovine pancreatic trypsin inhibitor (BPTI) were generated, starting from two different structures, each of which was held in two different charge states. The two starting structures were the X-ray crystal structure and one of the solution structures determined by NMR, and the charge states differed only in the ionization state of N terminus. Although it is evident that the MD simulations were too short to sample fully the equilibrium distribution of structures in each case, standard Poisson-Boltzmann titration state analysis of the resulting configurations shows general agreement between the overall titration behavior of the protein and the charge state assumed during MD simulation: at pH 7, the total net charge of the protein resulting from the titration analysis is consistently lower for the protein with the N terminus assumed to be neutral than for the protein with the N terminus assumed to be charged. For most of the ionizable residues, the differences in the calculated pKaS among the four trajectories are statistically negligible and remain in good agreement with the data obtained by crystal structure titration and by experiment. The exceptions include the N terminus, which responds directly to the change of its imposed charge; the C terminus, which in the NMR structure interacts strongly with the former; and a few other residues (Arg 1, Glu 7, Tyr 35, and Arg 42) whose pKaS reflect the initial structure and the limited trajectory lengths. This study illustrates the importance of the careful assignment of protonation states at the start of MD simulations and points to the need for simulation methods that allow for the variation of the protonation state in the calculation of equilibrium properties.     

Cooperativity and effects of ionic strength on the binding of sodium n-dodecyl sulphate to lysozyme

The binding of sodium n-dodecyl sulphate to lysozyme has been measured by equilibrium dialysis at 25°C and pH 3.2 over a range of ionic strengts from 0.0119 to 0.2119. Binding isotherms in the region corresponding to ionic binding between the surfactant anions and cationic amino acid residues on the protein have been interpreted in terms of the Hill equation and exhibit positive cooperativity with Hill coefficients in the region of 7–11. The Gibbs energies of binding have been calculated from the Hill binding constants and from the Wyman binding potentials. The stability of the surfactant-protein complexes is discussed in relation to the stability of surfactant micelles. Ionic binding of the surfactant is weakened and hydrophobic binding strengthened by increasing ionic strength.     

Holling's “hungry mantid” model for the invertebrate functional response considered as a Markov process. III. Stable satiation distribution

In this paper, we study an analytical model describing predatory behaviour. It is assumed that the parameter describing the predator's behaviour is its satiation. Using semigroup methods and compactness arguments we prove that a stable satiation distribution is reached if t. Furthermore, using a Trotter-Kato theorem we justify the transition to the much simpler problem that is obtained if the prey biomass tends to zero.     

Glucose-6-phosphate dehydrogenase porbandar: A new slow variant with slightly reduced activity in a South African family of Indian descent

A new variant of the red cell enzyme glucose-6-phosphate dehydrogenase has been detected in a South African male of Indian descent and in several of his relatives. The enzyme variant is characterized by slow electrophoretic mobility, low Michaelis constants for the substrates glucose-6-phosphate and NADP, and increased utilization of the substrate analogues 2-deoxyglucose-6-phosphate and deamino-NADP relative to the normal (B⁺) enzyme. There is no evidence that the enzyme variant, for which the name G6PD Porbandar is suggested, is associated with any hematological abnormality.The Atomic Energy Board and the South African Medical Research Council provided support for part of this work.     

Role of unstirred layer in intestinal absorption of phenylalanine in vivo

The appearance rate of l- and d-phenylalanine in the venous blood of rat jejunal loops in vivo is increased up to 60% if the intraluminal solution is mixed more efficiently by the simultaneous perfusion of air. The effect decreases as the luminal concentration is increased to 100 mmol/1. Thus, the apparent Michaelis constants are by 50% lower in the case of the reduced unstirred layer (26 to 17 for l- and 9 to 6 mmol/1 for d-phenylalanine).The enhancement of the absorption and the reduction of the Michaelis constants can be attributed to the reduction of the effective unstirred layer thickness by about 400–500 μm.