We utilized molecular dynamics (MD) simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) free energy calculations to investigate the specificity of two oligonucleotide probes, namely probe B and probe D, in detecting single-stranded DNA (ssDNA) within three bacteria families: Enterobacteriaceae, Pasteurellaceae, and Vibrionaceae. Due to the limited understanding of molecular mechanisms in the previous research, we have extended the discussion to focus specifically on investigating the binding process of bacteria-probe DNA duplexes, with an emphasis on analyzing the binding free energy. The role of electrostatic contributions in the specificity between the oligonucleotide probes and the bacterial ssDNAs was investigated and found to be crucial. Our calculations yielded results that were highly consistent with the experimental data. Through our study, we have successfully exhibited the benefits of utilizing in-silico approaches as a powerful virtual-screening tool, particularly in research areas that demand a thorough comprehension of molecular interactions.
Tyrosinase is a key enzyme in melanogenesis. Generally, mushroom tyrosinase from A. bisporus had been used as a model in skin-whitening agent tests employed in the cosmetic industry. The recently obtained crystal structure of bacterial tyrosinase from B. megaterium has high similarity (33.5%) to the human enzyme and thus it was used as a template for constructing of the human model. Binding of tyrosinase to a series of its inhibitors was simulated by automated docking calculations. Docking and MD simulation results suggested that N81, N260, H263, and M280 are involved in the binding of inhibitors to mushroom tyrosinase. E195 and H208 are important residues in bacterial tyrosinase, while E230, S245, N249, H252, V262, and S265 bind to inhibitors and are important in forming pi interaction in human tyrosinase.
Conformational search for the most stable geometry connection of 16 sets of polydopamine (PDA) tetramer subunits has been systematically investigated using density functional theory (DFT) calculations. Our results indicated that the more planar subunits are, the more stable they are. This finding is in good agreement with recent experimental observations, which have suggested that PDA are composed of the nearly planar subunits that appear to be stacked together via the π-π interactions to form graphite-like layered aggregates associated with the balance of the intramolecular hydrogen bonds and steric effects from the indole and catechol moieties. Molecular dynamics (MD) simulations of 16 spherical clusters of the tetramer subunits of PDA in the gas and aqueous phase were performed at 298 K and confirmed the stability of supramolecular tetramer aggregates. The complex formation and binding energy of all 16 clusters are very strong although the shapes of the clusters in aqueous solution are not spherical and are very much different from those in the gas phase. The aggregations of all 16 clusters in aqueous solution were also confirmed from the profiles of the Kratky plot and the radius of gyration of all clusters. Our MD results in both gas phase and aqueous solution pointed out that there are high possibilities of aggregations of the 16 kinds of tetramer subunits although the conformations of each tetramer subunit are not flat. In summary, this work brings an insight into the controversial structure of PDA tetramer units and explains some of the important structural features found in the aqueous phase in comparison to the gas phase.
The data is obtained from exploring the modulatory activities of bioflavonoids on P-glycoprotein function by ligand-based approaches. Multivariate Linear-QSAR models for predicting the induced/inhibitory activities of the flavonoids were created. Molecular descriptors were initially used as independent variables and a dependent variable was expressed as pFAR. The variables were then used in MLR analysis by stepwise regression calculation to build the linear QSAR data. The entire dataset consisted of 23 bioflavonoids was used as a training set. Regarding the obtained MLR QSAR model, R of 0.963, R (2)=0.927, [Formula: see text], SEE=0.197, F=33.849 and q (2)=0.927 were achieved. The true predictabilities of QSAR model were justified by evaluation with the external dataset (Table 4). The pFARs of representative flavonoids were predicted by MLR QSAR modelling. The data showed that internal and external validations may generate the same conclusion.
Extracellular signal-regulated kinases 1 and 2 (ERK1/2) play key roles in promoting cell survival and proliferation through the phosphorylation of various substrates. Remarkable antitumour activity is found in many inhibitors that act upstream of the ERK pathway. However, drug-resistant tumour cells invariably emerge after their use due to the reactivation of ERK1/2 signalling. ERK1/2 inhibitors have shown clinical efficacy as a therapeutic strategy for the treatment of tumours with mitogen-activated protein kinase (MAPK) upstream target mutations. These inhibitors may be used as a possible strategy to overcome acquired resistance to MAPK inhibitors. Here, we report a class of repeat proteins-designed ankyrin repeat protein (DARPin) macromolecules targeting ERK2 as inhibitors. The structural basis of ERK2-DARPin interactions based on molecular dynamics (MD) simulations was studied. The information was then used to predict stabilizing mutations employing a web-based algorithm, MAESTRO. To evaluate whether these design strategies were successfully deployed, we performed all-atom, explicit-solvent molecular dynamics (MD) simulations. Two mutations, Ala → Asp and Ser → Leu, were found to perform better than the original sequence (DARPin E40) based on the associated energy and key residues involved in protein-protein interaction. MD simulations and analysis of the data obtained on these mutations supported our predictions.
In this work, molecular docking, pharmacophore modeling and molecular dynamics (MD) simulation were rendered for the mouse P-glycoprotein (P-gp) (code: 4Q9H) and bioflavonoids; amorphigenin, chrysin, epigallocatechin, formononetin and rotenone including a positive control; verapamil to identify protein-ligand interaction features including binding affinities, interaction characteristics, hot-spot amino acid residues and complex stabilities. These flavonoids occupied the same binding site with high binding affinities and shared the same key residues for their binding interactions and the binding region of the flavonoids was revealed that overlapped the ATP binding region with hydrophobic and hydrophilic interactions suggesting a competitive inhibition mechanism of the compounds. Root mean square deviations (RMSDs) analysis of MD trajectories of the protein-ligand complexes and NBD2 residues, and ligands pointed out these residues were stable throughout the duration of MD simulations. Thus, the applied preliminary structure-based molecular modeling approach of interactions between NBD2 and flavonoids may be gainful to realize the intimate inhibition mechanism of P-gp at NBD2 level and on the basis of the obtained data, it can be concluded that these bioflavonoids have the potential to cause herb-drug interactions or be used as lead molecules for the inhibition of P-gp (as anti-multidrug resistance agents) via the NBD2 blocking mechanism in future.
Phosphodiesterase 1B (PDE1B) and PDE10A are dual-specificity PDEs that hydrolyse both cyclic adenosine monophosphate and cyclic guanosine monophosphate, and are highly expressed in the striatum. Several reports have suggested that PDE10A inhibitors may present a promising approach for the treatment of positive symptoms of schizophrenia, whereas PDE1B inhibitors may present a novel mechanism to modulate cognitive deficits. Previously, we have reported a novel dual inhibitor of PDE1B and PDE10A, compound 2 [(3-fluorophenyl)(2-methyl-2,3-dihydro-4H-benzo[b][1,4]oxazin-4-yl)methanone] which has shown inhibitory activity for human recombinant PDE1B and PDE10A in vitro. In the present study, the safety profile of compound 2 has been evaluated in rats in the acute oral toxicity study, as well as; the antipsychotic-like effects in the rat model of schizophrenia. Compound 2 was tolerated up to 1 g/kg when administered at a single oral dose. Additionally, compound 2 has strongly suppressed ketamine-induced hyperlocomotion, which presented a model for the positive symptoms of schizophrenia. It has also shown an ability to attenuate social isolation induced by chronic administration of ketamine and enhanced recognition memory of rats in the novel object recognition test. Altogether, our results suggest that compound 2 represents a promising therapy for the treatment of the three symptomatic domains of schizophrenia.
Dengue virus (DENV) infection causes mild to severe illness in humans that can lead to fatality in severe cases. Currently, no specific drug is available for the treatment of DENV infection. Thus, the development of an anti-DENV drug is urgently required. Cordycepin (3'-deoxyadenosine), which is a major bioactive compound in Cordyceps (ascomycete) fungus that has been used for centuries in Chinese traditional medicine, was reported to exhibit antiviral activity. However, the anti-DENV activity of cordycepin is unknown. We hypothesized that cordycepin exerts anti-DENV activity and that, as an adenosine derivative, it inhibits DENV replication. To test this hypothesis, we investigated the anti-DENV activity of cordycepin in DENV-infected Vero cells. Cordycepin treatment significantly decreased DENV protein at a half-maximal effective concentration (EC50) of 26.94 μM. Moreover, DENV RNA was dramatically decreased in cordycepin-treated Vero cells, indicating its effectiveness in inhibiting viral RNA replication. Via in silico molecular docking, the binding of cordycepin to DENV non-structural protein 5 (NS5), which is an important enzyme for RNA synthesis, at both the methyltransferase (MTase) and RNA-dependent RNA polymerase (RdRp) domains, was predicted. The results of this study demonstrate that cordycepin is able to inhibit DENV replication, which portends its potential as an anti-dengue therapy.
Protein-protein interaction plays an essential role in almost all cellular processes and biological functions. Coupling molecular dynamics (MD) simulations and nanoparticle tracking analysis (NTA) assay offered a simple, rapid, and direct approach in monitoring the protein-protein binding process and predicting the binding affinity. Our case study of designed ankyrin repeats proteins (DARPins)-AnkGAG1D4 and the single point mutated AnkGAG1D4-Y56A for HIV-1 capsid protein (CA) were investigated. As reported, AnkGAG1D4 bound with CA for inhibitory activity; however, it lost its inhibitory strength when tyrosine at residue 56 AnkGAG1D4, the most key residue was replaced by alanine (AnkGAG1D4-Y56A). Through NTA, the binding of DARPins and CA was measured by monitoring the increment of the hydrodynamic radius of the AnkGAG1D4-gold conjugated nanoparticles (AnkGAG1D4-GNP) and AnkGAG1D4-Y56A-GNP upon interaction with CA in buffer solution. The size of the AnkGAG1D4-GNP increased when it interacted with CA but not AnkGAG1D4-Y56A-GNP. In addition, a much higher binding free energy (∆GB) of AnkGAG1D4-Y56A (-31 kcal/mol) obtained from MD further suggested affinity for CA completely reduced compared to AnkGAG1D4 (-60 kcal/mol). The possible mechanism of the protein-protein binding was explored in detail by decomposing the binding free energy for crucial residues identification and hydrogen bond analysis.