A stitch in time: Efficient computation of genomic DNA melting bubbles

Background: It is of biological interest to make genome-wide predictions of the locations of DNA melting bubbles using statistical mechanics models. Computationally, this poses the challenge that a generic search through all combinations of bubble starts and ends is quadratic. Results: An efficient algorithm is described, which shows that the time complexity of the task is O(NlogN) rather than quadratic. The algorithm exploits that bubble lengths may be limited, but without a prior assumption of a maximal bubble length. No approximations, such as windowing, have been introduced to reduce the time complexity. More than just finding the bubbles, the algorithm produces a stitch profile, which is a probabilistic graphical model of bubbles and helical regions. The algorithm applies a probability peak finding method based on a hierarchical analysis of the energy barriers in the Poland-Scheraga model. Conclusions: Exact and fast computation of genomic stitch profiles is thus feasible. Sequences of several megabases have been computed, only limited by computer memory. Possible applications are the genome-wide comparisons of bubbles with promotors, TSS, viral integration sites, and other melting-related regions.Comment: 16 pages, 10 figure

A mixed methods approach to evaluating community drug distributor performance in the control of neglected tropical diseases

BACKGROUND: Trusted literate, or semi-literate, community drug distributors (CDDs) are the primary implementers in integrated preventive chemotherapy (IPC) programmes for Neglected Tropical Disease (NTD) control. The CDDs are responsible for safely distributing drugs and for galvanising communities to repeatedly, often over many years, receive annual treatment, create and update treatment registers, monitor for side-effects and compile treatment coverage reports. These individuals are 'volunteers' for the programmes and do not receive remuneration for their annual work commitment. METHODS: A mixed methods approach, which included pictorial diaries to prospectively record CDD use of time, structured interviews and focus group discussions, was used to triangulate data on how 58 CDDs allocated their time towards their routine family activities and to NTD Programme activities in Uganda. The opportunity costs of CDD time were valued, performance assessed by determining the relationship between time and programme coverage, and CDD motivation for participating in the programme was explored. RESULTS: Key findings showed approximately 2.5 working weeks (range 0.6-11.4 working weeks) were spent on NTD Programme activities per year. The amount of time on NTD control activities significantly increased between the one and three deliveries that were required within an IPC campaign. CDD time spent on NTD Programme activities significantly reduced time available for subsistence and income generating engagements. As CDDs took more time to complete NTD Programme activities, their treatment performance, in terms of validated coverage, significantly decreased. Motivation for the programme was reported as low and CDDs felt undervalued. CONCLUSIONS: CDDs contribute a considerable amount of opportunity cost to the overall economic cost of the NTD Programme in Uganda due to the commitment of their time. Nevertheless, programme coverage of at least 75 %, as required by the World Health Organisation, is not being achieved and vulnerable individuals may not have access to treatment as a consequence of sub-optimal performance by the CDDs due to workload and programmatic factors

Capturing the essence of folding and functions of biomolecules using Coarse-Grained Models

The distances over which biological molecules and their complexes can function range from a few nanometres, in the case of folded structures, to millimetres, for example during chromosome organization. Describing phenomena that cover such diverse length, and also time scales, requires models that capture the underlying physics for the particular length scale of interest. Theoretical ideas, in particular, concepts from polymer physics, have guided the development of coarse-grained models to study folding of DNA, RNA, and proteins. More recently, such models and their variants have been applied to the functions of biological nanomachines. Simulations using coarse-grained models are now poised to address a wide range of problems in biology.Comment: 37 pages, 8 figure

No effect of 24 h severe energy restriction on appetite regulation and ad libitum energy intake in overweight and obese males

Background/Objectives: Long-term success of weight loss diets might depend on how the appetite regulatory system responds to energy restriction (ER). This study determined the effect of 24 h severe ER on subjective and hormonal appetite regulation, subsequent ad libitum energy intake and metabolism. Subjects/Methods: In randomised order, eight overweight or obese males consumed a 24 h diet containing either 100% (12105 (1174 kJ; energy balance; EB) or 25% (3039 (295) kJ; ER) of estimated daily energy requirements (EER). An individualised standard breakfast containing 25% of EER (3216 (341) kJ) was consumed the following morning and resting energy expenditure, substrate utilisation and plasma concentrations of acylated ghrelin, glucagon-like peptide-1 (GLP-17–36), glucose-dependant insulinotropic peptide (GIP1–42), glucose, insulin and non-esterified fatty acid (NEFA) were determined for 4 h after breakfast. Ad libitum energy intake was assessed in the laboratory on day 2 and via food records on day 3. Subjective appetite was assessed throughout. Results: Energy intake was not different between trials for day 2 (EB: 14946 (1272) kJ; ER: 15251 (2114) kJ; P=0.623), day 3 (EB: 10580 (2457) kJ; 10812 (4357) kJ; P=0.832) or day 2 and 3 combined (P=0.693). Subjective appetite was increased during ER on day 1 (P0.381). Acylated ghrelin, GLP-17–36 and insulin were not different between trials (P>0.104). Post-breakfast area under the curve (AUC) for NEFA (P<0.05) and GIP1–42 (P<0.01) were greater during ER compared with EB. Fat oxidation was greater (P<0.01) and carbohydrate oxidation was lower (P<0.01) during ER, but energy expenditure was not different between trials (P=0.158). Conclusions: These results suggest that 24 h severe ER does not affect appetite regulation or energy intake in the subsequent 48 h. This style of dieting may be conducive to maintenance of a negative EB by limiting compensatory eating behaviour, and therefore may assist with weight loss

Prediction of Optimal Folding Routes of Proteins That Satisfy the Principle of Lowest Entropy Loss: Dynamic Contact Maps and Optimal Control

An optimization model is introduced in which proteins try to evade high energy regions of the folding landscape, and prefer low entropy loss routes during folding. We make use of the framework of optimal control whose convenient solution provides practical and useful insight into the sequence of events during folding. We assume that the native state is available. As the protein folds, it makes different set of contacts at different folding steps. The dynamic contact map is constructed from these contacts. The topology of the dynamic contact map changes during the course of folding and this information is utilized in the dynamic optimization model. The solution is obtained using the optimal control theory. We show that the optimal solution can be cast into the form of a Gaussian Network that governs the optimal folding dynamics. Simulation results on three examples (CI2, Sso7d and Villin) show that folding starts by the formation of local clusters. Non-local clusters generally require the formation of several local clusters. Non-local clusters form cooperatively and not sequentially. We also observe that the optimal controller prefers “zipping” or small loop closure steps during folding. The folding routes predicted by the proposed method bear strong resemblance to the results in the literature

Evolution of protein-coupled RNA dynamics during hierarchical assembly of ribosomal complexes

Assembly of 30S ribosomes involves the hierarchical addition of ribosomal proteins that progressively stabilize the folded 16S rRNA. Here, we use three-color single molecule FRET to show how combinations of ribosomal proteins uS4, uS17 and bS20 in the 16S 5&apos; domain enable the recruitment of protein bS16, the next protein to join the complex. Analysis of real-time bS16 binding events shows that bS16 binds both native and non-native forms of the rRNA. The native rRNA conformation is increasingly favored after bS16 binds, explaining how bS16 drives later steps of 30S assembly. Chemical footprinting and molecular dynamics simulations show that each ribosomal protein switches the 16S conformation and dampens fluctuations at the interface between rRNA subdomains where bS16 binds. The results suggest that specific protein-induced changes in the rRNA dynamics underlie the hierarchy of 30S assembly and simplify the search for the native ribosome structure

Inversion of the balance between hydrophobic and hydrogen bonding interactions in protein folding and aggregation.

Identifying the forces that drive proteins to misfold and aggregate, rather than to fold into their functional states, is fundamental to our understanding of living systems and to our ability to combat protein deposition disorders such as Alzheimer's disease and the spongiform encephalopathies. We report here the finding that the balance between hydrophobic and hydrogen bonding interactions is different for proteins in the processes of folding to their native states and misfolding to the alternative amyloid structures. We find that the minima of the protein free energy landscape for folding and misfolding tend to be respectively dominated by hydrophobic and by hydrogen bonding interactions. These results characterise the nature of the interactions that determine the competition between folding and misfolding of proteins by revealing that the stability of native proteins is primarily determined by hydrophobic interactions between side-chains, while the stability of amyloid fibrils depends more on backbone intermolecular hydrogen bonding interactions

A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry

We present here a review of the fundamental topics of Hartree-Fock theory in Quantum Chemistry. From the molecular Hamiltonian, using and discussing the Born-Oppenheimer approximation, we arrive to the Hartree and Hartree-Fock equations for the electronic problem. Special emphasis is placed in the most relevant mathematical aspects of the theoretical derivation of the final equations, as well as in the results regarding the existence and uniqueness of their solutions. All Hartree-Fock versions with different spin restrictions are systematically extracted from the general case, thus providing a unifying framework. Then, the discretization of the one-electron orbitals space is reviewed and the Roothaan-Hall formalism introduced. This leads to a exposition of the basic underlying concepts related to the construction and selection of Gaussian basis sets, focusing in algorithmic efficiency issues. Finally, we close the review with a section in which the most relevant modern developments (specially those related to the design of linear-scaling methods) are commented and linked to the issues discussed. The whole work is intentionally introductory and rather self-contained, so that it may be useful for non experts that aim to use quantum chemical methods in interdisciplinary applications. Moreover, much material that is found scattered in the literature has been put together here to facilitate comprehension and to serve as a handy reference.Comment: 64 pages, 3 figures, tMPH2e.cls style file, doublesp, mathbbol and subeqn package

Physiological and molecular responses to an acute bout of reduced-exertion high-intensity interval training (REHIT)

PurposeWe have previously shown that 6 weeks of reduced-exertion high-intensity interval training (REHIT) improves V˙O2V˙O2 max in sedentary men and women and insulin sensitivity in men. Here, we present two studies examining the acute physiological and molecular responses to REHIT.MethodsIn Study 1, five men and six women (age: 26 ± 7 year, BMI: 23 ± 3 kg m−2, V˙O2V˙O2 max: 51 ± 11 ml kg−1 min−1) performed a single 10-min REHIT cycling session (60 W and two 20-s ‘all-out’ sprints), with vastus lateralis biopsies taken before and 0, 30, and 180 min post-exercise for analysis of glycogen content, phosphorylation of AMPK, p38 MAPK and ACC, and gene expression of PGC1α and GLUT4. In Study 2, eight men (21 ± 2 year; 25 ± 4 kg·m−2; 39 ± 10 ml kg−1 min−1) performed three trials (REHIT, 30-min cycling at 50 % of V˙O2V˙O2 max, and a resting control condition) in a randomised cross-over design. Expired air, venous blood samples, and subjective measures of appetite and fatigue were collected before and 0, 15, 30, and 90 min post-exercise.ResultsAcutely, REHIT was associated with a decrease in muscle glycogen, increased ACC phosphorylation, and activation of PGC1α. When compared to aerobic exercise, changes in V˙O2V˙O2 , RER, plasma volume, and plasma lactate and ghrelin were significantly more pronounced with REHIT, whereas plasma glucose, NEFAs, PYY, and measures of appetite were unaffected.ConclusionsCollectively, these data demonstrate that REHIT is associated with a pronounced disturbance of physiological homeostasis and associated activation of signalling pathways, which together may help explain previously observed adaptations once considered exclusive to aerobic exercise

Outer membrane protein folding from an energy landscape perspective

The cell envelope is essential for the survival of Gram-negative bacteria. This specialised membrane is densely packed with outer membrane proteins (OMPs), which perform a variety of functions. How OMPs fold into this crowded environment remains an open question. Here, we review current knowledge about OFMP folding mechanisms in vitro and discuss how the need to fold to a stable native state has shaped their folding energy landscapes. We also highlight the role of chaperones and the β-barrel assembly machinery (BAM) in assisting OMP folding in vivo and discuss proposed mechanisms by which this fascinating machinery may catalyse OMP folding