Insights into conformational regulation of PfMATE transporter from Pyrococcus furiosus induced by alternating protonation state of Asp41 residue: A molecular dynamics simulation study.

BACKGROUND: Multidrug and toxic compound extrusion (MATE) family transporters induce multiple-drug resistance (MDR) of bacterial pathogens and cancer cells, thus causing critical reductions in the therapeutic efficacies of antibiotics and anti-cancer drugs. Unfortunately, to date, the details and intrinsic reason about conformational regulation mechanism of MATE transporters remain elusive. METHOD: In this work, molecular dynamics (MD) simulations were conducted to explore the conformational regulation mechanism of PfMATE transporter from Pyrococcus furiosus based on different protonation state of Asp41. Two (MD) simulation systems were investigated: a system with protonation of Asp41 and a system without protonation of Asp41, which were named by D184(H)D41(H) system and D184(H) system, respectively. RESULTS AND CONCLUSIONS: Firstly, MD simulation results indicate that conformational changes mainly happen in extracellular regions of PfMATE protein. Further analysis reveals that PfMATE protein experiences different motion mode and forms different conformation based on different protonation state of Asp41. In the D184(H)D41(H) system, PfMATE experiences an opening motion and forms a more outward-open conformation. As for the D184(H) system, the protein has an anticlockwise rotational motion with the channel axis of protein and the more outward-open conformation does not appear. It can be inferred that protonation of Asp41 is essential for conformational regulation of PfMATE during transporting substrates. GENERAL SIGNIFICANCE: These findings provide intrinsic information for understanding the conformational regulation mechanism of PfMATE and will be very meaningful to explore the MDR mechanism of PfMATE further.

Search for ß2 adrenergic receptor ligands by virtual screening via grid computing and investigation of binding modes by docking and molecular dynamics simulations.

We designed a program called MolGridCal that can be used to screen small molecule database in grid computing on basis of JPPF grid environment. Based on MolGridCal program, we proposed an integrated strategy for virtual screening and binding mode investigation by combining molecular docking, molecular dynamics (MD) simulations and free energy calculations. To test the effectiveness of MolGridCal, we screened potential ligands for ß2 adrenergic receptor (ß2AR) from a database containing 50,000 small molecules. MolGridCal can not only send tasks to the grid server automatically, but also can distribute tasks using the screensaver function. As for the results of virtual screening, the known agonist BI-167107 of ß2AR is ranked among the top 2% of the screened candidates, indicating MolGridCal program can give reasonable results. To further study the binding mode and refine the results of MolGridCal, more accurate docking and scoring methods are used to estimate the binding affinity for the top three molecules (agonist BI-167107, neutral antagonist alprenolol and inverse agonist ICI 118,551). The results indicate agonist BI-167107 has the best binding affinity. MD simulation and free energy calculation are employed to investigate the dynamic interaction mechanism between the ligands and ß2AR. The results show that the agonist BI-167107 also has the lowest binding free energy. This study can provide a new way to perform virtual screening effectively through integrating molecular docking based on grid computing, MD simulations and free energy calculations. The source codes of MolGridCal are freely available at http://molgridcal.codeplex.com.

Ligand induced change of ß2 adrenergic receptor from active to inactive conformation and its implication for the closed/open state of the water channel: insight from molecular dynamics simulation, free energy calculation and Markov state model analysis.

The reported crystal structures of ß2 adrenergic receptor (ß2AR) reveal that the open and closed states of the water channel are correlated with the inactive and active conformations of ß2AR. However, more details about the process by which the water channel states are affected by the active to inactive conformational change of ß2AR remain illusive. In this work, molecular dynamics simulations are performed to study the dynamical inactive and active conformational change of ß2AR induced by inverse agonist ICI 118,551. Markov state model analysis and free energy calculation are employed to explore the open and close states of the water channel. The simulation results show that inverse agonist ICI 118,551 can induce water channel opening during the conformational transition of ß2AR. Markov state model (MSM) analysis proves that the energy contour can be divided into seven states. States S1, S2 and S5, which represent the active conformation of ß2AR, show that the water channel is in the closed state, while states S4 and S6, which correspond to the intermediate state conformation of ß2AR, indicate the water channel opens gradually. State S7, which represents the inactive structure of ß2AR, corresponds to the full open state of the water channel. The opening mechanism of the water channel is involved in the ligand-induced conformational change of ß2AR. These results can provide useful information for understanding the opening mechanism of the water channel and will be useful for the rational design of potent inverse agonists of ß2AR.

Semisynthesis and quantitative structure-activity relationship (QSAR) study of some cholesterol-based hydrazone derivatives as insecticidal agents.

In continuation of our program aimed at the discovery and development of natural-product-based insecticidal agents, four series of novel cholesterol-based hydrazone derivatives were synthesized, and their insecticidal activity was tested against the pre-third-instar larvae of oriental armyworm, Mythimna separata (Walker) in vivo at 1mg/mL. All the derivatives showed the better insecticidal activity than their precursor cholesterol. Quantitative structure-activity relationship (QSAR) model demonstrated that six descriptors such as RDF085v, Mor06u, Mor11u, Dv, HATS0v and H-046, are likely to influence the insecticidal activity of these compounds. Among them, two important ones are the Mor06u and RDF085v.

Synthesis and quantitative structure-activity relationship (QSAR) study of novel N-arylsulfonyl-3-acylindole arylcarbonyl hydrazone derivatives as nematicidal agents.

In continuation of our program aimed at the discovery and development of natural-product-based pesticidal agents, 54 novel N-arylsulfonyl-3-acylindole arylcarbonyl hydrazone derivatives were prepared, and their structures were well characterized by ¹H NMR, ¹³C NMR, HRMS, ESI-MS, and mp. Their nematicidal activity was evaluated against that of the pine wood nematode, Bursaphelenchus xylophilus in vivo. Among all of the derivatives, especially V-12 and V-39 displayed the best promising nematicidal activity with LC50 values of 1.0969 and 1.2632 mg/L, respectively. This suggested that introduction of R¹ and R² together as the electron-withdrawing substituents, R³ as the methyl group, and R4 as the phenyl with the electron-donating substituents could be taken into account for further preparation of these kinds of compounds as nematicidal agents. Six selected descriptors are a WHIM descriptor (E1m), two GETAWAY descriptors (R1m+ and R3m+), a Burden eigenvalues descriptor (BEHm8), and two edge-adjacency index descriptors (EEig05x and EEig13d). Quantitative structure-activity relationship (QSAR) studies demonstrated that the structural factors, such as molecular mass (a negative correlation with the bioactivity) and molecular polarity (a positive correlation with bioactivity), are likely to govern the nematicidal activities of these compounds. For this model, the correlation coefficient (R²(training set)), the leave-one-out cross-validation correlation coefficient (Q²(LOO)), and the 7-fold cross-validation correlation coefficient (Q²(7-fold)) were 0.791, 0.701, and 0.715, respectively. The external cross-validation correlation coefficient (Q²ext) and the root-mean-square error for the test set (RMSE(test set)) were 0.774 and 3.412, respectively. This study will pave the way for future design, structural modification, and development of indole derivatives as nematicidal agents.

Synthesis and quantitative structure-activity relationship (QSAR) study of novel 4-acyloxypodophyllotoxin derivatives modified in the A and C rings as insecticidal agents.

In continuation of our program aimed at the discovery and development of natural-product-based insecticidal agents, we have synthesized three series of novel 4-acyloxy compounds derived from podophyllotoxin modified in the A and C rings, which is isolated as the main secondary metabolite from the roots and rhizomes of Podophyllum hexandrum . Their insecticidal activity was preliminarily evaluated against the pre-third-instar larvae of Mythimna separata in vivo. Compound 9g displayed the best promising insecticidal activity. It revealed that cleavage of the 6,7-methylenedioxy group of podophyllotoxin will lead to a less active compound and that the C-4 position of podophyllotoxin was the important modification location. A quantitative structure-activity relationship (QSAR) model was developed by genetic algorithm combined with multiple linear regression (GA-MLR). For this model, the squared correlation coefficient (R(2)) is 0.914, the leave-one-out cross-validation correlation coefficient (Q(2)(LOO)) is 0.881, and the root-mean-square error (RMSE) is 0.024. Five descriptors, BEHm2, Mor14v, Wap, G1v, and RDF020e, are likely to influence the biological activity of these compounds. Among them, two important ones are BEHm2 and Mor14v. This study will pave the way for further design, structural modification, and development of podophyllotoxin derivatives as insecticidal agents.

Synthesis and quantitative structure-activity relationship (QSAR) study of novel isoxazoline and oxime derivatives of podophyllotoxin as insecticidal agents.

In continuation of our program aimed at the discovery and development of natural-product-based insecticidal agents, 33 isoxazoline and oxime derivatives of podophyllotoxin modified in the C and D rings were synthesized and their structures were characterized by Proton nuclear magnetic resonance ((1)H NMR), high-resolution mass spectrometry (HRMS), electrospray ionization-mass spectrometry (ESI-MS), optical rotation, melting point (mp), and infrared (IR) spectroscopy. The stereochemical configurations of compounds 5e, 5f, and 9f were unambiguously determined by X-ray crystallography. Their insecticidal activity was evaluated against the pre-third-instar larvae of northern armyworm, Mythimna separata (Walker), in vivo. Compounds 5e, 9c, 11g, and 11h especially exhibited more promising insecticidal activity than toosendanin, a commercial botanical insecticide extracted from Melia azedarach . A genetic algorithm combined with multiple linear regression (GA-MLR) calculation is performed by the MOBY DIGS package. Five selected descriptors are as follows: one two-dimensional (2D) autocorrelation descriptor (GATS4e), one edge adjacency indice (EEig06x), one RDF descriptor (RDF080v), one three-dimensional (3D) MoRSE descriptor (Mor09v), and one atom-centered fragment (H-052) descriptor. Quantitative structure-activity relationship studies demonstrated that the insecticidal activity of these compounds was mainly influenced by many factors, such as electronic distribution, steric factors, etc. For this model, the standard deviation error in prediction (SDEP) is 0.0592, the correlation coefficient (R(2)) is 0.861, and the leave-one-out cross-validation correlation coefficient (Q(2)loo) is 0.797.