Using thermodynamic integration MD simulation to compute relative protein–ligand binding free energy of a GSK3β kinase inhibitor and its analogs

Hexaconazole exposure may lead to Parkinson via disrupting glucocerebrosidase and parkin: molecular interaction, dynamics, MMPBSA and DFT based in-silico predictive toxicology

Hexaconazole is a known fungicide for agricultural purposes. It has bioaccumulation ability which makes it important for its toxicological characterization. There are various neurological impacts of pollutants on human health. Therefore, in this study, we have done predictive analyses of the interaction mechanism of hexaconazole by molecular interaction analysis, molecular dynamics simulation, and Poisson-Boltzmann surface area (MM-PBSA) to assess hexaconazole's potency to disrupt the homeostasis of glucocerebrosidase (-7.9 kcal/mol) and parkin (-5.67 kcal/mol) proteins which have significant roles in the manifestation of Parkinson disease. The findings reveal that hexaconazole has the potency to form stable interactions with glucocerebrosidase and parkin. This research provides a molecular and atomic-level understanding of how hexaconazole exposure may disrupt the homeostasis of glucocerebrosidase and parkin. The root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration, and hydrogen bonding exhibited the potent molecular interactions of hexaconazole, which may lead to neurological manifestations such as Parkinson disease.

Potency of Hexaconazole to Disrupt Endocrine Function with Sex Hormone-Binding Globulin

Hexaconazole is widely used as a fungicide for agricultural purposes. However, the endocrine-disrupting potential of hexaconazole is still under investigation. In addition, an experimental study found that hexaconazole may disrupt the normal synthesis of steroidal hormones. The potency of hexaconazole to bind with sex hormone-binding globulin (SHBG), a plasma carrier protein that binds androgens and oestrogens, is unknown. In this study, we evaluated the efficacy of hexaconazole to bind with SHBG by molecular interaction, a molecular dynamics method. In addition, principal component analysis was performed to understand the dynamical behaviour of hexaconazole with SHBG in comparison with dihydrotestosterone and aminoglutethimide. The binding scores of hexaconazole, dihydrotestosterone, and aminoglutethimide with SHBG were found to be -7.12 kcal/mol, -11.41 kcal/mol, and -6.84 kcal/mol, respectively. With respect to stable molecular interaction, hexaconazole showed similar molecular dynamics patterns of root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), and hydrogen bonding. The solvent surface area (SASA) and principal component analysis (PCA) of hexaconazole exhibit similar patterns in comparison with dihydrotestosterone and aminoglutethimide. These results show that hexaconazole has a stable molecular interaction with SHBG, which may acquire the active site of the native ligand, resulting in significant endocrine disruption during agricultural work.

Sunscreen Ingredient Octocrylene's Potency to Disrupt Vitamin D Synthesis

Octocrylene is a widely used ingredient in sunscreen products, and it has been observed that the use of sunscreen has been increasing over the last few decades. In this paper, we investigated the way in which sunscreen's ingredient octocrylene may disrupt normal vitamin D synthesis pathway, resulting in an imbalance in vitamin D levels in the body. The key techniques used for this insilico investigation were molecular docking, molecular dynamic (MD) simulation, and MMPBSA-based assessment. Vitamin D abnormalities have become very common in human health. Unknown exposure to chemicals may be one of the important risk factors. In molecular docking analysis, octocrylene exhibited a binding energy of -11.52 kcal/mol with vitamin D binding protein (1KXP) and -11.71 for the calcitriol native ligand. Octocrylene had a binding potency of -11.152 kcal/mol with the vitamin D receptor (1DB1), and calcitriol had a binding potency of -8.73 kcal/mol. In addition, octocrylene has shown binding energy of -8.96 kcal/mol with CYP2R1, and the calcitriol binding energy was -10.36 kcal/mol. Regarding stability, the root-mean-square deviation (RMSD), the root-mean-square fluctuation (RMSF), the radius of gyration, hydrogen bonding, and the solvent-accessible surface area (SASA) exhibited that octocrylene has a stable binding pattern similar to calcitriol. These findings revealed that incessant exposure to octocrylene may disrupt normal vitamin D synthesis.

Dynamics of metal binding and mutation in yybP-ykoY riboswitch of Lactococcus lactis

Riboswitch is a regulatory segment of messenger RNA (mRNA), which by binding to various cellular metabolites regulates the activity of mRNA via modulating transcription, translation, alternative splicing, and stability of the mRNA. yybP–ykoY riboswitch of Lactococcus lactis, which is present upstream of the yoaB gene, functions as a Mn²⁺-specific genetic ON-switch, and modulates expression of proteins which are significant for Mn²⁺ homeostasis. The P1.1 switch helix of the aptamer domain of the riboswitch contains an intrinsic transcription terminator structure, which gets stabilized with Mn²⁺ binding and causes disruption of terminator structure and allows the continuation of transcription. The current research work involved the evaluation of structural and dynamical properties of the yybP-ykoY riboswitch of L. lactis in its Mn²⁺-free, Mn²⁺-bound (wild-type), and Mn²⁺-bound mutant (A41U) states by applying molecular dynamics simulations. Based on the simulations, the effects of Mn²⁺ absence and A41U mutation were evaluated on the structure and dynamics of the riboswitches followed by the computation of the free energy of metal binding in the wild-type and the mutant riboswitches. The simulation results provided insights into the properties of the riboswitch with the focus on the dynamics of the P1.1 switch helix, and the manganese binding site designated as M_B site, as well as the relative stability of the wild-type and the mutant riboswitches, which helped to understand the structural and dynamical role of the metal ion involved in the function of Mn²⁺-sensing riboswitch.

The current research work involved the evaluation of structural and dynamical properties of yybP–ykoY riboswitch of L. lactis in Mn²⁺-free, Mn²⁺-bound (wild-type), and Mn²⁺-bound mutant (A41U) states by applying molecular dynamics simulations.

Screening cyclooxygenase-2 inhibitors from Allium sativum L. compounds: in silico approach

Inflammation is a natural protective response toward various simulators, including tissue damage or pathogens. The cyclooxygenase-2 (COX-2) is a very important protein in triggering pain and inflammation. Previous studies have claimed that Allium sativum offers a wide range of anti-inflammatory therapeutics for human consumption. Drug discovery is a complicated process, though in silico methods can make this procedure simpler and more cost-effective. At the current study, we performed the virtual screening of eight Allium sativum-derived compounds via molecular docking with COX-2 enzyme and confirmed the binding energy by docking score estimate followed by ADMET and drug-likeness investigation. The resulting highest-docking scored compound was exposed to molecular dynamics simulation (MDS) for evaluating stability of the docked enzyme-ligand complex and to gauge the oscillation and conformational alterations for the time of enzyme-ligand interaction. The factors of RMSD, RMSF, hydrogen bond interactions, and Rg after 100 ns of MDS proved the stability of alliin in the active site of COX-2 in comparison with celecoxib (CEL) as the control. Moreover, we investigated the binding affinity analysis of all compounds via MM/PBSA method. The results from this study suggest that alliin (a sulfuric compound) exhibits a higher binding affinity for the COX-2 enzyme compared to the other compounds and CEL. Alliin showed to be a possible anti-inflammatory therapeutic candidate for managing the inflammatory conditions.

Calculation of Energy for RNA/RNA and DNA/RNA Duplex Formation by Molecular Dynamics Simulation

Abstract

The development of approaches for predictive calculation of hybridization properties of various nucleic acid (NA) derivatives is the basis for the rational design of the NA-based constructs. Modern advances in computer modeling methods provide the feasibility of these calculations. We have analyzed the possibility of calculating the energy of DNA/RNA and RNA/RNA duplex formation using representative sets of complexes (65 and 75 complexes, respectively). We used the classical molecular dynamics (MD) method, the MMPBSA or MMGBSA approaches to calculate the enthalpy (ΔH°) component, and the quasi-harmonic approximation (Q-Harm) or the normal mode analysis (NMA) methods to calculate the entropy (ΔS°) contribution to the Gibbs energy ($$\Delta G_{{37}}^{^\circ }$$ ) of the NA complex formation. We have found that the MMGBSA method in the analysis of the MD trajectory of only the NA duplex and the empirical linear approximation allow calculation of the enthalpy of formation of the DNA, RNA, and hybrid duplexes of various lengths and GC content with an accuracy of 8.6%. Within each type of complex, the combination of rather efficient MMGBSA and Q-Harm approaches being applied to the trajectory of only the bimolecular complex makes it possible to calculate the $$\Delta G_{{37}}^{^\circ }$$ of the duplex formation with an error value of 10%. The high accuracy of predictive calculation for different types of natural complexes (DNA/RNA, DNA/RNA, and RNA/RNA) indicates the possibility of extending the considered approach to analogs and derivatives of nucleic acids, which gives a fundamental opportunity in the future to perform rational design of new types of NA-targeted sequence-specific compounds.

Morusflavone, a New Therapeutic Candidate for Prostate Cancer by CYP17A1 Inhibition: Exhibited by Molecular Docking and Dynamics Simulation

Morusflavone, a flavonoid from Morus alba L., was evaluated for its interactive ability and stability with CYP17A1, in comparison with abiraterone, which is a Food and Drug Administration (FDA)-approved CYP17A1 inhibitor. CYP17A1 inhibition is an important therapeutic target for prostate cancer. The CHAMM36 force field was used to perform molecular dynamics (MD) simulations in this study. The results show that Morusflavone has significant interactive ability and stability for CYP17A1, in comparison with abiraterone. The final interaction energies for the Morusflavone–CYP17A1 and abiraterone–CYP17A1 complexes were −246.252 KJ/mol and −207.86 KJ/mol, respectively. Since there are only limited therapeutic agents available, such as abiraterone, galeterone, and seviteronel, which are being developed for prostate cancer, information on any potent natural anticancer compounds, such as vinca alkaloids, for prostate cancer treatment is limited. The results of this study show that CYP17A1 inhibition by Morusflavone could be an important therapeutic target for prostate cancer. Further preclinical and clinical evaluations of the lead compound Morusflavone are required to evaluate whether it can serve as a potential inhibitor of CYP17A1, which will be a new hope for prostate cancer treatment.

Molecular docking and dynamics simulation to screen interactive potency and stability of fungicide tebuconazole with thyroid and sex hormone-binding globulin: Implications of endocrine and reproductive interruptions

Tebuconazole is a widely used fungicide in agriculture, and it may easily enter in the human food chain. In addition, tebuconzaol skin permeation coefficient (Log Kp) is -5.55 cm/s and it does not violate Lipinski's rule. It may mimic as a ligand for various endocrine and reproductive receptor leading to toxicological response or disease manifestation. We studied interactive potential of tebuconazole with thyroid and sex hormone-binding globulin. The main methods for this in silico analyses are molecular docking and molecular dynamic (MD) simulation. This paper explores how agriculture fungicide tebuconzaol exposure can be a risk for endocrine and reprotoxicity due to its stable interactive potency with thyroid and sex hormone-binding globulin (2CEO and 1D2S). Thyroid impairment is one of the most common endocrine issues in human health. In molecular docking analyses, tebuconazole exhibited binding potency of -6.28 kcal/mol with 2CEO compared to its native ligand thyroxin and inhibitor propylthiouracil which had the binding potency of -9.9 and -4.49 kcal/mol, respectively. The binding score of tebuconzaol with 1D2S was found to be -7.54 kcal/mol compared to native ligand dihydrotestosteron and inhibitor aminoglutethimide which exhibited the binding score of -6.84 and -11.41 kcal/mol, respectively. Therefore, each complex was subjected to MD simulation for comparative assessment of physical movement. The root mean square deviation (RMSD), root mean square fluctuation (RMSF), Radius of Gyration and hydrogen bonding exhibited that fluconazole had stable binding pattern with 2CEO and 1D2S which was almost similar to native ligand and its inhibitor. Study revealed that tebuconazole may lead to potent endocrine and reproductive disruptions.

Synthesis of exfoliated multilayer graphene and its putative interactions with SARS-CoV-2 virus investigated through computational studies

Our work investigates the interaction of synthesized graphene with the SARS-CoV-2 virus using molecular docking and molecular dynamics (MD) simulation method. The layer dependent inhibitory effect of graphene nanosheets on spike receptor-binding domain of 6LZG, complexed with host receptor i.e. angiotensin-converting enzyme 2 (ACE2) of SARS-CoV-2 was investigated through computational study. Graphene sample was synthesized using mechanical exfoliation with shear stress and its mechanism of inhibition towards the SARS-CoV-2 virus was explored by molecular docking and molecular dynamics (MD) simulation method. The thermodynamics study for the free binding energy of graphene towards the SARS-CoV-2 virus was analyzed. The binding energy of graphene towards the virus increased with an increasing number of layers. It shows the highest affinity of -17.5 Kcal/mol in molecular docking while ΔG_binding is in the order of -28.01 ± 0.04 5 Kcal/mol for the seven-layers structure. The increase in carbon layers is associated with an increasing number of edge sp³ -type carbon, providing greater curvature, further increase the surface reactivity responsible for high binding efficiency. The MD simulation data reveals the high inhibition efficiency of the synthesized graphene towards SARS-CoV-2 virus which would help to design future in-vitro studies. The graphene system could find potential applications in personal protective equipment and diagnostic kits.Communicated by Ramaswamy H. Sarma.

Exploring pH Dependent Host/Guest Binding Affinities

When the electrostatic environment surrounding binding partners changes between unbound and bound states, the net uptake or release of a proton is possible by either binding partner. This process is pH-dependent in that the free energy required to uptake or release the proton varies with pH. This pH-dependence is typically not considered in conventional free energy methods where the use of fixed protonation states is the norm. In the present paper, we apply a simple two-step approach to calculate the pH-dependent binding free energy of a model cucubit[7]uril host/guest system. By use of λ-dynamics with an enhanced sampling protocol, adaptive landscape flattening, pK_a shifts and reference binding free energies upon complexation were determined. This information enables the construction of pH-dependent binding profiles that accurately capture the pK_a shifts and reproduce binding free energies at the different pH conditions that were observed experimentally. Our calculations illustrate a general framework for computing pH-dependent binding free energies but also point to some issues in modeling the molecular charge distributions within this series of molecules with CGenFF. However, by introducing some minor charge modifications to the CGenFF force field, we saw significant improvement in accuracy of the calculated pK_a shifts.

Evaluation of a sesquiterpene as possible drug lead against gelatinases via molecular dynamics simulations

Malignant tumors can be targeted by accounting for their metastatic capabilities. Matrix metalloproteinases (MMPs) are the key players in tumor metastasis facilitating through their proteolytic activities of angiogenesis and extracellular matrix components (ECM) degradation. MMP-2 and MMP-9 being the members of a distinguished class of MMPs more commonly known as gelatinases are the prominent enzymes which are involved in different cancer progression stages. Targeting these isoforms specifically has always been a challenging task due to highly similar structural and functional features among the other members of MMPs with well preserve active sites containing catalytic zinc atom that was the only reason that none of the MMP inhibitor has been successfully marketed for the tumor pathology up till now. Therefore, non-competitive inhibitors with different structural attributed are needed to be evaluated at the molecular level for further experiments. The present study deals with the application of molecular dynamics simulation for the investigation of an alternative pathway for the inhibition of MMP-2 and MMP-9 by a sesquiterpene isolated from Polygonum barbatum which demonstrates the characteristics binding to the S1' subsite of the enzymes followed by in vitro gene expression studies. The simulation results provide information on the possible binding profile producing inhibitory effects imposed by the inhibitor to these enzymes by acquiring different structural and dynamical features. Moreover, thermodynamic quantities based on the computationally intensive thermodynamic integration approach were also obtained in terms of inhibitor binding affinity computed for the inhibitor against MMP-2 and MMP-9 that completely augmented the experimental gene expression study.Communicated by Ramaswamy H. Sarma.

Molecular docking and dynamic approach to virtual screen inhibitors against Esbp of Candidatus Liberibacter asiaticus

Citrus greening (huanglongbing) is the most destructive disease of citrus worldwide caused by Candidatus Liberibacter asiaticus (CLA). Currently, no strategies have been developed to manage the Huanglongbing (HLB) disease and to stop the spreading of this disease to new citrus areas. Esbp is an extracellular solute-binding protein, involved in the uptake of iron in CLA. Thus, inhibiting this process may be a promising approach to design a drug against CLA. Thus, the present study focused on the identification of novel effective inhibitors which can inhibit the activity of CLas Esbp. A series of small molecules were screened against the CLas Esbp and the binding affinities were assessed using docking simulation studies. Top scored molecules were screened for different pharmacophore properties and Inhibitory Concentration 50 (IC50) values. Density functional theory was employed to check the chemical properties of the molecules. Further, Molecular Dynamics simulation analysis like RMSD, RMSF, Rg, SASA and MMPBSA results reveal that the identified molecules (ZINC03143779, ZINC05491830, ZINC19210425, ZINC08750867, and ZINC14671545) exhibit a good binding affinity for CLas Esbp and results in the formation of stable CLas Esbp-inhibitor(s) complex. The present study reported that these compounds appeared to be the suitable novel inhibitor of CLas Esbp and pave the way to further development of antimicrobial agents against CLA.

Rapid and accurate structure-based therapeutic peptide design using GPU accelerated thermodynamic integration

Peptide-based therapeutics are an alternative to small molecule drugs as they offer superior specificity, lower toxicity, and easy synthesis. Here we present an approach that leverages the dramatic performance increase afforded by the recent arrival of GPU accelerated thermodynamic integration (TI). GPU TI facilitates very fast, highly accurate binding affinity optimization of peptides against therapeutic targets. We benchmarked TI predictions using published peptide binding optimization studies. Prediction of mutations involving charged side-chains was found to be less accurate than for non-charged, and use of a more complex 3-step TI protocol was found to boost accuracy in these cases. Using the 3-step protocol for non-charged side-chains either had no effect or was detrimental. We use the benchmarked pipeline to optimize a peptide binding to our recently discovered cancer target: EME1. TI calculations predict beneficial mutations using both canonical and non-canonical amino acids. We validate these predictions using fluorescence polarization and confirm that binding affinity is increased. We further demonstrate that this increase translates to a significant reduction in pancreatic cancer cell viability.

Identification of novel phosphatidic acid binding domain on sphingosine kinase 1 of Arabidopsis thaliana

Phosphatidic acid (PA) is an important lipid signaling molecule which interacts with Arabidopsis thaliana Sphingosine kinase1 (AtSPHK1) during several abiotic stresses particularly drought stress as a result of Abscisic acid (ABA) signaling in guard cells. PA molecules respond by generating lipid signal and/or by binding and translocating target proteins to membrane. However, site of interaction and role of PA binding to AtSPHK1 is not clear yet. Owing to the importance of AtSPHK1 during stress signaling it is imperative to decipher the site of PA interaction with AtSPHK1. To identify the PA binding region of AtSPHK1, various deletion fragments from N-terminal and C-terminal region were prepared. Results from protein lipid overlay assay using various truncated proteins of AtSPHK1 suggested the involvement of N-terminal region, between 110 and 205 amino acids, in binding with PA. In-silico analyses performed to build homologous structure of AtSPHK1 revealed that PA docking occurs in the hydrophobic cavity of DAG-Kinase domain. Deletion of amino acids ¹⁸²VSGDGI¹⁸⁷ perturbed PA-AtSPHK1 binding, indicating an essential role of these six amino acids in PA-AtSPHK1 binding.

Fast and flexible gpu accelerated binding free energy calculations within the amber molecular dynamics package

Alchemical free energy (AFE) calculations based on molecular dynamics (MD) simulations are key tools in both improving our understanding of a wide variety of biological processes and accelerating the design and optimization of therapeutics for numerous diseases. Computing power and theory have, however, long been insufficient to enable AFE calculations to be routinely applied in early stage drug discovery. One of the major difficulties in performing AFE calculations is the length of time required for calculations to converge to an ensemble average. CPU implementations of MD-based free energy algorithms can effectively only reach tens of nanoseconds per day for systems on the order of 50,000 atoms, even running on massively parallel supercomputers. Therefore, converged free energy calculations on large numbers of potential lead compounds are often untenable, preventing researchers from gaining crucial insight into molecular recognition, potential druggability and other crucial areas of interest. Graphics Processing Units (GPUs) can help address this. We present here a seamless GPU implementation, within the PMEMD module of the AMBER molecular dynamics package, of thermodynamic integration (TI) capable of reaching speeds of >140 ns/day for a 44,907-atom system, with accuracy equivalent to the existing CPU implementation in AMBER. The implementation described here is currently part of the AMBER 18 beta code and will be an integral part of the upcoming version 18 release of AMBER. © 2018 Wiley Periodicals, Inc.

A computational approach for designing D-proteins with non-canonical amino acid optimised binding affinity

Redesigning protein surface topology to improve target binding holds great promise in the search for highly selective therapeutics. While significant binding improvements can be achieved using natural amino acids, the introduction of non-canonical residues vastly increases sequence space and thus the chance to significantly out-compete native partners. The potency of protein inhibitors can be further enhanced by synthesising mirror image, D-amino versions. This renders them non-immunogenic and makes them highly resistant to proteolytic degradation. Current experimental design methods often preclude the use of D-amino acids and non-canonical amino acids for a variety of reasons. To address this, we build an in silico pipeline for D-protein designs featuring non-canonical amino acids. For a test scaffold we use an existing D-protein inhibitor of VEGF: D-RFX001. We benchmark the approach by recapitulating previous experimental optimisation with canonical amino acids. Subsequent incorporation of non-canonical amino acids allows designs that are predicted to improve binding affinity by up to -7.18 kcal/mol.

A Molecular Docking and Dynamics Approach to Screen Potent Inhibitors Against Fosfomycin Resistant Enzyme in Clinical Klebsiella pneumoniae

Klebsiella pneumoniae, BA6753 was cultured from a patient in the Clinical Microbiology Laboratory of Christian Medical College. K. pneumoniae, BA6753 has a multidrug resistance plasmid encoding novel FosA variant-7, fosfomycin resistance enzyme. Minimal side effects and a wide range of bactericidal activity of fosfomycin have resulted in its expanded clinical use that prompts the rise of fosfomycin-resistant strains. At present, there are no effective inhibitors available to conflict the FosA-medicated fosfomycin resistance. To develop effective FosA inhibitors, it is crucial to understand the structural and dynamic properties of resistance enzymes. Hence, the present study focuses on the identification of potent inhibitors that can effectively bind to the fosfomycin resistance enzyme, thus predispose the target to inactivate by the second antibiotic. Initially, a series of active compounds were screened against the resistant enzyme, and the binding affinities were confirmed using docking simulation analysis. For efficient activity, the binding affinity of the resistance enzyme ought to be high with the inhibitor than the fosfomycin drug. Consequently, the enzyme-ligand complex which showed higher binding affinity than the fosfomycin was employed for subsequent analysis. The stability of the top scoring enzyme-ligand complex was further validated using molecular dynamics simulation studies. On the whole, we presume that the compound 19583672 demonstrates a higher binding affinity for the resistance enzyme comparing to other compounds and fosfomycin. We believe that further enhancement of the lead compound can serve as a potential inhibitor against resistance enzyme in drug discovery process. J. Cell. Biochem. 118: 4088-4094, 2017. © 2017 Wiley Periodicals, Inc.

Investigation of the bindings of a class of inhibitors with GSK3β kinase using thermodynamic integration MD simulation and kinase assay

GSK3β kinase is a noteworthy target for discovery of the drugs that will be used to treat several diseases. In the effort to identify a new inhibitor lead compound, we utilized thermodynamic integration (TI)-molecular dynamics (MD) simulation and kinase assay to investigate the bindings between GSK3β kinase and five compounds that were analogous to a known inhibitor with an available crystal structure. TI-MD simulations of the first two compounds (analogs 1 and 2) were used for calibration. The computed binding affinities of analogs 1 and 2 agreed well with the experimental results. The rest three compounds (analogs 3-5) were newly obtained from a database search, and their affinity data were newly measured in our labs. TI-MD simulations predicted the binding modes and the computed ΔΔG values have a reasonably good correlation with the experimental affinity data. These newly identified inhibitors appear to be new leads according to our survey of GSK3β inhibitors listed in recent review articles. The predicted binding modes of these compounds should aid in designing new derivatives of these compounds in the future.

A Molecular Docking and Dynamics Study to Screen Potent Anti-Staphylococcal Compounds Against Ceftaroline Resistant MRSA

World Health Organization reports that methicillin-resistant Staphylococcus aureus (MRSA) is the origin of higher proportion of hospital acquired infections. In order to combat the effect of MRSA infection, an ideal drug should stimulate the allosteric exposure of active site, prompting penicillin binding proteins (PBP2a) to bind with that particular compound. Ceftaroline shows high binding affinity towards PBP2a and also confers resistance against degrading enzymes. Recently, two amino acid alterations in the allosteric site of PBP2a, asparagine (N) to lysine (K) at position 146 and glutamic acid (E) to lysine at position 150 are reported to confer resistance against ceftaroline resulting in the rise of ceftaroline-resistant MRSA strains. The present study focuses on the identification of potential ligands that can effectively bind with allosteric site of PBP2a, that leads to the access of active site and entry of a β-lactam antibiotic for effective inhibition. The results obtained from our study will be useful for designing effective compounds with potential therapeutic effects against ceftaroline resistant MRSA strains.

Ethyl 4-(4-methylphenyl)-4-pentenoate from Vetiveria zizanioides Inhibits Dengue NS2B-NS3 Protease and Prevents Viral Assembly: A Computational Molecular Dynamics and Docking Study

Around 50 % of the world's population is at the risk of dengue, a viral infection. Presently, there are not many drugs and prophylactic measures available to control dengue viral infection, and hence, there is an urgent need to develop effective antidengue compound from natural sources. In the current study, we explored the antiviral properties of the medicinal plant Vetiveria zizanioides against dengue virus. Initially, the antiviral properties of active compounds were examined using docking analysis along with reference ligand. The enzyme-ligand complex which showed higher binding affinity than the reference ligand was employed for subsequent analysis. The stability of the top scoring enzyme-ligand complex was further validated using molecular simulation studies. On the whole, the study reveals that the compound Ethyl 4-(4-methylphenyl)-4-pentenoate has an effective antiviral property, which can serve as a potential lead molecule in drug discovery process.

Molecular dynamics simulation studies of GSK-3β ATP competitive inhibitors: understanding the factors contributing to selectivity

Glycogen synthase kinase-3 is a constitutively acting, multifunctional serine threonine kinase, the role of which has been implicated in several physiological pathways and has emerged as a promising target for the treatment of type-II diabetes and Alzheimer's disease. In order to provide a detailed understanding of the origin of selectivity determinants of ATP competitive inhibitors, molecular dynamics simulations in combination with MM-PBSA binding energy calculations were performed using crystal structures of GSK-3β and CDK-2 in complex with 12 ATP competitive inhibitors. Analysis of energy contributions indicate that electrostatic interaction energy dictates the selectivity of ATP competitive inhibitors against CDK-2. Key interactions as well as residues that potentially make a major contribution to the binding free energy were identified at the ATP binding site. This analysis stresses the need for the inhibitors to interact with Lys85, Thr138, and Arg141 in the binding site of GSK-3β to show selectivity. The residue-wise energy decomposition analysis further suggested the additional role of Gln185 in determining the selectivity of maleimides. The results obtained in this study can be utilized to design new selective GSK-3 ATP competitive inhibitors.

Investigation of differences in the binding affinities of two analogous ligands for untagged and tagged p38 kinase using thermodynamic integration MD simulation

Thermodynamic integration (TI) molecular dynamics (MD) simulations for the binding of a pair of a reference ("ref") ligand and an analogous ("analog") ligand to either tagged (with six extra residues at the N-terminus) or untagged p38 kinase proteins were carried out in order to probe how the binding affinity is influenced by the presence or absence of the peptide tag in p38 kinase. This possible effect of protein length on the binding affinity of a ligand-which is seldom addressed in the literature-is important because, even when two labs claim to have performed experiments with the same protein, they may actually have studied variants of the same protein with different lengths because they applied different protein expression conditions/procedures. Thus, if we wanted to compare ligand binding affinities measured in the two labs, it would be necessary to account for any variation in ligand binding affinity with protein length. The pair of ligand-p38 kinase complexes examined in this work (pdb codes: 3d7z and 3lhj, respectively) were ideal for investigating this effect. The experimentally determined binding energy for the ref ligand with the untagged p38 kinase was -10.9 kcal mol(-1), while that for the analog ligand with the tagged p38 kinase was -11.9 kcal mol(-1). The present TI-MD simulation of the mutation of the ref ligand into the analog ligand while the ligand is bound to the untagged p38 kinase predicted that the binding affinity of the analog ligand is 2.0 kcal mol(-1) greater than that of the ref ligand. A similar simulation also indicated that the same was true for ligand binding to the tagged protein, but in this case the binding affinity for the analog ligand is 2.5 kcal mol(-1) larger than that for the ref ligand. These results therefore suggest that the presence of the peptide tag on p38 kinase increased the difference in the binding energies of the ligands by a small amount of 0.5 kcal mol(-1). This result supports the assumption that the presence of a peptide tag has only a minor effect on ΔG values. The error bars in the computed ΔG values were then estimated via confidence interval analysis and a time autocorrelation function for the quantity dV/dλ. The estimated correlation time was ~0.5 ps and the error bar in the ΔG values estimated using nanosecond-scale simulations was ±0.3 kcal mol(-1) at a confidence level of 95%. These predicted results can be verified in future experiments and should prove useful in subsequent similar studies. Graphical Abstract Thermodynamic cycles for binding of two analogous ligands with untagged and tagged p38 kinases and associated Gibbs free energy.

Characterization of the Interaction between Gallic Acid and Lysozyme by Molecular Dynamics Simulation and Optical Spectroscopy

The binding interaction between gallic acid (GA) and lysozyme (LYS) was investigated and compared by molecular dynamics (MD) simulation and spectral techniques. The results from spectroscopy indicate that GA binds to LYS to generate a static complex. The binding constants and thermodynamic parameters were calculated. MD simulation revealed that the main driving forces for GA binding to LYS are hydrogen bonding and hydrophobic interactions. The root-mean-square deviation verified that GA and LYS bind to form a stable complex, while the root-mean-square fluctuation results showed that the stability of the GA-LYS complex at 298 K was higher than that at 310 K. The calculated free binding energies from the molecular mechanics/Poisson-Boltzmann surface area method showed that van der Waals forces and electrostatic interactions are the predominant intermolecular forces. The MD simulation was consistent with the spectral experiments. This study provides a reference for future study of the pharmacological mechanism of GA.