Molecular mechanisms of allosteric probe dependence in μ opioid receptor

Allostery is one of the most important features of proteins. It greatly contributes to the complexity of life, since it enables possibility of precise tuning of protein function, as well as performing more than one function per protein. Probe dependence is one of the unique features of allostery. It allows a protein to respond differently to the same allosteric modulator when different drugs or transmitters are bound. Unfortunately, allosteric mechanisms are difficult to investigate experimentally. Instead, they can be reproduced artificially in simulations. We simulated in silico a native-like cell membrane fragment with an active-state human μ opioid receptor (MOR) in order to investigate diverse effects of a receptor’s positive allosteric modulator on various agonists. Particular emphasis on native-likeness of the environment was put. We managed to reproduce the experimentally observed effects, which allowed us to take deeper insight into their underlying mechanisms. We found an allosteric pathway in the receptor, leading from the ligand binding site to the intracellular, effector site. We observed that the modulator affected the pathway, inducing different resultant responses for full and partial agonists.     

Molecular dynamics and combined docking studies for the identification of Zaire ebola virus inhibitors

Ebola virus (EBOV) is a lethal human pathogen with a risk of global spread of its zoonotic infections, and Ebolavirus Zaire specifically has the highest fatality rate amongst other species. There is a need for continuous effort towards having therapies, as a single licensed treatment to neutralize the EBOV is yet to come into reality. This present study virtually screened the MCULE database containing almost 36 million compounds against the structure of a Zaire Ebola viral protein (VP) 35 and a consensus scoring of both MCULE and CLCDDW docking programs remarked five compounds as potential hits. These compounds, with binding energies ranging from –7.9 to –8.9?kcal/mol, were assessed for predictions of their physicochemical and bioactivity properties, as well as absorption, distribution, metabolism, excretion, and toxicity (ADMET) criteria. The results of the 50?ns molecular dynamics simulations showed the presence of dynamic stability between ligand and protein complexes, and the structures remained significantly unchanged at the ligand-binding site throughout the simulation period. Both docking analysis and molecular dynamics simulation studies suggested strong binding affinity towards the receptor cavity and these selected compounds as potential inhibitors against the Zaire Ebola VP 35. With respect to inhibition constant values, bioavailability radar and other physicochemical properties, compound A (MCULE-1018045960-0-1) appeared to be the most promising hit compound. However, the ligand efficiency and ligand efficiency scale need improvement during optimization, and also validation via in vitro and in vivo studies are necessary to finally make a lead compound in treating Ebola virus diseases.

Communicated by Ramaswamy H. Sarma     

Diffusion of energetic compounds through biological membrane: Application of classical MD and COSMOmic approximations

Computational studies of the potential biological impact of several energetic compounds were performed. The most commonly used explosives were considered in the present studies: trinitrotoluene (TNT), 2,4-dinitrotoluene (2,4-DNT), 2,4-dinitroanisole (DNAN), and 5-Nitro-2,4-dihydro-3H-1,2,4-triazol-3-one (NTO). The effect of such factors as ionic strength and presence of DMSO in the water solution on the structure of the membrane were considered using the POPC lipid bilayer as an example. Molecular dynamics (MD) simulations revealed that, even on a short-time scale, the influence of those additives is noticeable, and therefore those factors should always be taken into account. The MD and the COSMOmic approaches were used to elucidate the ability of the energetic compounds to penetrate the living cell. Calculated free energy profiles and partitioning coefficients revealed distributions of the compounds in the lipid bilayer as well as an overall ability to enter the cell. MD in this case provides a better representation of the free energy profile, while the COSMOmic approach works better to predict log(K_lipw) values. The effect of the functional group was observed for the profiles that were obtained using the MD method.     

Molecular Modelling reveals the inhibition mechanism and structure–activity relationship of curcumin and its analogues to Staphylococcal aureus Sortase A

Previous studies found that the activity of Sortase A, a bacterial surface protein from Staphylococcus aureus, was inhibited by curcumin and its analogues. To explore this inhibitory mechanism, Sortase A and its inhibitors in complex systems were studied by molecular docking, molecular modelling, binding energy decomposition calculation and steered molecular dynamics simulations. Energy decomposition analysis indicated that PRO-163, LEU-169, GLN-172, ILE-182 and ILE-199 are key residues in Sortase A-inhibitor complexes. Furthermore, interactions between the methoxyl group on the benzene ring in the conjugated molecule (curcumin, demethoxycurcumin, bisdemethoxycurcumin) and VAL-168, LEU-169 and GLN-172 induce the inhibitory activity based on the energy decomposition and distance analyses between the whole residues and inhibitors. However, because of its coiled structure, the non-conjugated molecule, tetrahydrocurcumin, with key residues in the binding sites of Sortase A, interacted weakly with SrtA, leading to the loss of inhibitory activity. Based on these results, the methoxyl group on the benzene ring in the conjugated molecule largely influenced the inhibitory activity of the Sortase A inhibitors.     

Identification of isobavachalcone as a potential drug for rice blast disease caused by the fungus Magnaporthe grisea

The rice blast disease caused by the fungus Magnaporthe grisea is one of the most devastating rice diseases, but there is no effective fungicide toward chitinase which is a key enzyme of M. grisea. In this study, we observed that distortion and cell-wall damage of M. grisea hyphae were significantly under the scanning electron micrograph after a 24-h treatment with 10?mg/L isobavachalcone (IBC) extracted from Psoralea corylifolia L. To further explore the effect of IBC on the cell wall of M. grisea, we examined changes in enzymes associated with cell wall degradation by enzyme activity experiments, treated liquid culture mycelia with 10?mg/L IBC for 1?h. Results displayed that chitinase was obviously more active than control group. To illustrate the interactions between IBC and chitinase, the studies of homology modeling and molecular docking were carried out successively. The results revealed that IBC had hydrogen bonds with residues ASP267 and ARG276 of chitinase. Besides, these nonpolar residues TYR270, PRO271, VAL272, LEU310, PRO311, TYR316, and LEU317 were able to form strong hydrophobic interactions. Binding energies of the chitinase-IBC complexes were calculated by MM-GBSA showed that the ΔG_bind score of molecular dynamics had lower binding energy and more stable than docking complexes. All above, IBC owns significant agonistic activity in chitinase and would be a potent fungicide to inhibit the growth of M. grisea. We hope the above information provides an important insight for understanding the interactions between IBC and chitinase, which may be useful in the discovery of a novel potent agonist.

Communicated by Ramaswamy H. Sarma     

Simulation of protein diffusion: a sensitive probe of protein–solvent interactions

Aqueous solutions of Candida antarctica lipase B (CALB) were simulated considering three different water models (SPC/E, TIP3P, TIP4P) by a series of molecular dynamics (MD) simulations of three different box sizes (L = 9, 14, and 19 nm) to determine the diffusion coefficient, the water viscosity and the protein density. The protein–water systems were equilibrated for 500 ns, followed by 100 ns production runs which were analysed. The diffusional properties of CALB were characterized by the Stokes radius (R_S), which was derived from the diffusion coefficient and the viscosity. R_S was compared to the geometric radius (R_G) of CALB, which was derived from the protein density. R_S and R_G differed by 0.27 nm for SPC/E and by 0.40 and 0.39 nm for TIP3P and TIP4P, respectively, which characterizes the thickness of the diffusive hydration layer on the protein surface. The simulated hydration layer of CALB resulted in agreement with those experimentally determined for other seven different proteins of comparable size. By avoiding the most common pitfalls, protein diffusion can be reliably simulated: simulating different box sizes to account for the finite size effect, equilibrating the protein–water system sufficiently, and using the complete production run for the determination of the diffusion coefficient.     

Probing light chain mutation effects on thrombin via molecular dynamics simulations and machine learning

Thrombin is a key component for chemotherapeutic and antithrombotic therapy development. As the physiologic and pathologic roles of the light chain still remain vague, here, we continue previous efforts to understand the impacts of the disease-associated single deletion of LYS9 in the light chain. By combining supervised and unsupervised machine learning methodologies and more traditional structural analyses on data from 10 μs molecular dynamics simulations, we show that the conformational ensemble of the ΔK9 mutant is significantly perturbed. Our analyses consistently indicate that LYS9 deletion destabilizes both the catalytic cleft and regulatory functional regions and result in some conformational changes that occur in tens to hundreds of nanosecond scaled motions. We also reveal that the two forms of thrombin each prefer a distinct binding mode of a Na⁺ ion. We expand our understanding of previous experimental observations and shed light on the mechanisms of the LYS9 deletion associated bleeding disorder by providing consistent but more quantitative and detailed structural analyses than early studies in literature. With a novel application of supervised learning, i.e. the decision tree learning on the hydrogen bonding features in the wild-type and ΔK9 mutant forms of thrombin, we predict that seven pairs of critical hydrogen bonding interactions are significant for establishing distinct behaviors of wild-type thrombin and its ΔK9 mutant form. Our calculations indicate the LYS9 in the light chain has both localized and long-range allosteric effects on thrombin, supporting the opinion that light chain has an important role as an allosteric effector.     

Exploring the structural determinants of novel xanthine derivatives as A2B adenosine receptor antagonists: a computational study

Adenosine is a ubiquitous endogenous nucleoside that controls numerous physiological functions via interacting with its specific G-coupled receptors. Activation of adenosine receptors (AdoRs), particularly A_2B AdoRs promotes the release of inflammatory cytokines; reduces vascular permeabilization and induces angiogenesis, thereby making A_2B AdoR becomes a potentially pharmacological target for drug development. Presently, for investigating the structural determinants of 164 xanthine derivatives as A_2B AdoR antagonists, we performed an in silico study integrating with 3D-QSAR, docking and molecular dynamics (MD) simulation. The obtained optimal model shows strong predictability (Q²?=?0.647, R²_ncv?=?0.955, and R²_pred?=?0.848). Additionally, to explore the binding mode of the ligand with A_2B AdoR and to understand their binding mechanism, docking analysis, MD simulations (20?ns), and the calculation of binding free energy were also carried out. Finally, the structural determinants of these xanthine derivatives were identified and a total of 20 novel A_2B AdoR antagonists with improved potency were computationally designed, and their synthetic feasibility and selectivity were also evaluated. The information derived from the present study offers a better appreciation for exploring the interaction mechanism of the ligand with A_2B AdoR, which could be helpful for designing novel potent A_2B AdoR antagonists.

Communicated by Ramaswamy H. Sarma