Population and community-level compositional patterns shape the realized niche of the rare arctic-alpine species Carex lachenalii Schkuhr

The ability for vegetative growth and development of generative organs often reflects an adaptation to the environment and may be a suitable proxy for understanding population dynamics of rare relict species. An example of such a plant is Carex lachenalii Schkuhr, an arctic-alpine species, in the temperate zone of Europe only occurring in isolated localities of high-elevation mountain ranges. We aimed to assess whether there were relationships between flower production and clonal growth of C. lachenalii, both at the tuft and plot level, and how co-occurring vegetation could modify this relationship. In the study we focused on population-level traits of C. lachenalii, vegetation traits and components of functional diversity. At the tuft level we found that the proportion of flowering ramets of C. lachenalii decreased with increasing diameter of the tuft. At the plot level, in snowbed vegetation C. lachenalii produced more flowering ramets. We suggest this is due to higher environmental stress, expressed by high importance of habitat filtering (low functional dispersion) in shaping species composition of co-occurring vegetation. In granite grasslands and milder environment (expressed by higher functional dispersion), C. lachenalii produced more vegetative ramets, which we suggest is a result of a more competitive environment. While in snowbeds investment in flowering ramets could promote successful persistence of C. lachenalii, survival of subpopulations occurring in the highly competitive conditions of granite grasslands may be uncertain due to potentially weak adaptation to competition with graminoids and dwarf shrubs.     

Potential evolutionary trade-off between feeding and stability in Cambrian cinctan echinoderms

Reconstructing the function and behaviour of extinct groups of echinoderms is problematic because there are no modern analogues for their aberrant body plans. Cinctans, an enigmatic group of Cambrian echinoderms, exemplify this problem: their asymmetrical body plan differentiates them from all living species. Here, we used computational fluid dynamics to analyse the functional performance of cinctans without assuming an extant comparative model. Three-dimensional models of six species from across cinctan phylogeny were used in computer simulations of water flow. The results demonstrate that cinctans with strongly flattened bodies produced much less drag than species characterized by dorsal protuberances or swellings, suggesting the former were more stable on the seafloor. However, unlike the flattened forms, cinctans with high-relief bodies were able to passively direct flow towards the mouth and associated food grooves, indicating that they were capable of more efficient feeding on particles suspended in the water. This study provides evidence of a previously unknown evolutionary trade-off between feeding and stability in Cambrian cinctan echinoderms.     

Biotic and abiotic factors driving the diversification dynamics of Crocodylia

Species diversity patterns are governed by complex interactions among biotic and abiotic factors over time and space, but are essentially the result of the diversification dynamics (differential speciation and extinction rates) over the long-term evolutionary history of a clade. Previous studies have suggested that temporal variation in global temperature drove long-term diversity changes in Crocodylia, a monophyletic group of large ectothermic organisms. We use a large database of crocodylian fossil occurrences (192 spp.) and body mass estimations, under a taxic approach, to characterize the global diversification dynamics of crocodylians since the Cretaceous, and their correlation with multiple biotic and abiotic factors in a Bayesian framework. The diversification dynamic of crocodylians, which appears to have originated in the Turonian (c. 92.5 Ma), is characterized by several phases with high extinction and speciation rates within a predominantly low long-term mean rate. Our results reveal long-term diversification dynamics of Crocodylia to be a highly complex process driven by a combination of biotic and abiotic factors which influenced the speciation and extinction rates in dissimilar ways. Higher crocodylian extinction rates are related to low body mass disparity, indicating selective extinctions of taxa at both ends of the body mass spectrum. Speciation rate slowdowns are noted when the diversity of the clade is high and the warm temperate climatic belt is reduced. Our finding supports the idea that temporal variations of body mass disparity, self-diversity, and the warm climate belt size provided more direct mechanistic explanations for crocodylian diversification than do proxies of global temperature.     

Binding mechanism of spinosine and venenatine molecules with p300 HAT enzyme: Molecular screening,molecular dynamics and free-energy analysis

The chromatin modification is regulated by the histone acetyltransferase (HAT) and histone deacetyltransferase (HDAC) enzymes; abnormal function of these enzymes leads to several malignant diseases. The inhibition of these enzymes using natural ligand molecules is an emerging technique to cure these diseases. The in vitro analysis of natural molecules, venenatine, spinosine, palmatine and taxodione are giving the best inhibition rate against p300 HAT enzyme. However, the detailed understanding of binding and the stability of these molecules with p300 HAT is not yet known. The aim of the present study is focused to determine the binding strength of the molecules from molecular dynamics simulation analysis. The docking analysis confirms that, the venenatine (−6.97 kcal/mol - conformer 8), spinosine (−6.52 kcal/mol conformer −10), palmatine (−5.72 kcal/mol conformer-3) and taxodione (−4.99 kcal/mol conformer-4) molecules form strong hydrogen bonding interactions with the key amino acid residues (Arg1410, Thr1411 and Trp1466) present in the active site of p300. In the molecular dynamics (MD) simulation, the spinosine retain these key interactions with the active site amino acid residues (Arg1410, Thr1411, and Trp1466) than venenatine and are stable throughout the simulation. The RMSD value of spinosine (0.5 to 1.3 Å) and venenatine (0.3 to 1.3 Å) are almost equal during the MD simulation. However, during the MD simulation, the intermolecular interaction between venenatine and the active site amino acid residues (Arg1410, Thr1411, and Trp1466) decreased on comparing with the spinosine-p300 interaction. The binding free energy of the spinosine (−15.30 kcal/mol) is relatively higher than the venenatine (−11.8 kcal/mol); this increment is attributed to the strong hydrogen bonding interactions of spinosine molecule with the active site amino acid residues of p300.     

The mechanism of high affinity pentasaccharide binding to antithrombin,insights from Gaussian accelerated molecular dynamics simulations

Abstract

The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from ‘static’ structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a ‘non-activated’ AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. ‘Generalized correlation’ matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.

Communicated by Ramaswamy H. Sarma     

Inhibitory effect of α-ketoglutaric acid on α-glucosidase: integrating molecular dynamics simulation and inhibition kinetics

Abstract

The inhibition of α-glucosidase is used as a key clinical approach to treat type 2 diabetes mellitus and thus, we assessed the inhibitory effect of α-ketoglutaric acid (AKG) on α-glucosidase with both an enzyme kinetic assay and computational simulations. AKG bound to the active site and interacted with several key residues, including ASP68, PHE157, PHE177, PHE311, ARG312, TYR313, ASN412, ILE434 and ARG439, as detected by protein–ligand docking and molecular dynamics simulations. Subsequently, we confirmed the action of AKG on α-glucosidase as mixed-type inhibition with reversible and rapid binding. The relevant kinetic parameter IC₅₀ was measured (IC₅₀ = 1.738?±?0.041?mM), and the dissociation constant was determined (K_{i Slope} = 0.46?±?0.04?mM). Regarding the relationship between structure and activity, a high AKG concentration induced the slight modulation of the shape of the active site, as monitored by hydrophobic exposure. This tertiary conformational change was linked to AKG inhibition and mostly involved regional changes in the active site. Our study provides insight into the functional role of AKG due to its structural property of a hydroxyphenyl ring that interacts with the active site. We suggest that similar hydroxyphenyl ring-containing compounds targeting key residues in the active site might be potential α-glucosidase inhibitors. Abbreviations AKG alpha-ketoglutaric acid

pNPG 4-nitrophenyl-α-d-glucopyranoside

ANS 1-anilinonaphthalene-8-sulfonate

MD molecular dynamics.

Communicated by Ramaswamy H. Sarma     

The investigation of the G-quadruplex aptamer selectivity to Pb2+ ion: a joint molecular dynamics simulation and density functional theory study

Abstract

The aptamers with the ability to form a G-quadruplex structure can be stable in the presence of some ions. Hence, study of the interactions between such aptamers and ions can be beneficial to determine the highest selective aptamer toward an ion. In this article, molecular dynamics (MD) simulations and quantum mechanics (QM) calculations have been applied to investigate the selectivity of the T30695 aptamer toward Pb²⁺ in comparison with some ions. The Free Energy Landscape (FEL) analysis indicates that Pb²⁺ has remained inside the aptamer during the MD simulation, while the other ions have left it. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) binding energies prove that the conformational stability of the aptamer is the highest in the presence of Pb²⁺. According to the compaction parameters, the greatest compressed ion-aptamer complex, and hence, the highest ion-aptamer interaction have been induced in the presence of Pb²⁺. The contact maps clarify the closer contacts between the nucleotides of the aptamer in the presence of Pb²⁺. The density functional theory (DFT) results show that Pb²⁺ forms the most stable complex with the aptamer, which is consistent with the MD results. The QM calculations reveal that the N-H bonds and the O…H distances are the longest and the shortest, respectively, in the presence of Pb²⁺. The obtained results verify that the strongest hydrogen bonds (HBs), and hence, the most compressed aptamer structure are induced by Pb²⁺. Besides, atoms in molecules (AIM) and natural bond orbital (NBO) analyses confirm the results.

Communicated by Ramaswamy H. Sarma     

Simulated and experimental estimates of hydrodynamic drag from bio-logging tags

Drag force acting on swimming marine mammals is difficult to measure directly. Researchers often use simple modeling and kinematic measurements from animals, or computational fluid dynamics (CFD) simulations to estimate drag. However, studies that compare these methods are lacking. Here, computational simulation and physical experiments were used to estimate drag forces on gliding bottlenose dolphins (Tursiops truncatus). To facilitate comparison, variable drag loading (no-tag, tag, tag + 4, tag + 8) was used to increase force in both simulations and experiments. During the experiments, two dolphins were trained to perform controlled glides with variable loading. CFD simulations of dolphin/tag geometry in steady flow (1–6 m/s) were used to model drag forces. We expect both techniques will capture relative changes created by experimental conditions, but absolute forces predicted by the methods will differ. CFD estimates were within a calculated 90% confidence interval of the experimental results for all but the tag condition. Relative drag increase predicted by the simulation vs. experiment, respectively, differed by between 21% and 31%: tag, 4% vs. 33%; tag + 4, 47% vs. 68%; and tag + 8, 108% vs. 77%. The results from this work provide a direct comparison of computational and experimental estimates of drag, and provide a framework to quantify uncertainty.