Search

In our quest to discover effective inhibitors against severe acute respiratory syndrome coronavirus 2 helicase, a diverse set of more than 300 naturally occurring antiviral metabolites was investigated. Employing advanced computational techniques, we initiated the selection process by analyzing and comparing the co-crystallized ligand (VXG) of the severe acute respiratory syndrome coronavirus 2 helicase protein (PDB ID: 5RMM) to identify compounds with structurally similar features and potential for comparable binding. Through structural similarity and pharmacophore research, 13 compounds that shared important characteristics with VXG were pinpointed. Subsequently, these candidates were subjected to molecular docking to identify seven compounds that demonstrated favorable energy profiles and accurate binding to the severe acute respiratory syndrome coronavirus 2 helicase. Among these, mycophenolic acid emerged as the most promising candidate. To ensure the safety and viability of the selected compounds, we conducted ADMET tests, which confirmed the favorable characteristics of mycophenolic acid, and the safety of atropine and plumbagin. Building on these results, we performed additional analyses on mycophenolic acid, including various molecular dynamics simulations. These investigations demonstrated that mycophenolic acid exhibited optimal binding to the severe acute respiratory syndrome coronavirus 2 helicase, maintaining flawless dynamics throughout the simulations. Furthermore, the Molecular Mechanics Poisson–Boltzmann Surface Area tests provided strong evidence that mycophenolic acid successfully formed a stable connection with the severe acute respiratory syndrome coronavirus 2 helicase, with a calculated free energy value of −294 kJ mol−1. These encouraging findings provide a solid foundation for further research, including in vitro and in vivo studies, on the three identified compounds. The potential efficacy of these compounds as treatment options for coronavirus-19 warrants further exploration and may hold significant promise in the ongoing fight against the pandemic. [ABSTRACT FROM AUTHOR]

The detrimental effect of coronavirus disease 2019 (COVID-19) pandemic has manifested itself as a global crisis. Currently, no specific treatment options are available for COVID-19, so therapeutic interventions to tackle the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection must be urgently established. Therefore, cohesive and multidimensional efforts are required to identify new therapies or investigate the efficacy of small molecules and existing drugs against SARS-CoV-2. Since the RNA-dependent RNA Polymerase (RdRP) of SARS-CoV-2 is a promising therapeutic target, this study addresses the identification of antiviral molecules that can specifically target SARS-CoV-2 RdRP. The computational approach of drug development was used to screen the antiviral molecules from two antiviral libraries (Life Chemicals [LC] and ASINEX) against RdRP. Here, we report six antiviral molecules (F3407-4105, F6523-2250, F6559-0746 from LC and BDG 33693278, BDG 33693315, LAS 34156196 from ASINEX), which show substantial interactions with key amino acid residues of the active site of SARS-CoV-2 RdRP and exhibit higher binding affinity (>7.5 kcalmol−1) than Galidesivir, an Food and Drug Administration-approved inhibitor of the same. Further, molecular dynamics simulation and Molecular Mechanics Poisson-Boltzmann Surface Area results confirmed that identified molecules with RdRP formed higher stable RdRP-inhibitor(s) complex than RdRP-Galidesvir complex. Our findings suggest that these molecules could be potential inhibitors of SARS-CoV-2 RdRP. However, further in vitro and preclinical experiments would be required to validate these potential inhibitors of SARS-CoV-2 protein. [ABSTRACT FROM AUTHOR]

Both the inactive‐ and active‐state CB1 receptor crystal structures have now been solved, allowing their application in various structure‐based drug design methods. One potential method utilizing these crystal structures is the Molecular Mechanics/Poisson–Boltzmann Surface Area (MM/PBSA) method of predicting relative binding free [ABSTRACT FROM AUTHOR]

4-Hydroxyphenylpyruvate dioxygenase (HPPD) is a significant enzyme in the biosynthesis of plastoquinone and tocopherol. Moreover, it is also a potential target to develop new herbicide. The technology of computer-aided drug design (CADD) is a useful tool in the efficient discovery of new HPPD inhibitors. Forty-three compounds with known activities were used to generate comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models based on common framework and molecular docking. The structural contribution to the activity was determined, which provided further information for the design of novel inhibitors. Molecular docking was used to explain the changes in activity caused by the binding mode between ligand and protein. The molecular dynamics (MD) results indicated that the electrostatic energy was the major driving force for ligand–protein interaction and the Phe403 made the greatest contribution to the binding. The present work has provided useful information for the rational design of novel HPPD inhibitors with improved activity.

In the present work a combination of traditional and steered molecular dynamics based techniques were employed to identify potential inhibitors against the human BRD4 protein (BRD4- BD1); an established drug target for multiple illnesses including various malignancies. Quinoline derivatives that were synthesized in-house were tested for their potential as new BRD4-BD1 inhibitors. Initially molecular docking experiments were performed to determine the binding poses of BRD4-BD1 inhibitors. To learn more about the thermodynamics of inhibitor binding to the BRD4-BD1 active site, the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) free energy calculations were conducted afterwards. The findings of the MM-PBSA analysis were further reinforced by performing steered umbrella sampling simulations which revealed crucial details about the binding/unbinding process of the most potent quinoline derivatives at the BRD4-BD1 active site. We report a novel quinoline derivative which can be developed into a fully functional BRD4-BD1 inhibitor after experimental validation. The identified compound (4 g) shows better properties than the standard BRD4-BD1 inhibitors considered in the study. The study also highlights the crucial role of Gln78, Phe79, Trp81, Pro82, Phe83, Gln84, Gln85, Val87, Leu92, Leu94, Tyr97, Met105, Cys136, Asn140, Ile146 and Met149 in inhibitor binding. The study provides a possible lead candidate and key amino acids involved in inhibitor recognition and binding at the active site of BRD4-BD1 protein. The findings might be of significance to medicinal chemists involved in the development of potent BRD4-BD1 inhibitors., (© 2024 Wiley Periodicals LLC.)

The vasoactive proline-rich oligopeptide termed BPP-BrachyNH 2 (H-WPPPKVSP-NH 2 ) induces in vitro inhibitory activity of angiotensin I-converting enzyme (ACE) in rat blood serum. In the present study, the removal of N-terminal tryptophan or C-terminal proline from BPP-BrachyNH 2 was investigated in order to predict which structural components are important or required for interaction with ACE. Furthermore, the toxicological profile was assessed by in silico prediction and in vitro MTT assay. Two BPP-BrachyNH 2 analogues (des-Trp 1 -BPP-BrachyNH 2 and des-Pro 8 -BPP-BrachyNH 2 ) were synthesized, and in vitro and in silico ACE inhibitory activity and toxicological profile were assessed. The des-Trp 1 -BPP-BrachyNH 2 and des-Pro 8 -BPP-BrachyNH 2 were respectively 3.2- and 29.5-fold less active than the BPP-BrachyNH 2 -induced ACE inhibitory activity. Molecular Dynamic and Molecular Mechanics Poisson-Boltzmann Surface Area simulations (MM-PBSA) demonstrated that the ACE/BBP-BrachyNH 2 complex showed lower binding and van der Wall energies than the ACE/des-Pro 8 -BPP-BrachyNH 2 complex, therefore having better stability. The removal of the N-terminal tryptophan increased the in silico predicted toxicological effects and cytotoxicity when compared with BPP-BrachyNH 2 or des-Pro 8 -BPP-BrachyNH 2 . Otherwise, des-Pro 8 -BPP-BrachyNH 2 was 190-fold less cytotoxic than BPP-BrachyNH 2 . Thus, the removal of C-terminal proline residue was able to markedly decrease both the BPP-BrachyNH 2 -induced ACE inhibitory and cytotoxic effects assessed by in vitro and in silico approaches. In conclusion, the aminoacid sequence of BPP-BrachyNH 2 is essential for its ACE inhibitory activity and associated with an acceptable toxicological profile. The perspective of the interactions of BPP-BrachyNH 2 with ACE found in the present study can be used for development of drugs with differential therapeutic profile than current ACE inhibitors. [ABSTRACT FROM AUTHOR]

In this study, we estimated the optimum concentration of copper ions that are effective in the stability and the structural changes of human growth hormone (hGH) protein in the combination of different concentrations of these ions at the molecular level using molecular dynamics simulation by Gromacs 4.6.5 software. Moreover, to estimate the binding affinity of copper ions to hGH protein, binding free energies is calculated by the molecular mechanics Poisson-Boltzmann Surface Area (MM-PBSA). The analysis of molecular dynamics (MD) trajectories as dictionary of the secondary structure of protein (DSSP), solvent accessible surface area (SASA) and binding free energy calculations show that hGH protein structure is more stabilized by increasing a limited concentration of copper ions. These findings align with our previous experimental studies. [ABSTRACT FROM AUTHOR]

ABSTRACT Variation among crystal structures of the λ Cro dimer highlights conformational flexibility. The structures range from a wild type closed to a mutant fully open conformation, but it is unclear if each represents a stable solution state or if one may be the result of crystal packing. Here we use molecular dynamics (MD) simulation to investigate the energetics of crystal packing interfaces and the influence of site-directed mutagenesis on them in order to examine the effect of crystal packing on wild type and mutant Cro dimer conformation. Replica exchange MD of mutant Cro in solution shows that the observed conformational differences between the wild type and mutant protein are not the direct consequence of mutation. Instead, simulation of Cro in different crystal environments reveals that mutation affects the stability of crystal forms. Molecular Mechanics Poisson-Boltzmann Surface Area binding energy calculations reveal the detailed energetics of packing interfaces. Packing interfaces can have diverse properties in strength, energetic components, and some are stronger than the biological dimer interface. Further analysis shows that mutation can strengthen packing interfaces by as much as ∼5 kcal/mol in either crystal environment. Thus, in the case of Cro, mutation provides an additional energetic contribution during crystal formation that may stabilize a fully open higher energy state. Moreover, the effect of mutation in the lattice can extend to packing interfaces not involving mutation sites. Our results provide insight into possible models for the effect of crystallization on Cro conformational dynamics and emphasize careful consideration of protein crystal structures. Proteins 2014; 82:1128-1141. © 2013 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

ABSTRACT The ectodomain of the human epidermal growth factor receptor (hEGFR) controls input to several cell signalling networks via binding with extracellular growth factors. To gain insight into the dynamics and ligand binding of the ectodomain, the hEGFR monomer was subjected to molecular dynamics simulation. The monomer was found to be substantially more flexible than the ectodomain dimer studied previously. Simulations where the endogeneous ligand EGF binds to either Subdomain I or Subdomain III, or where hEGFR is unbound, show significant differences in dynamics. The molecular mechanics Poisson-Boltzmann surface area method has been used to derive relative free energies of ligand binding, and we find that the ligand is capable of binding either subdomain with a slight preference for III. Alanine-scanning calculations for the effect of selected ligand mutants on binding reproduce the trends of affinity measurements. Taken together, these results emphasize the possible role of the ectodomain monomer in the initial step of ligand binding, and add details to the static picture obtained from crystal structures. Proteins 2013; 81:1931-1943. © 2013 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

In the later stages of drug design projects, accurately predicting relative binding affinities of chemically similar compounds to a biomolecular target is of utmost importance for making decisions based on the ranking of such compounds. So far, the extensive application of binding free energy approaches has been hampered by the complex and time-consuming setup of such calculations. We introduce the free energy workflow (FEW) tool that facilitates setup and execution of binding free energy calculations with the AMBER suite for multiple ligands. FEW allows performing free energy calculations according to the implicit solvent molecular mechanics (MM-PB(GB)SA), the linear interaction energy, and the thermodynamic integration approaches. We describe the tool's architecture and functionality and demonstrate in a show case study on Factor Xa inhibitors that the time needed for the preparation and analysis of free energy calculations is considerably reduced with FEW compared to a fully manual procedure. © 2013 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

Piperlongumine (PPL) is an alkaloid extracted from several pepper species that exhibits anti-inflammatory and anti-carcinogenic properties. Nevertheless, the molecular mode of action of PPL that confers such powerful pharmacological properties remains unknown. From this perspective, spectroscopic methods aided by computational modeling were employed to characterize the interaction between PPL and nucleotide-binding domain of heat shock protein 70 (NBD/HSP70), which is involved in the pathogenesis of several diseases. Steady-state fluorescence spectroscopy along with time-resolved fluorescence revealed the complex formation based on a static quenching mechanism. Van't Hoff analyses showed that the binding of PPL toward NBD is driven by equivalent contributions of entropic and enthalpic factors. Furthermore, IDF and Scatchard methods applied to fluorescence intensities determined two cooperative binding sites with Kb of (6.3 ± 0.2) × 104 M−1. Circular dichroism determined the thermal stability of the NBD domain and showed that PPL caused minor changes in the protein secondary structure. Computational simulations elucidated the microenvironment of these interactions, showing that the binding sites are composed mainly of polar amino acids and the predominant interaction of PPL with NBD is Van der Waals in nature. [ABSTRACT FROM AUTHOR]

4-Hydroxyphenylpyruvate dioxygenase (HPPD) is a significant enzyme in the biosynthesis of plastoquinone and tocopherol. Moreover, it is also a potential target to develop new herbicide. The technology of computer-aided drug design (CADD) is a useful tool in the efficient discovery of new HPPD inhibitors. Forty-three compounds with known activities were used to generate comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models based on common framework and molecular docking. The structural contribution to the activity was determined, which provided further information for the design of novel inhibitors. Molecular docking was used to explain the changes in activity caused by the binding mode between ligand and protein. The molecular dynamics (MD) results indicated that the electrostatic energy was the major driving force for ligand-protein interaction and the Phe403 made the greatest contribution to the binding. The present work has provided useful information for the rational design of novel HPPD inhibitors with improved activity.

Tyrosine hydroxylase (TH) catalyzes the rate-limiting step in the synthesis of catecholamine neurotransmitters, and a reduction in TH activity is associated with several neurological diseases. Human TH is regulated, among other mechanisms, by Ser19-phosphorylation-dependent interaction with 14-3-3 proteins. The N-terminal sequence (residues 1-43), which corresponds to an extension to the TH regulatory domain, also interacts with negatively charged membranes. By using X-ray crystallography together with molecular dynamics simulations and structural bioinformatics analysis, we have probed the conformations of the Ser19-phosphorylated N-terminal peptide [THp-(1-43)] bound to 14-3-3γ, free in solution and bound to a phospholipid bilayer, and of the unphosphorylated peptide TH-(1-43) both free and bilayer bound. As seen in the crystal structure of THp-(1-43) complexed with 14-3-3γ, the region surrounding pSer19 adopts an extended conformation in the bound state, whereas THp-(1-43) adopts a bent conformation when free in solution, with higher content of secondary structure and higher number of internal hydrogen bonds. TH-(1-43) in solution presents the highest mobility and least defined structure of all forms studied, and it shows an energetically more favorable interaction with membranes relative to THp-(1-43). Cationic residues, notably Arg15 and Arg16, which are the recognition sites of the kinases phosphorylating at Ser19, are also contributing to the interaction with the membrane. Our results reveal the structural flexibility of this region of TH, in accordance with the functional versatility and conformational adaptation to different partners. Furthermore, this structural information has potential relevance for the development of therapeutics for neurodegenerative disorders, through modulation of TH-partner interactions., (© 2013.)