Generic reclassification and species boundaries in the rediscovered freshwater mussel ‘<Emphasis Type="Italic">Quadrula’ mitchelli</Emphasis> (Simpson in Dall, 1896)

The Central Texas endemic freshwater mussel, Quadrula mitchelli (Simpson in Dall, 1896), had been presumed extinct until relict populations were recently rediscovered. To help guide ongoing and future conservation efforts focused on Q. mitchelli we set out to resolve several uncertainties regarding its evolutionary history, specifically its unknown generic position and untested species boundaries. We designed a molecular matrix consisting of two loci (cytochrome c oxidase subunit I and internal transcribed spacer I) and 57 terminal taxa to test the generic position of Q. mitchelli using Bayesian inference and maximum likelihood phylogenetic reconstruction. We also employed two Bayesian species validation methods to test five a priori species models (i.e. hypotheses of species delimitation). Our study is the first to test the generic position of Q. mitchelli and we found robust support for its inclusion in the genus Fusconaia. Accordingly, we introduce the binomial, Fusconaia mitchelli comb. nov., to accurately represent the systematic position of the species. We resolved F. mitchelli individuals in two well supported and divergent clades that were generally distinguished as distinct species using Bayesian species validation methods, although alternative hypotheses of species delineation were also supported. Despite strong evidence of genetic isolation within F. mitchelli, we do not advocate for species-level status of the two clades as they are allopatrically distributed and no morphological, behavioral, or ecological characters are known to distinguish them. These results are discussed in the context of the systematics, distribution, and conservation of F. mitchelli.     

Landscape genetics of a tropical rescue pollinator

Pollination services are increasingly threatened by the loss and modification of natural habitats, posing a risk to the maintenance of both native plant biodiversity and agricultural production. In order to safeguard pollination services, it is essential to examine the impacts of habitat degradation on the population dynamics of key pollinators and identify potential “rescue pollinators” capable of persisting in these human-altered landscapes. Using a landscape genetic approach, we assessed the impact of landscape structure on genetic differentiation in the widely-distributed tropical stingless bee Trigona spinipes (Apidae: Meliponini) across agricultural landscape mosaics composed of coffee plantations and Atlantic forest fragments in southeastern Brazil. We genotyped 115 bees at 16 specific and highly polymorphic microsatellite loci, developed using next-generation sequencing. Our results reveal that T. spinipes is capable of dispersing across remarkably long distances, as we did not find genetic differentiation across a 200 km range, nor fine-scale spatial genetic structure. Furthermore, gene flow was not affected by forest cover, land cover, or elevation, indicating that reproductive individuals are able to disperse well through agricultural landscapes and across altitudinal gradients. We also found evidence of a recent population expansion, suggesting that this opportunistic stingless bee is capable of colonizing degraded habitats. Our results thus suggest that T. spinipes can persist in heavily-altered landscapes and can be regarded as a rescue pollinator, potentially compensating for the decline of other native pollinators in degraded tropical landscapes.     

Exploring ligand dissociation pathways from aminopeptidase N using random acceleration molecular dynamics simulation

Aminopeptidase N (APN) is a zinc-dependent ectopeptidase involved in cell proliferation, secretion, invasion, and angiogenesis, and is widely recognized as an important cancer target. However, the mechanisms whereby ligands leave the active site of APN remain unknown. Investigating ligand dissociation processes is quite difficult, both in classical simulation methods and in experimental approaches. In this study, random acceleration molecular dynamics (RAMD) simulation was used to investigate the potential dissociation pathways of ligand from APN. The results revealed three pathways (channels A, B and C) for ligand release. Channel A, which matches the hypothetical channel region, was the most preferred region for bestatin to dissociate from the enzyme, and is probably the major channel for the inner bound ligand. In addition, two alternative channels (channels B and C) were shown to be possible pathways for ligand egression. Meanwhile, we identified key residues controlling the dynamic features of APN channels. Identification of the dissociation routes will provide further mechanistic insights into APN, which will benefit the development of more promising APN inhibitors.

Open image in new window

Graphical Abstract The release pathways of bestatin inside active site of aminopeptidase N were simulated using RAMD simulation

     

Protonation and electronic structure of 2,6-dichlorophenolindophenolate during reduction. A theoretical study including explicit solvent

Protonation in the two-electron/two-proton reduction processes of 2,6-dichlorophenolindophenolate (DCIP) is investigated combining density functional theory (DFT) and molecular dynamics (MD) methods. DCIP (anion), DCIP^?– (radical anion), and DCIP^2? (dianion) are considered, including the electronic structure analysis from the prospective of quantum theory of atoms and molecules (QTAIM). It is shown that oxygen on the indophenolate moiety and nitrogen are the first and/or the second proton acceptor sites and their energetic order depends on the total charge of the system. MD simulations of differently charged species interacting with the solvent molecules have been performed for methanol, water, and oxonium cation (H₃O⁺). Methanol and water molecules are found to form only hydrogen bonds with the solute irrespective of its charge. The calculated pK_a values show that the imino group of DCIPH^? is a weaker acid than water. While in the case of DCIP (and DCIP^?–) plus oxonium cation, proton transfer from the solvent to the solute was evidenced for both aforementioned acceptor sites. In addition, MD simulations of bulks containing 15 and 43 molecules of water around the DCIP molecule have been performed, revealing the formation of 2–4 hydrogen bonds.

Open image in new window

Graphical Abstract 2,6-Dichlorophenolindophenolate interacts with solvent molecules (water, oxonium cation and methanol). Hydrogen transfer and electronic structure are studied by DFT and molecular dynamics methods

     

Modulating the strength of tetrel bonding through beryllium bonding

Quantum chemical calculations were performed to investigate the stability of the ternary complexes BeH₂···XMH₃···NH₃ (X?=?F, Cl, and Br; M?=?C, Si, and Ge) and the corresponding binary complexes at the atomic level. Our results reveal that the stability of the XMH₃···BeH₂ complexes is mainly due to both a strong beryllium bond and a weak tetrel–hydride interaction, while the XMH₃···NH₃ complexes are stabilized by a tetrel bond. The beryllium bond with a halogen atom as the electron donor has many features in common with a beryllium bond with an O or N atom as the electron donor, although they do exhibit some different characteristics. The stability of the XMH₃···NH₃ complex is dominated by the electrostatic interaction, while the orbital interaction also makes an important contribution. Interestingly, as the identities of the X and M atoms are varied, the strength of the tetrel bond fluctuates in an irregular manner, which can explained by changes in electrostatic potentials and orbital interactions. In the ternary systems, both the beryllium bond and the tetrel bond are enhanced, which is mainly ascribed to increased electrostatic potentials on the corresponding atoms and charge transfer. In particular, when compared to the strengths of the tetrel and beryllium bonds in the binary systems, in the ternary systems the tetrel bond is enhanced to a greater degree than the beryllium bond.

Open image in new window

Graphical Abstract A tetrel bond can be strengthened greatly by a beryllium bond

     

De novo design of VEGFR-2 tyrosine kinase inhibitors based on a linked-fragment approach

Vascular endothelial growth factor receptor-2 (VEGFR-2) tyrosine kinase inhibitors have been demonstrated to possess substantial antitumor activity. VEGFR-2 tyrosine kinase inhibitors are crucial for development of antitumor drugs. Based on the crystal structure of VEGFR-2 tyrosine kinase, a linked-fragment strategy was employed to design novel VEGFR-2 tyrosine kinase inhibitors, and 1000 compounds were generated in this process. Absorption, distribution, metabolism, excretion and toxicity (ADMET) were used to screen the 1000 compounds, and 59 compounds were acceptable. Scaffold hopping was then used for further screening, and only four compounds were obtained in this way. Then, the binding energy of the four molecules to VEGFR-2 tyrosine kinase was calculated using molecular docking, and their values were found to be lower than that of Sorafenib. Finally, molecular dynamics simulations were performed on the complex of the compound with the lowest binding energy with VEGFR-2 tyrosine kinase, and the binding model was analyzed. At the end, four chemical entities with novel structures were obtained, and were suggested for experimental testing in future studies.     

Insight into drug resistance mechanisms and discovery of potential inhibitors against wild-type and L1196M mutant ALK from FDA-approved drugs

Anaplastic lymphoma kinase (ALK) plays a crucial role in multiple malignant cancers. It is known as a well-established target for the treatment of ALK-dependent cancers. Even though substantial efforts have been made to develop ALK inhibitors, only crizotinib, ceritinib, and alectinib had been approved by the U.S. Food and Drug Administration for patients with ALK-positive non-small cell lung cancer (NSCLC). The secondary mutations with drug-resistance bring up difficulties to develop effective drugs for ALK-positive cancers. To give a comprehensive understanding of molecular mechanism underlying inhibitor response to ALK tyrosine kinase mutations, we established an accurate assessment for the extensive profile of drug against ALK mutations by means of computational approaches. The molecular mechanics-generalized Born surface area (MM-GBSA) method based on molecular dynamics (MD) simulation was carried out to calculate relative binding free energies for receptor-drug systems. In addition, the structure-based virtual screening was utilized to screen effective inhibitors targeting wild-type ALK and the gatekeeper mutation L1196M from 3180 approved drugs. Finally, the mechanism of drug resistance was discussed, several novel potential wild-type and L1196M mutant ALK inhibitors were successfully identified.     

The protonated 2-halogenated imidazolium cation as the noncovalent interaction donor: the σ-hole and π-hole interactions

The σ-hole and π-hole of the protonated 2-halogenated imidazolium cation (XC₃H₄N₂ ⁺; X = F, Cl, Br, I) were investigated and analyzed. The monomers of (CH₃)₃SiY(Y=F, Cl, Br, I), considered as the Lewis base, were combined with the σ-hole and π-hole of XC₃H₄N₂ ⁺ to form the σ-hole and π-hole interactions in the bimolecular complexes (CH₃)₃SiY?·?·?·?XC₃H₄N₂ ⁺ and (CH₃)₃SiY?·?·?·?C₃(X)H₄N₂ ⁺(X/Y=F, Cl, Br, I), respectively. For both the σ-hole and π-hole interactions, the equilibrium geometries of complexes show regular changes according to the sequence of heavy sequence of the noncovalent interaction acceptors and donors. The electrostatic energy is the main contribution in the formation of both kinds of interactions, it has linear relations with the V _S,max values of σ-hole and the V′ _S,max values of π-hole. Both the σ-hole and π-hole interactions belong to the closed-shell and noncovalent interactions. The π-hole interactions are stronger than the σ-hole interactions. For the π-hole interactions, the contribution percents of the dispersion energies are somewhat greater than those of the σ-hole interactions, while it is contrary for the polarization energy.

Open image in new window

Graphical Abstract The protonated 2-halogenated imidazolium cation as the noncovalent interaction donor: the σ-hole and π-hole interactions?

     

pH dependence of ligand-induced human epidermal growth factor receptor activation investigated by molecular dynamics simulations

The activation of human epidermal growth factor receptor (hEGFR) involves a large conformational change in its soluble extracellular domains (sECD, residues 1–620), from a tethered to an extended conformation upon binding of ligands, such as EGF. It has been reported that this dynamic process is pH-dependent, that is, hEGFR can be activated by EGF at high pH to form an extended dimer but remains as an inactive monomer at low pH. In this paper, we perform all-atom molecular dynamics (MD) simulations starting from the tethered conformation of sECD:EGF complex, at pH 5.0 and 8.5, respectively. Simulation results indicate that sECD:EGF shows different dynamic properties between the two pHs, and the complex may have a higher tendency of activation at pH 8.5. Twenty residues, including 13 histidines, in sECD:EGF have different protonation states between the two pHs (calculated by the H++ server). The charge distribution at pH 8.5 is more favorable for forming an extended conformation toward the active state of sECD than that at pH 5.0. Our study may shed light on the mechanism of pH dependence of hEGFR activation.

Open image in new window

Graphical abstract pH dependence of ligand-induced human epidermal growth factor receptor activation

     

Computational investigation on MB<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript> (M = Li-Cs,Be-Ba,Sc-La and Ti; <Emphasis Type="Italic">n</Emphasis> = 28 and 38)

Differing from the weakly antiaromatic B₈₀ buckyball, the medium-sized C ₁–B₂₈ and D _2h –B₃₈, as well as their mono- to tetra-anions, are highly aromatic, as indicated by the negative nucleus-independent chemical shifts (NICSs) at their cage centers. The interior cavities and high aromaticity of the B₂₈ and B₃₈ cages render them very promising hosts to accommodate diverse metal atoms. Accordingly, we carried out systematic density functional theory (DFT) computations on the structures, stabilities and electronic properties of metalloborofullerenes MB _n (M?=?Li, Na, K, Rb, Cs, Be, Mg, Ca, Sr, Ba, Sc, Y, La and Ti; n?=?28 and 38). Among them, besides the recently reported M@B₃₈(M?=?Sc, Y and Ti) [Lu et al. (2015) Phys Chem Chem Phys 17:20897–20902], Ti@B₂₈ and M@B₃₈ (M?=?Ca and La) also favor endohedral structures with large binding energies, and are suggested promising targets for experimental applications. Note that Ti@B₂₈ is the first endohedral derivative based on the new B₂₈ fullerene, and La@B₃₈ features the largest metal size inside a B₃₈ cage thus far. These endohedral derivatives, as exemplified by Ca@B₃₈, may exhibit σ and π double aromaticity over the whole cage surface, indicating their considerable stability. In contrast, the other metals prefer to reside at the exterior cage surface, due mainly to the mismatch of their sizes with the boron cages, though the size match is not the only factor to determine their doping form. Furthermore, the infrared absorption spectra and ¹¹B nuclear magnetic resonance spectra of the three new M@B _n complexes were computed to assist future experimental characterization.

Open image in new window

Graphical Abstract Putting more metals into medium-sized boron cages!