Use of the FACTS solvation model for protein-ligand docking calculations. Application to EADock.

Protein-ligand docking has made important progress during the last decade and has become a powerful tool for drug development, opening the way to virtual high throughput screening and in silico structure-based ligand design. Despite the flattering picture that has been drawn, recent publications have shown that the docking problem is far from being solved, and that more developments are still needed to achieve high successful prediction rates and accuracy. Introducing an accurate description of the solvation effect upon binding is thought to be essential to achieve this goal. In particular, EADock uses the Generalized Born Molecular Volume 2 (GBMV2) solvent model, which has been shown to reproduce accurately the desolvation energies calculated by solving the Poisson equation. Here, the implementation of the Fast Analytical Continuum Treatment of Solvation (FACTS) as an implicit solvation model in small molecules docking calculations has been assessed using the EADock docking program. Our results strongly support the use of FACTS for docking. The success rates of EADock/FACTS and EADock/GBMV2 are similar, i.e. around 75% for local docking and 65% for blind docking. However, these results come at a much lower computational cost: FACTS is 10 times faster than GBMV2 in calculating the total electrostatic energy, and allows a speed up of EADock by a factor of 4. This study also supports the EADock development strategy relying on the CHARMM package for energy calculations, which enables straightforward implementation and testing of the latest developments in the field of Molecular Modeling.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Binding free energy differences in a TCR-peptide-MHC complex induced by a peptide mutation: a simulation analysis.

Recognition by the T-cell receptor (TCR) of immunogenic peptides presented by class I major histocompatibility complexes (MHCs) is the determining event in the specific cellular immune response against virus-infected cells or tumor cells. It is of great interest, therefore, to elucidate the molecular principles upon which the selectivity of a TCR is based. These principles can in turn be used to design therapeutic approaches, such as peptide-based immunotherapies of cancer. In this study, free energy simulation methods are used to analyze the binding free energy difference of a particular TCR (A6) for a wild-type peptide (Tax) and a mutant peptide (Tax P6A), both presented in HLA A2. The computed free energy difference is 2.9 kcal/mol, in good agreement with the experimental value. This makes possible the use of the simulation results for obtaining an understanding of the origin of the free energy difference which was not available from the experimental results. A free energy component analysis makes possible the decomposition of the free energy difference between the binding of the wild-type and mutant peptide into its components. Of particular interest is the fact that better solvation of the mutant peptide when bound to the MHC molecule is an important contribution to the greater affinity of the TCR for the latter. The results make possible identification of the residues of the TCR which are important for the selectivity. This provides an understanding of the molecular principles that govern the recognition. The possibility of using free energy simulations in designing peptide derivatives for cancer immunotherapy is briefly discussed.

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.

New measurements of energy expenditure and physical activity in chronic kidney disease.

The accurate estimation of total daily energy expenditure (TEE) in chronic kidney patients is essential to allow the provision of nutritional requirements; however, it remains a challenge to collect actual physical activity and resting energy expenditure in maintenance dialysis patients. The direct measurement of TEE by direct calorimetry or doubly labeled water cannot be used easily so that, in clinical practice, TEE is usually estimated from resting energy expenditure and physical activity. Prediction equations may also be used to estimate resting energy expenditure; however, their use has been poorly documented in dialysis patients. Recently, a new system called SenseWear Armband (BodyMedia, Pittsburgh, PA) was developed to assess TEE, but so far no data have been published in chronic kidney disease patients. The aim of this review is to describe new measurements of energy expenditure and physical activity in chronic kidney disease patients.

Creation of an Age-Adjusted, Dual-Energy X-ray Absorptiometry-Derived Trabecular Bone Score Curve for the Lumbar Spine in Non-Hispanic US White Women.

The trabecular bone score (TBS, Med-Imaps, Pessac, France) is an index of bone microarchitecture texture extracted from anteroposterior dual-energy X-ray absorptiometry images of the spine. Previous studies have documented the ability of TBS of the spine to differentiate between women with and without fractures among age- and areal bone mineral density (aBMD)-matched controls, as well as to predict future fractures. In this cross-sectional analysis of data collected from 3 geographically dispersed facilities in the United States, we investigated age-related changes in the microarchitecture of lumbar vertebrae as assessed by TBS in a cohort of non-Hispanic US white American women. All subjects were 30 yr of age and older and had an L1-L4aBMDZ-score within ±2 SD of the population mean. Individuals were excluded if they had fractures, were on any osteoporosis treatment, or had any illness that would be expected to impact bone metabolism. All data were extracted from Prodigy dual-energy X-ray absorptiometry devices (GE-Lunar, Madison, WI). Cross-calibrations between the 3 participating centers were performed for TBS and aBMD. aBMD and TBS were evaluated for spine L1-L4 but also for all other possible vertebral combinations. To validate the cohort, a comparison between the aBMD normative data of our cohort and US non-Hispanic white Lunar data provided by the manufacturer was performed. A database of 619 non-Hispanic US white women, ages 30-90 yr, was created. aBMD normative data obtained from this cohort were not statistically different from the non-Hispanic US white Lunar normative data provided by the manufacturer (p = 0.30). This outcome thereby indirectly validates our cohort. TBS values at L1-L4 were weakly inversely correlated with body mass index (r = -0.17) and weight (r = -0.16) and not correlated with height. TBS values for all lumbar vertebral combinations decreased significantly with age. There was a linear decrease of 16.0% (-2.47 T-score) in TBS at L1-L4 between 45 and 90 yr of age (vs. -2.34 for aBMD). Microarchitectural loss rate increased after age 65 by 50% (-0.004 to -0.006). Similar results were obtained for other combinations of lumbar vertebra. TBS, an index of bone microarchitectural texture, decreases with advancing age in non-Hispanic US white women. Little change in TBS is observed between ages 30 and 45. Thereafter, a progressive decrease is observed with advancing age. The changes we observed in these American women are similar to that previously reported for a French population of white women (r(2) > 0.99). This reference database will facilitate the use of TBS to assess bone microarchitectural deterioration in clinical practice.

Enteral versus parenteral nutrition: comparison of energy metabolism in healthy subjects.

Continuous respiratory exchange measurements were performed on 10 healthy young women for 1 h before, 3 h during, and 3 h after either parenteral (iv) or intragastric (ig) administration of a nutrient mixture (52% glucose, 18% amino acid, and 30% lipid energy) infused at twice the postabsorptive resting energy expenditure (REE). REE rose from 0.98 +/- 0.02 (iv) and 0.99 +/- 0.02 kcal/min (ig) postabsorptively to 1.13 +/- 0.03 (iv) and 1.13 +/- 0.02 kcal/min (ig), resulting in nutrient-induced thermogenesis of 10 +/- 0.6 and 9.3 +/- 0.9%, respectively, when related to the metabolizable energy. The respiratory quotient rose from preinfusion values of 0.81 +/- 0.02 (iv) and 0.80 +/- 0.01 (ig) to 0.86 +/- 0.01 (iv) and 0.85 +/- 0.01 (ig). After nutrient administration the respiratory quotient fell significantly to below the preinfusion values. Plasma glucose and insulin concentrations rose during nutrient administration but were higher during the intravenous route. It is concluded that, although the response time to intragastric administration was delayed, the thermic effects and overall substrate oxidations were comparable during intravenous or intragastric administration, albeit, at lower plasma glucose and insulin concentrations via the intragastric route.

Macronutrients and energy balance in obesity.

Energy balance is the difference between metabolizable energy intake and total energy expenditure. Energy intake is difficult to measure accurately; changes in body weight, for example, are not a good measure of the adequacy of energy intake, because fluctuations in body weight are common even if the overall trend is toward weight loss. It is now customary to assess energy requirements indirectly from total energy expenditure. Total energy expenditure consists of basal metabolism, postprandial thermogenesis, and physical activity. Energy expenditure is related to both body weight and body composition. A reduction in total energy expenditure accompanies weight loss, because basal metabolic rate decreases with the loss of lean tissue mass. Similarly, with weight gain, there is an increase in basal metabolic rate, because lean tissue mass grows to support the increase in fat tissue mass. Excess energy intake over energy expenditure causes weight gain and an accompanying increase in total energy expenditure. Following a period of adaptation, total energy expenditure will match energy intake and body weight will stabilize at a higher level. This same relationship holds for weight loss. Respiratory quotient (measured in steady state) is an indication of the proportion of energy expenditure derived from fat and carbohydrate oxidation. Over long periods of time, fat balance is equivalent to energy balance, as an excess of fat intake over fat oxidation causes fat storage.

Social isolation causers mortality by disrupting energy homeostrasis in ants

Social deprivation can have negative effects on the lives of social animals, including humans, yet little is known about the mechanisms by which social withdrawal affects animal health. Here we show that in the carpenter ant Camponotus fellah, socially isolated workers have a greatly reduced life span relative to ants kept in groups of ten individuals. By using a new tracking system, we found that social isolation resulted in important behavioral changes and greatly increased locomotor activity. The higher activity of single ants and their increased propensity to leave the nest to move along the walls suggested that the increased mortality of isolated ants might stem from an imbalance of energy income and expenditure. This view was supported by the finding that while isolated ants ingested the same amount of food as grouped ants, they retained food in the crop, hence preventing its use as an energy source. Moreover, the difference in life span between single and grouped individuals vanished when ants were not fed. This study thus underlines the role of social interactions as key regulators of energy balance, which ultimately affects aging and health in a highly social organism.

Effects of a glucose meal on energy metabolism in patients with cirrhosis before and after liver transplantation.

HYPOTHESIS: Liver transplantation results in hepatic denervation. This may produce alterations of liver energy and substrate metabolism, which may contribute to weight gain after liver transplantation. DESIGN: Prospective clinical study. SETTING: Liver transplantation clinics in a university hospital. PATIENTS: Seven nondiabetic patients with cirrhosis were recruited while on a waiting list for liver transplantation. Seven healthy subjects were recruited as controls. INTERVENTION: Orthotopic liver transplantation. MAIN OUTCOME MEASURES: Evaluation of energy and substrate metabolism after ingestion of a glucose load with indirect calorimetry was performed before, 2 to 6 weeks after, and 5 to 19 months after transplantation. Whole-body glucose oxidation and storage and glucose-induced thermogenesis were calculated. RESULTS: Patients with cirrhosis had modestly elevated resting energy expenditure and normal glucose-induced thermogenesis and postprandial glucose oxidation and storage. These measures remained unchanged after liver transplantation despite a significant increase in postprandial glycemia. Patients, however, gained an average of 3 kg of body weight after 5 to 19 months compared with their weight before transplantation. CONCLUSION: Liver denervation secondary to transplantation does not lead to alterations of energy metabolism after ingestion of a glucose load.