Relating side-chain mobility in proteins to rotameric transitions: Insights from molecular dynamics simulations and NMR

The dynamic aspect of proteins is fundamental to understanding protein stability and function. One of the goals of NMR studies of side-chain dynamics in proteins is to relate spin relaxation rates to discrete conformational states and the timescales of interconversion between those states. Reported here is a physical analysis of side-chain dynamics that occur on a timescale commensurate with monitoring by ²H spin relaxation within methyl groups. Motivated by observations made from tens-of-nanoseconds long MD simulations on the small protein eglin c in explicit solvent, we propose a simple molecular mechanics-based model for the motions of side-chain methyl groups. By using a Boltzmann distribution within rotamers, and by considering the transitions between different rotamer states, the model semi-quantitatively correlates the population of rotamer states with ‘model-free’ order parameters typically fitted from NMR relaxation experiments. Two easy-to-use, analytical expressions are given for converting S²_axis’ values (order parameter for C–CH₃ bond) into side-chain rotamer populations. These predict that S²_axis’ values below 0.8 result from population of more than one rotameric state. The relations are shown to predict rotameric sampling with reasonable accuracy on the ps–ns timescale for eglin c and are validated for longer timescales on ubiquitin, for which side-chain residual dipolar coupling (RDC) data have been collected.     

Bias-corrected maximum likelihood estimator of the negative binomial dispersion parameter

We derive a first-order bias-corrected maximum likelihood estimator for the negative binomial dispersion parameter. This estimator is compared, in terms of bias and efficiency, with the maximum likelihood estimator investigated by Piegorsch (1990, Biometrics46, 863-867), the moment and the maximum extended quasi-likelihood estimators investigated by Clark and Perry (1989, Biometrics45, 309-316), and a double-extended quasi-likelihood estimator. The bias-corrected maximum likelihood estimator has superior bias and efficiency properties in most instances. For ease of comparison we give results for the two-parameter negative binomial model. However, an example involving negative binomial regression is given.     

Modeling late entry bias in survival analysis

In a failure time analysis, we sometimes observe additional study subjects who enter during the study period. These late entries are treated as left-truncated data in the statistical literature. However, with real data, there is a substantial possibility that the delayed entries may have extremely different hazards compared to the other standard subjects. We focus on a situation in which such entry bias might arise in the analysis of survival data. The purpose of the present article is to develop an appropriate methodology for making inference about data including late entries. We construct a model that includes parameters for the effect of delayed entry bias having no specification for the distribution of entry time. We also discuss likelihood inference based on this model and derive the asymptotic behavior of estimates. A simulation study is conducted for a finite sample size in order to compare the analysis results using our method with those using the standard method, where independence between entry time and failure time is assumed. We apply this method to mortality analysis among atomic bomb survivors defined in a geographical study region.     

Transient response of aerobic and anoxic activated sludge activities to sudden substrate concentration changes

The state-of-the-art understanding of activated sludge processes as summarized in activated sludge models (ASMs) predicts an instantaneous increase in the biomass activity (which is measured, e.g., by the corresponding respiration rate OUR, NUR, etc.) under sudden substrate concentration changes. Experimental data (e.g., short-term batch respiration experiments under aerobic or anoxic conditions) collected for the calibration of the dynamic models (ASMs) often exhibit a transient phenomenon while attaining maximum activity, which cannot be explained by the current understanding of the activated sludge process. That transient phenomenon exhibits itself immediately upon addition of a substrate source to an endogenously respiring activated sludge sample and it usually takes a few minutes until the activated sludge reaches its maximum possible rate under given environmental conditions. This discrepancy between the state-of-the-art model and the experimental data is addressed in detail in this investigation. It is shown that the discrepancy is not caused by an error in the experimental set-up/data but it is rather due to model inadequacy. Among the hypotheses proposed, it appears that this transient response of the activated sludge most likely results from the sequence of intracellular reactions involved in substrate degradation by the activated sludge. Results from studies performed elsewhere with pure cultures (S. cerevisae and E. coli) support the hypothesis. The transient phenomenon can be described by a dynamic metabolic network model or by a simple first-order model, as adopted in this study. The transient phenomenon occurring in short-term batch respiration experiments is shown to interfere severely with parameter estimation if not modeled properly (2.8%, 11.5%, and 16.8% relative errors [average of three experiments] on Y(H), micro(maxH), and K(S), respectively). Proper modeling of this transient phenomenon whose time constant is on the order of minutes (1 to 3 min) is expected to contribute fundamentally to a better understanding and modeling of Orbal, carousel, and SBR-type treatment plants with fast-alternating process conditions, although such studies are beyond the scope of this report.     

A weighted measure for the similarity analysis of DNA sequences

Here we propose a weighted measure for the similarity analysis of DNA sequences. It is based on LZ complexity and (0,1) characteristic sequences of DNA sequences. This weighted measure enables biologists to extract similarity information from biological sequences according to their requirements. For example, by this weighted measure, one can obtain either the full similarity information or a similarity analysis from a given biological aspect. Moreover, the length of DNA sequence is not problematic. The application of the weighted measure to the similarity analysis of β-globin genes from nine species shows its flexibility.     

Characteristics of soil fauna community in the Dongjiao coconut plantation ecosystem in Hainan, China

From July 2002 to July 2004, we investigated the soil fauna in the Dongjiao coconut plantation of Hainan Island. A total of 5,378 specimens were obtained. These species represented 4 phyla, 12 classes and 27 genera. The soil animal clusters are as follows: All major animal groups presented in a tropical rain forest are present in Dongjiao coconut plantation's soil community. The Dongjiao coconut plantation demonstrated the typical characteristics of a tropical soil animal community. The number of species and the diversity index (H) increase from the high latitude areas toward the equator; The dominance index (C) decreased from the high latitude areas toward the equator; the ratio of Acarina/Collembola and the percentage occurrence of increase along the declination of latitude; the occurrence of termites, a typical member of tropical community, varies from absent to present and a dramatic increase along the declination of latitude. Compared to primary tropical rainforests, the Dongjian coconut plantation community is relatively low in species diversity, and has a high dominance index and low diversity index. This may partially due to some characters of the plantation: singular tree species, monsoon climate, seashore location, high pH and salinity of the soil, soil moisture and other environmental factors. A seasonal change occurs in the community structure but is not obvious. Soil around human residences has a higher organic material content, and has higher counts of specimens and numbers of species, compared to other three sampling sites. Coconut production at locations around human residences is higher than at any other microhabitats. The high production is positively correlated with the richness of animal community in the soil.     

马铃薯吸附干燥特性及模型拟合

在实验的基础上研究了马铃薯的干燥特性,考察了温度、湿度、风速和粒子大小对马铃薯的干燥特性的影响。对马铃薯的干燥特性进行模型拟合结果表明:马铃薯吸附式低温干燥过程分为调整、恒速和降速3个阶段。含水率从76.88%变化到65.44%这段为调整阶段;从65.44%变化到47.34%这段为恒速阶段,主要是去除非结合水;从47.34%降到27.7%为降速阶段,主要脱除非结合水及部分结合水,干燥速度下降。温度是影响干燥的主要因素,温度越高,湿度越小,介质流速快,干燥速度越快。Page模型能够较好地拟合马铃薯干燥过程的特性。     

低温弱光胁迫及恢复对切花菊光合作用和叶绿素荧光参数的影响

以切花菊品种‘神马’为试材,在偏低温弱光(16℃/12℃,PFD100μmol.m-2.s-1)和临界低温弱光(12℃/8℃,PFD60μmol.m-2.s-1)下分别胁迫11d,然后转入正常条件(22℃/18℃,PFD450μmol.m-2.s-1)恢复11d,研究不同低温弱光强度及恢复对菊花光合作用和叶绿素荧光参数的影响.结果表明:低温弱光导致菊花叶片的净光合速率(Pn)和气孔限制值(Ls)下降,而胞间CO2浓度(Ci)上升.偏低温弱光胁迫下菊花叶片暗适应下最大光化学效率(Fv/Fm)和初始荧光(Fo)无明显变化,但光适应下最大光化学效率(Fv′/Fm′)在处理前期略有下降,后期则有所回升;而临界低温弱光处理的Fo明显升高,Fv/Fm和Fv′/Fm′显著降低.PSⅡ光合电子传递量子效率(ΦPSⅡ)、光化学猝灭系数(qP)和表观光合电子传递速率(ETR)均随着低温弱光胁迫程度的增加和时间的延长而降低;偏低温弱光处理植株在解除胁迫后能迅速恢复到对照水平,而临界低温弱光处理植株回升速度较慢;同时,低温弱光胁迫下吸收光强用于分配光化学反应部分(Prate)的比例减少,而天线热耗散(Drate)和反应中心的能量耗散(Ex)比例上升,但天线热耗散为过剩光能的主要分配途径.     

基于图论模型的蛋白质结构类型的研究

提出了基于图论模型的H系数分类蛋白质结构为H结型和NH结型的方法.论述了蛋白质结构中序列不相邻的C_α原子之间的空间距离与序列相邻的C_α原子之间空间距离的关系.用此方法对PDB的66个单链蛋白质结构进行分类,结果显示H结型占18.2%.H结在全α型中出现比例较高,在全β型中出现比例较小,所以H结倾向出现在含有α螺旋的蛋白质结构中.