Exploring the Binding Mechanism of PF-07321332 SARS-CoV-2 Protease Inhibitor through Molecular Dynamics and Binding Free Energy Simulations

The novel coronavirus disease, caused by severe acute respiratory coronavirus 2 (SARS-CoV-2), rapidly spreading around the world, poses a major threat to the global public health. Herein, we demonstrated the binding mechanism of PF-07321332, α-ketoamide, lopinavir, and ritonavir to the coronavirus 3-chymotrypsin-like-protease (3CL^pro) by means of docking and molecular dynamic (MD) simulations. The analysis of MD trajectories of 3CL^pro with PF-07321332, α-ketoamide, lopinavir, and ritonavir revealed that 3CL^pro–PF-07321332 and 3CL^pro–α-ketoamide complexes remained stable compared with 3CL^pro–ritonavir and 3CL^pro–lopinavir. Investigating the dynamic behavior of ligand–protein interaction, ligands PF-07321332 and α-ketoamide showed stronger bonding via making interactions with catalytic dyad residues His41–Cys145 of 3CL^pro. Lopinavir and ritonavir were unable to disrupt the catalytic dyad, as illustrated by increased bond length during the MD simulation. To decipher the ligand binding mode and affinity, ligand interactions with SARS-CoV-2 proteases and binding energy were calculated. The binding energy of the bespoke antiviral PF-07321332 clinical candidate was two times higher than that of α-ketoamide and three times than that of lopinavir and ritonavir. Our study elucidated in detail the binding mechanism of the potent PF-07321332 to 3CL^pro along with the low potency of lopinavir and ritonavir due to weak binding affinity demonstrated by the binding energy data. This study will be helpful for the development and optimization of more specific compounds to combat coronavirus disease.     

Comparative Molecular Dynamics Simulations of Mitogen-Activated Protein Kinase-Activated Protein Kinase 5

The mitogen-activated protein kinase-activated protein kinase MK5 is a substrate of the mitogen-activated protein kinases p38, ERK3 and ERK4. Cell culture and animal studies have demonstrated that MK5 is involved in tumour suppression and promotion, embryogenesis, anxiety, cell motility and cell cycle regulation. In the present study, homology models of MK5 were used for molecular dynamics (MD) simulations of: (1) MK5 alone; (2) MK5 in complex with an inhibitor; and (3) MK5 in complex with the interaction partner p38α. The calculations showed that the inhibitor occupied the active site and disrupted the intramolecular network of amino acids. However, intramolecular interactions consistent with an inactive protein kinase fold were not formed. MD with p38α showed that not only the p38 docking region, but also amino acids in the activation segment, αH helix, P-loop, regulatory phosphorylation region and the C-terminal of MK5 may be involved in forming a very stable MK5-p38α complex, and that p38α binding decreases the residual fluctuation of the MK5 model. Electrostatic Potential Surface (EPS) calculations of MK5 and p38α showed that electrostatic interactions are important for recognition and binding.     

Combined Modeling Study of the Binding Characteristics of Natural Compounds,Derived from Psoralea Fruits,to β-Amyloid Peptide Monomer

A dysfunctional protein aggregation in the nervous system can lead to several neurodegenerative disorders that result in intracellular inclusions or extracellular aggregates. An early critical event within the pathogenesis of Alzheimer’s disease is the accumulation of amyloid beta peptide within the brain. Natural compounds isolated from Psoralea Fructus (PF) have significant anti-Alzheimer effects as strong inhibitors of Aβ42 aggregation. Computer simulations provide a powerful means of linking experimental findings to nanoscale molecular events. As part of this research four prenylated compounds, the active ingredients of Psoralea Fructus (PF), were studied as Aβ42 accumulation inhibitors using molecular simulations modeling. In order to resolve the binding modes of the ligands and identify the main interactions of Aβ42 residues, we performed a 100 ns molecular dynamics simulation and binding free energy calculations starting from the model of the compounds obtained from the docking study. This study was able to pinpoint the key amino acid residues in the Aβ42 active site and provide useful information that could benefit the development of new Aβ42 accumulation inhibitors.     

Potential Activity of Fevicordin-A from Phaleria macrocarpa (Scheff) Boerl. Seeds as Estrogen Receptor Antagonist Based on Cytotoxicity and Molecular Modelling Studies

Fevicordin-A (FevA) isolated from Phaleria macrocarpa (Scheff) Boerl. seeds was evaluated for its potential anticancer activity by in vitro and in silico approaches. Cytotoxicity studies indicated that FevA was selective against cell lines of human breast adenocarcinoma (MCF-7) with an IC₅₀ value of 6.4 μM. At 11.2 μM, FevA resulted in 76.8% cell death of T-47D human breast cancer cell lines. Critical pharmacophore features amongst human Estrogen Receptor-α (hERα) antagonists were conserved in FevA with regard to a hypothesis that they could make notable contributions to its pharmacological activity. The binding stability as well as the dynamic behavior of FevA towards the hERα receptor in agonist and antagonist binding sites were probed using molecular dynamics (MD) simulation approach. Analysis of MD simulation suggested that the tail of FevA was accountable for the repulsion of the C-terminal of Helix-11 (H11) in both agonist and antagonist receptor forms. The flexibility of loop-534 indicated the ability to disrupt the hydrogen bond zipper network between H3 and H11 in hERα. In addition, MM/GBSA calculation from the molecular dynamic simulations also revealed a stronger binding affinity of FevA in antagonistic action as compared to that of agonistic action. Collectively, both the experimental and computational results indicated that FevA has potential as a candidate for an anticancer agent, which is worth promoting for further preclinical evaluation.     

In Silico Screening of Novel α1-GABAA Receptor PAMs towards Schizophrenia Based on Combined Modeling Studies of Imidazo [1,2-a]-Pyridines

The ionotropic GABA_A receptor (GABA_AR) has been proven to be an important target of atypical antipsychotics. A novel series of imidazo [1,2-a]-pyridine derivatives, as selective positive allosteric modulators (PAMs) of α1-containing GABA_ARs with potent antipsychotic activities, have been reported recently. To better clarify the pharmacological essentiality of these PAMs and explore novel antipsychotics hits, three-dimensional quantitative structure–activity relationships (3D-QSAR), molecular docking, pharmacophore modeling, and molecular dynamics (MD) were performed on 33 imidazo [1,2-a]-pyridines. The constructed 3D-QSAR models exhibited good predictive abilities. The dockings results and MD simulations demonstrated that hydrogen bonds, π–π stackings, and hydrophobic interactions play essential roles in the binding of these novel PAMs in the GABA_AR binding pocket. Four hit compounds (DS01–04) were then screened out by the combination of the constructed models and computations, including the pharmacophore model, Topomer Search, molecular dockings, ADME/T predictions, and MD simulations. The compounds DS03 and DS04, with higher docking scores and better predicted activities, were also found to be relatively stable in the binding pocket by MD simulations. These results might provide a significant theoretical direction or information for the rational design and development of novel α1-GABA_AR PAMs with antipsychotic activities.     

Screening of β1- and β2-Adrenergic Receptor Modulators through Advanced Pharmacoinformatics and Machine Learning Approaches

Cardiovascular diseases (CDs) are a major concern in the human race and one of the leading causes of death worldwide. β-Adrenergic receptors (β1-AR and β2-AR) play a crucial role in the overall regulation of cardiac function. In the present study, structure-based virtual screening, machine learning (ML), and a ligand-based similarity search were conducted for the PubChem database against both β1- and β2-AR. Initially, all docked molecules were screened using the threshold binding energy value. Molecules with a better binding affinity were further used for segregation as active and inactive through ML. The pharmacokinetic assessment was carried out on molecules retained in the above step. Further, similarity searching of the ChEMBL and DrugBank databases was performed. From detailed analysis of the above data, four compounds for each of β1- and β2-AR were found to be promising in nature. A number of critical ligand-binding amino acids formed potential hydrogen bonds and hydrophobic interactions. Finally, a molecular dynamics (MD) simulation study of each molecule bound with the respective target was performed. A number of parameters obtained from the MD simulation trajectories were calculated and substantiated the stability between the protein-ligand complex. Hence, it can be postulated that the final molecules might be crucial for CDs subjected to experimental validation.     

Combinatorial Virtual Library Screening Study of Transforming Growth Factor-β2–Chondroitin Sulfate System

Transforming growth factor-beta (TGF-β), a member of the TGF-β cytokine superfamily, is known to bind to sulfated glycosaminoglycans (GAGs), but the nature of this interaction remains unclear. In a recent study, we found that preterm human milk TGF-β2 is sequestered by chondroitin sulfate (CS) in its proteoglycan form. To understand the molecular basis of the TGF-β2–CS interaction, we utilized the computational combinatorial virtual library screening (CVLS) approach in tandem with molecular dynamics (MD) simulations. All possible CS oligosaccharides were generated in a combinatorial manner to give 24 di- (CS02), 192 tetra- (CS04), and 1536 hexa- (CS06) saccharides. This library of 1752 CS oligosaccharides was first screened against TGF-β2 using the dual filter CVLS algorithm in which the GOLDScore and root-mean-square-difference (RMSD) between the best bound poses were used as surrogate markers for in silico affinity and in silico specificity. CVLS predicted that both the chain length and level of sulfation are critical for the high affinity and high specificity recognition of TGF-β2. Interestingly, CVLS led to identification of two distinct sites of GAG binding on TGF-β2. CVLS also deduced the preferred composition of the high specificity hexasaccharides, which were further assessed in all-atom explicit solvent MD simulations. The MD results confirmed that both sites of binding form stable GAG–protein complexes. More specifically, the highly selective CS chains were found to engage the TGF-β2 monomer with high affinity. Overall, this work present key principles of recognition with regard to the TGF-β2–CS system. In the process, it led to the generation of the in silico library of all possible CS oligosaccharides, which can be used for advanced studies on other protein–CS systems. Finally, the study led to the identification of unique CS sequences that are predicted to selectively recognize TGF-β2 and may out-compete common natural CS biopolymers.     

Molecular Dynamics Simulations Based on 1-Phenyl-4-Benzoyl-1-Hydro-Triazole ERRα Inverse Agonists

Estrogen-related receptor α (ERRα), which is overexpressed in a variety of cancers has been considered as an effective target for anticancer therapy. ERRα inverse agonists have been proven to effectively inhibit the migration and invasion of cancer cells. As few crystalline complexes have been reported, molecular dynamics (MD) simulations were carried out in this study to deepen the understanding of the interaction mechanism between inverse agonists and ERRα. The binding free energy was analyzed by the MM-GBSA method. The results show that the total binding free energy was positively correlated with the biological activity of an inverse agonist. The interaction of the inverse agonist with the hydrophobic interlayer composed of Phe328 and Phe495 had an important impact on the biological activity of inverse agonists, which was confirmed by the decomposition of energy on residues. As Glu331 flipped and formed a hydrogen bond with Arg372 in the MD simulation process, the formation of hydrogen bond interaction with Glu331 was not a necessary condition for the compound to act as an inverse agonist. These rules provide guidance for the design of new inverse agonists.     

Covalent Complex of DNA and Bacterial Topoisomerase: Implications in Antibacterial Drug Development

A topoisomerase-DNA transient covalent complex can be a druggable target for novel topoisomerase poison inhibitors that represent a new class of antibacterial or anticancer drugs. Herein, we have investigated molecular features of the functionally important Escherichia coli topoisomerase I (EctopoI)-DNA covalent complex (EctopoIcc) for molecular simulations, which is very useful in the development of new antibacterial drugs. To demonstrate the usefulness of our approach, we used a model small molecule (SM), NSC76027, obtained from virtual screening. We examined the direct binding of NSC76027 to EctopoI as well as inhibition of EctopoI relaxation activity of this SM via experimental techniques. We then performed molecular dynamics (MD) simulations to investigate the dynamics and stability of EctopoIcc and EctopoI-NSC76027-DNA ternary complex. Our simulation results show that NSC76027 forms a stable ternary complex with EctopoIcc. EctopoI investigated here also serves as a model system for investigating a complex of topoisomerase and DNA in which DNA is covalently attached to the protein.     

Docking and Molecular Dynamic of Microalgae Compounds as Potential Inhibitors of Beta-Lactamase

Bacterial resistance is responsible for a wide variety of health problems, both in children and adults. The persistence of symptoms and infections are mainly treated with β-lactam antibiotics. The increasing resistance to those antibiotics by bacterial pathogens generated the emergence of extended-spectrum β-lactamases (ESBLs), an actual public health problem. This is due to rapid mutations of bacteria when exposed to antibiotics. In this case, β-lactamases are enzymes used by bacteria to hydrolyze the beta-lactam rings present in the antibiotics. Therefore, it was necessary to explore novel molecules as potential β-lactamases inhibitors to find antibacterial compounds against infection caused by ESBLs. A computational methodology based on molecular docking and molecular dynamic simulations was used to find new microalgae metabolites inhibitors of β-lactamase. Six 3D β-lactamase proteins were selected, and the molecular docking revealed that the metabolites belonging to the same structural families, such as phenylacridine (4-Ph), quercetin (Qn), and cryptophycin (Cryp), exhibit a better binding score and binding energy than commercial clinical medicine β-lactamase inhibitors, such as clavulanic acid, sulbactam, and tazobactam. These results indicate that 4-Ph, Qn, and Cryp molecules, homologous from microalgae metabolites, could be used, likely as novel β-lactamase inhibitors or as structural templates for new in-silico pharmaceutical designs, with the possibility of combatting β-lactam resistance     

Structural and Dynamical Insight into PPARγ Antagonism: In Silico Study of the Ligand-Receptor Interactions of Non-Covalent Antagonists

The structural and dynamical properties of the peroxisome proliferator-activated receptor γ (PPARγ) nuclear receptor have been broadly studied in its agonist state but little is known about the key features required for the receptor antagonistic activity. Here we report a series of molecular dynamics (MD) simulations in combination with free energy estimation of the recently discovered class of non-covalent PPARγ antagonists. Their binding modes and dynamical behavior are described in details. Two key interactions have been detected within the cavity between helices H3, H11 and the activation helix H12, as well as with H12. The strength of the ligand-amino acid residues interactions has been analyzed in relation to the specificity of the ligand dynamical and antagonistic features. According to our results, the PPARγ activation helix does not undergo dramatic conformational changes, as seen in other nuclear receptors, but rather perturbations that occur through a significant ligand-induced reshaping of the ligand-receptor and the receptor-coactivator binding pockets. The H12 residue Tyr473 and the charge clamp residue Glu471 play a central role for the receptor transformations. Our results also demonstrate that MD can be a helpful tool for the compound phenotype characterization (full agonists, partial agonists or antagonists) when insufficient experimental data are available.     

In Silico Investigation of Potential Applications of Gamma Carbonic Anhydrases as Catalysts of CO2 Biomineralization Processes: A Visit to the Thermophilic Bacteria Persephonella hydrogeniphila,Persephonella marina,Thermosulfidibacter takaii,and Thermus thermophilus

Carbonic anhydrases (CAs) have been identified as ideal catalysts for CO₂ sequestration. Here, we report the sequence and structural analyses as well as the molecular dynamics (MD) simulations of four γ-CAs from thermophilic bacteria. Three of these, Persephonella marina, Persephonella hydrogeniphila, and Thermosulfidibacter takaii originate from hydrothermal vents and one, Thermus thermophilus HB8, from hot springs. Protein sequences were retrieved and aligned with previously characterized γ-CAs, revealing differences in the catalytic pocket residues. Further analysis of the structures following homology modeling revealed a hydrophobic patch in the catalytic pocket, presumed important for CO₂ binding. Monitoring of proton shuttling residue His69 (P. marina γ-CA numbering) during MD simulations of P. hydrogeniphila and P. marina’s γ-CAs (γ-PhCA and γ-PmCA), showed a different behavior to that observed in the γ-CA of Escherichia coli, which periodically coordinates Zn²⁺. This work also involved the search for hotspot residues that contribute to interface stability. Some of these residues were further identified as key in protein communication via betweenness centrality metric of dynamic residue network analysis. T. takaii’s γ-CA showed marginally lower thermostability compared to the other three γ-CA proteins with an increase in conformations visited at high temperatures being observed. Hydrogen bond analysis revealed important interactions, some unique and others common in all γ-CAs, which contribute to interface formation and thermostability. The seemingly thermostable γ-CA from T. thermophilus strangely showed increased unsynchronized residue motions at 423 K. γ-PhCA and γ-PmCA were, however, preliminarily considered suitable as prospective thermostable CO₂ sequestration agents.     

Mechanistic Insights into the Polymorphic Associations and Cross-Seeding of Aβ and hIAPP in the Presence of Histidine Tautomerism: An All-Atom Molecular Dynamic Study

Hundreds of millions of people around the world have been affected by Type 2 diabetes (T2D) which is a metabolic disorder. Clinical research has revealed T2D as a possible risk factor for Alzheimer’s disease (AD) development (and vice versa). Amyloid-β (Aβ) and human islet amyloid polypeptide are the main pathological species in AD and T2D, respectively. However, the mechanisms by which these two amyloidogenic peptides co-aggregate are largely uninvestigated. Herein, for the first time, we present the cross-seeding between Amylin1-37 and Aβ40 considering the particular effect of the histidine tautomerism at atomic resolution applying the all-atom molecular dynamics (MD) simulations for heterodimeric complexes. The results via random seed MD simulations indicated that the Aβ40(δδδ) isomer in cross-talking with Islet(ε) and Islet(δ) isomers could retain or increase the β-sheet content in its structure that may make it more prone to further aggregation and exhibit higher toxicity. The other tautomeric isomers which initially did not have a β-sheet structure in their monomeric forms did not show any generated β-sheet, except for one seed of the Islet(ε) and Aβ40(εεε) heterodimers complex that displayed a small amount of formed β-sheet. This computational research may provide a different point of view to examine all possible parameters that may contribute to the development of AD and T2D and provide a better understanding of the pathological link between these two severe diseases.     

Revisiting the Proposition of Binding Pockets and Bioactive Poses for GSK-3β Allosteric Modulators Addressed to Neurodegenerative Diseases

Glycogen synthase kinase-3 beta (GSK-3β) is an enzyme pertinently linked to neurodegenerative diseases since it is associated with the regulation of key neuropathological features in the central nervous system. Among the different kinds of inhibitors of this kinase, the allosteric ones stand out due to their selective and subtle modulation, lowering the chance of producing side effects. The mechanism of GSK-3β allosteric modulators may be considered still vague in terms of elucidating a well-defined binding pocket and a bioactive pose for them. In this context, we propose to reinvestigate and reinforce such knowledge by the application of an extensive set of in silico methodologies, such as cavity detection, ligand 3D shape analysis and docking (with robust validation of corresponding protocols), and molecular dynamics. The results here obtained were consensually consistent in furnishing new structural data, in particular by providing a solid bioactive pose of one of the most representative GSK-3β allosteric modulators. We further applied this to the prospect for new compounds by ligand-based virtual screening and analyzed the potential of the two obtained virtual hits by quantum chemical calculations. All potential hits achieved will be subsequently tested by in vitro assays in order to validate our approaches as well as to unveil novel chemical entities as GSK-3β allosteric modulators.     

Structure and Dynamics of Meprin β in Complex with a Hydroxamate-Based Inhibitor

The astacin protease Meprin β represents an emerging target for drug development due to its potential involvement in disorders such as acute and chronic kidney injury and fibrosis. Here, we elaborate on the structural basis of inhibition by a specific Meprin β inhibitor. Our analysis of the crystal structure suggests different binding modes of the inhibitor to the active site. This flexibility is caused, at least in part, by movement of the C-terminal region of the protease domain (CTD). The CTD movement narrows the active site cleft upon inhibitor binding. Compared with other astacin proteases, among these the highly homologous isoenzyme Meprin α, differences in the subsites account for the unique selectivity of the inhibitor. Although the inhibitor shows substantial flexibility in orientation within the active site, the structural data as well as binding analyses, including molecular dynamics simulations, support a contribution of electrostatic interactions, presumably by arginine residues, to binding and specificity. Collectively, the results presented here and previously support an induced fit and substantial movement of the CTD upon ligand binding and, possibly, during catalysis. To the best of our knowledge, we here present the first structure of a Meprin β holoenzyme containing a zinc ion and a specific inhibitor bound to the active site. The structural data will guide rational drug design and the discovery of highly potent Meprin inhibitors.     

Molecular Dynamics Simulations of Acylpeptide Hydrolase Bound to Chlorpyrifosmethyl Oxon and Dichlorvos

Acylpeptide hydrolases (APHs) catalyze the removal of N-acylated amino acids from blocked peptides. Like other prolyloligopeptidase (POP) family members, APHs are believed to be important targets for drug design. To date, the binding pose of organophosphorus (OP) compounds of APH, as well as the different OP compounds binding and inducing conformational changes in two domains, namely, α/β hydrolase and β-propeller, remain poorly understood. We report a computational study of APH bound to chlorpyrifosmethyl oxon and dichlorvos. In our docking study, Val471 and Gly368 are important residues for chlorpyrifosmethyl oxon and dichlorvos binding. Molecular dynamics simulations were also performed to explore the conformational changes between the chlorpyrifosmethyl oxon and dichlorvos bound to APH, which indicated that the structural feature of chlorpyrifosmethyl oxon binding in APH permitted partial opening of the β-propeller fold and allowed the chlorpyrifosmethyl oxon to easily enter the catalytic site. These results may facilitate the design of APH-targeting drugs with improved efficacy.     

DAG/PKCδ and IP3/Ca2+/CaMK IIβ Operate in Parallel to Each Other in PLCγ1-Driven Cell Proliferation and Migration of Human Gastric Adenocarcinoma Cells,through Akt/mTOR/S6 Pathway

Phosphoinositide specific phospholipase Cγ (PLCγ) activates diacylglycerol (DAG)/protein kinase C (PKC) and inositol 1,4,5-trisphosphate (IP3)/Ca²⁺/calmodulin-dependent protein kinase II (CaMK II) axes to regulate import events in some cancer cells, including gastric adenocarcinoma cells. However, whether DAG/PKCδ and IP3/Ca²⁺/CaMK IIβ axes are simultaneously involved in PLCγ1-driven cell proliferation and migration of human gastric adenocarcinoma cells and the underlying mechanism are not elucidated. Here, we investigated the role of DAG/PKCδ or CaMK IIβ in PLCγ1-driven cell proliferation and migration of human gastric adenocarcinoma cells, using the BGC-823 cell line. The results indicated that the inhibition of PKCδ and CaMK IIβ could block cell proliferation and migration of BGC-823 cells as well as the effect of inhibiting PLCγ1, including the decrease of cell viability, the increase of apoptotic index, the down-regulation of matrix metalloproteinase (MMP) 9 expression level, and the decrease of cell migration rate. Both DAG/PKCδ and CaMK IIβ triggered protein kinase B (Akt)/mammalian target of rapamycin (mTOR)/S6 pathway to regulate protein synthesis. The data indicate that DAG/PKCδ and IP3/Ca²⁺/CaMK IIβ operate in parallel to each other in PLCγ1-driven cell proliferation and migration of human gastric adenocarcinoma cells through Akt/mTOR/S6 pathway, with important implication for validating PLCγ1 as a molecular biomarker in early gastric cancer diagnosis and disease surveillance.