Lack of coupling between secondary structure formation and collapse in a model polypeptide that mimics early folding intermediates, the F2 fragment of the Escherichia coli tryptophan-synthase beta chain.

The isolated, 101-residue long C-terminal (so called F2) fragment of the beta chain from Escherichia coli tryptophan synthase was shown previously to fold into an ensemble of conformations that are condensed, to contain large amounts of highly dynamic secondary structures, and to behave as a good model of structured intermediates that form at the very early stages of protein folding. Here, solvent perturbations were used to investigate the forces that are involved in stabilizing the secondary structure (monitored by far-UV CD) and the condensation of the polypeptide chain (monitored by dynamic light scattering) in isolated F2. It was observed that neither the ionic strength, nor the pH (between 7 and 10), nor salts of the Hofmeister series affected the global secondary structure contents of F2, whereas some of these salts affected the collapse slightly. Addition of trifluoroethanol resulted in a large increase in both the amount of secondary structure and the Stokes radius of F2. Conversely, F2 became more condensed upon raising the temperature from 4 to 60 degrees C, whereas in this temperature range, the secondary structure undergoes significant melting. These observations lead to the conclusion that, in isolated F2, there is no coupling between the hydrophobic collapse and the secondary structure. This finding will be discussed in terms of early events in protein folding.     

Coevolution: Mathematical analysis of host-parasite interactions

An SIS epidemic transmitted by two similar strains of parasite acting on a host population of three genotypes which differ in their reaction to the disease is modelled and analyzed. Singular perturbation techniques are used to reduce the original system of nine differential equations to a coupled system of two equations describing the slowtime coevolution of gene frequency and parasite strain frequency.Karen Christine Beck died June 25, 1983 at home.Born February 8, 1952 in Madison, Wisconsin, She received a B.A. degree in 1974 from Luther College, Decorah, Iowa and a Ph.D. in mathematics in 1980 from the University of Iowa. Since that time she has been an instructor in the Mathematics Department at the University of Utah. She was to become an Assistant Professor at the University of Texas, Arlington, beginning Autumn, 1983. Dr. Beck's areas of specialization in mathematics were Mathematical Analysis and Mathematical Biology. She published numerous research articles that resolved various problems in these areas.     

Heuristics for the Hodgkin–Huxley system

Hodgkin and Huxley (HH) discovered that voltages control ionic currents in nerve membranes. This led them to describe electrical activity in a neuronal membrane patch in terms of an electronic circuit whose characteristics were determined using empirical data. Due to the complexity of this model, a variety of heuristics, including relaxation oscillator circuits and integrate-and-fire models, have been used to investigate activity in neurons, and these simpler models have been successful in suggesting experiments and explaining observations. Connections between most of the simpler models had not been made clear until recently. Shown here are connections between these heuristics and the full HH model. In particular, we study a new model (Type III circuit): It includes the van der Pol-based models; it can be approximated by a simple integrate-and-fire model; and it creates voltages and currents that correspond, respectively, to the h and V components of the HH system.     

Optimization-driven identification of genetic perturbations accelerates the convergence of model parameters in ensemble modeling of metabolic networks

The ensemble modeling (EM) approach has shown promise in capturing kinetic and regulatory effects in the modeling of metabolic networks. Efficacy of the EM procedure relies on the identification of model parameterizations that adequately describe all observed metabolic phenotypes upon perturbation. In this study, we propose an optimization-based algorithm for the systematic identification of genetic/enzyme perturbations to maximally reduce the number of models retained in the ensemble after each round of model screening. The key premise here is to design perturbations that will maximally scatter the predicted steady-state fluxes over the ensemble parameterizations. We demonstrate the applicability of this procedure for an Escherichia coli metabolic model of central metabolism by successively identifying single, double, and triple enzyme perturbations that cause the maximum degree of flux separation between models in the ensemble. Results revealed that optimal perturbations are not always located close to reaction(s) whose fluxes are measured, especially when multiple perturbations are considered. In addition, there appears to be a maximum number of simultaneous perturbations beyond which no appreciable increase in the divergence of flux predictions is achieved. Overall, this study provides a systematic way of optimally designing genetic perturbations for populating the ensemble of models with relevant model parameterizations.     

Relationship between myocardial bridge compression severity and haemodynamic perturbations

Objectives: This study aims to examine the alteration in coronary haemodynamics with increasing the severity of vessel compression caused by myocardial bridging (MB).

Methods: Angiography and intravascular ultrasound were performed in 10 patients with MB with varying severities of systolic compression in the left anterior descending (LAD) artery. Computer models of MB were developed and transient computational fluid dynamics simulations were performed to derive distribution of blood residence time and shear stress.

Results: With increasing the severity of bridge compression, a decreasing trend was observed in the shear stress over proximal segment whereas an increasing trend was found in the shear stress over bridge segment. When patients were divided into 2 groups based on the average systolic vessel compression in the whole cohort (%CR_ave = 27.38), patients with bridges with major systolic compression (>%CR_ave) had smaller shear stress and higher residence time in the proximal segment compared to those with bridges with minor systolic compression (<%CR_ave) (0.37?±?0.23 vs 0.69?±?0.29?Pa and 0.0037?±?0.0069 vs 0.022?±?0.0094?s). In contrast, patients with bridges with major systolic compression had greater shear stress in the bridge segment compared to those with bridges with minor systolic compression (2.49?±?2.06 vs 1.13?±?0.89?Pa). No significant difference was found in the distal shear stress of patients with bridges with major and minor systolic compression.

Conclusion: Our findings revealed a direct relationship between the severity of systolic compression of MB and haemodynamic perturbations in the proximal segment such that the increased systolic vessel compression was associated with decreased shear stress and increased blood residence time.     

Pair potentials for protein folding: choice of reference states and sensitivity of predicted native states to variations in the interaction schemes

We examine the similarities and differences between two widely used knowledge-based potentials, which are expressed as contact matrices (consisting of 210 elements) that gives a scale for interaction energies between the naturally occurring amino acid residues. These are the Miyazawa-Jernigan contact interaction matrix M and the potential matrix S derived by Skolnick J et al., 1997, Protein Sci 6:676-688. Although the correlation between the two matrices is good, there is a relatively large dispersion between the elements. We show that when Thr is chosen as a reference solvent within the Miyazawa and Jernigan scheme, the dispersion between the M and S matrices is reduced. The resulting interaction matrix B gives hydrophobicities that are in very good agreement with experiment. The small dispersion between the S and B matrices, which arises due to differing reference states, is shown to have dramatic effect on the predicted native states of lattice models of proteins. These findings and other arguments are used to suggest that for reliable predictions of protein structures, pairwise additive potentials are not sufficient. We also establish that optimized protein sequences can tolerate relatively large random errors in the pair potentials. We conjecture that three body interaction may be needed to predict the folds of proteins in a reliable manner.     

Thresholds for macroparasite infections

We analyse here the equilibria of an infinite system of partial differential equations modelling the dynamics of a population infected by macroparasites. We find that it is possible to define a reproduction number R₀ that satisfies the intuitive definition, and yields a sharp threshold in the behaviour of the system: if R₀ < 1, the parasite-free equilibrium (PFE) is asymptotically stable and there are no endemic equilibria; if R₀ > 1, the PFE is unstable and there exists a unique endemic equilibrium. The results mainly confirm what had been obtained in simplified models, except for the fact that no backward bifurcation occurs in this model. The stability of equilibria is established by extending an abstract linearization principle and by analysing the spectra of appropriate operators.Revised version: 14 November 2003Supported in part by CNR under Grant n. 00.0142.ST74 Metodi e modelli matematici nello studio dei fenomeni biologici     

Space-time distribution patterns of shape Erica australis L. subsp. shape aragonensis (Willk) after experimental burning,cutting, and ploughing

The response of Erica australis to experimental burning, cutting and ploughing treatments was studied in two heathland communities where it was dominant. The treatments represent those most frequently originated by humans on these heathland communities throughout history. The response of this species in a community where it is dominant and there is no strong interspecific competition was also compared to that produced by it in a community where there is competition. It can be observed that the response to burning and cutting treatments is very similar with very fast spatial occupation in the first few years. From the point of view of time cover values increased in a pronounced manner during the first few years and this increase was stabilized from the fourth year. However, from this moment on a greater increase in this species' maximum height is evident. The response to ploughing is slower according to recovery mechanism (seedlings). Recovery is comparatively less in the area where there is no strong competition than in that where it exists between species. RésuméDans deux communautés de bruyère dominées par Erica australis on a étudié la réponse de cette espèce aux traitements expérimentaux de brûlage, coupe et labour. Ces traitements représentent les actions les plus fréquentes que l'homme a exercé sur ces communautés de bruyère tout au long de l'histoire. De la même façon, on a fait la comparaison de la réponse de cette espèce quand elle se trouve dans une communauté dominée principalement par elle, où elle n'y a pas une forte compétition entre les espèces, avec la réponse quand elle se situe dans une communauté où plusieurs espèces coexistent. On observe que la réponse aux traitements de brûlage et coupe est très semblable, en présentant une occupation spatiale très rapide pendant les premières années. Du point de vue temporel, leurs valeurs de couverture augmentent de façon plus prononcée pendant les premières années, et c'est à partir de la quatrième année que ce développement se stabilise. Cependant, c'est à partir de ce moment là qu'on remarque une augmentation plus grande en ce qui concerne la hauteur maximale de cette espèce. La réponse au traitement de labour est plus lente en relation avec la régénération par semences. L'effet de la compétition se traduit par une récupération qui suit un mode d'augmentation opposé à la situation originaire. Dans la surface où il n'existe pas une forte compétition, la récupération est comparativement plus faible que dans la station où il y existait une compétition entre les espèces.Mots clés: Maquis, Mécanismes reproductives, Perturbations expérimentales, Régénération, Rejets     

Cytochrome c1 exhibits two binding sites for cytochrome c in plants

In plants, channeling of cytochrome c molecules between complexes III and IV has been purported to shuttle electrons within the supercomplexes instead of carrying electrons by random diffusion across the intermembrane bulk phase. However, the mode plant cytochrome c behaves inside a supercomplex such as the respirasome, formed by complexes I, III and IV, remains obscure from a structural point of view. Here, we report ab-initio Brownian dynamics calculations and nuclear magnetic resonance-driven docking computations showing two binding sites for plant cytochrome c at the head soluble domain of plant cytochrome c₁, namely a non-productive (or distal) site with a long heme-to-heme distance and a functional (or proximal) site with the two heme groups close enough as to allow electron transfer. As inferred from isothermal titration calorimetry experiments, the two binding sites exhibit different equilibrium dissociation constants, for both reduced and oxidized species, that are all within the micromolar range, thus revealing the transient nature of such a respiratory complex. Although the docking of cytochrome c at the distal site occurs at the interface between cytochrome c₁ and the Rieske subunit, it is fully compatible with the complex III structure. In our model, the extra distal site in complex III could indeed facilitate the functional cytochrome c channeling towards complex IV by building a “floating boat bridge” of cytochrome c molecules (between complexes III and IV) in plant respirasome.     

Changes in dissolved organic material determine exposure of stream benthic communities to UV-B radiation and heavy metals: implications for climate change

Changes in regional climate in the Rocky Mountains over the next 100 years are expected to have significant effects on biogeochemical cycles and hydrological processes. In particular, decreased discharge and lower stream depth during summer when ultraviolet radiation (UVR) is the highest combined with greater photo-oxidation of dissolved organic materials (DOM) will significantly increase exposure of benthic communities to UVR. Communities in many Rocky Mountain streams are simultaneously exposed to elevated metals from abandoned mines, the toxicity and bioavailability of which are also determined by DOM. We integrated field surveys of 19 streams (21 sites) along a gradient of metal contamination with microcosm and field experiments conducted in Colorado, USA, and New Zealand to investigate the influence of DOM on bioavailability of heavy metals and exposure of benthic communities to UVR. Spatial and seasonal variation in DOM were closely related to stream discharge and significantly influenced heavy metal uptake in benthic organisms. Qualitative and quantitative changes in DOM resulting from exposure to sunlight increased UV-B (290–320 nm) penetration and toxicity of heavy metals. Results of microcosm experiments showed that benthic communities from a metal-polluted stream were tolerant of metals, but were more sensitive to UV-B than communities from a reference stream. We speculate that the greater sensitivity of these communities to UV-B resulted from costs associated with metal tolerance. Exclusion of UVR from 12 separate Colorado streams and from outdoor stream microcosms in New Zealand increased the abundance of benthic organisms (mayflies, stoneflies, and caddisflies) by 18% and 54%, respectively. Our findings demonstrate the importance of considering changes in regional climate and UV-B exposure when assessing the effects of local anthropogenic stressors.