Regularized numerical integration of multibody dynamics with the generalized alpha method*

This paper discusses the consistent regularization property of the generalized α method when applied as an integrator to an initial value high index and singular differential-algebraic equation model of a multibody system. The regularization comes from within the discretization itself and the discretization remains consistent over the range of values the regularization parameter may take. The regularization involves increase of the smallest singular values of the ill-conditioned Jacobian of the discretization and is different from Baumgarte and similar techniques which tend to be inconsistent for poor choice of regularization parameter. This regularization also helps where pre-conditioning the Jacobian by scaling is of limited effect, for example, when the scleronomic constraints contain multiple closed loops or singular configuration or when high index path constraints are present. The feed-forward control in Kane's equation models is additionally considered in the numerical examples to illustrate the effect of regularization. The discretization presented in this work is adopted to the first order DAE system (unlike the original method which is intended for second order systems) for its A-stability and same order of accuracy for positions and velocities.

Thermal activation parameters for low temperature deformation of alpha-hafnium

Abstract is not available.

Effect of beta-n-oxalyl-L-alpha, beta-diaminopropionic acid on glutamate uptake by synaptosomes

THE unusual amino acid beta-N-oxalyl-L-alpha, beta-diaminopropionic acid (ODAP), isolated from the seeds of Lathyrus sativus is a potent neurotoxin1−3. It produces biochemical changes in the brain typical of an excitant amino acid and is implicated in the aetiology of human neurolathyrism caused by eating the seeds of L. sativus 4−6. It may act as a glutamate antagonist: ODAP inhibits glutamate oxidation7 possibly by inhibiting glutamate uptake in bovine brain mitochondria; it also acts as a competitive inhibitor of glutamate uptake in certain strains of yeast8, and a similar process might occur at the synaptic level. Any effect of ODAP on glutamate uptake at synapses is significant in view of the neurotransmitter function of glutamate, which seems to be neuroexcitory as well as neurotoxic9−12. But Balcar and Johnston13 have shown with rat brain slices that ODAP does not inhibit the glutamate uptake by the high affinity system.

A Case Study of the Conformation of Poly(alpha-aminoisobutyric acid): alpha- or 310-Helix

The relative stabilities of a- and Blo-helical structures for polymers of a-aminoisobutyric acid (Aib) have been worked out, using the classical potential energy functions. To make a comparative study, we have used Buckingham "6-exp" and Kitaigorodsky's potential functions. Conformational analysis of the dipeptide segment with Aib residue indicates the necessity for nonplanar distortion of the peptide unit, which is a common feature in the observed crystal structures with Aib residues. In the range of Aw -10 to +loo studied, a-helical conformations are preferred in the region -3" < Aw < +loo, and Blo-helical conformations are preferred in the region -3" > Aw > -10'. Minimum energy conformations for right-handed structures are found in the +ue region of Aw and correspondingly for left-handed structures in the -ue region of Aw. For Aw - 6", a-helical structures have four- or near fourfold symmetry with h - 1.5 A. Such a helix with n = 4 and h = 1.5 A is termed an a'-helix. This structure is found to be consistent with the electron diffraction data of Malcolm3 and energetically more favorable than the standard 310-helix.

Mapping of assembled epitopic regions of human chorionic gonadotropin reveals proximity of CTP alpha to the determinant loop beta 93-100

Three overlapping assembled epitopes of beta hCG have been mapped using MAb probes and a single step solid phase radioimmunoassay. These epitopes have been shown to be at receptor binding region comprising of the loop region beta Cys93-Cys100. Importance of disulphide bonds in maintaining integrity of these epitopes is assessed. Two MAbs (INN 58 and INN 22) interact with the beta region as well as the alpha C-terminal peptide, while the other MAb INN 24 interacts with only the beta region. Cross-reactivity pattern with beta hCG and hLH as web as the reported crystal structure of hCG substantiates the epitope identification. The results demonstrate utility of MAbs as probes in investigations on three-dimensional structure of gonadatropins.

Characterisation of the DNA-polymerase-alpha-primase complex from the silk glands of Bombyx mori

Silk gland cells ofBombyx mori undergo chromosomal endoduplication throughout larval development. The DNA content of both posterior and middle silk gland nuclei increased by 300000 times the haploid genomic content, amounting to 18 rounds of replication. The DNA doubling time is approximately 48 h and 24 h during the fourth and fifth instars of larval development. However, DNA content does not change during the interim moult. Concomitant with DNA content, DNA polymerase activity also increases as development progressed. Enzyme activity is predominantly due to DNA polymerase with no detectable level of polymerase . DNA polymerase from silk gland extracts was purified to homogeneity (using a series of columns involving ionexchange, gel-filtration and affintiy chromatography), resulting in a 4000-fold increase in specific activity. The enzyme is a heterogeneous multimer of high molecular mass, and the catalytic (polymerase) activity is resident in the 180-kDa subunit. The enzyme shows a PI of 6.2 and theKm values for the dNTP vary over 5-16 . The polymerase is tightly associated with primase activity and initiates primer synthesis in the presence of ribonucleoside triphosphates on a single-stranded DNA template. The primase activity is resident in the 45-kDa subunit. The enzyme is devoid of any detectable exonuclease activity. The abundance of DNA polymerase α in silk glands and its strong association with the nuclear matrix suggest a role in the DNA endoduplication process.

Voltage gated Na+ channels contribute to membrane voltage fluctuation in \alpha T3-1 pituitary gonadotroph cells

\alpha T3-1 cells showed a slope resistance of 1.8 G\omega. The cell membrane surface was not smooth and a scanning electron micrograph showed a complex structure with blebs and microvilli like projections. The cells showed spontaneous fluctuations at zero current resting membrane potential and hyperpolarization increased the amplitude of membrane potential fluctuations. The amplitude of membrane potential fluctuations at hyperpolarized membrane potential was attenuated on application of TTX to the bath solution. The potential at which half steady state inactivation of isolated sodium current occurred, was at a very hyperpolarized potential (-95.4 mV). The study presented in this paper shows that the voltage gated sodium channels contribute to the increase in the amplitude of electrical noise with hyperpolarization in \alpha T3-1 cells.

Crystal structures and thermal analyses of alkyl 2-deoxy-alpha-D-arabino-hexopyranosides

The crystal structures of alkyl 2-deoxy-alpha-D-arabino-hexopyranosides, with the alkyl chain lengths from C-8 to C-18, are established by the single crystal X-ray structural determination. The even-alkyl chain length derivatives crystallized orthorhombic, with space group P2(1)2(1)2(1), whereas the odd-alkyl chain length derivatives crystallized monoclinic, with space group P2(1). The sugar moieties retained a C-4(1) chair conformation and the conformation of the alkyl chains was all-trans. The molecules formed a bilayer structure, in which alkyl chains were interdigitated.The hydrogen bonds, originating from the sugar moieties, were observed in adjacent layers and also within the same layer, resulting in the formation of infinite chains. The alkyl chains arranged parallel to each other and formed planar structures. The thermal properties of the alkyl 2-deoxy glucosides were analyzed further. It was observed that none of the derivatives exhibited mesomorphism. This study establishes that the absence of the hydroxyl group at C-2 of the sugar moiety results in a non-mesogenic nature of the alkyl 2-deoxy-alpha-D-glycosides, as opposed to the profound mesogenic nature of the normal alkyl glycosides.

A Remarkably Simple One-Step Procedure for the Preparation of alpha-Bromo-alpha,beta-Unsaturated Carbonyl Compounds

An easy and convenient one-step procedure for the conversion of alpha,beta-unsaturated carbonyl compounds into their corresponding bromo-enones using NBS-Et3N center dot 3HBr in the presence of potassium carbonate in dichloromethane at 0 degrees C to room temperature under very mild conditions in high yields and significantly shorter times, is reported.

Ruthenium(II)-CO complexes of N-[(2-pyridyl)methyliden]-alpha(or beta)-aminonaphthalene: Synthesis, spectral studies, crystal structure, redox properties and DFT calculation

The characterization and properties of trans-(X)-[RuX2(CO)(2)(alpha/beta-NaiPy)] (1, 2) (alpha-NaiPy (a), beta-NaiPy (b); X = Cl (1), I (2)) are described in this work. The structures are confirmed by single crystal X-ray diffraction studies. Reaction of these compounds with Me3NO in MeCN has isolated monocarbonyl trans-(X)-RuX2(CO)(MeCN)(alpha/beta-NaiPy)] (3, 4). The complexes show intense emission properties. Quantum yields of 1 and 2 (phi= 0.02-0.08) are higher than 3 and 4 (phi = 0.006-0.015). Voltammogram shows higher Ru(III)/Ru(II) (1.3-1.5 V) potential of 1 and 2 than that of 3 and 4 (0.8-0.9 V) that may be due to coordination of two pi-acidic CO groups in former. The electronic spectra and redox properties of the complexes are compared with the results obtained by density functional theory (DFT) and time-dependent density functional theory (TD-DFT) using polarizable continuum model (CPCM).

The fully implicit stochastic-alpha method for stiff stochastic differential equations

A fully implicit integration method for stochastic differential equations with significant multiplicative noise and stiffness in both the drift and diffusion coefficients has been constructed, analyzed and illustrated with numerical examples in this work. The method has strong order 1.0 consistency and has user-selectable parameters that allow the user to expand the stability region of the method to cover almost the entire drift-diffusion stability plane. The large stability region enables the method to take computationally efficient time steps. A system of chemical Langevin equations simulated with the method illustrates its computational efficiency.

Kinetic characterization of rat brain type IIA sodium channel alpha-subunit stably expressed in a somatic cell line

1. The rat brain type IIA Na+ channel alpha-subunit was stably expressed in Chinese hamster ovary (CHO) cells. Current through the expressed Na+ channels was studied using the whole-cell configuration of the patch clamp technique. The transient Na+ current was sensitive to TTX and showed a bell-shaped peak current vs. membrane potential relation. 2. Na+ current inactivation was better described by the sum of two exponentials in the potential range -30 to +40 mV, with. a dominating fast component and a small slower component. 3. The steady-state inactivation, h(infinity), was related to potential by a Boltzmann distribution, underlying thr ee states of the inactivation gate. 4. Recovery of the channels from inactivation at different potentials in the range -70 to -120 mV were characterized by al? initial delay which decreased with hyperpolarization. The time course was well fitted by the sum of two exponentials. In this case the slower exponential was the major component, and both time constants decreased with hyperpolarization. 5. For a working description of the Na+ channel inactivation in this preparation, with a minimal deviation from the Hodgkin-Huxley model, a three-state scheme of the form O reversible arrow I-1 reversible arrow I-2 was proposed, replacing the original two-state scheme of the Hodgkin-Huxley model, and the rate constants are reported. 6. The instantaneous current-voltage relationship showed marked deviation from linearity and was satisfactorily fitted by the constant-field equation. 7. The time course of activation was described by an m(x) model. However, the best-fitted value of x varied with the membrane potential and had a mean value of 2. 8. Effective gating charge was determined to be 4.7e from the slope of the activation plot, plotted on a logarithmic scale. 9. The rate constants of activation, alpha(m) and beta(m), were determined. Their functional dependence on the membrane potential was investigated.

Identification of alpha- and beta-Hydroxy Acid Containing Cyclodepsipeptides in Natural Peptide Mixtures Using Negative Ion Mass Spectrometry

Natural peptide libraries often contain cyclodepsipeptides containing alpha or beta hydroxy residues. Extracts of fungal hyphae of Isaria yield a microheterogenous cyclodepsipeptide mixture in which two classes of molecules can be identified by mass spectral fragmentation of negative ions. In the case of isaridins, which contain an alpha-hydroxy residue and a beta-amino acid residue, a characteristic product ion corresponding to a neutral loss of 72 Da is obtained. hi addition, neutral loss of water followed by a 72 Da loss is also observed. Two distinct modes of fragmentation rationalize the observed product ion distribution. The neutral loss of 72 Da has also been obtained for a roseotoxin component, which is also an alpha-hydroxy residue containing cyclodepsipeptide. In the case of isariins, which contain a beta-hydroxy acid residue, ring opening and subsequent loss of the terminal residue as an unsaturated ketene fragment, rationalizes the observed product ion formation. Fragmentation of negative ions provide characteristic neutral losses, which are diagnostic of the presence of alpha-hydroxy or beta-hydroxy residues.

The alpha method for solving differential algebraic inequality (DAI) systems

This paper describes an algorithm for ``direct numerical integration'' of the initial value Differential-Algebraic Inequalities (DAI) in a time stepping fashion using a sequential quadratic programming (SQP) method solver for detecting and satisfying active path constraints at each time step. The activation of a path constraint generally increases the condition number of the active discretized differential algebraic equation's (DAE) Jacobian and this difficulty is addressed by a regularization property of the alpha method. The algorithm is locally stable when index 1 and index 2 active path constraints and bounds are active. Subject to available regularization it is seen to be stable for active index 3 active path constraints in the numerical examples. For the high index active path constraints, the algorithm uses a user-selectable parameter to perturb the smaller singular values of the Jacobian with a view to reducing the condition number so that the simulation can proceed. The algorithm can be used as a relatively cheaper estimation tool for trajectory and control planning and in the context of model predictive control solutions. It can also be used to generate initial guess values of optimization variables used as input to inequality path constrained dynamic optimization problems. The method is illustrated with examples from space vehicle trajectory and robot path planning.

Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of C alpha Atoms

Background: The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the C alpha atoms which are modeled with modest accuracy. Methodology/Principal Findings: In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of C alpha atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of C alpha. We extend the method further to recognize potential protein-protein interface residues. Conclusion/Significance: Our approach to identify buried and exposed residues solely from the positions of C alpha atoms resulted in an accuracy of 84%, sensitivity of 83-89% and specificity of 67-94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70-96% and specificity of 58-94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from C alpha positions and all-atom models suggested that, recognition of interfacial residues using C alpha atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only C alpha positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.