Structural variability of C3larvin toxin. Intrinsic dynamics of the α/β fold of the C3-like group of mono-ADP-ribosyltransferase toxins

C3larvin toxin is a new member of the C3 class of the mono-ADP-ribosyltransferase toxin family. The C3 toxins are known to covalently modify small G-proteins, e.g. RhoA, impairing their function, and serving as virulence factors for an offending pathogen. A full-length X-ray structure of C3larvin (2.3 Å) revealed that the characteristic mixed α/β fold consists of a central β-core flanked by two helical regions. Topologically, the protein can be separated into N and C lobes, each formed by a β-sheet and an α-motif, and connected by exposed loops involved in the recognition, binding, and catalysis of the toxin/enzyme, i.e. the ADP-ribosylation turn–turn and phosphate–nicotinamide PN loops. Herein, we provide two new C3larvin X-ray structures and present a systematic study of the toxin dynamics by first analyzing the experimental variability of the X-ray data-set followed by contrasting those results with theoretical predictions based on Elastic Network Models (GNM and ANM). We identify residues that participate in the stability of the N-lobe, putative hinges at loop residues, and energy-favored deformation vectors compatible with conformational changes of the key loops and 3D-subdomains (N/C-lobes), among the X-ray structures. We analyze a larger ensemble of known C3bot1 conformations and conclude that the characteristic ‘crab-claw’ movement may be driven by the main intrinsic modes of motion. Finally, via computational simulations, we identify harmonic and anharmonic fluctuations that might define the C3larvin ‘native state.’ Implications for docking protocols are derived.     

Algal Population Growth Model Integrated with Toxicokinetics for Ecological Risk Assessment under Time-Varying Pesticide Exposure

An algal population growth model integrated with toxicokinetics was developed for assessing the effect of pesticides on population dynamics. This model is a simple one-compartment, first-order kinetic model in which toxicity (growth inhibition and mortality) depends on the intracellular effective concentration of a pesticide at a target site. The model's parameters were derived using an experimental study that investigated the effects of pretilachlor, bensulfuron-methyl, pentoxazone, and quinoclamine on the growth, mortality, and subsequent population recovery of the green alga Pseudokirchneriella subcapitata. Modeled and measured trajectories of algal population dynamics agreed well. The effect on population recovery was underestimated by the model that ignored the toxicokinetics. The four tested herbicides had a variety of toxicity characteristics and physicochemical properties, indicating the wide range of the model's applicability. Moreover, the developed model and the obtained model's parameters were extrapolated to predict long-term algal population dynamics under time-varying herbicide exposure. The calculated integral biomass lost compared with the control was considered a quantitative index of the population-level ecological risk. The model's prediction showed that the same exposure level (peak concentration is equivalent to EC50) indicated much different population-level effect depending on the herbicide.     

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.     

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.