Kim A. Sharp

The Common Feature of Leukemia-Associated IDH1 and IDH2 Mutations Is a Neomorphic Enzyme Activity Converting α-Ketoglutarate to 2-Hydroxyglutarate

The somatic mutations in cytosolic isocitrate dehydrogenase 1 (IDH1) observed in gliomas can lead to the production of 2-hydroxyglutarate (2HG). Here, we report that tumor 2HG is elevated in a high percentage of patients with cytogenetically normal acute myeloid leukemia (AML). Surprisingly, less than half of cases with elevated 2HG possessed IDH1 mutations. The remaining cases with elevated 2HG had mutations in IDH2, the mitochondrial homolog of IDH1. These data demonstrate that a shared feature of all cancer-associated IDH mutations is production of the oncometabolite 2HG. Furthermore, AML patients with IDH mutations display a significantly reduced number of other well characterized AML-associated mutations and/or associated chromosomal abnormalities, potentially implicating IDH mutation in a distinct mechanism of AML pathogenesis.

Entropy—enthalpy compensation: Fact or artifact?

The phenomenon of entropy–enthalpy (S‐H) compensation is widely invoked as an explanatory principle in thermodynamic analyses of proteins, ligands, and nucleic acids. It has been suggested that this compensation is an intrinsic property of either complex, fluctuating, or aqueous systems. The questions examined here are whether the observed compensation is extra‐thermodynamic (i.e., reflects anything more than the well‐known laws of statistical thermodynamics) and if so, what does it reveal about the system? Compensation is rather variably defined in the literature and different usages are discussed. The most precise and interesting one, which is considered here, is a linear relationship between ΔH and ΔS for some series of perturbations or changes in experimental variable. Some recent thermodynamic data on proteins purporting to show compensation is analyzed and shown to be better explained by other causes. A general statistical mechanical model of a complex system is analyzed to explore whether and under what conditions extra‐thermodynamic compensation can occur and what it reveals about the system. This model shows that the most likely behavior to be seen is linear S‐H compensation over a rather limited range of perturbations with a compensation temperature Tc = dΔH/dΔS within 20% of the experimental temperature. This behavior is insensitive to the details of the model, thus revealing little extra‐thermodynamic or causal information about the system. In addition, it will likely be difficult to distinguish this from more trivial forms of compensation in real experimental systems.

Theory of the electrokinetic behavior of human erythrocytes.

We develop a theory of electrophoresis of human erythrocytes that predicts mobilities significantly smaller than those based on the classical Smoluchowski relation. In the classical treatment the charge is assumed to be spread uniformly on the hydrodynamic surface. The present model takes into account that most of the charge, due mainly to sialic acid, is contained in the glycocalyx. The glycocalyx is modeled as a permeable layer of polyelectrolyte molecules anchored to the cell membrane. The charge is assumed to be uniformly distributed throughout this layer. The fluid flow in the layer is treated as being dominated by Stokes friction arising from idealized polymer segments. The Navier-Stokes equations are solved to give the dependence of electroosomotic velocity with distance from the cell surface. An expression for the electrophoretic mobility is obtained which contains two parameters (a) the thickness of the glycocalyx and (b) the mean polymer segment radius. The best fit to experimental data is obtained if these are given the values 75 A and 7 A, respectively. Deviation from experimental data at low ionic strength (less than 0.05 M) occurs. However, this deviation is in the direction one would expect if at low ionic strength the polyelectrolyte layer expands slightly due to decreased charge shielding.

The electrostatic potential of B‐DNA

Electrostatic potentials around DNA are obtained by solving the nonlinear Poisson‐Boltzmann (PB) equation. The detailed charge distribution of the DNA and the different polarizabilities of the macromolecule and solvent are included explicitly in the calculations. The PB equation is solved using extensions of a finite difference approach applied previously to proteins. Electrical potentials and ion concentrations are compared to those obtained with simpler models. It is found that the shape of the dielectric boundary between the macromolecule and solvent has significant effects on the calculated potentials near the surface, particularly in the grooves. Sequence‐specific patterns are found, the most surprising result being the existence of positive regions of potential near the bases in both the major and minor grooves. The effect of solvent and ionic atmosphere screening of phosphate‐phosphate repulsions is studied, and an effective dielectric function, appropriate for molecular mechanics simulations, is derived.

Entropy in protein folding and in protein-protein interactions.

The reduction of conformational entropy is a major barrier that has to be overcome in protein folding and binding. Changes in solvent entropy are also a major factor. Recent advances include clarification of the fundamental issues concerning the separation of entropy into components, the treatment of association entropy in binding, and the role of size and shape effects in solvation entropy. Advances in the application of entropy calculations include an emerging consensus for estimates of backbone and sidechain entropy loss in protein folding via use of numerically intensive methods for sampling, and use of the expanding protein-structure database.

Water structure changes induced by hydrophobic and polar solutes revealed by simulations and infrared spectroscopy

A combination of simulations and Fourier transform infrared spectroscopy was used to examine the effect of three ionic solutes (KCl, NaCl, and KSCN), the polar solute urea, and the osmolyte trimethylamine-N-oxide (TMAO) on a water structure. The ionic solutes increase the mean water–water H-bond angle in their first hydration shell concomitantly shifting the OH stretching mode to higher frequency, and shifting the HOH bending mode to lower frequency. TMAO decreases the mean water–water H-bond angle in its first hydration shell, shifts the OH stretching mode frequency down, and shifting the HOH bending mode frequency up. Urea has no effect on the mean H-bond angle, OH stretch, and HOH bend frequencies. These results can be explained in terms of changes in the relative proportions of two H-bond angle populations: Ionic solutes increase the population of more distorted (larger angle) H bonds relative to the less distorted population, TMAO has the reverse effect, while urea does not affect the H-bond angle probability distribution. The negligible effect of urea on water structure supports the direct binding model for urea-induced protein denaturation.A combination of simulations and Fourier transform infrared spectroscopy was used to examine the effect of three ionic solutes (KCl, NaCl, and KSCN), the polar solute urea, and the osmolyte trimethylamine-N-oxide (TMAO) on a water structure. The ionic solutes increase the mean water–water H-bond angle in their first hydration shell concomitantly shifting the OH stretching mode to higher frequency, and shifting the HOH bending mode to lower frequency. TMAO decreases the mean water–water H-bond angle in its first hydration shell, shifts the OH stretching mode frequency down, and shifting the HOH bending mode frequency up. Urea has no effect on the mean H-bond angle, OH stretch, and HOH bend frequencies. These results can be explained in terms of changes in the relative proportions of two H-bond angle populations: Ionic solutes increase the population of more distorted (larger angle) H bonds relative to the less distorted population, TMAO has the reverse effect, while urea does not affect the H-bond angle pr...

Identification of additional IDH mutations associated with oncometabolite R (−)-2-hydroxyglutarate production

Mutations in cytosolic isocitrate dehydrogenase 1 (IDH1) or its mitochondrial homolog IDH2 can lead to R(−)-2-hydroxyglutarate (2HG) production. To date, mutations in three active site arginine residues, IDH1 R132, IDH2 R172 and IDH2 R140, have been shown to result in the neomorphic production of 2HG. Here we report on three additional 2HG-producing IDH1 mutations: IDH1 R100, which is affected in adult glioma, IDH1 G97, which is mutated in colon cancer cell lines and pediatric glioblastoma, and IDH1 Y139. All these new mutants stereospecifically produced 2HGs (R) enantiomer. In contrast, we find that the IDH1 SNPs V71I and V178I, as well as a number of other single-sample reports of IDH non-synonymous mutation, did not elevate cellular 2HG levels in cells and retained the wild-type ability for isocitrate-dependent NADPH production. Finally, we report the existence of additional rare, but recurring mutations found in lymphoma and thyroid cancer, which while failing to elevate 2HG nonetheless displayed loss of function, indicating a possible tumorigenic mechanism for a non-2HG-producing subset of IDH mutations in some malignancies. These data broaden our understanding of how IDH mutations may contribute to cancer through either neomorphic R(−)-2HG production or reduced wild-type enzymatic activity, and highlight the potential value of metabolite screening in identifying IDH-mutated tumors associated with elevated oncometabolite levels.

On the calculation of absolute macromolecular binding free energies

The standard framework for calculating the absolute binding free energy of a macromolecular association reaction A + B → AB with an association constant KAB is to equate chemical potentials of the species on the left- and right-hand sides of this reaction and evaluate the chemical potentials from theory. This theory involves (usually hidden) assumptions about what constitutes the bound species, AB, and where the contribution of the solvent appears. We present here an alternative derivation that can be traced back to Bjerrum, in which the expectation value of KAB is obtained directly through the statistical mechanical method of evaluating its ensemble (Boltzmann-weighted) average. The generalized Bjerrum approach more clearly delineates: (i) the different contributions to binding; (ii) the origin of the much-discussed and somewhat controversial association entropy term; and (iii) where the solvent contribution appears. This approach also allows approximations required for practical evaluation of the binding constant in complex macromolecular systems, to be introduced in a well defined way. We provide an example, with application to test cases that illustrate a range of binding behavior.

Calculations of the electrostatic potential adjacent to model phospholipid bilayers.

We used the nonlinear Poisson-Boltzmann equation to calculate electrostatic potentials in the aqueous phase adjacent to model phospholipid bilayers containing mixtures of zwitterionic lipids (phosphatidylcholine) and acidic lipids (phosphatidylserine or phosphatidylglycerol). The aqueous phase (relative permittivity, epsilon r = 80) contains 0.1 M monovalent salt. When the bilayers contain < 11% acidic lipid, the -25 mV equipotential surfaces are discrete domes centered over the negatively charged lipids and are approximately twice the value calculated using Debye-Hückel theory. When the bilayers contain > 25% acidic lipid, the -25 mV equipotential profiles are essentially flat and agree well with the values calculated using Gouy-Chapman theory. When the bilayers contain 100% acidic lipid, all of the equipotential surfaces are flat and agree with Gouy-Chapman predictions (including the -100 mV surface, which is located only 1 A from the outermost atoms). Even our model bilayers are not simple systems: the charge on each lipid is distributed over several atoms, these partial charges are non-coplanar, there is a 2 A ion-exclusion region (epsilon r = 80) adjacent to the polar headgroups, and the molecular surface is rough. We investigated the effect of these four factors using smooth (or bumpy) epsilon r = 2 slabs with embedded point charges: these factors had only minor effects on the potential in the aqueous phase.

Calculation of the Electrophoretic Mobility of a Particle Bearing Bound Polyelectrolyte Using the Nonlinear Poisson-Boltzmann Equation

A numerical method for determining the electrophoretic mobility of a polyelectrolyte-coated particle is presented. The particle surface is modeled as having a permeable layer of polyelectrolyte molecules anchored to its surface. Fluid flow within the polyelectrolyte layer is subject to Stokes drag arising from the polyelectrolyte segments. The method allows arbitrary distribution of polymer segments and charge density normal to the surface to be used. The hydrodynamic plane of shear may also be varied. The potential profile is determined by a numerical solution to the nonlinearized Poisson-Boltzmann equation. The potential profile is then used in a numerical solution to the Navier-Stokes equation to give the required mobility. The use of the nonlinearized Poisson-Boltzmann equation extends the results to higher charge density/lower ionic strength conditions than previous treatments. The surface potentials and mobilities for three limiting charge distributions are compared for both the linear and nonlinear treatments to delimit the range of validity of the linear treatment. The utility of the numerical, nonlinear treatment is demonstrated by an improved fit to the electrophoretic mobility of human erythrocytes as a function of ionic strength in the range 10 to 150 mM.