A study of the stability of subcycling algorithms in structural dynamics

Algorithms for explicit integration of structural dynamics problems with multiple time steps (subcycling) are investigated. Only one such algorithm, due to Smolinski and Sleith has proved to be stable in a classical sense. A simplified version of this algorithm that retains its stability is presented. However, as with the original version, it can be shown to sacrifice accuracy to achieve stability. Another algorithm in use is shown to be only statistically stable, in that a probability of stability can be assigned if appropriate time step limits are observed. This probability improves rapidly with the number of degrees of freedom in a finite element model. The stability problems are shown to be a property of the central difference method itself, which is modified to give the subcycling algorithm. A related problem is shown to arise when a constraint equation in time is introduced into a time-continuous space-time finite element model. (C) 1998 Elsevier Science S.A.

The influence of joint position on the dynamics of perception-action coupling

Six right-handed subjects performed rhythmic flexion and extension movements of the index finger in time with an auditory metronome. On each block of trials, the wrist of the response hand was placed in a extended, neutral or flexed position. In the flex-on-the-beat condition, subjects were instructed to coordinate maximum excursion in the direction of finger flexion with each beat of the metronome. In the extend-on-the-beat condition, subjects were instructed to coordinate maximum excursion in the direction of finger extension with each beat of the metronome. The frequency of the metronome was increased from 2.00 Hz to 3.75 Hz in 8 steps (8 s epochs) of 0.25 Hz. During trials prepared in the extend-on-the-beat pattern, all subjects exhibited transitions to either a flex-on-the-beat pattern or to phase wandering as the frequency of pacing was increased. The time at which these transitions occurred was reliably influenced by the position of the wrist. Four subjects exhibited qualitative departures from the flex-on-the-beat pattern at pacing frequencies that were greater than those at which the extend-on-the-beat pattern could be maintained. The lime at which these departures occurred was not influenced by the position of the wrist. These results are discussed with reference to the constraints imposed on the coordination dynamics by the intrinsic properties of the neuromuscular-skeletal system.

Calcium-permeable AMPA receptors mediate long-term potentiation in interneurons in the amygdala

Fear conditioning is a paradigm that has been used as a model for emotional learning in animals'. The cellular correlate of fear conditioning is thought to be associative N-methyl-D-aspartate (NMDA) receptor-dependent synaptic plasticity within the amygdala(1-3). Here we show that glutamatergic synaptic transmission to inhibitory interneurons in the basolateral amygdala is mediated solely by alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors. In contrast to AMPA receptors at inputs to pyramidal neurons, these receptors have an inwardly rectifying current-voltage relationship, indicative of a high permeability to calcium(4 5), Tetanic stimulation of inputs to interneurons caused an immediate and sustained increase in the efficacy of these synapses. This potentiation required a rise in postsynaptic calcium, but was independent of NMDA receptor activation. The potentiation of excitatory inputs to interneurons was reflected as an increase in the amplitude of the GABAA-mediated inhibitory synaptic current in pyramidal neurons. These results demonstrate that excitatory synapses onto interneurons within a fear conditioning circuit show NMDA-receptor independent long-term potentiation. This plasticity might underlie the increased synchronization of activity between neurons in the basolateral amygdala after fear conditioning(6).

Folding, calcium binding, and structural characterization of a concatemer of the first and second ligand-binding modules of the low-density lipoprotein receptor

The ligand-binding domain of the low-density lipoprotein (LDL) receptor is comprised of seven tandemly repeated ligand-binding modules, each being approximately 40 amino acids long and containing six conserved cysteine residues. We have expressed and characterized a concatemer of the first two modules (LB1 and LB2) of the human LDL receptor. Oxidative folding of the recombinant concatemer (rLB(1-2)), in the presence of calcium ions, gave a single dominant isomer with six disulfide bonds. Peptic cleavage of the short Linker region that connects the last cysteine residue of LB1 and the first cysteine residue of LB2 yielded two discrete fragments, thus excluding the presence of intermodule disulfide bonds. The N-terminal module, LB1, reacted with a conformation-specific monoclonal antibody (IgG-C7) made to LB1 in the native LDL receptor. From this, we concluded that the first module was correctly folded, with the same set of disulfide bonds as LB1 of the LDL receptor. The disulfide bond connections of LB2 were identified from mass spectral analysis of fragments formed by digestion of the C-terminal peptic fragment with elastase. These data showed that the disulfide bonds of LB2 connected Cys(I) and Cys(III), Cys(II) and Cys(V), and Cys(IV) and Cys(VI). This pattern is identical to that found for recombinant LB1 and LB2. The concatemer has two high-affinity calcium-binding sites, one per module. An analysis of the secondary chemical shifts of C alpha protons shows that the conformations of LB1 and LB2 in the concatemer are very similar to those of the individual modules, with no evidence for strong interactions between the two modules.

Human vascular to cardiac tissue selectivity of L- and T-type calcium channel antagonists

1 Voltage-operated calcium channel (VOCC) antagonists are effective antihypertensive and antianginal agents but they also depress myocardial contractility. 2 We compared four L-type calcium channel antagonists, felodipine, nifedipine, amlodipine and verapamil and a relatively T-type selective calcium channel antagonist, mibefradil, on human and rat isolated tissue assays to determine their functional vascular to cardiac tissue selectivity (V/C) ratio. 3 The V/C ratio was calculated as the ratio of the IC50 value of the antagonist that reduced (by 50%) submaximally contracted (K+ 62 mM) human small arteries from the aortic vasa vasorum (vascular, V) mounted in a myograph and the IC50 value of the antagonist that reduced (-)-isoprenaline (6 nM) submaximally stimulated human right atrial trabeculae muscle (cardiac, C) mounted in organ chambers. 4 The average pIC(50) Values (-log IC50 M) for the human vascular preparations were felodipine 8.30, nifedipine 7.78, amlodipine 6.64, verapamil 6.26 and mibefradil 6.22. The average pIC(50) values for the cardiac muscle were felodipine 7.21, nifedipine 6.95, verapamil 6.91, amlodipine 5.94, and mibefradil 4.61. 5 The V/C ratio calculated as antilog [pIC(50)V-pIC(50)C] is thus mibefradil 41, felodipine 12, nifedipine 7, amlodipine 5 and verapamil 0.2. 6 In rat small mesenteric arteries the pIC(50) values for the five drugs were similar to the values for human vasa vasorum arteries contracted by K+ 62 mM. However for methoxamine (10 mu M) contraction in the rat arteries the pIC(50) values were lower for felodipine 7.24 and nifedipine 6.23, but similar for verapamil 6.13, amlodipine 6.28 and mibefradil 5.91. 7 In conclusion in the human tissue assays, the putative T-channel antagonist mibefradil shows the highest vascular to cardiac selectivity ratio; some 3 fold higher than the dihydropyridine, felodipine, and some 200 fold more vascular selective than the phenylalkylamine, verapamil. This favourable vascular to cardiac selectivity for mibefradil, from a new chemical class of VOCC antagonist, may be explained by its putative T-channel selectivity.

In vivo dynamics of axon pathfinding in the Drosophila CNS: A time-lapse study of an identified motoneuron

We developed a system for time-lapse observation of identified neurons in the central nervous system (CNS) of the Drosophila embryo. Using this system, we characterize the dynamics of filopodia and axon growth of the motorneuron RP2 as it navigates anteriorly through the CNS and then laterally along the intersegmental nerve (ISN) into the periphery. We find that both axonal extension and turning occur primarily through the process of filopodial dilation. In addition, we used the GAL4-UAS system to express the fusion protein Tau-GFP in a subset of neurons, allowing us to correlate RP2's patterns of growth with a subset of axons in its environment. In particular, we show that RP2's sharp lateral turn is coincident with the nascent ISN. (C) 1998 John Wiley & Sons, Inc.

Thermal diffusion in molecular dynamics simulations of thin film diamond deposition

Molecular dynamics simulations of carbon atom depositions are used to investigate energy diffusion from the impact zone. A modified Stillinger-Weber potential models the carbon interactions for both sp2 and sp3 bonding. Simulations were performed on 50 eV carbon atom depositions onto the (111) surface of a 3.8 x 3.4 x 1.0 nm diamond slab containing 2816 atoms in 11 layers of 256 atoms each. The bottom layer was thermostated to 300 K. At every 100th simulation time step (27 fs), the average local kinetic energy, and hence local temperature, is calculated. To do this the substrate is divided into a set of 15 concentric hemispherical zones, each of thickness one atomic diameter (0.14 nm) and centered on the impact point. A 50-eV incident atom heats the local impact zone above 10 000 K. After the initial large transient (200 fs) the impact zone has cooled below 3000 K, then near 1000 K by 1 ps. Thereafter the temperature profile decays approximately as described by diffusion theory, perturbed by atomic scale fluctuations. A continuum model of classical energy transfer is provided by the traditional thermal diffusion equation. The results show that continuum diffusion theory describes well energy diffusion in low energy atomic deposition processes, at distance and time scales larger than 1.5 nm and 1-2 ps, beyond which the energy decays essentially exponentially. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

Chaotic dynamics of cold atoms in far-off-resonant donut beams

We describe the classical two-dimensional nonlinear dynamics of cold atoms in far-off-resonant donut beams. We show that chaotic dynamics exists there for charge greater than unity, when the intensity of the beam is periodically modulated. The two-dimensional distributions of atoms in the (x,y) plant for charge 2 are simulated. We show that the atoms will accumulate on several ring regions when the system enters a regime of global chaos. [S1063-651X(99)03903-3].

Structure-activity relationships of omega-conotoxins MVIIA, MVIIC and 14 loop splice hybrids at N and P/Q-type calcium channels

The omega-conotoxins are a set of structurally related, four-loop, six cysteine containing peptides, that have a range of selectivities for different subtypes of the voltage-sensitive calcium channel (VSCC). To investigate the basis of the selectivity displayed by these peptides, we have studied the binding affinities of two naturally occurring omega-conotoxins, MVIIA and MVIIC and a series of 14 MVIIA/MVIIC loop hybrids using radioligand binding assays for N and P/Q-type Ca2+ channels in rat brain tissue. A selectivity profile was developed from the ratio of relative potencies at N-type VSCCs (using [I-125]GVIA radioligand binding assays) and P/Q-type VSCCs (using [I-125]MVIIC radioligand binding assays). in these peptides, loops 2 and 4 make the greatest contribution to VSCC subtype selectivity, while the effects of loops 1 and 3 are negligible. Peptides with homogenous combinations of loop 2 and 4 display clear selectivity preferences, while those with heterogeneous combinations of loops 2 and 4 are less discriminatory. H-1 NMR spectroscopy revealed that the global folds of MVIIA, MVIIC and the 14 loop hybrid peptides were similar; however, several differences in local structure were identified. Based on the binding data and the 3D structures of MVIIA, GVIA and MVIIC, we have developed a preliminary pharmacophore based on the omega-conotoxin residues most Likely to interact with the N-type VSCC. (C) 1999 Academic Press.

The plasma membrane calcium pump, its role and regulation: New complexities and possibilities

Significant progress has been achieved in elucidating the role of the plasma membrane Ca2+-ATPase in cellular Ca2+ homeostasis and physiology since the enzyme was first purified and physiology since the enzyme was first purified and cloned a number of years ago. The simple notion that the PM Ca2+-ATPase controls resting levels of [Ca2+](CYT) has been challenged by the complexity arising from the finding of four major isoforms and splice variants of the Ca2+ pump, and the finding that these are differentially localized in various organs and subcellular regions. Furthermore, the isoforms exhibit differential sensitivities to Ca2+, calmodulin, ATP, and kinase-mediated phosphorylation. The latter pathways of regulation can give rise to activation or inhibition of the Ca2+ pump activity, depending on the kinase and the particular Ca2+ pump isoform. Significant progress is being made in elucidating subtle and more profound roles of the PM Ca2+-ATPase in the control of cellular function. Further understanding of these roles awaits new studies in both transfected cells and intact organelles, a process that will be greatly aided by the development of new and selective Ca2+ pump inhibitors. (C) 1999 Elsevier Science Inc.

Conformations of cyclic octapeptides and the influence of heterocyclic ring constraints upon calcium binding

A comparison is made between the structures and calcium binding properties of four cyclic octapeptides that differ in the number of heterocyclic thiazole and oxazoline ring constraints. The conformations of the naturally occurring cyclic octapeptides ascidiacyclamide 1 and patellamide D 2, which each contain two oxazoline and two thiazole rings, are compared by H-1 NMR spectroscopy with the analogues cyclo(Thr-D-Val(Thz)-Ile)(2) 3 with just two thiazoles, and cyclo(Thr-D-Val-alpha Abu-Ile)(2) 4, with no 5-membered rings. The conformations observed in the solid state for ascidiacyclamide (saddle) and patellamide D (twisted figure of eight) were retained in solution, whilst peptide 3 was found to have a chair shape and peptide 4 displayed a range of conformations. The solid state structure of 4 revealed that the peptide takes a relatively planar conformation with a number of transannular hydrogen bonds, which are apparently retained in solution. Complexation studies utilising H-1 NMR and CD spectroscopy yielded 1∶1 calcium-peptide binding constants (log K) for the four peptides (2.9 (1), 2.8 (2), 4.0 (3) and 5.5 (4)) as well as a 1 : 2 metal-peptide binding constant for 3 (log K = 4.5). The affinity for Ca2+ thus decreases with increasing number of 5-membered ring constraints in the macrocycle (4 > 3 > 2 approximate to 1).

Calcium absorption in postmenopausal osteoporosis: Benefit of HRT plus calcitriol, but not HRT alone, in both malabsorbers and normal absorbers

In a randomized trial involving 71 postmenopausal osteoporotic women with vertebral compression fractures, radiocalcium absorption studies using the Ca-45 single isotope method (alpha) were performed at baseline and after 8 months of treatment with either continuous combined hormone replacement therapy (HRT, as piperazine estrone sulfate 0.625-0.937mg daily +/- medroxyprogesterone acetate 2.5 mg daily depending on uterine status) or HRT plus calcitriol 0.25 mu g twice daily. A calcium supplement of 600 mg nocte was given to only those women who had a daily calcium intake of less than 1 g per day at baseline, as assessed by recalled dietary intake. There was a significant decrease 0.74 (+/- 0.35 SD) to 0.58 (+/- 0.22), Delta alpha = -0.17 (+/- 0.26), p<0.0005] in alpha at 8 months compared with baseline in the HRT-treated group, but a significant increase [0.68 (+/- 0.31) to 0.84 (+/- 0.27), Delta alpha = +0.16 (+/- 0.30), p<0.003] in the HRT-plus-calcitriol treated patients, resulting in alpha being significantly higher after 8 months in the latter group than in the HRT-only group. Although 72% of the patients had been supplemented with calcium between the first and second studies, separate analyses revealed that the change in calcium intake had not affected the result. Further breakdown of the groups into baseline 'normal' absorbers (alpha greater than or equal to 0.55) and 'malabsorbers' (alpha <0.55) revealed that alpha decreased with HRT treatment only in the normal absorbers, and remained stable in the malabsorbers. Conversely, following HRT plus calcitriol treatment, alpha increased only in the malabsorbers, the normal absorbers in this group remaining unchanged. In conclusion, our data show that HRT, of the type and dose used in this study, did not produce an increase in absorption efficiency; it was in fact associated with a fall. increased absorption efficiency cannot be achieved unless calcitriol is used concurrently, and then only in patients with malabsorption. Calcitriol also had a significant effect in normal absorbers in that it prevented the decline in alpha seen with HRT alone, and thus should be considered in all patients with postmenopausal osteoporosis treated with HRT.

Neuromuscular-skeletal constraints upon the dynamics of unimanual and bimanual coordination

In the first of three experiments, 11 participants generated pronation and supination movements of the forearm, in time with an auditory metronome. The metronome frequency was increased in eight steps (0.25 Hz) from a base frequency of 1.75 Hz. On alternating trials, participants were required to coordinate either maximum pronation or maximum supination with each beat of the metronome. In each block of trials, the axis of rotation was either coincident with the long axis of the forearm, above this axis, or below this axis. The stability of the pronate-on-the-beat pattern, as indexed by the number of pattern changes, and the time of onset of pattern change, was greatest when the axis of rotation of the movement was below the long axis of the forearm. In contrast, the stability of the supinate-on-the-beat pattern was greatest when the axis of rotation of the movement was above the long axis of the forearm. In a second experiment, we examined how changes in the position of the axis of rotation alter the activation patterns of muscles that contribute to pronation and supination of the forearm. Variations in the relative dominance of the pronation and supination phases of the movement cycle across conditions were accounted for primarily by changes in the activation profile of flexor carpi radialis (FCR) and extensor carpi radialis longus (ECR). In the Final experiment we examined how these constraints impact upon the stability of bimanual coordination. Thirty-two participants were assigned at random to one of four conditions, each of which combined an axis of rotation configuration (bottom or top) for each limb. The participants generated both inphase (both limbs pronating simultaneously, and supinating simultaneously) and antiphase (left limb pronating and right limb supinating simultaneously, and vice versa) patterns of coordination. When the position of the axis of rotation was equivalent for the left and the right limb, transitions from antiphase to inphase patterns of coordination were Frequently observed. In marked contrast, when the position of the axis of rotation for the left and right limb was contradistinct, transitions From inphase to antiphase patterns of coordination occurred. The results demonstrated that when movements are performed in an appropriate mechanical context, inphase patterns of coordination are less stable than antiphase patterns.