A stochastic approximation algorithm with Markov chain Monte-Carlo method for incomplete data estimation problems

We propose a general procedure for solving incomplete data estimation problems. The procedure can be used to find the maximum likelihood estimate or to solve estimating equations in difficult cases such as estimation with the censored or truncated regression model, the nonlinear structural measurement error model, and the random effects model. The procedure is based on the general principle of stochastic approximation and the Markov chain Monte-Carlo method. Applying the theory on adaptive algorithms, we derive conditions under which the proposed procedure converges. Simulation studies also indicate that the proposed procedure consistently converges to the maximum likelihood estimate for the structural measurement error logistic regression model.

Guanidine hydrochloride blocks a critical step in the propagation of the prion-like determinant [PSI+] of Saccharomyces cerevisiae

The cytoplasmic heritable determinant [PSI+] of the yeast Saccharomyces cerevisiae reflects the prion-like properties of the chromosome-encoded protein Sup35p. This protein is known to be an essential eukaryote polypeptide release factor, namely eRF3. In a [PSI+] background, the prion conformer of Sup35p forms large oligomers, which results in the intracellular depletion of functional release factor and hence inefficient translation termination. We have investigated the process by which the [PSI+] determinant can be efficiently eliminated from strains, by growth in the presence of the protein denaturant guanidine hydrochloride (GuHCl). Strains are “cured” of [PSI+] by millimolar concentrations of GuHCl, well below that normally required for protein denaturation. Here we provide evidence indicating that the elimination of the [PSI+] determinant is not derived from the direct dissolution of self-replicating [PSI+] seeds by GuHCl. Although GuHCl does elicit a moderate stress response, the elimination of [PSI+] is not enhanced by stress, and furthermore, exhibits an absolute requirement for continued cell division. We propose that GuHCl inhibits a critical event in the propagation of the prion conformer and demonstrate that the kinetics of curing by GuHCl fit a random segregation model whereby the heritable [PSI+] element is diluted from a culture, after the total inhibition of prion replication by GuHCl.

An equilibrium statistical model for the spreading phase of open-ocean convection

A “most probable state” equilibrium statistical theory for random distributions of hetons in a closed basin is developed here in the context of two-layer quasigeostrophic models for the spreading phase of open-ocean convection. The theory depends only on bulk conserved quantities such as energy, circulation, and the range of values of potential vorticity in each layer. The simplest theory is formulated for a uniform cooling event over the entire basin that triggers a homogeneous random distribution of convective towers. For a small Rossby deformation radius typical for open-ocean convection sites, the most probable states that arise from this theory strongly resemble the saturated baroclinic states of the spreading phase of convection, with a stabilizing barotropic rim current and localized temperature anomaly.

A quantitative model for membrane fusion based on low-energy intermediates

The energetics of a fusion pathway is considered, starting from the contact site where two apposed membranes each locally protrude (as “nipples”) toward each other. The equilibrium distance between the tips of the two nipples is determined by a balance of physical forces: repulsion caused by hydration and attraction generated by fusion proteins. The energy to create the initial stalk, caused by bending of cis monolayer leaflets, is much less when the stalk forms between nipples rather than parallel flat membranes. The stalk cannot, however, expand by bending deformations alone, because this would necessitate the creation of a hydrophobic void of prohibitively high energy. But small movements of the lipids out of the plane of their monolayers allow transformation of the stalk into a modified stalk. This intermediate, not previously considered, is a low-energy structure that can reconfigure into a fusion pore via an additional intermediate, the prepore. The lipids of this latter structure are oriented as in a fusion pore, but the bilayer is locally compressed. All membrane rearrangements occur in a discrete local region without creation of an extended hemifusion diaphragm. Importantly, all steps of the proposed pathway are energetically feasible.