Torsional path integral Monte Carlo method for calculating the absolute quantum free energy of large molecules

A new technique for evaluating the absolute free energy of large molecules is presented. Quantum-mechanical contributions to the intramolecular torsions are included via the torsional path integral Monte Carlo (TPIMC) technique. Importance sampling schemes based on uncoupled free rotors and harmonic oscillators facilitate the use of the TPIMC technique for the direct evaluation of quantum partition functions. Absolute free energies are calculated for the molecules ethane, n-butane, n-octane, and enkephalin, and quantum contributions are found to be significant. Comparison of the TPIMC technique with the harmonic oscillator approximation and a variational technique is performed for the ethane molecule. For all molecules, the quantum contributions to free energy are found to be significant but slightly smaller than the quantum contributions to internal energy

Quantum free energies of the conformers of glycine on an ab initio potential energy surface

The torsional path integral Monte Carlo (TPIMC) technique is used to study the five lowest-energy conformers of glycine. The theoretical method provides an anharmonic and quantum-mechanical description of conformational free energy and is used for the first time with an ab initio potential energy surface. The 3-dimensional torsional potential energy surface of glycine was obtained at the MP2/6-311++G** level of theory and is optimized with respect to the non-torsional degrees of freedom. Calculated conformer populations compare well with those reported in recent matrix-isolation infrared spectroscopy experiments. An additional conformer, not yet observed, is predicted to be heavily populated in the thermal equilibria probed by experiment, and a new explanation for its elusiveness is provided. Quantum effects, such as zero point energy, are found to substantially alter conformer populations, and an algorithm for estimating the role of non-torsional vibrations in the conformational thermodynamics of a molecule is introduced

Immunocompetence in Hydra. Epithelial cells recognize self-nonself and react against it

The evolution of effective immunologic defense mechanisms in multicellular organisms involves the ability of host cells to distinguish betweeen self and nonself and to react appropriately to eliminate foreign tissue. By producing interspecies grafts we have obtained evidence that immunorecognition followed by incompatibility reactions occur in Hydra. Our results demonstrate that epithelial cells of Hydra recognize and phagocytose foreign hydra cells, indicating that they are the effector cells in the incompatibility reactions. This observation is consistent with the idea that immunocompetence appeared early in the evolution of multicellular organisms

Torsional path integral Monte Carlo method for the quantum simulation of large molecules

A molecular application is introduced for calculating quantum statistical mechanical expectation values of large molecules at nonzero temperatures. The Torsional Path Integral Monte Carlo (TPIMC) technique applies an uncoupled winding number formalism to the torsional degrees of freedom in molecular systems. The internal energy of the molecules ethane, n-butane, n-octane, and enkephalin are calculated at standard temperature using the TPIMC technique and compared to the expectation values obtained using the harmonic oscillator approximation and a variational technique. All studied molecules exhibited significant quantum mechanical contributions to their internal energy expectation values according to the TPIMC technique. The harmonic oscillator approximation approach to calculating the internal energy performs well for the molecules presented in this study but is limited by its neglect of both anharmonicity effects and the potential coupling of intramolecular torsion

Cloned interstitial stem cells grow as contiguous patches in hydra

The migration of interstitial cells was analyzed during the growth of stem cell clones in vivo. The spatial distribution of cloned cells was analyzed at a time by which extensive migration of interstitial cells could have occurred. All interstitial cell clones were found to form large contiguous patches of cells. The results indicate that there is little migration of large interstitial cells in undisturbed tissue during normal growth. This finding is surprising since numerous grafting experiments have shown extensive migration of these cells. The implications of finding nonrandomly distributed stem cells are discussed

Male and female stem cells and sex reversal in Hydra polyps

Single interstitial stem cells of male polyps of Hydra magnipapillata give rise to clones that differentiate either male or female gametes. To test the sexual stability of these clones, stem cells were recloned. The results indicate that stem cells from female clones are stable in their sexual differentiation capacity; male stem cells, by comparison, switch sexual phenotype at the rate of 10-2 per cell per generation. As a result, female polyps contain only female stem cells; male polyps contain a mixture of male and female stem cells. A model is presented in which the sexual phenotype of Hydra polyps is controlled by (i) the switching rate of male and female stem cells and (ii) the repression of female differentiation by male stem cells

Torsional anharmonicity in the conformational thermodynamics of flexible molecules

We present an algorithm for calculating the conformational thermodynamics of large, flexible molecules that combines ab initio electronic structure theory calculations with a torsional path integral Monte Carlo (TPIMC) simulation. The new algorithm overcomes the previous limitations of the TPIMC method by including the thermodynamic contributions of non-torsional vibrational modes and by affordably incorporating the ab initio calculation of conformer electronic energies, and it improves the conventional ab initio treatment of conformational thermodynamics by accounting for the anharmonicity of the torsional modes. Using previously published ab initio results and new TPIMC calculations, we apply the algorithm to the conformers of the adrenaline molecule

Growth regulation in Hydra. Relationship between epithelial cell cycle length and growth rate

The relationship between epithelial cell production and growth rate was investigated in Hydra attenuata under different feeding regimes. The increase of epithelial cell number was compared to the duration of the epithelial cell cycle using standard methods of cell cycle analysis. The results indicate that cell cycle changes accompanying changes in feeding regime are not sufficient to explain the altered growth rate. Under heavy feeding regimes, epithelial cell production equals tissue growth rate. At low feeding level or under starvation conditions the epithelial cell cycle lengthens and growth rate of epithelial cell population is slowed. However, the cell cycle changes are insufficient to account for the reduction in tissue growth and thus there is an effective overproduction of epithelial cells amounting to 10% per day. Evidence suggests that these excess cells are phagocytized by neighboring cells in the tissue. Thus phagocytosis is directly or indirectly involved in regulating the growth of hydra tissue

Reliability analysis of continuous fiber composite laminates

A composite lamina may be viewed as a homogeneous solid whose directional strengths are random variables. Calculation of the lamina reliability under a multi-axial stress state can be approached by either assuming that the strengths act separately (modal or independent action), or that they interact through a quadratic interaction criterion. The independent action reliability may be calculated in closed form, while interactive criteria require simulations; there is currently insufficient data to make a final determination of preference between them. Using independent action for illustration purposes, the lamina reliability may be plotted in either stress space or in a non-dimensional representation. For the typical laminated plate structure, the individual lamina reliabilities may be combined in order to produce formal upper and lower bounds of reliability for the laminate, similar in nature to the bounds on properties produced from variational elastic methods. These bounds are illustrated for a (0/plus or minus 15)sub s Graphite/Epoxy (GR/EP) laminate. And addition, simple physically plausible phenomenological rules are proposed for redistribution of load after a lamina has failed. These rules are illustrated by application to (0/plus or minus 15)sub s and (90/plus or minus 45/0)sub s GR/EP laminates and results are compared with respect to the proposed bounds