Plant protein-lipid interfaces studied by molecular dynamics simulations

The delineation of protein-lipid interfaces is essential for understanding the mechanisms of various membrane-associated processes crucial to plant development and growth, including signalling, trafficking, and membrane transport. Due to their highly dynamic nature, the precise characterization of lipid-protein interactions is challenging by experimental techniques. Molecular dynamics (MD) simulations provide a powerful computational alternative with a spatial-temporal resolution allowing the atomistic-level description. In this review, we aim to introduce plant scientists to the MD simulations. We describe different steps of performing the MD simulations and provide a broad survey of the MD studies investigating plant protein-lipid interfaces. Our aim is also to illustrate that combining the MD simulations with artificial intelligence-based protein structure determination opens unprecedented possibilities for future investigations of dynamic plant protein-lipid interfaces.

Misfolding of fukutin-related protein (FKRP) variants in congenital and limb girdle muscular dystrophies

Fukutin-related protein (FKRP, MIM ID 606596) variants cause a range of muscular dystrophies associated with hypo-glycosylation of the matrix receptor, α-dystroglycan. These disorders are almost exclusively caused by homozygous or compound heterozygous missense variants in the FKRP gene that encodes a ribitol phosphotransferase. To understand how seemingly diverse FKRP missense mutations may contribute to disease, we examined the synthesis, intracellular dynamics, and structural consequences of a panel of missense mutations that encompass the disease spectrum. Under non-reducing electrophoresis conditions, wild type FKRP appears to be monomeric whereas disease-causing FKRP mutants migrate as high molecular weight, disulfide-bonded aggregates. These results were recapitulated using cysteine-scanning mutagenesis suggesting that abnormal disulfide bonding may perturb FKRP folding. Using fluorescence recovery after photobleaching, we found that the intracellular mobility of most FKRP mutants in ATP-depleted cells is dramatically reduced but can, in most cases, be rescued with reducing agents. Mass spectrometry showed that wild type and mutant FKRP differentially associate with several endoplasmic reticulum (ER)-resident chaperones. Finally, structural modelling revealed that disease-associated FKRP missense variants affected the local environment of the protein in small but significant ways. These data demonstrate that protein misfolding contributes to the molecular pathophysiology of FKRP-deficient muscular dystrophies and suggest that molecules that rescue this folding defect could be used to treat these disorders.

Modelling 3D supramolecular structure from sparse single-molecule localisation microscopy data

Single-molecule localisation microscopy (SMLM) has the potential to reveal the underlying organisation of specific molecules within supramolecular complexes and their conformations, which is not possible with conventional microscope resolution. However, the detection efficiency for fluorescent molecules in cells can be limited in SMLM, even to below 1% in thick and dense samples. Segmentation of individual complexes can also be challenging. To overcome these problems, we have developed a software package termed PERPL: Pattern Extraction from Relative Positions of Localisations. This software assesses the relative likelihoods of models for underlying patterns behind incomplete SMLM data, based on the relative positions of pairs of localisations. We review its principles and demonstrate its use on the 3D lattice of Z-disk proteins in mammalian cardiomyocytes. We find known and novel features at ~20 nm with localisations of less than 1% of the target proteins, using mEos fluorescent protein constructs.

Evolution of anelloviruses from a circovirus-like ancestor through gradual augmentation of the jelly-roll capsid protein

Anelloviruses are highly prevalent in diverse mammals, including humans, but so far have not been linked to any disease and are considered to be part of the 'healthy virome'. These viruses have small circular single-stranded DNA (ssDNA) genomes and encode several proteins with no detectable sequence similarity to proteins of other known viruses. Thus, anelloviruses are the only family of eukaryotic ssDNA viruses currently not included in the realm Monodnaviria. To gain insights into the provenance of these enigmatic viruses, we sequenced more than 250 complete genomes of anelloviruses from nasal and vaginal swab samples of Weddell seal (Leptonychotes weddellii) from Antarctica and a fecal sample of grizzly bear (Ursus arctos horribilis) from the USA and performed a comprehensive family-wide analysis of the signature anellovirus protein ORF1. Using state-of-the-art remote sequence similarity detection approaches and structural modeling with AlphaFold2, we show that ORF1 orthologs from all Anelloviridae genera adopt a jelly-roll fold typical of viral capsid proteins (CPs), establishing an evolutionary link to other eukaryotic ssDNA viruses, specifically, circoviruses. However, unlike CPs of other ssDNA viruses, ORF1 encoded by anelloviruses from different genera display remarkable variation in size, due to insertions into the jelly-roll domain. In particular, the insertion between β-strands H and I forms a projection domain predicted to face away from the capsid surface and function at the interface of virus-host interactions. Consistent with this prediction and supported by recent experimental evidence, the outermost region of the projection domain is a mutational hotspot, where rapid evolution was likely precipitated by the host immune system. Collectively, our findings further expand the known diversity of anelloviruses and explain how anellovirus ORF1 proteins likely diverged from canonical jelly-roll CPs through gradual augmentation of the projection domain. We suggest assigning Anelloviridae to a new phylum, 'Commensaviricota', and including it into the kingdom Shotokuvirae (realm Monodnaviria), alongside Cressdnaviricota and Cossaviricota.

Isolation, purification and biochemical characterization of Conopomorpha cramerella farnesol dehydrogenase

In Southeast Asia, Conopomorpha cramerella (Snellen) which is commonly known as the cocoa pod borer (CPB) moth has been identified as the most detrimental pest of Theobroma cacao L. Apart from the various side effects on human health and non-target organisms, heavily relying on synthetic pyrethroid insecticides to control CPB infestations also increases the environmental contamination risks. Thus, developing biorational insecticides that minimally affect the non-target organism and environment by targeting the insect growth regulation process is needed to manage the pest population. In insects, juvenile hormones (JH) regulate critical biological events, especially metamorphosis, development and reproduction. Since the physiological roles of JH III vary among different organisms, the biochemical properties, especially substrate specificity and analogue inhibition, may also be different. Therefore, studies on the JH III biosynthetic pathway enzymes in both plants and insects are beneficial to discover more effective analogues. Bioinformatic analysis and biochemical characterization of a NADP⁺ -dependent farnesol dehydrogenase, an intermediate enzyme of the JH III pathway, from C. cramerella (CcFolDH), were described in this study. In addition, the farnesol analogues that may act as a potent analogue inhibitor for CcFolDH ware determined using in vitro enzymatic study. The phylogenetic analysis indicated that CcFolDH shared a close phylogenetic relationship to the honeybee's short-chain dehydrogenase/reductase. The 27 kDa CcFolDH has an NADP(H) binding domain with a typical Rossmann fold and is likely a homotetrameric protein in the solution. The enzyme had a greater preference for substrate trans, trans-farnesol and coenzyme NADP⁺ . In terms of analogue inhibitor inhibition, hexahydroxyfarnesyl acetone showed the highest inhibition (the lowest K_i ) compared to other farnesol analogues. Thus, hexahydroxyfarnesyl acetone would serve as the most potent active ingredient for future biorational pesticide management for C. cramerella infestation. Based on the bioinformatic analyses and biochemical characterizations conducted in this research, we proposed that rCcFolDH differs slightly from other reported farnesol dehydrogenases in terms of molecular weight, substrate preference, coenzymes utilization and analogue inhibitors selection.

Anastrozole-mediated modulation of mitochondrial activity by inhibition of mitochondrial permeability transition pore opening: an initial perspective

The mitochondrial permeability transition pore (mtPTP) plays a vital role in altering the structure and function of mitochondria. Cyclophilin D (CypD) is a mitochondrial protein that regulates mtPTP function and a known drug target for therapeutic studies involving mitochondria. While the effect of aromatase inhibition on the mtPTP has been studied previously, the effect of anastrozole on the mtPTP has not been completely elucidated. The role of anastrozole in modulating the mtPTP was evaluated by docking, molecular dynamics and network-guided studies using human CypD data. The peripheral blood mononuclear cells (PBMCs) of patients with mitochondrial disorders and healthy controls were treated with anastrozole and evaluated for mitochondrial permeability transition pore (mtPTP) function and apoptosis using a flow cytometer. Spectrophotometry was employed for estimating total ATP levels. The anastrozole-CypD complex is more stable than cyclosporin A (CsA)-CypD. Anastrozole performed better than cyclosporine in inhibiting mtPTP. Additional effects included inducing mitochondrial membrane depolarization and a reduction in mitochondrial swelling and superoxide generation, intrinsic caspase-3 activity and cellular apoptosis, along with an increase in ATP levels. Anastrozole may serve as a potential therapeutic agent for mitochondrial disorders and ameliorate the clinical phenotype by regulating the activity of mtPTP. However, further studies are required to substantiate our preliminary findings.Communicated by Ramaswamy H. Sarma.

Application of Transactional (Cross-lagged panel) Models in Mental Health Research: An Introduction and Review of Methodological Considerations

Transactional models employing cross-lagged panels have been used for over 40 years to examine the longitudinal relations and directional associations between variables of interest to child and adolescent mental health. Through a narrative synthesis of the literature, we provide an accessible overview of cross-lagged panels with attention to developing a research question, study design and assumptions, dynamic effects (including the random-intercept cross-lagged panel model), and reporting and interpretation of results. Implications and critical appraisal guidelines for readers are discussed throughout. Overall, several key points are highlighted, with particular emphasis on the intended level of inference, model and measure selection, and timing of assessments. Despite limitations in establishing causation, cross-lagged panel models are fundamental to non-experimental epidemiologic research in child mental health and development.

Structural modelling and in silico pharmacology of β-carboline alkaloids as potent 5-HT1A receptor antagonists and reuptake inhibitors

Serotonin (5-HT) antagonists and reuptake inhibitors (SARIs) are atypical antidepressants for managing major depressive disorder. They are oftentimes applied as adjuvants for ameliorating aftereffects of SSRI antidepressants including insomnia and sexual dysfunction. The few available candidates of this class including lorpiprazole and trazodone also display some daunting side effects, making a continuous search for improved alternatives essential. Natural β-carboline alkaloids (NβCs) are interestingly renowned with broad pharmacological spectrum against several neuropsychiatric disorders including depression. However, their potentials as SARIs remain underexplored. In this study, 982 NβCs retrieved from the Ambinter-Greenpharma (Amb) database were virtually screened for potent SARI alternatives using computational and biocheminformatics approaches: homology modelling of 5-HT1A receptor, Glide HTVS, SP and XP molecular docking, molecular dynamics (MD) simulation, ADMET and mutagenicity predictions. The homology receptor was validated as a good representative of human 5HT1A receptor using the RCSB structure validation and quality protocols. From the virtual screening against the 5-HT1A receptor, Amb ligands, Amb18709727 and Amb37857532 showed higher binding affinities by XP scores of -8.725 and -7.976 kcal/mol, and MMGBSA of -87.972 and -107.585 kcal/mol respectively compared to lorpiprazole, a reference SARI with XP score and MMGBSA of -6.512 and -62.788 kcal/mol respectively. They maintained ideal contacts with pharmacologically essential amino acid residues implicated in SARI mechanisms and expressed higher stability and compactness than lorpiprazole throughout the trajectories of 100 ns MD simulation. They also displayed interesting ADME, druggability, low toxicity and mutagenicity profiles, ideal for CNS drug prospects, thus, recommended as putative SARI candidates for further study.

Phase transitions and (p-T-X) behaviour of centrosymmetric perovskites: modelling with transformed crystallographic data

A reversible transformation of the unit-cell parameters and atomic coordinates of centrosymmetric perovskites ABX₃ into a Cartesian space is defined. Analytical expressions for the three vectors for the pseudocubic cell and three vectors for a BX₆ octahedron are derived for space groups Pbmn, Cmcm, Ibmm, P4/mbm, P4/nmc, I4/mcm and R3c. The following structural parameters may be derived from these vectors: up to six pseudocubic parameters defining octahedral geometry; length- and angle-based octahedral distortion parameters λ and σ; inclination angles of tilted octahedra, θ₁, θ₂ and θ₃; angles of tilt of octahedra; AX₁₂:BX₆ polyhedral volume ratio, V_A/V_B; parameters η_A and η_B defining the relative contraction of inner AX₈ polyhedra and expansion of BX₆ octahedra due to octahedral tilting. The application of these parameters is demonstrated by reference to published crystal structures. The variation of η_A and η_B with temperature in the compositional series Sr_xBa_1-xSnO₃ and Sr_xBa_1-xHfO₃, as well as the temperature series of BaPbO₃ and CaTiO₃, is related to the sequence of phases Pbmn → Ibmm→ Pm3m. Stabilization of the Cmcm phase is likewise interpreted in terms of these two parameters for NaTaO₃ and NaNbO₃. The pressure evolution of the structures of MgSiO₃, YAlO₃, (La_1-xNd_x)GaO₃ (0 ≤ x ≤ 1) and YAl_0.25Cr_0.75O₃ is modelled with the appropriate structural parameters, thereby also addressing the characteristics of the Pbmn → R3c transition. Simulation of MgSiO₃ up to 125 GPa and of YAlO₃ up to 52 GPa in space group Pbnm is carried out by using the Birch-Murnaghan equation of state. In both cases, full sets of oxygen coordinates assuming regular octahedra are generated. Octahedral distortion is also modelled in the latter system and predicted to have a key influence on structural evolution and the sequence of phase transitions. The core modelling procedures are made available as a Microsoft Excel file.

Topological analysis of a bacterial DedA protein associated with alkaline tolerance and antimicrobial resistance

Maintaining membrane integrity is of paramount importance to the survival of bacteria as the membrane is the site of multiple crucial cellular processes including energy generation, nutrient uptake and antimicrobial efflux. The DedA family of integral membrane proteins are widespread in bacteria and are associated with maintaining the integrity of the membrane. In addition, DedA proteins have been linked to resistance to multiple classes of antimicrobials in various microorganisms. Therefore, the DedA family are attractive targets for the development of new antibiotics. Despite DedA family members playing a key physiological role in many bacteria, their structure, function and physiological role remain unclear. To help illuminate the structure of the bacterial DedA proteins, we performed substituted cysteine accessibility method (SCAM) analysis on the most comprehensively characterized bacterial DedA protein, YqjA from Escherichia coli. By probing the accessibility of 15 cysteine residues across the length of YqjA using thiol reactive reagents, we mapped the topology of the protein. Using these data, we experimentally validated a structural model of YqjA generated using evolutionary covariance, which consists of an α-helical bundle with two re-entrant hairpin loops reminiscent of several secondary active transporters. In addition, our cysteine accessibility data suggest that YqjA forms an oligomer wherein the protomers are arranged in a parallel fashion. This experimentally verified model of YqjA lays the foundation for future work in understanding the function and mechanism of this interesting and important family.

Predicting the zoonotic capacity of mammals to transmit SARS-CoV-2

Back and forth transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) between humans and animals will establish wild reservoirs of virus that endanger long-term efforts to control COVID-19 in people and to protect vulnerable animal populations. Better targeting surveillance and laboratory experiments to validate zoonotic potential requires predicting high-risk host species. A major bottleneck to this effort is the few species with available sequences for angiotensin-converting enzyme 2 receptor, a key receptor required for viral cell entry. We overcome this bottleneck by combining species' ecological and biological traits with three-dimensional modelling of host-virus protein-protein interactions using machine learning. This approach enables predictions about the zoonotic capacity of SARS-CoV-2 for greater than 5000 mammals-an order of magnitude more species than previously possible. Our predictions are strongly corroborated by in vivo studies. The predicted zoonotic capacity and proximity to humans suggest enhanced transmission risk from several common mammals, and priority areas of geographic overlap between these species and global COVID-19 hotspots. With molecular data available for only a small fraction of potential animal hosts, linking data across biological scales offers a conceptual advance that may expand our predictive modelling capacity for zoonotic viruses with similarly unknown host ranges.

In Silico Analysis of the Enzymes Involved in Haloarchaeal Denitrification

During the last century, anthropogenic activities such as fertilization have led to an increase in pollution in many ecosystems by nitrogen compounds. Consequently, researchers aim to reduce nitrogen pollutants following different strategies. Some haloarchaea, owing to their denitrifier metabolism, have been proposed as good model organisms for the removal of not only nitrate, nitrite, and ammonium, but also (per)chlorates and bromate in brines and saline wastewater. Bacterial denitrification has been extensively described at the physiological, biochemical, and genetic levels. However, their haloarchaea counterparts remain poorly described. In previous work the model structure of nitric oxide reductase was analysed. In this study, a bioinformatic analysis of the sequences and the structural models of the nitrate, nitrite and nitrous oxide reductases has been described for the first time in the haloarchaeon model Haloferax mediterranei. The main residues involved in the catalytic mechanism and in the coordination of the metal centres have been explored to shed light on their structural characterization and classification. These results set the basis for understanding the molecular mechanism for haloarchaeal denitrification, necessary for the use and optimization of these microorganisms in bioremediation of saline environments among other potential applications including bioremediation of industrial waters.

ATSAS 3.0: expanded functionality and new tools for small-angle scattering data analysis

The ATSAS software suite encompasses a number of programs for the processing, visualization, analysis and modelling of small-angle scattering data, with a focus on the data measured from biological macromolecules. Here, new developments in the ATSAS 3.0 package are described. They include IMSIM, for simulating isotropic 2D scattering patterns; IMOP, to perform operations on 2D images and masks; DATRESAMPLE, a method for variance estimation of structural invariants through parametric resampling; DATFT, which computes the pair distance distribution function by a direct Fourier transform of the scattering data; PDDFFIT, to compute the scattering data from a pair distance distribution function, allowing comparison with the experimental data; a new module in DATMW for Bayesian consensus-based concentration-independent molecular weight estimation; DATMIF, an ab initio shape analysis method that optimizes the search model directly against the scattering data; DAMEMB, an application to set up the initial search volume for multiphase modelling of membrane proteins; ELLLIP, to perform quasi-atomistic modelling of liposomes with elliptical shapes; NMATOR, which models conformational changes in nucleic acid structures through normal mode analysis in torsion angle space; DAMMIX, which reconstructs the shape of an unknown intermediate in an evolving system; and LIPMIX and BILMIX, for modelling multilamellar and asymmetric lipid vesicles, respectively. In addition, technical updates were deployed to facilitate maintainability of the package, which include porting the PRIMUS graphical interface to Qt5, updating SASpy - a PyMOL plugin to run a subset of ATSAS tools - to be both Python 2 and 3 compatible, and adding utilities to facilitate mmCIF compatibility in future ATSAS releases. All these features are implemented in ATSAS 3.0, freely available for academic users at https://www.embl-hamburg.de/biosaxs/software.html.

Sources of Artifacts in SLODR Detection

Background

Spearman's law of diminishing returns (SLODR) states that intercorrelations between scores on tests of intellectual abilities were higher when the data set was comprised of subjects with lower intellectual abilities and vice versa. After almost a hundred years of research, this trend has only been detected on average.

Objective

To determine whether the very different results were obtained due to variations in scaling and the selection of subjects.

Design

We used three methods for SLODR detection based on moderated factor analysis (MFCA) to test real data and three sets of simulated data. Of the latter group, the first one simulated a real SLODR effect. The second one simulated the case of a different density of tasks of varying difficulty; it did not have a real SLODR effect. The third one simulated a skewed selection of respondents with different abilities and also did not have a real SLODR effect. We selected the simulation parameters so that the correlation matrix of the simulated data was similar to the matrix created from the real data, and all distributions had similar skewness parameters (about -0.3).

Results

The results of MFCA are contradictory and we cannot clearly distinguish by this method the dataset with real SLODR from datasets with similar correlation structure and skewness, but without a real SLODR effect. The results allow us to conclude that when effects like SLODR are very subtle and can be identified only with a large sample, then features of the psychometric scale become very important, because small variations of scale metrics may lead either to masking of real SLODR or to false identification of SLODR.

Bovine seminal plasma osteopontin: Structural modelling, recombinant expression and its relationship with semen quality

Osteopontin (OPN) is a multifunctional phosphoprotein that has been linked to fertility in bulls. However, the exact mechanism by which OPN contributes to fertilisation is yet unknown. The biotechnological use of OPN in bovine reproduction is promising but some gaps remain unfilled. The present work aimed: (a) to verify whether the seminal plasma OPN is associated with seminal traits and a standard breeding soundness exam; (b) to predict OPN interactions with integrins, CD44 and glycosaminoglycans through molecular docking; and (c) to develop a protocol for recombinant expression of OPN from vesicular gland cDNA. Ejaculates from top ranked bulls had higher amounts of seminal plasma OPN in comparison with bulls classified as questionable (p < .01). The structural modelling and molecular docking predictions indicated that bovine OPN binds to heparin disaccharide, hyaluronic acid and hyaluronan. In addition, docking studies described the binding complexes of OPN with CD44 and the integrin heterodimers α5β1 and αVβ3. Finally, expression of rOPN-6His was successfully obtained after 3 hr of induction with 0.5 mM IPTG at 37°C and a denaturing purification protocol resulted in efficiently purified recombinant OPN. The present results contribute to the development of biotechnological uses of OPN as a biomarker in bovine reproduction.

Structure-function analysis of ZAR1 immune receptor reveals key molecular interactions for activity

NLR (nucleotide-binding [NB] leucine-rich repeat [LRR] receptor) proteins are critical for inducing immune responses in response to pathogen proteins, and must be tightly modulated to prevent spurious activation in the absence of a pathogen. The ZAR1 NLR recognizes diverse effector proteins from Pseudomonas syringae, including HopZ1a, and Xanthomonas species. Receptor-like cytoplasmic kinases (RLCKs) such as ZED1, interact with ZAR1 and provide specificity for different effector proteins, such as HopZ1a. We previously developed a transient expression system in Nicotiana benthamiana that allowed us to demonstrate that ZAR1 function is conserved from the Brassicaceae to the Solanaceae. Here, we combined structural modelling of ZAR1, with molecular and functional assays in our transient system, to show that multiple intramolecular and intermolecular interactions modulate ZAR1 activity. We identified determinants required for the formation of the ZARCC oligomer and its activity. Lastly, we characterized intramolecular interactions between ZAR1 subdomains that participate in keeping ZAR1 immune complexes inactive. This work identifies molecular constraints on immune receptor function and activation.

A bidirectional path analysis model of smoking cessation self-efficacy and concurrent smoking status: impact on abstinence outcomes

Self-efficacy is routinely associated with abstinence in the addictions literature, and is a major component relapse-prevention models. The magnitude of this relationship has been brought into question following equivocal results in studies controlling for concurrent smoking status. The aim of our study was to clarify the relationship between cessation self-efficacy, smoking status, and cessation outcomes in a cohort of treatment-seeking smokers. Smokers participating in the FLEX trial, a randomized trial investigating the efficacy of three pharmacologic treatments for smoking cessation, completed questionnaires assessing cessation self-efficacy at baseline and at weeks 1, 3, 5 and 10 post-target quit date; smoking status was verified using expired carbon monoxide. Structural models were fit in order to ascertain the relationship between cessation self-efficacy and concurrent smoking at each time-point, and to assess the association between cessation self-efficacy, smoking and seven-day point prevalence smoking status at week 10. A total of 737 treatment-seeking smokers participated. In our path model, self-efficacy and smoking status at all time points were associated with week 10 abstinence (except week 3 self-efficacy), after controlling these values' previous time-points. All direct pathways between cessation self-efficacy and smoking were also significant, supporting a bidirectional relationship. Our results support a bidirectional and reciprocal relationship between cessation self-efficacy and concurrent smoking behavior; participants with higher confidence were more likely to be smoke-free, and concurrent smoking status predicted levels of confidence over the ensuing weeks. Both measures were associated with week 10 abstinence. Our results indicate that while correlated, both cessation self-efficacy and current smoking behavior during a cessation attempt are important independent markers of ultimate cessation success.

CD9 and CD81 Interactions and Their Structural Modelling in Sperm Prior to Fertilization

Proteins CD9 and CD81 are members of the tetraspanin superfamily and were detected in mammalian sperm, where they are suspected to form an active tetraspanin web and to participate in sperm–egg membrane fusion. The importance of these two proteins during the early stages of fertilization is supported by the complete sterility of CD9/CD81 double null female mice. In this study, the putative mechanism of CD9/CD81 involvement in tetraspanin web formation in sperm and its activity prior to fertilization was addressed. Confocal microscopy and colocalization assay was used to determine a mutual CD9/CD81 localization visualised in detail by super-resolution microscopy, and their interaction was address by co-immunoprecipitation. The species-specific traits in CD9 and CD81 distribution during sperm maturation were compared between mice and humans. A mutual position of CD9/CD81 is shown in human spermatozoa in the acrosomal cap, however in mice, CD9 and CD81 occupy a distinct area. During the acrosome reaction in human sperm, only CD9 is relocated, compared to the relocation of both proteins in mice. The structural modelling of CD9 and CD81 homologous and possibly heterologous network formation was used to propose their lateral Cis as well as Trans interactions within the sperm membrane and during sperm–egg membrane fusion.

Linking Compositional and Functional Predictions to Decipher the Biogeochemical Significance in DFAA Turnover of Abundant Bacterioplankton Lineages in the North Sea

Deciphering the ecological traits of abundant marine bacteria is a major challenge in marine microbial ecology. In the current study, we linked compositional and functional predictions to elucidate such traits for abundant bacterioplankton lineages in the North Sea. For this purpose, we investigated entire and active bacterioplankton composition along a transect ranging from the German Bight to the northern North Sea by pyrotag sequencing of bacterial 16S rRNA genes and transcripts. Functional profiles were inferred from 16S rRNA data using Tax4Fun. Bacterioplankton communities were dominated by well-known marine lineages including clusters/genera that are affiliated with the Roseobacter group and the Flavobacteria. Variations in community composition and function were significantly explained by measured environmental and microbial properties. Turnover of dissolved free amino acids (DFAA) showed the strongest correlation to community composition and function. We applied multinomial models, which enabled us to identify bacterial lineages involved in DFAA turnover. For instance, the genus Planktomarina was more abundant at higher DFAA turnover rates, suggesting its vital role in amino acid degradation. Functional predictions further indicated that Planktomarina is involved in leucine and isoleucine degradation. Overall, our results provide novel insights into the biogeochemical significance of abundant bacterioplankton lineages in the North Sea.

2017 publication guidelines for structural modelling of small-angle scattering data from biomolecules in solution: an update

In 2012, preliminary guidelines were published addressing sample quality, data acquisition and reduction, presentation of scattering data and validation, and modelling for biomolecular small-angle scattering (SAS) experiments. Biomolecular SAS has since continued to grow and authors have increasingly adopted the preliminary guidelines. In parallel, integrative/hybrid determination of biomolecular structures is a rapidly growing field that is expanding the scope of structural biology. For SAS to contribute maximally to this field, it is essential to ensure open access to the information required for evaluation of the quality of SAS samples and data, as well as the validity of SAS-based structural models. To this end, the preliminary guidelines for data presentation in a publication are reviewed and updated, and the deposition of data and associated models in a public archive is recommended. These guidelines and recommendations have been prepared in consultation with the members of the International Union of Crystallography (IUCr) Small-Angle Scattering and Journals Commissions, the Worldwide Protein Data Bank (wwPDB) Small-Angle Scattering Validation Task Force and additional experts in the field.

Progress in small-angle scattering from biological solutions at high-brilliance synchrotrons

Small-angle X-ray scattering (SAXS) is an established technique that provides low-resolution structural information on macromolecular solutions. Recent decades have witnessed significant progress in both experimental facilities and in novel data-analysis approaches, making SAXS a mainstream method for structural biology. The technique is routinely applied to directly reconstruct low-resolution shapes of proteins and to generate atomistic models of macromolecular assemblies using hybrid approaches. Very importantly, SAXS is capable of yielding structural information on systems with size and conformational polydispersity, including highly flexible objects. In addition, utilizing high-flux synchrotron facilities, time-resolved SAXS allows analysis of kinetic processes over time ranges from microseconds to hours. Dedicated bioSAXS beamlines now offer fully automated data-collection and analysis pipelines, where analysis and modelling is conducted on the fly. This enables SAXS to be employed as a high-throughput method to rapidly screen various sample conditions and additives. The growing SAXS user community is supported by developments in data and model archiving and quality criteria. This review illustrates the latest developments in SAXS, in particular highlighting time-resolved applications aimed at flexible and evolving systems.

Seismic and Restoration Assessment of Monumental Masonry Structures

Masonry structures are complex systems that require detailed knowledge and information regarding their response under seismic excitations. Appropriate modelling of a masonry structure is a prerequisite for a reliable earthquake-resistant design and/or assessment. However, modelling a real structure with a robust quantitative (mathematical) representation is a very difficult, complex and computationally-demanding task. The paper herein presents a new stochastic computational framework for earthquake-resistant design of masonry structural systems. The proposed framework is based on the probabilistic behavior of crucial parameters, such as material strength and seismic characteristics, and utilizes fragility analysis based on different failure criteria for the masonry material. The application of the proposed methodology is illustrated in the case of a historical and monumental masonry structure, namely the assessment of the seismic vulnerability of the Kaisariani Monastery, a byzantine church that was built in Athens, Greece, at the end of the 11th to the beginning of the 12th century. Useful conclusions are drawn regarding the effectiveness of the intervention techniques used for the reduction of the vulnerability of the case-study structure, by means of comparison of the results obtained.

ATSAS 2.8: a comprehensive data analysis suite for small-angle scattering from macromolecular solutions

Developments and improvements of the ATSAS software suite (versions 2.5–2.8) for analysis of small-angle scattering data of biological macromolecules or nanoparticles are described.

ATSAS is a comprehensive software suite for the analysis of small-angle scattering data from dilute solutions of biological macromolecules or nanoparticles. It contains applications for primary data processing and assessment, ab initio bead modelling, and model validation, as well as methods for the analysis of flexibility and mixtures. In addition, approaches are supported that utilize information from X-ray crystallography, nuclear magnetic resonance spectroscopy or atomistic homology modelling to construct hybrid models based on the scattering data. This article summarizes the progress made during the 2.5–2.8 ATSAS release series and highlights the latest developments. These include AMBIMETER, an assessment of the reconstruction ambiguity of experimental data; DATCLASS, a multiclass shape classification based on experimental data; SASRES, for estimating the resolution of ab initio model reconstructions; CHROMIXS, a convenient interface to analyse in-line size exclusion chromatography data; SHANUM, to evaluate the useful angular range in measured data; SREFLEX, to refine available high-resolution models using normal mode analysis; SUPALM for a rapid superposition of low- and high-resolution models; and SASPy, the ATSAS plugin for interactive modelling in PyMOL. All these features and other improvements are included in the ATSAS release 2.8, freely available for academic users from https://www.embl-hamburg.de/biosaxs/software.html.

Structural modelling and phylogenetic analyses of PgeIF4A2 (Eukaryotic translation initiation factor) from Pennisetum glaucum reveal signature motifs with a role in stress tolerance and development

Eukaryotic translation initiation factor 4A (eIF4A) is an indispensable component of the translation machinery and also play a role in developmental processes and stress alleviation in plants and animals. Different eIF4A isoforms are present in the cytosol of the cell, namely, eIF4A1, eIF4A2, and eIF4A3 and their expression is tightly regulated in cap-dependent translation. We revealed the structural model of PgeIF4A2 protein using the crystal structure of Homo sapiens eIF4A3 (PDB ID: 2J0S) as template by Modeller 9.12. The resultant PgeIF4A2 model structure was refined by PROCHECK, ProSA, Verify3D and RMSD that showed the model structure is reliable with 77 % amino acid sequence identity with template. Investigation revealed two conserved signatures for ATP-dependent RNA Helicase DEAD-box conserved site (VLDEADEML) and RNA helicase DEAD-box type, Q-motif in sheet-turn-helix and α-helical region respectively. All these conserved motifs are responsible for response during developmental stages and stress tolerance in plants.

Partial Least Squares with Structured Output for Modelling the Metabolomics Data Obtained from Complex Experimental Designs: A Study into the Y-Block Coding

Partial least squares (PLS) is one of the most commonly used supervised modelling approaches for analysing multivariate metabolomics data. PLS is typically employed as either a regression model (PLS-R) or a classification model (PLS-DA). However, in metabolomics studies it is common to investigate multiple, potentially interacting, factors simultaneously following a specific experimental design. Such data often cannot be considered as a “pure” regression or a classification problem. Nevertheless, these data have often still been treated as a regression or classification problem and this could lead to ambiguous results. In this study, we investigated the feasibility of designing a hybrid target matrix Y that better reflects the experimental design than simple regression or binary class membership coding commonly used in PLS modelling. The new design of Y coding was based on the same principle used by structural modelling in machine learning techniques. Two real metabolomics datasets were used as examples to illustrate how the new Y coding can improve the interpretability of the PLS model compared to classic regression/classification coding.

Redefining the structural motifs that determine RNA binding and RNA editing by pentatricopeptide repeat proteins in land plants

The pentatricopeptide repeat (PPR) proteins form one of the largest protein families in land plants. They are characterised by tandem 30-40 amino acid motifs that form an extended binding surface capable of sequence-specific recognition of RNA strands. Almost all of them are post-translationally targeted to plastids and mitochondria, where they play important roles in post-transcriptional processes including splicing, RNA editing and the initiation of translation. A code describing how PPR proteins recognise their RNA targets promises to accelerate research on these proteins, but making use of this code requires accurate definition and annotation of all of the various nucleotide-binding motifs in each protein. We have used a structural modelling approach to define 10 different variants of the PPR motif found in plant proteins, in addition to the putative deaminase motif that is found at the C-terminus of many RNA-editing factors. We show that the super-helical RNA-binding surface of RNA-editing factors is potentially longer than previously recognised. We used the redefined motifs to develop accurate and consistent annotations of PPR sequences from 109 genomes. We report a high error rate in PPR gene models in many public plant proteomes, due to gene fusions and insertions of spurious introns. These consistently annotated datasets across a wide range of species are valuable resources for future comparative genomics studies, and an essential pre-requisite for accurate large-scale computational predictions of PPR targets. We have created a web portal (http://www.plantppr.com) that provides open access to these resources for the community.

Structural insight into selectivity and resistance profiles of ROS1 tyrosine kinase inhibitors

Oncogenic ROS1 fusion proteins are molecular drivers in multiple malignancies, including a subset of non-small cell lung cancer (NSCLC). The phylogenetic proximity of the ROS1 and anaplastic lymphoma kinase (ALK) catalytic domains led to the clinical repurposing of the Food and Drug Administration (FDA)-approved ALK inhibitor crizotinib as a ROS1 inhibitor. Despite the antitumor activity of crizotinib observed in both ROS1- and ALK-rearranged NSCLC patients, resistance due to acquisition of ROS1 or ALK kinase domain mutations has been observed clinically, spurring the development of second-generation inhibitors. Here, we profile the sensitivity and selectivity of seven ROS1 and/or ALK inhibitors at various levels of clinical development. In contrast to crizotinib's dual ROS1/ALK activity, cabozantinib (XL-184) and its structural analog foretinib (XL-880) demonstrate a striking selectivity for ROS1 over ALK. Molecular dynamics simulation studies reveal structural features that distinguish the ROS1 and ALK kinase domains and contribute to differences in binding site and kinase selectivity of the inhibitors tested. Cell-based resistance profiling studies demonstrate that the ROS1-selective inhibitors retain efficacy against the recently reported CD74-ROS1(G2032R) mutant whereas the dual ROS1/ALK inhibitors are ineffective. Taken together, inhibitor profiling and stringent characterization of the structure-function differences between the ROS1 and ALK kinase domains will facilitate future rational drug design for ROS1- and ALK-driven NSCLC and other malignancies.

Ambiguity assessment of small-angle scattering curves from monodisperse systems

A novel approach is presented for an a priori assessment of the ambiguity associated with spherically averaged single-particle scattering. The approach is of broad interest to the structural biology community, allowing the rapid and model-independent assessment of the inherent non-uniqueness of three-dimensional shape reconstruction from scattering experiments on solutions of biological macromolecules. One-dimensional scattering curves recorded from monodisperse systems are nowadays routinely utilized to generate low-resolution particle shapes, but the potential ambiguity of such reconstructions remains a major issue. At present, the (non)uniqueness can only be assessed by a posteriori comparison and averaging of repetitive Monte Carlo-based shape-determination runs. The new a priori ambiguity measure is based on the number of distinct shape categories compatible with a given data set. For this purpose, a comprehensive library of over 14,000 shape topologies has been generated containing up to seven beads closely packed on a hexagonal grid. The computed scattering curves rescaled to keep only the shape topology rather than the overall size information provide a `scattering map' of this set of shapes. For a given scattering data set, one rapidly obtains the number of neighbours in the map and the associated shape topologies such that in addition to providing a quantitative ambiguity measure the algorithm may also serve as an alternative shape-analysis tool. The approach has been validated in model calculations on geometrical bodies and its usefulness is further demonstrated on a number of experimental X-ray scattering data sets from proteins in solution. A quantitative ambiguity score (a-score) is introduced to provide immediate and convenient guidance to the user on the uniqueness of the ab initio shape reconstruction from the given data set.

New developments in the ATSAS program package for small-angle scattering data analysis

New developments in the program package ATSAS (version 2.4) for the processing and analysis of isotropic small-angle X-ray and neutron scattering data are described. They include (i) multiplatform data manipulation and display tools, (ii) programs for automated data processing and calculation of overall parameters, (iii) improved usage of high- and low-resolution models from other structural methods, (iv) new algorithms to build three-dimensional models from weakly interacting oligomeric systems and complexes, and (v) enhanced tools to analyse data from mixtures and flexible systems. The new ATSAS release includes installers for current major platforms (Windows, Linux and Mac OSX) and provides improved indexed user documentation. The web-related developments, including a user discussion forum and a widened online access to run ATSAS programs, are also presented.

Biochemical characterization and structural modeling of human cathepsin E variant 2 in comparison to the wild-type protein

Cathepsin E splice variant 2 appears in a number of gastric carcinomas. Here we report detecting this variant in HeLa cells using polyclonal antibodies and biotinylated inhibitor pepstatin A. An overexpression of GFP fusion proteins of cathepsin E and its splice variant within HEK-293T cells was performed to show their localization. Their distribution under a fluorescence microscope showed that they are colocalized. We also expressed variants 1 and 2 of cathepsins E, with propeptide and without it, in Escherichia coli. After refolding from the inclusion bodies, the enzymatic activity and circular dichroism spectra of the splice variant 2 were compared to those of the wild-type mature active cathepsins E. While full-length cathepsin E variant 1 is activated at acid pH, the splice variant remains inactive. In contrast to the active cathepsin E, the splice variant 2 predominantly assumes β-sheet structure, prone to oligomerization, at least under in vitro conditions, as shown by atomic force microscopy as shallow disk-like particles. A comparative structure model of splice variant 2 was computed based on its alignment to the known structure of cathepsin E intermediate (Protein Data Bank code 1TZS) and used to rationalize its conformational properties and loss of activity.

The d subunit plays a central role in human vacuolar H(+)-ATPases

The multi-subunit vacuolar-type H(+)-ATPase consists of a V(1) domain (A-H subunits) catalyzing ATP hydrolysis and a V(0) domain (a, c, c', c", d, e) responsible for H(+) translocation. The mammalian V(0) d subunit is one of the least-well characterized, and its function and position within the pump are still unclear. It has two different forms encoded by separate genes, d1 being ubiquitous while d2 is predominantly expressed at the cell surface in kidney and osteoclast. To determine whether it forms part of the pump's central stalk as suggested by bacterial A-ATPase studies, or is peripheral as hypothesized from a yeast model, we investigated both human d subunit isoforms. In silico structural modelling demonstrated that human d1 and d2 are structural orthologues of bacterial subunit C, despite poor sequence identity. Expression studies of d1 and d2 showed that each can pull down the central stalk's D and F subunits from human kidney membrane, and in vitro studies using D and F further showed that the interactions between these proteins and the d subunit is direct. These data indicate that the d subunit in man is centrally located within the pump and is thus important in its rotary mechanism.