Investigation of the molecular interactions of triclocarban with human serum albumin using multispectroscopies and molecular modeling

Triclocarban (TCC), as a broad spectrum antibacterial agent widely used in personal care products, has recently been recognized as environmental pollutant with the potential of adversely affecting wildlife and human health. However, the behavior of TCC in blood circulatory system and the potential toxicity of TCC at the molecular level have been poorly investigated. In this study, the effect of TCC on human serum albumin (HSA) and the binding mechanism of TCC to HSA were examined using spectroscopic techniques and molecular modeling methods. The fluorescence results suggested that the fluorescence of HSA was quenched by TCC through a static quenching mechanism and nonradiation energy transfer, and TCC was bound to HSA with moderately strong binding affinity via hydrophobic interaction based on the analysis of the thermodynamic parameters. The site marker competitive experiments revealed that TCC bound into subdomain IIA (site I) of HSA. In addition, the results obtained from the circular dichroism, Fourier transform infrared (FT-IR), 8-anilino-1-naphthalenesulfonic acid fluorescence, synchronous fluorescence, three-dimensional fluorescence spectra and dynamic light scattering suggested the change in the microenvironment and conformation of HSA during the binding reaction. Finally, the best binding mode of TCC and specific interaction of TCC with amino acid residues were determined using molecular docking and molecular dynamics simulations. In a word, the present studies can provide a way to help us well understand the transport, distribution and toxicity effect of TCC when it diffused in the human body.

Communicated by Ramaswamy H. Sarma     

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.     

In vitro and in silico evaluation of fucosterol from Sargassum horridum as potential human acetylcholinesterase inhibitor

The fucosterol has been reported numerous biological activities. In this study, the activity in vitro of the fucosterol from Sargassum horridum as potential human acetylcholinesterase inhibitor was evaluated. The structural identification was obtained by nuclear magnetic resonance (NMR) spectroscopy and based on experimental data, we combined docking and molecular dynamics simulations coupled to the molecular-mechanics-generalized-born-surface-area approach to evaluating the structural and energetic basis for the molecular recognition of fucosterol and neostigmine at the binding site of acetylcholinesterase (AChE). In addition, the Lineweaver–Burk plot showed the nature of a non-competitive inhibition. The maximum velocity (V_max) and the constant of Michaelis–Menten (K_m) estimated for fucosterol (0.006 µM) were 0.015 1/V_o (ΔA/h and 6.399 1/[ACh] mM^?1, respectively. While, for neostigmine (0.14 µM), the V_max was 0.022 1/V_o (ΔA/h) and K_m of 6.726 1/[ACh] mM^?1, these results showed a more effective inhibition by fucosterol respect to neostigmine. Structural analysis revealed that neostigmine reaches the AChE binding site reported elsewhere, whereas fucosterol can act as a no-competitive and competitive acetylcholinesterase inhibitor, in agree with kinetic enzymatic experiments. Binding free energy calculations revealed that fucosterol reaches the acetylcholinesterase binding site with higher affinity than neostigmine, which is according to experimental results. Whereas the per-residue decomposition free energy analysis let us identify crucial residues involved in the molecular recognition of ligands by AChE. Results corroborate the ability of theoretical methods to provide crucial information at the atomic level about energetic and structural differences in the binding interaction and affinity from fucosterol with AChE.

Communicated by Ramaswamy H. Sarma     

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.     

A combination of bioactive and nonbioactive alkyl-peptides form a more stable nanofiber structure for differentiating neural stem cells: a molecular dynamics simulation survey

Self-assembling alkyl-peptides are important molecules due to their ability to construct nano-level structures such as nanofibers to be utilized as tissue engineering scaffolds. The bioactive epitope of FAQRVPP which acts as neural stem cells (NSCs) outgrowth inducing factor is used in nanofiber structures. Based on previous experimental studies the density and distribution pattern of the epitopes on the surface of the nanofibers plays an important role in the differentiation function efficiency. We decided to survey and compare the stability of two pre-constructed fiber structures in the forms of all-functionalized nanofiber (containing only bioactive alkyl-peptides) and distributed functionalized nanofiber (a combination of nonbioactive and bioactive alkyl-peptides with ratio 2:1). Our findings reveal that the all-functionalized fiber shows an unstable structure and is split into intermediate micelle-like structures to reduce compactness and steric hindrance of functional epitopes whereas the distributed functionalized fiber shows an integrated stable nanofiber with a more amount of beta sheets that are well-organized and oriented around the hydrophobic core. The hydrogen bonds and energy profiles of the structures indicate the role of hydrophobic interactions during the alkyl-chain core formation and the important role of electrostatic interactions and hydrogen bond network in the stability of the final structures. Finally, it seems that the possibility of the presence of intermediate structure is increased in the all-functionalized nanofiber environment, and it can reduce functional efficiency of the scaffolds. These findings can help to design more efficient nanofiber structures with different goals in scaffolds for tissue engineering. Abbreviations MD Molecular Dynamics

NSCs Neural Stem Cells

PME Particle mesh Ewald

RDF Radial Distribution Function

RG Radius of gyration

RASA Relative Accessible Surface Area

RMSD Root Mean Square Deviations

SASA Solvent Accessible Surface Area.

Communicated by Ramaswamy H. Sarma