scholarly journals Structural Differences In 3C-like protease (Mpro) From SARS-CoV and SARS-CoV-2: Molecular Insights For Drug Repurposing Against COVID-19 Revealed by Molecular Dynamics Simulations.

2021 ◽  
Author(s):  
Dhaval Patel ◽  
Meet Parmar ◽  
Ritik Thumar ◽  
Bhumi Patel ◽  
Mohd. Athar ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Binding Site ◽  
Active Site ◽  
Substrate Binding ◽  
Drug Repurposing ◽  
Molecular Networks ◽  
Substrate Binding Site ◽  
Sequence Identity ◽  
Dynamics Simulations

A recent fatal outbreak of novel coronavirus SARS-CoV-2, identified preliminary as a causative agent for series of unusual pneumonia cases in Wuhan city, China has infected more than 20 million individuals with more than 4 million mortalities. Since, the infection crossed geographical barriers, the WHO permanently named the causing disease as COVID-2019 by declaring it a pandemic situation. SARS-CoV-2 is an enveloped single-stranded RNA virus causing a wide range of pathological conditions from common cold symptoms to pneumonia and fatal severe respiratory syndrome. Genome sequencing of SARS-CoV-2 has revealed 96% identity to the bat coronavirus and 79.6% sequence identity to the previous SARS-CoV. The main protease (known as 3C-like proteinase/ Mpro) plays a vital role during the infection with the processing of replicase polyprotein thus offering an attractive target for therapeutic interventions. SARS-CoV and SARS-CoV-2 Mpro shares 97% sequence identity, with 12 variable residues but none of them are present in the catalytic and substrate binding site. With the high level of sequence and structural similarity and absence of any drug/vaccine against SARS-CoV-2, drug repurposing against Mpro is an effective strategy to combat COVID-19. Here, we report a detailed comparison of SARS-CoV-2 Mpro with SARS-CoV Mpro using molecular dynamics simulations to assess the impact of 12 divergent residues on the molecular microenvironment of Mpro. Structural comparison and analysis are made on how these variable residues affect the intra-molecular interactions between key residues in the monomer and biologically active dimer form of Mpro. The present MD simulations study concluded the change in microenvironment of active-site residues at the entrance (T25, T26, M49 and Q189), near the catalytic region (F140, H163, H164, M165 and H172) and other residues in substrate binding site (V35T, N65S, K88R and N180K) due to 12 mutation incorporated in the SARS-CoV-2 Mpro. It is also evident that SARS-CoV-2 dimer is more stable and less flexible state compared to monomer which may be due to these variable residues, mainly F140, E166 and H172 which are involved in dimerization. This also warrants a need for inhibitor design considering the more stable dimer form. The mutation accumulated in SARS-CoV-2 Mpro indirectly reconfigures the key molecular networks around the active site conferring a potential change in SARS-CoV-2, thus posing a challenge in drug repurposing SARS drugs for COVID-19. The new networks and changes in the microenvironment identified by our work might guide attempts needed for repurposing and identification of new Mpro inhibitors.

scholarly journals Molecular Dynamics Simulations of A27S and K120A Mutated PTP1B Reveals Selective Binding of the Bidentate Inhibitor

2019 ◽  
Vol 2019 ◽  
pp. 1-11
Author(s):  
Xi Chen ◽  
Xia Liu ◽  
Qiang Gan ◽  
Changgen Feng ◽  
Qian Zhang
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Binding Site ◽  
Protein Tyrosine Phosphatase ◽  
Active Site ◽  
Tyrosine Phosphatase ◽  
Selective Binding ◽  
Protein Tyrosine ◽  
Ptp1b Inhibitors ◽  
Dynamics Simulations

Protein tyrosine phosphatase 1B (PTP1B) is considered a potential target for the treatment of type II diabetes and obesity due to its critical negative role in the insulin signaling pathway. However, improving the selectivity of PTP1B inhibitors over the most closely related T-cell protein tyrosine phosphatase (TCPTP) remains a major challenge for inhibitor development. Lys120 at the active site and Ser27 at the second pTyr binding site are distinct in PTP1B and TCPTP, which may bring differences in binding affinity. To explore the determinant of selective binding of inhibitor, molecular dynamics simulations with binding free energy calculations were performed on K120A and A27S mutated PTP1B, and the internal changes induced by mutations were investigated. Results reveal that the presence of Lys120 induces a conformational change in the WPD-loop and YRD-motif and has a certain effect on the selective binding at the active site. Ser27 weakens the stability of the inhibitor at the second pTyr binding site by altering the orientation of the Arg24 and Arg254 side chains via hydrogen bonds. Further comparison of alanine scanning demonstrates that the reduction in the energy contribution of Arg254 caused by A27S mutation leads to a different inhibitory activity. These observations provide novel insights into the selective binding mechanism of PTP1B inhibitors to TCPTP.

Design of Inhibitors for Glyceraldehyde-3-phosphate Dehydrogenase (GAPDH) Enzyme of Leishmania mexicana

2020 ◽  
Vol 16 (6) ◽  
pp. 784-795
Author(s):  
Krisnna M.A. Alves ◽  
Fábio José Bonfim Cardoso ◽  
Kathia M. Honorio ◽  
Fábio A. de Molfetta
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Binding Site ◽  
Point Of View ◽  
New Drugs ◽  
Leishmania Mexicana ◽  
Receptor Complexes ◽  
Starting Point ◽  
Dynamics Simulations ◽  
Selection Of

Background:: Leishmaniosis is a neglected tropical disease and glyceraldehyde 3- phosphate dehydrogenase (GAPDH) is a key enzyme in the design of new drugs to fight this disease. Objective:: The present study aimed to evaluate potential inhibitors of GAPDH enzyme found in Leishmania mexicana (L. mexicana). Methods: A search for novel antileishmanial molecules was carried out based on similarities from the pharmacophoric point of view related to the binding site of the crystallographic enzyme using the ZINCPharmer server. The molecules selected in this screening were subjected to molecular docking and molecular dynamics simulations. Results:: Consensual analysis of the docking energy values was performed, resulting in the selection of ten compounds. These ligand-receptor complexes were visually inspected in order to analyze the main interactions and subjected to toxicophoric evaluation, culminating in the selection of three compounds, which were subsequently submitted to molecular dynamics simulations. The docking results showed that the selected compounds interacted with GAPDH from L. mexicana, especially by hydrogen bonds with Cys166, Arg249, His194, Thr167, and Thr226. From the results obtained from molecular dynamics, it was observed that one of the loop regions, corresponding to the residues 195-222, can be related to the fitting of the substrate at the binding site, assisting in the positioning and the molecular recognition via residues responsible for the catalytic activity. Conclusion:: he use of molecular modeling techniques enabled the identification of promising compounds as inhibitors of the GAPDH enzyme from L. mexicana, and the results obtained here can serve as a starting point to design new and more effective compounds than those currently available.

scholarly journals The Interplay of Cholesterol and Ligand Binding in hTSPO from Classical Molecular Dynamics Simulations

Molecules ◽  
2021 ◽  
Vol 26 (5) ◽  
pp. 1250
Author(s):  
Hien T. T. Lai ◽  
Alejandro Giorgetti ◽  
Giulia Rossetti ◽  
Toan T. Nguyen ◽  
Paolo Carloni ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Ligand Binding ◽  
Molecular Dynamics Simulations ◽  
Binding Site ◽  
Transmembrane Protein ◽  
Free Form ◽  
Experimental Investigations ◽  
Terminal Part ◽  
Dynamics Simulations ◽  
Cholesterol Binding

The translocator protein (TSPO) is a 18kDa transmembrane protein, ubiquitously present in human mitochondria. It is overexpressed in tumor cells and at the sites of neuroinflammation, thus representing an important biomarker, as well as a promising drug target. In mammalian TSPO, there are cholesterol–binding motifs, as well as a binding cavity able to accommodate different chemical compounds. Given the lack of structural information for the human protein, we built a model of human (h) TSPO in the apo state and in complex with PK11195, a molecule routinely used in positron emission tomography (PET) for imaging of neuroinflammatory sites. To better understand the interactions of PK11195 and cholesterol with this pharmacologically relevant protein, we ran molecular dynamics simulations of the apo and holo proteins embedded in a model membrane. We found that: (i) PK11195 stabilizes hTSPO structural fold; (ii) PK11195 might enter in the binding site through transmembrane helices I and II of hTSPO; (iii) PK11195 reduces the frequency of cholesterol binding to the lower, N–terminal part of hTSPO in the inner membrane leaflet, while this impact is less pronounced for the upper, C–terminal part in the outer membrane leaflet, where the ligand binding site is located; (iv) very interestingly, cholesterol most frequently binds simultaneously to the so-called CRAC and CARC regions in TM V in the free form (residues L150–X–Y152–X(3)–R156 and R135–X(2)–Y138–X(2)–L141, respectively). However, when the protein is in complex with PK11195, cholesterol binds equally frequently to the CRAC–resembling motif that we observed in TM I (residues L17–X(2)–F20–X(3)–R24) and to CRAC in TM V. We expect that the CRAC–like motif in TM I will be of interest in future experimental investigations. Thus, our MD simulations provide insight into the structural features of hTSPO and the previously unknown interplay between PK11195 and cholesterol interactions with this pharmacologically relevant protein.

scholarly journals Exploiting correlated molecular-dynamics networks to counteract enzyme activity–stability trade-off

2018 ◽  
Vol 115 (52) ◽  
pp. E12192-E12200 ◽  
Author(s):  
Haoran Yu ◽  
Paul A. Dalby
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Active Site ◽  
Dynamic Networks ◽  
Aromatic Aldehydes ◽  
Active Site Residues ◽  
Trade Off ◽  
Highly Correlated ◽  
Dynamics Simulations ◽  
The Impact

The directed evolution of enzymes for improved activity or substrate specificity commonly leads to a trade-off in stability. We have identified an activity–stability trade-off and a loss in unfolding cooperativity for a variant (3M) of Escherichia coli transketolase (TK) engineered to accept aromatic substrates. Molecular dynamics simulations of 3M revealed increased flexibility in several interconnected active-site regions that also form part of the dimer interface. Mutating the newly flexible active-site residues to regain stability risked losing the new activity. We hypothesized that stabilizing mutations could be targeted to residues outside of the active site, whose dynamics were correlated with the newly flexible active-site residues. We previously stabilized WT TK by targeting mutations to highly flexible regions. These regions were much less flexible in 3M and would not have been selected a priori as targets using the same strategy based on flexibility alone. However, their dynamics were highly correlated with the newly flexible active-site regions of 3M. Introducing the previous mutations into 3M reestablished the WT level of stability and unfolding cooperativity, giving a 10.8-fold improved half-life at 55 °C, and increased midpoint and aggregation onset temperatures by 3 °C and 4.3 °C, respectively. Even the activity toward aromatic aldehydes increased up to threefold. Molecular dynamics simulations confirmed that the mutations rigidified the active-site via the correlated network. This work provides insights into the impact of rigidifying mutations within highly correlated dynamic networks that could also be useful for developing improved computational protein engineering strategies.

Active site gate of M32 carboxypeptidases illuminated by crystal structure and molecular dynamics simulations

2017 ◽  
Vol 1865 (11) ◽  
pp. 1406-1415 ◽  
Author(s):  
Bhaskar Sharma ◽  
Sahayog N. Jamdar ◽  
Biplab Ghosh ◽  
Pooja Yadav ◽  
Ashwani Kumar ◽  
...  
Keyword(s):  
Crystal Structure ◽  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Active Site ◽  
Dynamics Simulations

A molecular dynamics study of the complete binding process of meropenem to New Delhi metallo-β-lactamase 1

2018 ◽  
Vol 20 (9) ◽  
pp. 6409-6420 ◽  
Author(s):  
Juan Duan ◽  
Chuncai Hu ◽  
Jiafan Guo ◽  
Lianxian Guo ◽  
Jia Sun ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Substrate Binding ◽  
New Delhi ◽  
Binding Process ◽  
Dynamics Simulations

We have investigated the substrate-binding pathways of NDM-1 via unbiased molecular dynamics simulations and metadynamics.

scholarly journals IDENTIFICATION OF POTENTIAL SALMONELLA TYPHI BETA-LACTAMASE TEM 1 INHIBITORS USING PEPTIDOMIMETICS, VIRTUAL SCREENING, AND MOLECULAR DYNAMICS SIMULATIONS

2018 ◽  
Vol 10 (1) ◽  
pp. 91
Author(s):  
Rakesh K. R. Pandit ◽  
Dinesh Gupta ◽  
Tapan K. Mukherjee
Keyword(s):  
Molecular Dynamics ◽  
Antibiotic Resistance ◽  
Virtual Screening ◽  
Molecular Dynamics Simulations ◽  
Active Site ◽  
Beta Lactamase ◽  
Ramachandran Plot ◽  
Small Peptide ◽  
In Silico Docking ◽  
Dynamics Simulations

Objective: The purpose of this study was to identify a potential peptidomimetic S. typhi Beta-lactamase TEM 1 inhibitor to tackle the antibiotic resistance among S. typhi.Methods: The potential peptidomimetic inhibitor was identified by in silico docking of the small peptide WFRKQLKW with S. typhi Beta-lactamase TEM 1. The 3D coordinate geometry of the residues of small peptide interacting with the active site of the receptor was generated and mimics were identified using PEP: MMs: MIMIC server. All the identified mimics were docked at the active site of the receptor using Autodock 4.2 and the best-docked complex was selected on the basis of binding energy and number of H-bonds. The complex was then subjected to molecular dynamics simulations of 30 ns using AMBER 12 software package. The stereochemical stability of the Beta-lactamase TEM 1-WFRKQLKW complex was estimated with the help of Ramachandran plot using PROCHECK tool.Results: In the present study, a new potential peptidomimetic inhibitor (ZINC05839264) of Beta-lactamase TEM 1 has been identified based on antimicrobial peptide WFRKQLKW by virtual screening of the MMsINC database. The docking and molecular simulation studies revealed that the mimic binds more tightly to the active site of the receptor than the peptide. The Ramachandran plot also shows that the Beta-lactamase TEM 1-mimic complex is stereo chemically more stable than Beta-lactamase TEM 1-WFRKQLKW complex as more number of residues (93.6%) are falling under the core region of the plot in case of the former.Conclusion: The study shows that the peptidomimetic compound can act as a potential inhibitor of S. typhi Beta-lactamase TEM 1 and further it can be developed into more effective therapeutic to tackle the problem of antibiotic resistance.

Sign in / Sign up

Export Citation Format

Share Document