skip to main content

Explore Research Products in the NSF-PAR

Title: Dynamic coupling of residues within proteins as a mechanistic foundation of many enigmatic pathogenic missense variants
Many pathogenic missense mutations are found in protein positions that are neither well-conserved nor fall in any known functional domains. Consequently, we lack any mechanistic underpinning of dysfunction caused by such mutations. We explored the disruption of allosteric dynamic coupling between these positions and the known functional sites as a possible mechanism for pathogenesis. In this study, we present an analysis of 591 pathogenic missense variants in 144 human enzymes that suggests that allosteric dynamic coupling of mutated positions with known active sites is a plausible biophysical mechanism and evidence of their functional importance. We illustrate this mechanism in a case study of β-Glucocerebrosidase (GCase) in which a vast majority of 94 sites harboring Gaucher disease-associated missense variants are located some distance away from the active site. An analysis of the conformational dynamics of GCase suggests that mutations on these distal sites cause changes in the flexibility of active site residues despite their distance, indicating a dynamic communication network throughout the protein. The disruption of the long-distance dynamic coupling caused by missense mutations may provide a plausible general mechanistic explanation for biological dysfunction and disease.  more » « less
Award ID(s):
1934848
NSF-PAR ID:
10354817
Author(s) / Creator(s):
; ; ; ; ; ;
Editor(s):
Wallqvist, Anders
Date Published:
Journal Name:
PLOS Computational Biology
Volume:
18
Issue:
4
ISSN:
1553-7358
Page Range / eLocation ID:
e1010006
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ; ; ; et al ( , Circulation)
    null (Ed.)
    Background: Hypertrophic cardiomyopathy (HCM) is caused by pathogenic variants in sarcomere protein genes that evoke hypercontractility, poor relaxation, and increased energy consumption by the heart and increased patient risks for arrhythmias and heart failure. Recent studies show that pathogenic missense variants in myosin, the molecular motor of the sarcomere, are clustered in residues that participate in dynamic conformational states of sarcomere proteins. We hypothesized that these conformations are essential to adapt contractile output for energy conservation and that pathophysiology of HCM results from destabilization of these conformations. Methods: We assayed myosin ATP binding to define the proportion of myosins in the super relaxed state (SRX) conformation or the disordered relaxed state (DRX) conformation in healthy rodent and human hearts, at baseline and in response to reduced hemodynamic demands of hibernation or pathogenic HCM variants. To determine the relationships between myosin conformations, sarcomere function, and cell biology, we assessed contractility, relaxation, and cardiomyocyte morphology and metabolism, with and without an allosteric modulator of myosin ATPase activity. We then tested whether the positions of myosin variants of unknown clinical significance that were identified in patients with HCM, predicted functional consequences and associations with heart failure and arrhythmias. Results: Myosins undergo physiological shifts between the SRX conformation that maximizes energy conservation and the DRX conformation that enables cross-bridge formation with greater ATP consumption. Systemic hemodynamic requirements, pharmacological modulators of myosin, and pathogenic myosin missense mutations influenced the proportions of these conformations. Hibernation increased the proportion of myosins in the SRX conformation, whereas pathogenic variants destabilized these and increased the proportion of myosins in the DRX conformation, which enhanced cardiomyocyte contractility, but impaired relaxation and evoked hypertrophic remodeling with increased energetic stress. Using structural locations to stratify variants of unknown clinical significance, we showed that the variants that destabilized myosin conformations were associated with higher rates of heart failure and arrhythmias in patients with HCM. Conclusions: Myosin conformations establish work-energy equipoise that is essential for life-long cellular homeostasis and heart function. Destabilization of myosin energy-conserving states promotes contractile abnormalities, morphological and metabolic remodeling, and adverse clinical outcomes in patients with HCM. Therapeutic restabilization corrects cellular contractile and metabolic phenotypes and may limit these adverse clinical outcomes in patients with HCM. 
    more » « less
  2. ; ; ; ; ; ; ( , eLife)

    We develop integrated co-evolution and dynamic coupling (ICDC) approach to identify, mutate, and assess distal sites to modulate function. We validate the approach first by analyzing the existing mutational fitness data of TEM-1 β-lactamase and show that allosteric positions co-evolved and dynamically coupled with the active site significantly modulate function. We further apply ICDC approach to identify positions and their mutations that can modulate binding affinity in a lectin, cyanovirin-N (CV-N), that selectively binds to dimannose, and predict binding energies of its variants through Adaptive BP-Dock. Computational and experimental analyses reveal that binding enhancing mutants identified by ICDC impact the dynamics of the binding pocket, and show that rigidification of the binding residues compensates for the entropic cost of binding. This work suggests a mechanism by which distal mutations modulate function through dynamic allostery and provides a blueprint to identify candidates for mutagenesis in order to optimize protein function.

     
    more » « less
  3. ; ( , Entropy)
    Understanding the underlying mechanisms behind protein allostery and non-additivity of substitution outcomes (i.e., epistasis) is critical when attempting to predict the functional impact of mutations, particularly at non-conserved sites. In an effort to model these two biological properties, we extend the framework of our metric to calculate dynamic coupling between residues, the Dynamic Coupling Index (DCI) to two new metrics: (i) EpiScore, which quantifies the difference between the residue fluctuation response of a functional site when two other positions are perturbed with random Brownian kicks simultaneously versus individually to capture the degree of cooperativity of these two other positions in modulating the dynamics of the functional site and (ii) DCIasym, which measures the degree of asymmetry between the residue fluctuation response of two sites when one or the other is perturbed with a random force. Applied to four independent systems, we successfully show that EpiScore and DCIasym can capture important biophysical properties in dual mutant substitution outcomes. We propose that allosteric regulation and the mechanisms underlying non-additive amino acid substitution outcomes (i.e., epistasis) can be understood as emergent properties of an anisotropic network of interactions where the inclusion of the full network of interactions is critical for accurate modeling. Consequently, mutations which drive towards a new function may require a fine balance between functional site asymmetry and strength of dynamic coupling with the functional sites. These two tools will provide mechanistic insight into both understanding and predicting the outcome of dual mutations. 
    more » « less
  4. ; ; ; ; ; ; ; ( , PNAS Nexus)
    Wand, Josh (Ed.)
    Abstract Hypertrophic cardiomyopathy (HCM) is an inherited disorder often caused by mutations to sarcomeric genes. Many different HCM-associated TPM1 mutations have been identified but they vary in their degrees of severity, prevalence, and rate of disease progression. The pathogenicity of many TPM1 variants detected in the clinical population remains unknown. Our objective was to employ a computational modeling pipeline to assess pathogenicity of one such variant of unknown significance, TPM1 S215L, and validate predictions using experimental methods. Molecular dynamic simulations of tropomyosin on actin suggest that the S215L significantly destabilizes the blocked regulatory state while increasing flexibility of the tropomyosin chain. These changes were quantitatively represented in a Markov model of thin-filament activation to infer the impacts of S215L on myofilament function. Simulations of in vitro motility and isometric twitch force predicted that the mutation would increase Ca2+ sensitivity and twitch force while slowing twitch relaxation. In vitro motility experiments with thin filaments containing TPM1 S215L revealed higher Ca2+ sensitivity compared with wild type. Three-dimensional genetically engineered heart tissues expressing TPM1 S215L exhibited hypercontractility, upregulation of hypertrophic gene markers, and diastolic dysfunction. These data form a mechanistic description of TPM1 S215L pathogenicity that starts with disruption of the mechanical and regulatory properties of tropomyosin, leading thereafter to hypercontractility and finally induction of a hypertrophic phenotype. These simulations and experiments support the classification of S215L as a pathogenic mutation and support the hypothesis that an inability to adequately inhibit actomyosin interactions is the mechanism whereby thin-filament mutations cause HCM. 
    more » « less
  5. ; ; ; ; ( , Communications Biology)
    Abstract Presenilin-1 (PS1) is the catalytic subunit of γ-secretase which cleaves within the transmembrane domain of over 150 peptide substrates. Dominant missense mutations in PS1 cause early-onset familial Alzheimer’s disease (FAD); however, the exact pathogenic mechanism remains unknown. Here we combined Gaussian accelerated molecular dynamics (GaMD) simulations and biochemical experiments to determine the effects of six representative PS1 FAD mutations (P117L, I143T, L166P, G384A, L435F, and L286V) on the enzyme-substrate interactions between γ-secretase and amyloid precursor protein (APP). Biochemical experiments showed that all six PS1 FAD mutations rendered γ-secretase less active for the endoproteolytic (ε) cleavage of APP. Distinct low-energy conformational states were identified from the free energy profiles of wildtype and PS1 FAD-mutant γ-secretase. The P117L and L286V FAD mutants could still sample the “Active” state for substrate cleavage, but with noticeably reduced conformational space compared with the wildtype. The other mutants hardly visited the “Active” state. The PS1 FAD mutants were found to reduce γ-secretase proteolytic activity by hindering APP residue L49 from proper orientation in the active site and/or disrupting the distance between the catalytic aspartates. Therefore, our findings provide mechanistic insights into how PS1 FAD mutations affect structural dynamics and enzyme-substrate interactions of γ-secretase and APP. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email