1
|
J AR, D SP, Arumainathan S. Digital nets conformational sampling (DNCS) - an enhanced sampling technique to explore the conformational space of intrinsically disordered peptides. Phys Chem Chem Phys 2024. [PMID: 39158517 DOI: 10.1039/d4cp01891e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/20/2024]
Abstract
We propose digital nets conformational sampling (DNCS) - an enhanced sampling technique to explore the conformational ensembles of peptides, especially intrinsically disordered peptides (IDPs). The DNCS algorithm relies on generating history-dependent samples of dihedral variables using bitwise XOR operations and binary angle measurements (BAM). The algorithm was initially studied using met-enkephalin, a highly elusive neuropeptide. The DNCS method predicted near-native structures and the energy landscape of met-enkephalin was observed to be in direct correlation with earlier studies on the neuropeptide. Clustering analysis revealed that there are only 24 low-lying conformations of the molecule. The DNCS method has then been tested for predicting optimal conformations of 42 oligopeptides of length varying from 3 to 8 residues. The closest-to-native structures of 86% of cases are near-native and 24% of them have a root mean square deviation of less than 1.00 Å with respect to their crystal structures. The results obtained reveal that the DNCS method performs well, that too in less computational time.
Collapse
Affiliation(s)
- Abraham Rebairo J
- Department of Nuclear Physics, University of Madras, Chennai, Tamil Nadu, India.
| | - Sam Paul D
- Centre of Advanced Study in Crystallography and Biophysics, University of Madras, Chennai, Tamil Nadu, India
| | - Stephen Arumainathan
- Department of Nuclear Physics, University of Madras, Chennai, Tamil Nadu, India.
- Department of Materials Science, University of Madras, Chennai, Tamil Nadu, India
| |
Collapse
|
2
|
Jung J, Yagi K, Tan C, Oshima H, Mori T, Yu I, Matsunaga Y, Kobayashi C, Ito S, Ugarte La Torre D, Sugita Y. GENESIS 2.1: High-Performance Molecular Dynamics Software for Enhanced Sampling and Free-Energy Calculations for Atomistic, Coarse-Grained, and Quantum Mechanics/Molecular Mechanics Models. J Phys Chem B 2024; 128:6028-6048. [PMID: 38876465 PMCID: PMC11215777 DOI: 10.1021/acs.jpcb.4c02096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2024] [Revised: 05/15/2024] [Accepted: 05/21/2024] [Indexed: 06/16/2024]
Abstract
GENeralized-Ensemble SImulation System (GENESIS) is a molecular dynamics (MD) software developed to simulate the conformational dynamics of a single biomolecule, as well as molecular interactions in large biomolecular assemblies and between multiple biomolecules in cellular environments. To achieve the latter purpose, the earlier versions of GENESIS emphasized high performance in atomistic MD simulations on massively parallel supercomputers, with or without graphics processing units (GPUs). Here, we implemented multiscale MD simulations that include atomistic, coarse-grained, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. They demonstrate high performance and are integrated with enhanced conformational sampling algorithms and free-energy calculations without using external programs except for the QM programs. In this article, we review new functions, molecular models, and other essential features in GENESIS version 2.1 and discuss ongoing developments for future releases.
Collapse
Affiliation(s)
- Jaewoon Jung
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
| | - Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
| | - Cheng Tan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Hiraku Oshima
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
- Graduate
School of Life Science, University of Hyogo, Harima Science Park City, Hyogo 678-1297, Japan
| | - Takaharu Mori
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
- Department
of Chemistry, Tokyo University of Science, Shinjuku-ku, Tokyo 162-8601, Japan
| | - Isseki Yu
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
- Department
of Bioinformatics, Maebashi Institute of
Technology, Maebashi, Gunma 371-0816, Japan
| | - Yasuhiro Matsunaga
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
- Graduate
School of Science and Engineering, Saitama
University, Saitama 338-8570, Japan
| | - Chigusa Kobayashi
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Shingo Ito
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
| | - Diego Ugarte La Torre
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
3
|
Tolokh IS, Folescu DE, Onufriev AV. Inclusion of Water Multipoles into the Implicit Solvation Framework Leads to Accuracy Gains. J Phys Chem B 2024; 128:5855-5873. [PMID: 38860842 PMCID: PMC11194828 DOI: 10.1021/acs.jpcb.4c00254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2024] [Revised: 05/28/2024] [Accepted: 05/29/2024] [Indexed: 06/12/2024]
Abstract
The current practical "workhorses" of the atomistic implicit solvation─the Poisson-Boltzmann (PB) and generalized Born (GB) models─face fundamental accuracy limitations. Here, we propose a computationally efficient implicit solvation framework, the Implicit Water Multipole GB (IWM-GB) model, that systematically incorporates the effects of multipole moments of water molecules in the first hydration shell of a solute, beyond the dipole water polarization already present at the PB/GB level. The framework explicitly accounts for coupling between polar and nonpolar contributions to the total solvation energy, which is missing from many implicit solvation models. An implementation of the framework, utilizing the GAFF force field and AM1-BCC atomic partial charges model, is parametrized and tested against the experimental hydration free energies of small molecules from the FreeSolv database. The resulting accuracy on the test set (RMSE ∼ 0.9 kcal/mol) is 12% better than that of the explicit solvation (TIP3P) treatment, which is orders of magnitude slower. We also find that the coupling between polar and nonpolar parts of the solvation free energy is essential to ensuring that several features of the IWM-GB model are physically meaningful, including the sign of the nonpolar contributions.
Collapse
Affiliation(s)
- Igor S. Tolokh
- Department
of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Dan E. Folescu
- Department
of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department
of Mathematics, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Alexey V. Onufriev
- Department
of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department
of Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
- Center
for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
4
|
Ghosh D, Biswas A, Radhakrishna M. Advanced computational approaches to understand protein aggregation. BIOPHYSICS REVIEWS 2024; 5:021302. [PMID: 38681860 PMCID: PMC11045254 DOI: 10.1063/5.0180691] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 03/18/2024] [Indexed: 05/01/2024]
Abstract
Protein aggregation is a widespread phenomenon implicated in debilitating diseases like Alzheimer's, Parkinson's, and cataracts, presenting complex hurdles for the field of molecular biology. In this review, we explore the evolving realm of computational methods and bioinformatics tools that have revolutionized our comprehension of protein aggregation. Beginning with a discussion of the multifaceted challenges associated with understanding this process and emphasizing the critical need for precise predictive tools, we highlight how computational techniques have become indispensable for understanding protein aggregation. We focus on molecular simulations, notably molecular dynamics (MD) simulations, spanning from atomistic to coarse-grained levels, which have emerged as pivotal tools in unraveling the complex dynamics governing protein aggregation in diseases such as cataracts, Alzheimer's, and Parkinson's. MD simulations provide microscopic insights into protein interactions and the subtleties of aggregation pathways, with advanced techniques like replica exchange molecular dynamics, Metadynamics (MetaD), and umbrella sampling enhancing our understanding by probing intricate energy landscapes and transition states. We delve into specific applications of MD simulations, elucidating the chaperone mechanism underlying cataract formation using Markov state modeling and the intricate pathways and interactions driving the toxic aggregate formation in Alzheimer's and Parkinson's disease. Transitioning we highlight how computational techniques, including bioinformatics, sequence analysis, structural data, machine learning algorithms, and artificial intelligence have become indispensable for predicting protein aggregation propensity and locating aggregation-prone regions within protein sequences. Throughout our exploration, we underscore the symbiotic relationship between computational approaches and empirical data, which has paved the way for potential therapeutic strategies against protein aggregation-related diseases. In conclusion, this review offers a comprehensive overview of advanced computational methodologies and bioinformatics tools that have catalyzed breakthroughs in unraveling the molecular basis of protein aggregation, with significant implications for clinical interventions, standing at the intersection of computational biology and experimental research.
Collapse
Affiliation(s)
- Deepshikha Ghosh
- Department of Biological Sciences and Engineering, Indian Institute of Technology (IIT) Gandhinagar, Palaj, Gujarat 382355, India
| | - Anushka Biswas
- Department of Chemical Engineering, Indian Institute of Technology (IIT) Gandhinagar, Palaj, Gujarat 382355, India
| | | |
Collapse
|
5
|
Mavroeidi P, Zorba LP, Tzouras NV, Neofotistos SP, Georgiou N, Sahin K, Şentürk M, Durdagi S, Vougioukalakis GC, Mavromoustakos T. Are Terminal Alkynes Necessary for MAO-A/MAO-B Inhibition? A New Scaffold Is Revealed. Molecules 2024; 29:2486. [PMID: 38893361 PMCID: PMC11174103 DOI: 10.3390/molecules29112486] [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/03/2024] [Revised: 05/15/2024] [Accepted: 05/22/2024] [Indexed: 06/21/2024] Open
Abstract
A versatile family of quaternary propargylamines was synthesized employing the KA2 multicomponent reaction, through the single-step coupling of a number of amines, ketones, and terminal alkynes. Sustainable synthetic procedures using transition metal catalysts were employed in all cases. The inhibitory activity of these molecules was evaluated against human monoaminoxidase (hMAO)-A and hMAO-B enzymes and was found to be significant. The IC50 values for hMAO-B range from 152.1 to 164.7 nM while the IC50 values for hMAO-A range from 765.6 to 861.6 nM. Furthermore, these compounds comply with Lipinski's rule of five and exhibit no predicted toxicity. To understand their binding properties with the two target enzymes, key interactions were studied using molecular docking, all-atom molecular dynamics (MD) simulations, and MM/GBSA binding free energy calculations. Overall, herein, the reported family of propargylamines exhibits promise as potential treatments for neurodegenerative disorders, such as Parkinson's disease. Interestingly, this is the first time a propargylamine scaffold bearing an internal alkyne has been reported to show activity against monoaminoxidases.
Collapse
Affiliation(s)
- Panagiou Mavroeidi
- Laboratory of Organic Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, 15771 Athens, Greece; (P.M.); (L.P.Z.); (N.V.T.); (S.P.N.); (N.G.)
| | - Leandros P. Zorba
- Laboratory of Organic Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, 15771 Athens, Greece; (P.M.); (L.P.Z.); (N.V.T.); (S.P.N.); (N.G.)
| | - Nikolaos V. Tzouras
- Laboratory of Organic Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, 15771 Athens, Greece; (P.M.); (L.P.Z.); (N.V.T.); (S.P.N.); (N.G.)
| | - Stavros P. Neofotistos
- Laboratory of Organic Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, 15771 Athens, Greece; (P.M.); (L.P.Z.); (N.V.T.); (S.P.N.); (N.G.)
| | - Nikitas Georgiou
- Laboratory of Organic Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, 15771 Athens, Greece; (P.M.); (L.P.Z.); (N.V.T.); (S.P.N.); (N.G.)
| | - Kader Sahin
- Department of Analytical Chemistry, School of Pharmacy, Bahcesehir University, 34349 Istanbul, Turkey;
| | - Murat Şentürk
- Department of Biochemistry, Faculty of Pharmacy, Agri Ibrahim Cecen University, 04100 Agri, Turkey;
| | - Serdar Durdagi
- Molecular Therapy Laboratory, Department of Pharmaceutical Chemistry, School of Pharmacy, Bahcesehir University, 34349 Istanbul, Turkey
- Computational Biology and Molecular Simulations Laboratory, Department of Biophysics, School of Medicine, Bahcesehir University, 34349 Istanbul, Turkey
- Laboratory for Innovative Drugs (Lab4IND), Computational Drug Design Center (HİTMER), Bahcesehir University, 34349 Istanbul, Turkey
| | - Georgios C. Vougioukalakis
- Laboratory of Organic Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, 15771 Athens, Greece; (P.M.); (L.P.Z.); (N.V.T.); (S.P.N.); (N.G.)
| | - Thomas Mavromoustakos
- Laboratory of Organic Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, 15771 Athens, Greece; (P.M.); (L.P.Z.); (N.V.T.); (S.P.N.); (N.G.)
| |
Collapse
|
6
|
Rasul HO, Thomas NV, Ghafour DD, Aziz BK, Salgado M G, Mendoza-Huizar LH, Candia LG. Searching possible SARS-CoV-2 main protease inhibitors in constituents from herbal medicines using in silico studies. J Biomol Struct Dyn 2024; 42:4234-4248. [PMID: 37349945 DOI: 10.1080/07391102.2023.2220040] [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/01/2023] [Accepted: 05/23/2023] [Indexed: 06/24/2023]
Abstract
The largest threat to civilization since the Second World War is the spread of the new coronavirus disease (COVID-19). Therefore, there is an urgent need for innovative therapeutic medicines to treat COVID-19. Reusing bio-actives is a workable and efficient strategy in the battle against new epidemics because the process of developing new drugs is time-consuming. This research aimed to identify which herbal remedies had the highest affinity for the receptor and assess a variety of them for potential targets to suppress the SARS-CoV-2 Mpro. The use of AutoDock Vina for structure-based virtual screening was done first due to the importance of protein interactions in the development of drugs. Molecular docking was used in the comparative study to assess 89 different chemicals from medicinal herbs. To anticipate their effectiveness against the primary protease of SARS-CoV-2, more analysis was done on the ADMET profile, drug-likeness, and Lipinski's rule of five. The next step involved three replicas of 100 ns-long molecular dynamics simulations on the potential candidates, which were preceded by calculations of the binding free energy of MM-GBSA. The outcomes showed that Achyrodimer A, Cinchonain Ib, Symphonone F, and Lupeol acetate all performed well and had the highest 6LU7 binding affinities. Using RMSD, RMSF, and protein-ligand interactions, the stability of the protein-ligand complex was assessed. The studies indicate that bioactive substances obtained from herbal medicines may function as a COVID-19 therapeutic agent, necessitating additional wet lab research to confirm their therapeutic potential, efficacy, and pharmacological capacity against the condition.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Hezha O Rasul
- Department of Pharmaceutical Chemistry, College of Science, Charmo University, Chamchamal, Sulaimani, Iraq
| | - Noel Vinay Thomas
- Department of BioMedical Science, College of Science, Komar University of Science and Technology, Sulaimani, Iraq
| | - Dlzar D Ghafour
- Department of Medical Laboratory Science, College of Science, Komar University of Science and Technology, Sulaimani, Iraq
- Department of Chemistry, College of Science, University of Sulaimani, Sulaimani, Iraq
| | - Bakhtyar K Aziz
- Department of Nanoscience and Applied Chemistry, College of Science, Charmo University, Chamchamal, Sulaimani, Iraq
| | | | - L H Mendoza-Huizar
- Academic Area of Chemistry, Mineral de la Reforma, Autonomous University of Hidalgo State, Hidalgo, México
| | - Lorena Gerli Candia
- Departamento de Química Ambiental, Facultad de Ciencias, Universidad Católica de la Santísima Concepción, Concepción, Chile
| |
Collapse
|
7
|
Elhadi A, Zhao D, Ali N, Sun F, Zhong S. Multi-method computational evaluation of the inhibitors against leucine-rich repeat kinase 2 G2019S mutant for Parkinson's disease. Mol Divers 2024:10.1007/s11030-024-10808-w. [PMID: 38396210 DOI: 10.1007/s11030-024-10808-w] [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: 09/05/2023] [Accepted: 01/07/2024] [Indexed: 02/25/2024]
Abstract
Leucine-rich repeat kinase 2 G2019S mutant (LRRK2 G2019S) is a potential target for Parkinson's disease therapy. In this work, the computational evaluation of the LRRK2 G2019S inhibitors was conducted via a combined approach which contains a preliminary screening of a large database of compounds via similarity and pharmacophore, a secondary selection via structure-based affinity prediction and molecular docking, and a rescoring treatment for the final selection. MD simulations and MM/GBSA calculations were performed to check the agreement between different prediction methods for these inhibitors. 331 experimental ligands were collected, and 170 were used to build the structure-activity relationship. Eight representative ligand structural models were employed in similarity searching and pharmacophore screening over 14 million compounds. The process for selecting proper molecular descriptors provides a successful sample which can be used as a general strategy in QSAR modelling. The rescoring used in this work presents an alternative useful treatment for ranking and selection.
Collapse
Affiliation(s)
- Ahmed Elhadi
- School of Bioengineering, Dalian University of Technology, Dalian, 116024, Liaoning, People's Republic of China
| | - Dan Zhao
- School of Bioengineering, Dalian University of Technology, Dalian, 116024, Liaoning, People's Republic of China
| | - Noman Ali
- School of Bioengineering, Dalian University of Technology, Dalian, 116024, Liaoning, People's Republic of China
| | - Fusheng Sun
- School of Bioengineering, Dalian University of Technology, Dalian, 116024, Liaoning, People's Republic of China
| | - Shijun Zhong
- School of Bioengineering, Dalian University of Technology, Dalian, 116024, Liaoning, People's Republic of China.
| |
Collapse
|
8
|
Xiao T, Zhou Y, Li B. Energy-Scaled Debye-Hückel Theory for the Electrostatic Solvation Free Energy in Size-Asymmetric Electrolyte Solutions. J Phys Chem B 2024; 128:1029-1039. [PMID: 38235680 DOI: 10.1021/acs.jpcb.3c07233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
In this report, an energy-scaled Debye-Hückel theory is developed for fast and accurate evaluation of the electrostatic solvation free energy in size-asymmetric electrolyte solutions. A size-asymmetric electrolyte solution is mapped to a dielectric continuum medium with Debye-Hückel-like response. Based on the scaling relation of the electrostatic energy of a spherical ion in the small and large size limits, a Padé polynomial is used to interpolate the electrostatic energy at finite size. The Padé polynomial is further interpreted as the electrostatic energy of an effective Debye-Hückel mean field model, depicted by a modified Debye parameter and a surface charge density due to the size asymmetry of the solvent ions. This theory can distinguish the electrostatic energies and the electrostatic solvation free energies of solutes with the same size but opposite charges. Application to charged hard and charged soft spheres demonstrates the accuracy of our approach.
Collapse
Affiliation(s)
- Tiejun Xiao
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, Guiyang 550018, People's Republic of China
| | - Yun Zhou
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, Guiyang 550018, People's Republic of China
| | - Bo Li
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, Guiyang 550018, People's Republic of China
| |
Collapse
|
9
|
Liu A, Zhang H, Zheng Q, Wang S. The Potential of Cyclodextrins as Inhibitors for the BM2 Protein: An In Silico Investigation. Molecules 2024; 29:620. [PMID: 38338365 PMCID: PMC10856705 DOI: 10.3390/molecules29030620] [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: 12/19/2023] [Revised: 01/23/2024] [Accepted: 01/26/2024] [Indexed: 02/12/2024] Open
Abstract
The influenza BM2 transmembrane domain (BM2TM), an acid-activated proton channel, is an attractive antiviral target due to its essential roles during influenza virus replication, whereas no effective inhibitors have been reported for BM2. In this study, we draw inspiration from the properties of cyclodextrins (CDs) and hypothesize that CDs of appropriate sizes may possess the potential to act as inhibitors of the BM2TM proton channel. To explore this possibility, molecular dynamics simulations were employed to assess their inhibitory capabilities. Our findings reveal that CD4, CD5, and CD6 are capable of binding to the BM2TM proton channel, resulting in disrupted water networks and reduced hydrogen bond occupancy between H19 and the solvent within the BM2TM channel necessary for proton conduction. Notably, CD4 completely obstructs the BM2TM water channel. Based on these observations, we propose that CD4, CD5, and CD6 individually contribute to diminishing the proton transfer efficiency of the BM2 protein, and CD4 demonstrates promising potential as an inhibitor for the BM2 proton channel.
Collapse
Affiliation(s)
- Aijun Liu
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China; (A.L.); (H.Z.)
| | - Hao Zhang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China; (A.L.); (H.Z.)
| | - Qingchuan Zheng
- School of Pharmaceutical Sciences, Jilin University, Changchun, 130021, China
| | - Song Wang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China; (A.L.); (H.Z.)
| |
Collapse
|
10
|
Xiao T, Song X. A Gaussian field approach to the solvation of spherical ions in electrolyte solutions. J Chem Phys 2024; 160:034102. [PMID: 38226821 DOI: 10.1063/5.0187141] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2023] [Accepted: 12/25/2023] [Indexed: 01/17/2024] Open
Abstract
In this work, the electrostatic response of an electrolyte solution to a spherical ion is studied with a Gaussian field theory. In order to capture the ionic correlation effect in concentrated solutions, the bulk dielectric response function is described by a two-Yukawa response function. The modified response function of the solution is solved analytically in the spherical geometry, from which the induced charge density and the electrostatic energy are also derived analytically. Comparisons with results for small ions in electrolyte solutions from the hyper-netted chain theory demonstrate the validity of the Gaussian field theory.
Collapse
Affiliation(s)
- Tiejun Xiao
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, Guiyang People's Republic of China
| | - Xueyu Song
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| |
Collapse
|
11
|
Bass L, Elder LH, Folescu DE, Forouzesh N, Tolokh IS, Karpatne A, Onufriev AV. Improving the Accuracy of Physics-Based Hydration-Free Energy Predictions by Machine Learning the Remaining Error Relative to the Experiment. J Chem Theory Comput 2024; 20:396-410. [PMID: 38149593 PMCID: PMC10950260 DOI: 10.1021/acs.jctc.3c00981] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2023]
Abstract
The accuracy of computational models of water is key to atomistic simulations of biomolecules. We propose a computationally efficient way to improve the accuracy of the prediction of hydration-free energies (HFEs) of small molecules: the remaining errors of the physics-based models relative to the experiment are predicted and mitigated by machine learning (ML) as a postprocessing step. Specifically, the trained graph convolutional neural network attempts to identify the "blind spots" in the physics-based model predictions, where the complex physics of aqueous solvation is poorly accounted for, and partially corrects for them. The strategy is explored for five classical solvent models representing various accuracy/speed trade-offs, from the fast analytical generalized Born (GB) to the popular TIP3P explicit solvent model; experimental HFEs of small neutral molecules from the FreeSolv set are used for the training and testing. For all of the models, the ML correction reduces the resulting root-mean-square error relative to the experiment for HFEs of small molecules, without significant overfitting and with negligible computational overhead. For example, on the test set, the relative accuracy improvement is 47% for the fast analytical GB, making it, after the ML correction, almost as accurate as uncorrected TIP3P. For the TIP3P model, the accuracy improvement is about 39%, bringing the ML-corrected model's accuracy below the 1 kcal/mol threshold. In general, the relative benefit of the ML corrections is smaller for more accurate physics-based models, reaching the lower limit of about 20% relative accuracy gain compared with that of the physics-based treatment alone. The proposed strategy of using ML to learn the remaining error of physics-based models offers a distinct advantage over training ML alone directly on reference HFEs: it preserves the correct overall trend, even well outside of the training set.
Collapse
Affiliation(s)
- Lewis Bass
- Department of Computer Engineering, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Luke H Elder
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Dan E Folescu
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department of Mathematics, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Negin Forouzesh
- Department of Computer Science, California State University, Los Angeles, California 90032, United States
| | - Igor S Tolokh
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Anuj Karpatne
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Alexey V Onufriev
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department of Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
- Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
12
|
Tettevi EJ, Kuevi DNO, Sumabe BK, Simpong DL, Maina MB, Dongdem JT, Osei-Atweneboana MY, Ocloo A. In Silico Identification of a Potential TNF-Alpha Binder Using a Structural Similarity: A Potential Drug Repurposing Approach to the Management of Alzheimer's Disease. BIOMED RESEARCH INTERNATIONAL 2024; 2024:9985719. [PMID: 38221912 PMCID: PMC10787656 DOI: 10.1155/2024/9985719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/02/2023] [Revised: 11/25/2023] [Accepted: 12/13/2023] [Indexed: 01/16/2024]
Abstract
Introduction Alzheimer's disease (AD) is a neurodegenerative disorder with no conclusive remedy. Yohimbine, found in Rauwolfia vomitoria, may reduce brain inflammation by targeting tumour necrosis factor-alpha (TNFα), implicated in AD pathogenesis. Metoserpate, a synthetic compound, may inhibit TNFα. The study is aimed at assessing the potential utility of repurposing metoserpate for TNFα inhibition to reduce neuronal damage and inflammation in AD. The development of safe and effective treatments for AD is crucial to address the growing burden of the disease, which is projected to double over the next two decades. Methods Our study repurposed an FDA-approved drug as TNFα inhibitor for AD management using structural similarity studies, molecular docking, and molecular dynamics simulations. Yohimbine was used as a reference compound. Molecular docking used SeeSAR, and molecular dynamics simulation used GROMACS. Results Metoserpate was selected from 10 compounds similar to yohimbine based on pharmacokinetic properties and FDA approval status. Molecular docking and simulation studies showed a stable interaction between metoserpate and TNFα over 100 ns (100000 ps). This suggests a reliable and robust interaction between the protein and ligand, supporting the potential utility of repurposing metoserpate for TNFα inhibition in AD treatment. Conclusion Our study has identified metoserpate, a previously FDA-approved antihypertensive agent, as a promising candidate for inhibiting TNFα in the management of AD.
Collapse
Affiliation(s)
- Edward Jenner Tettevi
- Department of Biochemistry, Cell and Molecular Biology, School of Biological Science, University of Ghana, Legon, Accra, P.O. Box LG 25, Ghana
- West African Centre for Cell Biology of Infectious Pathogens, School of Biological Science, University of Ghana, Legon, Accra, P.O. Box LG 25, Ghana
- Biomedical and Public Health Research Unit, Council for Scientific and Industrial Research-Water Research Institute, Accra, P.O. Box M 32, Ghana
| | - Deryl Nii Okantey Kuevi
- Biomedical and Public Health Research Unit, Council for Scientific and Industrial Research-Water Research Institute, Accra, P.O. Box M 32, Ghana
| | - Balagra Kasim Sumabe
- Biomedical and Public Health Research Unit, Council for Scientific and Industrial Research-Water Research Institute, Accra, P.O. Box M 32, Ghana
| | - David Larbi Simpong
- Department of Medical Laboratory Sciences, College of Health and Allied Sciences, University of Cape Coast, Cape Coast, Ghana
| | - Mahmoud B. Maina
- Serpell Laboratory, Sussex Neuroscience, School of Life Sciences, University of Sussex, UK
- Biomedical Science Research and Training Centre, College of Medical Sciences, Yobe State University, Damaturu, Nigeria
| | - Julius T. Dongdem
- Department of Biochemistry and Molecular Medicine, School of Medicine, University for Development Studies, Tamale Campus, Ghana
| | - Mike Y. Osei-Atweneboana
- Biomedical and Public Health Research Unit, Council for Scientific and Industrial Research-Water Research Institute, Accra, P.O. Box M 32, Ghana
- CSIR-College of Science and Technology, 2nd CSIR Close, Airport Residential Area, Behind Golden Tulip Hotel, Greater Accra Region, Ghana
| | - Augustine Ocloo
- Department of Biochemistry, Cell and Molecular Biology, School of Biological Science, University of Ghana, Legon, Accra, P.O. Box LG 25, Ghana
| |
Collapse
|
13
|
Odunitan TT, Saibu OA, Apanisile BT, Omoboyowa DA, Balogun TA, Awe AV, Ajayi TM, Olagunju GV, Mahmoud FM, Akinboade M, Adeniji CB, Abdulazeez WO. Integrating biocomputational techniques for Breast cancer drug discovery via the HER-2, BCRA, VEGF and ER protein targets. Comput Biol Med 2024; 168:107737. [PMID: 38000249 DOI: 10.1016/j.compbiomed.2023.107737] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2023] [Revised: 11/03/2023] [Accepted: 11/16/2023] [Indexed: 11/26/2023]
Abstract
Computational modelling remains an indispensable technique in drug discovery. With myriad of high computing resources, and improved modelling algorithms, there has been a high-speed in the drug development cycle with promising success rate compared to the traditional route. For example, lapatinib; a well-known anticancer drug with clinical applications was discovered with computational drug design techniques. Similarly, molecular modelling has been applied to various disease areas ranging from cancer to neurodegenerative diseases. The techniques ranges from high-throughput virtual screening, molecular mechanics with generalized Born and surface area solvation (MM/GBSA) to molecular dynamics simulation. This review focuses on the application of computational modelling tools in the identification of drug candidates for Breast cancer. First, we begin with a succinct overview of molecular modelling in the drug discovery process. Next, we take note of special efforts on the developments and applications of combining these techniques with particular emphasis on possible breast cancer therapeutic targets such as estrogen receptor (ER), human epidermal growth factor receptor 2 (HER2), vascular endothelial growth factor (VEGF), breast cancer gene 1 (BRCA1), and breast cancer gene 2 (BRCA2). Finally, we discussed the search for covalent inhibitors against these receptors using computational techniques, advances, pitfalls, possible solutions, and future perspectives.
Collapse
Affiliation(s)
- Tope T Odunitan
- Department of Biochemistry, Ladoke Akintola University of Technology, Ogbomoso, Oyo State, Nigeria; Genomics Unit, Helix Biogen Institute, Ogbomoso, Oyo State, Nigeria
| | - Oluwatosin A Saibu
- Department of Chemistry and Biochemistry, New Mexico State University, Las Cruces, NM, USA.
| | - Boluwatife T Apanisile
- Department of Nutrition and Dietetics, Ladoke Akintola University of Technology, Ogbomoso, Oyo State, Nigeria
| | - Damilola A Omoboyowa
- Department of Biochemistry, Adekunle Ajasin University, Akungba-Akoko, Oyo State, Nigeria
| | - Toheeb A Balogun
- Department of Biological Sciences, University of California, San Diego, CA, USA
| | - Adeyoola V Awe
- Department of Medical Laboratory Science, Lead City, University, Ibadan, Oyo State, Nigeria
| | - Temitope M Ajayi
- Department of Biochemistry, Ladoke Akintola University of Technology, Ogbomoso, Oyo State, Nigeria
| | - Grace V Olagunju
- Department of Molecular Biology, New Mexico State University, Las Cruces, NM, USA
| | - Fatimah M Mahmoud
- Department of Molecular Biology, New Mexico State University, Las Cruces, NM, USA
| | - Modinat Akinboade
- Department of Biochemistry, Ladoke Akintola University of Technology, Ogbomoso, Oyo State, Nigeria
| | - Catherine B Adeniji
- Department of Environmental Management and Toxicology, Lead City University, Ibadan, Oyo State, Nigeria
| | - Waliu O Abdulazeez
- Department of Biochemistry, Ladoke Akintola University of Technology, Ogbomoso, Oyo State, Nigeria
| |
Collapse
|
14
|
Fu Y, Christov PP, Kingsley PJ, Richie-Jannetta RM, Marnett LJ, Stone MP. Base-Displaced Intercalated Structure of the 3-(2-Deoxy-β-D-erythropentofuranosyl)-pyrimido[1,2- f]purine-6,10(3 H,5 H)-dione (6-oxo-M 1dG) Lesion in DNA. Chem Res Toxicol 2023; 36:1947-1960. [PMID: 37989274 PMCID: PMC10731638 DOI: 10.1021/acs.chemrestox.3c00226] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Revised: 09/21/2023] [Accepted: 09/25/2023] [Indexed: 11/23/2023]
Abstract
The genotoxic 3-(2-deoxy-β-D-erythro-pentofuranosyl)pyrimido[1,2-α]purin-10(3H)-one (M1dG) DNA lesion arises from endogenous exposures to base propenals generated by oxidative damage and from exposures to malondialdehyde (MDA), produced by lipid peroxidation. Once formed, M1dG may oxidize, in vivo, to 3-(2-deoxy-β-D-erythropentofuranosyl)-pyrimido[1,2-f]purine-6,10(3H,5H)-dione (6-oxo-M1dG). The latter blocks DNA replication and is a substrate for error-prone mutagenic bypass by the Y-family DNA polymerase hpol η. To examine structural consequences of 6-oxo-M1dG damage in DNA, we conducted NMR studies of 6-oxo-M1dG incorporated site-specifically into 5' -d(C1A2T3X4A5T6G7A8C9G10C11T12)-3':5'-d(A13G14C15G16T17C18A19T20C21A22T23G24)-3' (X = 6-oxo-M1dG). NMR spectra afforded detailed resonance assignments. Chemical shift analyses revealed that nucleobase C21, complementary to 6-oxo-M1dG, was deshielded compared with the unmodified duplex. Sequential NOEs between 6-oxo-M1dG and A5 were disrupted, as well as NOEs between T20 and C21 in the complementary strand. The structure of the 6-oxo-M1dG modified DNA duplex was refined by using molecular dynamics (rMD) calculations restrained by NOE data. It revealed that 6-oxo-M1dG intercalated into the duplex and remained in the anti-conformation about the glycosyl bond. The complementary cytosine C21 extruded into the major groove, accommodating the intercalated 6-oxo-M1dG. The 6-oxo-M1dG H7 and H8 protons faced toward the major groove, while the 6-oxo-M1dG imidazole proton H2 faced into the major groove. Structural perturbations to dsDNA were limited to the 6-oxo-M1dG damaged base pair and the flanking T3:A22 and A5:T20 base pairs. Both neighboring base pairs remained within the Watson-Crick hydrogen bonding contact. The 6-oxo-M1dG did not stack well with the 5'-neighboring base pair T3:A22 but showed improved stacking with the 3'-neighboring base pair A5:T20. Overall, the base-displaced intercalated structure was consistent with thermal destabilization of the 6-oxo-M1dG damaged DNA duplex; thermal melting temperature data showed a 15 °C decrease in Tm compared to the unmodified duplex. The structural consequences of 6-oxo-M1dG formation in DNA are evaluated in the context of the chemical biology of this lesion.
Collapse
Affiliation(s)
- Yizhi Fu
- Departments of Chemistry and Biochemistry,
and the Vanderbilt-Ingram Cancer Center, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Plamen P. Christov
- Departments of Chemistry and Biochemistry,
and the Vanderbilt-Ingram Cancer Center, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Philip J. Kingsley
- Departments of Chemistry and Biochemistry,
and the Vanderbilt-Ingram Cancer Center, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Robyn M. Richie-Jannetta
- Departments of Chemistry and Biochemistry,
and the Vanderbilt-Ingram Cancer Center, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Lawrence J. Marnett
- Departments of Chemistry and Biochemistry,
and the Vanderbilt-Ingram Cancer Center, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Michael P. Stone
- Departments of Chemistry and Biochemistry,
and the Vanderbilt-Ingram Cancer Center, Vanderbilt University, Nashville, Tennessee 37235, United States
| |
Collapse
|
15
|
Aoki E, Manabe N, Ohno S, Aoki T, Furukawa JI, Togayachi A, Aoki-Kinoshita K, Inokuchi JI, Kurosawa K, Kaname T, Yamaguchi Y, Nishihara S. Predicting the pathogenicity of missense variants based on protein instability to support diagnosis of patients with novel variants of ARSL. Mol Genet Metab Rep 2023; 37:101016. [PMID: 38053926 PMCID: PMC10694752 DOI: 10.1016/j.ymgmr.2023.101016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Revised: 10/20/2023] [Accepted: 10/21/2023] [Indexed: 12/07/2023] Open
Abstract
Rare diseases are estimated to affect 3.5%-5.9% of the population worldwide and are difficult to diagnose. Genome analysis is useful for diagnosis. However, since some variants, especially missense variants, are also difficult to interpret, tools to accurately predict the effect of missense variants are very important and needed. Here we developed a method, "VarMeter", to predict whether a missense variant is damaging based on Gibbs free energy and solvent-accessible surface area calculated from the AlphaFold 3D protein model. We applied this method to the whole-exome sequencing data of 900 individuals with rare or undiagnosed disease in our in-house database, and identified four who were hemizygous for missense variants of arylsulfatase L (ARSL; known as the genetic cause of chondrodysplasia punctata 1, CPDX1). Two individuals had a novel Ser89 to Asn (Ser89Asn) or Arg469 to Trp (Arg469Trp) substitution, respectively predicted as "damaging" or "benign"; the other two had an Arg111 to His (Arg111His) or Gly117 to Arg (Gly117Arg) substitution, respectively predicted as "damaging" or "possibly damaging" and previously reported in patients showing clinical manifestations of CDPX1. Expression and analysis of the missense variant proteins showed that the predicted pathogenic variants (Ser89Asn, Arg111His, and Gly117Arg) had complete loss of sulfatase activity and reduced protease resistance due to destabilization of protein structure, while the predicted benign variant (Arg469Trp) had activity and protease resistance comparable to those of wild-type ARSL. The individual with the novel pathogenic Ser89Asn variant exhibited characteristics of CDPX1, while the individual with the benign Arg469Trp variant exhibited no such characteristics. These findings demonstrate that VarMeter may be used to predict the deleteriousness of variants found in genome sequencing data and thereby support disease diagnosis.
Collapse
Affiliation(s)
- Eriko Aoki
- Glycan & Life Systems Integration Center (GaLSIC), Soka University, Hachioji 192-8577, Japan
| | - Noriyoshi Manabe
- Division of Structural Glycobiology, Institute of Molecular Biomembrane and Glycobiology, Tohoku Medical and Pharmaceutical University, Sendai 981-8558, Japan
| | - Shiho Ohno
- Division of Structural Glycobiology, Institute of Molecular Biomembrane and Glycobiology, Tohoku Medical and Pharmaceutical University, Sendai 981-8558, Japan
| | - Taiga Aoki
- Department of Genome Medicine, National Center for Child Health and Development, Tokyo 157-8535, Japan
| | - Jun-Ichi Furukawa
- Institute for Glyco-Core Research (iGCORE), Nagoya University, Nagoya 466-8550, Japan
| | - Akira Togayachi
- Glycan & Life Systems Integration Center (GaLSIC), Soka University, Hachioji 192-8577, Japan
| | - Kiyoko Aoki-Kinoshita
- Glycan & Life Systems Integration Center (GaLSIC), Soka University, Hachioji 192-8577, Japan
| | - Jin-Ichi Inokuchi
- Division of Structural Glycobiology, Institute of Molecular Biomembrane and Glycobiology, Tohoku Medical and Pharmaceutical University, Sendai 981-8558, Japan
| | - Kenji Kurosawa
- Division of Medical Genetics, Kanagawa Children's Medical Center, Yokohama 232-8555, Japan
| | - Tadashi Kaname
- Department of Genome Medicine, National Center for Child Health and Development, Tokyo 157-8535, Japan
| | - Yoshiki Yamaguchi
- Division of Structural Glycobiology, Institute of Molecular Biomembrane and Glycobiology, Tohoku Medical and Pharmaceutical University, Sendai 981-8558, Japan
| | - Shoko Nishihara
- Glycan & Life Systems Integration Center (GaLSIC), Soka University, Hachioji 192-8577, Japan
| |
Collapse
|
16
|
Awan B, Khan MA, Ahmad I, Masood A, Raza A, Khaliq S, Ullah F, Ahmed J, Khan MR. Norfloxacin derivatives as DNA gyrase and urease inhibitors: synthesis, biological evaluation and molecular docking. Future Med Chem 2023; 15:2181-2194. [PMID: 37997685 DOI: 10.4155/fmc-2023-0225] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Accepted: 10/16/2023] [Indexed: 11/25/2023] Open
Abstract
Background: DNA gyrase and urease enzymes are important targets for the treatment of gastroenteritis, appendicitis, tuberculosis, urinary tract infections and Crohn's disease. Materials & methods: Esterification of norfloxacin was performed to enhance DNA gyrase and urease enzyme inhibition potential. Structure elucidation and chemical characterization were done through spectral (1H NMR, Fourier transform IR, 13C NMR) and carbon, hydrogen, nitrogen and sulfur analysis along with molecular docking. Results & conclusion: The majority of derivatives exhibited significant results but the 3e derivative showed maximum bactericidal, DPPH scavenging (96%), DNA gyrase and urease enzyme inhibitory activity with IC50 of 0.15 ± 0.24 and 1.14 ± 0.11 μM respectively which was further supported by molecular docking studies. So, the active derivatives can serve as a lead compound for the treatment of various pathological conditions.
Collapse
Affiliation(s)
- Breena Awan
- Department of Pharmaceutical chemistry, Faculty of Pharmacy, The Islamia University of Bahawalpur, Bahawalpur, 63100, Pakistan
| | - Mohsin Abbas Khan
- Department of Pharmaceutical chemistry, Faculty of Pharmacy, The Islamia University of Bahawalpur, Bahawalpur, 63100, Pakistan
- Institute of Pharmaceutical Science, Faculty of Life Science & Medicine, King's College, London, SE1 9NH, UK
| | - Irshad Ahmad
- Department of Pharmaceutical chemistry, Faculty of Pharmacy, The Islamia University of Bahawalpur, Bahawalpur, 63100, Pakistan
| | - Anum Masood
- Department of Pharmaceutical chemistry, Faculty of Pharmacy, The Islamia University of Bahawalpur, Bahawalpur, 63100, Pakistan
| | - Asim Raza
- Department of Pharmaceutical chemistry, Faculty of Pharmacy, The Islamia University of Bahawalpur, Bahawalpur, 63100, Pakistan
| | - Saharish Khaliq
- Department of Pharmaceutical chemistry, Faculty of Pharmacy, The Islamia University of Bahawalpur, Bahawalpur, 63100, Pakistan
| | - Farhat Ullah
- Department of Pharmaceutical chemistry, Faculty of Pharmacy, The Islamia University of Bahawalpur, Bahawalpur, 63100, Pakistan
| | - Javed Ahmed
- Department of Pharmaceutical chemistry, Faculty of Pharmacy, The Islamia University of Bahawalpur, Bahawalpur, 63100, Pakistan
| | - Muhammad Rizwan Khan
- Department of Pharmaceutical chemistry, Faculty of Pharmacy, The Islamia University of Bahawalpur, Bahawalpur, 63100, Pakistan
| |
Collapse
|
17
|
Gallardo A, Dutagaci B. Binding of small molecule inhibitors to RNA polymerase-Spt5 complex impacts RNA and DNA stability. J Comput Aided Mol Des 2023; 38:1. [PMID: 37987925 PMCID: PMC10663202 DOI: 10.1007/s10822-023-00543-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Accepted: 11/08/2023] [Indexed: 11/22/2023]
Abstract
Spt5 is an elongation factor that associates with RNA polymerase II (Pol II) during transcription and has important functions in promoter-proximal pausing and elongation processivity. Spt5 was also recognized for its roles in the transcription of expanded-repeat genes that are related to neurodegenerative diseases. Recently, a set of Spt5-Pol II small molecule inhibitors (SPIs) were reported, which selectively inhibit mutant huntingtin gene transcription. Inhibition mechanisms as well as interaction sites of these SPIs with Pol II and Spt5 are not entirely known. In this study, we predicted the binding sites of three selected SPIs at the Pol II-Spt5 interface by docking and molecular dynamics simulations. Two molecules out of three demonstrated strong binding with Spt5 and Pol II, while the other molecule was more loosely bound and sampled multiple binding sites. Strongly bound SPIs indirectly affected RNA and DNA dynamics at the exit site as DNA became more flexible while RNA was stabilized by increased interactions with Spt5. Our results suggest that the transcription inhibition mechanism induced by SPIs can be related to Spt5-nucleic acid interactions, which were altered to some extent with strong binding of SPIs.
Collapse
Affiliation(s)
- Adan Gallardo
- Department of Molecular and Cell Biology, University of California Merced, 5200 North Lake Rd, Merced, CA, 95343, USA
| | - Bercem Dutagaci
- Department of Molecular and Cell Biology, University of California Merced, 5200 North Lake Rd, Merced, CA, 95343, USA.
| |
Collapse
|
18
|
Keiff F, Jacques Dit Lapierre TJW, Bernal FA, Kloss F. Design and synthesis of benzofuran- and naphthalene-fused thiazinones as antimycobacterial agents. Arch Pharm (Weinheim) 2023; 356:e2300356. [PMID: 37667452 DOI: 10.1002/ardp.202300356] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2023] [Revised: 08/16/2023] [Accepted: 08/17/2023] [Indexed: 09/06/2023]
Abstract
Benzothiazinones (BTZs) have widely inspired medicinal chemistry and translational research due to their remarkable antitubercular potency and clinical potential. While most structure-activity relationship campaigns have largely focused on lateral chain modifications and substituents on the BTZ core, scaffold hopping strategies have been rarely investigated previously. In this work, we report the first example of ring expansion of the BTZ core toward benzofuran- and naphthalene-fused thiazinones. In vitro testing showed micromolar activity for both compounds, and molecular docking simulations provided insights into their reduced inhibitory capacity toward the enzymatic target (DprE1). Calculated electrochemical potentials revealed a lower susceptibility to reduction as opposed to BTZ drug candidates, in line with the mechanistic requirement for covalent binding.
Collapse
Affiliation(s)
- François Keiff
- Transfer Group Anti-infectives, Leibniz Institute for Natural Product Research and Infection Biology-Leibniz-HKI, Jena, Germany
| | | | - Freddy A Bernal
- Transfer Group Anti-infectives, Leibniz Institute for Natural Product Research and Infection Biology-Leibniz-HKI, Jena, Germany
| | - Florian Kloss
- Transfer Group Anti-infectives, Leibniz Institute for Natural Product Research and Infection Biology-Leibniz-HKI, Jena, Germany
| |
Collapse
|
19
|
Fernandez-Lopez L, Roda S, Robles-Martín A, Muñoz-Tafalla R, Almendral D, Ferrer M, Guallar V. Enhancing the Hydrolytic Activity of a Lipase towards Larger Triglycerides through Lid Domain Engineering. Int J Mol Sci 2023; 24:13768. [PMID: 37762071 PMCID: PMC10530837 DOI: 10.3390/ijms241813768] [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/03/2023] [Revised: 08/23/2023] [Accepted: 09/01/2023] [Indexed: 09/29/2023] Open
Abstract
Lipases have valuable potential for industrial use, particularly those mostly active against water-insoluble substrates, such as triglycerides composed of long-carbon chain fatty acids. However, in most cases, engineered variants often need to be constructed to achieve optimal performance for such substrates. Protein engineering techniques have been reported as strategies for improving lipase characteristics by introducing specific mutations in the cap domain of esterases or in the lid domain of lipases or through lid domain swapping. Here, we improved the lipase activity of a lipase (WP_075743487.1, or LipMRD) retrieved from the Marine Metagenomics MarRef Database and assigned to the Actinoalloteichus genus. The improvement was achieved through site-directed mutagenesis and by substituting its lid domain (FRGTEITQIKDWLTDA) with that of Rhizopus delemar lipase (previously R. oryzae; UniProt accession number, I1BGQ3) (FRGTNSFRSAITDIVF). The results demonstrated that the redesigned mutants gain activity against bulkier triglycerides, such as glyceryl tridecanoate and tridodecanoate, olive oil, coconut oil, and palm oil. Residue W89 (LipMRD numbering) appears to be key to the increase in lipase activity, an increase that was also achieved with lid swapping. This study reinforces the importance of the lid domains and their amino acid compositions in determining the substrate specificity of lipases, but the generalization of the lid domain swapping between lipases or the introduction of specific mutations in the lid domain to improve lipase activity may require further investigation.
Collapse
Affiliation(s)
- Laura Fernandez-Lopez
- Instituto de Catalisis y Petroleoquimica (ICP), CSIC, 28049 Madrid, Spain; (L.F.-L.); (D.A.)
| | - Sergi Roda
- Department of Life Sciences, Barcelona Supercomputing Center (BSC), 08034 Barcelona, Spain; (S.R.); (A.R.-M.); (R.M.-T.)
| | - Ana Robles-Martín
- Department of Life Sciences, Barcelona Supercomputing Center (BSC), 08034 Barcelona, Spain; (S.R.); (A.R.-M.); (R.M.-T.)
- PhD Programme, Faculty of Pharmacy and Food Science, Universitat de Barcelona (UB), 08007 Barcelona, Spain
| | - Rubén Muñoz-Tafalla
- Department of Life Sciences, Barcelona Supercomputing Center (BSC), 08034 Barcelona, Spain; (S.R.); (A.R.-M.); (R.M.-T.)
- PhD Programme, Faculty of Pharmacy and Food Science, Universitat de Barcelona (UB), 08007 Barcelona, Spain
| | - David Almendral
- Instituto de Catalisis y Petroleoquimica (ICP), CSIC, 28049 Madrid, Spain; (L.F.-L.); (D.A.)
| | - Manuel Ferrer
- Instituto de Catalisis y Petroleoquimica (ICP), CSIC, 28049 Madrid, Spain; (L.F.-L.); (D.A.)
| | - Víctor Guallar
- Department of Life Sciences, Barcelona Supercomputing Center (BSC), 08034 Barcelona, Spain; (S.R.); (A.R.-M.); (R.M.-T.)
- Institution for Research and Advanced Studies (ICREA), 08010 Barcelona, Spain
| |
Collapse
|
20
|
Li Y, Zhang Z, Wang R. HydraMap v.2: Prediction of Hydration Sites and Desolvation Energy with Refined Statistical Potentials. J Chem Inf Model 2023; 63:4749-4761. [PMID: 37433022 DOI: 10.1021/acs.jcim.3c00408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2023]
Abstract
The complex network of water molecules within the binding pocket of a target protein undergoes alterations upon ligand binding, presenting a significant challenge for conventional molecular modeling methods to accurately characterize and compute the associated energy changes. We have previously developed an empirical method, HydraMap (J. Chem. Inf. Model. 2020, 60, 4359-4375), which employs statistical potentials to predict hydration sites and compute desolvation energy, achieving a reasonable balance between accuracy and speed. In this work, we present its improved version, namely, HydraMap v.2. We updated the statistical potentials for protein-water interactions through an analysis of 17 042 crystal protein structures. We also introduced a new feature to evaluate ligand-water interactions by incorporating statistical potentials derived from the solvated structures of 9878 small organic molecules produced by molecular dynamics simulations. By combining these potentials, HydraMap v.2 can predict and compare the hydration sites in a binding pocket before and after ligand binding, identifying key water molecules involved in the binding process, such as those forming bridging hydrogen bonds and unstable ones that can be replaced. We demonstrated the application of HydraMap v.2 in explaining the structure-activity relationship of a panel of MCL-1 inhibitors. The desolvation energies calculated by summing the energy change of each hydration site before and after ligand binding showed good correlation with known ligand binding affinities on six target proteins. In conclusion, HydraMap v.2 offers a cost-effective solution for estimating the desolvation energy during protein-ligand binding and also is practical in guiding lead optimization in structure-based drug discovery.
Collapse
Affiliation(s)
- Yan Li
- Department of Medicinal Chemistry, School of Pharmacy, Fudan University, 826 Zhangheng Road, Shanghai 201203, People's Republic of China
| | - Zhe Zhang
- Department of Medicinal Chemistry, School of Pharmacy, Fudan University, 826 Zhangheng Road, Shanghai 201203, People's Republic of China
| | - Renxiao Wang
- Department of Medicinal Chemistry, School of Pharmacy, Fudan University, 826 Zhangheng Road, Shanghai 201203, People's Republic of China
| |
Collapse
|
21
|
Ye S, Han Y, Wei Z, Li J. Binding Affinity and Mechanisms of Potential Antidepressants Targeting Human NMDA Receptors. Molecules 2023; 28:4346. [PMID: 37298821 PMCID: PMC10254814 DOI: 10.3390/molecules28114346] [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/17/2023] [Revised: 05/12/2023] [Accepted: 05/22/2023] [Indexed: 06/12/2023] Open
Abstract
Depression, a mental disorder that plagues the world, is a burden on many families. There is a great need for new, fast-acting antidepressants to be developed. N-methyl-D-aspartic acid (NMDA) is an ionotropic glutamate receptor that plays an important role in learning and memory processes and its TMD region is considered as a potential target to treat depression. However, due to the unclear binding sites and pathways, the mechanism of drug binding lacks basic explanation, which brings great complexity to the development of new drugs. In this study, we investigated the binding affinity and mechanisms of an FDA-approved antidepressant (S-ketamine) and seven potential antidepressants (R-ketamine, memantine, lanicemine, dextromethorphan, Ro 25-6981, ifenprodil, and traxoprodil) targeting the NMDA receptor by ligand-protein docking and molecular dynamics simulations. The results indicated that Ro 25-6981 has the strongest binding affinity to the TMD region of the NMDA receptor among the eight selected drugs, suggesting its potential effective inhibitory effect. We also calculated the critical binding-site residues at the active site and found that residues Leu124 and Met63 contributed the most to the binding energy by decomposing the free energy contributions on a per-residue basis. We further compared S-ketamine and its chiral molecule, R-ketamine, and found that R-ketamine had a stronger binding capacity to the NMDA receptor. This study provides a computational reference for the treatment of depression targeting NMDA receptors, and the proposed results will provide potential strategies for further antidepressant development and is a useful resource for the future discovery of fast-acting antidepressant candidates.
Collapse
Affiliation(s)
- Simin Ye
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai Institute of Maternal-Fetal Medicine and Gynecologic Oncology, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China;
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai 200240, China;
| | - Yanqiang Han
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai 200240, China;
| | - Zhiyun Wei
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai Institute of Maternal-Fetal Medicine and Gynecologic Oncology, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China;
| | - Jinjin Li
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai Institute of Maternal-Fetal Medicine and Gynecologic Oncology, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China;
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai 200240, China;
| |
Collapse
|
22
|
Molteni C, Forni D, Cagliani R, Arrigoni F, Pozzoli U, De Gioia L, Sironi M. Selective events at individual sites underlie the evolution of monkeypox virus clades. Virus Evol 2023; 9:vead031. [PMID: 37305708 PMCID: PMC10256197 DOI: 10.1093/ve/vead031] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Revised: 03/31/2023] [Accepted: 05/12/2023] [Indexed: 06/13/2023] Open
Abstract
In endemic regions (West Africa and the Congo Basin), the genetic diversity of monkeypox virus (MPXV) is geographically structured into two major clades (Clades I and II) that differ in virulence and host associations. Clade IIb is closely related to the B.1 lineage, which is dominating a worldwide outbreak initiated in 2022. Lineage B.1 has however accumulated mutations of unknown significance that most likely result from apolipoprotein B mRNA editing catalytic polypeptide-like 3 (APOBEC3) editing. We applied a population genetics-phylogenetics approach to investigate the evolution of MPXV during historical viral spread in Africa and to infer the distribution of fitness effects. We observed a high preponderance of codons evolving under strong purifying selection, particularly in viral genes involved in morphogenesis and replication or transcription. However, signals of positive selection were also detected and were enriched in genes involved in immunomodulation and/or virulence. In particular, several genes showing evidence of positive selection were found to hijack different steps of the cellular pathway that senses cytosolic DNA. Also, a few selected sites in genes that are not directly involved in immunomodulation are suggestive of antibody escape or other immune-mediated pressures. Because orthopoxvirus host range is primarily determined by the interaction with the host immune system, we suggest that the positive selection signals represent signatures of host adaptation and contribute to the different virulence of Clade I and II MPXVs. We also used the calculated selection coefficients to infer the effects of mutations that define the predominant human MPXV1 (hMPXV1) lineage B.1, as well as the changes that have been accumulating during the worldwide outbreak. Results indicated that a proportion of deleterious mutations were purged from the predominant outbreak lineage, whose spread was not driven by the presence of beneficial changes. Polymorphic mutations with a predicted beneficial effect on fitness are few and have a low frequency. It remains to be determined whether they have any significance for ongoing virus evolution.
Collapse
Affiliation(s)
| | | | | | - Federica Arrigoni
- Department of Biotechnology and Biosciences, University of Milan-Bicocca, Piazza della scienza, Milan 20126, Italy
| | - Uberto Pozzoli
- Scientific Institute IRCCS E. MEDEA, Bioinformatics, Via don Luigi Monza, Bosisio Parini 23842, Italy
| | - Luca De Gioia
- Department of Biotechnology and Biosciences, University of Milan-Bicocca, Piazza della scienza, Milan 20126, Italy
| | - Manuela Sironi
- Scientific Institute IRCCS E. MEDEA, Bioinformatics, Via don Luigi Monza, Bosisio Parini 23842, Italy
| |
Collapse
|
23
|
Zuzic L, Marzinek JK, Anand GS, Warwicker J, Bond PJ. A pH-dependent cluster of charges in a conserved cryptic pocket on flaviviral envelopes. eLife 2023; 12:82447. [PMID: 37144875 PMCID: PMC10162804 DOI: 10.7554/elife.82447] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2022] [Accepted: 04/18/2023] [Indexed: 05/06/2023] Open
Abstract
Flaviviruses are enveloped viruses which include human pathogens that are predominantly transmitted by mosquitoes and ticks. Some, such as dengue virus, exhibit the phenomenon of antibody-dependent enhancement (ADE) of disease, making vaccine-based routes of fighting infections problematic. The pH-dependent conformational change of the envelope (E) protein required for fusion between the viral and endosomal membranes is an attractive point of inhibition by antivirals as it has the potential to diminish the effects of ADE. We examined six flaviviruses by employing large-scale molecular dynamics (MD) simulations of raft systems that represent a substantial portion of the flaviviral envelope. We utilised a benzene-mapping approach that led to a discovery of shared hotspots and conserved cryptic sites. A cryptic pocket previously shown to bind a detergent molecule exhibited strain-specific characteristics. An alternative conserved cryptic site at the E protein domain interfaces showed a consistent dynamic behaviour across flaviviruses and contained a conserved cluster of ionisable residues. Constant-pH simulations revealed cluster and domain-interface disruption under low pH conditions. Based on this, we propose a cluster-dependent mechanism that addresses inconsistencies in the histidine-switch hypothesis and highlights the role of cluster protonation in orchestrating the domain dissociation pivotal for the formation of the fusogenic trimer.
Collapse
Affiliation(s)
- Lorena Zuzic
- Bioinformatics Institute (A*STAR), Singapore, Singapore
- Department of Chemistry, Manchester Institute of Biotechnology, The University of Manchester, Manchester, United Kingdom
| | | | - Ganesh S Anand
- Department of Biological Sciences, 16 Science Drive 4, National University of Singapore, Singapore, Singapore
- Department of Chemistry, The Pennsylvania State University, University Park, United States
| | - Jim Warwicker
- School of Biological Sciences, Faculty of Biology, Medicine and Health, Manchester Institute of Biotechnology, The University of Manchester, Manchester, United Kingdom
| | - Peter J Bond
- Bioinformatics Institute (A*STAR), Singapore, Singapore
- Department of Biological Sciences, 16 Science Drive 4, National University of Singapore, Singapore, Singapore
| |
Collapse
|
24
|
Blatter M, Meylan C, Cléry A, Giambruno R, Nikolaev Y, Heidecker M, Solanki JA, Diaz MO, Gabellini D, Allain FHT. RNA binding induces an allosteric switch in Cyp33 to repress MLL1-mediated transcription. SCIENCE ADVANCES 2023; 9:eadf5330. [PMID: 37075125 PMCID: PMC10115415 DOI: 10.1126/sciadv.adf5330] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Mixed-lineage leukemia 1 (MLL1) is a transcription activator of the HOX family, which binds to specific epigenetic marks on histone H3 through its third plant homeodomain (PHD3) domain. Through an unknown mechanism, MLL1 activity is repressed by cyclophilin 33 (Cyp33), which binds to MLL1 PHD3. We determined solution structures of Cyp33 RNA recognition motif (RRM) free, bound to RNA, to MLL1 PHD3, and to both MLL1 and the histone H3 lysine N6-trimethylated. We found that a conserved α helix, amino-terminal to the RRM domain, adopts three different positions facilitating a cascade of binding events. These conformational changes are triggered by Cyp33 RNA binding and ultimately lead to MLL1 release from the histone mark. Together, our mechanistic findings rationalize how Cyp33 binding to MLL1 can switch chromatin to a transcriptional repressive state triggered by RNA binding as a negative feedback loop.
Collapse
Affiliation(s)
- Markus Blatter
- Department of Biology, Institute of Biochemistry, ETH Zurich, 8093 Zurich, Switzerland
- Corresponding author. (F.H.-T.A.); (M.B.)
| | - Charlotte Meylan
- Department of Biology, Institute of Biochemistry, ETH Zurich, 8093 Zurich, Switzerland
| | - Antoine Cléry
- Department of Biology, Institute of Biochemistry, ETH Zurich, 8093 Zurich, Switzerland
| | - Roberto Giambruno
- Gene Expression and Muscular Dystrophy Unit, Division of Genetics and Cell Biology, IRCCS San Raffaele Scientific Institute, Milan 20132, Italy
| | - Yaroslav Nikolaev
- Department of Biology, Institute of Biochemistry, ETH Zurich, 8093 Zurich, Switzerland
| | - Michel Heidecker
- Department of Biology, Institute of Biochemistry, ETH Zurich, 8093 Zurich, Switzerland
| | - Jessica Arvindbhai Solanki
- Department of Microbiology and Immunology, Stritch School of Medicine, Loyola University of Chicago Medical Center, University of Chicago, Chicago, IL, USA
| | - Manuel O. Diaz
- Department of Microbiology and Immunology, Stritch School of Medicine, Loyola University of Chicago Medical Center, University of Chicago, Chicago, IL, USA
| | - Davide Gabellini
- Gene Expression and Muscular Dystrophy Unit, Division of Genetics and Cell Biology, IRCCS San Raffaele Scientific Institute, Milan 20132, Italy
| | - Frédéric H.-T. Allain
- Department of Biology, Institute of Biochemistry, ETH Zurich, 8093 Zurich, Switzerland
- Corresponding author. (F.H.-T.A.); (M.B.)
| |
Collapse
|
25
|
Roda S, Terholsen H, Meyer JRH, Cañellas-Solé A, Guallar V, Bornscheuer U, Kazemi M. AsiteDesign: a Semirational Algorithm for an Automated Enzyme Design. J Phys Chem B 2023; 127:2661-2670. [PMID: 36944360 PMCID: PMC10068746 DOI: 10.1021/acs.jpcb.2c07091] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/23/2023]
Abstract
With advances in protein structure predictions, the number of available high-quality structures has increased dramatically. In light of these advances, structure-based enzyme engineering is expected to become increasingly important for optimizing biocatalysts for industrial processes. Here, we present AsiteDesign, a Monte Carlo-based protocol for structure-based engineering of active sites. AsiteDesign provides a framework for introducing new catalytic residues in a given binding pocket to either create a new catalytic activity or alter the existing one. AsiteDesign is implemented using pyRosetta and incorporates enhanced sampling techniques to efficiently explore the search space. The protocol was tested by designing an alternative catalytic triad in the active site of Pseudomonas fluorescens esterase (PFE). The designed variant was experimentally verified to be active, demonstrating that AsiteDesign can find alternative catalytic triads. Additionally, the AsiteDesign protocol was employed to enhance the hydrolysis of a bulky chiral substrate (1-phenyl-2-pentyl acetate) by PFE. The experimental verification of the designed variants demonstrated that F158L/F198A and F125A/F158L mutations increased the hydrolysis of 1-phenyl-2-pentyl acetate from 8.9 to 66.7 and 23.4%, respectively, and reversed the enantioselectivity of the enzyme from (R) to (S)-enantiopreference, with 32 and 55% enantiomeric excess (ee), respectively.
Collapse
Affiliation(s)
- Sergi Roda
- Barcelona Supercomputing Center (BSC), Plaça d'Eusebi Güell, 1-3, Barcelona 08034, Spain
| | - Henrik Terholsen
- Department of Biotechnology & Enzyme Catalysis, Institute of Biochemistry, University of Greifswald, Felix-Hausdorff-Str. 4, D-17487 Greifswald, Germany
| | - Jule Ruth Heike Meyer
- Department of Biotechnology & Enzyme Catalysis, Institute of Biochemistry, University of Greifswald, Felix-Hausdorff-Str. 4, D-17487 Greifswald, Germany
| | - Albert Cañellas-Solé
- Barcelona Supercomputing Center (BSC), Plaça d'Eusebi Güell, 1-3, Barcelona 08034, Spain
| | - Victor Guallar
- Barcelona Supercomputing Center (BSC), Plaça d'Eusebi Güell, 1-3, Barcelona 08034, Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Passeig de Lluís Companys, 23, Barcelona 08010, Spain
| | - Uwe Bornscheuer
- Department of Biotechnology & Enzyme Catalysis, Institute of Biochemistry, University of Greifswald, Felix-Hausdorff-Str. 4, D-17487 Greifswald, Germany
| | - Masoud Kazemi
- Barcelona Supercomputing Center (BSC), Plaça d'Eusebi Güell, 1-3, Barcelona 08034, Spain
- Biomatter Designs, Žirmu̅n̨ g. 139A, Vilnius 09120, Lithuania
| |
Collapse
|
26
|
Parkin D, Takano M. The Intrinsic Radius as a Key Parameter in the Generalized Born Model to Adjust Protein-Protein Electrostatic Interaction. Int J Mol Sci 2023; 24:ijms24054700. [PMID: 36902130 PMCID: PMC10003186 DOI: 10.3390/ijms24054700] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2023] [Revised: 01/30/2023] [Accepted: 02/19/2023] [Indexed: 03/05/2023] Open
Abstract
The generalized Born (GB) model is an extension of the continuum dielectric theory of Born solvation energy and is a powerful method for accelerating the molecular dynamic (MD) simulations of charged biological molecules in water. While the effective dielectric constant of water that varies as a function of the separation distance between solute molecules is incorporated into the GB model, adjustment of the parameters is indispensable for accurate calculation of the Coulomb (electrostatic) energy. One of the key parameters is the lower limit of the spatial integral of the energy density of the electric field around a charged atom, known as the intrinsic radius ρ. Although ad hoc adjustment of ρ has been conducted to improve the Coulombic (ionic) bond stability, the physical mechanism by which ρ affects the Coulomb energy remains unclear. Via energetic analysis of three differently sized systems, here, we clarify that the Coulomb bond stability increases with increasing ρ and that the increased stability is caused by the interaction energy term, not by the self-energy (desolvation energy) term, as was supposed previously. Our results suggest that the use of larger values for the intrinsic radii of hydrogen and oxygen atoms, together with the use of a relatively small value for the spatial integration cutoff in the GB model, can better reproduce the Coulombic attraction between protein molecules.
Collapse
Affiliation(s)
- Dan Parkin
- Research Institute for Science and Engineering, Waseda University, Okubo 3-4-1, Sinjuku-ku, Tokyo 169-8555, Japan
| | - Mitsunori Takano
- Research Institute for Science and Engineering, Waseda University, Okubo 3-4-1, Sinjuku-ku, Tokyo 169-8555, Japan
- Department of Pure and Applied Physics, Waseda University, Okubo 3-4-1, Sinjuku-ku, Tokyo 169-8555, Japan
- Correspondence: ; Tel.: +81-3-5286-3512
| |
Collapse
|
27
|
Synthesis, Characterization, and Biological Evaluation of 2-(N-((2'-(2H-tetrazole-5-yl)-[1,1'-biphenyl]-4yl)-methyl)-pentanamido)-3-methyl Butanoic Acid Derivatives. Molecules 2023; 28:molecules28041908. [PMID: 36838897 PMCID: PMC9959676 DOI: 10.3390/molecules28041908] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Revised: 02/10/2023] [Accepted: 02/14/2023] [Indexed: 02/19/2023] Open
Abstract
This study aimed to evaluate 2-(N-((2'-(2H-tetrazole-5-yl)-[1,1'-biphenyl]-4yl)-methyl)-pentanamido)-3-methyl butanoic acid-based ester derivatives as a new class of angiotensin-II receptor antagonists. For this purpose, a series of compounds were synthesized using a variety of phenols. Their chemical characterization was established by FTIR, 1HNMR, and 13CNMR techniques. The biological activities including antioxidant potentials using the DPPH assay, the antihypertensive assay, the urease enzyme inhibition assay, and the antibacterial assay using agar well diffusion methods were performed. All the new compounds showed significant free radical scavenging potentials more than the parent drug while retaining antihypertensive potentials along with urease inhibition properties. However, the AV2 test compound was found to be the most potent against hypertension. Most of the synthesized analogs showed urease inhibitory actions. Molecular docking studies were performed for all the active analogs to decode the binding detail of the ligands with receptors of the enzyme's active site.
Collapse
|
28
|
Kunche L, Natarajan U. Conformations and Solvation of Synthetic Polymers in Water by Generalized Born Implicit-Solvent Molecular Dynamics Simulations: Stereoisomers of Poly(acrylic acid) and Poly(methacrylic acid). J Phys Chem B 2023; 127:1244-1253. [PMID: 36705523 DOI: 10.1021/acs.jpcb.2c06658] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
We present the GB-OBC model as an approach for implicit-solvent MD simulations of a synthetic macromolecule in water. The model is tested and found to be successful in reproducing the chain dimensions and predicting the free energy of solvation of carboxylic acid vinyl polymers. The influence of stereochemistry and the hydrophobic nature of the polymer was investigated as a function of chain length (20 < N < 600) for poly(acrylic acid) (PAA) and poly(methacrylic acid) (PMA). The dimensionless parameters of the GB-OBC model were parameterized to be applicable to PAA and PMA. Scaling relations for chain dimensions obtained using implicit-solvent MD simulations in this study are in good agreement with those from experiments, theory of solvated chains in good solvents, and all-atom MD simulations in explicit water. Results show that ⟨Rg2⟩/NL2 is greater for the atactic chain as compared to the isotactic chain, for PAA as well as PMA. ⟨Rg2⟩/NL2 values of chains attain constancy in water for N = 200, with the values being greater for PMA. The PMA chain is conformationally more perturbed than the PAA chain, for both isotactic and atactic stereochemistry. The solvation free energy ΔGhyd of PAA and PMA in water is negative for all chain lengths (N = 20-600) and becomes more favorable with an increase in molecular weight. The ΔGhyd values for isotactic and atactic chains are identical at lower values of N but differ slightly for N > 300. Irrespective of the hydrophobic nature of the polymer, the atactic chain is thermodynamically more soluble in water as compared to the isotactic chain. The isotactic chain is less hydrophilic as compared to the atactic chain due to the closer proximity of the COOH groups along the backbone. This implicit solvent method is an effective way to accurately simulate the configurational properties and solvation of synthetic polymers in water.
Collapse
Affiliation(s)
- Lakshmikumar Kunche
- Macromolecular Modeling and Simulation Lab, Department of Chemical Engineering, Indian Institute of Technology (IIT) Madras, Chennai600036, India
| | - Upendra Natarajan
- Macromolecular Modeling and Simulation Lab, Department of Chemical Engineering, Indian Institute of Technology (IIT) Madras, Chennai600036, India
| |
Collapse
|
29
|
Chen L, Huang X, Li Y, Zhao B, Liang M, Wang R. Structural and energetic basis of interaction between human estrogen-related receptor γ and environmental endocrine disruptors from multiple molecular dynamics simulations and free energy predictions. JOURNAL OF HAZARDOUS MATERIALS 2023; 443:130174. [PMID: 36265380 DOI: 10.1016/j.jhazmat.2022.130174] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/24/2022] [Revised: 10/04/2022] [Accepted: 10/09/2022] [Indexed: 06/16/2023]
Abstract
Environmental endocrine disruptors (EEDs), a class of molecules that are widespread in our environment, may adversely affect the endocrine system. Exploring the interactions between these compounds and their potential targets is important for assessing their role in the organism. Focused on the human estrogen-related receptor γ (hERRγ) with BPA, BPB, HPTE, BPE, BP(2,2)(Et), and BP(2,2)(MeO) complexes, respectively, we groped for the mechanisms of conformational changes and interactions of hERRγ when binding to these six EEDs by combining multiple molecular dynamics (MD) simulations with energy prediction (MM-PBSA and solvated interaction energy (SIE)). Dynamics analysis results revealed these six EEDs have different effects on the internal dynamics of hERRγ, resulting in significant changes in the interaction of key residues around Leu268, Val313, Leu345, and Phe435 with EEDs, and thus affected its binding energy with these EEDs. The energy calculations further demonstrated that van der Waals interactions are critical for these EEDs binding to hERRγ. These results present detailed molecular insight into the interaction features between EEDs and hERRγ and help guide the search for safer alternatives to BPA.
Collapse
Affiliation(s)
- Lin Chen
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China.
| | - Xu Huang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China
| | - Yufei Li
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China
| | - Bing Zhao
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China; Heilongjiang Provincial Key Laboratory of Surface Active Agent and Auxiliary, Qiqihar 161006, PR China
| | - Min Liang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China
| | - Ruige Wang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China
| |
Collapse
|
30
|
Wang R, Sun H, Chen W, Zhao B, Chen L. Molecular basis of ssDNA recognition by RBM45 protein of neurodegenerative disease from multiple molecular dynamics simulations and energy predictions. J Mol Graph Model 2023; 118:108377. [PMID: 36435031 DOI: 10.1016/j.jmgm.2022.108377] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2022] [Revised: 11/09/2022] [Accepted: 11/14/2022] [Indexed: 11/21/2022]
Abstract
Neurodegenerative diseases (NDD) are a group of cognitive and behavioral disorders characterized by progressive loss of neuronal structure and function. As the population ages, the incidence is getting higher and higher, but there is currently no effective treatment. The details of RNA/DNA recognition by the RNA-binding protein RBM45 closely related to neurodegenerative diseases through its two tandem RNA-recognition domains at its N-terminus have important implications for structure-based drug discovery against degenerative diseases. To explore the key characteristics of interaction between ssDNA and RBM45, we performed multiple molecular dynamics (MD) simulations along with MM-PBSA energy prediction on the complexes of wild type (WT) and three mutant RBM45s (K100A, F124A/Y165A, and F29A/F70A/F124A/Y165A) with ssDNA, respectively. The findings suggest that these mutated residues of RBM45 modify the interaction of their surrounding residues with ssDNA, thereby affecting RBM45 protein binding to ssDNA. In contrast with WT RBM45 protein, variations in van der Waals and electrostatic interactions with ssDNA caused by these three RBM45 mutants are critical to affect binding between them. In addition, energy analysis showed that RBM45 is a specific ssDNA-binding protein. The results of our work provide valuable theoretical guidance for the design effective drugs of NDD.
Collapse
Affiliation(s)
- Ruige Wang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, PR China
| | - Han Sun
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, PR China
| | - Wei Chen
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, PR China.
| | - Bing Zhao
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, PR China; Heilongjiang Provincial Key Laboratory of Surface Active Agent and Auxiliary, Qiqihar, 161006, PR China
| | - Lin Chen
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, PR China.
| |
Collapse
|
31
|
Chang Y, Hawkins BA, Du JJ, Groundwater PW, Hibbs DE, Lai F. A Guide to In Silico Drug Design. Pharmaceutics 2022; 15:pharmaceutics15010049. [PMID: 36678678 PMCID: PMC9867171 DOI: 10.3390/pharmaceutics15010049] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2022] [Revised: 12/16/2022] [Accepted: 12/17/2022] [Indexed: 12/28/2022] Open
Abstract
The drug discovery process is a rocky path that is full of challenges, with the result that very few candidates progress from hit compound to a commercially available product, often due to factors, such as poor binding affinity, off-target effects, or physicochemical properties, such as solubility or stability. This process is further complicated by high research and development costs and time requirements. It is thus important to optimise every step of the process in order to maximise the chances of success. As a result of the recent advancements in computer power and technology, computer-aided drug design (CADD) has become an integral part of modern drug discovery to guide and accelerate the process. In this review, we present an overview of the important CADD methods and applications, such as in silico structure prediction, refinement, modelling and target validation, that are commonly used in this area.
Collapse
Affiliation(s)
- Yiqun Chang
- Sydney Pharmacy School, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia
| | - Bryson A. Hawkins
- Sydney Pharmacy School, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia
| | - Jonathan J. Du
- Department of Biochemistry, Emory University School of Medicine, Atlanta, GA 30322, USA
| | - Paul W. Groundwater
- Sydney Pharmacy School, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia
| | - David E. Hibbs
- Sydney Pharmacy School, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia
| | - Felcia Lai
- Sydney Pharmacy School, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia
- Correspondence:
| |
Collapse
|
32
|
Doğan N, Yavuz SÇ, Sahin K, Orhan MD, Muhammed HK, Calis S, Küp FÖ, Avsar T, Akkoc S, Tapera M, Sahin O, Kilic T, Durdagi S, Saripinar E. Synthesis, Characterization, Biological Activity and Molecular Modeling Studies of Novel Aminoguanidine Derivatives. ChemistrySelect 2022. [DOI: 10.1002/slct.202202819] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Affiliation(s)
- Nuriye Doğan
- Department of Chemistry Erciyes University Kayseri Türkiye
| | - Sevtap Çağlar Yavuz
- Department of Medical Services and Technicians İliç Dursun Yıldırım Vocational School Erzincan Binali Yıldırım University 24700 Erzincan Türkiye
| | - Kader Sahin
- Computational Biology and Molecular Simulations Laboratory Department of Biophysics School of Medicine Bahcesehir University Istanbul Türkiye
- Department of Analytical Chemistry School of Pharmacy Bahcesehir University Istanbul Türkiye
| | - Muge Didem Orhan
- Department of Medical Biology School of Medicine Bahcesehir University Istanbul Türkiye
| | | | - Seyma Calis
- Department of Medical Biology School of Medicine Bahcesehir University Istanbul Türkiye
| | - Fatma Öztürk Küp
- Department of Biology Faculty of Science Erciyes University Kayseri Türkiye
| | - Timucin Avsar
- Department of Medical Biology School of Medicine Bahcesehir University Istanbul Türkiye
| | - Senem Akkoc
- Department of Basic Pharmaceutical Sciences Faculty of Pharmacy Suleyman Demirel University Isparta Türkiye
| | - Michael Tapera
- Department of Chemistry Erciyes University Kayseri Türkiye
| | - Onur Sahin
- Scientific and Technological Research Application and Research Center Sinop University Sinop Türkiye
| | - Turker Kilic
- Deperment of Neurosurgery School of Medicine Bahcesehir University Istanbul Türkiye
| | - Serdar Durdagi
- Computational Biology and Molecular Simulations Laboratory Department of Biophysics School of Medicine Bahcesehir University Istanbul Türkiye
- Department of Pharmaceutical Chemistry School of Pharmacy Bahcesehir University Istanbul Türkiye
| | - Emin Saripinar
- Department of Chemistry Erciyes University Kayseri Türkiye
| |
Collapse
|
33
|
Adhikari RS, Parambathu AV, Chapman WG, Asthagiri DN. Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory. J Phys Chem B 2022; 126:9607-9616. [PMID: 36354351 DOI: 10.1021/acs.jpcb.2c05725] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)─an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration─physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.
Collapse
Affiliation(s)
- Rohan S Adhikari
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas77005, United States
| | - Arjun Valiya Parambathu
- Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, Delaware19711, United States
| | - Walter G Chapman
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas77005, United States
| | | |
Collapse
|
34
|
Goodman J, Attwood D, Kiely J, Coladas Mato P, Luxton R. Modeling Peptide Nucleic Acid Binding Enthalpies Using MM-GBSA. J Phys Chem B 2022; 126:9528-9538. [PMID: 36375178 PMCID: PMC9706560 DOI: 10.1021/acs.jpcb.2c05547] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
The binding enthalpies of peptide nucleic acid (PNA) homoduplexes were predicted using a molecular mechanics generalized Born surface area approach. Using the nucleic acid nearest-neighbor model, these were decomposed into sequence parameters which could replicate the enthalpies from thermal melting experiments with a mean error of 8.7%. These results present the first systematic computational investigation into the relationship between sequence and binding energy for PNA homoduplexes and identified a stabilizing helix initiation enthalpy not observed for nucleic acids with phosphoribose backbones.
Collapse
Affiliation(s)
- Jack Goodman
- University
of the West of England, BristolBS16 1QY, U.K.,
| | - David Attwood
- University
of the West of England, BristolBS16 1QY, U.K.
| | - Janice Kiely
- University
of the West of England, BristolBS16 1QY, U.K.
| | | | - Richard Luxton
- University
of the West of England, BristolBS16 1QY, U.K.
| |
Collapse
|
35
|
Fernandez-Lopez L, Roda S, Gonzalez-Alfonso JL, Plou FJ, Guallar V, Ferrer M. Design and Characterization of In-One Protease-Esterase PluriZyme. Int J Mol Sci 2022; 23:13337. [PMID: 36362119 PMCID: PMC9655419 DOI: 10.3390/ijms232113337] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2022] [Revised: 10/25/2022] [Accepted: 10/29/2022] [Indexed: 10/14/2023] Open
Abstract
Proteases are abundant in prokaryotic genomes (~10 per genome), but their recovery encounters expression problems, as only 1% can be produced at high levels; this value differs from that of similarly abundant esterases (1-15 per genome), 50% of which can be expressed at good levels. Here, we design a catalytically efficient artificial protease that can be easily produced. The PluriZyme EH1AB1 with two active sites supporting the esterase activity was employed. A Leu24Cys mutation in EH1AB1, remodelled one of the esterase sites into a proteolytic one through the incorporation of a catalytic dyad (Cys24 and His214). The resulting artificial enzyme, EH1AB1C, efficiently hydrolysed (azo)casein at pH 6.5-8.0 and 60-70 °C. The presence of both esterase and protease activities in the same scaffold allowed the one-pot cascade synthesis (55.0 ± 0.6% conversion, 24 h) of L-histidine methyl ester from the dipeptide L-carnosine in the presence of methanol. This study demonstrates that active sites supporting proteolytic activity can be artificially introduced into an esterase scaffold to design easy-to-produce in-one protease-esterase PluriZymes for cascade reactions, namely, the synthesis of amino acid esters from dipeptides. It is also possible to design artificial proteases with good production yields, in contrast to natural proteases that are difficult to express.
Collapse
Affiliation(s)
| | - Sergi Roda
- Department of Life Sciences, Barcelona Supercomputing Center (BSC), 08034 Barcelona, Spain
| | | | | | - Víctor Guallar
- Department of Life Sciences, Barcelona Supercomputing Center (BSC), 08034 Barcelona, Spain
- Institution for Research and Advanced Studies (ICREA), 08010 Barcelona, Spain
| | - Manuel Ferrer
- Department of Applied Biocatalysis, ICP, CSIC, 28049 Madrid, Spain
| |
Collapse
|
36
|
Williams-Noonan BJ, Speer MN, Le TC, Sadek MM, Thompson PE, Norton RS, Yuriev E, Barlow N, Chalmers DK, Yarovsky I. Membrane Permeating Macrocycles: Design Guidelines from Machine Learning. J Chem Inf Model 2022; 62:4605-4619. [PMID: 36178379 DOI: 10.1021/acs.jcim.2c00809] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The ability to predict cell-permeable candidate molecules has great potential to assist drug discovery projects. Large molecules that lie beyond the Rule of Five (bRo5) are increasingly important as drug candidates and tool molecules for chemical biology. However, such large molecules usually do not cross cell membranes and cannot access intracellular targets or be developed as orally bioavailable drugs. Here, we describe a random forest (RF) machine learning model for the prediction of passive membrane permeation rates developed using a set of over 1000 bRo5 macrocyclic compounds. The model is based on easily calculated chemical features/descriptors as independent variables. Our random forest (RF) model substantially outperforms a multiple linear regression model based on the same features and achieves better performance metrics than previously reported models using the same underlying data. These features include: (1) polar surface area in water, (2) the octanol-water partitioning coefficient, (3) the number of hydrogen-bond donors, (4) the sum of the topological distances between nitrogen atoms, (5) the sum of the topological distances between nitrogen and oxygen atoms, and (6) the multiple molecular path count of order 2. The last three features represent molecular flexibility, the ability of the molecule to adopt different conformations in the aqueous and membrane interior phases, and the molecular "chameleonicity." Guided by the model, we propose design guidelines for membrane-permeating macrocycles. It is anticipated that this model will be useful in guiding the design of large, bioactive molecules for medicinal chemistry and chemical biology applications.
Collapse
Affiliation(s)
- Billy J Williams-Noonan
- School of Engineering, RMIT University, Melbourne3001, Australia.,Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville3052, Australia
| | - Melissa N Speer
- University of Melbourne, Faculty of Engineering and Information Technology, Carlton3053, Australia
| | - Tu C Le
- School of Engineering, RMIT University, Melbourne3001, Australia
| | - Maiada M Sadek
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville3052, Australia
| | - Philip E Thompson
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville3052, Australia
| | - Raymond S Norton
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville3052, Australia.,ARC Centre for Fragment-Based Design, Monash University, Parkville, 3052, Australia
| | - Elizabeth Yuriev
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville3052, Australia
| | - Nicholas Barlow
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville3052, Australia
| | - David K Chalmers
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville3052, Australia
| | - Irene Yarovsky
- School of Engineering, RMIT University, Melbourne3001, Australia
| |
Collapse
|
37
|
Argikar U, Blatter M, Bednarczyk D, Chen Z, Cho YS, Doré M, Dumouchel JL, Ho S, Hoegenauer K, Kawanami T, Mathieu S, Meredith E, Möbitz H, Murphy SK, Parthasarathy S, Soldermann CP, Santos J, Silver S, Skolnik S, Stojanovic A. Paradoxical Increase of Permeability and Lipophilicity with the Increasing Topological Polar Surface Area within a Series of PRMT5 Inhibitors. J Med Chem 2022; 65:12386-12402. [PMID: 36069672 DOI: 10.1021/acs.jmedchem.2c01068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An imidazolone → triazolone replacement addressed the limited passive permeability of a series of protein arginine methyl transferase 5 (PRMT5) inhibitors. This increase in passive permeability was unexpected given the increase in the hydrogen bond acceptor (HBA) count and topological polar surface area (TPSA), two descriptors that are typically inversely correlated with permeability. Quantum mechanics (QM) calculations revealed that this unusual effect was due to an electronically driven disconnect between TPSA and 3D-PSA, which manifests in a reduction in overall HBA strength as indicated by the HBA moment descriptor from COSMO-RS (conductor-like screening model for real solvation). HBA moment was subsequently deployed as a design parameter leading to the discovery of inhibitors with not only improved passive permeability but also reduced P-glycoprotein (P-gp) transport. Our case study suggests that hidden polarity as quantified by TPSA-3DPSA can be rationally designed through QM calculations.
Collapse
Affiliation(s)
- Upendra Argikar
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | - Markus Blatter
- Novartis Institutes for BioMedical Research, Basel 4002, Switzerland
| | - Dallas Bednarczyk
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | - Zhuoliang Chen
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | - Young Shin Cho
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | - Michaël Doré
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | - Jennifer L Dumouchel
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | - Samuel Ho
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | | | - Toshio Kawanami
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | - Simon Mathieu
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | - Erik Meredith
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | - Henrik Möbitz
- Novartis Institutes for BioMedical Research, Basel 4002, Switzerland
| | - Stephen K Murphy
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | | | | | - Jobette Santos
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | - Serena Silver
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | - Suzanne Skolnik
- Novartis Institutes for BioMedical Research, Cambridge, Massachusetts 02139, United States
| | | |
Collapse
|
38
|
Sinha S, Tam B, Wang SM. Applications of Molecular Dynamics Simulation in Protein Study. MEMBRANES 2022; 12:membranes12090844. [PMID: 36135863 PMCID: PMC9505860 DOI: 10.3390/membranes12090844] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/14/2022] [Revised: 08/24/2022] [Accepted: 08/25/2022] [Indexed: 05/29/2023]
Abstract
Molecular Dynamics (MD) Simulations is increasingly used as a powerful tool to study protein structure-related questions. Starting from the early simulation study on the photoisomerization in rhodopsin in 1976, MD Simulations has been used to study protein function, protein stability, protein-protein interaction, enzymatic reactions and drug-protein interactions, and membrane proteins. In this review, we provide a brief review for the history of MD Simulations application and the current status of MD Simulations applications in protein studies.
Collapse
|
39
|
Folescu DE, Onufriev AV. A Closed-Form, Analytical Approximation for Apparent Surface Charge and Electric Field of Molecules. ACS OMEGA 2022; 7:26123-26136. [PMID: 35936397 PMCID: PMC9352323 DOI: 10.1021/acsomega.2c01484] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/12/2022] [Accepted: 06/10/2022] [Indexed: 06/15/2023]
Abstract
Closed-form, analytical approximations for electrostatic properties of molecules are of unique value as these can provide computational speed, versatility, and physical insight. Here, we have derived a simple, closed-form formula for the apparent surface charge (ASC) as well as for the electric field generated by a molecular charge distribution in aqueous solution. The approximation, with no fitted parameters, was tested against numerical solutions of the Poisson equation, where it has produced a significant speed-up. For neutral small molecules, the hydration free energies estimated from the closed-form ASC formula are within 0.8 kcal/mol RMSD from the numerical Poisson reference; the electric field at the surface is in quantitative agreement with the reference. Performance of the approximation was also tested on larger structures, including a protein, a DNA fragment, and a viral receptor-target complex. For all structures tested, a near-quantitative agreement with the numerical Poisson reference was achieved, except in regions of high negative curvature, where the new approximation is still qualitatively correct. A unique efficiency feature of the proposed "source-based″ closed-form approximation is that the ASC and electric field can be estimated individually at any point or surface patch, without the need to obtain the full global solution. An open-source software implementation of the method is available: https://people.cs.vt.edu/~onufriev/CODES/aasc.zip.
Collapse
Affiliation(s)
- Dan E. Folescu
- Department
of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Alexey V. Onufriev
- Department
of Computer Science, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department
of Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
- Center
for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
40
|
Dragan P, Atzei A, Sanmukh SG, Latek D. Computational and experimental approaches to probe GPCR activation and signaling. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2022; 193:1-36. [PMID: 36357073 DOI: 10.1016/bs.pmbts.2022.06.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
G protein-coupled receptors (GPCRs) regulate different physiological functions, e.g., sensation, growth, digestion, reproductivity, nervous and immune systems response, and many others. In eukaryotes, they are also responsible for intercellular communication in response to pathogens. The major primary messengers binding to these cell-surface receptors constitute small-molecule or peptide hormones and neurotransmitters, nucleotides, lipids as well as small proteins. The simplicity of the way how GPCR signaling can be regulated by their endogenous agonists prompted the usage of GPCRs as major drug targets in modern pharmacology. Drugs targeting GPCRs inhibit pathological processes at the very beginning. This enables to significantly reduce the occurrence of morphological changes caused by diseases. Until recently, X-ray crystallography was the method of the first choice to obtain high-resolution structural information about GPCRs. Following X-ray crystallography, cryo-EM gained attention in GPCR studies as a quick and low-cost alternative. FRET microscopy is also widely used for GPCRs in the analysis of protein-protein interactions (PPIs) in intact cells as well as for screening purposes. Regarding computational methods, molecular dynamics (MD) for many years has proven its usefulness in studying the GPCR activation. MODELLER and Rosetta were widely used to generate preliminary homology models of GPCRs for MD simulation systems. Apart from the conventional all-atom approach with explicitly defined solvent, also other techniques have been applied to GPCRs, e.g., MARTINI or hybrid methods involving the coarse-grained representation, less demanding regarding computational resources, and thus offering much larger simulation timescales.
Collapse
Affiliation(s)
- Paulina Dragan
- Faculty of Chemistry, University of Warsaw, Warsaw, Poland
| | | | | | - Dorota Latek
- Faculty of Chemistry, University of Warsaw, Warsaw, Poland.
| |
Collapse
|
41
|
Palladium(II) Complexes of Substituted Salicylaldehydes: Synthesis, Characterization and Investigation of Their Biological Profile. Pharmaceuticals (Basel) 2022; 15:ph15070886. [PMID: 35890184 PMCID: PMC9323974 DOI: 10.3390/ph15070886] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Revised: 07/06/2022] [Accepted: 07/14/2022] [Indexed: 02/04/2023] Open
Abstract
Five palladium(II) complexes of substituted salicylaldehydes (X-saloH, X = 4-Et2N (for 1), 3,5-diBr (for 2), 3,5-diCl (for 3), 5-F (for 4) or 4-OMe (for 5)) bearing the general formula [Pd(X-salo)2] were synthesized and structurally characterized. The crystal structure of complex [Pd(4-Et2N-salo)2] was determined by single-crystal X-ray crystallography. The complexes can scavenge 1,1-diphenyl-picrylhydrazyl and 2,2′-azinobis(3-ethylbenzothiazoline-6-sulfonic acid) radicals and reduce H2O2. They are active against two Gram-positive (Staphylococcus aureus and Bacillus subtilis) and two Gram-negative (Escherichia coli and Xanthomonas campestris) bacterial strains. The complexes interact strongly with calf-thymus DNA via intercalation, as deduced by diverse techniques and via the determination of their binding constants. Complexes interact reversibly with bovine and human serum albumin. Complementary insights into their possible mechanisms of bioactivity at the molecular level were provided by molecular docking calculations, exploring in silico their ability to bind to calf-thymus DNA, Escherichia coli and Staphylococcus aureus DNA-gyrase, 5-lipoxygenase, and membrane transport lipid protein 5-lipoxygenase-activating protein, contributing to the understanding of the role complexes 1–5 can play both as antioxidant and antibacterial agents. Furthermore, in silico predictive tools have been employed to study the chemical reactivity, molecular properties and drug-likeness of the complexes, and also the drug-induced changes of gene expression profile (as protein- and mRNA-based prediction results), the sites of metabolism, the substrate/metabolite specificity, the cytotoxicity for cancer and non-cancer cell lines, the acute rat toxicity, the rodent organ-specific carcinogenicity, the anti-target interaction profiles, the environmental ecotoxicity, and finally the activity spectra profile of the compounds.
Collapse
|
42
|
Gao K, Wang R, Chen J, Cheng L, Frishcosy J, Huzumi Y, Qiu Y, Schluckbier T, Wei X, Wei GW. Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2. Chem Rev 2022; 122:11287-11368. [PMID: 35594413 PMCID: PMC9159519 DOI: 10.1021/acs.chemrev.1c00965] [Citation(s) in RCA: 31] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.
Collapse
Affiliation(s)
- Kaifu Gao
- Department
of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Rui Wang
- Department
of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Jiahui Chen
- Department
of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Limei Cheng
- Clinical
Pharmacology and Pharmacometrics, Bristol
Myers Squibb, Princeton, New Jersey 08536, United States
| | - Jaclyn Frishcosy
- Department
of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Yuta Huzumi
- Department
of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Yuchi Qiu
- Department
of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Tom Schluckbier
- Department
of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Xiaoqi Wei
- Department
of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Guo-Wei Wei
- Department
of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
- Department
of Electrical and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, United States
- Department
of Biochemistry and Molecular Biology, Michigan
State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
43
|
Cong Y, Li M, Qi Y, Zhang JZH. A fast-slow method to treat solute dynamics in explicit solvent. Phys Chem Chem Phys 2022; 24:14498-14510. [PMID: 35665790 DOI: 10.1039/d2cp00732k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Aiming to reduce the computational cost in the current explicit solvent molecular dynamics (MD) simulation, this paper proposes a fast-slow method for the fast MD simulation of biomolecules in explicit solvent. This fast-slow method divides the entire system into two parts: a core layer (typically solute or biomolecule) and a peripheral layer (typically solvent molecules). The core layer is treated using standard MD method but the peripheral layer is treated by a slower dynamics method to reduce the computational cost. We compared four different simulation models in testing calculations for several small proteins. These include gas-phase, implicit solvent, fast-slow explicit solvent and standard explicit solvent MD simulations. Our study shows that gas-phase and implicit solvent models do not provide a realistic solvent environment and fail to correctly produce reliable dynamic structures of proteins. On the other hand, the fast-slow method can essentially reproduce the same solvent effect as the standard explicit solvent model while gaining an order of magnitude in efficiency. This fast-slow method thus provides an efficient approach for accelerating the MD simulation of biomolecules in explicit solvent.
Collapse
Affiliation(s)
- Yalong Cong
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University at Shanghai, 200062, China.
| | - Mengxin Li
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University at Shanghai, 200062, China.
| | - Yifei Qi
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University at Shanghai, 200062, China.
| | - John Z H Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University at Shanghai, 200062, China. .,CAS Key Laboratory of Quantitative Engineering Biology, Shenzhen Institute of Synthetic Biology, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, Guangdong, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Department of Chemistry, New York University, NY, NY 10003, USA.,Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi, 030006, China
| |
Collapse
|
44
|
Wang A, Levi M, Mohanty U, Whitford PC. Diffuse Ions Coordinate Dynamics in a Ribonucleoprotein Assembly. J Am Chem Soc 2022; 144:9510-9522. [PMID: 35593477 DOI: 10.1021/jacs.2c04082] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Proper ionic concentrations are required for the functional dynamics of RNA and ribonucleoprotein (RNP) assemblies. While experimental and computational techniques have provided many insights into the properties of chelated ions, less is known about the energetic contributions of diffuse ions to large-scale conformational rearrangements. To address this, we present a model that is designed to quantify the influence of diffuse monovalent and divalent ions on the dynamics of biomolecular assemblies. This model employs all-atom (non-H) resolution and explicit ions, where effective potentials account for hydration effects. We first show that the model accurately predicts the number of excess Mg2+ ions for prototypical RNA systems, at a level comparable to modern coarse-grained models. We then apply the model to a complete ribosome and show how the balance between diffuse Mg2+ and K+ ions can control the dynamics of tRNA molecules during translation. The model predicts differential effects of diffuse ions on the free-energy barrier associated with tRNA entry and the energy of tRNA binding to the ribosome. Together, this analysis reveals the direct impact of diffuse ions on the dynamics of an RNP assembly.
Collapse
Affiliation(s)
- Ailun Wang
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States.,Center for Theoretical Biological Physics, Northeastern University, 360 Huntington Avenue, Boston, Massachusetts 02115, United States
| | - Mariana Levi
- Department of Physics, Northeastern University, Dana Research Center 111, 360 Huntington Avenue, Boston, Massachusetts 02115, United States
| | - Udayan Mohanty
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Paul C Whitford
- Center for Theoretical Biological Physics, Northeastern University, 360 Huntington Avenue, Boston, Massachusetts 02115, United States.,Department of Physics, Northeastern University, Dana Research Center 111, 360 Huntington Avenue, Boston, Massachusetts 02115, United States
| |
Collapse
|
45
|
Dobler L, Oliveira RR. Automated Search For The Low-lying Energy Isomers of Rhamnolipids and Related Organometallic Complexes. Chemphyschem 2022; 23:e202200111. [PMID: 35588462 DOI: 10.1002/cphc.202200111] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 05/19/2022] [Indexed: 11/10/2022]
Abstract
Rhamnolipids (RMLs) are a widely studied biosurfactant due to their high biodegradability and environmentally friendly pro duction. However, the knowledge of the structure-property relationship of RMLs is imperative for the design of highly efficient applications. Aiming to a better understanding of it at a molecular level, we performed an automated search for low energy structures of the most abundant RMLs, namely, Rha-C 10 , Rha-C 10 -C 10 , Rha-Rha-C 10 and Rha-Rha-C 10 -C 10 and their respective C 2 -congeners. Besides that, selected neutral metal complexes were also considered. We found that several low-energy congeners have internal hydrogen bonds. Moreover, geometries in "closed" conformation were always more stable than "open" ones. Fi nally, the energy differences between open and closed conformations of K + , Ni 2 + , Cu 2 + and Zn 2 + complexes were found to be 23.5 kcal mol -1 , 62.8 kcal mol -1 , 24.3 kcal mol -1 and 41.6 kcal mol -1 , respectively, indicating a huge structural reorganization after the complex formation.
Collapse
Affiliation(s)
- Leticia Dobler
- Universidade Federal do Rio de Janeiro, Chemistry Institute, BRAZIL
| | - Ricardo Rodrigues Oliveira
- Universidade Federal do Rio de Janeiro, Physical Chemistry, Av. Athos da Silveira Ramos, Technological Center, Block A, 304, University City, 21941-590, Rio de Janeiro, BRAZIL
| |
Collapse
|
46
|
Panday S, Alexov E. Protein-Protein Binding Free Energy Predictions with the MM/PBSA Approach Complemented with the Gaussian-Based Method for Entropy Estimation. ACS OMEGA 2022; 7:11057-11067. [PMID: 35415339 PMCID: PMC8991903 DOI: 10.1021/acsomega.1c07037] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Accepted: 03/10/2022] [Indexed: 06/14/2023]
Abstract
Here, we present a Gaussian-based method for estimation of protein-protein binding entropy to augment the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method for computational prediction of binding free energy (ΔG). The method is termed f5-MM/PBSA/E, where "E" stands for entropy and f5 for five adjustable parameters. The enthalpy components of ΔG (molecular mechanics, polar and non-polar solvation energies) are computed from a single implicit solvent generalized Born (GB) energy minimized structure of a protein-protein complex, while the binding entropy is computed using independently GB energy minimized unbound and bound structures. It should be emphasized that the f5-MM/PBSA/E method does not use snapshots, just energy minimized structures, and is thus very fast and computationally efficient. The method is trained and benchmarked in 5-fold validation test over a data set consisting of 46 protein-protein binding cases with experimentally determined dissociation constant K d values. This data set has been used for benchmarking in recently published protein-protein binding studies that apply conventional MM/PBSA and MM/PBSA with an enhanced sampling method. The f5-MM/PBSA/E tested on the same data set achieves similar or better performance than these computationally demanding approaches, making it an excellent choice for high throughput protein-protein binding affinity prediction studies.
Collapse
|
47
|
Wang R, Zheng Q. Multiple Molecular Dynamics Simulations and Free-Energy Predictions Uncover the Susceptibility of Variants of HIV-1 Protease against Inhibitors Darunavir and KNI-1657. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:14407-14418. [PMID: 34851643 DOI: 10.1021/acs.langmuir.1c02348] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
HIV-1 protease (PR) is considered to be the main targets of anti-AIDS drug design because of its role in the proteolytic processing of viral polyproteins. However, the emergence of drug-resistant HIV has become a major problem in the therapy of HIV-1-infected patients. Focused on the complexes of wild type (WT) PR and two mutant PRs (V32I/L33F/I54M/V82I and V32I/L33F/I54M/I84 V) with inhibitors Darunavir (DRV) and KNI-1657 (KNI), respectively, we have conducted research on the conformational dynamics and the resistance mechanism caused by residue mutations through multiple molecular dynamics (MD) simulations combined with an energy (MM-PBSA and solvated interaction energy (SIE)) prediction. The results indicate that mutated residues of PR alter the distance between flap regions and catalytic sites, the volume of the inner catalytic site, and the curling degree of the flap tips, thereby affecting DRV and KNI inhibitor binding to PR. These mutated residues reduced the binding affinity of the two mutant PRs to DRV, resulting in drug resistance, whereas the two mutant PRs increase the binding affinity with KNI, indicating they enhance the sensitivity to KNI. Compared with the WT PR, the changes in van der Waals interaction and electrostatic interaction in the two variant PRs play a vital part in the binding of PR with DRV and KNI. These results may supply valuable guidance for the design of anti-AIDS drugs targeting PR.
Collapse
|
48
|
Discovery of inhibitors targeting protein tyrosine phosphatase 1B using a combined virtual screening approach. Mol Divers 2021; 26:2159-2174. [PMID: 34655403 DOI: 10.1007/s11030-021-10323-2] [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/14/2021] [Accepted: 09/21/2021] [Indexed: 10/20/2022]
Abstract
Protein tyrosine phosphatase 1B (PTP1B) acts as a therapeutic target for type 2 diabetes. However, the major challenges of PTP1B drug discovery are the poor selectivity and the weak oral bioavailability. In this study, we performed a combined virtual screening approach including multicomplex pharmacophore, molecular docking-based screening, van der Waals energy normalization, pose scaling factor, ADMET evaluation, and molecular dynamics simulation to select PTP1B inhibitors from three databases (PubChem, ChEMBL, and ZINC). We identified three potential PTP1B inhibitors, compounds 1, 4, and 5, with favorable binding energy and good oral bioavailability. The energetic and geometrical analyses show that the three compounds are stably bound to PTP1B, via occupying both the catalytic site (site A) and the proximal noncatalytic site (site B or C). Such occupancy may improve the selectivity. This work not only provided a feasible virtual screening protocol, but also suggested three potential PTP1B inhibitors for the treatment of type 2 diabetes.
Collapse
|
49
|
Vassetti D, Oǧuz IC, Labat F. Generalizing Continuum Solvation in Crystal to Nonaqueous Solvents: Implementation, Parametrization, and Application to Molecules and Surfaces. J Chem Theory Comput 2021; 17:6432-6448. [PMID: 34488338 DOI: 10.1021/acs.jctc.1c00611] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present an extension of a generalized finite-difference Poisson-Boltzmann (FDPB) continuum solvation model based on a self-consistent reaction field treatment to nonaqueous solvents. Implementation and reparametrization of the cavitation, dispersion, and structural (CDS) effects nonelectrostatic model are presented in CRYSTAL, with applications to both finite and infinite periodic systems. For neutral finite systems, computed errors with respect to available experimental data on free energies of solvation of 2523 solutes in 91 solvents, as well as 144 transfer energies from water to 14 organic solvents are on par with the reference SM12 solvation model for which the CDS parameters have been developed. Calculations performed on a TiO2 anatase surface and compared to VASPsol data revealed an overall very good agreement of computed solvation energies, surface energies, as well as band structure changes upon solvation in three different solvents, validating the general applicability of the reparametrized FDPB approach to neutral nonperiodic and periodic solutes in aqueous and nonaqueous solvents. For ionic species, while the reparametrized CDS model led to large errors on free energies of solvation of anions, addition of a corrective term based on Abraham's acidity of the solvent significantly improved the accuracy of the proposed continuum solvation model, leading to errors on aqueous pKa of a test set of 83 solutes divided by a factor of 4 compared to the reference solvation model based on density (SMD). Overall, therefore, these encouraging results demonstrate that the generalized FDPB continuum solvation model can be applied to a broad range of solutes in various solvents, ranging from finite neutral or charged solutes to extended periodic surfaces.
Collapse
Affiliation(s)
- Dario Vassetti
- Chemical Theory and Modelling Group, Institute of Chemistry for Life and Health Sciences, Chimie ParisTech, PSL University, CNRS, F-75005 Paris, France
| | - Ismail Can Oǧuz
- Chemical Theory and Modelling Group, Institute of Chemistry for Life and Health Sciences, Chimie ParisTech, PSL University, CNRS, F-75005 Paris, France
| | - Frédéric Labat
- Chemical Theory and Modelling Group, Institute of Chemistry for Life and Health Sciences, Chimie ParisTech, PSL University, CNRS, F-75005 Paris, France
| |
Collapse
|
50
|
Zhou S, Zhang Y, Cheng LT, Li B. Prediction of multiple dry-wet transition pathways with a mesoscale variational approach. J Chem Phys 2021; 155:124110. [PMID: 34598586 DOI: 10.1063/5.0061773] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Water fluctuates in a hydrophobic confinement, forming multiple dry and wet hydration states through evaporation and condensation. Transitions between such states are critical to both thermodynamics and kinetics of solute molecular processes, such as protein folding and protein-ligand binding and unbinding. To efficiently predict such dry-wet transition paths, we develop a hybrid approach that combines a variational implicit solvation model, a generalized string method for minimum free-energy paths, and the level-set numerical implementation. This approach is applied to three molecular systems: two hydrophobic plates, a carbon nanotube, and a synthetic host molecule Cucurbit[7]uril. Without an explicit description of individual water molecules, our mesoscale approach effectively captures multiple dry and wet hydration states, multiple dry-wet transition paths, such as those geometrically symmetric and asymmetric paths, and transition states, providing activation energy barriers between different states. Further analysis shows that energy barriers depend on mesoscopic lengths, such as the separation distance between the two plates and the cross section diameter of the nanotube, and that the electrostatic interactions strongly influence the dry-wet transitions. With the inclusion of solute atomic motion, general collective variables as reaction coordinates, and the finite-temperature string method, together with an improved treatment of continuum electrostatics, our approach can be further developed to sample an ensemble of transition paths, providing more accurate predictions of the transition kinetics.
Collapse
Affiliation(s)
- Shenggao Zhou
- School of Mathematical Sciences and MOE-LSC, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Yanan Zhang
- School of Mathematical Sciences, Soochow University, Suzhou 215006, China
| | - Li-Tien Cheng
- Department of Mathematics, University of California San Diego, La Jolla, California 92093-0112, USA
| | - Bo Li
- Department of Mathematics, University of California San Diego, La Jolla, California 92093-0112, USA
| |
Collapse
|