Flocculation in brewing yeasts: a computer simulation study

A new simulator, INDISIM-FLOC, based on the individual-based simulator INDISIM, is used to examine the predictions of two different models of yeast flocculation. The first, proposed by Calleja is known as the "addition" model. The second, proposed by Stratford is known as the "cascade" model. The simulations show that the latter exhibits a better qualitative agreement with available experimental data.     

肿瘤坏死因子与其受体相互作用的计算机模拟研究

利用同源模建方法,以肿瘤坏死因子(TNF)R55受体胞外区的晶体结构为参考模板,预测了TNFR75受体胞外区Cys18-Phe147片段的三维结构。根据R55受体胞外区与LT相结合的复合物的晶体结构,预测了TNF与R55及R75胞外区的复合物的三维结构,模拟了TNF与受体之间的相互作用。由于TNF与受体的作用形式是三聚体对三聚体,因此在模拟TNF与受体相互作用时选择了包括一个非对称的TNF三聚体和一个受体(R55或R75)单体的模拟系统。结合已有的突变体实验结果,利用计算机分析手段,发现了一些TNF突变体之所以具有受体选择性的三维结构基础和发挥了关键作用的氨基酸残基以及这些残基之间的主要作用形式。研究深化了对已有的突变体实验结果的认识,建立了不同的实验结果之间的内在关联,为以后有目的的新型突变体设计和实验研究打下了基础。     

A mathematical model and computer simulation study of insulin receptor regulation

A homeomorphic mathematical model of cell surface insulin receptor regulation is developed. The overall structure of the model is based on molecular mechanisms suggested by in vivo and in vitro experimental evidence from many different cell types. Model parameters correspond to cellular processes which are constrained by known boundry value conditions. As an example, computer simulation results are compared with published data from BC3H-1 myocytes in culture. With appropriate parameter choice, this model is able to simulate data from other cell types. Cellular processes which are explicitly represented in the model include: bound and unbound receptor endocytosis, receptor recycling, intracellular receptor degradation, and state-dependent receptor synthesis. Most of these processes are represented as first-order events. Using more complex representations of the model structure with higher order rate constants or saturable pathways does not qualitatively improve simulation results. Simulations are able to reproduce ligand-induced down and up regulation of receptors as well as the initial spontaneous display of surface insulin receptors. To demonstrate the behavior of our model and illustrate its utility for explaining insulin receptor regulation for a variety of conditions, simulations for which experimental data is unavailable for direct comparison are also shown. We believe the structure of our model is sufficient to explain insulin receptor regulation in a wide variety of cell types. In addition our model may aid in understanding the receptor component of insulin resistance (decreased sensitivity or responsiveness to insulin) seen in pathological states such as obesity and diabetes mellitus. Finally, this model may be applicable to the study of the regulation of other polypeptide hormone receptors.     

Solvation free energies of glutamate and its metal complexes: A computer simulation study

Most biological ion channels demonstrate a high degree of selectivity for one type of ion more than others, and in many cases, how they control attaining this is still not clear. So we have studied on some metal ion compounds of glutamate. The Glutamate and its meal ion compounds (Ca²⁺, Na⁺, K⁺ and Li⁺) were first modeled by ab initio calculations and then Monte Carlo simulation was used to calculate solvation free energies and also the complexes free energies for the related structures. The results indicated that Glutamate-Ca²⁺ have more stability in water than other metal ion. Also, it was found out that the more movement in ions; less stability of the structure would result. This trend can be seen both in gas and liquid phase.     

Reaction kinetics and mechanisms of penicillin amidase: A comparative study by computer simulation

Based on experimental results, two different kinetic models and reaction mechanisms of penicillin amidase (penicillin amidohydrolase EC 3.5.1.11) have been studied and analysed. The enzyme from Escherichia coli shows an ordered uni-bi reaction mechanism while the enzyme from Bacillus megaterium shows kinetics of double inhibition by the reaction products. The difference in the reaction mechanism is elucidated by two possible mechanistic models on a theoretical basis. Also suggested is the analytical method by which two different reaction mechanisms can be tested and confirmed.     

Statistical comparison of power spectra: A computer program

A general computer program, of interest to neurobiologists, for comparison of uncorrelated power spectra is presented. The program, ALPHAF, is written in Fortran V language and is easily compiled and executed in any computer with a Fortran software package. The purpose of the program is to read power spectral density (PSD) values into a computer from digital tape and then perform the following functions: (1) compute and report power; (2) compute and report averages of the computed power values over selected frequency bands and (3) compute and report statistical standardised normal distributions for a frequency spectrum of interest. The standardised normal distribution statistic is derived by using equation 7.20, page 250 of Random data: Analysis and measurement procedures (Bendat and Piersol, 1971).     

Infanticide and fertility among Eskimos: A computer simulation

Until recently, certain Eskimo groups were reported to practice female infanticide in the belief that the time spent suckling a girl would delay the mother's next opportunity to bear a son, males being preferred to females because of their future role as providers in a hunting economy. From sex ratios in census data, rates of female infanticide of up to 66% for some groups have been inferred, leading some ethnographers to conclude that these groups were headed for extinction. Eskimo beliefs regarding the effects of infanticide on fertility, however, are in accord with the results of research on the relation of fertility and lactation: The cessation of lactation following infanticide would significantly shortern the expected interval until the next birth. Given this fact and available field data regarding the parameters of Eskimo population growth, the present computer simulation indicates that Eskimo populations could sustain a rate of 30% female infanticide and still survive. Higher reported rates are explained as the combined result of female infanticide plus the tendency of ethnographers to overestimate to overestimate the ages of juvenile females relative to juvenile males.     

The effect of amino acid substitution on protein-folding and -unfolding transition studied by computer simulation

Protein-folding and -unfolding transitions were studied by the method of computer simulation. The protein was modeled as a two-dimensional lattice polymer. Various energy terms were assumed to be operative between units composing the polymer. But hydrophobic interactions were neglected explicitly. Both thermodynamic and kinetic quantities were obtained from the simulation, and from their temperature dependence in the transition zone characteristics of the conformational transition of proteins were discussed. Two amino acid substituted models, differing in the location of substitution, were studied and compared with the original in order to clarify the effect of substitution on conformational transition of proteins. The following conclusions were reached in this study: (1) The relaxation time of the slow mode, which reflects the overall folding and unfolding processes, shows a peak near the transition temperature, while that of the fast mode is almost independent of temperature. The peak of the slow mode occurs at a slightly lower temperature than the transition temperature. (2) The dependence of the logarithm of the rate constants on the inverse of temperature (Arrhenius plot) is linear. Therefore, the plot of the free energy of activation vs temperature is linear. (3) The values of kinetic parameters obtained suggest that in the activated state the intramolecular interactions are half broken, while the state is close to the native state on the entropy axis. (4) The amino acid substitution, which is modeled as having slightly unfavorable short-range interactions, causes the substituted ones to be slightly unstable. Moreover, it causes the folding transition to slow. From the analysis of the way slowing down is observed in the two substituted models, we conclude that a structure, designed to model a β-sheet, is formed before it gets assembled with other structures, which are designed to model α-helices. The process of assembly occurs nearly at the activated state of the folding and unfolding transition. (5) It is suggested from this study that the maximum of folding rate constant in the Arrhenius plot that has been observed experimentally in real proteins is likely due to hydrophobic interactions.     

Differential growth and plant tropisms: a study assisted by computer simulation

Tropisms and other movements of a plant organ result from alterations in local rates of cell elongation and a consequent development of a growth differential between its opposite sides. Relative elemental rates of elongation (RELELs) are useful to characterize the pattern of growth along and round an organ. We assume that the value of the RELEL at a given point is dependent on distance from the tip and that the distribution of values along the organ surface can be characterized in terms of the spread and the position of the maximum value. A computer model is described which accommodates these parameters and simulates tropic curvatures due to differential growth. Additional regulatory functions help to return the simulated organ to its original orientation. Particular attention is given to the simulation of root gravitropism because here not only do each of the various growth and regulatory parameters have a known biological counterpart, but some can also be given an actual quantitative value. The growth characteristics relate to the biophysical properties of cells in the elongation zone of the root, while the regulatory functions relate to aspects of the graviperception and transmission systems. We believe that, given a suitably flexible model, computer simulation is a powerful means of characterizing, in a quantitative way, the contribution of each parameter to the elongation of plant organs in general and their tropisms in particular.     

A computer simulation study of the hydrated proton in a synthetic proton channel

Classical molecular dynamics simulations using the multistate empirical valence bond model for aqueous proton transport were performed to characterize the hydration structure of an excess proton inside a leucine-serine synthetic ion channel, LS2. For such a nonuniform pore size ion channel, it is found that the Zundel ion (H(5)O(2)(+)) solvation structure is generally more stable in narrow channel regions than in wider channel regions, which is in agreement with a recent study on idealized hydrophobic proton channels. However, considerable diversity in the relative stability of the Zundel to Eigen cation (H(9)O(4)(+)) was observed. Three of the five wide channel regions, one located at the channel's center and the other two located near the channel mouths, are found to show extraordinary preference for the Eigen solvation structure. This implies that proton hopping is inhibited in these regions and therefore suggests that these regions may behave as barriers in the proton conducting pathway inside the channel. The proton solvation is also greatly influenced by the local molecular environment of the protein. In particular, the polar side chains of the Ser residues, which are intimately involved in the solvation structure, can greatly influence proton solvation. However, no preference of the influence by the various Ser side chains was found; they can either promote or prevent the formation of certain solvation structures.     

Relationship between the reproductive capability of a subdivided population and its genetic structure: computer simulation, the simplest multilocus model

An original computer model, simulating joint genetic and demographic dynamics of subdivided populations, is proposed. The model accounts for the reverse effect of the genetic structure on the reproductive capability of a population, which is based on a postulated limited set of biallelic loci, controlling variation in an adaptive quantitative trait. The model allows to simulate spreading of the population, originating from a single small colony, which involves genetic and demographic equilibrium (a normal population-genetic process) and reorganization of the genetic structure of the subdivided population upon anthropogenic pressure, associated with a decrease in its reproductive capability (an adverse population-genetic process).     

A mathematical model and computer simulation study of insulin sensitive glucose transporter regulation

A mathematical model of insulin sensitive glucose transporter regulation is developed. Model structure is based on experimental evidence from adipocytes and myocytes. Model parameters correspond with known cellular processes. As an example, computer simulation results are compared with data from rat adipocytes. Cellular processes explicitly represented in the model include state-dependent glucose transporter synthesis and degradation rates, insulin sensitive glucose transporter translocation rates, and a glucose transporter endocytosis rate. Most of these processes are represented as first-order events. Using more complex representations of the model structure (e.g. higher order rate constants or saturable pathways) or alternative structures did not result in qualitatively better results. The model is able to accurately simulate the insulin sensitive, insulin concentration dependent, reversible translocation of glucose transporters observed in normal adipocytes. The model is also able to accurately simulate the changes in regulation of glucose transporter translocation observed with increases in cell surface area. Finally, the model can simulate pathogenic states which induce impairment of glucose transporter regulation (e.g. altered glucose transporter regulation in adipocytes from rats on high fat diets, rats with streptozotocin induced diabetes, and fasted rats). Since the structure of our model is sufficient to explain glucose transporter regulation in both normal and pathological states, it may aid in understanding the post-receptor components of insulin resistance (decreased sensitivity or responsiveness to insulin) seen in pathological states such as obesity and diabetes mellitus.     

Comparison between effects of various partial inferior turbinectomy options on nasal airflow: a computer simulation study

Partial inferior turbinectomy is typically performed on patients suffering from chronic nasal obstruction due to hypertrophy of inferior turbinates and is refractory to other more conservative treatments. The effects of the various options of incision performed on the inferior turbinate in terms of the resulting nasal airflow pattern are examined using computational fluid mechanics. The pressure drops across the severely blocked nose and healthy nose models were found to be 32.3 and 12.3 Pa, respectively, whereas the pressure drops across the nasal cavity following one-third turbinate resection, total turbinate resection and front-end resection were obtained as 5.8, 6.1 and 30.5 Pa correspondingly. Based on the total pressure drop results, the one-third resection option seems to be better than the front-end surgery and the total turbinate resection.     

The structure of percolated polymer systems: a computer simulation study

We studied the percolation process in a system consisting of long flexible polymer chains and solvent molecules. The polymer chains were approximated by linear sequences of beads on a two-dimensional triangular lattice. The system was athermal and the excluded volume was the only potential. The properties of the model system across the entire range of polymer concentrations were determined by Monte Carlo simulations employing a cooperative motion algorithm (CMA). The scaling behavior and the structure of the percolation clusters are presented and discussed.     

Transcription regulation by steroid hormones: a computer simulation study.

Three-dimensional structures of complexes of 66 amino acid-DNA binding domains of human progesterone (hPR), estrogen (hER) and glucocorticoid (hGR) receptors (proteins), with ten base pair DNA duplexes: d(AGGTCATGCT).d(AGCATGACCT) and d(AGAACATGCT).d(AGCATGTTCT) were obtained using computer modeling and molecular mechanics techniques. Cartesian coordinates for the proteins were obtained from: 1) structural data of hER and hGR by NMR spectroscopy; 2) steric constraints imposed by tetrahedral coordination of the zinc ion to Cys residues, and 3) energy minimization in torsional and cartesian space. The proteins were made to interact with DNA (in B-form) in major groove through alpha-helical linker between the two zinc fingers. The geometry of the complexes was obtained by allowing them to slide, glide, penetrate in to and out of the groove, and to rotate about the helical axis. The complexes were energy minimized. Also maximized was the number of H-bonds between proteins and DNA. The complex structures were refined by molecular mechanics using AMBER 3.0. Structural parameters of DNA were analyzed in each complex and compared with those of native DNA optimized separately. The stereochemical differences of the complexes are discussed.     

Reciprocal interaction between IK1 and If in biological pacemakers: A simulation study

Pacemaking dysfunction (PD) may result in heart rhythm disorders, syncope or even death. Current treatment of PD using implanted electronic pacemakers has some limitations, such as finite battery life and the risk of repeated surgery. As such, the biological pacemaker has been proposed as a potential alternative to the electronic pacemaker for PD treatment. Experimentally and computationally, it has been shown that bio-engineered pacemaker cells can be generated from non-rhythmic ventricular myocytes (VMs) by knocking out genes related to the inward rectifier potassium channel current (I_K1) or by overexpressing hyperpolarization-activated cyclic nucleotide gated channel genes responsible for the “funny” current (I_f). However, it is unclear if a bio-engineered pacemaker based on the modification of I_K1- and I_f-related channels simultaneously would enhance the ability and stability of bio-engineered pacemaking action potentials. In this study, the possible mechanism(s) responsible for VMs to generate spontaneous pacemaking activity by regulating I_K1 and I_f density were investigated by a computational approach. Our results showed that there was a reciprocal interaction between I_K1 and I_f in ventricular pacemaker model. The effect of I_K1 depression on generating ventricular pacemaker was mono-phasic while that of I_f augmentation was bi-phasic. A moderate increase of I_f promoted pacemaking activity but excessive increase of I_f resulted in a slowdown in the pacemaking rate and even an unstable pacemaking state. The dedicated interplay between I_K1 and I_f in generating stable pacemaking and dysrhythmias was evaluated. Finally, a theoretical analysis in the I_K1/I_f parameter space for generating pacemaking action potentials in different states was provided. In conclusion, to the best of our knowledge, this study provides a wide theoretical insight into understandings for generating stable and robust pacemaker cells from non-pacemaking VMs by the interplay of I_K1 and I_f, which may be helpful in designing engineered biological pacemakers for application purposes.