1
|
Barron MP, Vilseck JZ. A λ-Dynamics Investigation of Insulin Wakayama and Other A3 Variant Binding Affinities to the Insulin Receptor. J Chem Inf Model 2024; 64:5657-5670. [PMID: 38963805 PMCID: PMC11268370 DOI: 10.1021/acs.jcim.4c00662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2024] [Revised: 06/21/2024] [Accepted: 06/24/2024] [Indexed: 07/06/2024]
Abstract
Insulin Wakayama is a clinical insulin variant where a conserved valine at the third residue on insulin's A chain (ValA3) is replaced with a leucine (LeuA3), weakening insulin receptor (IR) binding by 140-500-fold. This severe impact on binding from a subtle modification has posed an intriguing problem for decades. Although experimental investigations of natural and unnatural A3 mutations have highlighted the sensitivity of insulin-IR binding at this site, atomistic explanations of these binding trends have remained elusive. We investigate this problem computationally using λ-dynamics free energy calculations to model structural changes in response to perturbations of the ValA3 side chain and to calculate associated relative changes in binding free energy (ΔΔGbind). The Wakayama LeuA3 mutation and seven other A3 substitutions were studied in this work. The calculated ΔΔGbind results showed high agreement compared to experimental binding potencies with a Pearson correlation of 0.88 and a mean unsigned error of 0.68 kcal/mol. Extensive structural analyses of λ-dynamics trajectories revealed that critical interactions were disrupted between insulin and the insulin receptor as a result of the A3 mutations. This investigation also quantifies the effect that adding an A3 Cδ atom or losing an A3 Cγ atom has on insulin's binding affinity to the IR. Thus, λ-dynamics was able to successfully model the effects of mutations to insulin's A3 side chain on its protein-protein interactions with the IR and shed new light on a decades-old mystery: the exquisite sensitivity of hormone-receptor binding to a subtle modification of an invariant insulin residue.
Collapse
Affiliation(s)
- Monica P Barron
- Department of Biochemistry and Molecular Biology, Indiana University School of Medicine, Indianapolis, Indiana 46202, United States
- Center for Computational Biology and Bioinformatics, Indiana University School of Medicine, Indianapolis, Indiana 46202, United States
| | - Jonah Z Vilseck
- Department of Biochemistry and Molecular Biology, Indiana University School of Medicine, Indianapolis, Indiana 46202, United States
- Center for Computational Biology and Bioinformatics, Indiana University School of Medicine, Indianapolis, Indiana 46202, United States
| |
Collapse
|
2
|
Barron MP, Vilseck JZ. A λ-dynamics investigation of insulin Wakayama and other A3 variant binding affinities to the insulin receptor. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.03.15.585233. [PMID: 38559010 PMCID: PMC10979964 DOI: 10.1101/2024.03.15.585233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Insulin Wakayama is a clinical insulin variant where a conserved valine at the third residue on insulin's A chain (ValA3) is replaced with a leucine (LeuA3), impairing insulin receptor (IR) binding by 140-500 fold. This severe impact on binding from such a subtle modification has posed an intriguing problem for decades. Although experimental investigations of natural and unnatural A3 mutations have highlighted the sensitivity of insulin-IR binding to minor changes at this site, an atomistic explanation of these binding trends has remained elusive. We investigate this problem computationally using λ-dynamics free energy calculations to model structural changes in response to perturbations of the ValA3 side chain and to calculate associated relative changes in binding free energy (ΔΔGbind). The Wakayama LeuA3 mutation and seven other A3 substitutions were studied in this work. The calculated ΔΔGbind results showed high agreement compared to experimental binding potencies with a Pearson correlation of 0.88 and a mean unsigned error of 0.68 kcal/mol. Extensive structural analyses of λ-dynamics trajectories revealed that critical interactions were disrupted between insulin and the insulin receptor as a result of the A3 mutations. This investigation also quantifies the effect that adding an A3 Cδ atom or losing an A3 Cγ atom has on insulin's binding affinity to the IR. Thus, λ-dynamics was able to successfully model the effects of subtle modifications to insulin's A3 side chain on its protein-protein interactions with the IR and shed new light on a decades-old mystery: the exquisite sensitivity of hormone-receptor binding to a subtle modification of an invariant insulin residue.
Collapse
Affiliation(s)
- Monica P. Barron
- Department of Biochemistry and Molecular Biology, Indiana University School of Medicine, Indianapolis, Indiana 46202, United States
- Center for Computational Biology and Bioinformatics, Indiana University School of Medicine, Indianapolis, Indiana 46202, United States
| | - Jonah Z. Vilseck
- Department of Biochemistry and Molecular Biology, Indiana University School of Medicine, Indianapolis, Indiana 46202, United States
- Center for Computational Biology and Bioinformatics, Indiana University School of Medicine, Indianapolis, Indiana 46202, United States
| |
Collapse
|
3
|
Pathophysiological significance of protein hydrophobic interactions: An emerging hypothesis. Med Hypotheses 2017; 110:15-22. [PMID: 29317059 DOI: 10.1016/j.mehy.2017.10.021] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2017] [Revised: 10/09/2017] [Accepted: 10/24/2017] [Indexed: 12/17/2022]
Abstract
Fibrinogen is a unique protein that is converted into an insoluble fibrin in a single enzymatic event, which is a characteristic feature of fibrinogen due to its susceptibility to fibrinolytic degradation and dissolution. Although thrombosis is a result of activated blood coagulation, no explanation is being offered for the persistent presence of fibrin deposits in the affected organs. A classic example is stroke, in which the thrombolytic therapy is effective only during the first 3-4 h after the onset of thrombosis. This phenomenon can now be explained in terms of the modification of fibrinogen structure induced by hydroxyl radicals generated during the period of ischemia caused, in turn, by the blocking of the blood flow within the obstructed vessels. Fibrinogen modification involves intra-to intermolecular disulfide rearrangement induced by the reductive power of hydroxyl radicals that result in the exposition of buried hydrophobic epitopes. Such epitopes react readily with each other forming linkages stronger than the peptide covalent bonds, thus rendering them resistant to the proteolytic degradation. Also, limited reduction of human serum albumin (HSA) generates hydrophobic polymers that form huge insoluble complexes with fibrinogen. Consequently, such insoluble copolymers can be deposited within the circulation of various organs leading to their dysfunction. In conclusion, the study of protein hydrophobic interactions induced by a variety of nutritional and/or environmental factors can provide a rational explanation for a number of pathologic conditions including cardiovascular, neurologic, and other degenerative diseases including cancer.
Collapse
|
4
|
Singh R, Bansal R, Rathore AS, Goel G. Equilibrium Ensembles for Insulin Folding from Bias-Exchange Metadynamics. Biophys J 2017; 112:1571-1585. [PMID: 28445749 DOI: 10.1016/j.bpj.2017.03.015] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2016] [Revised: 03/03/2017] [Accepted: 03/20/2017] [Indexed: 12/29/2022] Open
Abstract
Earliest events in the aggregation process, such as single molecule reconfiguration, are extremely important and the most difficult to characterize in experiments. To this end, we have used well-tempered bias exchange metadynamics simulations to determine the equilibrium ensembles of an insulin molecule under amyloidogenic conditions of low pH and high temperature. A bin-based clustering method that uses statistics accumulated in bias exchange metadynamics trajectories was employed to construct a detailed thermodynamic and kinetic model of insulin folding. The highest lifetime, lowest free-energy ensemble identified consisted of native conformations adopted by a folded insulin monomer in solution, namely, the R-, the Rf-, and the T-states of insulin. The lowest free-energy structure had a root mean square deviation of only 0.15 nm from native x-ray structure. The second longest-lived metastable state was an unfolded, compact monomer with little similarity to the native structure. We have identified three additional long-lived, metastable states from the bin-based model. We then carried out an exhaustive structural characterization of metastable states on the basis of tertiary contact maps and per-residue accessible surface areas. We have also determined the lowest free-energy path between two longest-lived metastable states and confirm earlier findings of non-two-state folding for insulin through a folding intermediate. The ensemble containing the monomeric intermediate retained 58% of native hydrophobic contacts, however, accompanied by a complete loss of native secondary structure. We have discussed the relative importance of nativelike versus nonnative tertiary contacts for the folding transition. We also provide a simple measure to determine the importance of an individual residue for folding transition. Finally, we have compared and contrasted this intermediate with experimental data obtained in spectroscopic, crystallographic, and calorimetric measurements during early stages of insulin aggregation. We have also determined stability of monomeric insulin by incubation at a very low concentration to isolate protein-protein interaction effects.
Collapse
Affiliation(s)
- Richa Singh
- Department of Chemical Engineering, Indian Institute of Technology Delhi, New Delhi, India
| | - Rohit Bansal
- Department of Chemical Engineering, Indian Institute of Technology Delhi, New Delhi, India
| | - Anurag Singh Rathore
- Department of Chemical Engineering, Indian Institute of Technology Delhi, New Delhi, India
| | - Gaurav Goel
- Department of Chemical Engineering, Indian Institute of Technology Delhi, New Delhi, India.
| |
Collapse
|
5
|
Computational study of the activity, dynamics, energetics and conformations of insulin analogues using molecular dynamics simulations: Application to hyperinsulinemia and the critical residue B26. Biochem Biophys Rep 2017; 11:182-190. [PMID: 28955783 PMCID: PMC5614686 DOI: 10.1016/j.bbrep.2017.04.006] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2017] [Accepted: 04/13/2017] [Indexed: 12/15/2022] Open
Abstract
Due to the increasing prevalence of diabetes, finding therapeutic analogues for insulin has become an urgent issue. While many experimental studies have been performed towards this end, they have limited scope to examine all aspects of the effect of a mutation. Computational studies can help to overcome these limitations, however, relatively few studies that focus on insulin analogues have been performed to date. Here, we present a comprehensive computational study of insulin analogues-three mutant insulins that have been identified with hyperinsulinemia and three mutations on the critical B26 residue that exhibit similar binding affinity to the insulin receptor-using molecular dynamics simulations with the aim of predicting how mutations of insulin affect its activity, dynamics, energetics and conformations. The time evolution of the conformers is studied in long simulations. The probability density function and potential of mean force calculations are performed on each insulin analogue to unravel the effect of mutations on the dynamics and energetics of insulin activation. Our conformational study can decrypt the key features and molecular mechanisms that are responsible for an enhanced or reduced activity of an insulin analogue. We find two key results: 1) hyperinsulinemia may be due to the drastically reduced activity (and binding affinity) of the mutant insulins. 2) Y26BS and Y26BE are promising therapeutic candidates for insulin as they are more active than WT-insulin. The analysis in this work can be readily applied to any set of mutations on insulin to guide development of more effective therapeutic analogues.
Collapse
|
6
|
Affiliation(s)
- Pierre De Meyts
- Department of Diabetes Biology; Novo Nordisk A/S; Måløv Denmark
- De Meyts R&D Consulting; Kraainem; Belgium
| |
Collapse
|
7
|
Abstract
Some mutations of the insulin gene cause hyperinsulinemia or hyperproinsulinemia. Replacement of biologically important amino acid leads to defective receptor binding, longer half-life and hyperinsulinemia. Three mutant insulins have been identified: (i) insulin Chicago (F49L or PheB25Leu); (ii) insulin Los Angeles (F48S or PheB24Ser); (iii) and insulin Wakayama (V92L or ValA3Leu). Replacement of amino acid is necessary for proinsulin processing results in hyperproinsulinemia. Four types have been identified: (i) proinsulin Providence (H34D); (ii) proinsulin Tokyo (R89H); (iii) proinsulin Kyoto (R89L); and (iv) proinsulin Oxford (R89P). Three of these are processing site mutations. The mutation of proinsulin Providence, in contrast, is thought to cause sorting abnormality. Compared with normal proinsulin, a significant amount of proinsulin Providence enters the constitutive pathway where processing does not occur. These insulin gene mutations with hyper(pro)insulinemia were very rare, showed only mild diabetes or glucose intolerance, and hyper(pro)insulinemia was the key for their diagnosis. However, this situation changed dramatically after the identification of insulin gene mutations as a cause of neonatal diabetes. This class of insulin gene mutations does not show hyper(pro)insulinemia. Mutations at the cysteine residue or creating a new cysteine will disturb the correct disulfide bonding and proper conformation, and finally will lead to misfolded proinsulin accumulation, endoplasmic reticulum stress and apoptosis of pancreatic β-cells. Maturity-onset diabetes of the young (MODY) or an autoantibody-negative type 1-like phenotype has also been reported. Very recently, recessive mutations with reduced insulin biosynthesis have been reported. The importance of insulin gene mutation in the pathogenesis of diabetes will increase a great deal and give us a new understanding of β-cell biology and diabetes. (J Diabetes Invest, doi: 10.1111/j.2040-1124.2011.00100.x, 2011).
Collapse
Affiliation(s)
- Masahiro Nishi
- Department of Metabolism and Clinical Nutrition, Wakayama Medical University
| | - Kishio Nanjo
- Research Center of Rural Medicine, Nachi‐Katsuura Spa Hospital, Wakayama, Japan
| |
Collapse
|
8
|
Støy J, Steiner DF, Park SY, Ye H, Philipson LH, Bell GI. Clinical and molecular genetics of neonatal diabetes due to mutations in the insulin gene. Rev Endocr Metab Disord 2010; 11:205-15. [PMID: 20938745 PMCID: PMC2974937 DOI: 10.1007/s11154-010-9151-3] [Citation(s) in RCA: 100] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/01/2022]
Abstract
Over the last decade our insight into the causes of neonatal diabetes has greatly expanded. Neonatal diabetes was once considered a variant of type 1 diabetes that presented early in life. Recent advances in our understanding of this disorder have established that neonatal diabetes is not an autoimmune disease, but rather is a monogenic form of diabetes resulting from mutations in a number of different genes encoding proteins that play a key role in the normal function of the pancreatic beta-cell. Moreover, a correct genetic diagnosis can affect treatment and clinical outcome. This is especially true for patients with mutations in the genes KCNJ11 or ABCC8 that encode the two protein subunits (Kir6.2 and SUR1, respectively) of the ATP-sensitive potassium channel. These patients can be treated with oral sulfonylurea drugs with better glycemic control and quality of life. Recently, mutations in the insulin gene (INS) itself have been identified as another cause of neonatal diabetes. In this article, we review the role of INS mutations in the pathophysiology of neonatal diabetes.
Collapse
Affiliation(s)
- Julie Støy
- Department of Internal Medicine and Endocrinology, Aarhus University Hospital, Nørrebrogade 44, 8000, Aarhus C, Denmark.
| | | | | | | | | | | |
Collapse
|
9
|
Yano H, Kitano N, Morimoto M, Polonsky KS, Imura H, Seino Y. A novel point mutation in the human insulin gene giving rise to hyperproinsulinemia (proinsulin Kyoto). J Clin Invest 1992; 89:1902-7. [PMID: 1601997 PMCID: PMC295889 DOI: 10.1172/jci115795] [Citation(s) in RCA: 39] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
We have identified a 65-yr-old nonobese Japanese man with diabetes mellitus, fasting hyperinsulinemia (150-300 pM), and a reduced fasting C-peptide/insulin molar ratio of 2.5-3.0. Fasting hyperinsulinemia was also found in his son and daughter. Analysis of insulin isolated from the serum of the proband and his son by reverse-phase high performance liquid chromatography revealed a minor peak coeluting with human insulin and a major peak of proinsulin-like materials. The insulin gene of the patient was amplified by the polymerase chain reaction and the products were sequenced. A novel point mutation was identified in which guanine was replaced by thymine. The substitution gives rise to a new HindIII recognition site and results in the amino acid replacement of leucine for arginine at position 65. These results indicate that the amino-acid replacement prevents recognition of the C-peptide-A chain dibasic protease and results in an elevation of proinsulin-like materials in the circulation. Furthermore, in this family the proinsulin-like materials is due to a biosynthetic defect, inherited as an autosomal dominant trait. Rapid detection of this mutation can be accomplished by HindIII restriction enzyme mapping of polymerase chain reaction-generated DNA, which enables us to facilitate the diagnosis and screening.
Collapse
Affiliation(s)
- H Yano
- Department of Metabolism and Clinical Nutrition, Kyoto University Faculty of Medicine, Japan
| | | | | | | | | | | |
Collapse
|
10
|
Awata T, Iwamoto Y, Matsuda A, Kuzuya T. Hypervariable region 5'-flanking [Leu A3]insulin gene of insulin Tochigi is different from those of insulin Wakayama I,II. Diabetes Res Clin Pract 1990; 8:183-6. [PMID: 2187661 DOI: 10.1016/0168-8227(90)90115-a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Three families with abnormal insulinemia have been reported in Japan and sequencing analysis revealed that they had the same point mutation in one allele of the insulin genes causing [Leu A3]insulin. To estimate whether or not this same mutation came from a common ancestor we determined the sequence of the hypervariable region 5'-flanking the third [Leu A3]insulin allele (insulin Tochigi). This region is composed of 42 tandem repeating oligonucleotides, is 599 base pairs long and the sequence is 5' cdi jfa faa aba baa aaa fab aaa caa aac aca cba aaf ccb 3' (abbreviated as a = ACAGGGGTGTGGGG; b = ACAGGGGTCTGGGG; c = ACAGGGGTCCTGGGG; d = ACAGGGGTCCGGGG; f = ACAGGGGTCCCGGGG; i = ACAGGGTCCTGGGG; j = ACAGGGGTGTGAGG). The length of this region is different from those of the first and second [Leu A3]insulin alleles (insulin Wakayama I,II). This difference suggests either that insulin Tochigi and insulin Wakayama I,II are not of the same origin, or that three cases of [Leu A3]insulin in Japan have the same ancestor but recombination has occurred in this region at some point in the past.
Collapse
Affiliation(s)
- T Awata
- Division of Endocrinology and Metabolism, Jichi Medical School, Tochigi-ken, Japan
| | | | | | | |
Collapse
|