PubMed HealthSearch

SEARCH · PubMed Health

Results for “Structural equation modeling”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5Linked to original sources

Formation of topographic maps and columnar microstructures in nerve fields.

Topographic connections are found in many parts of the vertebrate nervous systems, known for example as retinotopy. The self-organizing ability of Hebb type modifiable synapses plays an important role in forming, at least in refining, the topographic connections. We present a mathematical analysis of a revised version of the Willshaw-Malsburg model of topographic formation, solving the equations of synaptic self-organization coupled with the field equation of neural excitations. The equilibrium solutions are obtained and their stability is studied. It is proved that two cases exist depending on parameters. In one case, the smooth topographic organization is obtained as a stable equilibrium of the equations. In the other case, this solution becomes unstable, and instead the topographic organization with columnar microstructures appears. This might explain the columnar structures in the cerebrum. The theory is confirmed by computer simulated experiments.

Computers

Representation of Nonepistatic selection models and analysis of multilocus Hardy-Weinberg Equilibrium configurations.

The paper develops conditions for the existence and the stability of central equilibria emanating from selection recombination interaction with generalized nonepistatic selection forms operating in multilocus multiallele systems. The selection structure admits a natural representation as simple sums of Kronecker products based on a common set of marginal selection components. A flexible parametrization of the recombination process is introduced leading to a canonical derivation of the transformation equations connecting gamete frequency states over successive generations. Conditions for the existence and stability of multilocus Hardy-Weinberg (H.W.) type equilibria are elaborated for the classical nonepistatic models (multiplicative and additive viability effects across loci) as well as for generalized nonepistatic selection expressions. It is established that the range of recombination distributions maintaining a stable H.W. polymorphic equilibrium is confined to loose linkage in the pure multiplicative case, but is not restricted in the additive model. In the bisexual case we ascertain for the generalized nonepistatic model the stability conditions of a common H.W polymorphism.

Alleles

The [18F]fluorodeoxyglucose method for the measurement of local cerebral glucose utilization in man.

A method has been developed to measure local glucose consumption in the various structures of the brain in man with three-dimensional resolution. [18F]-2-deoxy-2-fluoro-D-glucose is used as a tracer for the exchange of glucose between plasma and brain and its phosphorylation by hexokinase in the tissue. A mathematical model and derived operational equation are used which enable local cerebral glucose consumption to be calculated in terms of the following measurable variables. An intravenous bolus of [18F]-2-deoxy-2-fluoro-D-glucose is given and the arterial specific activity monitored for a predetermined period of from 30 to 120 minutes. Starting at 30 minutes, the activity in a series of sections through the brain is determined with three-dimensional resolution by an emission tomographic scanner. The method was used to measure local cerebral glucose consumption in two normal volunteers. The values in gray matter structures range from 5.79 mg/100 g per minute in the cerebellar cortex to 10.27 in the visual cortex, whereas, in white matter structures, the values range from 3.64 mg/100 g per minute in the corpus callosum to 4.22 in the occipital lobe. Average values for gray matter, white matter, and whole brain metabolic rates, calculated as a weighted average based on the approximate volume of each structure, are 8.05, 3.80, and 5.90 mg/100 g per minute, respectively. The value of 5.9 mg/100 g per minute compares favorably with values previously reported.

Adult

[Calculations of airway parameters from expiratory concentration volume values based on the model of stationary bronchoalveolar gas exchange (author's transl)].

Approximating the cumulative cross section of the terminal airway structure by a single parabolic exponential term with the unknown parameters initial cross section a, exponential steepness b and external dead space VDE, calculations of the stationary concentration front F along airway length L, raising during constant flow inspiration, were done, solving the gas transport equation for stationary conditions numerically and integrating the distribution of dF/dL. The conversion to the expiratory F-V profile was done neglecting pulmonary inhomogeneity and expiratory front deformation. Calculations were applied on a scale of different a and b enclosing normal and some pathological lung models. Applying different inspiratory flows, there is a linear correlation between the reciprocal steepness dV/dF of expiratory F-V diagrams and the flow diffusion ratio V/D*, where D* is the effective diffusivity of gases. Compared to normal, the model of nonobstructive emphysema shows higher values of dV/dF and the general obstruction lower ones; all three have the same increase of dV/dF along V/D*. The peripheral obstruction type has the highest slope of regression along V/D* and the obstructive emphysema the lowest one. Mathematical approximations of the regression between dV/dF of the expiratory F-V diagram at certain intervals of F and V/D* were used to present a procedure calculating airway parameters a, b, VDE and the effective diffusivity D* from a sequence of F-V diagrams applying different inspiratory flows.

Airway Obstruction

Light-scattering investigations of the subunit dissociation of human hemoglobin A. Effects of various neutral salts.

The effectiveness of various salts of the Hofmeister series as dissociating agents for human hemoglobin A tetramers has been investigated by light-scattering molecular-weight measurements. Dissociation of hemoglobin to half-molecules of alpha beta dimers follows the order of the series dictated predominantly by the sequence of the anions F- less than Cl- less than Br- less than ClO4- less than SCN-, I-, with the cations Na+ and K+ having relatively little effect on the observed dissociation. The use of equations derived for predicting the effects of dissociating reagents on the structure of subunit proteins [Herskovits, T. T., and Ibanez, V. S. (1976), Biochemistry 15, 5715] together with Setschenow constants based on the model amino acid data of Nandi and Robinson were found to give a satisfactory account of the dissociation behavior observed with many of the salts, giving reasonable estimates of the number of amino acids that form the smaller contact area of the alpha beta subunits of hemoglobin shown by the Perutz crystallographic model. The analysis of the dissociation data also extends the utility of the Setschenow constants tested for the characterization of the dissociation behavior of other subunit proteins.

Anions

Neuronal excitability modulation over the sleep cycle: a structural and mathematical model.

A model for control of the desynchronized phase of the sleep cycle postulates reciprocal interaction between cells in the pontine gigantocellular tegmental field (FTG cells) and cells in the nucleus locus coeruleus and nucleus subcoeruleus (LC cells). This physiological model leads to equations of the Lotka-Volterra type; the time course of activity predicted by the model is in good agreement with actual long-term recordings of FTG cells and single-cycle data for LC cells.

Animals

Simple mathematical models with very complicated dynamics.

First-order difference equations arise in many contexts in the biological, economic and social sciences. Such equations, even though simple and deterministic, can exhibit a surprising array of dynamical behaviour, from stable points, to a bifurcating hiearchy of stable cycles, to apparently random fluctuations. There are consequently many fascinating problems, some concerned with delicate mathematical aspects of the fine structure of the trajectories, and some concerned with the practical implications and applications. This is an interpretive review of them.

Diffusion

A model of pattern formation in insect embryogenesis.

A model is proposed in which the interaction of an autocatalytic substance with a short diffusion range--the activator--and its more diffusible antagonist--the inhibitor--leads to a local high concentration of activator at the posterior pole of the egg. The inhibitor, which is then produced mainly in this activated region, diffuses into the rest of the egg, where it acts as a 'morphogen', forming a concentration gradient which supplies positional information. This model can account quantitatively for the patterns resulting from a large number of different experiments performed during early insect development, including ligation, u.v.-irradiation and microsurgical manipulations. The formation of additional posterior structures is interpreted as the result of the appearance of a new activator peak. Omission of segments after ligation of the egg is explained as the result of accumulation of morphogen (the inhibitor) at the posterior side of the ligation and a decrease of morphogen on the anterior side. In order to account for certain quantitative features of the ligation experiments it is necessary to assume that determination in response to the morphogen gradient is a slow, stepwise process, in which the nuclei or cells first pass through determination stages characteristic for more anterior structures until they ultimately form a given structure. The mutual interactions of activator and inhibitor are expressed as a set of partial differential equations. The individual experiments are simulated by solving these equations by use of a computer.

Animals

Human adipocyte cholesterol. Concentration, localization, synthesis, and turnover.

By analysis of 124 specimens in 16 different patients, isolated human adipocyte cholesterol concentration is highly correlated with fat cell size but not with plasma cholesterol concentration. Less than 6 percent of total cholesterol is esterified; after subcellular fractionation, 88 percent of the cholesterol is recovered in the triglyceride-rich supernatant oil. This latter finding supports the observation that fat cell cholesterol is determined by triglyceride content, and hence by fat cell size. After intravenous administrtion of radioactive cholesterol, the sum of a three-exponential equation was fit simultaneously to both the plasma and adipocyte specific activity time curves in six patients. In five of the six, a slowly turning over pool (pool 3) closely fit the adipocyte data. Two model structures, mammillary and catenary, were fitted to the data. There was no synthesis in pool 3 using a mammillary model but a mean 5.3 percent of the total body production rate was found in compartment 3 if a catenary model was assumed. Although a catenary model is biologically unlikely, it could not be excluded. Obesity is associated with an increased cholesterol synthetic rate equal to 20 mg/day for each kilogram of body fat. To test (by an independent method) if this synthesis might be occurring in adipose tissue, human fat cells were obtained under a wide variety of dietary conditions and incubated in vitro with radioactive glucose or acetate. Incorportation of these precursors into sterol could account for no more than 1 mg cholesterol synthesis/kg fat per day. These in vitro data taken together with the in vivo mammillary compartmental analysis data are compatible with the possiblity that the excess cholesterol synthesis of obesity occurs in pool 1, most likely from hepatic or intestinal sites.

Adipose Tissue

A continuous selective model for an X-linked locus.

Neglecting age-structure, but taking into account matings with differential fertility in Mendelian reproduction, a continuous selective model is formulated for a single X-linked locus with an arbitrary number of alleles. Without restricting the mating system, differential equations are derived for the genotypic and allelic frequencies. Assuming random mating, no selection, and constant fertilities and mortalities, these differential equations are solved explicitly. For this case, in contrast to the corresponding phenomenon in the usual model with discrete, non-overlapping generation, the difference between the frequencies of any allele in males and females approaches zero without oscillation.

Age Factors

Linear relationships between plasma binding and lipophilicity of disopyramide derivatives.

The extent of plasma binding and the partition coefficient of disopyramide and 20 disopyramide derivatives were determined. Structural variations on the four functional groups around the tetrahedral carbon in the disopyramide molecule was found to influence both parameters to varying degrees. Three linear equations were developed to correlate the observed effects, depending on the type of chemical modification. The linear correlation between drug-plasma interaction and lipophilic character was analyzed theoretically. A simple model was derived to relate quantitatively the variation in the extent of plasma binding to the change in lipophilicity of disopyramide derivatives.

Arrhythmias, Cardiac

Magnetorheography; calculation of blood flow from surface induced potentials.

This research is an attempt to measure blood flow without doing any hurt to the body. When the body is placed in a high magnetic field the blood-flow-induced electromotive force (EMF) can be detected with the surface electrodes using the same principle as an electromagnetic flow-meter. This blood flow measurement is named magnetorheography (its recording: MRG). A theoretical analysis of the potential in the tissue shows that blood flow can be calculated from the surface-induced EMF (MRG) without measuring the radius and the depth of the vessel from the skin surface, when the surface is fairly flat in the vicinity of the vessel. A model experiment verified the theory. In order to apply the theory to the measurement in an in situ situation, a structural coefficient was introduced which is related to the external shape, internal tissue compositions and their impedance of the body. Using such a coefficient the flow calculated from MRG, by an equation including the coefficient, agreed with the flow actually recorded. In experiments on the thighs of 8 dogs weighing 7 to 17 Kg, MRG proved to be proportional to arterial flow and the coefficients were approximately constant. These results indicate that blood flow may possibly be estimated by the calculation from MRG.

Animals

Nonlinear regression methods in design of experiments and mathematical modelling. Applications to the analysis of the steady-state kinetics of glutathione reductase.

A branching reaction pathway involving a ping pong and a sequential loop has been proposed for glutathione reductase (Biochem. Biophys. Res. Commun. 53 (1973) 1151). In the present investigation nonlinear regression methods have been applied in the fitting of rate equations to experimental data to test the validity of the model proposed and to discriminate between alternative mathematical models (cf. FEBS Lett. 26 (1972) 252). In the best rate law, some of the parameters were numerically redundant. Therefore, a feature-wise analysis of the rate equation was carried out by varying one substrate concentration at a time. The overall strategy used was a cyclic procedure involving: experimentation - analysis of data - modelling - design of experiments - new experimentation etc. Consideration was given to the experimental error structure and to the importance of weighting in the regression analysis. In the design of experiments for discrimination between rival models, a previously defined discrimination function was used. The results of the analysis support the branching reaction scheme proposed for glutathione reductase.

Glutathione

Analysis of afferent responses from isolated semicircular canal of the guitarfish using rotational acceleration white-noise inputs. II. Estimation of linear system parameters and gain and phase spectra.

Quantitative estimates were computed for exponential coefficients and rate constants contributing to afferent unit impulse responses obtained from bundles innervating specific regions of the semicircular canal. The grouping of these estimates into specific response classes provided quantitative correlations with specific anatomical regions of innervation of the crista. Linear system gain and phase spectra were computed also, by applying Fourier transformations to unit impulse responses, for purposes of comparison with previous studies employing frequency domain analyses. Responses fitted by third-order linear system equations were specific to afferents innervating the crest and transition regions of the crista; whereas those fitted by overdamped, second-order equations were specific to afferents innervating the slopes and transition crista regions. It was concluded that strictly mechanical models of the transduction process are inadequate to account for the diverse and spatially distributed classes of observed responses and, moreover, structural features such as different hair cell types or efferent innervation effects could be excluded as inoperative in this preparation. The alternative hypothesis was suggested that certain of the observed subcomponents could be direct reflections of the initial mechanical stimulus, but that other subcomponents were reflections of more complex filtering mechanisms operating at the cellular or synaptic levels.

Animals

Acquisition of the covalent quaternary structure of an immunoglobulin G molecule. Theoretical reoxidation models.

A theoretical format suitable for analyzing diverse complex kinetic systems where reaction pathways may exhibit cooperativity, is developed to account for the in vitro kinetics of air reoxidation of a human IgGlkappa immunoglobulin, in which the four interchain disulfide bonds have been reduced (Sears, D.W., et al. (1977), Biochemistry 16 (preceding paper in this issue)). The equations relate experimentally determined concentrations of the reactants, product, and macroscopic intermediates to the probability of occurrence of any of 12 distinct microscopic states. The concentrations of the macroscopic species--light chains (L), heavy chains (H), covalently associated intermediates HL, H2, H2L, and product H2L2--are also explicitly related to an observable function of those concentrations, the sulfhydryl titer, r. Since values of r can be calculated for any microscopic state of the system, the theory is in principle capable of a complete description of the experimental reaction course. In practice, limited experimental information precludes a unique solution of the equations at present, and it was not judged worthwhile to attempt curve fitting or approximation of probability terms with adjustable parameters of unknown physical significance. Certain special cases of the theory are, however, readily and exactly solved. These include models in which reoxidation is a random process and others in which the probabilities for inter-HL and inter-HH bond formation are different but independent of one another. The experimental results in the case of the IgGlkappa studied here clearly depart form the predicted behavior in either of these models. The initial probability for formation of a bond between a heavy and light chain is 1.5 to 2 times greater than for a bond between heavy chains. From the fact that this ratio changes as the reaction proceeds, and from the pattern of variation in concentration of intermediates during the reaction, it is concluded that the reoxidation process is not random, and that the bonds do not form independently, exhibiting instead kinetic cooperativity. The results are discussed in terms of assembly pathways in this and related systems. A novel feature of the theory is that it eliminates time as an explicit variable in the treatment of the kinetic process. This makes it especially useful for discerning whether the formation of one bond influences the reaction probability of another in a system where several similar or intrinsically identical reactions occur, and where kinetic order is difficult to establish. While the theory here is formulated in terms of the reoxidation reactions, it is, for example, equally applicable to the reduction process described (Sears, D.W., et al. (1977b), Biochemistry 16 (following paper in this issue)).

Disulfides

Dynamic properties of polyelectrolyte calcium membranes.

Shashoua observed spontaneous oscillations in a polyelectrolyte membrane formed by interfacial precipitates of polyacid and polybase. We have here undertaken experimental and theoretical studies of polyglutamic acid-Ca++ membrane in order to clarify the processes involved in this dynamic behavior. We find a region of distinct hysteresis in the voltage current curve for this system. A sharp transition from a state of low membrane resistance to one of high resistance occurs at a current density different from that of inverse transition. This membrane system is modeled as a two layer structure: a negatively charged layer alpha made of ionized polyelectrolyte in series with a neutral region beta in which the polymeric ionic sites are masked by calcium ion. This structure results in a difference in the transference number for the mobile ions, causing salt accumulation at the interfacial region during a current flow in the alpha to beta direction. This altered salt concentration induces a change of polymeric conformation, which in turn affects the membrane permeability and the rate of accumulation. Based upon nonequilibrium thermodynamic flow equations, and a two-state representation of membrane macromolecular conformation, this model displays a region of hysteresis in the current range of experimental observations.

Calcium

Evolution of enzyme catalytic power. Characteristics of optimal catalysis evaluated for the simplest plausible kinetic model.

1. Evolutionary changes in the structure of an enzyme that provide an increase in its K(m) value are considered. Provided that K(m) increases as a result of increases in the forward rate constants of the catalysis relative to the reverse rate constants, the enzyme catalyses the conversion of a fixed concentration of its substrate more rapidly when its structure provides that K(m)>[S] than when K(m)<[S]. 2. Catalytic efficiency of enzymes is discussed in terms of the simplest plausible model, the Haldane [(1930) Enzymes, Longmans, London] reversible three-step model: [Formula: see text] The rate equation for the forward reaction of this model (formation of P) may be written in the simple form: [Formula: see text] K(eq.) is the equilibrium constant (=[P](eq.)/[S](eq.)), and k(cat.)=V/[E](T), where [E](T) is the total enzyme concentration. 3. To assess the effectiveness of an enzyme, it is necessary only to determine the extent to which the constraints of a particular kinetic mechanism permit v(2) (v when K(m)>>[S]) to approach v(d) (the diffusion-limited rate). 4. The value of the optimal rate of catalysis (v(opt.), the maximal value of v(2)) is dictated by the equilibrium constant for the reaction, K(eq.); v(2)=v(d)/a, where [Formula: see text] when k(+1) is assumed equal to k(-3), and v(opt.)=v(d)/a(min.). When K(eq.)>/=1, it is necessary that k(+2)>>k(-1) for a to take its minimum value, a(min.); when K(eq.)<<1, it is necessary only that k(+2)>>K(eq.).k(-1), i.e. a can equal a(min.) even if k(+2) >1, v(opt.)=v(d); when K(eq.)=1, v(opt.)=v(d)/2, and when K(eq.)<<1, v(opt.)=K(eq.).v(d). 5. The analysis, together with predicted effects of evolutionary pressure, suggests that in practice the rates of the fastest enzyme-catalysed freely reversible reactions might be expected to be lower than the value of k(+1)[E](T)[S] by about an order of magnitude, particularly if K(eq.)<1. 6. The existing literature suggests that, in general, appropriate values of K(m) have evolved for the provision of high rates of catalysis but that many values of k(cat.) are not large enough to provide optimal rates of catalysis unless the value of k(+1)in vivo is lower than its value in free solution.

Biological Evolution

Multiple inhibition of glutathione S-transferase A from rat liver by glutathione derivatives: kinetic analysis supporting a steady-state random sequential mechanism.

Glutathione derivatives inhibit glutathione S-transferase A [cf. Biochem. J. (1975) 147, 513--522]. The steady-state kinetics of this inhibition have been investigated in detail by using S-octyglutathione, glutathione disulphide and S-(2-chloro-4-nitrophenyl)glutathione: the last compound is a product of the enzyme-catalused reaction. Interpreted in terms of generalized denotations of inhibition patterns, the compounds were found to be competitive with the substrate glutathione. Double-inhibition experiments involving simultaneous use of two inhibitors indicated exclusive binding of the inhibitors to the enzyme. The discrimination between alternative rate equations has been based on the results of weighted non-linear regression analysis. The experimental error was determined by replicate measurements and was found to increase with velocity. The established error structure was used as a basis for weighting in the regression and to construct confidence levels for the judgement of goodness-of-fit of rate equations fitted to experimental data. The results obtained support a steady-state random model for the mechanism of action of glutathione S-transferase A and exclude a number of simple kinetic models.

Animals