Attracting cavities for docking. Replacing the rough energy landscape of the protein by a smooth attracting landscape.

Molecular docking is a computational approach for predicting the most probable position of ligands in the binding sites of macromolecules and constitutes the cornerstone of structure-based computer-aided drug design. Here, we present a new algorithm called Attracting Cavities that allows molecular docking to be performed by simple energy minimizations only. The approach consists in transiently replacing the rough potential energy hypersurface of the protein by a smooth attracting potential driving the ligands into protein cavities. The actual protein energy landscape is reintroduced in a second step to refine the ligand position. The scoring function of Attracting Cavities is based on the CHARMM force field and the FACTS solvation model. The approach was tested on the 85 experimental ligand-protein structures included in the Astex diverse set and achieved a success rate of 80% in reproducing the experimental binding mode starting from a completely randomized ligand conformer. The algorithm thus compares favorably with current state-of-the-art docking programs. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

NO rebinding to myoglobin: a reactive molecular dynamics study.

The rebinding of NO to myoglobin after photolysis is studied using the 'reactive molecular dynamics' method. In this approach the energy of the system is evaluated on two potential energy surfaces that include the heme-ligand interactions which change between liganded and unliganded myoglobin. This makes it possible to take into account in a simple way, the high dimensionality of the transition seam connecting the reactant and product states. The dynamics of the dissociated NO molecules are examined, and the geometrical and energetic properties of the transition seam are studied. Analysis of the frequency of recrossing shows that the height of the effective rebinding barrier is dependent on the time after photodissociation. This effect is due mainly to protein relaxation and may contribute to the experimentally observed non-exponential rebinding rate of NO, as has been suggested previously.

Local optima networks of hard combinatorial landscapes

Combinatorial optimization involves finding an optimal solution in a finite set of options; many everyday life problems are of this kind. However, the number of options grows exponentially with the size of the problem, such that an exhaustive search for the best solution is practically infeasible beyond a certain problem size. When efficient algorithms are not available, a practical approach to obtain an approximate solution to the problem at hand, is to start with an educated guess and gradually refine it until we have a good-enough solution. Roughly speaking, this is how local search heuristics work. These stochastic algorithms navigate the problem search space by iteratively turning the current solution into new candidate solutions, guiding the search towards better solutions. The search performance, therefore, depends on structural aspects of the search space, which in turn depend on the move operator being used to modify solutions. A common way to characterize the search space of a problem is through the study of its fitness landscape, a mathematical object comprising the space of all possible solutions, their value with respect to the optimization objective, and a relationship of neighborhood defined by the move operator. The landscape metaphor is used to explain the search dynamics as a sort of potential function. The concept is indeed similar to that of potential energy surfaces in physical chemistry. Borrowing ideas from that field, we propose to extend to combinatorial landscapes the notion of the inherent network formed by energy minima in energy landscapes. In our case, energy minima are the local optima of the combinatorial problem, and we explore several definitions for the network edges. At first, we perform an exhaustive sampling of local optima basins of attraction, and define weighted transitions between basins by accounting for all the possible ways of crossing the basins frontier via one random move. Then, we reduce the computational burden by only counting the chances of escaping a given basin via random kick moves that start at the local optimum. Finally, we approximate network edges from the search trajectory of simple search heuristics, mining the frequency and inter-arrival time with which the heuristic visits local optima. Through these methodologies, we build a weighted directed graph that provides a synthetic view of the whole landscape, and that we can characterize using the tools of complex networks science. We argue that the network characterization can advance our understanding of the structural and dynamical properties of hard combinatorial landscapes. We apply our approach to prototypical problems such as the Quadratic Assignment Problem, the NK model of rugged landscapes, and the Permutation Flow-shop Scheduling Problem. We show that some network metrics can differentiate problem classes, correlate with problem non-linearity, and predict problem hardness as measured from the performances of trajectory-based local search heuristics.

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility.

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.

Oligomerization of the alpha 1a- and alpha 1b-adrenergic receptor subtypes. Potential implications in receptor internalization.

We combined biophysical, biochemical, and pharmacological approaches to investigate the ability of the alpha 1a- and alpha 1b-adrenergic receptor (AR) subtypes to form homo- and hetero-oligomers. Receptors tagged with different epitopes (hemagglutinin and Myc) or fluorescent proteins (cyan and green fluorescent proteins) were transiently expressed in HEK-293 cells either individually or in different combinations. Fluorescence resonance energy transfer measurements provided evidence that both the alpha 1a- and alpha 1b-AR can form homo-oligomers with similar transfer efficiency of approximately 0.10. Hetero-oligomers could also be observed between the alpha 1b- and the alpha 1a-AR subtypes but not between the alpha 1b-AR and the beta2-AR, the NK1 tachykinin, or the CCR5 chemokine receptors. Oligomerization of the alpha 1b-AR did not require the integrity of its C-tail, of two glycophorin motifs, or of the N-linked glycosylation sites at its N terminus. In contrast, helix I and, to a lesser extent, helix VII were found to play a role in the alpha 1b-AR homo-oligomerization. Receptor oligomerization was not influenced by the agonist epinephrine or by the inverse agonist prazosin. A constitutively active (A293E) as well as a signaling-deficient (R143E) mutant displayed oligomerization features similar to those of the wild type alpha 1b-AR. Confocal imaging revealed that oligomerization of the alpha1-AR subtypes correlated with their ability to co-internalize upon exposure to the agonist. The alpha 1a-selective agonist oxymetazoline induced the co-internalization of the alpha 1a- and alpha 1b-AR, whereas the alpha 1b-AR could not co-internalize with the NK1 tachykinin or CCR5 chemokine receptors. Oligomerization might therefore represent an additional mechanism regulating the physiological responses mediated by the alpha 1a- and alpha 1b-AR subtypes.

The pollutant diethylhexyl phthalate regulates hepatic energy metabolism via species-specific PPARalpha-dependent mechanisms.

Background: The modulation of energetic homeostasis by pollutants has recently emerged as a potential contributor to the onset of metabolic disorders. Diethylhexyl phthalate (DEHP) is a widely used industrial plasticizer to which humans are widely exposed. Phthalates can activate the three peroxisome proliferatoractivated receptor (PPAR) isotypes on cellular models and induce peroxisome proliferation in rodents.Objectives: In this study, we aimed to evaluate the systemic and metabolic consequences of DEHP exposure that have remained so far unexplored and to characterize the underlying molecular mechanisms of action.Methods: As a proof of concept and mechanism, genetically engineered mouse models of PPARs were exposed to high doses of DEHP, followed by metabolic and molecular analyses.Results: DEHP-treated mice were protected from diet-induced obesity via PPARalpha-dependent activation of hepatic fatty acid catabolism, whereas the activity of neither PPARbeta nor PPARgamma was affected. However, the lean phenotype observed in response to DEHP in wild-type mice was surprisingly abolished in PPARalpha-humanized mice. These species differences are associated with a different pattern of coregulator recruitment.Conclusion: These results demonstrate that DEHP exerts species-specific metabolic actions that rely to a large extent on PPARalpha signaling and highlight the metabolic importance of the species-specific activation of PPARalpha by xenobiotic compounds. Editor's SummaryDiethylhexyl phthalate (DEHP) is an industrial plasticizer used in cosmetics, medical devices, food packaging, and other applications. Evidence that DEHP metabolites can activate peroxisome proliferatoractivated receptors (PPARs) involved in fatty acid oxidation (PPARalpha and PPARbeta) and adiposite function and insulin resistance (PPARgamma) has raised concerns about potential effects of DEHP on metabolic homeostasis. In rodents, PPARalpha activation also induces hepatic peroxisome proliferation, but this response to PPARalpha activation is not observed in humans. Feige et al. (p. 234) evaluated systemic and metabolic consequences of high-dose oral DEHP in combination with a high-fat diet in wild-type mice and genetically engineered mouse PPAR models. The authors report that mice exposed to DEHP gained less weight than controls, without modifying their feeding behavior; they also exhibited lower triglyceride levels, smaller adipocytes, and improved glucose tolerance compared with controls. These effects, which were observed in mice fed both high-fat and standard diets, appeared to be mediated by PPARalpha-dependent activation of hepatic fatty acid catabolism without apparent involvement of PPARbeta or PPARgamma. However, mouse models that expressed human (versus mouse) PPARalpha tended to gain more weight on a high-fat diet than their DHEP-unexposed counterparts. The authors conclude that findings support species-specific metabolic effects of DEHP mediated by PPARalpha activation.

The role of energy and nutritional support in the intensive care unit.

Malnutrition is common in critically ill, hospitalized patients and so represents a major problem for intensive care. Nutritional support can be beneficial in such cases and may help preserve vital organ and immune function. Energy requirements, route of delivery and potential complications of nutritional support are discussed in this paper.

Prevention and treatment of protein energy wasting in chronic kidney disease patients: a consensus statement by the International Society of Renal Nutrition and Metabolism.

Protein energy wasting (PEW) is common in patients with chronic kidney disease (CKD) and is associated with adverse clinical outcomes, especially in individuals receiving maintenance dialysis therapy. A multitude of factors can affect the nutritional and metabolic status of CKD patients requiring a combination of therapeutic maneuvers to prevent or reverse protein and energy depletion. These include optimizing dietary nutrient intake, appropriate treatment of metabolic disturbances such as metabolic acidosis, systemic inflammation, and hormonal deficiencies, and prescribing optimized dialytic regimens. In patients where oral dietary intake from regular meals cannot maintain adequate nutritional status, nutritional supplementation, administered orally, enterally, or parenterally, is shown to be effective in replenishing protein and energy stores. In clinical practice, the advantages of oral nutritional supplements include proven efficacy, safety, and compliance. Anabolic strategies such as anabolic steroids, growth hormone, and exercise, in combination with nutritional supplementation or alone, have been shown to improve protein stores and represent potential additional approaches for the treatment of PEW. Appetite stimulants, anti-inflammatory interventions, and newer anabolic agents are emerging as novel therapies. While numerous epidemiological data suggest that an improvement in biomarkers of nutritional status is associated with improved survival, there are no large randomized clinical trials that have tested the effectiveness of nutritional interventions on mortality and morbidity.

Sex differences in obesity and the regulation of energy homeostasis.

Obesity prevalence is generally higher in women than in men, and there is also a sex difference in body fat distribution. Sex differences in obesity can be explained in part by the influence of gonadal steroids on body composition and appetite; however, behavioural, socio-cultural and chromosomal factors may also play a role. This review, which evolved from the 2008 Stock Conference on sex differences in obesity, summarizes current research and recommendations related to hormonal and neuroendocrine influences on energy balance and fat distribution. A number of important gaps in the research are identified, including a need for more studies on chromosomal sex effects on energy balance, the role of socio-cultural (i.e. gender) factors in obesity and the potential deleterious effects of high-fat diets during pregnancy on the foetus. Furthermore, there is a paucity of clinical trials examining sex-specific approaches and outcomes of obesity treatment (lifestyle-based or pharmacological), and research is urgently needed to determine whether current weight loss programmes, largely developed and tested on women, are appropriate for men. Last, it is important that both animal and clinical research on obesity be designed and analysed in such a way that data can be separately examined in both men and women.

Determinants of brain cell metabolic phenotypes and energy substrate utilization unraveled with a modeling approach.

Although all brain cells bear in principle a comparable potential in terms of energetics, in reality they exhibit different metabolic profiles. The specific biochemical characteristics explaining such disparities and their relative importance are largely unknown. Using a modeling approach, we show that modifying the kinetic parameters of pyruvate dehydrogenase and mitochondrial NADH shuttling within a realistic interval can yield a striking switch in lactate flux direction. In this context, cells having essentially an oxidative profile exhibit pronounced extracellular lactate uptake and consumption. However, they can be turned into cells with prominent aerobic glycolysis by selectively reducing the aforementioned parameters. In the case of primarily oxidative cells, we also examined the role of glycolysis and lactate transport in providing pyruvate to mitochondria in order to sustain oxidative phosphorylation. The results show that changes in lactate transport capacity and extracellular lactate concentration within the range described experimentally can sustain enhanced oxidative metabolism upon activation. Such a demonstration provides key elements to understand why certain brain cell types constitutively adopt a particular metabolic profile and how specific features can be altered under different physiological and pathological conditions in order to face evolving energy demands.

Evaluation of the potential use of trabecular bone score to complement bone mineral density in the diagnosis of osteoporosis: a preliminary spine BMD-matched, case-control study.

The trabecular bone score (TBS) is a new parameter that is determined from gray-level analysis of dual-energy X-ray absorptiometry (DXA) images. It relies on the mean thickness and volume fraction of trabecular bone microarchitecture. This was a preliminary case-control study to evaluate the potential diagnostic value of TBS as a complement to bone mineral density (BMD), by comparing postmenopausal women with and without fractures. The sample consisted of 45 women with osteoporotic fractures (5 hip fractures, 20 vertebral fractures, and 20 other types of fracture) and 155 women without a fracture. Stratification was performed, taking into account each type of fracture (except hip), and women with and without fractures were matched for age and spine BMD. BMD and TBS were measured at the total spine. TBS measured at the total spine revealed a significant difference between the fracture and age- and spine BMD-matched nonfracture group, when considering all types of fractures and vertebral fractures. In these cases, the diagnostic value of the combination of BMD and TBS likely will be higher compared with that of BMD alone. TBS, as evaluated from standard DXA scans directly, potentially complements BMD in the detection of osteoporotic fractures. Prospective studies are necessary to fully evaluate the potential role of TBS as a complementary risk factor for fracture.

Hypertriglyceridemia: a potential side effect of propofol sedation in critical illness.

PURPOSE: Hypertriglyceridemia (hyperTG) is common among intensive care unit (ICU) patients, but knowledge about hyperTG risk factors is scarce. The present study aims to identify risk factors favoring its development in patients requiring prolonged ICU treatment. METHODS: Prospective observational study in the medicosurgical ICU of a university teaching hospital. All consecutive patients staying ≥4 days were enrolled. Potential risk factors were recorded: pathology, energy intake, amount and type of nutritional lipids, intake of propofol, glucose intake, laboratory parameters, and drugs. Triglyceride (TG) levels were assessed three times weekly. Statistics was based on two-way analysis of variance (ANOVA) and linear regression with potential risk factors. RESULTS: Out of 1,301 consecutive admissions, 220 patients were eligible, of whom 99 (45 %) presented hyperTG (triglycerides >2 mmol/L). HyperTG patients were younger, heavier, with more brain injury and multiple trauma. Intake of propofol (mg/kg/h) and lipids' propofol had the highest correlation with plasma TG (r (2) = 0.28 and 0.26, respectively, both p < 0.001). Infection and inflammation were associated with development of hyperTG [C-reactive protein (CRP), r (2) = 0.19, p = 0.004]. No strong association could be found with nutritional lipids or other risk factors. Outcome was similar in normo- and hyperTG patients. CONCLUSIONS: HyperTG is frequent in the ICU but is not associated with adverse outcome. Propofol and accompanying lipid emulsion are the strongest risk factors. Our results suggest that plasma TG should be monitored at least twice weekly in patients on propofol. The clinical consequences of propofol-related hyperTG should be investigated in further studies.

Verbal labels selectively bias brain responses to high-energy foods.

The influence of external factors on food preferences and choices is poorly understood. Knowing which and how food-external cues impact the sensory processing and cognitive valuation of food would provide a strong benefit toward a more integrative understanding of food intake behavior and potential means of interfering with deviant eating patterns to avoid detrimental health consequences for individuals in the long run. We investigated whether written labels with positive and negative (as opposed to 'neutral') valence differentially modulate the spatio-temporal brain dynamics in response to the subsequent viewing of high- and low-energetic food images. Electrical neuroimaging analyses were applied to visual evoked potentials (VEPs) from 20 normal-weight participants. VEPs and source estimations in response to high- and low- energy foods were differentially affected by the valence of preceding word labels over the ~260-300 ms post-stimulus period. These effects were only observed when high-energy foods were preceded by labels with positive valence. Neural sources in occipital as well as posterior, frontal, insular and cingulate regions were down-regulated. These findings favor cognitive-affective influences especially on the visual responses to high-energetic food cues, potentially indicating decreases in cognitive control and goal-adaptive behavior. Inverse correlations between insular activity and effectiveness in food classification further indicate that this down-regulation directly impacts food-related behavior.