PubMed Health⌕ Search

Biomedical subjects

Marylyn D Ritchie

Publications and source records attributed to Marylyn D Ritchie.

At least 19 recordsLinked to original sources

CASTER-DTA: Equivariant Graph Neural Networks for Predicting Drug-Target Affinity.

Accurately determining the binding affinity of a ligand with a protein is important for drug design, development, and screening. With the advent of accessible protein structure prediction methods such as AlphaFold, predicted protein 3D structures are readily available; however, methods for predicting binding affinity currently do not take full advantage of 3D protein information. Here, we present CASTER-DTA (Cross-Attention with Structural Target Equivariant Representations for Drug-Target Affinity), which uses an equivariant graph neural network to learn more robust protein representations alongside a standard graph neural network to learn molecular representations to predict drug-target affinity. We augment these representations by incorporating an attention-based mechanism between protein residues and drug atoms to improve interpretability. We show that CASTER-DTA represents a state-of-the-art improvement on multiple benchmarks for predicting drug-target affinity and that it generates novel insights for several related tasks. We then apply CASTER-DTA to create a large resource of the binding affinities of every FDA-approved drug against every protein in the human proteome and make these predictions freely available for download. We also make available a web server for researchers to apply a pretrained CASTER-DTA model for predicting binding affinities between arbitrary proteins and drugs.

deep learning↗

Integrating Imaging-Derived Clinical Endotypes with Plasma Proteomics and External Polygenic Risk Scores Enhances Coronary Microvascular Disease Risk Prediction.

Coronary microvascular disease (CMVD) is an underdiagnosed but significant contributor to the burden of ischemic heart disease, characterized by angina and myocardial infarction. The development of risk prediction models such as polygenic risk scores (PRS) for CMVD has been limited by a lack of large-scale genome-wide association studies (GWAS). However, there is significant overlap between CMVD and enrollment criteria for coronary artery disease (CAD) GWAS. In this study, we developed CMVD PRS models by selecting variants identified in a CMVD GWAS and applying weights from an external CAD GWAS, using CMVD-associated loci as proxies for the genetic risk. We integrated plasma proteomics, clinical measures from perfusion PET imaging, and PRS to evaluate their contributions to CMVD risk prediction in comprehensive machine and deep learning models. We then developed a novel unsupervised endotyping framework for CMVD from perfusion PET-derived myocardial blood flow data, revealing distinct patient subgroups beyond traditional case-control definitions. This imaging-based stratification substantially improved classification performance alongside plasma proteomics and PRS, achieving AUROCs between 0.65 and 0.73 per class, significantly outperforming binary classifiers and existing clinical models, highlighting the potential of this stratification approach to enable more precise and personalized diagnosis by capturing the underlying heterogeneity of CMVD. This work represents the first application of imaging-based endotyping and the integration of genetic and proteomic data for CMVD risk prediction, establishing a framework for multimodal modeling in complex diseases.

Cardiovascular Disease↗

Discovery of ancestry-specific variants associated with clopidogrel response among Caribbean Hispanics.

High on-treatment platelet reactivity (HTPR) with clopidogrel predicts ischemic events in adults with coronary artery disease, and while HTPR varies by ethnicity, no genome-wide association study (GWAS) of clopidogrel response has been conducted in Caribbean Hispanics. This study aimed to identify genetic predictors of HTPR in a cohort of 511 Puerto Rican cardiovascular patients treated with clopidogrel, stratified by P2Y12 reaction units (PRU) into responders and non-responders (HTPR). Local ancestry inference (LAI) and traditional GWAS identified variants in the CYP2C19 region associated with HTPR, primarily in individuals with European ancestry. Three variants (OSBPL10 rs1376606, DERL3 rs5030613, RGS6 rs9323567) showed suggestive significance, and a variant in UNC5C was linked to increased HTPR risk. These findings highlight the unique genetic landscape of Caribbean Hispanics and challenge the significance of CYP2C19*2 in predicting clopidogrel response in patients with high non-European ancestry. Further studies are needed to replicate these results in other diverse cohorts.

Journal Article↗

Bidirectional Risk Modulator and Modifier Variant of Dilated and Hypertrophic Cardiomyopathy in BAG3.

IMPORTANCE: The genetic factors that modulate the reduced penetrance and variable expressivity of heritable dilated cardiomyopathy (DCM) are largely unknown. BAG3 genetic variants have been implicated in both DCM and hypertrophic cardiomyopathy (HCM), nominating BAG3 as a gene that harbors potential modifier variants in DCM. OBJECTIVE: To interrogate the clinical traits and diseases associated with BAG3 coding variation. DESIGN, SETTING, AND PARTICIPANTS: This was a cross-sectional study in the Penn Medicine BioBank (PMBB) enrolling patients of the University of Pennsylvania Health System's clinical practice sites from 2014 to 2023. Whole-exome sequencing (WES) was linked to electronic health record (EHR) data to associate BAG3 coding variants with EHR phenotypes. This was a health care population-based study including individuals of European and African genetic ancestry in the PMBB with WES linked to EHR phenotypes, with replication studies in BioVU, UK Biobank, MyCode, and DCM Precision Medicine Study. EXPOSURES: Carrier status for BAG3 coding variants. MAIN OUTCOMES AND MEASURES: Association of BAG3 coding variation with clinical diagnoses, echocardiographic traits, and longitudinal outcomes. RESULTS: In PMBB (n = 43 731; median [IQR] age, 65 [50-76] years; 21 907 female [50.1%]), among 30 324 European and 11 198 African individuals, the common C151R variant was associated with decreased risk for DCM (odds ratio [OR], 0.85; 95% CI, 0.78-0.92) and simultaneous increased risk for HCM (OR, 1.59; 95% CI, 1.25-2.02), which was confirmed in the replication cohorts. C151R carriers exhibited improved longitudinal outcomes compared with noncarriers as assessed by age at death (hazard ratio [HR], 0.85; 95% CI, 0.74-0.96; median [IQR] age, 71.8 [63.1-80.7] in carriers and 70.3 [61.6-79.2] in noncarriers) and heart transplant (HR, 0.81; 95% CI, 0.66-0.99; median [IQR] age, 56.7 [46.1-63.1] in carriers and 55.6 [45.2-62.9] in noncarriers). C151R was associated with reduced risk of DCM (OR, 0.42; 95% CI, 0.24-0.74) and heart failure (OR, 0.27; 95% CI, 0.14-0.50) among individuals harboring truncating TTN variants in exons with high cardiac expression (n = 358). CONCLUSIONS AND RELEVANCE: BAG3 C151R was identified as a bidirectional modulator of risk along the DCM-HCM spectrum, as well as an important genetic modifier variant in TTN-mediated DCM. This work expands on the understanding of the etiology and penetrance of DCM, suggesting that BAG3 C151R is an important genetic modifier variant contributing to the variable expressivity of DCM, warranting further exploration of its mechanisms and of genetic modifiers in DCM more broadly.

Humans↗

Genetic variation in the mitochondrial enzyme carbamyl-phosphate synthetase I predisposes children to increased pulmonary artery pressure following surgical repair of congenital heart defects: a validated genetic association study.

Increased pulmonary artery pressure (PAP) can complicate the postoperative care of children undergoing surgical repair of congenital heart defects. Endogenous NO regulates PAP and is derived from arginine supplied by the urea cycle. The rate-limiting step in the urea cycle is catalyzed by a mitochondrial enzyme, carbamoyl-phosphate synthetase I (CPSI). A well-characterized polymorphism in the gene encoding CPSI (T1405N) has previously been implicated in neonatal pulmonary hypertension. A consecutive modeling cohort of children (N=131) with congenital heart defects requiring surgery was prospectively evaluated to determine key factors associated with increased postoperative PAP, defined as a mean PAP>20 mmHg for at least 1h during the 48h following surgery measured by an indwelling pulmonary artery catheter. Multiple dimensionality reduction (MDR) was used to both internally validate observations and develop optimal two-variable through five-variable models that were tested prospectively in a validation cohort (N=41). Unconditional logistic regression analysis of the modeling cohort revealed that age (OR=0.92, p=0.01), CPSI T1405N genotype (AC vs. AA: OR=4.08, p=0.04, CC vs. AA: OR=5.96, p=0.01), and Down syndrome (OR=5.25, p=0.04) were independent predictors of this complex phenotype. MDR predicted that the best two-variable model consisted of age and CPSI T1405N genotype (p<0.001). This two-variable model correctly predicted 73% of the outcomes from the validation cohort. A five-variable model that added race, gender and Down's syndrome was not significantly better than the two-variable model. In conclusion, the CPSI T1405N genotype appears to be an important new factor in predicting susceptibility to increased PAP following surgical repair of congenital cardiac defects in children.

Carbamoyl-Phosphate Synthase (Ammonia)↗

Renin-angiotensin system gene polymorphisms and coronary artery disease in a large angiographic cohort: detection of high order gene-gene interaction.

There have been many reports regarding the association between renin-angiotensin system (RAS) gene polymorphisms and coronary artery disease (CAD) or acute myocardial infarction (AMI), but the results are inconsistent. In the present study, we used several new approaches with multilocus data to reappraise this issue in a large and relatively homogeneous Taiwanese population. A total of 1254 consecutive patients who underwent cardiac catheterization (735 with documented coronary artery disease and 519 without) between 1996 and 2003 were recruited. Angiotensin-converting enzyme gene insertion/deletion (I/D) polymorphism; T174M, M235T, G-6A, A-20C, G-152A and G-217A polymorphisms of the angiotensinogen gene; and A1166C polymorphism of the angiotensin II type I receptor gene were genotyped. In single-locus analyses, no locus was associated with CAD, history of AMI and three-vessel CAD, either with or without adjustment for conventional CAD risk factors. For multilocus analyses, we recreated a balanced population, with the controls individually matched to the cases regarding the conventional CAD risk factors. We found that the angiotensinogen gene haplotype profile was significantly different between the cases and controls (chi2=31.6, P=0.030) in haplotype analyses. Furthermore, significant three-locus (G-217A, M235T and I/D) gene-gene interactions were detected by multifactor-dimensionality reduction method (highest cross-validation consistency 10.0, lowest prediction error 40.56%, P=0.017) and many even higher order gene-gene interactions by multilocus genotype disequilibrium tests (16 genotype disequilibria exclusively found in the controls, all of which included at least two genes among AGT, ACE and AT1R genes). Our study is the first to demonstrate epistatic, high-order, gene-gene interactions between RAS gene polymorphisms and CAD. These results are compatible with the concept of multilocus and multi-gene effects in complex diseases that would be missed with conventional approaches.

Aged↗

Complex gene-gene interactions in multiple sclerosis: a multifactorial approach reveals associations with inflammatory genes.

The complex inheritance involved in multiple sclerosis (MS) risk has been extensively investigated, but our understanding of MS genetics remains rudimentary. In this study, we explore 51 single nucleotide polymorphisms (SNPs) in 36 candidate genes from the inflammatory pathway and test for gene-gene interactions using complementary case-control, discordant sibling pair, and trio family study designs. We used a sample of 421 carefully diagnosed MS cases and 96 unrelated, healthy controls; discordant sibling pairs from 146 multiplex families; and 275 trio families. We used multifactor dimensionality reduction to explore gene-gene interactions. Based on our analyses, we have identified several statistically significant models including both main effect models and two-locus, three-locus, and four-locus epistasis models that predict MS disease risk with between approximately 61% and 85% accuracy. These results suggest that significant epistasis, or gene-gene interactions, may exist even in the absence of statistically significant individual main effects.

Case-Control Studies↗

Immunogenetics of CD4 lymphocyte count recovery during antiretroviral therapy: An AIDS Clinical Trials Group study.

During antiretroviral therapy, CD4 lymphocyte count increases are modest in some patients despite virologic control. We explored whether polymorphisms in genes important for T cell expansion, survival, and apoptosis are associated with the magnitude of CD4 lymphocyte count recovery during antiretroviral therapy. We studied treatment-naive individuals who achieved sustained control of plasma viremia (<400 HIV-1 RNA copies/mL) for at least 48 weeks after initiation of antiretroviral therapy and compared genotypes among individuals who had an increase of either <200 or > or =200 CD4 cells/mm3 from baseline. A total of 137 single-nucleotide polymorphisms across 17 genes were characterized in 873 study participants. In multivariate analyses that controlled for clinical variables, polymorphisms in genes encoding tumor necrosis factor (TNF)-related apoptosis-inducing ligand (TRAIL), TNF- alpha , Bcl-2-interacting molecule (Bim), interleukin (IL)-15, and IL-15 receptor alpha chain (IL-15R alpha ) were associated with the magnitude of the increase in CD4 lymphocyte count, as were haplotypes in genes encoding interferon- alpha , IL-2, and IL-15R alpha (P < .05, for each). Multifactor dimensionality reduction identified a gene-gene interaction between IL-2/IL-15 receptor common beta chain and IL-2/IL-7/IL-15 receptor common gamma chain. Immune recovery during antiretroviral therapy is a complex phenotype that is influenced by multiple genetic variants. Future studies should validate these tentative associations and define underlying mechanisms.

Adult↗

Drug transporter and metabolizing enzyme gene variants and nonnucleoside reverse-transcriptase inhibitor hepatotoxicity.

This nested case-control study examined relationships between MDR1, CYP2B6, and CYP3A4 variants and hepatotoxicity during antiretroviral therapy with either efavirenz- or nevirapine-containing regimens. Decreased risk of hepatotoxicity was associated with MDR1 3435C-->T (odds ratio, 0.254; P=.021). An interaction between MDR1 and hepatitis B surface antigen status predicted risk with 82% accuracy (P<.001).

Adult↗

Hemochromatosis (HFE) gene mutations and peripheral neuropathy during antiretroviral therapy.

OBJECTIVE: Peripheral neuropathy (PN) often complicates nucleoside reverse transcriptase inhibitor (NRTI) therapy of HIV infection and may involve mitochondrial dysfunction. Since iron deficiency is associated with some types of PN, and iron is essential for mitochondrial function, we tested the hypothesis that hemochromatosis (HFE) gene mutations influence susceptibility to NRTI-induced PN. DESIGN: Case-control study involving multicenter, AIDS Clinical Trials Group (ACTG) protocol 384 and ACTG Human DNA Repository specimens. METHODS: Study participants were randomized to receive three- or four-drug antiretroviral therapy with didanosine (ddI) plus stavudine (d4T) or zidovudine plus lamivudine, given with efavirenz, nelfinavir, or both, with up to three years of follow-up. PN was ascertained based on signs and symptoms. HFE C282Y and H63D genotypes were determined, and associations with PN were assessed using logistic regression. RESULTS: : Of 509 participants, 147 (29%) developed PN, 73% of whom had been randomized to receive ddI plus d4T. Among ddI/d4T-ever-treated individuals, HFE C282Y heterozygotes developed PN on ddI/d4T significantly less often than C282Y non-carriers, adjusting for age, CD4 lymphocyte count and viral load at baseline, and concomitant antiretroviral drugs [6% vs. 35%, respectively, in whites; adjusted odds ratio (OR), 0.17; 95% confidence interval (CI) 0.03-0.83; P = 0.021]. Regardless of race/ethnicity, ddI/d4T-associated PN was uncommon in C282Y heterozygotes [race-adjusted OR, 0.30; 95% CI 0.09-0.96); P = 0.042]. CONCLUSIONS: Iron-loading HFE mutations such as C282Y are associated with a decreased risk of PN during antiretroviral therapy. This finding has potential implications for the prediction and prevention of NRTI-associated PN, particularly in populations at risk of iron deficiency.

Adolescent↗

Parallel multifactor dimensionality reduction: a tool for the large-scale analysis of gene-gene interactions.

UNLABELLED: Parallel multifactor dimensionality reduction is a tool for large-scale analysis of gene-gene and gene-environment interactions. The MDR algorithm was redesigned to allow an unlimited number of study subjects, total variables and variable states, and to remove restrictions on the order of interactions being analyzed. In addition, the algorithm is markedly more efficient, with approximately 150-fold decrease in runtime for equivalent analyses. To facilitate the processing of large datasets, the algorithm was made parallel. AVAILABILITY: Parallel MDR is freely available for non-commercial research institutions. For full details see http://chgr.mc.vanderbilt.edu/ritchielab/pMDR. An open-source version of MDR software is available at http://www.epistasis.org.

Algorithms↗

Single-nucleotide polymorphisms for diagnosis of salt-sensitive hypertension.

BACKGROUND: Salt-sensitive (SS) hypertension affects >30 million Americans and is often associated with low plasma renin activity. We tested the diagnostic validity of several candidate genes for SS and low-renin hypertension. METHODS: In Japanese patients with newly diagnosed, untreated hypertension (n = 184), we studied polymorphisms in 10 genes, including G protein-coupled receptor kinase type 4 (GRK4), some variations of which are associated with hypertension and impair D1 receptor (D1R)-inhibited renal sodium transport. We used the multifactor dimensionality reduction method to determine the genotype associated with salt sensitivity (> or =10% increase in blood pressure with high sodium intake) or low renin. To determine whether the GRK4 genotype is associated with impaired D1R function, we tested the natriuretic effect of docarpamine, a dopamine prodrug, in normotensive individuals with or without GRK4 polymorphisms (n = 18). RESULTS: A genetic model based on GRK4 R65L, GRK4 A142V, and GRK4 A486V was 94.4% predictive of SS hypertension, whereas the single-locus model with only GRK4 A142V was 78.4% predictive, and a 2-locus model of GRK4 A142V and CYP11B2 C-344T was 77.8% predictive of low-renin hypertension. Sodium excretion was inversely related to the number of GRK4 variants in hypertensive persons, and the natriuretic response to dopaminergic stimulation was impaired in normotensive persons having > or =3 GRK4 gene variants. CONCLUSIONS: GRK4 gene variants are associated with SS and low-renin hypertension. However, the genetic model predicting SS hypertension is different from the model for low renin, suggesting genetic differences in these 2 phenotypes. Like low-renin testing, screening for GRK4 variants may be a useful diagnostic adjunct for detection of SS hypertension.

Asian People↗

GPNN: power studies and applications of a neural network method for detecting gene-gene interactions in studies of human disease.

BACKGROUND: The identification and characterization of genes that influence the risk of common, complex multifactorial disease primarily through interactions with other genes and environmental factors remains a statistical and computational challenge in genetic epidemiology. We have previously introduced a genetic programming optimized neural network (GPNN) as a method for optimizing the architecture of a neural network to improve the identification of gene combinations associated with disease risk. The goal of this study was to evaluate the power of GPNN for identifying high-order gene-gene interactions. We were also interested in applying GPNN to a real data analysis in Parkinson's disease. RESULTS: We show that GPNN has high power to detect even relatively small genetic effects (2-3% heritability) in simulated data models involving two and three locus interactions. The limits of detection were reached under conditions with very small heritability (<1%) or when interactions involved more than three loci. We tested GPNN on a real dataset comprised of Parkinson's disease cases and controls and found a two locus interaction between the DLST gene and sex. CONCLUSION: These results indicate that GPNN may be a useful pattern recognition approach for detecting gene-gene and gene-environment interactions.

Algorithms↗

The effect of reduction in cross-validation intervals on the performance of multifactor dimensionality reduction.

Multifactor Dimensionality Reduction (MDR) was developed to detect genetic polymorphisms that present an increased risk of disease. Cross-validation (CV) is an important part of the MDR algorithm, as it prevents over-fitting and allows the predictive ability of a model to be evaluated. CV is a computationally intensive step in the MDR algorithm. Traditionally, MDR has been implemented using 10-fold CV. In order to reduce computation time and therefore allow MDR analysis to be applied to larger datasets, we evaluated the possibility of eliminating or reducing the number of CV intervals used for analysis. We found that eliminating CV made final model selection impossible, but that reducing the number of CV intervals from ten to five caused no loss of power, thereby reducing the computation time of the algorithm by half. The validity of this reduction was confirmed with data from an Alzheimer's disease (AD) study.

Alzheimer Disease↗

Different contributions of polymorphisms in VKORC1 and CYP2C9 to intra- and inter-population differences in maintenance dose of warfarin in Japanese, Caucasians and African-Americans.

OBJECTIVE: To investigate pharmacokinetic and pharmacodynamic factors associated with population differences in warfarin doses needed to achieve anticoagulation, in particular the possible involvement of genetic variability in vitamin K epoxide reductase (VKOR) and CYP2C9. METHODS: Warfarin maintenance dose, unbound plasma S-warfarin concentration [Cu(S)] and INR were determined in 157 Caucasians, 172 Japanese, and 36 African-Americans stably anticoagulated patients. In a subset (n = 166), fully carboxylated plasma normal prothrombin levels (NPT) were also measured. Genotyping for seven CYP2C9 (CYP2C9*1 through 6 and *11) and seven VKORC1 variants were performed in 115 Caucasians and 64 Japanese patients and 66 healthy African-Americans. Multivariate analysis was performed to identify covariates associated with warfarin requirement. RESULTS: The relationship between NPT and Cu(S) indicated Japanese are more susceptible to inhibition of NPT production by S-warfarin than the other two populations. VKORC1 1173 C > T had a greater frequency in Japanese (89.1%) than Caucasians (42.2%) and African-Americans (8.6%). CYP2C9 variants with reduced metabolizing ability were less frequent in Japanese compared to the other two populations. The median warfarin dose was significantly higher in Caucasians than Japanese patients (5.5 versus 3.5 mg/day), however, when matched for CYP2C9*1 homozygosity, no difference in dose was observed between VKORC1 genotype-matched groups. Furthermore, VKORC1 1173C > T and CYP2C9 (*2/*3/*11) genotypes, age and weight were identified as independent covariates contributing to interpatient variability in warfarin dosage. CONCLUSIONS: Both VKORC1 and CYP2C9 polymorphisms contribute to inter-population difference in warfarin doses among the three populations, but their contribution to intra-population variability may differ within each population.

Aryl Hydrocarbon Hydroxylases↗

Multilocus genetic interactions and response to efavirenz-containing regimens: an adult AIDS clinical trials group study.

OBJECTIVE: For the HIV-1 reverse transcriptase inhibitor efavirenz, variant drug transporter gene ABCB1 may predict virologic response but not plasma efavirenz exposure. Conversely, variant drug metabolizing enzyme gene CYP2B6 predicts greater plasma efavirenz exposure but not virologic response. We examined whether long-term responses to efavirenz, and/or plasma efavirenz exposure, are better predicted by multilocus genetic interactions than by individual polymorphisms. MATERIALS AND METHODS: We studied antiretroviral-naïve study participants randomized to receive efavirenz (with or without nelfinavir) plus two nucleoside analogues in study ACTG 384, and who had DNA available for analysis. Participants were followed up for up to 3 years. Nine single nucleotide polymorphisms in ABCB1, CYP2B6, CYP3A4, CYP3A5 and CYP2C19 were identified. Gene-gene interactions were identified using multifactor dimensionality reduction. RESULTS: Among 340 efavirenz recipients, higher efavirenz AUC24 h values were associated with a single locus model involving CYP2B6 516G>T (73% accuracy; P<0.001). This was also the best model among blacks (69% accuracy; P<0.001), whereas among whites the best model involved a gene-gene interaction between CYP2B6 516G>T and ABCB1 2677G>T (82% accuracy, P<0.001). Among 155 participants who received efavirenz without nelfinavir, virologic failure was associated with a two-locus interaction between ABCB1 2677G>T and CYP2B6 516G>T (65% accuracy, P<0.001). Toxicity failure was best predicted by an interaction between ABCB1 2677G>T and ABCB1 3435C>T (71% accuracy, P<0.001). CONCLUSIONS: Multilocus genetic interactions between variant drug metabolism and transporter genes may predict efavirenz pharmacokinetics and treatment responses. This finding may have implications for better individualizing antiretroviral therapy.

Acquired Immunodeficiency Syndrome↗

Estrogens, enzyme variants, and breast cancer: a risk model.

Oxidative metabolites of estrogens have been implicated in the development of breast cancer, yet relatively little is known about the metabolism of estrogens in the normal breast. We developed a mathematical model of mammary estrogen metabolism based on the conversion of 17beta-estradiol (E(2)) by the enzymes cytochrome P450 (CYP) 1A1 and CYP1B1, catechol-O-methyltransferase (COMT), and glutathione S-transferase P1 into eight metabolites [i.e., two catechol estrogens, 2-hydroxyestradiol (2-OHE(2)) and 4-hydroxyestradiol (4-OHE(2)); three methoxyestrogens, 2-methoxyestradiol, 2-hydroxy-3-methoxyestradiol, and 4-methoxyestradiol; and three glutathione (SG)-estrogen conjugates, 2-OHE(2)-1-SG, 2-OHE(2)-4-SG, and 4-OHE(2)-2-SG]. When used with experimentally determined rate constants with purified enzymes, the model provides for a kinetic analysis of the entire metabolic pathway. The predicted concentration of each metabolite during a 30-minute reaction agreed well with the experimentally derived results. The model also enables simulation for the transient quinones, E(2)-2,3-quinone (E(2)-2,3-Q) and E(2)-3,4-quinone (E(2)-3,4-Q), which are not amenable to direct quantitation. Using experimentally derived rate constants for genetic variants of CYP1A1, CYP1B1, and COMT, we used the model to simulate the kinetic effect of enzyme polymorphisms on the pathway and identified those haplotypes generating the largest amounts of catechols and quinones. Application of the model to a breast cancer case-control population identified a subset of women with an increased risk of breast cancer based on their enzyme haplotypes and consequent E(2)-3,4-Q production. This in silico model integrates both kinetic and genomic data to yield a comprehensive view of estrogen metabolomics in the breast. The model offers the opportunity to combine metabolic, genetic, and lifetime exposure data in assessing estrogens as a breast cancer risk factor.

Area Under Curve↗

Multifactor dimensionality reduction: an analysis strategy for modelling and detecting gene-gene interactions in human genetics and pharmacogenomics studies.

The detection of gene-gene and gene-environment interactions associated with complex human disease or pharmacogenomic endpoints is a difficult challenge for human geneticists. Unlike rare, Mendelian diseases that are associated with a single gene, most common diseases are caused by the non-linear interaction of numerous genetic and environmental variables. The dimensionality involved in the evaluation of combinations of many such variables quickly diminishes the usefulness of traditional, parametric statistical methods. Multifactor dimensionality reduction (MDR) is a novel and powerful statistical tool for detecting and modelling epistasis. MDR is a non-parametric and model-free approach that has been shown to have reasonable power to detect epistasis in both theoretical and empirical studies. MDR has detected interactions in diseases such as sporadic breast cancer, multiple sclerosis and essential hypertension. As this method is more frequently applied, and was gained acceptance in the study of human disease and pharmacogenomics, it is becoming increasingly important that the implementation of the MDR approach is properly understood. As with all statistical methods, MDR is only powerful and useful when implemented correctly. Concerns regarding dataset structure, configuration parameters and the proper execution of permutation testing in reference to a particular dataset and configuration are essential to the method's effectiveness. The detection, characterisation and interpretation of gene-gene and gene-environment interactions are expected to improve the diagnosis, prevention and treatment of common human diseases. MDR can be a powerful tool in reaching these goals when used appropriately.

Environment↗