RosettaDDGPrediction for high-throughput mutational scans: From stability to binding

Reliable prediction of free energy changes upon amino acid substitutions (ΔΔGs) is crucial to investigate their impact on protein stability and protein–protein interaction. Advances in experimental mutational scans allow high-throughput studies thanks to multiplex techniques. On the other hand, genomics initiatives provide a large amount of data on disease-related variants that can benefit from analyses with structure-based methods. Therefore, the computational field should keep the same pace and provide new tools for fast and accurate high-throughput ΔΔG calculations. In this context, the Rosetta modeling suite implements effective approaches to predict folding/unfolding ΔΔGs in a protein monomer upon amino acid substitutions and calculate the changes in binding free energy in protein complexes. However, their application can be challenging to users without extensive experience with Rosetta. Furthermore, Rosetta protocols for ΔΔG prediction are designed considering one variant at a time, making the setup of high-throughput screenings cumbersome. For these reasons, we devised RosettaDDGPrediction, a customizable Python wrapper designed to run free energy calculations on a set of amino acid substitutions using Rosetta protocols with little intervention from the user. Moreover, RosettaDDGPrediction assists with checking completed runs and aggregates raw data for multiple variants, as well as generates publication-ready graphics. We showed the potential of the tool in four case studies, including variants of uncertain significance in childhood cancer, proteins with known experimental unfolding ΔΔGs values, interactions between target proteins and disordered motifs, and phosphomimetics. RosettaDDGPrediction is available, free of charge and under GNU General Public License v3.0, at https://github.com/ELELAB/RosettaDDGPrediction .     

Molecular analysis of a microaerobically induced operon required for hydrogenase synthesis in Rhizobium leguminosarum biovar viciae

The nucleotide sequence (6138 bp) of a microaerobically inducible region (hupV/VI) from the Rhizobium leguminosarum bv. viciae hydrogenase gene cluster has been determined. Six genes, arranged as a single operon, were identified, and designated hypA, B, F, C, D and E based on the sequence similarities of all of them, except hypF, to genes from the hydrogenase pleiotropic operon (hyp) from Escherichia coli. The gene products from hypBFCDE were identified by in vivo expression analysis in E. coli, and their molecular sizes were consistent with those predicted from the nucleotide sequence. Transposon Tn5 insertions into hypB, hypF, hypD and hypE resulted in R. leguminosarum mutants that lacked any hydrogenase activity in symbiosis with peas, but still were able to synthesize the polypeptide for the hydrogenase large subunit. The gene products HypA, HypB, HypF and HypD contained CX₂C motifs characteristic of metal-binding proteins. In addition, HypB bore a long histidine-rich stretch of amino acids near the N-terminus, suggesting a possible role in nickel binding for this protein. The gene product HypF, which was translationally coupled to HypB, presented two cysteine motifs (CX₂CX₈₁CX₂C) with a capacity to form zinc finger-like structures in the N-terminal third of the protein. A role in nickel metabolism in relation to hydrogenase synthesis is postulated for proteins HypB and HypF.     

Novel Miscanthus hybrids: Modelling productivity on marginal land in Europe using dynamics of canopy development determined by light interception

New biomass crop hybrids for bioeconomic expansion require yield projections to determine their potential for strategic land use planning in the face of global challenges. Our biomass growth simulation incorporates radiation interception and conversion efficiency. Models often use leaf area to predict interception which is demanding to determine accurately, so instead we use low-cost rapid light interception measurements using a simple laboratory-made line ceptometer and relate the dynamics of canopy closure to thermal time, and to measurements of biomass. We apply the model to project the European biomass potentials of new market-ready hybrids for 2020–2030. Field measurements are easier to collect, the calibration is seasonally dynamic and reduces influence of weather variation between field sites. The model obtained is conservative, being calibrated by crops of varying establishment and varying maturity on less productive (marginal) land. This results in conservative projections of miscanthus hybrids for 2020–2030 based on 10% land use conversion of the least (productive) grassland and arable for farm diversification, which show a European potential of 80.7–89.7 Mt year⁻¹ biomass, with potential for 1.2–1.3 EJ year⁻¹ energy and 36.3–40.3 Mt year⁻¹ carbon capture, with seeded Miscanthus sacchariflorus × sinensis displaying highest yield potential. Simulated biomass projections must be viewed in light of the field measurements on less productive land with high soil water deficits. We are attempting to model the results from an ambitious and novel project combining new hybrids across Europe with agronomy which has not been perfected on less productive sites. Nevertheless, at the time of energy sourcing issues, seed-propagated miscanthus hybrids for the upscaled provision of bioenergy offer an alternative source of renewable energy. If European countries provide incentives for growers to invest, seeded hybrids can improve product availability and biomass yields over the current commercial miscanthus variety.     

Poleward increase in feeding efficiency of leafminer Stigmella lapponica (Lepidoptera: Nepticulidae) in a latitudinal gradient crossing a boreal forest zone

Damage to plant communities imposed by insect herbivores generally decreases from low to high latitudes. This decrease is routinely attributed to declines in herbivore abundance and/or diversity, whereas latitudinal changes in per capita food consumption remain virtually unknown. Here, we tested the hypothesis that the lifetime food consumption by a herbivore individual decreases from low to high latitudes due to a temperature-driven decrease in metabolic expenses. From 2016 to 2019, we explored latitudinal changes in multiple characteristics of linear (gallery) mines made by larvae of the pygmy moth, Stigmella lapponica, in leaves of downy birch, Betula pubescens. The mined leaves were larger than intact leaves at the southern end of our latitudinal gradient (at 60°N) but smaller than intact leaves at its northern end (at 69°N), suggesting that female oviposition preference changes with latitude. No latitudinal changes were observed in larval size, mine length or area, and in per capita food consumption, but the larval feeding efficiency (quantified as the ratio between larval size and mine size) increased with latitude. Consequently, S. lapponica larvae consumed less foliar biomass at higher latitudes than at lower latitudes to reach the same size. Based on space-for-time substitution, we suggest that climate warming will increase metabolic expenses of insect herbivores with uncertain consequences for plant–herbivore interactions.     

Low genetic diversity and population connectivity fuel vulnerability to climate change for the Tertiary relict pine Pinus bungeana

Endemic species are important components of regional biodiversity and hold the key to understanding local adaptation and evolutionary processes that shape species distributions. This study investigated the biogeographic history of a relict conifer Pinus bungeana Zucc. ex Endl. confined to central China. We examined genetic diversity in P. bungeana using genotyping-by-sequencing and chloroplast and mitochondrial DNA markers. We performed spatial and temporal inference of recent genetic and demographic changes, and dissected the impacts of geography and environmental gradients on population differentiation. We then projected P. bungeana's risk of decline under future climates. We found extremely low nucleotide diversity (average π 0.0014), and strong population structure (global F_ST 0.234) even at regional scales, reflecting long-term isolation in small populations. The species experienced severe bottlenecks in the early Pliocene and continued to decline in the Pleistocene in the western distribution, whereas the east expanded recently. Local adaptation played a small (8%) but significant role in population diversity. Low genetic diversity in fragmented populations makes the species highly vulnerable to climate change, particularly in marginal and relict populations. We suggest that conservation efforts should focus on enhancing gene pool and population growth through assisted migration within each genetic cluster to reduce the risk of further genetic drift and extinction.     

Systematics of Mukdenia and Oresitrophe (Saxifragaceae): Insights from genome skimming data

Oresitrophe and Mukdenia (Saxifragaceae) are epilithic sister genera used in traditional Chinese medicine. The taxonomy of Mukdenia, especially of M. acanthifolia, has been controversial. To address this, we produced plastid and mitochondrial data using genome skimming for Mukdenia acanthifolia and Mukdenia rossii, including three individuals of each species. We assembled complete plastomes, mitochondrial CDS and nuclear ribosomal ETS/ITS sequences using these data. Comparative analysis shows that the plastomes of Mukdenia and Oresitrophe are relatively conservative in terms of genome size, structure, gene content, RNA editing sites and codon usage. Five plastid regions that represent hotspots of change (trnH-psbA, psbC-trnS, trnM-atpE, petA-psbJ and ccsA-ndhD) are identified within Mukdenia, and six regions (trnH-psbA, petN-psbM, trnM-atpE, rps16-trnQ, ycf1 and ndhF) contain a higher number of species-specific parsimony-informative sites that may serve as potential DNA barcodes for species identification. To infer phylogenetic relationships between Mukdenia and Oresitrophe, we combined our data with published data based on three different datasets. The monophyly of each species (Oresitrophe rupifraga, M. acanthifolia and M. rossii) and the inferred topology ((M. rossii, M. acanthifolia), O. rupifraga) are well supported in trees reconstructed using the complete plastome sequences, but M. acanthifolia and M. rossii did not form a separate clade in the trees based on ETS + ITS data, while the mitochondrial CDS trees are not well-resolved. We found low recovery of genes in the Angiosperms353 target enrichment panel from our unenriched genome skimming data. Hybridization or incomplete lineage sorting may be the cause of discordance between trees reconstructed from organellar and nuclear data. Considering its morphological distinctiveness and our molecular phylogenetic results, we strongly recommend that M. acanthifolia be treated as a distinct species.     

First chemical synthesis of a scorpion alpha-toxin affecting sodium channels: the Aah I toxin of Androctonus australis hector.

Aah I is a 63-residue alpha-toxin isolated from the venom of the Buthidae scorpion Androctonus australis hector, which is considered to be the most dangerous species. We report here the first chemical synthesis of Aah I by the solid-phase method, using a Fmoc strategy. The synthetic toxin I (sAah I) was renatured in DMSO-Tris buffer, purified and subjected to thorough analysis and comparison with the natural toxin. The sAah I showed physico-chemical (CD spectrum, molecular mass, HPLC elution), biochemical (amino-acid composition, sequence), immunochemical and pharmacological properties similar to those of the natural toxin. The synthetic toxin was recognized by a conformation-dependent monoclonal anti-Aah I antibody, with an IC50 value close to that for the natural toxin. Following intracerebroventricular injection, the synthetic and the natural toxins were similarly lethal to mice. In voltage-clamp experiments, Na(v) 1.2 sodium channel inactivation was inhibited by the application of sAah I or of the natural toxin in a similar way. This work describes a simple protocol for the chemical synthesis of a scorpion alpha-toxin, making it possible to produce structural analogues in time.     

Studies on the catalytic mechanism of H2-forming methylenetetrahydromethanopterin dehydrogenase: para-ortho H2 conversion rates in H2O and D2O

H₂–forming N ⁵,N ¹⁰ –methylenetetrahydromethanopterin dehydrogenase is a novel type of hydrogenase that contains neither nickel nor iron-sulfur clusters. Evidence has been presented that the reaction mechanism catalyzed by the enzyme is very similar to that of the formation of carbocations and H₂ from alkanes under superacidic conditions. We present here further results in support of this mechanism. It was found that the purified enzyme per se did not catalyze the conversion of para H₂ to ortho H₂, a reaction catalyzed by all other hydrogenases known to date. However, it catalyzed the conversion in the presence of the substrate N ⁵,N ¹⁰ –methenyltetrahydromethanopterin (CH≡H₄MPT⁺), indicating that for heterolytic cleavage of H₂ the enzyme-CH≡H₄MPT⁺ complex is required. In D₂O, the formation of HD and D₂ from H₂ rather than a para–ortho H₂ conversion was observed, indicating that after heterolytic cleavage of H₂ the dissociation of the proton from the enzyme-substrate complex is fast relative to the re-formation of free H₂.