In silico substrate-binding profiling for SARS-CoV-2 Main protease (Mpro) using Hexapeptide substrates
- Zabo, Sophakama, Lobb, Kevin A
- Authors: Zabo, Sophakama , Lobb, Kevin A
- Date: 2023
- Subjects: To be catalogued
- Language: English
- Type: text , article
- Identifier: http://hdl.handle.net/10962/452711 , vital:75164 , xlink:href="https://doi.org/10.3390/v15071480"
- Description: The SARS-CoV-2 main protease (Mpro) is essential for the life cycle of the COVID-19 virus. It cleaves the two polyproteins at 11 positions to generate mature proteins for virion formation. The cleavage site on these polyproteins is known to be Leu-Gln↓(Ser/Ala/Gly). A range of hexapeptides that follow the known sequence for recognition and cleavage was constructed using RDKit libraries and complexed with the crystal structure of Mpro (PDB ID 6XHM) through extensive molecular docking calculations. A subset of 131 of these complexes underwent 20 ns molecular dynamics simulations. The analyses of the trajectories from molecular dynamics included principal component analysis (PCA), and a method to compare PCA plots from separate trajectories was developed in terms of encoding PCA progression during the simulations. The hexapeptides formed stable complexes as expected, with reproducible molecular docking of the substrates given the extensiveness of the procedure. Only Lys-Leu-Gln*** (KLQ***) sequence complexes were studied for molecular dynamics. In this subset of complexes, the PCA analysis identified four classifications of protein motions across these sequences. KLQ*** complexes illustrated the effect of changes in substrate on the active site, with implications for understanding the substrate recognition of Mpro and informing the development of small molecule inhibitors.
- Full Text:
- Date Issued: 2023
- Authors: Zabo, Sophakama , Lobb, Kevin A
- Date: 2023
- Subjects: To be catalogued
- Language: English
- Type: text , article
- Identifier: http://hdl.handle.net/10962/452711 , vital:75164 , xlink:href="https://doi.org/10.3390/v15071480"
- Description: The SARS-CoV-2 main protease (Mpro) is essential for the life cycle of the COVID-19 virus. It cleaves the two polyproteins at 11 positions to generate mature proteins for virion formation. The cleavage site on these polyproteins is known to be Leu-Gln↓(Ser/Ala/Gly). A range of hexapeptides that follow the known sequence for recognition and cleavage was constructed using RDKit libraries and complexed with the crystal structure of Mpro (PDB ID 6XHM) through extensive molecular docking calculations. A subset of 131 of these complexes underwent 20 ns molecular dynamics simulations. The analyses of the trajectories from molecular dynamics included principal component analysis (PCA), and a method to compare PCA plots from separate trajectories was developed in terms of encoding PCA progression during the simulations. The hexapeptides formed stable complexes as expected, with reproducible molecular docking of the substrates given the extensiveness of the procedure. Only Lys-Leu-Gln*** (KLQ***) sequence complexes were studied for molecular dynamics. In this subset of complexes, the PCA analysis identified four classifications of protein motions across these sequences. KLQ*** complexes illustrated the effect of changes in substrate on the active site, with implications for understanding the substrate recognition of Mpro and informing the development of small molecule inhibitors.
- Full Text:
- Date Issued: 2023
Insights into the Dynamics and Binding of Two Polyprotein Substrate Cleavage Points in the Context of the SARS-CoV-2 Main and Papain-like Proteases
- Sanusi, Zainab K, Lobb, Kevin A
- Authors: Sanusi, Zainab K , Lobb, Kevin A
- Date: 2022
- Subjects: To be catalogued
- Language: English
- Type: text , article
- Identifier: http://hdl.handle.net/10962/452773 , vital:75169 , xlink:href="https://doi.org/10.3390/molecules27238251"
- Description: It is well known that vital enzymes in the replication process of the coronavirus are the SARS-CoV-2 PLpro and SARS-CoV-2 3CLpro, both of which are important targets in the search for anti-coronavirus agents. These two enzymes are responsible for cleavage at various polyprotein sites in the SARS-CoV-2 lifecycle. Herein, the dynamics of the polyprotein cleavage sequences for the boundary between non-structural proteins Nsp1 and Nsp2 (CS1) and between Nsp2 and Nsp3 (CS2) in complex with both the papain-like protein PLpro and the main protease 3CLpro were explored using computational methods. The post dynamics analysis reveals that CS1 and CS2 both have greater stability when complexed with PLpro. Of these two, greater stability is observed for the CS1–PLpro complex, while destabilization resulting in loss of CS2 from the PLpro active site is observed for CS2-PLpro, suggesting the rate of exchange by the papain-like protease is faster for CS2 compared to CS1. On the other hand, the 3CLpro main protease also reveals stability for CS1 suggesting that the main protease could also play a potential role in the cleavage at point CS1. However, destabilization occurs early in the simulation for the complex CLpro–CS2 suggesting a poor interaction and non-plausible protease cleavage of the polyprotein at CS2 by the main protease. These findings could be used as a guide in the development and design of potent COVID-19 antiviral inhibitors that mimic the CS1 cleavage site.
- Full Text:
- Date Issued: 2022
- Authors: Sanusi, Zainab K , Lobb, Kevin A
- Date: 2022
- Subjects: To be catalogued
- Language: English
- Type: text , article
- Identifier: http://hdl.handle.net/10962/452773 , vital:75169 , xlink:href="https://doi.org/10.3390/molecules27238251"
- Description: It is well known that vital enzymes in the replication process of the coronavirus are the SARS-CoV-2 PLpro and SARS-CoV-2 3CLpro, both of which are important targets in the search for anti-coronavirus agents. These two enzymes are responsible for cleavage at various polyprotein sites in the SARS-CoV-2 lifecycle. Herein, the dynamics of the polyprotein cleavage sequences for the boundary between non-structural proteins Nsp1 and Nsp2 (CS1) and between Nsp2 and Nsp3 (CS2) in complex with both the papain-like protein PLpro and the main protease 3CLpro were explored using computational methods. The post dynamics analysis reveals that CS1 and CS2 both have greater stability when complexed with PLpro. Of these two, greater stability is observed for the CS1–PLpro complex, while destabilization resulting in loss of CS2 from the PLpro active site is observed for CS2-PLpro, suggesting the rate of exchange by the papain-like protease is faster for CS2 compared to CS1. On the other hand, the 3CLpro main protease also reveals stability for CS1 suggesting that the main protease could also play a potential role in the cleavage at point CS1. However, destabilization occurs early in the simulation for the complex CLpro–CS2 suggesting a poor interaction and non-plausible protease cleavage of the polyprotein at CS2 by the main protease. These findings could be used as a guide in the development and design of potent COVID-19 antiviral inhibitors that mimic the CS1 cleavage site.
- Full Text:
- Date Issued: 2022
Force Field Parameters for Fe2+ 4S2− 4 Clusters of Dihydropyrimidine Dehydrogenase, the 5-Fluorouracil Cancer Drug Deactivation Protein: A Step towards In Silico Pharmacogenomics Studies
- Tendwa, Maureen B, Chebon-Bore, Lorna, Lobb, Kevin A, Musyoka, Thommas M, Taştan Bishop, Özlem
- Authors: Tendwa, Maureen B , Chebon-Bore, Lorna , Lobb, Kevin A , Musyoka, Thommas M , Taştan Bishop, Özlem
- Date: 2021
- Subjects: To be catalogued
- Language: English
- Type: text , article
- Identifier: http://hdl.handle.net/10962/451078 , vital:75016 , xlink:href="https://doi.org/10.3390/molecules26102929 "
- Description: The dimeric dihydropyrimidine dehydrogenase (DPD), metalloenzyme, an adjunct anti-cancer drug target, contains highly specialized 4 × Fe2+4S2−4 clusters per chain. These clusters facilitate the catalysis of the rate-limiting step in the pyrimidine degradation pathway through a harmonized electron transfer cascade that triggers a redox catabolic reaction. In the process, the bulk of the administered 5-fluorouracil (5-FU) cancer drug is inactivated, while a small proportion is activated to nucleic acid antimetabolites. The occurrence of missense mutations in DPD protein within the general population, including those of African descent, has adverse toxicity effects due to altered 5-FU metabolism. Thus, deciphering mutation effects on protein structure and function is vital, especially for precision medicine purposes. We previously proposed combining molecular dynamics (MD) and dynamic residue network (DRN) analysis to decipher the molecular mechanisms of missense mutations in other proteins. However, the presence of Fe2+4S2−4 clusters in DPD poses a challenge for such in silico studies. The existing AMBER force field parameters cannot accurately describe the Fe2+ center coordination exhibited by this enzyme. Therefore, this study aimed to derive AMBER force field parameters for DPD enzyme Fe2+ centers, using the original Seminario method and the collation features Visual Force Field Derivation Toolkit as a supportive approach. All-atom MD simulations were performed to validate the results. Both approaches generated similar force field parameters, which accurately described the human DPD protein Fe2+4S2−4 cluster architecture. This information is crucial and opens new avenues for in silico cancer pharmacogenomics and drug discovery related research on 5-FU drug efficacy and toxicity issues.
- Full Text:
- Date Issued: 2021
- Authors: Tendwa, Maureen B , Chebon-Bore, Lorna , Lobb, Kevin A , Musyoka, Thommas M , Taştan Bishop, Özlem
- Date: 2021
- Subjects: To be catalogued
- Language: English
- Type: text , article
- Identifier: http://hdl.handle.net/10962/451078 , vital:75016 , xlink:href="https://doi.org/10.3390/molecules26102929 "
- Description: The dimeric dihydropyrimidine dehydrogenase (DPD), metalloenzyme, an adjunct anti-cancer drug target, contains highly specialized 4 × Fe2+4S2−4 clusters per chain. These clusters facilitate the catalysis of the rate-limiting step in the pyrimidine degradation pathway through a harmonized electron transfer cascade that triggers a redox catabolic reaction. In the process, the bulk of the administered 5-fluorouracil (5-FU) cancer drug is inactivated, while a small proportion is activated to nucleic acid antimetabolites. The occurrence of missense mutations in DPD protein within the general population, including those of African descent, has adverse toxicity effects due to altered 5-FU metabolism. Thus, deciphering mutation effects on protein structure and function is vital, especially for precision medicine purposes. We previously proposed combining molecular dynamics (MD) and dynamic residue network (DRN) analysis to decipher the molecular mechanisms of missense mutations in other proteins. However, the presence of Fe2+4S2−4 clusters in DPD poses a challenge for such in silico studies. The existing AMBER force field parameters cannot accurately describe the Fe2+ center coordination exhibited by this enzyme. Therefore, this study aimed to derive AMBER force field parameters for DPD enzyme Fe2+ centers, using the original Seminario method and the collation features Visual Force Field Derivation Toolkit as a supportive approach. All-atom MD simulations were performed to validate the results. Both approaches generated similar force field parameters, which accurately described the human DPD protein Fe2+4S2−4 cluster architecture. This information is crucial and opens new avenues for in silico cancer pharmacogenomics and drug discovery related research on 5-FU drug efficacy and toxicity issues.
- Full Text:
- Date Issued: 2021
The Formation of 2, 2, 4-Trimethyl-2, 3-dihydro-1 H-1, 5-Benzodiazepine from 1, 2-Diaminobenzene in the Presence of Acetone
- Odame, Felix, Kleyi, Phumelele, Hosten, Eric C, Betz, Richard, Lobb, Kevin A, Tshentu, Zenixole R
- Authors: Odame, Felix , Kleyi, Phumelele , Hosten, Eric C , Betz, Richard , Lobb, Kevin A , Tshentu, Zenixole R
- Date: 2013
- Subjects: To be catalogued
- Language: English
- Type: text , article
- Identifier: http://hdl.handle.net/10962/448952 , vital:74773 , xlink:href="https://doi.org/10.3390/molecules181114293"
- Description: In an attempt to synthesize a 2-substituted benzimidazole from the reaction of o-phenylenediamine and isophthalic acid in the presence of acetone and ethanol under microwave irradiation, a salt of the isophthalate ion and 2,2,4-trimethyl-2,3-dihydro-1H-1,5-benzodiazepin-5-ium ion was obtained. The condensation of two moles of acetone with the amine groups resulted in the formation of the benzodiazepine which crystallized as an iminium cation forming a salt with the isophthalate anion. The formation of benzodiazepine was also confirmed by performing the reaction of o-phenylenediamine with excess acetone in ethanol under conventional heating conditions. The compounds were characterized by NMR, FTIR, HRMS and microanalysis as well as X-ray crystallography. The reaction mechanism leading to the formation of benzodiazepine is also discussed.
- Full Text:
- Date Issued: 2013
- Authors: Odame, Felix , Kleyi, Phumelele , Hosten, Eric C , Betz, Richard , Lobb, Kevin A , Tshentu, Zenixole R
- Date: 2013
- Subjects: To be catalogued
- Language: English
- Type: text , article
- Identifier: http://hdl.handle.net/10962/448952 , vital:74773 , xlink:href="https://doi.org/10.3390/molecules181114293"
- Description: In an attempt to synthesize a 2-substituted benzimidazole from the reaction of o-phenylenediamine and isophthalic acid in the presence of acetone and ethanol under microwave irradiation, a salt of the isophthalate ion and 2,2,4-trimethyl-2,3-dihydro-1H-1,5-benzodiazepin-5-ium ion was obtained. The condensation of two moles of acetone with the amine groups resulted in the formation of the benzodiazepine which crystallized as an iminium cation forming a salt with the isophthalate anion. The formation of benzodiazepine was also confirmed by performing the reaction of o-phenylenediamine with excess acetone in ethanol under conventional heating conditions. The compounds were characterized by NMR, FTIR, HRMS and microanalysis as well as X-ray crystallography. The reaction mechanism leading to the formation of benzodiazepine is also discussed.
- Full Text:
- Date Issued: 2013
- «
- ‹
- 1
- ›
- »