Homology modeling and explicit membrane molecular dynamics simulation to delineate the mode of binding of thiazolidinediones into FFAR1 and the mechanism of receptor activation

Free fatty acid receptor 1 (FFAR1) is a member of a previously characterized cluster of orphan G protein-coupled receptors (GPCRs). Later, this orphan receptor was identified as a target of medium- to long-chain free fatty acids in β-cells of the pancreas. Administration of FFAR1 agonists has been proved to potentiate glucose-stimulated insulin secretion from pancreatic β-cells. It was reported that some thiazolidinediones (TZDs), the best studied PPARγ agonists, are also able to stimulate FFAR1 in a dose-dependent manner. In the present study, a homology model of the human FFAR1 was constructed and inserted into a pre-equilibrated DPPC/TIP3P membrane system. This system was then simulated for 20 ns in complex with the FFAR1 agonist GW9085, as well as rosiglitazone and pioglitazone. We noticed that the salt bridge between Glu172 and Arg258 and the H bond between Glu145 and His153 could be responsible for the stabilization of the receptor in the inactive state. Moreover, we described for the first time the binding mode of TZDs in the binding site of FFAR1. The thiazolidinedione head forms a hydrogen bonding network with the critical polar residues in the binding site, Arg258 and Asn244, while the rest of the molecule is embedded into the receptor hydrophobic pocket. Based on this modeling study, we arrived at a proposal of the pharmacophore required for binding to both PPARγ and FFAR1. Insights gained from this investigation should provide future directions for the design of novel dual acting antidiabetic agents.     

Structural investigations into the binding mode of novel neolignans Cmp10 and Cmp19 microtubule stabilizers by in silico molecular docking,molecular dynamics,and binding free energy calculations

Microtubule stabilizers provide an important mode of treatment via mitotic cell arrest of cancer cells. Recently, we reported two novel neolignans derivatives Cmp10 and Cmp19 showing anticancer activity and working as microtubule stabilizers at micromolar concentrations. In this study, we have explored the binding site, mode of binding, and stabilization by two novel microtubule stabilizers Cmp10 and Cmp19 using in silico molecular docking, molecular dynamics (MD) simulation, and binding free energy calculations. Molecular docking studies were performed to explore the β-tubulin binding site of Cmp10 and Cmp19. Further, MD simulations were used to probe the β-tubulin stabilization mechanism by Cmp10 and Cmp19. Binding affinity was also compared for Cmp10 and Cmp19 using binding free energy calculations. Our docking results revealed that both the compounds bind at Ptxl binding site in β-tubulin. MD simulation studies showed that Cmp10 and Cmp19 binding stabilizes M-loop (Phe272-Val288) residues of β-tubulin and prevent its dynamics, leading to a better packing between α and β subunits from adjacent tubulin dimers. In addition, His229, Ser280 and Gln281, and Arg278, Thr276, and Ser232 were found to be the key amino acid residues forming H-bonds with Cmp10 and Cmp19, respectively. Consequently, binding free energy calculations indicated that Cmp10 (?113.655 kJ/mol) had better binding compared to Cmp19 (?95.216 kJ/mol). This study provides useful insight for better understanding of the binding mechanism of Cmp10 and Cmp19 and will be helpful in designing novel microtubule stabilizers.     

Influence of genetic architecture on contemporary local evolution in the soapberry bug,Jadera haematoloma: artificial selection on beak length

Little is known about the influence of genetic architecture on local adaptation. We investigated the genetic architecture of the rapid contemporary evolution of mouthparts, the flight polymorphism and life history traits in the soapberry bug Jadera haematoloma (Hemiptera) using laboratory selection. The mouthparts of these seed‐feeding bugs have adapted in 40–50 years by decreasing in length following novel natural selection induced by a host switch to the seeds of an introduced tree with smaller fruits than those of the native host vine. Laboratory selection on beak length in both an ancestral population feeding on the native host and a derived population feeding on the introduced host reveals genetic variance allowing a rapid response (heritabilities of 0.51–0.87) to selection for either longer or shorter beaks. This selection resulted in reverse evolution by restoring long beaks in the derived population and forward evolution by re‐creating short beaks in the ancestral bugs. There were strong genetic correlations (0.68–0.84) in both populations between beak lengths and the frequency of flight morphs, with short beaks associated with short wings. The results reveal a genetically interrelated set of adaptive multivariate traits including both beak length and flight morph. This suite of traits reflects host plant patchiness and seeding phenology. Weaker evidence suggests that egg mass and early egg production may be elements of the same suite. Reversible or forward evolution thus may occur in a broad set of genetically correlated multivariate traits undergoing rapid contemporary adaptation to altered local environments.     

Disentangling invasion processes in a dynamic shipping-boating network

The relative importance of multiple vectors to the initial establishment, spread and population dynamics of invasive species remains poorly understood. This study used molecular methods to clarify the roles of commercial shipping and recreational boating in the invasion by the cosmopolitan tunicate, Botryllus schlosseri. We evaluated (i) single vs. multiple introduction scenarios, (ii) the relative importance of shipping and boating to primary introductions, (iii) the interaction between these vectors for spread (i.e. the presence of a shipping-boating network) and (iv) the role of boating in determining population similarity. Tunicates were sampled from 26 populations along the Nova Scotia, Canada, coast that were exposed to either shipping (i.e. ports) or boating (i.e. marinas) activities. A total of 874 individuals (c. 30 per population) from five ports and 21 marinas was collected and analysed using both mitochondrial cytochrome c oxidase subunit I gene (COI) and 10 nuclear microsatellite markers. The geographical location of multiple hotspot populations indicates that multiple invasions have occurred in Nova Scotia. A loss of genetic diversity from port to marina populations suggests a stronger influence of ships than recreational boats on primary coastal introductions. Population genetic similarity analysis reveals a dependence of marina populations on those that had been previously established in ports. Empirical data on marina connectivity because of boating better explains patterns in population similarities than does natural spread. We conclude that frequent primary introductions arise by ships and that secondary spread occurs gradually thereafter around individual ports, facilitated by recreational boating.     

Between cheap and costly signals: the evolution of partially honest communication

Costly signalling theory has become a common explanation for honest communication when interests conflict. In this paper, we provide an alternative explanation for partially honest communication that does not require significant signal costs. We show that this alternative is at least as plausible as traditional costly signalling, and we suggest a number of experiments that might be used to distinguish the two theories.     

Dicoumarol derivatives: Green synthesis and molecular modelling studies of their anti-LOX activity

Dicoumarol derivatives were synthesized in the InCl₃ catalyzed pseudo three-component reactions of 4-hydroxycoumarin with aromatic aldehydes in excellent yields. The reactions were performed in water under microwave irradiation. All synthesized compounds were characterized using NMR, IR, and UV–Vis spectroscopy, as well as with TD-DFT. Obtained dicoumarols were subjected to evaluation of their in vitro lipid peroxidation and soybean lipoxygenase inhibition activities. It was shown that five of ten examined compounds (3e, 3h, 3b, 3d, 3f) possess significant potential of antilipid peroxidation (84–97%), and that compounds 3b, 3e, 3h provided the highest soybean lipoxygenase (LOX-Ib) inhibition (IC₅₀ = 52.5 µM) and 3i somewhat lower activity (IC₅₀ = 55.5 µM). The bioactive conformations of the best LOX-Ib inhibitors were obtained by means of molecular docking and molecular dynamics. It was shown that, within the bioactive conformations interior to LOX-Ib active site, the most active compounds form the pyramidal structure made of two 4-hydroxycoumarin cores and a central phenyl substituent. This form serves as a spatial barrier which prevents LOX-Ib Fe²⁺/Fe³⁺ ion activity to generate the coordinative bond with the C13 hydroxyl group of the α-linoleate. It is worth pointing out that the most active compounds 3b, 3e, 3h and 3i can be candidates for further examination of their in vitro and in vivo anti-inflammatory activity and that molecular modeling study results provide possibility to screen bioactive conformations and elucidate the mechanism of dicoumarols anti-LOX activity.     

Dissipative particle dynamics simulation of a gold nanoparticle system

Dissipative particle dynamics (DPD) was carried out to study systems containing gold atoms, organic ether (oligohydroquinonyl ether terminated with a thiol group) and organic solvents. The components in the simulated system are very different in size and chemical nature. Our simulation showed that the reproduction of the macroscopic experimental phase separation, properly dividing the polymeric molecule into beads, selecting the size of gold bead, and choosing the appropriate interaction parameters between beads are crucial. In addition, the solvent effect was the dominant factor for the formation of spherical aggregates of Au atoms and organic ether molecules. We report the interaction strengths between the solvent and gold clusters. Our work has demonstrated that DPD methods can be applied to the study of complex meso-scale systems.