Modeling the precision and robustness of Hunchback border during Drosophila embryonic development

During anterior-posterior axis specification in the Drosophila embryo, the Hunchback (Hb) protein forms a sharp boundary at the mid-point of the embryo with great positional precision. While Bicoid (Bcd) is a known upstream regulator for hb expression, there is evidence to suggest that Hb effectively filters out “noisy” data received from varied Bcd gradients. We use mathematical models to explore simple regulatory networks which filter out such noise to produce a precise Hb boundary. We find that in addition to Bcd and Hb, at least one freely evolving protein is necessary. An automated search yields a number of examples of three-protein networks exhibiting the desired precision. In all such networks, Hb diffuses much slower than the third protein. In addition, the action of Hb on the third protein is the opposite of the action of the third protein on hb (i.e. if Hb activates the third protein, then the third protein inhibits hb expression, and vice versa). Most of the discovered systems satisfy the known biological properties, that Bcd activates hb, and that Hb activates its own expression. We find that all network topologies satisfying these constraints arise among the networks exhibiting the desired precision. Investigating the dynamics of these networks, we find that under a general class of non-uniform initial conditions, Bcd can be eliminated from the system and the spatiotemporal evolution of these two proteins alone is sufficient to recapture the dynamics. We hypothesize that Bcd is needed only to spatially disturb the gradient of the third protein, and then becomes unnecessary in the further evolution of the Hb border. This provides a possible explanation as to why the Hb dynamics are robust under perturbations of the Bcd gradient. Under this hypothesis, other proteins would be able to assume the role of Bcd in our simulations (possibly in the case of evolutionary divergences or a redundancy in the process), with the only constraint that they act to positively regulate hb.     

Collective decision-making and behavioral polymorphism in group living organisms

Collective foraging in group living animal populations displaying behavioral polymorphism is considered. Using mathematical modeling it is shown that symmetric, spatially homogeneous (food sources are used equally) and asymmetric, spatially inhomogeneous (only one food source is used) regimes can coexist, as a result of differential amplification of choice depending on behavioral type. The model accounts for recent experimental results on social caterpillars not only confirming this coexistence, but also showing the relationship between the two types of regime and the ratio of active to inactive individuals.     

苯那普利对糖尿病大鼠肾小球MMP-9影响的实验研究

目的探讨苯那普利对糖尿病大鼠肾小球基质金属蛋白酶-9（MMP-9）表达的影响。方法将雄性SD大鼠30只,随机抽出10只动物作为正常对照组（C组）,其余20只采用腹腔注射STZ制备糖尿病大鼠模型。将这20只造模成功的糖尿病大鼠模型随机分为2组：糖尿病组（DM组）和药物处理组（DB组）。其中药物处理组即DB组大鼠每天应用苯那普利进行治疗性灌胃,C组和DM组用等量自来水灌胃。分别于治疗4、8周后,测定皿糖和内生肌苷清除率,用代谢笼收集24h尿量。利用免疫组织化学染色及图像分析方法,对MMP-9在肾组织不同时间点的含量变化进行定性分析。结果苯那普利治疗组血糖水平和内生肌苷清除率明显降低,与糖尿病模型组相比有显著性差异。免疫组织化学染色显示：MMP-9在正常肾小球有阳性表达;DM组大鼠病程第8周时,MMP-9在肾小球表达明显减弱;DB组大鼠病程第8周时,经苯那普利治疗后,与同期糖尿病组比较表达明显上调。图像分析方法定量分析免疫组织化学染色阳性面积百分比,苯那普利治疗后第4周和第8周糖尿病肾病的大鼠肾小球MMP-9的阳性面积百分比分别为25．55％±3．44％和20．10％±2．11％,与糖尿病组对比有明显的提高（P〈0．01）。结论MMP-9在糖尿病肾小球中的表达随病程进展逐渐减弱,苯那普利可能部分通过上调MMP-9在肾小球中的表达起到保护肾脏的作用。     

计算机辅助设计可溶性BLyS受体, eBCMA-Fc 融合蛋白及其在大肠杆菌中的表达

B 细胞成熟抗原 (BCMA)是 B淋巴细胞刺激因子(BLyS)的受体之一．它的胞外区与人IgG1 Fc的融合蛋白eBCMA-Fc,又称为诱饵受体,具有拮抗BLyS的活性．为了设计新的拮抗肽,基于BCMA和Fc的晶体结构,通过计算机图形学技术、分子模拟方法,建立了eBCMA-Fc融合蛋白的三维理论结构．利用均方根位移（root mean square distance, RMSD）对eBCMA-Fc融合蛋白与单体eBCMA、Fc构象差异进行分析．融合蛋白eBCMA-Fc中的eBCMA段与单体eBCMA的主链碳原子间RMSD值为0.036 nm,Fc段与单体Fc的主链碳原子间RMSD值为0.064 nm．结果表明,对比单体,融合蛋白eBCMA-Fc并未因eBCMA与Fc直接连接而发生构象的变化．分子对接方法显示,融合蛋白eBCMA-Fc中的BCMA与BLyS作用,而Fc扮演着稳定BCMA构象的支架作用．为进一步验证上述理论分析,构建eBCMA-Fc融合基因,并将载有eBCMA-Fc融合基因的原核表达质粒转化BL21 (DE3)菌、在细菌中表达．目的蛋白经蛋白A亲和柱纯化大约为36 kD,与理论预测值34 kD相近．免疫印迹表明抗人IgG抗体能够识别eBCMA-Fc融合蛋白．ELISA证实,eBCMA-Fc融合蛋白能够结合BLyS．随着eBCMA-Fc融合蛋白增加,结合BLyS的融合蛋白也相应增加．而对照人IgG,即使在高浓度条件下,也不结合BLyS．此外,eBCMA-Fc 融合蛋白能够抑制BLyS对B细胞肿瘤Daudi细胞的作用．这些研究为下一步设计和筛选BLyS拮抗肽提供了实验基础．     

Overview of Mathematical Approaches Used to Model Bacterial Chemotaxis II: Bacterial Populations

We review the application of mathematical modeling to understanding the behavior of populations of chemotactic bacteria. The application of continuum mathematical models, in particular generalized Keller–Segel models, is discussed along with attempts to incorporate the microscale (individual) behavior on the macroscale, modeling the interaction between different species of bacteria, the interaction of bacteria with their environment, and methods used to obtain experimentally verified parameter values. We allude briefly to the role of modeling pattern formation in understanding collective behavior within bacterial populations. Various aspects of each model are discussed and areas for possible future research are postulated.     

Myostatin三维结构模建及分子进化分析

Myostatin(MST)为肌肉生长负调节因子,其功能受抑制可导致肌肉量增加.对MST核酸序列进行序列比对,构建进化树;采用同源模建方法首次模建MST成熟肽生物活性二聚体的四级结构,并预测MST与其受体ActRIIB的相互作用模式.进化树将肌肉生长抑制素基因(MSTN)分成4个亚家族:哺乳动物MSTN,鸟类MSTN以及鱼类MSTN 1和2.MST受纯化选择作用,在不同物种的直系同源基因具有较高的刚源性,其中哺乳动物、鸟类MST C端活性肽氨基酸序列高度保守.表明哺乳动物、鸟类MST的结构、功能类似,且信号传导路径可能一致;而鱼类MST的调控机制可能存在较大差异.MST结构及其表面静电势和疏水氨基酸分布表明静电力和疏水相互作用在MST与其受体结合过程中可能起到十分重要的作用.     

中国萤火虫云南窗萤荧光素酶cDNA的克隆表达和序列分析(英文）

从一种来自中国日行性萤火虫（云南窗萤）发光器官mRNA中克隆、测序并表达了有功能的荧光素酶。云南窗萤荧光素酶的cDNA序列有1 647个碱基,编码548个氨基酸残基。从推测得到的氨基酸序列的比对分析得出：云南窗萤的荧光素酶与来自Lampyris noctiluca, L. turkestanicus和Nyctophila cf. caucasica三种萤火虫的荧光素酶有97.8%的序列一致性。从推测得出的氨基酸序列进行系统发育分析,其结果表明：云南窗萤和Lampyris+Nyctophila聚在一起, 与同属的发光强夜行性的萤火虫不形成的单系。云南窗萤荧光素酶在大肠杆菌中表达的条带大约70 kDa,并且在有荧光素存在时发出黄绿色荧光。对荧光素酶的结构模拟和分析表明,云南窗萤荧光素酶基因的氨基端和羧基端结构域之间的裂沟处存在这5个多肽环,这正是从其他荧光素酶推测得到的催化荧光反应时的底物结合位点。云南窗萤和窗萤属的其他3种萤火虫的荧光素酶相比,有13个不同氨基酸位点,位于模拟分子结构的表面。对于这些多肽环、不同氨基酸残基和晶体结构的进一步研究有利于解释日行和夜行性萤火虫荧光素酶的差异。     

The ecological effect of phenotypic plasticity — Analyzing complex interaction networks (COIN) with agent-based models

Analyzing complex dynamics of ecological systems is complicated by two important facts: First, phenotypic plasticity allows individual organisms to adapt their reaction norms in terms of morphology, anatomy, physiology and behavior to changing local environmental conditions and trophic relationships. Secondly, individual reactions and ecological dynamics are often determined by indirect interactions through reaction chains and networks involving feedback processes.

We present an agent-based modeling framework which allows to represent and analyze ecological systems that include phenotypic changes in individual performances and indirect interactions within heterogeneous and temporal changing environments. We denote this structure of interacting components as COmplex Interaction Network (COIN).

Three examples illustrate the potential of the system to analyze complex ecological processes that incorporate changing phenotypes on the individual level:

• A model on fish population dynamics of roach (Rutilus rutilus) leads to a differentiation in fish length resulting in a conspicuous distribution that influences reproduction capability and thus indirectly the fitness.

• Modeling the reproduction phase of the passerine bird Erithacus rubecula (European Robin) illustrates variation in the behavior of higher organisms in dependence of environmental factors. Changes in reproduction success and in the proportion of different activities are the results.

• The morphological reaction of plants to changes in fundamental environmental parameters is illustrated by the black alder (Alnus glutinosa) model. Specification of physiological processes and the interaction structure on the level of modules allow to represent the reaction to changes in irradiance and temperature accurately.

Applying the COIN-approach, individual plasticity emerges as a structural and functional implication in a self-organized manner. The examples illustrate the potential to integrate existing approaches to represent detailed and complex traits for higher order organisms and to combine ecological and evolutionary aspects.

Mathematical Model of Growth and Heterologous Hantavirus Protein Production of the Recombinant Yeast Saccharomyces cerevisiae

A segregated mathematical model was developed for the analysis and interpretation of cultivation data of growth of the recombinant yeast Saccharomyces cerevisiae on multiple substrates (glucose, maltose, pyruvate, ethanol, acetate, and galactose). The model accounts for substrate consumption, plasmid stability, and production level of a model protein, a modified nucleocapsid protein of the Puumala virus. Recombinant nucleocapsid proteins from different Hantaviruses have previously been demonstrated as suitable antigens for diagnostics as well as for sero‐epidemiological studies. The model is based on a system of 10 nonlinear ordinary differential equations and accounts for the influence of various factors, e.g., selective pressure for enhancing plasmid stability by formaldehyde or the toxic effects of the intracellular accumulation of the heterologous protein on cell growth and product yield. The model allows the growth of two populations of cells to be simulated: plasmid‐bearing and plasmid‐free yeast cells, which have lost the plasmid during cultivation. Based on the model, sensitivity studies in respect to parameter changes were performed. These enabled, for example, the evaluation of the impact of an increase in the initial concentration of nutrients and growth factors (e.g., vitamins, microelements, etc.) on the biomass yield and the heterologous protein production level. As expected, the productivity of the heterologous protein in S. cerevisiae is closely correlated with plasmid stability. The 25 free model parameters, including the yield coefficients for different growth stages and dynamic constants, were estimated by nonlinear techniques, and the model was validated against a data set not used for parameter estimation. The simulation results were found to be in good agreement with the experimental data.     

The effects of the false vocal fold gaps on intralaryngeal pressure distributions and their effects on phonation

Human phonation does not merely depend on the vibration of the vocal folds. Research by clinical and computer simulations has demonstrated that the false vocal fold (FVF) is an important laryngeal con-striction that plays a vital role during human voice production. This study explored the effects of the FVF gaps using both the three-dimensional Plexiglas model and the numerical computation methods. Twelve FVF gaps (ranging from 0.02 to 2.06 cm) were used in this study at three glottal angles (uniform and convergent/divergent 40°), two minimal glottal diameters (Dg) (0.04 cm and 0.06 cm) separately, and the constant subglottal pressure (8 cm H2O). The results suggested that (1) the intralaryngeal pressure was the lowest and the flow was the highest (least flow resistance) when the FVF gap was 1.5-2 times greater than Dg; (2) the divergent glottal angle gave lower pressure and greater flow than the conver-gent and uniform glottal angle as there were no FVF conditions; (3) the presence of the FVF decreased the effects of the glottal angle to a certain extent; and more importantly, (4) the presence of the FVF also moved the separation points downstream, straightened the glottal jet for a longer distance, decreased the overall laryngeal resistance, and reduced the energy dissipation, suggesting the significance of FVF in efficient voice production. These results may be incorporated in the phonatory models (physical or computational) for better understanding of vocal mechanics. The results might also be helpful in exploring the surgical and rehabilitative intervention of related voice problems.