Calcium-independent contraction in lysed cell models of teleost retinal cones: activation by unregulated myosin light chain kinase or high magnesium and loss of cAMP inhibition.

The retinal cones of teleost fish contract at dawn and elongate at dusk. We have previously reported that we can selectively induce detergent-lysed models of cones to undergo either reactivated contraction or reactivated elongation, with rates and morphology comparable to those observed in vivo. Reactivated contraction is ATP dependent, activated by Ca2+, and inhibited by cAMP. In addition, reactivated cone contraction exhibits several properties that suggest that myosin phosphorylation plays a role in mediating Ca2+-activation (Porrello, K., and B. Burnside, 1984, J. Cell Biol., 98:2230-2238). We report here that lysed cone models can be induced to contract in the absence of Ca2+ by incubation with trypsin-digested, unregulated myosin light chain kinase (MLCK) obtained from smooth muscle. This observation provides further evidence that MLCK plays a role in regulating cone contraction. We also report here that lysed cone models can be induced to contract in the absence of Ca2+ by incubation with high concentrations of MgCl2 (10-20 mM). Mg2+-induced reactivated contraction is supported by inosine triphosphate (ITP) just as well as by ATP. Because ITP will not serve as a substrate for MLCK, this finding suggests that Mg2+-activation of contraction does not require myosin phosphorylation. Although Ca2+-induced contraction is completely blocked by cAMP at concentrations less than 10 microM, cAMP has no effect on cone contraction activated by unregulated MLCK or by high Mg2+ in the absence of Ca2+. Because trypsin digestion of MLCK cleaves off not only the Ca2+/calmodulin-binding site but also the site phosphorylated by cAMP-dependent protein kinase, and because Mg2+ activation of cone contraction circumvents MLCK action altogether, both these observations would be expected if cAMP inhibits reactivated cone contraction by catalyzing the phosphorylation of MLCK and thus reducing its affinity for Ca2+, as has been described for smooth muscle. Together our results suggest that in lysed cone models, myosin phosphorylation is sufficient for activating cone contraction, even in the absence of other Ca2+-mediated events, that cAMP inhibition of contraction is mediated by cAMP-dependent phosphorylation of MLCK, and that 10-20 mM Mg2+ can activate actin-myosin interaction to produce contraction in the absence of myosin phosphorylation

The paramagnetic and glass transitions in sudoku

We study the statistical mechanics of a model glassy system based on a familiar and popular mathematical puzzle. Sudoku puzzles provide a very rare example of a class of frustrated systems with a unique groundstate without symmetry. Here, the puzzle is recast as thermodynamic system where the number of violated rules defines the energy. We use Monte Carlo simulation to show that the "Sudoku Hamiltonian" exhibits two transitions as a function of temperature, a paramagnetic and a glass transition. Of these, the intermediate condensed phase is the only one which visits the ground state (i.e. it solves the puzzle, though this is not the purpose of the study). Both transitions are associated with an entropy change, paramagnetism measured from the dynamics of the Monte Carlo run, showing a peak in specific heat, while the residual glass entropy is determined by finding multiple instances of the glass by repeated annealing. There are relatively few such simple models for frustrated or glassy systems which exhibit both ordering and glass transitions, sudoku puzzles are unique for the ease with which they can be obtained with the proof of the existence of a unique ground state via the satisfiability of all constraints. Simulations suggest that in the glass phase there is an increase in information entropy with lowering temperature. In fact, we have shown that sudoku have the type of rugged energy landscape with multiple minima which typifies glasses in many physical systems, and this puzzling result is a manifestation of the paradox of the residual glass entropy. These readily-available puzzles can now be used as solvable model Hamiltonian systems for studying the glass transition

The MOLDY short-range molecular dynamics package

We describe a parallelised version of the MOLDY molecular dynamics program. This Fortran code is aimed at systems which may be described by short-range potentials and specifically those which may be addressed with the embedded atom method. This includes a wide range of transition metals and alloys. MOLDY provides a range of options in terms of the molecular dynamics ensemble used and the boundary conditions which may be applied. A number of standard potentials are provided, and the modular structure of the code allows new potentials to be added easily. The code is parallelised using OpenMP and can therefore be run on shared memory systems, including modern multicore processors. Particular attention is paid to the updates required in the main force loop, where synchronisation is often required in OpenMP implementations of molecular dynamics. We examine the performance of the parallel code in detail and give some examples of applications to realistic problems, including the dynamic compression of copper and carbon migration in an iron-carbon alloy

Ab initio studies of structural instabilities in magnesium silicate perovskite

Density-functional simulations are used to calculate structural properties and high-symmetry phonons of the hypothetical cubic phase, the stable orthorhombic phase and an intermediate tetragonal phase of magnesium silicate perovskite. We show that the structure of the stable phase is well described by freezing in a small number of unstable phonons into the cubic phase. We use the frequencies of these unstable modes to estimate transition temperatures for cubic--tetragonal and tetragonal--orthorhombic phase transitions. These are investigated further to find that the coupling with the strain suggests that phonons give a better representation than rigid unit modes. The phonons of an intermediate tetragonal phase were found to be stable except for two rotational modes. The eigenvectors of the most unstable mode of each of the cubic and tetragonal phases account for all the positional parameters of the orthorhombic phase. The phase boundary for the orthorhombic--tetragonal transition intersects possible mantle geotherms, suggesting that the tetragonal phase may be present in the lower mantle.Comment: 16 pages, REVTEX, 7 postscript figures (Fig 1 very large, contact Authors if required); submitted to Physics and Chemistry of Mineral