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 .     

Biotechnical Characteristics of Root Systems of Typical Mediterranean Species

Quantifying patterns of fine root dynamics is crucial to the understanding of ecosystem structure and function, and in predicting how ecosystems respond to disturbance. Part of this understanding involves consideration of the carbon lost through root turnover. In the context of the rainfall pattern in the tropics, it was hypothesised that rainfall would strongly influence fine root biomass and longevity. A field study was conducted to determine root biomass, elemental composition and the influence of rainfall on longevity of fine roots in a tropical lowland evergreen rainforest at Danum Valley, Sabah, Malaysia. A combination of root coring, elemental analysis and rhizotron observation methods were used. Fine (less than 2 mm diameter) root biomass was relatively low (1700 kg ha −1) compared with previously described rainforest data. Standing root biomass was positively correlated with preceding rainfall, and the low fine root biomass in the dry season contained higher concentrations of N and lower concentrations of P and K than at other times. Observations on rhizotrons demonstrated that the decrease in fine root biomass in the dry season was a product of both a decrease in fine root length appearance and an increase in fine root length disappearance. Fitting an overall model to root survival time showed significant effects of rainfall preceding root disappearance, with the hazard of root disappearance decreasing by 8 for each 1 mm increase in the average daily (30 day) rainfall preceding root disappearance. While it is acknowledged that other factors have a part to play, this work demonstrates the importance of rainfall and soil moisture in influencing root biomass and root disappearance in this tropical rainforest.     

Direct inhibition of hexokinase activity by metformin at least partially impairs glucose metabolism and tumor growth in experimental breast cancer

Emerging evidence suggests that metformin, a widely used anti-diabetic drug, may be useful in the prevention and treatment of different cancers. In the present study, we demonstrate that metformin directly inhibits the enzymatic function of hexokinase (HK) I and II in a cell line of triple-negative breast cancer (MDA-MB-231). The inhibition is selective for these isoforms, as documented by experiments with purified HK I and II as well as with cell lysates. Measurements of ¹⁸F-fluoro-deoxyglycose uptake document that it is dose- and time-dependent and powerful enough to virtually abolish glucose consumption despite unchanged availability of membrane glucose transporters. The profound energetic imbalance activates phosphorylation and is subsequently followed by cell death. More importantly, the “in vivo” relevance of this effect is confirmed by studies of orthotopic xenografts of MDA-MB-231 cells in athymic (nu/nu) mice. Administration of high drug doses after tumor development caused an evident tumor necrosis in a time as short as 48 h. On the other hand, 1 mo metformin treatment markedly reduced cancer glucose consumption and growth. Taken together, our results strongly suggest that HK inhibition contributes to metformin therapeutic and preventive potential in breast cancer.     

A structural modeling approach for the understanding of initiation and elongation of ALS-linked superoxide dismutase fibrils

Familial amyotrophic lateral sclerosis caused by mutations in copper-zinc superoxide dismutase (SOD1) is characterized by the presence of SOD1-rich inclusions in spinal cords. It has been shown that a reduced intra-subunit disulfide bridge apo-SOD1 can rapidly initiate fibrillation forming an inter-subunits disulfide under mild, physiologically accessible conditions. Once initiated, elongation can proceed via recruitment of either apo or partially metallated disulfide-intact SOD1 and the presence of copper, but not zinc, ions inhibit fibrillation. We propose a structural model, refined through molecular dynamics simulations, that, taking into account these experimental findings, provides a molecular explanation for the initiation and the elongation of SOD1 fibrils in physiological conditions. The model indicates the occurrence of a new dimeric unit, prone to interact one with the other due to the presence of a wide hydrophobic surface and specific electrostatic interactions. The model has dimensions consistent with the SOD1 fibril size observed through electron microscopy and provides a structural basis for the understanding of SOD1 fibrillation.

Improvement in cyst recovery and molecular detection of Giardia duodenalis from stool samples

Molecular Biology Reports - Molecular detection of Giardia duodenalis by polymerase chain reaction (PCR) is difficult in faecal samples due to inhibitors that contaminate DNA preparations, or due...     

Molecular dynamics simulations show how the FMRP Ile304Asn mutation destabilizes the KH2 domain structure and affects its function

Mutations or deletions of FMRP, involved in the regulation of mRNA metabolism in brain, lead to the Fragile X syndrome (FXS), the most frequent form of inherited intellectual disability. A severe manifestation of the disease has been associated with the Ile304Asn mutation, located on the KH2 domain of the protein. Several hypotheses have been proposed to explain the possible molecular mechanism responsible for the drastic effect of this mutation in humans. Here, we performed a molecular dynamics simulation and show that the Ile304Asn mutation destabilizes the hydrophobic core producing a partial unfolding of two α-helices and a displacement of a third one. The affected regions show increased residue flexibility and motion. Molecular docking analysis revealed strongly reduced binding to a model single-stranded nucleic acid in agreement with known data that the two partially unfolded helices form the RNA-binding surface. The third helix, which we show here to be also affected, is involved in the PAK1 protein interaction. These two functional binding sites on the KH2 domain do not overlap spatially, and therefore, they can simultaneously bind their targets. Since the Ile304Asn mutation affects both binding sites, this may justify the severe clinical manifestation observed in the patient in which both mRNA metabolism activity and cytoskeleton remodeling would be affected.     

Seasonal and environmental influences on the calling behaviour of Eurasian Scops Owls

Capsule Spontaneous calls and replies to playback of Scops Owls were significantly more frequent during the breeding period, when paired birds defend territories. The probability of spontaneous calls varied with moon phase, with a peak occurring during nights with a full moon. In some cases, during the winter, Scops Owls responded to playback calls of Little Owls, thus suggesting possible intra-guild competition between these species.