Expanding molecular modeling and design tools to non-natural sidechains.

Protein-protein interactions encode the wiring diagram of cellular signaling pathways and their deregulations underlie a variety of diseases, such as cancer. Inhibiting protein-protein interactions with peptide derivatives is a promising way to develop new biological and therapeutic tools. Here, we develop a general framework to computationally handle hundreds of non-natural amino acid sidechains and predict the effect of inserting them into peptides or proteins. We first generate all structural files (pdb and mol2), as well as parameters and topologies for standard molecular mechanics software (CHARMM and Gromacs). Accurate predictions of rotamer probabilities are provided using a novel combined knowledge and physics based strategy. Non-natural sidechains are useful to increase peptide ligand binding affinity. Our results obtained on non-natural mutants of a BCL9 peptide targeting beta-catenin show very good correlation between predicted and experimental binding free-energies, indicating that such predictions can be used to design new inhibitors. Data generated in this work, as well as PyMOL and UCSF Chimera plug-ins for user-friendly visualization of non-natural sidechains, are all available at http://www.swisssidechain.ch. Our results enable researchers to rapidly and efficiently work with hundreds of non-natural sidechains.

iLOGP: A Simple, Robust, and Efficient Description of n-Octanol/Water Partition Coefficient for Drug Design Using the GB/SA Approach.

The n-octanol/water partition coefficient (log Po/w) is a key physicochemical parameter for drug discovery, design, and development. Here, we present a physics-based approach that shows a strong linear correlation between the computed solvation free energy in implicit solvents and the experimental log Po/w on a cleansed data set of more than 17,500 molecules. After internal validation by five-fold cross-validation and data randomization, the predictive power of the most interesting multiple linear model, based on two GB/SA parameters solely, was tested on two different external sets of molecules. On the Martel druglike test set, the predictive power of the best model (N = 706, r = 0.64, MAE = 1.18, and RMSE = 1.40) is similar to six well-established empirical methods. On the 17-drug test set, our model outperformed all compared empirical methodologies (N = 17, r = 0.94, MAE = 0.38, and RMSE = 0.52). The physical basis of our original GB/SA approach together with its predictive capacity, computational efficiency (1 to 2 s per molecule), and tridimensional molecular graphics capability lay the foundations for a promising predictor, the implicit log P method (iLOGP), to complement the portfolio of drug design tools developed and provided by the SIB Swiss Institute of Bioinformatics.

Comment on recent modeling studies of astrocyte-neuron metabolic interactions.

Recent years have seen a surge in mathematical modeling of the various aspects of neuron-astrocyte interactions, and the field of brain energy metabolism is no exception in that regard. Despite the advent of biophysical models in the field, the long-lasting debate on the role of lactate in brain energy metabolism is still unresolved. Quite the contrary, it has been ported to the world of differential equations. Here, we summarize the present state of this discussion from the modeler's point of view and bring some crucial points to the attention of the non-mathematically proficient reader.

Modeling complex wells with the multi-scale finite-volume method

In this paper, an extension of the multi-scale finite-volume (MSFV) method is devised, which allows to Simulate flow and transport in reservoirs with complex well configurations. The new framework fits nicely into the data Structure of the original MSFV method,and has the important property that large patches covering the whole well are not required. For each well. an additional degree of freedom is introduced. While the treatment of pressure-constraint wells is trivial (the well-bore reference pressure is explicitly specified), additional equations have to be solved to obtain the unknown well-bore pressure of rate-constraint wells. Numerical Simulations of test cases with multiple complex wells demonstrate the ability of the new algorithm to capture the interference between the various wells and the reservoir accurately. (c) 2008 Elsevier Inc. All rights reserved.

The predictive value of trabecular bone score (TBS) on whole lumbar vertebrae mechanics: an ex vivo study.

We investigated the association of trabecular bone score (TBS) with microarchitecture and mechanical behavior of human lumbar vertebrae. We found that TBS reflects vertebral trabecular microarchitecture and is an independent predictor of vertebral mechanics. However, the addition of TBS to areal BMD (aBMD) did not significantly improve prediction of vertebral strength. INTRODUCTION: The trabecular bone score (TBS) is a gray-level measure of texture using a modified experimental variogram which can be extracted from dual-energy X-ray absorptiometry (DXA) images. The current study aimed to confirm whether TBS is associated with trabecular microarchitecture and mechanics of human lumbar vertebrae, and if its combination with BMD improves prediction of fracture risk. METHODS: Lumbar vertebrae (L3) were harvested fresh from 16 donors. The anteroposterior and lateral bone mineral content (BMC) and areal BMD (aBMD) of the vertebral body were measured using DXA; then, the TBS was extracted using TBS iNsight software (Medimaps SA, France). The trabecular bone volume (Tb.BV/tissue volume, TV), trabecular thickness (Tb.Th), degree of anisotropy, and structure model index (SMI) were measured using microcomputed tomography. Quasi-static uniaxial compressive testing was performed on L3 vertebral bodies to assess failure load and stiffness. RESULTS: The TBS was significantly correlated to Tb.BV/TV and SMI (râeuro0/00=âeuro0/000.58 and -0.62; pâeuro0/00=âeuro0/000.02, 0.01), but not related to BMC and BMD. TBS was significantly correlated with stiffness (râeuro0/00=âeuro0/000.64; pâeuro0/00=âeuro0/000.007), independently of bone mass. Using stepwise multiple regression models, we failed to demonstrate that the combination of BMD and TBS was better at explaining mechanical behavior than either variable alone. However, the combination TBS, Tb.Th, and BMC did perform better than each parameter alone, explaining 79 % of the variability in stiffness. CONCLUSIONS: In our study, TBS was associated with microarchitecture parameters and with vertebral mechanical behavior, but TBS did not improve prediction of vertebral biomechanical properties in addition to aBMD.

The predictive value of preoperative B-type natriuretic peptide in children undergooing heart surgery.

Introduction: Low cardiac output syndrome is frequent in childrenafter heart surgery for congenital heart disease and may result in pooroutcome and increased morbidity. In the adult population, preoperativebrain natriuretic peptide (BNP) was shown to be predictive of postoperative complications. In children, the value of preoperative BNP onpostoperative outcome is not so clear. The aim of this study was todetermine the predictive value of preoperative BNP on postoperativeoutcome and low cardiac output syndrome in children after heartsurgery for congenital heart disease.Methods: We examined, retrospectively, the postoperative course of97 pediatric patients (mean age 3.7 years, range 0-14 years old) whounderwent heart surgery in a tertiary care pediatric intensive caresetting. NTproBNP was measured preoperatively in all patients(median 412 pg/ml, range 12-35'000 pg/ml). Patients were divided intothree groups according to their NTproBNP levels (group 1: 0-300 pg/ml, group 2: 300-600 pg/ml, group 3: >600 pg/ml) and then,correlations with postoperative outcomes were examined.Results: We found that patients with a high preoperative BNP requiredmore frequently prolonged (>2 days) mechanical ventilation (33%vs 40% vs 61%, p = 0.045) and stayed more frequently longer than6 days in the intensive care unit (42% vs 50% vs 71%, p = 0.03).However, high preoperative BNP was not correlated with occurrenceof low cardiac output syndrome.Conclusion: Preoperative BNP cannot be used, in children, as areliable and sole predictor of postoperative low cardiac outputsyndrome. However it may help identify, before surgery, those patientsat risk of having a difficult postoperative course.

Carbon and oxygen isotope study of hydrothermal carbonates in the zinc-lead deposits of the San Vicente district, central Peru: A quantitative modeling on mixing processes and CO2 degassing

Carbon and oxygen isotope studies of the host and gangue carbonates of Mississippi Valley-type zinc-lead deposits in the San Vicente District hosted in the Upper Triassic to Lower Jurassic dolostones of the Pucara basin (central Peru) were used to constrain models of the ore formation. A mixing model between an incoming hot saline slightly acidic radiogenic (Pb, Sr) fluid and the native formation water explains the overall isotopic variation (delta(13)C = - 11.5 to + 2.5 parts per thousand relative to PDB and delta(18)O = + 18.0 to + 24.3 parts per thousand relative to SMOW) of the carbonate generations. The dolomites formed during the main ore stage show a narrower range (delta(13)C = - 0.1 to + 1.7 parts per thousand and delta(18)O = + 18.7 to + 23.4 parts per thousand) which is explained by exchange between the mineralizing fluids and the host carbonates combined with changes in temperature and pressure. This model of fluid-rock interaction explains the pervasive alteration of the host dolomite I and precipitation of sphalerite I. The open-space filling hydrothermal white sparry dolomite and the coexisting sphalerite II formed by prolonged fluid-host dolomite interaction and limited CO2 degassing. Late void-filling dolomite III (or calcite) and the associated sphalerite III formed as the consequence of CO2 degassing and concomitant pH increase of a slightly acidic ore fluid. Widespread brecciation is associated to CO2 outgassing. Consequently, pressure variability plays a major role in the ore precipitation during the late hydrothermal events in San Vicente. The presence of native sulfur associated with extremely carbon-light calcites replacing evaporitic sulfates (e.g., delta(13)C = - 11.5 parts per thousand), altered native organic matter and heavier hydrothermal bitumen (from - 27.0 to - 23.0 parts per thousand delta(13)C) points to thermochemical reduction of sulfate and/or thiosulfate. The delta(13)C- and delta(18)O-values of the altered host dolostone and hydrothermal carbonates, and the carbon isotope composition of the associated organic matter show a strong regional homogeneity. These results coupled with the strong mineralogical and petrographic similarities of the different MVT occurrences perhaps reflects the fact that the mineralizing processes were similar in the whole San Vicente belt, suggesting the existence of a common regional mineralizing hydrothermal system with interconnected plumbing.

Identification of optimal structural connectivity using functional connectivity and neural modeling.

The complex network dynamics that arise from the interaction of the brain's structural and functional architectures give rise to mental function. Theoretical models demonstrate that the structure-function relation is maximal when the global network dynamics operate at a critical point of state transition. In the present work, we used a dynamic mean-field neural model to fit empirical structural connectivity (SC) and functional connectivity (FC) data acquired in humans and macaques and developed a new iterative-fitting algorithm to optimize the SC matrix based on the FC matrix. A dramatic improvement of the fitting of the matrices was obtained with the addition of a small number of anatomical links, particularly cross-hemispheric connections, and reweighting of existing connections. We suggest that the notion of a critical working point, where the structure-function interplay is maximal, may provide a new way to link behavior and cognition, and a new perspective to understand recovery of function in clinical conditions.

Physiologically based pharmacokinetic modeling of arsenic in the mouse

A remarkable feature of the carcinogenicity of inorganic arsenic is that while human exposures to high concentrations of inorganic arsenic in drinking water are associated with increases in skin, lung, and bladder cancer, inorganic arsenic has not typically caused tumors in standard laboratory animal test protocols. Inorganic arsenic administered for periods of up to 2 yr to various strains of laboratory mice, including the Swiss CD-1, Swiss CR:NIH(S), C57Bl/6p53(+/-), and C57Bl/6p53(+/+), has not resulted in significant increases in tumor incidence. However, Ng et al. (1999) have reported a 40% tumor incidence in C57Bl/6J mice exposed to arsenic in their drinking water throughout their lifetime, with no tumors reported in controls. In order to investigate the potential role of tissue dosimetry in differential susceptibility to arsenic carcinogenicity, a physiologically based pharmacokinetic (PBPK) model for inorganic arsenic in the rat, hamster, monkey, and human (Mann et al., 1996a, 1996b) was extended to describe the kinetics in the mouse. The PBPK model was parameterized in the mouse using published data from acute exposures of B6C3F1 mice to arsenate, arsenite, monomethylarsonic acid (MMA), and dimethylarsinic acid (DMA) and validated using data from acute exposures of C57Black mice. Predictions of the acute model were then compared with data from chronic exposures. There was no evidence of changes in the apparent volume of distribution or in the tissue-plasma concentration ratios between acute and chronic exposure that might support the possibility of inducible arsenite efflux. The PBPK model was also used to project tissue dosimetry in the C57Bl/6J study, in comparison with tissue levels in studies having shorter duration but higher arsenic treatment concentrations. The model evaluation indicates that pharmacokinetic factors do not provide an explanation for the difference in outcomes across the various mouse bioassays. Other possible explanations may relate to strain-specific differences, or to the different durations of dosing in each of the mouse studies, given the evidence that inorganic arsenic is likely to be active in the later stages of the carcinogenic process. [Authors]

Predictive value of the MGMT promoter methylation status in metastatic melanoma patients receiving first-line temozolomide plus bevacizumab in the trial SAKK 50/07.

The O6-methylguanine-DNA-methyltransferase (MGMT) promoter methylation status is a predictive parameter for the response of malignant gliomas to alkylating agents such as temozolomide. First clinical trials with temozolomide plus bevacizumab therapy in metastatic melanoma patients are ongoing, although the predictive value of the MGMT promoter methylation status in this setting remains unclear. We assessed MGMT promoter methylation in formalin-fixed, primary tumor tissue of metastatic melanoma patients treated with first-line temozolomide and bevacizumab from the trial SAKK 50/07 by methylation-specific polymerase chain reaction. In addition, the MGMT expression levels were also analyzed by MGMT immunohistochemistry. Eleven of 42 primary melanomas (26%) revealed a methylated MGMT promoter. Promoter methylation was significantly associated with response rates CR + PR versus SD + PD according to RECIST (response evaluation criteria in solid tumors) (p<0.05) with a trend to prolonged median progression-free survival (8.1 versus 3.4 months, p>0.05). Immunohistochemically different protein expression patterns with heterogeneous and homogeneous nuclear MGMT expression were identified. Negative MGMT expression levels were associated with overall disease stabilization CR+PR+SD versus PD (p=0.05). There was only a poor correlation between MGMT methylation and lack of MGMT expression. A significant proportion of melanomas have a methylated MGMT promoter. The MGMT promoter methylation status may be a promising predictive marker for temozolomide therapy in metastatic melanoma patients. Larger sample sizes may help to validate significant differences in survival type endpoints.

Modeling the impact of lesions in the human brain.

Lesions of anatomical brain networks result in functional disturbances of brain systems and behavior which depend sensitively, often unpredictably, on the lesion site. The availability of whole-brain maps of structural connections within the human cerebrum and our increased understanding of the physiology and large-scale dynamics of cortical networks allow us to investigate the functional consequences of focal brain lesions in a computational model. We simulate the dynamic effects of lesions placed in different regions of the cerebral cortex by recording changes in the pattern of endogenous ("resting-state") neural activity. We find that lesions produce specific patterns of altered functional connectivity among distant regions of cortex, often affecting both cortical hemispheres. The magnitude of these dynamic effects depends on the lesion location and is partly predicted by structural network properties of the lesion site. In the model, lesions along the cortical midline and in the vicinity of the temporo-parietal junction result in large and widely distributed changes in functional connectivity, while lesions of primary sensory or motor regions remain more localized. The model suggests that dynamic lesion effects can be predicted on the basis of specific network measures of structural brain networks and that these effects may be related to known behavioral and cognitive consequences of brain lesions.

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.