1
|
Barraza GA, Castro-Guijarro AC, de la Fuente Hoffmann V, Bolívar Ávila SJ, Flamini MI, Sanchez AM. Drug repositioning for rosacea disease: Biological TARGET identification, molecular docking, pharmacophore mapping, and molecular dynamics analysis. Comput Biol Med 2024; 181:108988. [PMID: 39168013 DOI: 10.1016/j.compbiomed.2024.108988] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2024] [Revised: 07/16/2024] [Accepted: 08/02/2024] [Indexed: 08/23/2024]
Abstract
Rosacea is a chronic dermatological condition that currently lacks a clear treatment approach due to an uncomprehensive knowledge of its pathogenesis. The main obstacle lies in understanding its etiology and the mode of action of the different drugs used. This study aims to clarify these aspects by employing drug repositioning. Using an in silico approach, we performed a transcriptomic analysis comparing samples from individuals with diverse types of rosacea to those from healthy controls to identify genes deregulated in this disease. Subsequently, we realized molecular docking and molecular dynamics studies to assess the binding affinity of drugs currently used to treat rosacea and drugs that target proteins interacting with, and thus affecting, proteins deregulated in rosacea. Our findings revealed that the downregulation of SKAP2 and upregulation of S100A7A in rosacea, could be involved in the pathogenesis of the disease. Furthermore, considering the drugs currently used for rosacea management, we demonstrated stable interactions between isotretinoin and BFH772 with SKAP2, and permethrin and PAC-14028 with S100A7A. Similarly, considering drugs targeting SKAP2 and S100A7A interactome proteins, we found that pitavastatin and dasatinib exert stable interactions with SKAP2, and lovastatin and tirbanibulin with S100A7A. In addition, we determine that the types of bonds involved in the interactions were different in SKAP2 from S100A7A. The drug-SKAP2 interactions are hydrogen bonds, whereas the drug-S100A7A interactions are of the hydrophobic type. In conclusion, our study provides evidence for the possible contribution of SKAP2 and S100A7A to rosacea pathology. Furthermore, it provides significant information on the molecular interactions between drugs and these proteins, highlighting the importance of considering structural features and binding interactions in the design of targeted therapies for skin disorders such as rosacea.
Collapse
Affiliation(s)
- Gustavo Adolfo Barraza
- Laboratorio de Transducción de Señales y Movimiento Celular, Instituto de Medicina y Biología Experimental de Cuyo (IMBECU), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo, Mendoza, Argentina
| | - Ana Carla Castro-Guijarro
- Laboratorio de Biología Tumoral, Instituto de Medicina y Biología Experimental de Cuyo (IMBECU), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo, Mendoza, Argentina
| | - Valentina de la Fuente Hoffmann
- Laboratorio de Transducción de Señales y Movimiento Celular, Instituto de Medicina y Biología Experimental de Cuyo (IMBECU), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo, Mendoza, Argentina
| | - Santiago Junior Bolívar Ávila
- Instituto de Química Rosario (IQUIR), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), and Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, Santa Fe, Argentina
| | - Marina Inés Flamini
- Laboratorio de Biología Tumoral, Instituto de Medicina y Biología Experimental de Cuyo (IMBECU), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo, Mendoza, Argentina.
| | - Angel Matias Sanchez
- Laboratorio de Transducción de Señales y Movimiento Celular, Instituto de Medicina y Biología Experimental de Cuyo (IMBECU), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo, Mendoza, Argentina.
| |
Collapse
|
2
|
Dziadek ŁJ, Sieradzan AK, Czaplewski C, Zalewski M, Banaś F, Toczek M, Nisterenko W, Grudinin S, Liwo A, Giełdoń A. Assessment of Four Theoretical Approaches to Predict Protein Flexibility in the Crystal Phase and Solution. J Chem Theory Comput 2024; 20:7667-7681. [PMID: 39171852 DOI: 10.1021/acs.jctc.4c00754] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
In this paper, we evaluated the ability of four coarse-grained methods to predict protein flexible regions with potential biological importance, UNRES-flex, UNRES-DSSP-flex (based on the united residue model of polypeptide chains without and with secondary structure restraints, respectively), CABS-flex (based on the C-α, C-β, and side chain model), and nonlinear rigid block normal mode analysis (NOLB) with a set of 100 protein structures determined by NMR spectroscopy or X-ray crystallography, with all secondary structure types. End regions with high fluctuations were excluded from analysis. The Pearson and Spearman correlation coefficients were used to quantify the conformity between the calculated and experimental fluctuation profiles, the latter determined from NMR ensembles and X-ray B-factors, respectively. For X-ray structures (corresponding to proteins in a crowded environment), NOLB resulted in the best agreement between the predicted and experimental fluctuation profiles, while for NMR structures (corresponding to proteins in solution), the ranking of performance is CABS-flex > UNRES-DSSP-flex > UNRES-flex > NOLB; however, CABS-flex sometimes exaggerated the extent of small fluctuations, as opposed to UNRES-DSSP-flex.
Collapse
Affiliation(s)
- Ł J Dziadek
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - A K Sieradzan
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - C Czaplewski
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
- School of Computational Sciences, Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-gu, Seoul 02455, Republic of Korea
| | - M Zalewski
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - F Banaś
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - M Toczek
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - W Nisterenko
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - S Grudinin
- LJK, University Grenoble Alpes, CNRS, Grenoble INP, F-38000 Grenoble, France
| | - A Liwo
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| | - A Giełdoń
- Faculty of Chemistry, University of Gdansk, ul. Wita-Stwosza 63, 80-308 Gdańsk, Poland
| |
Collapse
|
3
|
Abbass J, Parisi C. Machine learning-based prediction of proteins' architecture using sequences of amino acids and structural alphabets. J Biomol Struct Dyn 2024:1-16. [PMID: 38505995 DOI: 10.1080/07391102.2024.2328736] [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: 11/28/2023] [Accepted: 03/05/2024] [Indexed: 03/21/2024]
Abstract
In addition to the growth of protein structures generated through wet laboratory experiments and deposited in the PDB repository, AlphaFold predictions have significantly contributed to the creation of a much larger database of protein structures. Annotating such a vast number of structures has become an increasingly challenging task. CATH is widely recognized as one the most common platforms for addressing this challenge, as it classifies proteins based on their structural and evolutionary relationships, offering the scientific community an invaluable resource for uncovering various properties, including functional annotations. While CATH annotation involves - to some extent - human intervention, keeping up with the classification of the rapidly expanding repositories of protein structures has become exceedingly difficult. Therefore, there is a pressing need for a fully automated approach. On the other hand, the abundance of protein sequences stemming from next generation sequencing technologies, lacking structural annotations, presents an additional challenge to the scientific community. Consequently, 'pre-annotating' protein sequences with structural features, ensuring a high level of precision, could prove highly advantageous. In this paper, after a thorough investigation, we introduce a novel machine-learning model capable of classifying any protein domain, whether it has a known structure or not, into one of the 40 main CATH Architectures. We achieve an F1 Score of 0.92 using only the amino acid sequence and a score of 0.94 using both the sequence of amino acids and the sequence of structural alphabets.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Jad Abbass
- School of Computer Science and Mathematics, Kingston University, London, UK
| | - Charles Parisi
- School of Computer Science and Mathematics, Kingston University, London, UK
- Telecom Physique Strasbourg, Strasbourg University, Strasbourg, France
| |
Collapse
|
4
|
General Trends of the Camelidae Antibody V HHs Domain Dynamics. Int J Mol Sci 2023; 24:ijms24054511. [PMID: 36901942 PMCID: PMC10003728 DOI: 10.3390/ijms24054511] [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: 02/09/2023] [Revised: 02/22/2023] [Accepted: 02/23/2023] [Indexed: 03/03/2023] Open
Abstract
Conformational flexibility plays an essential role in antibodies' functional and structural stability. They facilitate and determine the strength of antigen-antibody interactions. Camelidae express an interesting subtype of single-chain antibody, named Heavy Chain only Antibody. They have only one N-terminal Variable domain (VHH) per chain, composed of Frameworks (FRs) and Complementarity Determining regions (CDRs) like their VH and VL counterparts in IgG. Even when expressed independently, VHH domains display excellent solubility and (thermo)stability, which helps them to retain their impressive interaction capabilities. Sequence and structural features of VHH domains contributing to these abilities have already been studied compared to classical antibodies. To have the broadest view and understand the changes in dynamics of these macromolecules, large-scale molecular dynamics simulations for a large number of non-redundant VHH structures have been performed for the first time. This analysis reveals the most prevalent movements in these domains. It reveals the four main classes of VHHs dynamics. Diverse local changes were observed in CDRs with various intensities. Similarly, different types of constraints were observed in CDRs, while FRs close to CDRs were sometimes primarily impacted. This study sheds light on the changes in flexibility in different regions of VHH that may impact their in silico design.
Collapse
|
5
|
Alanazi AM, Bakheit AH, Attwa MW, Abdelhameed AS. Spectroscopic, molecular docking and dynamic simulation studies of binding between the new anticancer agent olmutinib and human serum albumin. J Biomol Struct Dyn 2022; 40:14236-14246. [PMID: 34766879 DOI: 10.1080/07391102.2021.2001380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
In this study, the interaction between human serum albumin (HSA), which is the key bio-distributor of exogenous and endogenous compounds in the human bloodstream, and HM61713 (Olmutinib; OMB), which is used as an anticancer drug, is examined by multiple spectroscopic techniques (steady-state fluorescence, UV spectrophotometry, synchronous, and 3 D fluorescence) combined with molecular docking and molecular dynamic simulation investigations. The fluorescence results clearly demonstrated quenching in HSA fluorescence in the existence of OMB indicating the formation of complex and have also shown that the interaction is static. Fluorescence spectroscopy was used to obtain the binding constant values that revealed a strong interaction between the HSA and OMB at 298 K with a binding constant of 7.39x104 M-1 suggesting strong interaction. OMB binds to HSA at site I (IIA). Electrostatic forces and H-bonding were the main binding forces of main bonding between HSA and OMB as revealed by docking and thermodynamic results.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Amer M Alanazi
- Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Ahmed H Bakheit
- Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Mohamed W Attwa
- Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Ali Saber Abdelhameed
- Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| |
Collapse
|
6
|
Reis CRCD, Souza HVCMD, Leme RPP, Castelo-Branco FS, Fernandes TVA, Boechat N, Dias LRS, Hoelz LVB. Study of the dynamic behavior of the cruzain enzyme in free and complexed forms with competitive and noncovalent benzimidazole inhibitors. J Biomol Struct Dyn 2022:1-15. [PMID: 35475501 DOI: 10.1080/07391102.2022.2067238] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
There are only two drugs for the treatment of Chagas disease, namely, nifurtimox and benznidazole, that can cause several adverse effects. Despite the effectiveness of these drugs in the disease's acute phase, they are not recognized as curative in the chronic phase, establishing the need for more effective treatment in all stages of the disease. Cruzain is an enzyme that plays a vital role in the life cycle of the etiologic agent, the protozoan Trypanosoma cruzi, being relevant as a therapeutic target in the planning of new drugs. Using molecular docking and dynamics simulations, we have investigated the structural and dynamic factors that can be involved in the enzyme inhibition process at the atomic-molecular level by benzimidazole compounds that are potent cruzain inhibitors with in vitro trypanocidal activity. The study suggests that these inhibitors bind cruzain through steric and hydrogen bonding interactions without altering its secondary structure content and protein compaction. Besides, we observed that these inhibitors decrease the correlation of movements between Cα-atoms of cruzain, increasing the number of atomic communities, mainly in the α-helix that presents the catalytic Cys25 residue. As expected, we also observed a correlation between the inhibitory activity of each inhibitor and their respective binding-free energies, reinforcing that the affinity of the complexes seems to be a relevant factor for enzymatic inhibition. Hence, the results presented in this work contribute to a better understanding of the cruzain enzyme inhibition mechanism through competitive and non-covalent inhibitors.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Caroline Rodrigues Chaves Dos Reis
- Programa de Pós-Graduação em Ciências Aplicadas a Produtos Para Saúde (PPG-CAPS), Faculdade de Farmácia, Universidade Federal Fluminense, Rua Mário Viana, Niterói, RJ, Brazil.,Faculdade de Farmácia, Laboratório de Química Medicinal - LQMed, Universidade Federal Fluminense (UFF), Santa Rosa, Niterói, RJ, Brazil.,Fundação Oswaldo Cruz, Instituto de Tecnologia em Fármacos, Farmanguinhos - Fiocruz, Laboratório de Síntese Fármacos - LASFAR, Manguinhos Rio de Janeiro, RJ, Brazil
| | - Hellen Valério Chaves Moura de Souza
- Fundação Oswaldo Cruz, Instituto de Tecnologia em Fármacos, Farmanguinhos - Fiocruz, Laboratório de Síntese Fármacos - LASFAR, Manguinhos Rio de Janeiro, RJ, Brazil
| | - Rennan Papaleo Paes Leme
- Fundação Oswaldo Cruz, Instituto de Tecnologia em Fármacos, Farmanguinhos - Fiocruz, Laboratório de Síntese Fármacos - LASFAR, Manguinhos Rio de Janeiro, RJ, Brazil.,Programa de Pós-Graduação em Farmacologia e Química Medicinal - Instituto de Ciências Biomédicas, Universidade Federal do Rio de Janeiro, UFRJ, Centro de Ciências da Saúde (CCS), Bloco J - Ilha do Fundão, Rio de Janeiro, RJ, Brazil
| | - Frederico Silva Castelo-Branco
- Fundação Oswaldo Cruz, Instituto de Tecnologia em Fármacos, Farmanguinhos - Fiocruz, Laboratório de Síntese Fármacos - LASFAR, Manguinhos Rio de Janeiro, RJ, Brazil
| | - Tácio Vinício Amorim Fernandes
- Fundação Oswaldo Cruz, Instituto de Tecnologia em Fármacos, Farmanguinhos - Fiocruz, Laboratório de Síntese Fármacos - LASFAR, Manguinhos Rio de Janeiro, RJ, Brazil
| | - Núbia Boechat
- Fundação Oswaldo Cruz, Instituto de Tecnologia em Fármacos, Farmanguinhos - Fiocruz, Laboratório de Síntese Fármacos - LASFAR, Manguinhos Rio de Janeiro, RJ, Brazil
| | - Luiza Rosaria Sousa Dias
- Faculdade de Farmácia, Laboratório de Química Medicinal - LQMed, Universidade Federal Fluminense (UFF), Santa Rosa, Niterói, RJ, Brazil
| | - Lucas Villas Bôas Hoelz
- Fundação Oswaldo Cruz, Instituto de Tecnologia em Fármacos, Farmanguinhos - Fiocruz, Laboratório de Síntese Fármacos - LASFAR, Manguinhos Rio de Janeiro, RJ, Brazil
| |
Collapse
|
7
|
Analysis of Integrin α IIb Subunit Dynamics Reveals Long-Range Effects of Missense Mutations on Calf Domains. Int J Mol Sci 2022; 23:ijms23020858. [PMID: 35055046 PMCID: PMC8776176 DOI: 10.3390/ijms23020858] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2021] [Revised: 12/23/2021] [Accepted: 12/30/2021] [Indexed: 11/17/2022] Open
Abstract
Integrin αIIbβ3, a glycoprotein complex expressed at the platelet surface, is involved in platelet aggregation and contributes to primary haemostasis. Several integrin αIIbβ3 polymorphisms prevent the aggregation that causes haemorrhagic syndromes, such as Glanzmann thrombasthenia (GT). Access to 3D structure allows understanding the structural effects of polymorphisms related to GT. In a previous analysis using Molecular Dynamics (MD) simulations of αIIbCalf-1 domain structure, it was observed that GT associated with single amino acid variation affects distant loops, but not the mutated position. In this study, experiments are extended to Calf-1, Thigh, and Calf-2 domains. Two loops in Calf-2 are unstructured and therefore are modelled expertly using biophysical restraints. Surprisingly, MD revealed the presence of rigid zones in these loops. Detailed analysis with structural alphabet, the Proteins Blocks (PBs), allowed observing local changes in highly flexible regions. The variant P741R located at C-terminal of Calf-1 revealed that the Calf-2 presence did not affect the results obtained with isolated Calf-1 domain. Simulations for Calf-1 + Calf-2, and Thigh + Calf-1 variant systems are designed to comprehend the impact of five single amino acid variations in these domains. Distant conformational changes are observed, thus highlighting the potential role of allostery in the structural basis of GT.
Collapse
|
8
|
Schwarz D, Georges G, Kelm S, Shi J, Vangone A, Deane CM. Co-evolutionary distance predictions contain flexibility information. Bioinformatics 2021; 38:65-72. [PMID: 34383892 DOI: 10.1093/bioinformatics/btab562] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Revised: 06/19/2021] [Accepted: 08/10/2021] [Indexed: 02/03/2023] Open
Abstract
MOTIVATION Co-evolution analysis can be used to accurately predict residue-residue contacts from multiple sequence alignments. The introduction of machine-learning techniques has enabled substantial improvements in precision and a shift from predicting binary contacts to predict distances between pairs of residues. These developments have significantly improved the accuracy of de novo prediction of static protein structures. With AlphaFold2 lifting the accuracy of some predicted protein models close to experimental levels, structure prediction research will move on to other challenges. One of those areas is the prediction of more than one conformation of a protein. Here, we examine the potential of residue-residue distance predictions to be informative of protein flexibility rather than simply static structure. RESULTS We used DMPfold to predict distance distributions for every residue pair in a set of proteins that showed both rigid and flexible behaviour. Residue pairs that were in contact in at least one reference structure were classified as rigid, flexible or neither. The predicted distance distribution of each residue pair was analysed for local maxima of probability indicating the most likely distance or distances between a pair of residues. We found that rigid residue pairs tended to have only a single local maximum in their predicted distance distributions while flexible residue pairs more often had multiple local maxima. These results suggest that the shape of predicted distance distributions contains information on the rigidity or flexibility of a protein and its constituent residues. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Dominik Schwarz
- Department of Statistics, University of Oxford, Oxford OX1 3LB, UK
| | - Guy Georges
- Department of Computational Engineering and Data Science, Large Molecule Research, Penzberg 82377, Germany
| | - Sebastian Kelm
- Computer-Aided Drug Design, UCB Pharma, Slough SL1 3WE, UK
| | - Jiye Shi
- Computer-Aided Drug Design, UCB Pharma, Slough SL1 3WE, UK
| | - Anna Vangone
- Department of Computational Engineering and Data Science, Large Molecule Research, Penzberg 82377, Germany
| | | |
Collapse
|
9
|
Peng W, Wang T, Liang XR, Yang YS, Wang QZ, Cheng HF, Peng YK, Ding F. Characterizing the potentially neuronal acetylcholinesterase reactivity toward chiral pyraclofos: Enantioselective insights from spectroscopy, in silico docking, molecular dynamics simulation and per-residue energy decomposition studies. J Mol Graph Model 2021; 110:108069. [PMID: 34773872 DOI: 10.1016/j.jmgm.2021.108069] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2021] [Revised: 10/27/2021] [Accepted: 10/28/2021] [Indexed: 10/20/2022]
Abstract
Chiral organophosphorus agents are distributed ubiquitously in the environment, but the neuroactivity of these asymmetric chemicals to humans remains uncertain. This scenario was to explore the stereoselective neurobiological response of human acetylcholinesterase (AChE) to chiral pyraclofos at the enantiomeric scale, and then decipher the microscopic basis of enantioselective neurotoxicity of pyraclofos enantiomers. The results indicated that (R)-/(S)-pyraclofos can form the bioconjugates with AChE with a stoichiometric ratio of 1:1, but the neuronal affinity of (R)-pyraclofos (K = 6.31 × 104 M-1) with AChE was larger than that of (S)-pyraclofos (K = 1.86 × 104 M-1), and significant enantioselectivity was existed in the biochemical reaction. The modes of neurobiological action revealed that pyraclofos enantiomers were situated at the substrate binding domain, and the strength of the overall noncovalent bonds between (S)-pyraclofos and the residues was weaker than that of (R)-pyraclofos, resulting in the high inhibitory effect of (R)-pyraclofos toward the activity of AChE. Dynamic enantioselective biointeractions illustrated that the intervention of inherent conformational flexibility in the AChE-(R)-pyraclofos was greater than that of the AChE-(S)-pyraclofos, which arises from the big spatial displacement and the conformational flip of the binding domain composed of the residues Thr-64~Asn-89, Gly-122~Asp-134, and Thr-436~Tyr-449. Energy decomposition exhibited that the Gibbs free energies of the AChE-(R)-/(S)-pyraclofos were ΔG° = -37.4/-30.2 kJ mol-1, respectively, and the disparity comes from the electrostatic energy during the stereoselective neurochemical reactions. Quantitative conformational analysis further confirmed the atomic-scale computational chemistry conclusions, and the perturbation of (S)-pyraclofos on the AChE's ordered conformation was lower than that of (R)-pyraclofos, which is germane to the interaction energies of the crucial residues, e.g. Tyr-124, Tyr-337, Asp-74, Trp-86, and Tyr-119. Evidently, this attempt will contribute mechanistic information to uncovering the neurobiological effects of chiral organophosphates on the body.
Collapse
Affiliation(s)
- Wei Peng
- School of Water and Environment, Chang'an University, Xi'an, 710054, China; Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education, Chang'an University, Xi'an, 710054, China; State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, 361005, China
| | - Tao Wang
- School of Water and Environment, Chang'an University, Xi'an, 710054, China; Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education, Chang'an University, Xi'an, 710054, China
| | - Xiang-Rong Liang
- School of Water and Environment, Chang'an University, Xi'an, 710054, China
| | - Yu-Sen Yang
- School of Water and Environment, Chang'an University, Xi'an, 710054, China
| | - Qi-Zhao Wang
- School of Water and Environment, Chang'an University, Xi'an, 710054, China; Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education, Chang'an University, Xi'an, 710054, China
| | - Hong-Fei Cheng
- School of Earth Science and Resources, Chang'an University, Xi'an, 710054, China
| | - Yu-Kui Peng
- Xining Center for Agricultural Product Quality and Safety Testing, Xining, 810016, China
| | - Fei Ding
- School of Water and Environment, Chang'an University, Xi'an, 710054, China; Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education, Chang'an University, Xi'an, 710054, China; Department of Agricultural Chemistry, Qingdao Agricultural University, Qingdao, 266109, China.
| |
Collapse
|
10
|
El Jahrani N, Cretin G, de Brevern AG. CALR-ETdb, the database of calreticulin variants diversity in essential thrombocythemia. Platelets 2021; 33:157-167. [PMID: 33444113 DOI: 10.1080/09537104.2020.1869712] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Abstract
Essential thrombocythemia (ET) is a blood cancer defined by a strong increase of platelet numbers. A quarter of patients suffering from ET show mutations in the last exon of calreticulin (CALR) gene. Two variants named type 1 and type 2 represent 85% of these patients. However, a large number of other variants have been determined. In this study, we have compiled variants taken from COSMIC database and literature leading to 155 different variants. This large number of variants allowed redefining 5 new classes extending the classification of type 1-like and type 2-like to a finer description. These analyses showed that last class, named E, corresponding to more than 10% of CALR variants seemed not attached to ET. Structural properties analyzed showed that CALR variants associated to ET have common features. All the compiled and refined information had been included into a freely dedicated database CALR-ETdb (https://www.dsimb.inserm.fr/CALR-ET).
Collapse
Affiliation(s)
- Nora El Jahrani
- Université de Paris, UMR_S 1134, Université De La Réunion, Université Des Antilles, Paris, France.,INSERM, U 1134, DSIMB, Paris, France.,Institut National De La Transfusion Sanguine (INTS), Paris, France.,Laboratoire d'Excellence GR-Ex, Paris, France
| | - Gabriel Cretin
- Université de Paris, UMR_S 1134, Université De La Réunion, Université Des Antilles, Paris, France.,INSERM, U 1134, DSIMB, Paris, France.,Institut National De La Transfusion Sanguine (INTS), Paris, France.,Laboratoire d'Excellence GR-Ex, Paris, France
| | - Alexandre G de Brevern
- Université de Paris, UMR_S 1134, Université De La Réunion, Université Des Antilles, Paris, France.,INSERM, U 1134, DSIMB, Paris, France.,Institut National De La Transfusion Sanguine (INTS), Paris, France.,Laboratoire d'Excellence GR-Ex, Paris, France
| |
Collapse
|
11
|
de Brevern AG. Impact of protein dynamics on secondary structure prediction. Biochimie 2020; 179:14-22. [PMID: 32946990 DOI: 10.1016/j.biochi.2020.09.006] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2020] [Revised: 09/04/2020] [Accepted: 09/10/2020] [Indexed: 02/08/2023]
Abstract
Protein 3D structures support their biological functions. As the number of protein structures is negligible in regards to the number of available protein sequences, prediction methodologies relying only on protein sequences are essential tools. In this field, protein secondary structure prediction (PSSPs) is a mature area, and is considered to have reached a plateau. Nonetheless, proteins are highly dynamical macromolecules, a property that could impact the PSSP methods. Indeed, in a previous study, the stability of local protein conformations was evaluated demonstrating that some regions easily changed to another type of secondary structure. The protein sequences of this dataset were used by PSSPs and their results compared to molecular dynamics to investigate their potential impact on the quality of the secondary structure prediction. Interestingly, a direct link is observed between the quality of the prediction and the stability of the assignment to the secondary structure state. The more stable a local protein conformation is, the better the prediction will be. The secondary structure assignment not taken from the crystallized structures but from the conformations observed during the dynamics slightly increase the quality of the secondary structure prediction. These results show that evaluation of PSSPs can be done differently, but also that the notion of dynamics can be included in development of PSSPs and other approaches such as de novo approaches.
Collapse
Affiliation(s)
- Alexandre G de Brevern
- Biologie Intégrée Du Globule Rouge UMR_S1134, Inserm, Université de Paris, Univ. de la Réunion, Univ. des Antilles, F-75739, Paris, France; Laboratoire D'Excellence GR-Ex, F-75739, Paris, France; Institut National de la Transfusion Sanguine (INTS), F-75739, Paris, France; IBL, F-75015, Paris, France.
| |
Collapse
|
12
|
de Brevern AG. Analysis of Protein Disorder Predictions in the Light of a Protein Structural Alphabet. Biomolecules 2020; 10:biom10071080. [PMID: 32698546 PMCID: PMC7408373 DOI: 10.3390/biom10071080] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Revised: 07/14/2020] [Accepted: 07/18/2020] [Indexed: 12/30/2022] Open
Abstract
Intrinsically-disordered protein (IDP) characterization was an amazing change of paradigm in our classical sequence-structure-function theory. Moreover, IDPs are over-represented in major disease pathways and are now often targeted using small molecules for therapeutic purposes. This has had created a complex continuum from order-that encompasses rigid and flexible regions-to disorder regions; the latter being not accessible through classical crystallographic methodologies. In X-ray structures, the notion of order is dictated by access to resolved atom positions, providing rigidity and flexibility information with low and high experimental B-factors, while disorder is associated with the missing (non-resolved) residues. Nonetheless, some rigid regions can be found in disorder regions. Using ensembles of IDPs, their local conformations were analyzed in the light of a structural alphabet. An entropy index derived from this structural alphabet allowed us to propose a continuum of states from rigidity to flexibility and finally disorder. In this study, the analysis was extended to comparing these results to disorder predictions, underlying a limited correlation, and so opening new ideas to characterize and predict disorder.
Collapse
Affiliation(s)
- Alexandre G de Brevern
- INSERM, UMR_S 1134, DSIMB, Univ Paris, INTS, Laboratoire d'Excellence GR-Ex, 75015 Paris, France
| |
Collapse
|