Experimental tests of villin subdomain folding simulations

We have used laser temperature-jump to investigate the kinetics and mechanism of folding the 35 residue subdomain of the villin headpiece. The relaxation kinetics are biphasic with a sub-microsecond phase corresponding to a helix-coil transition and a slower microsecond phase corresponding to overall unfolding/refolding. At 300 K, the folding time is 4.3(+/-0.6) micros, making it the fastest folding, naturally occurring protein, with a rate close to the theoretical speed limit. This time is in remarkable agreement with the prediction of 5 (+11,-3) micros by Zagrovic et al. from atomistic molecular dynamics simulations using an implicit solvent model. We test their prediction that replacement of the C-terminal phenylalanine residue with alanine will increase the folding rate by removing a transient non-native interaction. We find that the alanine substitution has no effect on the folding rate or on the equilibrium constant. Implications of this result for the validity of the simulated folding mechanism are discussed.     

Bayesian survival analysis using a MARS model

A Bayesian multivariate adaptive regression spline fitting approach is used to model univariate and multivariate survival data with censoring. The possible models contain the proportional hazards model as a subclass and automatically detect departures from this. A reversible jump Markov chain Monte Carlo algorithm is described to obtain the estimate of the hazard function as well as the survival curve.     

On the Bayesian estimation of a closed population size in the presence of heterogeneity and model uncertainty

Summary .   We consider the estimation of the size of a closed population, often of interest for wild animal populations, using a capture–recapture study. The estimate of the total population size can be very sensitive to the choice of model used to fit to the data. We consider a Bayesian approach, in which we consider all eight plausible models initially described by Otis et al. (1978, Wildlife Monographs 62, 1–135) within a single framework, including models containing an individual heterogeneity component. We show how we are able to obtain a model-averaged estimate of the total population, incorporating both parameter and model uncertainty. To illustrate the methodology we initially perform a simulation study and analyze two datasets where the population size is known, before considering a real example relating to a population of dolphins off northeast Scotland.     

Bayesian meta-analysis for longitudinal data models using multivariate mixture priors

We propose a class of longitudinal data models with random effects that generalizes currently used models in two important ways. First, the random-effects model is a flexible mixture of multivariate normals, accommodating population heterogeneity, outliers, and nonlinearity in the regression on subject-specific covariates. Second, the model includes a hierarchical extension to allow for meta-analysis over related studies. The random-effects distributions are decomposed into one part that is common across all related studies (common measure), and one part that is specific to each study and that captures the variability intrinsic between patients within the same study. Both the common measure and the study-specific measures are parameterized as mixture-of-normals models. We carry out inference using reversible jump posterior simulation to allow a random number of terms in the mixtures. The sampler takes advantage of the small number of entertained models. The motivating application is the analysis of two studies carried out by the Cancer and Leukemia Group B (CALGB). In both studies, we record for each patient white blood cell counts (WBC) over time to characterize the toxic effects of treatment. The WBCs are modeled through a nonlinear hierarchical model that gathers the information from both studies.     

Gender difference in the relationship of performance in the handgrip and standing long jump tests to lean limb volume in young adults

Groups of young, adult males and females performed the handgrip and standing long jump tests. Their total forearm and leg volumes were calculated from a series of circumference and length measurements, and the lean volumes (bone + muscle) calculated by taking the skinfold thickness into consideration. In the handgrip, the mean female performance was 298 N compared with 496 N for the males. In the standing long jump, mean performance expressed as distance x body mass was 87.3 kg.m for females compared with 137.7 kg.m for males. These superior performances of males could simply reflect their greater muscle mass, as the mean lean volumes of female and male limbs respectively were 0.54 l and 0.89 l for forearms, and 11.82 l and 14.82 l for the two legs. However, when the performances of males and females were grouped by lean limb volume, it was found that while in both tests there were linear relationships, males and females did not share a common line. In both tests the male relationship was at a higher level than the female; therefore, for a given lean volume, the male performance was significantly superior to that of the female. The gender difference found in this study has not been seen in other studies in which the performance of skeletal muscle has been related to the cross-sectional area of the active muscles and the possible reasons for the differences are considered.     

Proton transfer reactions in aqueous solutions of pyridoxamine phosphate

UV-visible and ¹³C NMR measurements described in the literature and our ³¹P NMR measurements support the following mechanism of proton transfer reactions in aqueous solutions of pyridoxamine phosphate: Only the tautomeric equilibrium between neutral form, A _N, and zwitterion, A _Z, which is analogous to the tautomeric equilibrium of 3-hydroxypyridine in aqueous solution, is important, and that equilibrium does not change upon the dissociation of the second phosphate proton. With these simplifying assumption, we have simulated the relaxation spectrum of the proton transfer reactions of pyridoxamine phosphate in water using parameters from analogous reactions and compared it with our ultrasound and temperature jump measurements. We have found that the relaxation process measured by the temperature jump experiment is mainly caused by the overall reaction A _N=A _Z (or A _N ^- =A _Z ^- ) and the ultrasound absorption at the isoelectric point between pK₂ and pK₃ is mainly caused by the overall reaction .     

Physicochemical aspects of carbohydrate binding to some plant lectins with binding preference forN-acetylgalactosamine and galactose

This contribution illustrates the advantages of some chromophoric and fluorophoric carbohydrate derivatives such asp-nitrophenyl (pNO₂Phe) or 4-methylumbelliferyl (MeUmb) glycosides andN-dansylgalactosamine in studies of the binding equilibrium and kinetics with some plant lectins. The methods used involve continuous titrations of changes in ligand or protein absorption and ligand fluorescence, including substitution titrations as well as stopped-flow, temperature-jump or pressure-jump relaxation kinetics. When monitored by temperature-jump relaxation, binding of MeUmbαGal to the bloodgroup A specific lectin GSAI-A₄ fromGriffonia simplicifolia is a simple bimolecular association with parametersk ₊ = 9.4 × 10⁴ M^-1 s^-1 andk _-1 = 5.3 s^-1 at 23°C, but binding to the GSAI-B₄ lectin is biphasic. The complementarity of the peanut agglutinin binding site with Galβ1 → 3GalNAc that occurs in manyO-glycoproteins follows from enthalpic considerations and also from the value of the dissociation-rate parameterk _-1 = 0.24 s^-1 of the MeUmbβGalβl → 3GalNAc.lectin complex. This value, obtained by stopped-flow kinetics is 100 times smaller than for other mono-and disaccharides investigated. The binding mechanism is simple and the derivatisation of Galβ1 → 3GalNAc does not affect the affinity to a considerable degree. The binding preference of tetravalentsoybean agglutinin for MeαGalNAc over MeαGal by a factor of 25 is mainly of enthalpic origin with an additional 7 kJ mol^-1; the NAc group causes perturbation of a tryptophanyl residue, evidenced by protein difference absorption spectrometry. In the glycosides, a large aglycon likeβpNO₂ Phe orβMeUmb hardly affects the affinity of SBA but a largeN-dansyl group increases the affinity by a factor 20 as compared to GalNAc. The 10-fold increase in carbohydrate-specificN-dansylgalactosamine fluorescence, together with a very favourable entropic contribution point at the presence of a hydrophobic region in the vicinity of the carbohydrate-binding site. The dissociation-rate parameter of the MeUmbβGalNAc SBA complex is slower than for any reported monosaccharide-lectin complex: 0.4 s^-1. The divalent lectin fromErythrina cristagalli preferentially binds the Galβ1 → 4GlcNAc structure that occurs in manyN-glycoproteins. The combining site was mapped thermodynamically with carbohydrates ranging from mono-to pentasaccharides as derived fromN-glycoproteins. Here, N-dansylgalactosamine was used as a fluorescent indicator ligand in substitution titrations. When Galβ1 → 4GlcNAc was linkedα1 → 2 orα1 → 6 to Man, the binding enthalpy and entropy remained practically constant. Application of stopped flow kinetics and pressure-jump relaxation withN-dansylgalactosamine gave mono-exponential signal changes with a concentration dependence corresponding tok ₊ = 4.8 x 10⁴ M^-1 s^-1 k _- = 0.4 to 0.66 s^-1 and a change in reaction volume of+7ml/mol.     

Charge transfer in P-type ATPases investigated on planar membranes

Planar lipid bilayers, e.g., black lipid membranes (BLM) and solid supported membranes (SSM), have been employed to investigate charge movements during the reaction cycle of P-type ATPases. The BLM/SSM method allows a direct measurement of the electrical currents generated by the cation transporter following chemical activation by a substrate concentration jump. The electrical current transients provides information about the reaction mechanism of the enzyme. In particular, the BLM/SSM technique allows identification of electrogenic steps which in turn may be used to localize ion translocation during the reaction cycle of the pump. In addition, using the high time resolution of the technique, especially when rapid activation via caged ATP is employed, rate constants of electrogenic and electroneutral steps can be determined. In the present review, we will discuss the main results obtained by the BLM and SSM methods and how they have contributed to unravel the transport mechanism of P-type ATPases.     

A Survey of λ Repressor Fragments from Two-State to Downhill Folding

We survey the two-state to downhill folding transition by examining 20 λ_6-85^? mutants that cover a wide range of stabilities and folding rates. We investigated four new λ_6-85^? mutants designed to fold especially rapidly. Two were engineered using the core remodeling of Lim and Sauer, and two were engineered using Ferreiro et al.'s frustratometer. These proteins have probe-dependent melting temperatures as high as 80 °C and exhibit a fast molecular phase with the characteristic temperature dependence of the amplitude expected for downhill folding. The survey reveals a correlation between melting temperature and downhill folding previously observed for the β-sheet protein WW domain. A simple model explains this correlation and predicts the melting temperature at which downhill folding becomes possible. An X-ray crystal structure with a 1.64-Å resolution of a fast-folding mutant fragment shows regions of enhanced rigidity compared to the full wild-type protein.