1
|
Hasan A, Alonazi WB, Ibrahim M, Bin L. Immunoinformatics and Reverse Vaccinology Approach for the Identification of Potential Vaccine Candidates against Vandammella animalimors. Microorganisms 2024; 12:1270. [PMID: 39065039 PMCID: PMC11278545 DOI: 10.3390/microorganisms12071270] [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: 04/29/2024] [Revised: 06/17/2024] [Accepted: 06/19/2024] [Indexed: 07/28/2024] Open
Abstract
Vandammella animalimorsus is a Gram-negative and non-motile bacterium typically transmitted to humans through direct contact with the saliva of infected animals, primarily through biting, scratches, or licks on fractured skin. The absence of a confirmed post-exposure treatment of V. animalimorsus bacterium highlights the imperative for developing an effective vaccine. We intended to determine potential vaccine candidates and paradigm a chimeric vaccine against V. animalimorsus by accessible public data analysis of the strain by utilizing reverse vaccinology. By subtractive genomics, five outer membranes were prioritized as potential vaccine candidates out of 2590 proteins. Based on the instability index and transmembrane helices, a multidrug transporter protein with locus ID A0A2A2AHJ4 was designated as a potential candidate for vaccine construct. Sixteen immunodominant epitopes were retrieved by utilizing the Immune Epitope Database. The epitope encodes the strong binding affinity, nonallergenic properties, non-toxicity, high antigenicity scores, and high solubility revealing the more appropriate vaccine construct. By utilizing appropriate linkers and adjuvants alongside a suitable adjuvant molecule, the epitopes were integrated into a chimeric vaccine to enhance immunogenicity, successfully eliciting both adaptive and innate immune responses. Moreover, the promising physicochemical features, the binding confirmation of the vaccine to the major innate immune receptor TLR-4, and molecular dynamics simulations of the designed vaccine have revealed the promising potential of the selected candidate. The integration of computational methods and omics data has demonstrated significant advantages in discovering novel vaccine targets and mitigating vaccine failure rates during clinical trials in recent years.
Collapse
Affiliation(s)
- Ahmad Hasan
- State Key Laboratory of Rice Biology and Breeding, Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Biotechnology, Zhejiang University, Hangzhou 310058, China; (A.H.); (M.I.)
| | - Wadi B. Alonazi
- Health Administration Department, College of Business Administration, King Saud University, Riyadh 11421, Saudi Arabia;
| | - Muhammad Ibrahim
- State Key Laboratory of Rice Biology and Breeding, Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Biotechnology, Zhejiang University, Hangzhou 310058, China; (A.H.); (M.I.)
| | - Li Bin
- State Key Laboratory of Rice Biology and Breeding, Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Biotechnology, Zhejiang University, Hangzhou 310058, China; (A.H.); (M.I.)
| |
Collapse
|
2
|
Subhadarshini S, Tandon H, Srinivasan N, Sowdhamini R. Normal Mode Analysis Elicits Conformational Shifts in Proteins at Both Proximal and Distal Regions to the Phosphosite Stemming from Single-Site Phosphorylation. ACS OMEGA 2024; 9:24520-24537. [PMID: 38882086 PMCID: PMC11170700 DOI: 10.1021/acsomega.4c00523] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/16/2024] [Revised: 04/29/2024] [Accepted: 05/13/2024] [Indexed: 06/18/2024]
Abstract
Phosphorylation, a fundamental biochemical switch, intricately regulates protein function and signaling pathways. Our study employs extensive computational structural analyses on a curated data set of phosphorylated and unphosphorylated protein pairs to explore the multifaceted impact of phosphorylation on protein conformation. Using normal mode analysis (NMA), we investigated changes in protein flexibility post-phosphorylation, highlighting an enhanced level of structural dynamism. Our findings reveal that phosphorylation induces not only local changes at the phosphorylation site but also extensive alterations in distant regions, showcasing its far-reaching influence on protein structure-dynamics. Through in-depth case studies on polyubiquitin B and glycogen synthase kinase-3 beta, we elucidate how phosphorylation at distinct sites leads to variable structural and dynamic modifications, potentially dictating functional outcomes. While phosphorylation largely preserves the residue motion correlation, it significantly disrupts low-frequency global modes, presenting a dualistic impact on protein dynamics. We also explored alterations in the total accessible surface area (ASA), emphasizing region-specific changes around phosphorylation sites. This study sheds light on phosphorylation-induced conformational changes, dynamic modulation, and surface accessibility alterations, leveraging an integrated computational approach with RMSD, NMA, and ASA, thereby contributing to a comprehensive understanding of cellular regulation and suggesting promising avenues for therapeutic interventions.
Collapse
Affiliation(s)
| | - Himani Tandon
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore 560012, India
| | | | - Ramanathan Sowdhamini
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore 560012, India
- Computational Approaches to Protein Science, National Centre for Biological Sciences, Bangalore 560065, India
- Computational Biology, Institute of Bioinformatics and Applied Biotechnology, Bangalore 560100, India
| |
Collapse
|
3
|
Kolossváry I. A Fresh Look at the Normal Mode Analysis of Proteins: Introducing Allosteric Co-Vibrational Modes. JACS AU 2024; 4:1303-1309. [PMID: 38665643 PMCID: PMC11040550 DOI: 10.1021/jacsau.4c00109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/04/2024] [Revised: 03/28/2024] [Accepted: 03/29/2024] [Indexed: 04/28/2024]
Abstract
We propose a new way of utilizing normal modes to study protein conformational transitions. Instead of considering individual modes independently, we show that a weighted mixture of low-frequency vibrational modes can reveal dynamic information about the conformational mechanism in more detail than any single mode can. The weights in the mixed mode, termed the allosteric covibrational mode, are determined using a simple model where the conformational transition is viewed as a perturbation of the coupled harmonic oscillator associated with either of the two conformations. We demonstrate our theory in a biologically relevant example of high pharmaceutical interest involving the V617F mutation of Janus 2 tyrosine kinase (JAK2).
Collapse
Affiliation(s)
- István Kolossváry
- Department of Biomedical
Engineering, Boston University, Boston, Massachusetts 02215, United States
| |
Collapse
|
4
|
Ghaffar SA, Tahir H, Muhammad S, Shahid M, Naqqash T, Faisal M, Albekairi TH, Alshammari A, Albekairi NA, Manzoor I. Designing of a multi-epitopes based vaccine against Haemophilius parainfluenzae and its validation through integrated computational approaches. Front Immunol 2024; 15:1380732. [PMID: 38690283 PMCID: PMC11058264 DOI: 10.3389/fimmu.2024.1380732] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Accepted: 04/01/2024] [Indexed: 05/02/2024] Open
Abstract
Haemophilus parainfluenzae is a Gram-negative opportunist pathogen within the mucus of the nose and mouth without significant symptoms and has an ability to cause various infections ranging from ear, eye, and sinus to pneumonia. A concerning development is the increasing resistance of H. parainfluenzae to beta-lactam antibiotics, with the potential to cause dental infections or abscesses. The principal objective of this investigation is to utilize bioinformatics and immuno-informatic methodologies in the development of a candidate multi-epitope Vaccine. The investigation focuses on identifying potential epitopes for both B cells (B lymphocytes) and T cells (helper T lymphocytes and cytotoxic T lymphocytes) based on high non-toxic and non-allergenic characteristics. The selection process involves identifying human leukocyte antigen alleles demonstrating strong associations with recognized antigenic and overlapping epitopes. Notably, the chosen alleles aim to provide coverage for 90% of the global population. Multi-epitope constructs were designed by using suitable linker sequences. To enhance the immunological potential, an adjuvant sequence was incorporated using the EAAAK linker. The final vaccine construct, comprising 344 amino acids, was achieved after the addition of adjuvants and linkers. This multi-epitope Vaccine demonstrates notable antigenicity and possesses favorable physiochemical characteristics. The three-dimensional conformation underwent modeling and refinement, validated through in-silico methods. Additionally, a protein-protein molecular docking analysis was conducted to predict effective binding poses between the multi-epitope Vaccine and the Toll-like receptor 4 protein. The Molecular Dynamics (MD) investigation of the docked TLR4-vaccine complex demonstrated consistent stability over the simulation period, primarily attributed to electrostatic energy. The docked complex displayed minimal deformation and enhanced rigidity in the motion of residues during the dynamic simulation. Furthermore, codon translational optimization and computational cloning was performed to ensure the reliability and proper expression of the multi-Epitope Vaccine. It is crucial to emphasize that despite these computational validations, experimental research in the laboratory is imperative to demonstrate the immunogenicity and protective efficacy of the developed vaccine. This would involve practical assessments to ascertain the real-world effectiveness of the multi-epitope Vaccine.
Collapse
Affiliation(s)
- Sana Abdul Ghaffar
- Department of Bioinformatics and Biotechnology, Government College University, Faisalabad, Pakistan
| | - Haneen Tahir
- Department of Bioinformatics and Biotechnology, Government College University, Faisalabad, Pakistan
| | - Sher Muhammad
- Department of Bioinformatics and Biotechnology, Government College University, Faisalabad, Pakistan
| | - Muhammad Shahid
- Department of Bioinformatics and Biotechnology, Government College University, Faisalabad, Pakistan
| | - Tahir Naqqash
- Institute of Molecular Biology and Biotechnology, Bahauddin Zakariya University, Multan, Pakistan
| | | | - Thamer H. Albekairi
- Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Abdulrahman Alshammari
- Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Norah A. Albekairi
- Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Irfan Manzoor
- Department of Bioinformatics and Biotechnology, Government College University, Faisalabad, Pakistan
- Department of Biology, Indiana University, Bloomington, IN, United States
| |
Collapse
|
5
|
Araki K, Watanabe-Nakayama T, Sasaki D, Sasaki YC, Mio K. Molecular Dynamics Mappings of the CCT/TRiC Complex-Mediated Protein Folding Cycle Using Diffracted X-ray Tracking. Int J Mol Sci 2023; 24:14850. [PMID: 37834298 PMCID: PMC10573753 DOI: 10.3390/ijms241914850] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2023] [Revised: 09/27/2023] [Accepted: 09/29/2023] [Indexed: 10/15/2023] Open
Abstract
The CCT/TRiC complex is a type II chaperonin that undergoes ATP-driven conformational changes during its functional cycle. Structural studies have provided valuable insights into the mechanism of this process, but real-time dynamics analyses of mammalian type II chaperonins are still scarce. We used diffracted X-ray tracking (DXT) to investigate the intramolecular dynamics of the CCT complex. We focused on three surface-exposed loop regions of the CCT1 subunit: the loop regions of the equatorial domain (E domain), the E and intermediate domain (I domain) juncture near the ATP-binding region, and the apical domain (A domain). Our results showed that the CCT1 subunit predominantly displayed rotational motion, with larger mean square displacement (MSD) values for twist (χ) angles compared with tilt (θ) angles. Nucleotide binding had a significant impact on the dynamics. In the absence of nucleotides, the region between the E and I domain juncture could act as a pivotal axis, allowing for greater motion of the E domain and A domain. In the presence of nucleotides, the nucleotides could wedge into the ATP-binding region, weakening the role of the region between the E and I domain juncture as the rotational axis and causing the CCT complex to adopt a more compact structure. This led to less expanded MSD curves for the E domain and A domain compared with nucleotide-absent conditions. This change may help to stabilize the functional conformation during substrate binding. This study is the first to use DXT to probe the real-time molecular dynamics of mammalian type II chaperonins at the millisecond level. Our findings provide new insights into the complex dynamics of chaperonins and their role in the functional folding cycle.
Collapse
Affiliation(s)
- Kazutaka Araki
- AIST-UTokyo Advanced Operando-Measurement Technology Open Innovation Laboratory (OPERANDO-OIL), National Institute of Advanced Industrial Science and Technology (AIST), 6-2-3 Kashiwanoha, Chiba 277-0882, Japan;
| | - Takahiro Watanabe-Nakayama
- WPI Nano Life Science Institute (WPI-NanoLSI), Kanazawa University, Kakuma-machi, Kanazawa 920-1192, Japan;
| | - Daisuke Sasaki
- Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5 Kashiwanoha, Chiba 277-8561, Japan (Y.C.S.)
| | - Yuji C. Sasaki
- Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5 Kashiwanoha, Chiba 277-8561, Japan (Y.C.S.)
| | - Kazuhiro Mio
- AIST-UTokyo Advanced Operando-Measurement Technology Open Innovation Laboratory (OPERANDO-OIL), National Institute of Advanced Industrial Science and Technology (AIST), 6-2-3 Kashiwanoha, Chiba 277-0882, Japan;
| |
Collapse
|
6
|
Khan MW, Murali A. Normal mode analysis and comparative study of intrinsic dynamics of alcohol oxidase enzymes from GMC protein family. J Biomol Struct Dyn 2023:1-16. [PMID: 37676256 DOI: 10.1080/07391102.2023.2255275] [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: 06/09/2023] [Accepted: 08/30/2023] [Indexed: 09/08/2023]
Abstract
Glucose-Methanol-Choline (GMC) family enzymes are very important in catalyzing the oxidation of a wide range of structurally diverse substrates. Enzymes that constitute the GMC family, share a common tertiary fold but < 25% sequence identity. Cofactor FAD, FAD binding signature motif, and similar structural scaffold of the active site are common features of oxidoreductase enzymes of the GMC family. Protein functionality mainly depends on protein three-dimensional structures and dynamics. In this study, we used the normal mode analysis method to search the intrinsic dynamics of GMC family enzymes. We have explored the dynamical behavior of enzymes with unique substrate catabolism and active site characteristics from different classes of the GMC family. Analysis of individual enzymes and comparative ensemble analysis of enzymes from different classes has shown conserved dynamic motion at FAD binding sites. The present study revealed that GMC enzymes share a strong dynamic similarity (Bhattacharyya coefficient >90% and root mean squared inner product >52%) despite low sequence identity across the GMC family enzymes. The study predicts that local deformation energy between atoms of the enzyme may be responsible for the catalysis of different substrates. This study may help that intrinsic dynamics can be used to make meaningful classifications of proteins or enzymes from different organisms.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Mohammad Wahab Khan
- Department of Bioinformatics, School of Life Science, Pondicherry University, Puducherry, India
| | - Ayaluru Murali
- Department of Bioinformatics, School of Life Science, Pondicherry University, Puducherry, India
| |
Collapse
|
7
|
Motamedi H, Alvandi A, Fathollahi M, Ari MM, Moradi S, Moradi J, Abiri R. In silico designing and immunoinformatics analysis of a novel peptide vaccine against metallo-beta-lactamase (VIM and IMP) variants. PLoS One 2023; 18:e0275237. [PMID: 37471423 PMCID: PMC10358925 DOI: 10.1371/journal.pone.0275237] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2022] [Accepted: 09/12/2022] [Indexed: 07/22/2023] Open
Abstract
The rapid spread of acquired metallo-beta-lactamases (MBLs) among gram negative pathogens is becoming a global concern. Improper use of broad-spectrum antibiotics can trigger the colonization and spread of resistant strains which lead to increased mortality and significant economic loss. In the present study, diverse immunoinformatic approaches are applied to design a potential epitope-based vaccine against VIM and IMP MBLs. The amino acid sequences of VIM and IMP variants were retrieved from the GenBank database. ABCpred and BCPred online Web servers were used to analyze linear B cell epitopes, while IEDB was used to determine the dominant T cell epitopes. Sequence validation, allergenicity, toxicity and physiochemical analysis were performed using web servers. Seven sequences were identified for linear B cell dominant epitopes and 4 sequences were considered as dominant CD4+ T cell epitopes, and the predicted epitopes were joined by KK and GPGPG linkers. Stabilized multi-epitope protein structure was obtained using molecular dynamics simulation. Molecular docking showed that the designed vaccine exhibited sustainable and strong binding interactions with Toll-like receptor 4 (TLR4). Finally, codon adaptation and in silico cloning studies were performed to design an effective vaccine production strategy. Immune simulation significantly provided high levels of immunoglobulins, T helper cells, T-cytotoxic cells and INF-γ. Even though the introduced vaccine candidate demonstrates a very potent immunogenic potential, but wet-lab validation is required to further assessment of the effectiveness of this proposed vaccine candidate.
Collapse
Affiliation(s)
- Hamid Motamedi
- Department of Microbiology, School of Medicine, Kermanshah University of Medical Sciences, Kermanshah, Iran
- Student Research Committee, School of Medicine, Kermanshah University of Medical Sciences, Kermanshah, Iran
| | - Amirhoushang Alvandi
- Medical Technology Research Center, Health Technology Institute, Kermanshah University of Medical Sciences, Kermanshah, Iran
| | - Matin Fathollahi
- Department of Microbiology, School of Medicine, Kermanshah University of Medical Sciences, Kermanshah, Iran
- Student Research Committee, School of Medicine, Kermanshah University of Medical Sciences, Kermanshah, Iran
| | - Marzie Mahdizade Ari
- Department of Microbiology, School of Medicine, Iran University of Medical Sciences, Tehran, Iran
- Microbial Biotechnology Research Centre, Iran University of Medical Sciences, Tehran, Iran
| | - Sajad Moradi
- Nano Drug Delivery Research Center, Health Technology Institute, Kermanshah University of Medical Sciences, Kermanshah, Iran
| | - Jale Moradi
- Department of Microbiology, School of Medicine, Kermanshah University of Medical Sciences, Kermanshah, Iran
| | - Ramin Abiri
- Department of Microbiology, School of Medicine, Kermanshah University of Medical Sciences, Kermanshah, Iran
- Fertility and Infertility Research Center, Health Technology Institute, Kermanshah University of Medical Sciences, Kermanshah, Iran
| |
Collapse
|
8
|
Gaur NK, Ghosh B, Goyal VD, Kulkarni K, Makde RD. Evolutionary conservation of protein dynamics: insights from all-atom molecular dynamics simulations of 'peptidase' domain of Spt16. J Biomol Struct Dyn 2023; 41:1445-1457. [PMID: 34971347 DOI: 10.1080/07391102.2021.2021990] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Protein function is encoded in its sequence, manifested in its three-dimensional structure, and facilitated by its dynamics. Studies have suggested that protein structures with higher sequence similarity could have more similar patterns of dynamics. However, such studies of protein dynamics within and across protein families typically rely on coarse-grained models, or approximate metrics like crystallographic B-factors. This study uses µs scale molecular dynamics (MD) simulations to explore the conservation of dynamics among homologs of ∼50 kDa N-terminal module of Spt16 (Spt16N). Spt16N from Saccharomyces cerevisiae (Sc-Spt16N) and three of its homologs with 30-40% sequence identities were available in the PDB. To make our data-set more comprehensive, the crystal structure of an additional homolog (62% sequence identity with Sc-Spt16N) was solved at 1.7 Å resolution. Cumulative MD simulations of 6 µs were carried out on these Spt16N structures and on two additional protein structures with varying degrees of similarity to it. The simulations revealed that correlation in patterns of backbone fluctuations vary linearly with sequence identity. This trend could not be inferred using crystallographic B-factors. Further, normal mode analysis suggested a similar pattern of inter-domain (inter-lobe) motions not only among Spt16N homologs, but also in the M24 peptidase structure. On the other hand, MD simulation results highlighted conserved motions that were found unique for Spt16N protein, this along with electrostatics trends shed light on functional aspects of Spt16N.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Neeraj K Gaur
- Beamline Development and Application Section, Bhabha Atomic Research Centre, Mumbai, India.,Division of Biochemical Sciences, CSIR-National Chemical Laboratory, Pune, India.,Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
| | - Biplab Ghosh
- Beamline Development and Application Section, Bhabha Atomic Research Centre, Mumbai, India
| | - Venuka Durani Goyal
- Beamline Development and Application Section, Bhabha Atomic Research Centre, Mumbai, India
| | - Kiran Kulkarni
- Division of Biochemical Sciences, CSIR-National Chemical Laboratory, Pune, India.,Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
| | - Ravindra D Makde
- Beamline Development and Application Section, Bhabha Atomic Research Centre, Mumbai, India
| |
Collapse
|
9
|
Chandrasekhar G, Chandra Sekar P, Srinivasan E, Amarnath A, Pengyong H, Rajasekaran R. Molecular simulation unravels the amyloidogenic misfolding of nascent ApoA1 protein, driven by deleterious point mutations occurring in between 170-178 hotspot region. J Biomol Struct Dyn 2022; 40:13278-13290. [PMID: 34613891 DOI: 10.1080/07391102.2021.1986134] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Protein ApoA1 is extensively studied for its role in lipid metabolism. Its seedy dark side of amyloid formulation remains relatively understudied yet. Due to genetic mutations, the protein pathologically misshapes into its amyloid form that gets accumulated in various organs, including the heart. To contrive effective therapeutics against this debilitating congenital disorder, it is imperative to comprehend the structural ramifications induced by mutations in APoA1's dynamic conformation. Till now, several point mutations have been implicated in ApoA1's amyloidosis, although only a handful has been examined considerably. Especially, the single nucleotide polymorphisms (SNPs) that occur in-between 170-178 mutation hotspot site of APoA1 needs to be investigated, since most of them are culpable of amyloid deposition in the heart. To that effect, in the present study, we have computationally quantified and studied the ApoA1's biomolecular modifications fostered by SNPs in the 170-178 mutation hotspot. Findings from discrete molecular dynamics simulation studies indicate that the SNPs have noticeably steered the ApoA1's behaviour from its native structural dynamics. Analysis of protein's secondary structural changes exhibits a considerable change upon mutations. Further, subjecting the protein structures to simulated thermal denaturation shows increased resistance to denaturation among mutants when compared to native. Further, normal mode analysis of protein's dynamic motion also shows discrepancy in its dynamic structural change upon SNP. These structural digressions induced by SNPs can very well be the biomolecular incendiary that drives ApoA1 into its amyloidogenesis. And, understanding these structural modifications initiates a better understanding of SNP's amyloidogenic pathology on APoA1.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- G Chandrasekhar
- Bioinformatics Lab, Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology (Deemed to be University), Vellore, Tamil Nadu, India
| | - P Chandra Sekar
- Bioinformatics Lab, Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology (Deemed to be University), Vellore, Tamil Nadu, India
| | - E Srinivasan
- Bioinformatics Lab, Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology (Deemed to be University), Vellore, Tamil Nadu, India
| | - A Amarnath
- Bioinformatics Lab, Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology (Deemed to be University), Vellore, Tamil Nadu, India
| | - H Pengyong
- Central Lab, Changzhi Medical College, Changzhi, China
| | - R Rajasekaran
- Bioinformatics Lab, Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology (Deemed to be University), Vellore, Tamil Nadu, India
| |
Collapse
|
10
|
Mertens HDT. Computational methods for the analysis of solution small-angle X-ray scattering of biomolecules: ATSAS. Methods Enzymol 2022; 678:193-236. [PMID: 36641208 DOI: 10.1016/bs.mie.2022.09.019] [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: 01/04/2023]
Abstract
The ATSAS software suite provides a comprehensive set of programs for the processing, analysis and modeling of small-angle scattering data, tailored for but not limited to data acquired on biological macromolecules. In this review the major components and developments in the ATSAS package are described, with a focus on user driven application. Data reduction, analysis and modeling approaches and strategies will be introduced and discussed. At the time of writing the latest package, ATSAS 3.1, is freely available for academic users at: https://www.embl-hamburg.de/biosaxs/software.html.
Collapse
|
11
|
Periwal N, Rathod SB, Sarma S, Johar GS, Jain A, Barnwal RP, Srivastava KR, Kaur B, Arora P, Sood V. Time Series Analysis of SARS-CoV-2 Genomes and Correlations among Highly Prevalent Mutations. Microbiol Spectr 2022; 10:e0121922. [PMID: 36069583 PMCID: PMC9603882 DOI: 10.1128/spectrum.01219-22] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Accepted: 08/03/2022] [Indexed: 12/30/2022] Open
Abstract
The efforts of the scientific community to tame the recent pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) seem to have been diluted by the emergence of new viral strains. Therefore, it is imperative to understand the effect of mutations on viral evolution. We performed a time series analysis on 59,541 SARS-CoV-2 genomic sequences from around the world to gain insights into the kinetics of the mutations arising in the viral genomes. These 59,541 genomes were grouped according to month (January 2020 to March 2021) based on the collection date. Meta-analysis of these data led us to identify significant mutations in viral genomes. Pearson correlation of these mutations led us to the identification of 16 comutations. Among these comutations, some of the individual mutations have been shown to contribute to viral replication and fitness, suggesting a possible role of other unexplored mutations in viral evolution. We observed that the mutations 241C>T in the 5' untranslated region (UTR), 3037C>T in nsp3, 14408C>T in the RNA-dependent RNA polymerase (RdRp), and 23403A>G in spike are correlated with each other and were grouped in a single cluster by hierarchical clustering. These mutations have replaced the wild-type nucleotides in SARS-CoV-2 sequences. Additionally, we employed a suite of computational tools to investigate the effects of T85I (1059C>T), P323L (14408C>T), and Q57H (25563G>T) mutations in nsp2, RdRp, and the ORF3a protein of SARS-CoV-2, respectively. We observed that the mutations T85I and Q57H tend to be deleterious and destabilize the respective wild-type protein, whereas P323L in RdRp tends to be neutral and has a stabilizing effect. IMPORTANCE We performed a meta-analysis on SARS-CoV-2 genomes categorized by collection month and identified several significant mutations. Pearson correlation analysis of these significant mutations identified 16 comutations having absolute correlation coefficients of >0.4 and a frequency of >30% in the genomes used in this study. The correlation results were further validated by another statistical tool called hierarchical clustering, where mutations were grouped in clusters on the basis of their similarity. We identified several positive and negative correlations among comutations in SARS-CoV-2 isolates from around the world which might contribute to viral pathogenesis. The negative correlations among some of the mutations in SARS-CoV-2 identified in this study warrant further investigations. Further analysis of mutations such as T85I in nsp2 and Q57H in ORF3a protein revealed that these mutations tend to destabilize the protein relative to the wild type, whereas P323L in RdRp is neutral and has a stabilizing effect. Thus, we have identified several comutations which can be further characterized to gain insights into SARS-CoV-2 evolution.
Collapse
Affiliation(s)
- Neha Periwal
- Department of Biochemistry, SCLS, Jamia Hamdard, New Delhi, India
| | - Shravan B. Rathod
- Department of Chemistry, Smt. S. M. Panchal Science College, Talod, Gujarat, India
| | - Sankritya Sarma
- Department of Zoology, Hansraj College, University of Delhi, New Delhi, India
| | | | - Avantika Jain
- Department of Biochemistry, SCLS, Jamia Hamdard, New Delhi, India
- Delhi Institute of Pharmaceutical Sciences and Research, New Delhi, Delhi, India
| | - Ravi P. Barnwal
- Department of Biophysics, Panjab University, Chandigarh, India
| | | | - Baljeet Kaur
- Department of Computer Science, Hansraj College, University of Delhi, New Delhi, India
| | - Pooja Arora
- Department of Zoology, Hansraj College, University of Delhi, New Delhi, India
| | - Vikas Sood
- Department of Biochemistry, SCLS, Jamia Hamdard, New Delhi, India
| |
Collapse
|
12
|
Tang QY, Ren W, Wang J, Kaneko K. The Statistical Trends of Protein Evolution: A Lesson from AlphaFold Database. Mol Biol Evol 2022; 39:msac197. [PMID: 36108094 PMCID: PMC9550990 DOI: 10.1093/molbev/msac197] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The recent development of artificial intelligence provides us with new and powerful tools for studying the mysterious relationship between organism evolution and protein evolution. In this work, based on the AlphaFold Protein Structure Database (AlphaFold DB), we perform comparative analyses of the proteins of different organisms. The statistics of AlphaFold-predicted structures show that, for organisms with higher complexity, their constituent proteins will have larger radii of gyration, higher coil fractions, and slower vibrations, statistically. By conducting normal mode analysis and scaling analyses, we demonstrate that higher organismal complexity correlates with lower fractal dimensions in both the structure and dynamics of the constituent proteins, suggesting that higher functional specialization is associated with higher organismal complexity. We also uncover the topology and sequence bases of these correlations. As the organismal complexity increases, the residue contact networks of the constituent proteins will be more assortative, and these proteins will have a higher degree of hydrophilic-hydrophobic segregation in the sequences. Furthermore, by comparing the statistical structural proximity across the proteomes with the phylogenetic tree of homologous proteins, we show that, statistical structural proximity across the proteomes may indirectly reflect the phylogenetic proximity, indicating a statistical trend of protein evolution in parallel with organism evolution. This study provides new insights into how the diversity in the functionality of proteins increases and how the dimensionality of the manifold of protein dynamics reduces during evolution, contributing to the understanding of the origin and evolution of lives.
Collapse
Affiliation(s)
- Qian-Yuan Tang
- Laboratory for Neural Computation and Adaptation, RIKEN Center for Brain Science, 2-1 Hirosawa, Wako, Saitama 351-0106, Japan
| | - Weitong Ren
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Jun Wang
- School of Physics, National Laboratory of Solid State Microstructure, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Kunihiko Kaneko
- Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Komaba, Meguro, Tokyo 153-8902, Japan
- The Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, Copenhagen 2100-DK, Denmark
| |
Collapse
|
13
|
Chakraborty C, Bhattacharya M, Sharma AR, Mallik B. Omicron (B.1.1.529) - A new heavily mutated variant: Mapped location and probable properties of its mutations with an emphasis on S-glycoprotein. Int J Biol Macromol 2022; 219:980-997. [PMID: 35952818 PMCID: PMC9359758 DOI: 10.1016/j.ijbiomac.2022.07.254] [Citation(s) in RCA: 33] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2022] [Revised: 05/23/2022] [Accepted: 07/31/2022] [Indexed: 12/17/2022]
Abstract
Omicron, another SARS-CoV-2 variant, has been recorded and reported as a VoC. It has already spread across >30 countries and is a highly mutated variant. We tried to understand the role of mutations in the investigated variants by comparison with previous characterized VoC. We have mapped the mutations in Omicron S-glycoprotein's secondary and tertiary structure landscape using bioinformatics tools and statistical software and developed different models. In addition, we analyzed the effect of diverse mutations in antibody binding regions of the S-glycoprotein on the binding affinity of the investigated antibodies. This study has chosen eight significant mutations in Omicron (D614G, E484A, N501Y, Q493K, K417N, S477N, Y505H G496S), and seven of them are located in the RBD region. We also performed a comparative analysis of the ΔΔG score of these mutations to understand the stabilizing or destabilizing properties of the investigated mutations. The analysis outcome shows that D614G, Q493K, and S477N mutations are stable mutations with ΔΔG scores of 0.351 kcal/mol, 0.470 kcal/mol, and 0.628 kcal/mol, respectively, according to DynaMut estimations. While other mutations (E484A, N501Y, K417N, Y505H, G496S) showed destabilizing results. The D614G, E484A, N501Y, K417N, Y505H, and G496S mutations increased the molecular flexibility of S-glycoprotein to interact with the ACE2 receptor, increasing the variant's infectivity. Our study will contribute to research on the SARS-CoV-2 variant, Omicron, by providing information on the mutational pattern and exciting properties of these eight significant mutations, such as antibody escape and infectivity quotient (stabilizing or destabilizing; increased or decreased molecular flexibility of S-glycoprotein to interact with the human ACE2 receptor).
Collapse
Affiliation(s)
- Chiranjib Chakraborty
- Department of Biotechnology, School of Life Science and Biotechnology, Adamas University, Kolkata, West Bengal 700126, India.
| | - Manojit Bhattacharya
- Department of Zoology, Fakir Mohan University, Vyasa Vihar, Balasore 756020, Odisha, India
| | - Ashish Ranjan Sharma
- Institute for Skeletal Aging & Orthopaedic Surgery, Hallym University-Chuncheon Sacred Heart Hospital, Chuncheon-si 24252, Gangwon-do, South Korea
| | - Bidyut Mallik
- Department of Applied Science, Galgotias College of Engineering and Technology, Knowledge Park-II, Greater Noida, Uttar Pradesh 201306, India
| |
Collapse
|
14
|
Laurents DV. AlphaFold 2 and NMR Spectroscopy: Partners to Understand Protein Structure, Dynamics and Function. Front Mol Biosci 2022; 9:906437. [PMID: 35655760 PMCID: PMC9152297 DOI: 10.3389/fmolb.2022.906437] [Citation(s) in RCA: 25] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 04/25/2022] [Indexed: 11/29/2022] Open
Abstract
The artificial intelligence program AlphaFold 2 is revolutionizing the field of protein structure determination as it accurately predicts the 3D structure of two thirds of the human proteome. Its predictions can be used directly as structural models or indirectly as aids for experimental structure determination using X-ray crystallography, CryoEM or NMR spectroscopy. Nevertheless, AlphaFold 2 can neither afford insight into how proteins fold, nor can it determine protein stability or dynamics. Rare folds or minor alternative conformations are also not predicted by AlphaFold 2 and the program does not forecast the impact of post translational modifications, mutations or ligand binding. The remaining third of human proteome which is poorly predicted largely corresponds to intrinsically disordered regions of proteins. Key to regulation and signaling networks, these disordered regions often form biomolecular condensates or amyloids. Fortunately, the limitations of AlphaFold 2 are largely complemented by NMR spectroscopy. This experimental approach provides information on protein folding and dynamics as well as biomolecular condensates and amyloids and their modulation by experimental conditions, small molecules, post translational modifications, mutations, flanking sequence, interactions with other proteins, RNA and virus. Together, NMR spectroscopy and AlphaFold 2 can collaborate to advance our comprehension of proteins.
Collapse
|
15
|
Zimmermann MT. Molecular Modeling is an Enabling Approach to Complement and Enhance Channelopathy Research. Compr Physiol 2022; 12:3141-3166. [PMID: 35578963 DOI: 10.1002/cphy.c190047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Hundreds of human membrane proteins form channels that transport necessary ions and compounds, including drugs and metabolites, yet details of their normal function or how function is altered by genetic variants to cause diseases are often unknown. Without this knowledge, researchers are less equipped to develop approaches to diagnose and treat channelopathies. High-resolution computational approaches such as molecular modeling enable researchers to investigate channelopathy protein function, facilitate detailed hypothesis generation, and produce data that is difficult to gather experimentally. Molecular modeling can be tailored to each physiologic context that a protein may act within, some of which may currently be difficult or impossible to assay experimentally. Because many genomic variants are observed in channelopathy proteins from high-throughput sequencing studies, methods with mechanistic value are needed to interpret their effects. The eminent field of structural bioinformatics integrates techniques from multiple disciplines including molecular modeling, computational chemistry, biophysics, and biochemistry, to develop mechanistic hypotheses and enhance the information available for understanding function. Molecular modeling and simulation access 3D and time-dependent information, not currently predictable from sequence. Thus, molecular modeling is valuable for increasing the resolution with which the natural function of protein channels can be investigated, and for interpreting how genomic variants alter them to produce physiologic changes that manifest as channelopathies. © 2022 American Physiological Society. Compr Physiol 12:3141-3166, 2022.
Collapse
Affiliation(s)
- Michael T Zimmermann
- Bioinformatics Research and Development Laboratory, Genomic Sciences and Precision Medicine Center, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Clinical and Translational Sciences Institute, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Department of Biochemistry, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| |
Collapse
|
16
|
BEHZADI PAYAM, GAJDÁCS MÁRIÓ. Worldwide Protein Data Bank (wwPDB): A virtual treasure for research in biotechnology. Eur J Microbiol Immunol (Bp) 2021; 11:77-86. [PMID: 34908533 PMCID: PMC8830413 DOI: 10.1556/1886.2021.00020] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2021] [Accepted: 11/23/2021] [Indexed: 12/25/2022] Open
Abstract
The Research Collaboratory for Structural Bioinformatics Protein Data Bank (RSCB PDB) provides a wide range of digital data regarding biology and biomedicine. This huge internet resource involves a wide range of important biological data, obtained from experiments around the globe by different scientists. The Worldwide Protein Data Bank (wwPDB) represents a brilliant collection of 3D structure data associated with important and vital biomolecules including nucleic acids (RNAs and DNAs) and proteins. Moreover, this database accumulates knowledge regarding function and evolution of biomacromolecules which supports different disciplines such as biotechnology. 3D structure, functional characteristics and phylogenetic properties of biomacromolecules give a deep understanding of the biomolecules' characteristics. An important advantage of the wwPDB database is the data updating time, which is done every week. This updating process helps users to have the newest data and information for their projects. The data and information in wwPDB can be a great support to have an accurate imagination and illustrations of the biomacromolecules in biotechnology. As demonstrated by the SARS-CoV-2 pandemic, rapidly reliable and accessible biological data for microbiology, immunology, vaccinology, and drug development are critical to address many healthcare-related challenges that are facing humanity. The aim of this paper is to introduce the readers to wwPDB, and to highlight the importance of this database in biotechnology, with the expectation that the number of scientists interested in the utilization of Protein Data Bank's resources will increase substantially in the coming years.
Collapse
Affiliation(s)
- PAYAM BEHZADI
- Department of Microbiology, College of Basic Sciences, Shahr-e-Qods Branch, Islamic Azad University, Tehran, 37541-374, Iran
| | - MÁRIÓ GAJDÁCS
- Department of Oral Biology and Experimental Dental Research, Faculty of Dentistry, University of Szeged, 6720, Szeged, Hungary,*Corresponding author. Tel.: +36-62-342-532. E-mail:
| |
Collapse
|
17
|
Modeling coronavirus spike protein dynamics: implications for immunogenicity and immune escape. Biophys J 2021; 120:5592-5618. [PMID: 34767789 PMCID: PMC8577870 DOI: 10.1016/j.bpj.2021.11.009] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2020] [Revised: 03/19/2021] [Accepted: 11/04/2021] [Indexed: 12/23/2022] Open
Abstract
The ongoing COVID-19 pandemic is a global public health emergency requiring urgent development of efficacious vaccines. While concentrated research efforts have focused primarily on antibody-based vaccines that neutralize SARS-CoV-2, and several first-generation vaccines have either been approved or received emergency use authorization, it is forecasted that COVID-19 will become an endemic disease requiring updated second-generation vaccines. The SARS-CoV-2 surface spike (S) glycoprotein represents a prime target for vaccine development because antibodies that block viral attachment and entry, i.e., neutralizing antibodies, bind almost exclusively to the receptor-binding domain. Here, we develop computational models for a large subset of S proteins associated with SARS-CoV-2, implemented through coarse-grained elastic network models and normal mode analysis. We then analyze local protein domain dynamics of the S protein systems and their thermal stability to characterize structural and dynamical variability among them. These results are compared against existing experimental data and used to elucidate the impact and mechanisms of SARS-CoV-2 S protein mutations and their associated antibody binding behavior. We construct a SARS-CoV-2 antigenic map and offer predictions about the neutralization capabilities of antibody and S mutant combinations based on protein dynamic signatures. We then compare SARS-CoV-2 S protein dynamics to SARS-CoV and MERS-CoV S proteins to investigate differing antibody binding and cellular fusion mechanisms that may explain the high transmissibility of SARS-CoV-2. The outbreaks associated with SARS-CoV, MERS-CoV, and SARS-CoV-2 over the last two decades suggest that the threat presented by coronaviruses is ever-changing and long term. Our results provide insights into the dynamics-driven mechanisms of immunogenicity associated with coronavirus S proteins and present a new, to our knowledge, approach to characterize and screen potential mutant candidates for immunogen design, as well as to characterize emerging natural variants that may escape vaccine-induced antibody responses.
Collapse
|
18
|
Ni D, Chai Z, Wang Y, Li M, Yu Z, Liu Y, Lu S, Zhang J. Along the allostery stream: Recent advances in computational methods for allosteric drug discovery. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1585] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Duan Ni
- College of Pharmacy Ningxia Medical University Yinchuan China
- The Charles Perkins Centre University of Sydney Sydney New South Wales Australia
| | - Zongtao Chai
- Department of Hepatic Surgery VI, Eastern Hepatobiliary Surgery Hospital Second Military Medical University Shanghai China
| | - Ying Wang
- State Key Laboratory of Oncogenes and Related Genes, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education Shanghai Jiao Tong University School of Medicine Shanghai China
| | - Mingyu Li
- State Key Laboratory of Oncogenes and Related Genes, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education Shanghai Jiao Tong University School of Medicine Shanghai China
| | | | - Yaqin Liu
- Medicinal Chemistry and Bioinformatics Center Shanghai Jiao Tong University School of Medicine Shanghai China
| | - Shaoyong Lu
- College of Pharmacy Ningxia Medical University Yinchuan China
- State Key Laboratory of Oncogenes and Related Genes, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education Shanghai Jiao Tong University School of Medicine Shanghai China
- Medicinal Chemistry and Bioinformatics Center Shanghai Jiao Tong University School of Medicine Shanghai China
| | - Jian Zhang
- College of Pharmacy Ningxia Medical University Yinchuan China
- State Key Laboratory of Oncogenes and Related Genes, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education Shanghai Jiao Tong University School of Medicine Shanghai China
- Medicinal Chemistry and Bioinformatics Center Shanghai Jiao Tong University School of Medicine Shanghai China
- School of Pharmaceutical Sciences Zhengzhou University Zhengzhou China
| |
Collapse
|
19
|
Kanapeckaitė A, Beaurivage C, Jančorienė L, Mažeikienė A. In silico drug discovery for a complex immunotherapeutic target - human c-Rel protein. Biophys Chem 2021; 276:106593. [PMID: 34087524 DOI: 10.1016/j.bpc.2021.106593] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Revised: 03/28/2021] [Accepted: 04/12/2021] [Indexed: 12/22/2022]
Abstract
Target evaluation and rational drug design rely on identifying and characterising small-molecule binding sites on therapeutically relevant target proteins. Immunotherapeutics development is especially challenging because of complex disease etiology and heterogenous nature of targets. c-Rel protein, a promising target in many human inflammatory and cancer pathologies, was selected as a case study for an effective in silico screening platform development since this transcription factor currently has no successful therapeutic inhibitors or modulators. This study introduces a novel in silico screening approach to probe binding sites using structural validation sets, molecular modelling and describes a method of a computer-aided drug design when a crystal structure is not available for the target of interest. In addition, we showed that binding sites can be analysed with the machine learning as well as molecular simulation approaches to help assess and systematically analyse how drug candidates can exert their mode of action. Finally, this cutting-edge approach was subjected to a high through-put virtual screen of selected 34 M drug-like compounds filtered from a library of 659 M compounds by identifying the most promising structures and proposing potential action mechanisms for the future development of highly selective human c-Rel inhibitors and/or modulators.
Collapse
Affiliation(s)
| | | | - Ligita Jančorienė
- Vilnius University Medical Faculty InsTtute of Clinical Medicine, Clinic of InfecTous Diseases and Dermatovenerology, Santariškių str. 14, 08406 Vilnius, Lithuania
| | - Asta Mažeikienė
- Department of Physiology, Biochemistry, Microbiology and Laboratory Medicine, Institute of Biomedical Sciences, Faculty of Medicine, Vilnius University, M. K. Čiurlionio g. 21, LT-03101, Vilnius, Lithuania
| |
Collapse
|
20
|
Tang QY, Kaneko K. Dynamics-Evolution Correspondence in Protein Structures. PHYSICAL REVIEW LETTERS 2021; 127:098103. [PMID: 34506164 DOI: 10.1103/physrevlett.127.098103] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Accepted: 07/28/2021] [Indexed: 06/13/2023]
Abstract
The genotype-phenotype mapping of proteins is a fundamental question in structural biology. In this Letter, with the analysis of a large dataset of proteins from hundreds of protein families, we quantitatively demonstrate the correlations between the noise-induced protein dynamics and mutation-induced variations of native structures, indicating the dynamics-evolution correspondence of proteins. Based on the investigations of the linear responses of native proteins, the origin of such a correspondence is elucidated. It is essential that the noise- and mutation-induced deformations of the proteins are restricted on a common low-dimensional subspace, as confirmed from the data. These results suggest an evolutionary mechanism of the proteins gaining both dynamical flexibility and evolutionary structural variability.
Collapse
Affiliation(s)
- Qian-Yuan Tang
- Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Komaba 3-8-1, Meguro-ku, Tokyo 153-8902, Japan
- Lab for Neural Computation and Adaptation, RIKEN Center for Brain Science, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Kunihiko Kaneko
- Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Komaba 3-8-1, Meguro-ku, Tokyo 153-8902, Japan
| |
Collapse
|
21
|
Teruel N, Mailhot O, Najmanovich RJ. Modelling conformational state dynamics and its role on infection for SARS-CoV-2 Spike protein variants. PLoS Comput Biol 2021; 17:e1009286. [PMID: 34351895 PMCID: PMC8384204 DOI: 10.1371/journal.pcbi.1009286] [Citation(s) in RCA: 54] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Revised: 08/24/2021] [Accepted: 07/17/2021] [Indexed: 01/21/2023] Open
Abstract
The SARS-CoV-2 Spike protein needs to be in an open-state conformation to interact with ACE2 to initiate viral entry. We utilise coarse-grained normal mode analysis to model the dynamics of Spike and calculate transition probabilities between states for 17081 variants including experimentally observed variants. Our results correctly model an increase in open-state occupancy for the more infectious D614G via an increase in flexibility of the closed-state and decrease of flexibility of the open-state. We predict the same effect for several mutations on glycine residues (404, 416, 504, 252) as well as residues K417, D467 and N501, including the N501Y mutation recently observed within the B.1.1.7, 501.V2 and P1 strains. This is, to our knowledge, the first use of normal mode analysis to model conformational state transitions and the effect of mutations on such transitions. The specific mutations of Spike identified here may guide future studies to increase our understanding of SARS-CoV-2 infection mechanisms and guide public health in their surveillance efforts. The present work explores the molecular mechanisms underlying and potentially helping new strains of SARS-CoV-2 to gain an evolutionary advantage during the ongoing COVID-19 pandemics. We show how a computational method called normal mode analysis that treats protein dynamics in a simplified manner is capable to predict the higher propensity of the Spike protein to be in the open state in which it is capable to interact with the human ACE2 receptor and thus facilitate cell entry. Because the simulation of the simplified computational model is relatively less demanding on resources than alternative methods, we were able to simulate over 17000 mutations in the SARS-CoV-2 Spike protein to identify multiple mutations that if they were to appear as the virus continues to evolve, could confer an evolutionary advantage. As a matter of fact, our predictions foresaw the emergence of particular mutations such as N501Y that appeared in several variants of concern. Our results can inform public health regarding new variants and serves as a proof of concept for the application of normal mode analysis to study the effect of mutations on both, protein dynamics and conformational transitions in a high-throughput manner.
Collapse
Affiliation(s)
- Natália Teruel
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Olivier Mailhot
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
- Institute for Research in Immunology and Cancer (IRIC), Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Rafael J. Najmanovich
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
- * E-mail:
| |
Collapse
|
22
|
Majumder S, Chaudhuri D, Datta J, Giri K. Exploring the intrinsic dynamics of SARS-CoV-2, SARS-CoV and MERS-CoV spike glycoprotein through normal mode analysis using anisotropic network model. J Mol Graph Model 2021; 102:107778. [PMID: 33099199 PMCID: PMC7567490 DOI: 10.1016/j.jmgm.2020.107778] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Revised: 10/05/2020] [Accepted: 10/13/2020] [Indexed: 12/19/2022]
Abstract
COVID-19 caused by SARS-CoV-2 have become a global pandemic with serious rate of fatalities. SARS-CoV and MERS-CoV have also caused serious outbreak previously but the intensity was much lower than the ongoing SARS-CoV-2. The main infectivity factor of all the three viruses is the spike glycoprotein. In this study we have examined the intrinsic dynamics of the prefusion, lying state of trimeric S protein of these viruses through Normal Mode Analysis using Anisotropic Network Model. The dynamic modes of the S proteins of the aforementioned viruses were compared by root mean square inner product (RMSIP), spectral overlap and cosine correlation matrix. S proteins show homogenous correlated or anticorrelated motions among their domains but direction of Cα atom among the spike proteins show less similarity. SARS-CoV-2 spike shows high vertically upward motion of the receptor binding motif implying its propensity for binding with the receptor even in the lying state. MERS-CoV spike shows unique dynamical motion compared to the other two S protein indicated by low RMSIP, spectral overlap and cosine correlation value. This study will guide in developing common potential inhibitor molecules against closed state of spike protein of these viruses to prevent conformational switching from lying to standing state.
Collapse
Affiliation(s)
| | | | - Joyeeta Datta
- Department of Life Sciences, Presidency University, Kolkata, India
| | - Kalyan Giri
- Department of Life Sciences, Presidency University, Kolkata, India.
| |
Collapse
|
23
|
Cirauqui Diaz N, Frezza E, Martin J. Using normal mode analysis on protein structural models. How far can we go on our predictions? Proteins 2020; 89:531-543. [PMID: 33349977 DOI: 10.1002/prot.26037] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2020] [Accepted: 12/12/2020] [Indexed: 01/01/2023]
Abstract
Normal mode analysis (NMA) is a fast and inexpensive approach that is largely used to gain insight into functional protein motions, and more recently to create conformations for further computational studies. However, when the protein structure is unknown, the use of computational models is necessary. Here, we analyze the capacity of NMA in internal coordinate space to predict protein motion, its intrinsic flexibility, and atomic displacements, using protein models instead of native structures, and the possibility to use it for model refinement. Our results show that NMA is quite insensitive to modeling errors, but that calculations are strictly reliable only for very accurate models. Our study also suggests that internal NMA is a more suitable tool for the improvement of structural models, and for integrating them with experimental data or in other computational techniques, such as protein docking or more refined molecular dynamics simulations.
Collapse
Affiliation(s)
- Nuria Cirauqui Diaz
- CNRS, UMR 5086 Molecular Microbiology and Structural Biochemistry, Université de Lyon, Lyon, France
| | - Elisa Frezza
- CiTCoM, CNRS, Université de Paris, Paris, France
| | - Juliette Martin
- CNRS, UMR 5086 Molecular Microbiology and Structural Biochemistry, Université de Lyon, Lyon, France
| |
Collapse
|
24
|
Expression and characterization of l-arabinose isomerase from Geobacillus stearothermophilus for improved activity under acidic condition. Protein Expr Purif 2020; 175:105692. [DOI: 10.1016/j.pep.2020.105692] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2019] [Revised: 04/11/2020] [Accepted: 06/24/2020] [Indexed: 11/21/2022]
|
25
|
Abstract
Iodothyronine deiodinases (Dios) are important selenoproteins that control the concentration of the active thyroid hormone (TH) triiodothyronine through regioselective deiodination. The X-ray structure of a truncated monomer of Type III Dio (Dio3), which deiodinates TH inner rings through a selenocysteine (Sec) residue, revealed a thioredoxin-fold catalytic domain supplemented with an unstructured Ω-loop. Loop dynamics are driven by interactions of the conserved Trp207 with solvent in multi-microsecond molecular dynamics simulations of the Dio3 thioredoxin(Trx)-fold domain. Hydrogen bonding interactions of Glu200 with residues conserved across the Dio family anchor the loop’s N-terminus to the active site Ser-Cys-Thr-Sec sequence. A key long-lived loop conformation coincides with the opening of a cryptic pocket that accommodates thyroxine (T4) through an I⋯Se halogen bond to Sec170 and the amino acid group with a polar cleft. The Dio3-T4 complex is stabilized by an I⋯O halogen bond between an outer ring iodine and Asp211, consistent with Dio3 selectivity for inner ring deiodination. Non-conservation of residues, such as Asp211, in other Dio types in the flexible portion of the loop sequence suggests a mechanism for regioselectivity through Dio type-specific loop conformations. Cys168 is proposed to attack the selenenyl iodide intermediate to regenerate Dio3 based upon structural comparison with related Trx-fold proteins.
Collapse
|
26
|
Grant BJ, Skjaerven L, Yao XQ. The Bio3D packages for structural bioinformatics. Protein Sci 2020; 30:20-30. [PMID: 32734663 DOI: 10.1002/pro.3923] [Citation(s) in RCA: 172] [Impact Index Per Article: 43.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2020] [Revised: 07/20/2020] [Accepted: 07/23/2020] [Indexed: 12/21/2022]
Abstract
Bio3D is a family of R packages for the analysis of biomolecular sequence, structure, and dynamics. Major functionality includes biomolecular database searching and retrieval, sequence and structure conservation analysis, ensemble normal mode analysis, protein structure and correlation network analysis, principal component, and related multivariate analysis methods. Here, we review recent package developments, including a new underlying segregation into separate packages for distinct analysis, and introduce a new method for structure analysis named ensemble difference distance matrix analysis (eDDM). The eDDM approach calculates and compares atomic distance matrices across large sets of homologous atomic structures to help identify the residue wise determinants underlying specific functional processes. An eDDM workflow is detailed along with an example application to a large protein family. As a new member of the Bio3D family, the Bio3D-eddm package supports both experimental and theoretical simulation-generated structures, is integrated with other methods for dissecting sequence-structure-function relationships, and can be used in a highly automated and reproducible manner. Bio3D is distributed as an integrated set of platform independent open source R packages available from: http://thegrantlab.org/bio3d/.
Collapse
Affiliation(s)
- Barry J Grant
- Division of Biological Sciences, Section of Molecular Biology, University of California, San Diego, La Jolla, California, USA
| | - Lars Skjaerven
- Division of Biological Sciences, Section of Molecular Biology, University of California, San Diego, La Jolla, California, USA
| | - Xin-Qiu Yao
- Department of Chemistry, Georgia State University, Atlanta, Georgia, USA
| |
Collapse
|
27
|
Levoin N, Jean M, Legembre P. CD95 Structure, Aggregation and Cell Signaling. Front Cell Dev Biol 2020; 8:314. [PMID: 32432115 PMCID: PMC7214685 DOI: 10.3389/fcell.2020.00314] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2020] [Accepted: 04/08/2020] [Indexed: 01/16/2023] Open
Abstract
CD95 is a pre-ligand-associated transmembrane (TM) receptor. The interaction with its ligand CD95L brings to a next level its aggregation and triggers different signaling pathways, leading to cell motility, differentiation or cell death. This diversity of biological responses associated with a unique receptor devoid of enzymatic property raises the question of whether different ligands exist, or whether the fine-tuned control of CD95 aggregation and conformation, its distribution within certain plasma membrane sub-domains or the pattern of post-translational modifications account for this such broad-range of cell signaling. Herein, we review how the different domains of CD95 and their post-translational modifications or the different forms of CD95L can participate in the receptor aggregation and induction of cell signaling. Understanding how CD95 response goes from cell death to cell proliferation, differentiation and motility is a prerequisite to reveal novel therapeutic options to treat chronic inflammatory disorders and cancers.
Collapse
Affiliation(s)
| | - Mickael Jean
- Univ Rennes, CNRS, ISCR-UMR 6226, Rennes, France
| | | |
Collapse
|
28
|
Tang QY, Kaneko K. Long-range correlation in protein dynamics: Confirmation by structural data and normal mode analysis. PLoS Comput Biol 2020; 16:e1007670. [PMID: 32053592 PMCID: PMC7043781 DOI: 10.1371/journal.pcbi.1007670] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Revised: 02/26/2020] [Accepted: 01/21/2020] [Indexed: 11/18/2022] Open
Abstract
Proteins in cellular environments are highly susceptible. Local perturbations to any residue can be sensed by other spatially distal residues in the protein molecule, showing long-range correlations in the native dynamics of proteins. The long-range correlations of proteins contribute to many biological processes such as allostery, catalysis, and transportation. Revealing the structural origin of such long-range correlations is of great significance in understanding the design principle of biologically functional proteins. In this work, based on a large set of globular proteins determined by X-ray crystallography, by conducting normal mode analysis with the elastic network models, we demonstrate that such long-range correlations are encoded in the native topology of the proteins. To understand how native topology defines the structure and the dynamics of the proteins, we conduct scaling analysis on the size dependence of the slowest vibration mode, average path length, and modularity. Our results quantitatively describe how native proteins balance between order and disorder, showing both dense packing and fractal topology. It is suggested that the balance between stability and flexibility acts as an evolutionary constraint for proteins at different sizes. Overall, our result not only gives a new perspective bridging the protein structure and its dynamics but also reveals a universal principle in the evolution of proteins at all different sizes. The long-range correlated fluctuations are closely related to many biological processes of the proteins, such as catalysis, ligand binding, biomolecular recognition, and transportation. In this paper, we elucidate the structural origin of the long-range correlation and describe how native contact topology defines the slow-mode dynamics of the native proteins. Our result suggests an evolutionary constraint for proteins at different sizes, which may shed light on solving many biophysical problems such as structure prediction, multi-scale molecular simulations, and the design of molecular machines. Moreover, in statistical physics, as the long-range correlations are notable signs of the critical point, unveiling the origin of such criticality can extend our understanding of the organizing principle of a large variety of complex systems.
Collapse
Affiliation(s)
- Qian-Yuan Tang
- Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Tokyo, Japan
- * E-mail:
| | - Kunihiko Kaneko
- Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Tokyo, Japan
| |
Collapse
|
29
|
Arai S, Shibazaki C, Shimizu R, Adachi M, Ishibashi M, Tokunaga H, Tokunaga M. Catalytic mechanism and evolutionary characteristics of thioredoxin from Halobacterium salinarum NRC-1. Acta Crystallogr D Struct Biol 2020; 76:73-84. [PMID: 31909745 DOI: 10.1107/s2059798319015894] [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: 08/01/2019] [Accepted: 11/25/2019] [Indexed: 01/08/2023] Open
Abstract
Thioredoxin (TRX) is an important antioxidant against oxidative stress. TRX from the extremely halophilic archaeon Halobacterium salinarum NRC-1 (HsTRX-A), which has the highest acidic residue content [(Asp + Glu)/(Arg + Lys + His) = 9.0] among known TRXs, was chosen to elucidate the catalytic mechanism and evolutionary characteristics associated with haloadaptation. X-ray crystallographic analysis revealed that the main-chain structure of HsTRX-A is similar to those of homologous TRXs; for example, the root-mean-square deviations on Cα atoms were <2.3 Å for extant archaeal TRXs and <1.5 Å for resurrected Precambrian TRXs. A unique water network was located near the active-site residues (Cys45 and Cys48) in HsTRX-A, which may enhance the proton transfer required for the reduction of substrates under a high-salt environment. The high density of negative charges on the molecular surface (3.6 × 10-3 e Å-2) should improve the solubility and haloadaptivity. Moreover, circular-dichroism measurements and enzymatic assays using a mutant HsTRX-A with deletion of the long flexible N-terminal region (Ala2-Pro17) revealed that Ala2-Pro17 improves the structural stability and the enzymatic activity of HsTRX-A under high-salt environments (>2 M NaCl). The elongation of the N-terminal region in HsTRX-A accompanies the increased hydrophilicity and acidic residue content but does not affect the structure of the active site. These observations offer insights into molecular evolution for haloadaptation and potential applications in halophilic protein-related biotechnology.
Collapse
Affiliation(s)
- Shigeki Arai
- Institute for Quantum Life Science, National Institutes for Quantum and Radiological Science and Technology, 2-4 Shirakata, Tokai, Ibaraki 319-1106, Japan
| | - Chie Shibazaki
- Institute for Quantum Life Science, National Institutes for Quantum and Radiological Science and Technology, 2-4 Shirakata, Tokai, Ibaraki 319-1106, Japan
| | - Rumi Shimizu
- Institute for Quantum Life Science, National Institutes for Quantum and Radiological Science and Technology, 2-4 Shirakata, Tokai, Ibaraki 319-1106, Japan
| | - Motoyasu Adachi
- Institute for Quantum Life Science, National Institutes for Quantum and Radiological Science and Technology, 2-4 Shirakata, Tokai, Ibaraki 319-1106, Japan
| | - Matsujiro Ishibashi
- Applied and Molecular Microbiology, Faculty of Agriculture, Kagoshima University, 1-21-24 Korimoto, Kagoshima 890-0065, Japan
| | - Hiroko Tokunaga
- Applied and Molecular Microbiology, Faculty of Agriculture, Kagoshima University, 1-21-24 Korimoto, Kagoshima 890-0065, Japan
| | - Masao Tokunaga
- Applied and Molecular Microbiology, Faculty of Agriculture, Kagoshima University, 1-21-24 Korimoto, Kagoshima 890-0065, Japan
| |
Collapse
|
30
|
Wako H, Endo S. Dynamic properties of oligomers that characterize low-frequency normal modes. Biophys Physicobiol 2019; 16:220-231. [PMID: 31984175 PMCID: PMC6976002 DOI: 10.2142/biophysico.16.0_220] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2019] [Accepted: 07/09/2019] [Indexed: 01/09/2023] Open
Abstract
Dynamics of oligomeric proteins (one trimer, two tetramers, and one hexamer) were studied by elastic network model-based normal mode analysis to characterize their large-scale concerted motions. First, the oligomer motions were simplified by considering rigid-body motions of individual subunits. The subunit motions were resolved into three components in a cylindrical coordinate system: radial, tangential, and axial ones. Single component is dominant in certain normal modes. However, more than one component is mixed in others. The subunits move symmetrically in certain normal modes and as a standing wave with several wave nodes in others. Secondly, special attention was paid to atoms on inter-subunit interfaces. Their displacement vectors were decomposed into intra-subunit deformative (internal) and rigid-body (external) motions in individual subunits. The fact that most of the cosines of the internal and external motion vectors were negative for the atoms on the inter-subunit interfaces, indicated their opposing movements. Finally, a structural network of residues defined for each normal mode was investigated; the network was constructed by connecting two residues in contact and moving coherently. The centrality measure “betweenness” of each residue was calculated for the networks. Several residues with significantly high betweenness were observed on the inter-subunit interfaces. The results indicate that these residues are responsible for oligomer dynamics. It was also observed that amino acid residues with significantly high betweenness were more conservative. This supports that the betweenness is an effective characteristic for identifying an important residue in protein dynamics.
Collapse
Affiliation(s)
- Hiroshi Wako
- School of Social Sciences, Waseda University, Shinjuku-ku, Tokyo 169-8050, Japan
| | - Shigeru Endo
- Department of Physics, School of Science, Kitasato University, Sagamihara, Kanagawa 252-0373, Japan
| |
Collapse
|
31
|
Yamato T, Laprévote O. Normal mode analysis and beyond. Biophys Physicobiol 2019; 16:322-327. [PMID: 31984187 PMCID: PMC6976091 DOI: 10.2142/biophysico.16.0_322] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2019] [Accepted: 08/02/2019] [Indexed: 01/05/2023] Open
Abstract
Normal mode analysis provides a powerful tool in biophysical computations. Particularly, we shed light on its application to protein properties because they directly lead to biological functions. As a result of normal mode analysis, the protein motion is represented as a linear combination of mutually independent normal mode vectors. It has been widely accepted that the large amplitude motions throughout the entire protein molecule can be well described with a few low-frequency normal modes. Furthermore, it is possible to represent the effect of external perturbations, e.g., ligand binding, hydrostatic pressure, as the shifts of normal mode variables. Making use of this advantage, we are able to explore mechanical properties of proteins such as Young's modulus and compressibility. Within thermally fluctuating protein molecules under physiological conditions, tightly packed amino acid residues interact with each other through heat and energy exchanges. Since the structure and dynamics of protein molecules are highly anisotropic, the flow of energy and heat should also be anisotropic. Based on the harmonic approximation of the heat current operator, it is possible to analyze the communication map of a protein molecule. By using this method, the energy transfer pathways of photoactive yellow protein were calculated. It turned out that these pathways are similar to those obtained via the Green-Kubo formalism with equilibrium molecular dynamics simulations, indicating that normal mode analysis captures the intrinsic nature of the transport properties of proteins.
Collapse
Affiliation(s)
- Takahisa Yamato
- Graduate School of Science, Nagoya University, Nagoya, Aichi 464-8602, Japan
| | - Olivier Laprévote
- Graduate School of Science, Nagoya University, Nagoya, Aichi 464-8602, Japan
- École supérieure de biotechnologie Strasbourg, 10413 – F-67412, Illkirsh France
| |
Collapse
|
32
|
Nojima H, Kiyota Y, Terashi G, Takeda-Shitaka M, Matsubara H. Geometrical Conversion of the EGFR Extracellular Domain by Adiabatic Mapping Combining Normal Mode Analysis of the Elastic Network Model and Energy Optimization. Chem Pharm Bull (Tokyo) 2019; 67:1061-1071. [PMID: 31582626 DOI: 10.1248/cpb.c19-00205] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
The activation of epidermal growth factor receptor (EGFR) involves the geometrical conversion of the extracellular domain (ECD) from the tethered to the extended forms with the dynamic rearrangement of the relative positions of four subdomains (SDs); however, this conversion process has not yet been thoroughly understood. We compare the two different forms of the X-ray crystal structures of ECD and simulate the ECD conversion process using adiabatic mapping that combines normal mode analysis of the elastic network model (ENM-NMA) and energy optimization. A comparison of the crystal structures reveals the rigidity of the intradomain geometry of the SD-I and -III backbone regardless of the form. The forward mapping from the tethered to the extended forms retains the intradomain geometry of the SD-I and -III backbone and reveals the trends to rearrange the relative positions of SD-I and -III and to dissociate the C-terminal tail of SD-IV from the hairpin loop in SD-II. The reverse mapping from the extended to the tethered forms complements the promotion of ECD conversion in the presence of epidermal growth factor (EGF).
Collapse
|
33
|
Aydınkal RM, Serçinoğlu O, Ozbek P. ProSNEx: a web-based application for exploration and analysis of protein structures using network formalism. Nucleic Acids Res 2019; 47:W471-W476. [PMID: 31114881 PMCID: PMC6602423 DOI: 10.1093/nar/gkz390] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2019] [Revised: 04/17/2019] [Accepted: 05/09/2019] [Indexed: 01/14/2023] Open
Abstract
ProSNEx (Protein Structure Network Explorer) is a web service for construction and analysis of Protein Structure Networks (PSNs) alongside amino acid flexibility, sequence conservation and annotation features. ProSNEx constructs a PSN by adding nodes to represent residues and edges between these nodes using user-specified interaction distance cutoffs for either carbon-alpha, carbon-beta or atom-pair contact networks. Different types of weighted networks can also be constructed by using either (i) the residue-residue interaction energies in the format returned by gRINN, resulting in a Protein Energy Network (PEN); (ii) the dynamical cross correlations from a coarse-grained Normal Mode Analysis (NMA) of the protein structure; (iii) interaction strength. Upon construction of the network, common network metrics (such as node centralities) as well as shortest paths between nodes and k-cliques are calculated. Moreover, additional features of each residue in the form of conservation scores and mutation/natural variant information are included in the analysis. By this way, tool offers an enhanced and direct comparison of network-based residue metrics with other types of biological information. ProSNEx is free and open to all users without login requirement at http://prosnex-tool.com.
Collapse
Affiliation(s)
- Rasim Murat Aydınkal
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
- Ali Nihat Gokyigit Foundation, Etiler, Istanbul 34340, Turkey
| | - Onur Serçinoğlu
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
- Department of Bioengineering, Faculty of Engineering, Recep Tayyip Erdoğan University, Rize 53100, Turkey
| | - Pemra Ozbek
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
| |
Collapse
|
34
|
Liang Z, Verkhivker GM, Hu G. Integration of network models and evolutionary analysis into high-throughput modeling of protein dynamics and allosteric regulation: theory, tools and applications. Brief Bioinform 2019; 21:815-835. [DOI: 10.1093/bib/bbz029] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2018] [Revised: 02/04/2019] [Accepted: 02/21/2019] [Indexed: 12/24/2022] Open
Abstract
Abstract
Proteins are dynamical entities that undergo a plethora of conformational changes, accomplishing their biological functions. Molecular dynamics simulation and normal mode analysis methods have become the gold standard for studying protein dynamics, analyzing molecular mechanism and allosteric regulation of biological systems. The enormous amount of the ensemble-based experimental and computational data on protein structure and dynamics has presented a major challenge for the high-throughput modeling of protein regulation and molecular mechanisms. In parallel, bioinformatics and systems biology approaches including genomic analysis, coevolution and network-based modeling have provided an array of powerful tools that complemented and enriched biophysical insights by enabling high-throughput analysis of biological data and dissection of global molecular signatures underlying mechanisms of protein function and interactions in the cellular environment. These developments have provided a powerful interdisciplinary framework for quantifying the relationships between protein dynamics and allosteric regulation, allowing for high-throughput modeling and engineering of molecular mechanisms. Here, we review fundamental advances in protein dynamics, network theory and coevolutionary analysis that have provided foundation for rapidly growing computational tools for modeling of allosteric regulation. We discuss recent developments in these interdisciplinary areas bridging computational biophysics and network biology, focusing on promising applications in allosteric regulations, including the investigation of allosteric communication pathways, protein–DNA/RNA interactions and disease mutations in genomic medicine. We conclude by formulating and discussing future directions and potential challenges facing quantitative computational investigations of allosteric regulatory mechanisms in protein systems.
Collapse
Affiliation(s)
- Zhongjie Liang
- School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Gennady M Verkhivker
- Department of Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA, USA
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA, USA
| | - Guang Hu
- School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| |
Collapse
|
35
|
Sorzano COS, Jiménez A, Mota J, Vilas JL, Maluenda D, Martínez M, Ramírez-Aportela E, Majtner T, Segura J, Sánchez-García R, Rancel Y, del Caño L, Conesa P, Melero R, Jonic S, Vargas J, Cazals F, Freyberg Z, Krieger J, Bahar I, Marabini R, Carazo JM. Survey of the analysis of continuous conformational variability of biological macromolecules by electron microscopy. Acta Crystallogr F Struct Biol Commun 2019; 75:19-32. [PMID: 30605122 PMCID: PMC6317454 DOI: 10.1107/s2053230x18015108] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2018] [Accepted: 10/26/2018] [Indexed: 11/10/2022] Open
Abstract
Single-particle analysis by electron microscopy is a well established technique for analyzing the three-dimensional structures of biological macromolecules. Besides its ability to produce high-resolution structures, it also provides insights into the dynamic behavior of the structures by elucidating their conformational variability. Here, the different image-processing methods currently available to study continuous conformational changes are reviewed.
Collapse
Affiliation(s)
| | - A. Jiménez
- National Center of Biotechnology (CSIC), Spain
| | - J. Mota
- National Center of Biotechnology (CSIC), Spain
| | - J. L. Vilas
- National Center of Biotechnology (CSIC), Spain
| | - D. Maluenda
- National Center of Biotechnology (CSIC), Spain
| | - M. Martínez
- National Center of Biotechnology (CSIC), Spain
| | | | - T. Majtner
- National Center of Biotechnology (CSIC), Spain
| | - J. Segura
- National Center of Biotechnology (CSIC), Spain
| | | | - Y. Rancel
- National Center of Biotechnology (CSIC), Spain
| | - L. del Caño
- National Center of Biotechnology (CSIC), Spain
| | - P. Conesa
- National Center of Biotechnology (CSIC), Spain
| | - R. Melero
- National Center of Biotechnology (CSIC), Spain
| | - S. Jonic
- Sorbonne Université, UMR CNRS 7590, Muséum National d’Histoire Naturelle, IRD, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, IMPMC, Paris, France
| | | | - F. Cazals
- Inria Sophia Antipolis – Méditerranée, France
| | | | | | | | | | | |
Collapse
|
36
|
Togashi Y, Flechsig H. Coarse-Grained Protein Dynamics Studies Using Elastic Network Models. Int J Mol Sci 2018; 19:ijms19123899. [PMID: 30563146 PMCID: PMC6320916 DOI: 10.3390/ijms19123899] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2018] [Revised: 11/28/2018] [Accepted: 12/03/2018] [Indexed: 01/03/2023] Open
Abstract
Elastic networks have been used as simple models of proteins to study their slow structural dynamics. They consist of point-like particles connected by linear Hookean springs and hence are convenient for linear normal mode analysis around a given reference structure. Furthermore, dynamic simulations using these models can provide new insights. As the computational cost associated with these models is considerably lower compared to that of all-atom models, they are also convenient for comparative studies between multiple protein structures. In this review, we introduce examples of coarse-grained molecular dynamics studies using elastic network models and their derivatives, focusing on the nonlinear phenomena, and discuss their applicability to large-scale macromolecular assemblies.
Collapse
Affiliation(s)
- Yuichi Togashi
- Research Center for the Mathematics on Chromatin Live Dynamics (RcMcD), Department of Mathematical and Life Sciences, Graduate School of Science, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Japan.
- RIKEN Center for Biosystems Dynamics Research (BDR), 6-2-3 Furuedai, Suita, Osaka 565-0874, Japan.
- Cybermedia Center, Osaka University, 5-1 Mihogaoka, Ibaraki, Osaka 567-0047, Japan.
| | - Holger Flechsig
- Nano Life Science Institute (WPI-NanoLSI), Kanazawa University, Kakuma-machi, Kanazawa, Ishikawa 920-1192, Japan.
| |
Collapse
|
37
|
Na H, Ben-Avraham D, Tirion MM. Slow normal modes of proteins are accurately reproduced across different platforms. Phys Biol 2018; 16:016003. [PMID: 30238928 DOI: 10.1088/1478-3975/aae333] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
The Protein data bank (PDB) (Berman et al 2000 Nucl. Acids Res. 28 235-42) contains the atomic structures of over 105 biomolecules with better than 2.8 Å resolution. The listing of the identities and coordinates of the atoms comprising each macromolecule permits an analysis of the slow-time vibrational response of these large systems to minor perturbations. 3D video animations of individual modes of oscillation demonstrate how regions interdigitate to create cohesive collective motions, providing a comprehensive framework for and familiarity with the overall 3D architecture. Furthermore, the isolation and representation of the softest, slowest deformation coordinates provide opportunities for the development of mechanical models of enzyme function. The eigenvector decomposition, therefore, must be accurate, reliable as well as rapid to be generally reported upon. We obtain the eigenmodes of a 1.2 Å 34 kDa PDB entry using either exclusively heavy atoms or partly or fully reduced atomic sets; Cartesian or internal coordinates; interatomic force fields derived either from a full Cartesian potential, a reduced atomic potential or a Gaussian distance-dependent potential; and independently developed software. These varied technologies are similar in that each maintains proper stereochemistry either by use of dihedral degrees of freedom which freezes bond lengths and bond angles, or by use of a full atomic potential that includes realistic bond length and angle restraints. We find that the shapes of the slowest eigenvectors are nearly identical, not merely similar.
Collapse
Affiliation(s)
- Hyuntae Na
- Computer Science, Penn State Harrisburg, Middletown, PA, United States of America
| | | | | |
Collapse
|
38
|
Laksmi FA, Arai S, Tsurumaru H, Nakamura Y, Saksono B, Tokunaga M, Ishibashi M. Improved substrate specificity for D-galactose of L-arabinose isomerase for industrial application. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2018; 1866:1084-1091. [PMID: 30282606 DOI: 10.1016/j.bbapap.2018.09.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/02/2018] [Revised: 08/22/2018] [Accepted: 09/10/2018] [Indexed: 10/28/2022]
Abstract
L-Arabinose isomerase isolated from Geobacillus stearothermophilus (GSAI) was modified to improve its substrate specificity for D-galactose for the production of D-tagatose, a potential reduced-energy sweetener. Among the selected residues, mutation at residue 18 produced a mutant strain, H18T, which exhibited increased activity for D-galactose compared with the wild-type (WT) enzyme. Analysis of the substrate specificity of H18T showed a 45.4% improvement for D-galactose. Replacing histidine with threonine at residue 18 resulted in approximately 2.7-fold and 1.8-fold higher substrate binding and catalytic efficiency, respectively, for D-galactose. Further enhancement of the specific activity and catalytic efficiency of H18T for D-galactose by up to 2.7-fold and 4.3-fold, respectively, was achieved by adding borate during L-arabinose isomerase catalysis. Moreover, H18T showed thermostability and no destabilization was detected, which is promising for the industrial production of D-tagatose.
Collapse
Affiliation(s)
- Fina Amreta Laksmi
- The United Graduate School of Agricultural Sciences, Kagoshima University, 1-21-24, Korimoto, Kagoshima 890-0065, Japan; Research Center for Biotechnology, Indonesian Institute of Sciences (LIPI), Jalan Raya Bogor Km. 46, Cibinong 16911, Indonesia
| | - Shigeki Arai
- National Institutes for Quantum and Radiological Science and Technology, 2-4 Shirakata, Tokai, Ibaraki 319-1106, Japan
| | - Hirohito Tsurumaru
- The United Graduate School of Agricultural Sciences, Kagoshima University, 1-21-24, Korimoto, Kagoshima 890-0065, Japan
| | - Yoshitaka Nakamura
- The United Graduate School of Agricultural Sciences, Kagoshima University, 1-21-24, Korimoto, Kagoshima 890-0065, Japan
| | - Budi Saksono
- Research Center for Biotechnology, Indonesian Institute of Sciences (LIPI), Jalan Raya Bogor Km. 46, Cibinong 16911, Indonesia
| | - Masao Tokunaga
- The United Graduate School of Agricultural Sciences, Kagoshima University, 1-21-24, Korimoto, Kagoshima 890-0065, Japan
| | - Matsujiro Ishibashi
- The United Graduate School of Agricultural Sciences, Kagoshima University, 1-21-24, Korimoto, Kagoshima 890-0065, Japan.
| |
Collapse
|
39
|
Bacterial flagellar switching: a molecular mechanism directed by the logic of an electric motor. J Mol Model 2018; 24:280. [DOI: 10.1007/s00894-018-3819-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2018] [Accepted: 08/30/2018] [Indexed: 11/27/2022]
|
40
|
Krüger A, Zimbres FM, Kronenberger T, Wrenger C. Molecular Modeling Applied to Nucleic Acid-Based Molecule Development. Biomolecules 2018; 8:E83. [PMID: 30150587 PMCID: PMC6163985 DOI: 10.3390/biom8030083] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2018] [Revised: 08/12/2018] [Accepted: 08/16/2018] [Indexed: 12/15/2022] Open
Abstract
Molecular modeling by means of docking and molecular dynamics (MD) has become an integral part of early drug discovery projects, enabling the screening and enrichment of large libraries of small molecules. In the past decades, special emphasis was drawn to nucleic acid (NA)-based molecules in the fields of therapy, diagnosis, and drug delivery. Research has increased dramatically with the advent of the SELEX (systematic evolution of ligands by exponential enrichment) technique, which results in single-stranded DNA or RNA sequences that bind with high affinity and specificity to their targets. Herein, we discuss the role and contribution of docking and MD to the development and optimization of new nucleic acid-based molecules. This review focuses on the different approaches currently available for molecular modeling applied to NA interaction with proteins. We discuss topics ranging from structure prediction to docking and MD, highlighting their main advantages and limitations and the influence of flexibility on their calculations.
Collapse
Affiliation(s)
- Arne Krüger
- Unit for Drug Discovery, Department of Parasitology, Institute of Biomedical Sciences, University of São Paulo, São Paulo, SP 05508-000, Brazil.
| | - Flávia M Zimbres
- Department of Biochemistry and Molecular Biology and Center for Tropical and Emerging Global Diseases, University of Georgia, Athens, GA 30602, USA.
| | - Thales Kronenberger
- Department of Internal Medicine VIII, University Hospital of Tübingen, 72076 Tübingen, Germany.
| | - Carsten Wrenger
- Unit for Drug Discovery, Department of Parasitology, Institute of Biomedical Sciences, University of São Paulo, São Paulo, SP 05508-000, Brazil.
| |
Collapse
|