Development and validation of a computational method for assessment of missense variants in hypertrophic cardiomyopathy

Assessing the significance of novel genetic variants revealed by DNA sequencing is a major challenge to the integration of genomic techniques with medical practice. Many variants remain difficult to classify by traditional genetic methods. Computational methods have been developed that could contribute to classifying these variants, but they have not been properly validated and are generally not considered mature enough to be used effectively in a clinical setting. We developed a computational method for predicting the effects of missense variants detected in patients with hypertrophic cardiomyopathy (HCM). We used a curated clinical data set of 74 missense variants in six genes associated with HCM to train and validate an automated predictor. The predictor is based on support vector regression and uses phylogenetic and structural features specific to genes involved in HCM. Ten-fold cross validation estimated our predictor's sensitivity at 94% (95% confidence interval: 83%-98%) and specificity at 89% (95% confidence interval: 72%-100%). This corresponds to an odds ratio of 10 for a prediction of pathogenic (95% confidence interval: 4.0-infinity), or an odds ratio of 9.9 for a prediction of benign (95% confidence interval: 4.6-21). Coverage (proportion of variants for which a prediction was made) was 57% (95% confidence interval: 49%-64%). This performance exceeds that of existing methods that are not specifically designed for HCM. The accuracy of this predictor provides support for the clinical use of automated predictions alongside family segregation and population frequency data in the interpretation of new missense variants and suggests future development of similar tools for other diseases.     

Simultaneous gas chromatography-tandem mass spectrometry quantification of symmetric and asymmetric dimethylarginine in human urine

Previously, we demonstrated the utility of a gas chromatography–tandem mass spectrometry (GC–MS/MS) method for the quantitative determination of asymmetric dimethylarginine (ADMA) in biological samples. Here we report the extension of this method to symmetric dimethylarginine (SDMA) in human urine. SDMA and ADMA were simultaneously quantitated in urine by using their in situ prepared trideuteromethyl esters as internal standards. The GC–MS/MS method was validated for SDMA and ADMA in spot urine samples of 19 healthy adults. In these samples, the creatinine-corrected excretion rate was 3.23 ± 0.63 μmol/mmol for SDMA and 3.14 ± 0.98 μmol/mmol for ADMA.     

m-Calpain activation in vitro does not require autolysis or subunit dissociation

Calpains are Ca(2+)-dependent, intracellular cysteine proteases involved in many physiological functions. How calpains are activated in the cell is unknown because the average intracellular concentration of Ca(2+) is orders of magnitude lower than that needed for half-maximal activation of the enzyme in vitro. Two of the proposed mechanisms by which calpains can overcome this Ca(2+) concentration differential are autoproteolysis (autolysis) and subunit dissociation, both of which could release constraints on the core by breaking the link between the anchor helix and the small subunit to allow the active site to form. By measuring the rate of autolysis at different sites in calpain, we show that while the anchor helix is one of the first targets to be cut, this occurs in the same time-frame as several potentially inactivating cleavages in Domain III. Thus autolytic activation would overlap with inactivation. We also show that the small subunit does not dissociate from the large subunit, but is proteolyzed to a 40-45k heterodimer of Domains IV and VI. It is likely that this autolysis-generated heterodimer has previously been misidentified as the small subunit homodimer produced by subunit dissociation. We propose a model for m-calpain activation that does not involve either autolysis or subunit dissociation.     

Kinetic mechanism of an aldehyde reductase of Saccharomyces cerevisiae that relieves toxicity of furfural and 5-hydroxymethylfurfural

An effective means of relieving the toxicity of furan aldehydes, furfural (FFA) and 5-hydroxymethylfurfural (HMF), on fermenting organisms is essential for achieving efficient fermentation of lignocellulosic biomass to ethanol and other products. Ari1p, an aldehyde reductase from Saccharomyces cerevisiae, has been shown to mitigate the toxicity of FFA and HMF by catalyzing the NADPH-dependent conversion to corresponding alcohols, furfuryl alcohol (FFOH) and 5-hydroxymethylfurfuryl alcohol (HMFOH). At pH 7.0 and 25°C, purified Ari1p catalyzes the NADPH-dependent reduction of substrates with the following values (k(cat) (s(-1)), k(cat)/K(m) (s(-1)mM(-1)), K(m) (mM)): FFA (23.3, 1.82, 12.8), HMF (4.08, 0.173, 23.6), and dl-glyceraldehyde (2.40, 0.0650, 37.0). When acting on HMF and dl-glyceraldehyde, the enzyme operates through an equilibrium ordered kinetic mechanism. In the physiological direction of the reaction, NADPH binds first and NADP(+) dissociates from the enzyme last, demonstrated by k(cat) of HMF and dl-glyceraldehyde that are independent of [NADPH] and (K(ia)(NADPH)/k(cat)) that extrapolate to zero at saturating HMF or dl-glyceraldehyde concentration. Microscopic kinetic parameters were determined for the HMF reaction (HMF+NADPH?HMFOH+NADP(+)), by applying steady-state, presteady-state, kinetic isotope effects, and dynamic modeling methods. Release of products, HMFOH and NADP(+), is 84% rate limiting to k(cat) in the forward direction. Equilibrium constants, [NADP(+)][FFOH]/[NADPH][FFA][H(+)]=5600×10(7)M(-1) and [NADP(+)][HMFOH]/[NADPH][HMF][H(+)]=4200×10(7)M(-1), favor the physiological direction mirrored by the slowness of hydride transfer in the non-physiological direction, NADP(+)-dependent oxidation of alcohols (k(cat) (s(-1)), k(cat)/K(m) (s(-1)mM(-1)), K(m) (mM)): FFOH (0.221, 0.00158, 140) and HMFOH (0.0105, 0.000104, 101).     

Thermodynamics and solvent effects on substrate and cofactor binding in Escherichia coli chromosomal dihydrofolate reductase

Chromosomal dihydrofolate reductase from Escherichia coli catalyzes the reduction of dihydrofolate to tetrahydrofolate using NADPH as a cofactor. The thermodynamics of ligand binding were examined using an isothermal titration calorimetry approach. Using buffers with different heats of ionization, zero to a small, fractional proton release was observed for dihydrofolate binding, while a proton was released upon NADP(+) binding. The role of water in binding was additionally monitored using a number of different osmolytes. Binding of NADP(+) is accompanied by the net release of ～5-24 water molecules, with a dependence on the identity of the osmolyte. In contrast, binding of dihydrofolate is weakened in the presence of osmolytes, consistent with "water uptake". Different effects are observed depending on the identity of the osmolyte. The net uptake of water upon dihydrofolate binding was previously observed in the nonhomologous R67-encoded dihydrofolate reductase (dfrB or type II enzyme) [Chopra, S., et al. (2008) J. Biol. Chem. 283, 4690-4698]. As R67 dihydrofolate reductase possesses a nonhomologous sequence and forms a tetrameric structure with a single active site pore, the observation of weaker DHF binding in the presence of osmolytes in both enzymes implicates cosolvent effects on free dihydrofolate. Consistent with this analysis, stopped flow experiments find betaine mostly affects DHF binding via changes in k(on), while betaine mostly affects NADPH binding via changes in k(off). Finally, nonadditive enthalpy terms when binary and ternary cofactor binding events are compared suggest the presence of long-lived conformational transitions that are not included in a simple thermodynamic cycle.