1
|
Bayarsaikhan B, Zsidó BZ, Börzsei R, Hetényi C. Efficient Refinement of Complex Structures of Flexible Histone Peptides Using Post-Docking Molecular Dynamics Protocols. Int J Mol Sci 2024; 25:5945. [PMID: 38892133 PMCID: PMC11172440 DOI: 10.3390/ijms25115945] [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: 04/24/2024] [Revised: 05/26/2024] [Accepted: 05/27/2024] [Indexed: 06/21/2024] Open
Abstract
Histones are keys to many epigenetic events and their complexes have therapeutic and diagnostic importance. The determination of the structures of histone complexes is fundamental in the design of new drugs. Computational molecular docking is widely used for the prediction of target-ligand complexes. Large, linear peptides like the tail regions of histones are challenging ligands for docking due to their large conformational flexibility, extensive hydration, and weak interactions with the shallow binding pockets of their reader proteins. Thus, fast docking methods often fail to produce complex structures of such peptide ligands at a level appropriate for drug design. To address this challenge, and improve the structural quality of the docked complexes, post-docking refinement has been applied using various molecular dynamics (MD) approaches. However, a final consensus has not been reached on the desired MD refinement protocol. In this present study, MD refinement strategies were systematically explored on a set of problematic complexes of histone peptide ligands with relatively large errors in their docked geometries. Six protocols were compared that differ in their MD simulation parameters. In all cases, pre-MD hydration of the complex interface regions was applied to avoid the unwanted presence of empty cavities. The best-performing protocol achieved a median of 32% improvement over the docked structures in terms of the change in root mean squared deviations from the experimental references. The influence of structural factors and explicit hydration on the performance of post-docking MD refinements are also discussed to help with their implementation in future methods and applications.
Collapse
Affiliation(s)
- Bayartsetseg Bayarsaikhan
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, H-7624 Pécs, Hungary; (B.B.); (B.Z.Z.); (R.B.)
| | - Balázs Zoltán Zsidó
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, H-7624 Pécs, Hungary; (B.B.); (B.Z.Z.); (R.B.)
| | - Rita Börzsei
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, H-7624 Pécs, Hungary; (B.B.); (B.Z.Z.); (R.B.)
| | - Csaba Hetényi
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, H-7624 Pécs, Hungary; (B.B.); (B.Z.Z.); (R.B.)
- National Laboratory for Drug Research and Development, Magyar tudósok krt. 2, H-1117 Budapest, Hungary
| |
Collapse
|
2
|
Krokidis MG, Dimitrakopoulos GN, Vrahatis AG, Exarchos TP, Vlamos P. Challenges and limitations in computational prediction of protein misfolding in neurodegenerative diseases. Front Comput Neurosci 2024; 17:1323182. [PMID: 38250244 PMCID: PMC10796696 DOI: 10.3389/fncom.2023.1323182] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2023] [Accepted: 12/19/2023] [Indexed: 01/23/2024] Open
Affiliation(s)
| | | | | | | | - Panagiotis Vlamos
- Bioinformatics and Human Electrophysiology Laboratory, Department of Informatics, Ionian University, Corfu, Greece
| |
Collapse
|
3
|
King JE, Koes DR. Interpreting forces as deep learning gradients improves quality of predicted protein structures. Biophys J 2023:S0006-3495(23)04149-8. [PMID: 38104241 DOI: 10.1016/j.bpj.2023.12.011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2023] [Revised: 10/30/2023] [Accepted: 12/12/2023] [Indexed: 12/19/2023] Open
Abstract
Protein structure predictions from deep learning models like AlphaFold2, despite their remarkable accuracy, are likely insufficient for direct use in downstream tasks like molecular docking. The functionality of such models could be improved with a combination of increased accuracy and physical intuition. We propose a new method to train deep learning protein structure prediction models using molecular dynamics force fields to work toward these goals. Our custom PyTorch loss function, OpenMM-Loss, represents the potential energy of a predicted structure. OpenMM-Loss can be applied to any all-atom representation of a protein structure capable of mapping into our software package, SidechainNet. We demonstrate our method's efficacy by finetuning OpenFold. We show that subsequently predicted protein structures, both before and after a relaxation procedure, exhibit comparable accuracy while displaying lower potential energy and improved structural quality as assessed by MolProbity metrics.
Collapse
Affiliation(s)
- Jonathan Edward King
- Joint PhD Program in Computational Biology, Carnegie Mellon University-University of Pittsburgh, Pittsburgh, Pennsylvania
| | - David Ryan Koes
- Computational & Systems Biology, University of Pittsburgh, Pittsburgh, Pennsylvania.
| |
Collapse
|
4
|
Mi J, Zhou X, Sun R, Han J. Disabling spidroin N-terminal homologs' reverse reaction unveils why its intermolecular disulfide bonds have not evolved for 380 million years. Int J Biol Macromol 2023; 249:125974. [PMID: 37499718 DOI: 10.1016/j.ijbiomac.2023.125974] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Revised: 07/21/2023] [Accepted: 07/22/2023] [Indexed: 07/29/2023]
Abstract
Spiders, ubiquitous predators known for their powerful silks, rely on spidroins that self-assemble from high-concentration solutions stored in silk glands, which are mediated by the NT and CT domains. CT homodimers containing intermolecular disulfide bonds enhance silk performance, promoting spider survival and reproduction. However, no NT capable of forming such disulfide bonds has been identified. Our study reveals that NT homodimers with sulfur substitution can form under alkaline conditions, shedding light on why spiders have not evolved intermolecular disulfide bonds in the NT module during their 380 million years of evolution. This discovery significantly advances our comprehension of spider evolution and silk spinning mechanisms, while also providing novel insights into protein storage, assembly, as well as the mechanisms and therapeutic strategies for neurodegenerative diseases associated with protein aggregation.
Collapse
Affiliation(s)
- Junpeng Mi
- College of Biological Science and Medical Engineering, Donghua University, Shanghai 201620, China
| | - Xingping Zhou
- College of Biological Science and Medical Engineering, Donghua University, Shanghai 201620, China
| | - Rou Sun
- College of Biological Science and Medical Engineering, Donghua University, Shanghai 201620, China
| | - Jiaojiao Han
- Department of Clinical Hematology and osology, Shanghai center for clinical laboratory, Shanghai 200126, China.
| |
Collapse
|
5
|
Kagaya Y, Zhang Z, Ibtehaz N, Wang X, Nakamura T, Huang D, Kihara D. NuFold: A Novel Tertiary RNA Structure Prediction Method Using Deep Learning with Flexible Nucleobase Center Representation. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.09.20.558715. [PMID: 37790488 PMCID: PMC10542152 DOI: 10.1101/2023.09.20.558715] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
RNA is not only playing a core role in the central dogma as mRNA between DNA and protein, but also many non-coding RNAs have been discovered to have unique and diverse biological functions. As genome sequences become increasingly available and our knowledge of RNA sequences grows, the study of RNA's structure and function has become more demanding. However, experimental determination of three-dimensional RNA structures is both costly and time-consuming, resulting in a substantial disparity between RNA sequence data and structural insights. In response to this challenge, we propose a novel computational approach that harnesses state-of-the-art deep learning architecture NuFold to accurately predict RNA tertiary structures. This approach aims to offer a cost-effective and efficient means of bridging the gap between RNA sequence information and structural comprehension. NuFold implements a nucleobase center representation, which allows it to reproduce all possible nucleotide conformations accurately.
Collapse
Affiliation(s)
- Yuki Kagaya
- Department of Biological Sciences, Purdue University, West Lafayette, Indiana, 47907, USA
| | - Zicong Zhang
- Department of Computer Science, Purdue University, West Lafayette, Indiana, 47907, USA
| | - Nabil Ibtehaz
- Department of Computer Science, Purdue University, West Lafayette, Indiana, 47907, USA
| | - Xiao Wang
- Department of Computer Science, Purdue University, West Lafayette, Indiana, 47907, USA
| | - Tsukasa Nakamura
- Department of Biological Sciences, Purdue University, West Lafayette, Indiana, 47907, USA
| | - David Huang
- Department of Biological Sciences, Purdue University, West Lafayette, Indiana, 47907, USA
| | - Daisuke Kihara
- Department of Biological Sciences, Purdue University, West Lafayette, Indiana, 47907, USA
- Department of Computer Science, Purdue University, West Lafayette, Indiana, 47907, USA
| |
Collapse
|