1
|
Pan S, Hale AT, Lemieux ME, Raval DK, Garton TP, Sadler B, Mahaney KB, Strahle JM. Iron homeostasis and post-hemorrhagic hydrocephalus: a review. Front Neurol 2024; 14:1287559. [PMID: 38283681 PMCID: PMC10811254 DOI: 10.3389/fneur.2023.1287559] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2023] [Accepted: 11/21/2023] [Indexed: 01/30/2024] Open
Abstract
Iron physiology is regulated by a complex interplay of extracellular transport systems, coordinated transcriptional responses, and iron efflux mechanisms. Dysregulation of iron metabolism can result in defects in myelination, neurotransmitter synthesis, and neuronal maturation. In neonates, germinal matrix-intraventricular hemorrhage (GMH-IVH) causes iron overload as a result of blood breakdown in the ventricles and brain parenchyma which can lead to post-hemorrhagic hydrocephalus (PHH). However, the precise mechanisms by which GMH-IVH results in PHH remain elusive. Understanding the molecular determinants of iron homeostasis in the developing brain may lead to improved therapies. This manuscript reviews the various roles iron has in brain development, characterizes our understanding of iron transport in the developing brain, and describes potential mechanisms by which iron overload may cause PHH and brain injury. We also review novel preclinical treatments for IVH that specifically target iron. Understanding iron handling within the brain and central nervous system may provide a basis for preventative, targeted treatments for iron-mediated pathogenesis of GMH-IVH and PHH.
Collapse
Affiliation(s)
- Shelei Pan
- Department of Neurosurgery, Washington University School of Medicine, Washington University in St. Louis, St. Louis, MO, United States
| | - Andrew T. Hale
- Department of Neurosurgery, University of Alabama at Birmingham School of Medicine, University of Alabama at Birmingham, Birmingham, AL, United States
| | - Mackenzie E. Lemieux
- Department of Neurosurgery, Washington University School of Medicine, Washington University in St. Louis, St. Louis, MO, United States
| | - Dhvanii K. Raval
- Department of Neurosurgery, Washington University School of Medicine, Washington University in St. Louis, St. Louis, MO, United States
| | - Thomas P. Garton
- Department of Neurology, Johns Hopkins University School of Medicine, Johns Hopkins University, Baltimore, MD, United States
| | - Brooke Sadler
- Department of Pediatrics, Washington University School of Medicine, Washington University in St. Louis, St. Louis, MO, United States
- Department of Hematology and Oncology, Washington University School of Medicine, Washington University in St. Louis, St. Louis, MO, United States
| | - Kelly B. Mahaney
- Department of Neurosurgery, Stanford University School of Medicine, Stanford University, Palo Alto, CA, United States
| | - Jennifer M. Strahle
- Department of Neurosurgery, Washington University School of Medicine, Washington University in St. Louis, St. Louis, MO, United States
- Department of Pediatrics, Washington University School of Medicine, Washington University in St. Louis, St. Louis, MO, United States
- Department of Orthopedic Surgery, Washington University School of Medicine, Washington University in St. Louis, St. Louis, MO, United States
| |
Collapse
|
2
|
Ma C, Chung DJ, Abramson D, Langley DR, Thayer KM. Mutagenic Activation of Glutathione Peroxidase-4: Approaches toward Rational Design of Allosteric Drugs. ACS OMEGA 2022; 7:29587-29597. [PMID: 36061715 PMCID: PMC9434792 DOI: 10.1021/acsomega.2c01289] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/03/2022] [Accepted: 07/22/2022] [Indexed: 06/15/2023]
Abstract
Glutathione peroxidase 4 (GPX4) reduces lipid hydroperoxides in lipid membranes, effectively inhibiting iron-dependent cell death or ferroptosis. The upregulation of the enzyme by the mutations at residues D21 and D23 has been suggested to be associated with higher protein activity, which confers more protection against neurodegenerative diseases such as Alzheimer's, Parkinson's, and Huntington's diseases. Therefore, it has become an attractive target for treating and preventing neurodegenerative diseases. However, identifying means of mimicking the beneficial effects of these mutations distant from the active site constitutes a formidable challenge in moving toward therapeutics. In this study, we explore using molecular dynamics simulations to computationally map the conformational and energetic landscape of the wild-type GPX4 protein and three mutant variants to identify the allosteric networks of the enzyme. We present the conformational dynamic profile providing the desired signature behavior of the enzyme. We also discuss the implications of these findings for drug design efforts.
Collapse
Affiliation(s)
- Chunyue Ma
- Department
of Mathematics & Computer Science, Wesleyan
University, Middletown, Connecticut 06459, United States
| | - Daniel J. Chung
- Department
of Chemistry, Wesleyan University, Middletown, Connecticut 06459, United States
- Molecular
Biophysics Program, Wesleyan University, Middletown, Connecticut 06459, United States
| | - Dylan Abramson
- Department
of Mathematics & Computer Science, Wesleyan
University, Middletown, Connecticut 06459, United States
| | - David R. Langley
- Department
of Chemistry, Wesleyan University, Middletown, Connecticut 06459, United States
- Molecular
Biophysics Program, Wesleyan University, Middletown, Connecticut 06459, United States
- Arvinas
Inc., New Haven, Connecticut 06511, United States
| | - Kelly M. Thayer
- Department
of Mathematics & Computer Science, Wesleyan
University, Middletown, Connecticut 06459, United States
- Department
of Chemistry, Wesleyan University, Middletown, Connecticut 06459, United States
- Molecular
Biophysics Program, Wesleyan University, Middletown, Connecticut 06459, United States
| |
Collapse
|
3
|
do Bomfim MR, Barbosa DB, de Carvalho PB, da Silva AM, de Oliveira TA, Taranto AG, Leite FHA. Identification of potential human beta-secretase 1 inhibitors by hierarchical virtual screening and molecular dynamics. J Biomol Struct Dyn 2022:1-15. [DOI: 10.1080/07391102.2022.2069155] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Mayra Ramos do Bomfim
- Programa de Pós-Graduação em Ciências Farmacêuticas, Universidade Estadual de Feira de Santana, Feira de Santana, Brazil
| | - Deyse Brito Barbosa
- Programa de Pós-Graduação em Ciências Farmacêuticas, Universidade Estadual de Feira de Santana, Feira de Santana, Brazil
| | | | - Alisson Marques da Silva
- Departamento de Informática, Gestão e Design, Centro Federal de Educação Tecnológica de Minas Gerais, Divinópolis, Brazil
| | - Tiago Alves de Oliveira
- Departamento de Informática, Gestão e Design, Centro Federal de Educação Tecnológica de Minas Gerais, Divinópolis, Brazil
- Departamento de Bioengenharia, Universidade Federal de São João del-Rei, São João del-Rei, Brazil
| | - Alex Gutterres Taranto
- Departamento de Bioengenharia, Universidade Federal de São João del-Rei, São João del-Rei, Brazil
- Faculty of Computing, University of Latvia (UL), Riga, Latvia
| | - Franco Henrique Andrade Leite
- Programa de Pós-Graduação em Ciências Farmacêuticas, Universidade Estadual de Feira de Santana, Feira de Santana, Brazil
| |
Collapse
|
4
|
Gueto-Tettay C, Martinez-Consuegra A, Pelaez-Bedoya L, Drosos-Ramirez JC. G-score: A function to solve the puzzle of modeling the protonation states of β-secretase binding pocket. J Mol Graph Model 2018; 85:1-12. [PMID: 30053756 DOI: 10.1016/j.jmgm.2018.07.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2018] [Revised: 07/13/2018] [Accepted: 07/16/2018] [Indexed: 10/28/2022]
Abstract
The population density concept has emerged as a proposal for the analysis of molecular dynamics results, the key characteristic of population density is the evaluation of the simultaneous occurrence of a set of relevant parameters for a system. However, despite its statistical strength, selection of the tolerance level for the comparison of different models may appear as arbitrary. This work introduces the G-score, a function which summarizes and categorizes the results of population density analysis. Additionally, it incorporates parameters based on rmsd and dihedral angles, besides the protein-protein and protein-ligand interatomic distances conventionally used, which complement each other to provide a better description of the behavior of the system. These newly-proposed tools were applied to determine the most probable protonation state of the aspartic dyad of BACE1, Asp93 and Asp289, in the presence of three types of transition state inhibitors namely: reduced amides, tertiary carbinamines and hydroxyethylamines. The results show a full agreement between G-score values and population density charts, with the advantage of allowing a quick and direct comparison among all the considered models. We anticipate that the simplicity of calculating the parameters employed in this study will permit the extensive use of population density and the G-score for other molecular systems.
Collapse
Affiliation(s)
- Carlos Gueto-Tettay
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Campus San Pablo, 130015, Colombia; Lund University, Faculty of Medicine, Department of Clinical Sciences Lund, Division of Infection Medicine, Lund, Sweden.
| | - Alejandro Martinez-Consuegra
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Campus San Pablo, 130015, Colombia
| | - Luis Pelaez-Bedoya
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Campus San Pablo, 130015, Colombia
| | - Juan Carlos Drosos-Ramirez
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Campus San Pablo, 130015, Colombia.
| |
Collapse
|
5
|
Gueto-Tettay C, Martinez-Consuegra A, Zuchniarz J, Gueto-Tettay LR, Drosos-Ramírez JC. A PM7 dynamic residue-ligand interactions energy landscape of the BACE1 inhibitory pathway by hydroxyethylamine compounds. Part I: The flap closure process. J Mol Graph Model 2017; 76:274-288. [PMID: 28746905 DOI: 10.1016/j.jmgm.2017.07.010] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2016] [Revised: 07/10/2017] [Accepted: 07/11/2017] [Indexed: 01/08/2023]
Abstract
BACE1 is an enzyme of scientific interest because it participates in the progression of Alzheimer's disease. Hydroxyethylamines (HEAs) are a family of compounds which exhibit inhibitory activity toward BACE1 at a nanomolar level, favorable pharmacokinetic properties and oral bioavailability. The first step in the inhibition of BACE1 by HEAs consists of their entrance into the protease active site and the resultant conformational change in the protein, from Apo to closed form. These two conformations differ in the position of an antiparallel loop (called the flap) which covers the entrance to the catalytic site. For BACE1, closure of this flap is vital to its catalytic activity and to inhibition of the enzyme due to the new interactions thereby formed with the ligand. In the present study a dynamic energy landscape of residue-ligand interaction energies (ReLIE) measured for 112 amino acids in the BACE1 active site and its immediate vicinity during the closure of the flap induced by 8 HEAs of different inhibitory power is presented. A total of 6.272 million ReLIE calculations, based on the PM7 semiempirical method, provided a deep and quantitative view of the first step in the inhibition of the aspartyl protease. The information suggests that residues Asp93, Asp289, Thr292, Thr293, Asn294 and Arg296 are anchor points for the ligand, accounting for approximately 45% of the total protein-ligand interaction. Additionally, flap closure improved the BACE1-HEA interaction by around 25%. Furthermore, the inhibitory activity of HEAs could be related to the capacity of these ligands to form said anchor point interactions and maintain them over time: the lack of some of these anchor interactions delayed flap closure or impeded it completely, or even caused the flap to reopen. The methodology employed here could be used as a tool to evaluate future structural modifications which lead to improvements in the favorability and stability of BACE1-HEA ReLIEs, aiding in the design of better inhibitors.
Collapse
Affiliation(s)
- Carlos Gueto-Tettay
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Cartagena, Colombia.
| | - Alejandro Martinez-Consuegra
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Cartagena, Colombia
| | - Joshua Zuchniarz
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Cartagena, Colombia
| | - Luis Roberto Gueto-Tettay
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Cartagena, Colombia
| | - Juan Carlos Drosos-Ramírez
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Cartagena, Colombia.
| |
Collapse
|
6
|
Hernández-Rodríguez M, Correa-Basurto J, Gutiérrez A, Vitorica J, Rosales-Hernández MC. Asp32 and Asp228 determine the selective inhibition of BACE1 as shown by docking and molecular dynamics simulations. Eur J Med Chem 2016; 124:1142-1154. [DOI: 10.1016/j.ejmech.2016.08.028] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2015] [Revised: 08/10/2016] [Accepted: 08/13/2016] [Indexed: 11/28/2022]
|
7
|
Gueto-Tettay C, Zuchniarz J, Fortich-Seca Y, Gueto-Tettay LR, Drosos-Ramirez JC. A molecular dynamics study of the BACE1 conformational change from Apo to closed form induced by hydroxyethylamine derived compounds. J Mol Graph Model 2016; 70:181-195. [PMID: 27750187 DOI: 10.1016/j.jmgm.2016.10.006] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2016] [Revised: 10/05/2016] [Accepted: 10/06/2016] [Indexed: 12/11/2022]
Abstract
BACE1 is an aspartyl protease which is a therapeutic target for Alzheimer's disease (AD) because of its participation in the rate-limiting step in the production of Aβ-peptide, the accumulation of which produces senile plaques and, in turn, the neurodegenerative effects associated with AD. The active site of this protease is composed in part by two aspartic residues (Asp93 and Asp289). Additionally, the catalytic site has been found to be covered by an antiparallel hairpin loop called the flap. The dynamics of this flap are fundamental to the catalytic function of the enzyme. When BACE1 is inactive (Apo), the flap adopts an open conformation, allowing a substrate or inhibitor to access the active site. Subsequent interaction with the ligand induces flap closure and the stabilization of the macromolecular complex. Further, the protonation state of the aspartic dyad is affected by the chemical nature of the species entering the active site, so that appropriate selection of protonation states for the ligand and the catalytic residues will permit the elucidation of the inhibitory pathway for BACE1. In the present study, comparative analysis of different combinations of protonation states for the BACE1-hydroxyethylamine (HEA) system is reported. HEAs are potent inhibitors of BACE1 with favorable pharmacological and kinetic properties, as well as oral bioavailability. The results of Molecular Dynamics (MD) simulations and population density calculations using 8 different parameters demonstrate that the LnAsp289 configuration (HEA with a neutral amine and the Asp289 residue protonated) is the only one which permits the expected conformational change in BACE1, from apo to closed form, after flap closure. Additionally, differences in their capacities to establish and maintain interactions with residues such as Asp93, Gly95, Thr133, Asp289, Gly291, and Asn294 during this step allow differentiation among the inhibitory activities of the HEAs. The results and methodology here reported will serve to elucidate the inhibitory pathway of other families of compounds that act as BACE1 inhibitors, as well as the design of better leader compounds for the treatment of AD.
Collapse
Affiliation(s)
- Carlos Gueto-Tettay
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Colombia
| | - Joshua Zuchniarz
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Colombia
| | - Yeyson Fortich-Seca
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Colombia
| | - Luis Roberto Gueto-Tettay
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Colombia
| | - Juan Carlos Drosos-Ramirez
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Colombia.
| |
Collapse
|
8
|
Gueto-Tettay C, Pestana-Nobles R, Drosos-Ramirez JC. Determination of the protonation state for the catalytic dyad in β-secretase when bound to hydroxyethylamine transition state analogue inhibitors: A molecular dynamics simulation study. J Mol Graph Model 2016; 66:155-67. [PMID: 27111489 DOI: 10.1016/j.jmgm.2016.04.003] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2016] [Revised: 04/09/2016] [Accepted: 04/12/2016] [Indexed: 01/12/2023]
Abstract
BACE1 is an aspartyl protease of pharmacological interest for its direct participation in Alzheimer's disease (AD) through β-amyloid peptide production. Two aspartic acid residues are present in the BACE1 catalytic region which can adopt multiple protonation states depending on the chemical nature of its inhibitors, i.e., monoprotonated, diprotonated and di-deprotonated states. In the present study a series of protein-ligand molecular dynamics (MD) simulations was carried out to identify the most feasible protonation state adopted by the catalytic dyad in the presence of hydroxyethylamine transition state analogue inhibitors. The MD trajectories revealed that the di-deprotonated state is most prefered in the presence of hydroxyethilamine (HEA) family inhibitors. This appears as a result after evaluating, for all 9 protonation state configurations during the simulation time, the deviations of a set of distances and dihedral angles measured on the ligand, protein and protein-ligand complex with reference to an X-ray experimental BACE1/HEA crystallographic structure. These results will help to clarify the phenomena related to the HEAs inhibitory pathway, and improve HEAs databases' virtual screening and ligand design processes targeting β-secretase protein.
Collapse
Affiliation(s)
- Carlos Gueto-Tettay
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Colombia
| | - Roberto Pestana-Nobles
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Colombia
| | - Juan Carlos Drosos-Ramirez
- Grupo de Química Bioorgánica, Facultad de Ciencias Exactas y Naturales, Universidad de Cartagena, Colombia.
| |
Collapse
|
9
|
Communication routes in ARID domains between distal residues in helix 5 and the DNA-binding loops. PLoS Comput Biol 2014; 10:e1003744. [PMID: 25187961 PMCID: PMC4154638 DOI: 10.1371/journal.pcbi.1003744] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2014] [Accepted: 06/12/2014] [Indexed: 11/19/2022] Open
Abstract
ARID is a DNA-binding domain involved in several transcriptional regulatory processes, including cell-cycle regulation and embryonic development. ARID domains are also targets of the Human Cancer Protein Interaction Network. Little is known about the molecular mechanisms related to conformational changes in the family of ARID domains. Thus, we have examined their structural dynamics to enrich the knowledge on this important family of regulatory proteins. In particular, we used an approach that integrates atomistic simulations and methods inspired by graph theory. To relate these properties to protein function we studied both the free and DNA-bound forms. The interaction with DNA not only stabilizes the conformations of the DNA-binding loops, but also strengthens pre-existing paths in the native ARID ensemble for long-range communication to those loops. Residues in helix 5 are identified as critical mediators for intramolecular communication to the DNA-binding regions. In particular, we identified a distal tyrosine that plays a key role in long-range communication to the DNA-binding loops and that is experimentally known to impair DNA-binding. Mutations at this tyrosine and in other residues of helix 5 are also demonstrated, by our approach, to affect the paths of communication to the DNA-binding loops and alter their native dynamics. Overall, our results are in agreement with a scenario in which ARID domains exist as an ensemble of substates, which are shifted by external perturbation, such as the interaction with DNA. Conformational changes at the DNA-binding loops are transmitted long-range by intramolecular paths, which have their heart in helix 5.
Collapse
|
10
|
Molecular Dynamics Simulations of Bromodomains Reveal Binding-Site Flexibility and Multiple Binding Modes of the Natural Ligand Acetyl-Lysine. Isr J Chem 2014. [DOI: 10.1002/ijch.201400009] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
|
11
|
Jia X, Zeng J, Zhang JZH, Mei Y. Accessing the applicability of polarized protein-specific charge in linear interaction energy analysis. J Comput Chem 2014; 35:737-47. [PMID: 24500844 DOI: 10.1002/jcc.23547] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2013] [Revised: 11/15/2013] [Accepted: 01/05/2014] [Indexed: 12/12/2022]
Abstract
The reliability of the linear interaction energy (LIE) depends on the atomic charge model used to delineate the Coulomb interaction between the ligand and its environment. In this work, the polarized protein-specific charge (PPC) implementing a recently proposed fitting scheme has been examined in the LIE calculations of the binding affinities for avidin and β-secretase binding complexes. This charge fitting scheme, termed delta restrained electrostatic potential, bypasses the prevalent numerical difficulty of rank deficiency in electrostatic-potential-based charge fitting methods via a dual-step fitting strategy. A remarkable consistency between the predicted binding affinities and the experimental measurement has been observed. This work serves as a direct evidence of PPC's applicability in rational drug design.
Collapse
Affiliation(s)
- Xiangyu Jia
- State Key Laboratory of Precision Spectroscopy, Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, 200062, China
| | | | | | | |
Collapse
|
12
|
Dhanjal JK, Goyal S, Sharma S, Hamid R, Grover A. Mechanistic insights into mode of action of potent natural antagonists of BACE-1 for checking Alzheimer's plaque pathology. Biochem Biophys Res Commun 2013; 443:1054-9. [PMID: 24365147 DOI: 10.1016/j.bbrc.2013.12.088] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2013] [Accepted: 12/17/2013] [Indexed: 11/16/2022]
Abstract
Alzheimer's is a neurodegenerative disorder resulting in memory loss and decline in cognitive abilities. Accumulation of extracellular beta amyloidal plaques is one of the major pathology associated with this disease. β-Secretase or BACE-1 performs the initial and rate limiting step of amyloidic pathway in which 37-43 amino acid long peptides are generated which aggregate to form plaques. Inhibition of this enzyme offers a viable prospect to check the growth of these plaques. Numerous efforts have been made in recent years for the generation of BACE-1 inhibitors but many of them failed during the preclinical or clinical trials due to drug related or drug induced toxicity. In the present work, we have used computational methods to screen a large dataset of natural compounds to search for small molecules having BACE-1 inhibitory activity with low toxicity to normal cells. Molecular dynamics simulations were performed to analyze molecular interactions between the screened compounds and the active residues of the enzyme. Herein, we report two natural compounds of inhibitory nature active against β-secretase enzyme of amyloidic pathway and are potent lead molecules against Alzheimer's disease.
Collapse
Affiliation(s)
| | - Sukriti Goyal
- Apaji Institute of Mathematics & Applied Computer Technology, Banasthali University, Tonk 304022, Rajasthan, India
| | - Sudhanshu Sharma
- Department of Biotechnology, Delhi Technological University, New Delhi 110042, India
| | - Rabia Hamid
- Department of Biochemistry, University of Kashmir, Srinagar 190006, India
| | - Abhinav Grover
- School of Biotechnology, Jawaharlal Nehru University, New Delhi 110067, India.
| |
Collapse
|
13
|
Buonfiglio R, Ferraro M, Falchi F, Cavalli A, Masetti M, Recanatini M. Collecting and assessing human lactate dehydrogenase-A conformations for structure-based virtual screening. J Chem Inf Model 2013; 53:2792-7. [PMID: 24138094 DOI: 10.1021/ci400543y] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Human lactate dehydrogenase-A (LDHA) is emerging as a promising anticancer target. Up to now, structure-based investigations for identifying inhibitors of this enzyme have not explicitly accounted for active site flexibility. In the present study, by combining replica exchange molecular dynamics with network and cluster analyses, we identified reliable LDHA conformations for structure-based ligand design. The selected conformations were challenged and validated by retrospective virtual screening simulations.
Collapse
Affiliation(s)
- Rosa Buonfiglio
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum-Università di Bologna , via Belmeloro 6, 40126 Bologna, Italy
| | | | | | | | | | | |
Collapse
|
14
|
Ahlstrom LS, Baker JL, Ehrlich K, Campbell ZT, Patel S, Vorontsov II, Tama F, Miyashita O. Network visualization of conformational sampling during molecular dynamics simulation. J Mol Graph Model 2013; 46:140-9. [PMID: 24211466 DOI: 10.1016/j.jmgm.2013.10.003] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2013] [Revised: 09/06/2013] [Accepted: 10/03/2013] [Indexed: 02/01/2023]
Abstract
Effective data reduction methods are necessary for uncovering the inherent conformational relationships present in large molecular dynamics (MD) trajectories. Clustering algorithms provide a means to interpret the conformational sampling of molecules during simulation by grouping trajectory snapshots into a few subgroups, or clusters, but the relationships between the individual clusters may not be readily understood. Here we show that network analysis can be used to visualize the dominant conformational states explored during simulation as well as the connectivity between them, providing a more coherent description of conformational space than traditional clustering techniques alone. We compare the results of network visualization against 11 clustering algorithms and principal component conformer plots. Several MD simulations of proteins undergoing different conformational changes demonstrate the effectiveness of networks in reaching functional conclusions.
Collapse
Affiliation(s)
- Logan S Ahlstrom
- Department of Chemistry and Biochemistry, University of Arizona, Tucson, AZ 85721, USA
| | | | | | | | | | | | | | | |
Collapse
|
15
|
LI QIANG, CHEN MIN, LIU HONGMIN, YANG LIQUN, YANG GUIYING. Expression of APP, BACE1, AChE and ChAT in an AD model in rats and the effect of donepezil hydrochloride treatment. Mol Med Rep 2012; 6:1450-4. [DOI: 10.3892/mmr.2012.1102] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2012] [Accepted: 09/14/2012] [Indexed: 11/06/2022] Open
|
16
|
Roy A, Post CB. Detection of long-range concerted motions in protein by a distance covariance. J Chem Theory Comput 2012; 8:3009-3014. [PMID: 23610564 PMCID: PMC3630994 DOI: 10.1021/ct300565f] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We asses the ability of a distance correlation coefficient (DiCC), calculated from distance covariance, for detecting long-range concerted motion in proteins. We establish a set of criteria for ideal correlation coefficient values based on the coefficient of determination in multi-dimension, R2. We compare in detail DiCC and conventional coefficients against these criteria. We demonstrate that in contrast to conventional correlation coefficients, which capture long-distance correlation adequately only with certain restrictions in multi-dimension, DiCC reflects appropriate correlation in both one- and multi-dimension. Finally we demonstrate the usefulness of DiCC for assessing long-distance correlated fluctuation in protein dynamics.
Collapse
Affiliation(s)
- Amitava Roy
- Medicinal Chemistry and Molecular Pharmacology, Purdue University, West Lafayette, USA
| | - Carol Beth Post
- Medicinal Chemistry and Molecular Pharmacology, Purdue University, West Lafayette, USA
| |
Collapse
|
17
|
Liu S, Fu R, Cheng X, Chen SP, Zhou LH. Exploring the binding of BACE-1 inhibitors using comparative binding energy analysis (COMBINE). BMC STRUCTURAL BIOLOGY 2012; 12:21. [PMID: 22925713 PMCID: PMC3533579 DOI: 10.1186/1472-6807-12-21] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/24/2011] [Accepted: 08/03/2012] [Indexed: 01/14/2023]
Abstract
BACKGROUND The inhibition of the activity of β-secretase (BACE-1) is a potentially important approach for the treatment of Alzheimer disease. To explore the mechanism of inhibition, we describe the use of 46 X-ray crystallographic BACE-1/inhibitor complexes to derive quantitative structure-activity relationship (QSAR) models. The inhibitors were aligned by superimposing 46 X-ray crystallographic BACE-1/inhibitor complexes, and gCOMBINE software was used to perform COMparative BINding Energy (COMBINE) analysis on these 46 minimized BACE-1/inhibitor complexes. The major advantage of the COMBINE analysis is that it can quantitatively extract key residues involved in binding the ligand and identify the nature of the interactions between the ligand and receptor. RESULTS By considering the contributions of the protein residues to the electrostatic and van der Waals intermolecular interaction energies, two predictive and robust COMBINE models were developed: (i) the 3-PC distance-dependent dielectric constant model (built from a single X-ray crystal structure) with a q2 value of 0.74 and an SDEC value of 0.521; and (ii) the 5-PC sigmoidal electrostatic model (built from the actual complexes present in the Brookhaven Protein Data Bank) with a q2 value of 0.79 and an SDEC value of 0.41. CONCLUSIONS These QSAR models and the information describing the inhibition provide useful insights into the design of novel inhibitors via the optimization of the interactions between ligands and those key residues of BACE-1.
Collapse
Affiliation(s)
- Shu Liu
- Guangdong Province Key Laboratory of Functional Molecules in Oceanic Microorganism, Zhong Shan School of Medicine, Sun Yat-Sen University, Guangzhou, 510080, People’s Republic of China
- Department of Anatomy, Zhong Shan School of Medicine, Sun Yat-Sen University, Guangzhou, 510080, People’s Republic of China
| | - Rao Fu
- Department of Anatomy, Zhong Shan School of Medicine, Sun Yat-Sen University, Guangzhou, 510080, People’s Republic of China
| | - Xiao Cheng
- Department of Anatomy, Zhong Shan School of Medicine, Sun Yat-Sen University, Guangzhou, 510080, People’s Republic of China
| | - Sheng-Ping Chen
- Guangdong Province Key Laboratory of Functional Molecules in Oceanic Microorganism, Zhong Shan School of Medicine, Sun Yat-Sen University, Guangzhou, 510080, People’s Republic of China
| | - Li-Hua Zhou
- Department of Anatomy, Zhong Shan School of Medicine, Sun Yat-Sen University, Guangzhou, 510080, People’s Republic of China
| |
Collapse
|
18
|
Zhou T, Caflisch A. Distribution of Reciprocal of Interatomic Distances: A Fast Structural Metric. J Chem Theory Comput 2012; 8:2930-7. [PMID: 26592131 DOI: 10.1021/ct3003145] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
We present a structural metric based on the Distribution of Reciprocal of Interatomic Distances (DRID) for evaluating geometrical similarity between two conformations of a molecule. A molecular conformation is described by a vector of 3N orientation-independent DRID descriptors where N is the number of molecular centroids, for example, the non-hydrogen atoms in all nonsymmetric groups of a peptide. For two real-world applications (pairwise comparison of snapshots from an explicit solvent simulation of a protease/peptide substrate complex and implicit solvent simulations of reversible folding of a 20-residue β-sheet peptide), the DRID-based metric is shown to be about 5 and 15 times faster than coordinate root-mean-square deviation (cRMSD) and distance root-mean-square deviation (dRMSD), respectively. A single core of a mainstream processor can perform about 10(8) DRID comparisons in one-half a minute. Importantly, the DRID metric has closer similarity to kinetic distance than does either cRMSD or dRMSD. Therefore, DRID is suitable for clustering molecular dynamics trajectories for kinetic analysis, for example, by Markov state models. Moreover, conformational space networks and free energy profiles derived by DRID-based clustering preserve the kinetic information.
Collapse
Affiliation(s)
- Ting Zhou
- Department of Biochemistry, University of Zurich , CH-8057 Zurich, Switzerland
| | - Amedeo Caflisch
- Department of Biochemistry, University of Zurich , CH-8057 Zurich, Switzerland
| |
Collapse
|
19
|
Brás NF, Ramos MJ, Fernandes PA. The catalytic mechanism of mouse renin studied with QM/MM calculations. Phys Chem Chem Phys 2012; 14:12605-13. [PMID: 22796659 DOI: 10.1039/c2cp41422h] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Hypertension is a chronic condition that affects nearly 25% of adults worldwide. As the Renin-Angiotensin-Aldosterone System is implicated in the control of blood pressure and body fluid homeostasis, its combined blockage is an attractive therapeutic strategy currently in use for the treatment of several cardiovascular conditions. We have performed QM/MM calculations to study the mouse renin catalytic mechanism in atomistic detail, using the N-terminal His6-Asn14 segment of angiotensinogen as substrate. The enzymatic reaction (hydrolysis of the peptidic bond between residues in the 10th and 11th positions) occurs through a general acid/base mechanism and, surprisingly, it is characterized by three mechanistic steps: it begins with the creation of a first very stable tetrahedral gem-diol intermediate, followed by protonation of the peptidic bond nitrogen, giving rise to a second intermediate. In a final step the peptidic bond is completely cleaved and both gem-diol hydroxyl protons are transferred to the catalytic dyad (Asp32 and Asp215). The final reaction products are two separate peptides with carboxylic acid and amine extremities. The activation energy for the formation of the gem-diol intermediate was calculated as 23.68 kcal mol(-1), whereas for the other steps the values were 15.51 kcal mol(-1) and 14.40 kcal mol(-1), respectively. The rate limiting states were the reactants and the first transition state. The associated barrier (23.68 kcal mol(-1)) is close to the experimental values for the angiotensinogen substrate (19.6 kcal mol(-1)). We have also tested the influence of the density functional on the activation and reaction energies. All eight density functionals tested (B3LYP, B3LYP-D3, X3LYP, M06, B1B95, BMK, mPWB1K and B2PLYP) gave very similar results.
Collapse
Affiliation(s)
- Natércia F Brás
- REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre s/n, 4169-007 Porto, Portugal
| | | | | |
Collapse
|