Lagrangian molecular dynamics using selected conformational degrees of freedom,with application to the pseudorotation dynamics of furanose rings

Using internal conformational degrees of freedom for biopolymers as natural variables, and introducing a Lagrangian dynamics approach, one can simulate time-dependent processes over a much longer time scale than in classical Newtonian molecular dynamics (MD) techniques. Two factors contribute to this: a substantial reduction in the number of degrees of freedom and a very large increase in the size of the time step. We present the Lagrangian equations of motion for repuckering transitions in model furanose ( F ), ribose ( R ), and 2′-deoxyribose ( dR ) ring systems using the pseudorotation phase angle as the single dynamic variable. As in most Lagrangian analyses, the effective masses for the R and dR models are dependent on conformation, and we test the behavior of this variable mass (VM) model. Since the variation in effective mass is small, the VM model is compared with a simplified constant mass (CM) model, which is shown to be an excellent approximation. The equations of motion for the CM and VM models are integrated with the leapfrog and the iterative leapfrog algorithms, respectively. The Lagrangian dynamics approach reduces the number of degrees of freedom from about 40 to 1, and allows the use of time steps on the order of 20 fs, about an order of magnitude greater than is used in conventional MD simulations. © 1994 John Wiley & Sons, Inc.     

Simulations of the erythrocyte cytoskeleton at large deformation. II. Micropipette aspiration.

Coarse-grained molecular models of the erythrocyte membrane's spectrin cytoskeleton are presented in Monte Carlo simulations of whole cells in micropipette aspiration. The nonlinear chain elasticity and sterics revealed in more microscopic cytoskeleton models (developed in a companion paper; Boey et al., 1998. Biophys. J. 75:1573-1583) are faithfully represented here by two- and three-body effective potentials. The number of degrees of freedom of the system are thereby reduced to a range that is computationally tractable. Three effective models for the triangulated cytoskeleton are developed: two models in which the cytoskeleton is stress-free and does or does not have internal attractive interactions, and a third model in which the cytoskeleton is prestressed in situ. These are employed in direct, finite-temperature simulations of erythrocyte deformation in a micropipette. All three models show reasonable agreement with aspiration measurements made on flaccid human erythrocytes, but the prestressed model alone yields optimal agreement with fluorescence imaging experiments. Ensemble-averaging of nonaxisymmetrical, deformed structures exhibiting anisotropic strain are thus shown to provide an answer to the basic question of how a triangulated mesh such as that of the red cell cytoskeleton deforms in experiment.     

Freedom of choice and the tobacco endgame

Endgame proposals strive for a tobacco-free (or at least cigarette-free) society. Some endgame proposals are radical and include, for example, a complete ban on cigarettes. Setting aside empirical worries, one worry is ethical: would such proposals not go too far in interfering with individual freedom? I argue that concerns around freedom do not speak against endgame proposals, including strong proposals such as a ban on cigarettes. I first argue that when balancing freedom with public health goals in tobacco control, the latter win out. But I also argue that, in principle, a concern with freedom itself already justifies endgame measures. First, such measures can increase people's lifetime freedom, that is, the freedom they have across their entire lives. Second, such measures can facilitate a better interpersonal distribution of freedom by increasing aggregate societal freedom and by reducing inequalities. Overall, freedom does not preclude strict tobacco control but supports it.     

Polytomies and phylogenetically independent contrasts: examination of the bounded degrees of freedom approach

We examined the effect of soft polytomies on the performance (Type I error rate and bias) of Felsenstein's (1985; Am. Nat. 125:1-15) method of phylogenetically independent contrasts for estimating a bivariate correlation. We specifically tested the adequacy of bounding degrees of freedom, as suggested by Purvis and Garland (1993; Syst. Biol. 42:569-575). We simulated bivariate character evolution under Brownian motion (assumed by independent contrasts) and eight other models on five phylogenetic trees. For non-Brownian motion simulations, the adequacy of branch-length standardization was checked with a simple diagnostic (Garland et al., 1992; Syst. Biol. 41:18-32), and transformations were applied as indicated. Surprisingly, soft polytomies tended to have negligible effects on Type I error rates when models other than Brownian motion were used. Overall, and irrespective of evolutionary model, degrees of freedom were appropriately bounded for hypothesis testing, and unbiased estimates of the correlation coefficient were obtained. Our results, along with those of previous simulation studies, suggest that independent contrasts can reliably be applied to real data, even with phylogenetic uncertainty.     

Designs for Discrimination Between Bivariate Binary Response Models

We propose a test statistic for discrimination between alternative bivariate binary response models and the optimal design procedure which is an extension of T-optimality. Under certain conditions we prove that the maximum value of the power can be-obtained when the degrees of freedom of the test statistic is one. The conclusion is the same as that in discrimination between alternative univariate separate models. However the test statistics are different.     

Conditional transitions in gaze dynamics: role of vestibular nuclei in eye-only and eye/head gaze behaviors

 The gaze control system governs distinct gaze behaviors, including visual fixation and gaze reorientations. Transitions between these gaze behaviors are frequent and smooth in healthy individuals. This study models these gaze-behavior transitions for different numbers of gaze degrees of freedom. Eye/head gaze behaviors have twice the number of degrees of freedom as eye-only gaze behaviors. Each gaze behavior is observable in the system dynamics and is correlated with neuronal behaviors in several, coordinated neural centers, including the vestibular nuclei. The coordination among the neural centers establishes a sensorimotor state which maintains each gaze behavior. This study develops a mathematical framework for synthesizing the coordination among neural centers in gaze sensorimotor states and focuses on the role of vestibular nuclei neurons in gaze sensorimotor state transitions. Received: 17 December 1999 / Accepted in revised form: 3 May 2001     

An empirical Bayes approach to inferring large-scale gene association networks

MOTIVATION: Genetic networks are often described statistically using graphical models (e.g. Bayesian networks). However, inferring the network structure offers a serious challenge in microarray analysis where the sample size is small compared to the number of considered genes. This renders many standard algorithms for graphical models inapplicable, and inferring genetic networks an 'ill-posed' inverse problem. METHODS: We introduce a novel framework for small-sample inference of graphical models from gene expression data. Specifically, we focus on the so-called graphical Gaussian models (GGMs) that are now frequently used to describe gene association networks and to detect conditionally dependent genes. Our new approach is based on (1) improved (regularized) small-sample point estimates of partial correlation, (2) an exact test of edge inclusion with adaptive estimation of the degree of freedom and (3) a heuristic network search based on false discovery rate multiple testing. Steps (2) and (3) correspond to an empirical Bayes estimate of the network topology. RESULTS: Using computer simulations, we investigate the sensitivity (power) and specificity (true negative rate) of the proposed framework to estimate GGMs from microarray data. This shows that it is possible to recover the true network topology with high accuracy even for small-sample datasets. Subsequently, we analyze gene expression data from a breast cancer tumor study and illustrate our approach by inferring a corresponding large-scale gene association network for 3883 genes.     

Completion and refinement of 3-D homology models with restricted molecular dynamics: application to targets 47, 58, and 111 in the CASP modeling competition and posterior analysis

A method is presented to refine models built by homology by the use of restricted molecular dynamics (MD) techniques. The basic idea behind this method is the use of structure validation software to determine for each residue the likelihood that it is modeled correctly. This information is used to determine constraints and restraints in an MD simulation including explicit solvent molecules, which is used for model refinement. The procedure is based on the idea that residues that the validation software identifies as correctly positioned should be strongly constrained or restrained in the MD simulations, whereas residues that are likely to be positioned wrongly should move freely. Two different protocols are compared: one (applied to CASP3 target T58) using full structural constraints with separate optimization of each short fragment and the other (applied to T47) allowing some freedom using harmonic restraining potentials, with automatic optimization of the whole molecule. Structures along the MD trajectory that scored best in structural checks were selected for the construction of models that appeared to be successful in the CASP3 competition. Model refinement with MD in general leads to a model that is less like the experimental structure (Levitt et al. Nature Struct Biol 1999;6:108-111). Actually, refined T47 was slightly improved compared to the starting model; changes in model T58 led not to further enhancement. After the X-ray structure of the modeled proteins became known, the procedure was evaluated for two targets (T47 and the CASP4 target T111) by comparing a long simulation in water with the experimental target structures. It was found that structural improvements could be obtained on a nanosecond time scale by allowing appropriate freedom in the simulation. Structural checks applied to fast fluctuations do not appear to be informative for the correctness of the structure. However, both a simple hydrogen bond count and a simple compactness measure, if averaged over times of typically 300 ps, correlate well with structural correctness and we suggest that criteria based on these properties may be used in computational folding strategies.     

Probing protein fold space with a simplified model

We probe the stability and near-native energy landscape of protein fold space using powerful conformational sampling methods together with simple reduced models and statistical potentials. Fold space is represented by a set of 280 protein domains spanning all topological classes and having a wide range of lengths (33-300 residues) amino acid composition and number of secondary structural elements. The degrees of freedom are taken as the loop torsion angles. This choice preserves the native secondary structure but allows the tertiary structure to change. The proteins are represented by three-point per residue, three-dimensional models with statistical potentials derived from a knowledge-based study of known protein structures. When this space is sampled by a combination of parallel tempering and equi-energy Monte Carlo, we find that the three-point model captures the known stability of protein native structures with stable energy basins that are near-native (all α: 4.77 Å, all β: 2.93 Å, α/β: 3.09 Å, α+β: 4.89 Å on average and within 6 Å for 71.41%, 92.85%, 94.29% and 64.28% for all-α, all-β, α/β and α+β, classes, respectively). Denatured structures also occur and these have interesting structural properties that shed light on the different landscape characteristics of α and β folds. We find that α/β proteins with alternating α and β segments (such as the β-barrel) are more stable than proteins in other fold classes.     

Haemodynamic assessment of human coronary arteries is affected by degree of freedom of artery movement

Abnormal haemodynamic parameters are associated with atheroma plaque progression and instability in coronary arteries. Flow recirculation, shear stress and pressure gradient are understood to be important pathogenic mediators in coronary disease. The effect of freedom of coronary artery movement on these parameters is still unknown. Fluid–structure interaction (FSI) simulations were carried out in 25 coronary artery models derived from authentic human coronaries in order to investigate the effect of degree of freedom of movement of the coronary arteries on flow recirculation, wall shear stress (WSS) and wall pressure gradient (WPG). Each FSI model had distinctive supports placed upon it. The quantitative and qualitative differences in flow recirculation, maximum wall shear stress (MWSS), areas of low wall shear stress (ALWSS) and maximum wall pressure gradient (MWPG) for each model were determined. The results showed that greater freedom of movement was associated with lower MWSS, smaller ALWSS, smaller flow recirculation zones and lower MWPG. With increasing percentage diameter stenosis (%DS), the effect of degree of freedom on flow recirculation and WSS diminished. Freedom of movement is an important variable to be considered for computational modelling of human coronary arteries, especially in the setting of mild to moderate stenosis.

Abbreviations: 3D: Three-dimensional; 3DR: Three-dimensional Reconstruction; 3D-QCA: Three-dimensional quantitative coronary angiography; ALWSS: Areas of low wall shear stress; CAD: Coronary artery disease; CFD: Computational fluid dynamics; %DS: Diameter stenosis percentage; EPCS: End point of counter-rotating streamlines; FSI: Fluid–structure interaction; IVUS: Intravascular ultrasound; LAD: Left anterior descending; MWSS: Maximum wall shear stress; SST: Shear stress transport; TAWSS: Time-averaged wall shear stress; WSS: wall shear stress; WPG: Wall pressure gradient; MWPG: Maximum wall pressure gradient; FFR: Fractional flow reserve; iFR: Instantaneous wave-free ratio     

Bayes factor between Student t and Gaussian mixed models within an animal breeding context

The implementation of Student t mixed models in animal breeding has been suggested as a useful statistical tool to effectively mute the impact of preferential treatment or other sources of outliers in field data. Nevertheless, these additional sources of variation are undeclared and we do not know whether a Student t mixed model is required or if a standard, and less parameterized, Gaussian mixed model would be sufficient to serve the intended purpose. Within this context, our aim was to develop the Bayes factor between two nested models that only differed in a bounded variable in order to easily compare a Student t and a Gaussian mixed model. It is important to highlight that the Student t density converges to a Gaussian process when degrees of freedom tend to infinity. The twomodels can then be viewed as nested models that differ in terms of degrees of freedom. The Bayes factor can be easily calculated from the output of a Markov chain Monte Carlo sampling of the complex model (Student t mixed model). The performance of this Bayes factor was tested under simulation and on a real dataset, using the deviation information criterion (DIC) as the standard reference criterion. The two statistical tools showed similar trends along the parameter space, although the Bayes factor appeared to be the more conservative. There was considerable evidence favoring the Student t mixed model for data sets simulated under Student t processes with limited degrees of freedom, and moderate advantages associated with using the Gaussian mixed model when working with datasets simulated with 50 or more degrees of freedom. For the analysis of real data (weight of Pietrain pigs at six months), both the Bayes factor and DIC slightly favored the Student t mixed model, with there being a reduced incidence of outlier individuals in this population.     

Sample reproducibility of genetic association using different multimarker TDTs in genome-wide association studies: characterization and a new approach

Multimarker Transmission/Disequilibrium Tests (TDTs) are very robust association tests to population admixture and structure which may be used to identify susceptibility loci in genome-wide association studies. Multimarker TDTs using several markers may increase power by capturing high-degree associations. However, there is also a risk of spurious associations and power reduction due to the increase in degrees of freedom. In this study we show that associations found by tests built on simple null hypotheses are highly reproducible in a second independent data set regardless the number of markers. As a test exhibiting this feature to its maximum, we introduce the multimarker 2-Groups TDT (mTDT(2G)), a test which under the hypothesis of no linkage, asymptotically follows a χ2 distribution with 1 degree of freedom regardless the number of markers. The statistic requires the division of parental haplotypes into two groups: disease susceptibility and disease protective haplotype groups. We assessed the test behavior by performing an extensive simulation study as well as a real-data study using several data sets of two complex diseases. We show that mTDT(2G) test is highly efficient and it achieves the highest power among all the tests used, even when the null hypothesis is tested in a second independent data set. Therefore, mTDT(2G) turns out to be a very promising multimarker TDT to perform genome-wide searches for disease susceptibility loci that may be used as a preprocessing step in the construction of more accurate genetic models to predict individual susceptibility to complex diseases.     

Phenomenological and molecular models of biological proton transfers

In this paper, the main models describing the transfer of a proton in a molecular system are presented. Models valid when the intersite coupling is weak (non-adiabatic and electronically adiabatic regimes) and strong (adiabatic regime) are described. We distinguish molecular models in which the rate constant is obtained by considering explicitly various degrees of freedom of the system and simpler, phenomenological models built to account for the kinetic isotope effect. The relations between the various models are discussed. Their application to specific systems is illustrated by several studies reported in the literature, with a special emphasis on biological systems.     

Two-arm trajectory planning in a manipulation task

This paper addresses the trajectory planning problem for a task which requires positioning and orienting an object firmly grasped by two hands at a visually specified goal configuration in the horizontal plane. The motor task involves three degrees of freedom (two translational and one rotational), and the motions of the arms are constrained by the physical coupling through the held object. Experimentally measured trajectories of two arms in the coordinated positioning/orienting task are presented. The hypothesis that the rotational and translational components of motions are decoupled and independently planned is tested. Two explicit mathematical models to account for the kinematic features of the two-arm motions are formulated, and the predictions of the models are compared with the experimental data. Both models extend the minimum-jerk model to the two-arm coordinated motions case. The trajectories predicted by the models were found to be in qualitative agreement with the experimental data. However, neither model could account for the observed configuration dependence of the motions, nor for some of the properties of the measured velocity components of the motions. Our findings support the idea that the rotational and translational components of two-arm motions in the positioning/orienting task are independently planned in extra-personal space, and are further combined in a hierarchical fashion to produce the observed motions. The tested models may serve as a basis for further investigations of issues pertinent to the generation of two-arm trajectories. Received: 27 March 1995 / Accepted in revised form: 17 September 1996     

An explanation of the relationships between mass, metabolic rate and characteristic skeletal length for birds and mammals

A theory explaining the relationships between metabolic rate and body mass for birds and mammals is developed in terms of the mitochondrion theory of aerobic metabolism, the sliding filament theory of muscle contraction, simple models of vertebrate anatomy and activity, and other propositions. Both body mass and metabolic rates are shown to result from a homogeneous vertebrate design quantitatively expressed in terms of a set of nearly invariant parameters with five degrees of freedom: propulsion technique, mitochondrion capability, non-skeletal muscle mass exponent, characteristic skeletal length, and sturdiness factor. The first three of the degrees of freedom are phylogenetic group specific. The last two vary considerably even within a single species. The theory is shown to agree satisfactorily with placental mammal, marsupial mammal, passerine bird and nonpasserine bird data. Algorithms for determining metabolic rates and body mass as functions of skeletal characteristic length are developed and quantitative estimates of the constants occurring in the algorithms are given.     

Quantitative evaluation of complex surveillance systems for pest and disease detection

Surveillance for disease detection is used primarily for early detection of incursions and to support assertions of freedom from disease.Analytical techniques used to evaluate the efficacy of such surveillance are equally applicable across the domains of invasive species (both plant and animal), diseases and pests of agriculture crops, livestock, and of wildlife. Scenario tree models of surveillance activities may be used to estimate their diagnostic sensitivities, or the probability that the target organism will be detected given that it is present at a defined level. This paper will outline techniques for estimating the sensitivity of both targeted and general surveillance activities, and for the surveillance system as a whole. Probability of freedom from the target organism may be estimated from the surveillance sensitivity, and this Bayesian approach may be extended to estimate current probability of freedom from appropriate use of historical and ongoing surveillance evidence.     

A comparison of three maximum likelihood models for stage-frequency data

A comparison has been made between the estimates obtained from maximum likelihood estimation of gamma, inverse normal, and normal distribution models for stage-frequency data. Results have been compared for six of sets of test data, and from many sets of simulated data. It is concluded that (1) some estimates may differ substantially between the models, (2) estimates from the correct model have little bias, and estimated standard errors are generally close to theoretical values, (3) there are problems in determining degrees of freedom for chi-squared goodness of fit tests, so that it is best to compare test statistics with simulated distributions, and (4) goodness of fit tests may not discriminate well between the three models.     

A stereochemical model of the transpeptidation complex

Molecular models are proposed to describe the relative arrangement of aminoacyl and peptidyl tRNAs when bound to their respective A and P sites on the ribosome. The crystallographically determined structures of tRNAasp and tRNAphe have served as the models for these bound structures, while the imposed steric constraints for the model complexes were based on the results of published experimental data. The constructed models satisfy the stereochemical requirements needed for codon-anticodon interaction and for peptide bond formation. In this paper, the results of the complex containing tRNAphe as the A and P site bound transfer RNAs, is compared to a similarly constructed model which uses tRNAasp as the ribosome-bound transfer RNAs. The models have the following three major features: 1) the aminoacyl and peptidyl transfer RNAs assume an angle of approximately 45 degrees relative to each other; 2) in providing the proper stereochemistry for peptide bond condensation, a significant kink must be present in the messenger RNA between the A site and P site codons; and 3) a comparison of the two model complexes indicates that structural variations between the tRNAs or any allosteric transitions of the transfer RNAs associated with codon-anticodon recognition may be accommodated in the model by way of freedom of rotation about the phosphate backbone bonds in the mRNA between consecutive codons.     

The maximum entropy formalism and the idiosyncratic theory of biodiversity

Why does the neutral theory, which is based on unrealistic assumptions, predict diversity patterns so accurately? Answering questions like this requires a radical change in the way we tackle them. The large number of degrees of freedom of ecosystems pose a fundamental obstacle to mechanistic modelling. However, there are tools of statistical physics, such as the maximum entropy formalism (MaxEnt), that allow transcending particular models to simultaneously work with immense families of models with different rules and parameters, sharing only well‐established features. We applied MaxEnt allowing species to be ecologically idiosyncratic, instead of constraining them to be equivalent as the neutral theory does. The answer we found is that neutral models are just a subset of the majority of plausible models that lead to the same patterns. Small variations in these patterns naturally lead to the main classical species abundance distributions, which are thus unified in a single framework.