Generalized solvent boundary potential for computer simulations

Copyright 2001 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in The Journal of Chemical Physics and may be found at http://dx.doi.org/10.1063/1.1336570.A general approach has been developed to allow accurate simulations of a small region part of a large macromolecular system while incorporating the influence of the remaining distant atoms with an effective boundary potential. The method is called the Generalized Solvent Boundary Potential (GSBP). By representing the surrounding solvent as a continuum dielectric, both the solvent-shielded static field from the distant atoms of the macromolecule and the reaction field from the dielectricsolvent acting on the atoms in the region of interest are included. The static field is calculated once, using the finite-difference Poisson–Boltzmann (PB) equation, and the result is stored on a discrete grid for efficient simulations. The solventreaction field is developed using a basis-set expansion whose coefficients correspond to generalized electrostatic multipoles. A matrix representing the reaction field Green’s function between those generalized multipoles is calculated only once using the PB equation and stored for efficient simulations. In the present work, the formalism is applied to both spherical and orthorhombic simulation regions for which orthonormal basis-sets exist based on spherical harmonics or cartesian Legendre polynomials. The GSBP method is also tested and illustrated with simple model systems and two detailed atomic systems: the active site region of aspartyl-tRNA synthetase (spherical region) and the interior of the KcsA potassium channel (orthorhombic region). Comparison with numerical finite-difference PB calculations shows that GSBP can accurately describe all long-range electrostatic interactions and remain computationally inexpensive

Electrostatic free energy calculations using the generalized solvent boundary potential method

Copyright 2002 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in The Journal of Chemical Physics and may be found at http://dx.doi.org/10.1063/1.1507108.Free energyperturbation (FEP) calculations using all-atom molecular dynamics simulations with a large number of explicit solvent molecules are a powerful approach to study ligand–macromolecule association processes at the atomic level. One strategy to carry out FEP calculations efficiently and reduce computational time is to consider the explicit dynamics of only a small number of atoms in a localized region around the ligand. Such an approximation is motivated by the observation that the factors governing binding specificity are dominated by interactions in the vicinity of the ligand. However, a straightforward truncation of the system may yield inaccurate results as the influence exerted by the remote regions of the macromolecule and the surrounding solvent through long-range electrostatic effects may be significant. To obtain meaningful results, it is important to incorporate the influence of the remote regions of the ligand–macromolecule complex implicitly using some effective potential. The generalized solvent boundary potential (GSBP) that was developed recently [W. Im, S. Bernèche, and B. Roux, J. Chem. Phys. 114, 2924 (2001)] is an efficient computational method to represent the long-range electrostaticinteractions arising from remote (outer) regions in simulations of a localized (inner) region with a small number of explicit atoms. In the present work, FEP calculations combined with GSBP are used to illustrate the importance of these long-range electrostatic factors in estimation of the charging free energy of an aspartate ligand bound to the aspartyl-tRNA synthetase. Calculations with explicit spherical simulation inner regions of different radii are used to test the accuracy of the GSBP method and also illustrate the importance of explicit protein and solvent dynamics in the free energy estimation. The influence of the represented outer region is tested using separate simulations in which the reaction field and/or the protein static field are excluded. Both components are shown to be essential to obtain quantitatively meaningful results. The ability of implicitly treating the influence of protein fluctuations in the outer region using a protein dielectric constant is examined. It is shown that accurate charging free energy calculations can be performed for this system with a spherical region of 15 to 20 Å radius, which roughly corresponds to 1500–3500 moving atoms. The results indicate that GSBP in combination with FEP calculations is a precise and efficient approach to include long-range electrostatic effects in the study of ligand binding to large macromolecules

GÖÇ VE YAŞAM

Uluslararası Bakalorya Programı, A1 dersi Türk Dili ve Edebiyatı alanında ele alınan bu tezde, Orhan Kemal'in Gurbet Kuşları adlı yapıtında göç olgusu nedenleri ve sonuçlarıyla beraber incelenmiştir. Göç olgusuyla değişen toplumsal yapı, ekonomik ve kültürel farklılıklar çerçevesinde değerlendirilmiştir. Bu tezin amacı, göç olgusunun toplumsal yapıda alt sınıf ve üst sınıflardaki bireyler üzerindeki etkilerini ortaya koymaktır. Üç ana bölümden oluşan tezin ilk bölümünde yapıta adını veren Gurbet Kuşları kavramı üzerinde durulmuştur. Köylülerin aidiyetsizliği ve uyum sorunu bu bölümde aktarılmıştır. Tezin ikinci bölümünde ise köylülerin köyden kente göç sürecinde yaşadıkları kadın ve erkek figürler üzerinden neden ve sonuçlarıyla işlenmiştir. Tezin üçüncü bölümünde şehirliler başlığı altından genel olarak şehirde – İstanbul – yaşayan insanların göç sürecinde köylülerle yaşadıkları uyumsuzluk ve çatışmalara yer verilmektedir. Çalışmada göç sürecinde şehre yerleşen figürlerin şehirlilerle aralarındaki ekonomik ve kültürel farklılıkların sınıflar arasında geçişe olanak tanımadığı sonucuna varılmıştır

All-d-Enantiomer of β-Amyloid Peptide Forms Ion Channels in Lipid Bilayers

Alzheimer’s disease (AD) is the most common type of senile dementia in aging populations. Amyloid β (Aβ)-mediated dysregulation of ionic homeostasis is the prevailing underlying mechanism leading to synaptic degeneration and neuronal death. Aβ-dependent ionic dysregulation most likely occurs either directly via unregulated ionic transport through the membrane or indirectly via Aβ binding to cell membrane receptors and subsequent opening of existing ion channels or transporters. Receptor binding is expected to involve a high degree of stereospecificity. Here, we investigated whether an Aβ peptide enantiomer, whose entire sequence consists of d-amino acids, can form ion-conducting channels; these channels can directly mediate Aβ effects even in the absence of receptor–peptide interactions. Using complementary approaches of planar lipid bilayer (PLB) electrophysiological recordings and molecular dynamics (MD) simulations, we show that the d-Aβ isomer exhibits ion conductance behavior in the bilayer indistinguishable from that described earlier for the l-Aβ isomer. The d isomer forms channel-like pores with heterogeneous ionic conductance similar to the l-Aβ isomer channels, and the d-isomer channel conductance is blocked by Zn2+, a known blocker of l-Aβ isomer channels. MD simulations further verify formation of β-barrel-like Aβ channels with d- and l-isomers, illustrating that both d- and l-Aβ barrels can conduct cations. The calculated values of the single-channel conductance are approximately in the range of the experimental values. These findings are in agreement with amyloids forming Ca2+ leaking, unregulated channels in AD, and suggest that Aβ toxicity is mediated through a receptor-independent, nonstereoselective mechanism

Allostery in Its Many Disguises: From Theory to Applications.

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts

Reorganization Energy for Internal Electron Transfer in Multicopper Oxidases.

We have calculated the reorganization energy for the intramolecular electron transfer between the reduced type 1 copper site and the peroxy intermediate of the trinuclear cluster in the multicopper oxidase CueO. The calculations are performed at the combined quantum mechanics and molecular mechanics (QM/MM) level, based on molecular dynamics simulations with tailored potentials for the two copper sites. We obtain a reorganization energy of 91-133 kJ/mol, depending on the theoretical treatment. The two Cu sites contribute by 12 and 22 kJ/mol to this energy, whereas the solvent contribution is 34 kJ/mol. The rest comes from the protein, involving small contributions from many residues. We have also estimated the energy difference between the two electron-transfer states and show that the reduction of the peroxy intermediate is exergonic by 43-87 kJ/mol, depending on the theoretical method. Both the solvent and the protein contribute to this energy difference, especially charged residues close to the two Cu sites. We compare these estimates with energies obtained from QM/MM optimizations and QM calculations in a vacuum and discuss differences between the results obtained at various levels of theory

Structural basis for hijacking of cellular LxxLL motifs by papillomavirus E6 oncoproteins

E6 viral oncoproteins are key players in epithelial tumors induced by papillomaviruses in vertebrates, including cervical cancer in humans. E6 proteins target many host proteins by specifically interacting with acidic LxxLL motifs. We solved the crystal structures of bovine (BPV1) and human (HPV16) papillomavirus E6 proteins bound to LxxLL peptides from the focal adhesion protein paxillin and the ubiquitin ligase E6AP, respectively. In both E6 proteins, two zinc domains and a linker helix form a basic-hydrophobic pocket, which captures helical LxxLL motifs in a way compatible with other interaction modes. Mutational inactivation of the LxxLL binding pocket disrupts the oncogenic activities of both E6 proteins. This work reveals the structural basis of both the multifunctionality and the oncogenicity of E6 proteins

Coffee and tea consumption in relation to inflammation and basal glucose metabolism in a multi-ethnic Asian population: a cross-sectional study

<p>Abstract</p> <p>Background</p> <p>Higher coffee consumption has been associated with a lower risk of type 2 diabetes in cohort studies, but the physiological pathways through which coffee affects glucose metabolism are not fully understood. The aim of this study was to evaluate the associations between habitual coffee and tea consumption and glucose metabolism in a multi-ethnic Asian population and possible mediation by inflammation.</p> <p>Methods</p> <p>We cross-sectionally examined the association between coffee, green tea, black tea and Oolong tea consumption and glycemic (fasting plasma glucose, HOMA-IR, HOMA-beta, plasma HbA1c) and inflammatory (plasma adiponectin and C-reactive protein) markers in a multi-ethnic Asian population (N = 4139).</p> <p>Results</p> <p>After adjusting for multiple confounders, we observed inverse associations between coffee and HOMA-IR (percent difference: - 8.8% for ≥ 3 cups/day versus rarely or never; <it>P_trend</it>= 0.007), but no significant associations between coffee and inflammatory markers. Tea consumption was not associated with glycemic markers, but green tea was inversely associated with plasma C-reactive protein concentrations (percent difference: - 12.2% for ≥ 1 cup/day versus < 1 cup/week; <it>P_trend</it>= 0.042).</p> <p>Conclusions</p> <p>These data provide additional evidence for a beneficial effect of habitual caffeinated coffee consumption on insulin sensitivity, and suggest that this effect is unlikely to be mediated by anti-inflammatory mechanisms.</p

Estimates of adherence and error analysis of physical activity data collected via accelerometry in a large study of free-living adults

<p>Abstract</p> <p>Background</p> <p>Activity monitors (AM) are small, electronic devices used to quantify the amount and intensity of physical activity (PA). Unfortunately, it has been demonstrated that data loss that occurs when AMs are not worn by subjects (removals during sleeping and waking hours) tend to result in biased estimates of PA and total energy expenditure (TEE). No study has reported the degree of data loss in a large study of adults, and/or the degree to which the estimates of PA and TEE are affected. Also, no study in adults has proposed a methodology to minimize the effects of AM removals.</p> <p>Methods</p> <p>Adherence estimates were generated from a pool of 524 women and men that wore AMs for 13 – 15 consecutive days. To simulate the effect of data loss due to AM removal, a reference dataset was first compiled from a subset consisting of 35 highly adherent subjects (24 HR; minimum of 20 hrs/day for seven consecutive days). AM removals were then simulated during sleep and between one and ten waking hours using this 24 HR dataset. Differences in the mean values for PA and TEE between the 24 HR reference dataset and the different simulations were compared using paired <it>t</it>-tests and/or coefficients of variation.</p> <p>Results</p> <p>The estimated average adherence of the pool of 524 subjects was 15.8 ± 3.4 hrs/day for approximately 11.7 ± 2.0 days. Simulated data loss due to AM removals during sleeping hours in the 24 HR database (n = 35), resulted in biased estimates of PA (p < 0.05), but not TEE. Losing as little as one hour of data from the 24 HR dataset during waking hours results in significant biases (p < 0.0001) and variability (coefficients of variation between 7 and 21%) in the estimates of PA. Inserting a constant value for sleep and imputing estimates for missing data during waking hours significantly improved the estimates of PA.</p> <p>Conclusion</p> <p>Although estimated adherence was good, measurements of PA can be improved by relatively simple imputation of missing AM data.</p