ConKit: a python interface to contact predictions.

SUMMARY: Recent advances in protein residue contact prediction algorithms have led to the emergence of many new methods and a variety of file formats. We present ConKit , an open source, modular and extensible Python interface which allows facile conversion between formats and provides an interface to analyses of sequence alignments and sets of contact predictions. AVAILABILITY AND IMPLEMENTATION: ConKit is available via the Python Package Index. The documentation can be found at http://www.conkit.org . ConKit is licensed under the BSD 3-Clause. CONTACT: hlfsimko@liverpool.ac.uk or drigden@liverpool.ac.uk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Crystal structure of the nipah virus phosphoprotein tetramerization domain.

The Nipah virus phosphoprotein (P) is multimeric and tethers the viral polymerase to the nucleocapsid. We present the crystal structure of the multimerization domain of Nipah virus P: a long, parallel, tetrameric, coiled coil with a small, α-helical cap structure. Across the paramyxoviruses, these domains share little sequence identity yet are similar in length and structural organization, suggesting a common requirement for scaffolding or spatial organization of the functions of P in the virus life cycle.

Aquaponics in schools: Hands-on learning about healthy eating and a healthy planet.

Our food system is giving rise to a growing social, health and environmental crisis. Much of the food consumed in the United Kingdom is cheap, nutrient-poor and highly processed, leading to under-consumption of essential foods such as grains, beans, vegetables and fruit. This has contributed to a rise in diet-related diseases, with approximately 22% of primary school leavers being overweight or obese. Food production is unsustainable with agriculture responsible for 10% of the UK's greenhouse gas emissions and intensive farming practices have led to a significant loss of soil carbon and a decline in biodiversity. COVID-19 increased inequalities in our food system. Therefore, there is an urgent need for interventions to counteract these adverse social, health and environmental impacts. Education can play a crucial role as an intervention to address challenges in the food system. We tested an innovative school initiative using portable aquaponic pods and aligned to the national curriculum, to engage pupils in food production and foster learning about sustainability, climate change and healthy eating. The evaluation, based on teacher surveys, aquapod chart data, student blogs and postcards and feedback from the development team, revealed positive impacts on students' environmental awareness, as well as sustainability and practical food production knowledge. However, the programme encountered logistical challenges and we therefore highlight future improvements to produce a curriculum programme that can be delivered at scale to enhance food education and empower pupils to drive the agenda on tackling food sustainability and climate change.

MrParse: finding homologues in the PDB and the EBI AlphaFold database for molecular replacement and more.

Crystallographers have an array of search-model options for structure solution by molecular replacement (MR). The well established options of homologous experimental structures and regular secondary-structure elements or motifs are increasingly supplemented by computational modelling. Such modelling may be carried out locally or may use pre-calculated predictions retrieved from databases such as the EBI AlphaFold database. MrParse is a new pipeline to help to streamline the decision process in MR by consolidating bioinformatic predictions in one place. When reflection data are provided, MrParse can rank any experimental homologues found using eLLG, which indicates the likelihood that a given search model will work in MR. Inbuilt displays of predicted secondary structure, coiled-coil and transmembrane regions further inform the choice of MR protocol. MrParse can also identify and rank homologues in the EBI AlphaFold database, a function that will also interest other structural biologists and bioinformaticians.

CCP4 Cloud for structure determination and project management in macromolecular crystallography.

Nowadays, progress in the determination of three-dimensional macromolecular structures from diffraction images is achieved partly at the cost of increasing data volumes. This is due to the deployment of modern high-speed, high-resolution detectors, the increased complexity and variety of crystallographic software, the use of extensive databases and high-performance computing. This limits what can be accomplished with personal, offline, computing equipment in terms of both productivity and maintainability. There is also an issue of long-term data maintenance and availability of structure-solution projects as the links between experimental observations and the final results deposited in the PDB. In this article, CCP4 Cloud, a new front-end of the CCP4 software suite, is presented which mitigates these effects by providing an online, cloud-based environment for crystallographic computation. CCP4 Cloud was developed for the efficient delivery of computing power, database services and seamless integration with web resources. It provides a rich graphical user interface that allows project sharing and long-term storage for structure-solution projects, and can be linked to data-producing facilities. The system is distributed with the CCP4 software suite version 7.1 and higher, and an online publicly available instance of CCP4 Cloud is provided by CCP4.

Extending the scope of coiled-coil crystal structure solution by AMPLE through improved ab initio modelling.

The phase problem remains a major barrier to overcome in protein structure solution by X-ray crystallography. In recent years, new molecular-replacement approaches using ab initio models and ideal secondary-structure components have greatly contributed to the solution of novel structures in the absence of clear homologues in the PDB or experimental phasing information. This has been particularly successful for highly α-helical structures, and especially coiled-coils, in which the relatively rigid α-helices provide very useful molecular-replacement fragments. This has been seen within the program AMPLE, which uses clustered and truncated ensembles of numerous ab initio models in structure solution, and is already accomplished for α-helical and coiled-coil structures. Here, an expansion in the scope of coiled-coil structure solution by AMPLE is reported, which has been achieved through general improvements in the pipeline, the removal of tNCS correction in molecular replacement and two improved methods for ab initio modelling. Of the latter improvements, enforcing the modelling of elongated helices overcame the bias towards globular folds and provided a rapid method (equivalent to the time requirements of the existing modelling procedures in AMPLE) for enhanced solution. Further, the modelling of two-, three- and four-helical oligomeric coiled-coils, and the use of full/partial oligomers in molecular replacement, provided additional success in difficult and lower resolution cases. Together, these approaches have enabled the solution of a number of parallel/antiparallel dimeric, trimeric and tetrameric coiled-coils at resolutions as low as 3.3 Å, and have thus overcome previous limitations in AMPLE and provided a new functionality in coiled-coil structure solution at lower resolutions. These new approaches have been incorporated into a new release of AMPLE in which automated elongated monomer and oligomer modelling may be activated by selecting `coiled-coil' mode.

Using Phaser and ensembles to improve the performance of SIMBAD.

The conventional approach to search-model identification in molecular replacement (MR) is to screen a database of known structures using the target sequence. However, this strategy is not always effective, for example when the relationship between sequence and structural similarity fails or when the crystal contents are not those expected. An alternative approach is to identify suitable search models directly from the experimental data. SIMBAD is a sequence-independent MR pipeline that uses either a crystal lattice search or MR functions to directly locate suitable search models from databases. The previous version of SIMBAD used the fast AMoRe rotation-function search. Here, a new version of SIMBAD which makes use of Phaser and its likelihood scoring to improve the sensitivity of the pipeline is presented. It is shown that the additional compute time potentially required by the more sophisticated scoring is counterbalanced by the greater sensitivity, allowing more cases to trigger early-termination criteria, rather than running to completion. Using Phaser solved 17 out of 25 test cases in comparison to the ten solved with AMoRe, and it is shown that use of ensemble search models produces additional performance benefits.

Molecular replacement using structure predictions from databases.

Molecular replacement (MR) is the predominant route to solution of the phase problem in macromolecular crystallography. Where the lack of a suitable homologue precludes conventional MR, one option is to predict the target structure using bioinformatics. Such modelling, in the absence of homologous templates, is called ab initio or de novo modelling. Recently, the accuracy of such models has improved significantly as a result of the availability, in many cases, of residue-contact predictions derived from evolutionary covariance analysis. Covariance-assisted ab initio models representing structurally uncharacterized Pfam families are now available on a large scale in databases, potentially representing a valuable and easily accessible supplement to the PDB as a source of search models. Here, the unconventional MR pipeline AMPLE is employed to explore the value of structure predictions in the GREMLIN and PconsFam databases. It was tested whether these deposited predictions, processed in various ways, could solve the structures of PDB entries that were subsequently deposited. The results were encouraging: nine of 27 GREMLIN cases were solved, covering target lengths of 109-355 residues and a resolution range of 1.4-2.9 Å, and with target-model shared sequence identity as low as 20%. The cluster-and-truncate approach in AMPLE proved to be essential for most successes. For the overall lower quality structure predictions in the PconsFam database, remodelling with Rosetta within the AMPLE pipeline proved to be the best approach, generating ensemble search models from single-structure deposits. Finally, it is shown that the AMPLE-obtained search models deriving from GREMLIN deposits are of sufficiently high quality to be selected by the sequence-independent MR pipeline SIMBAD. Overall, the results help to point the way towards the optimal use of the expanding databases of ab initio structure predictions.

Electron delocalization in the metallabenzenes: a computational analysis of ring currents.

Many metallabenzene complexes appear to exhibit an enhanced thermodynamic stability which has been attributed to the concept of aromaticity. Analysis of the ring currents induced by a magnetic field, either by direct visualization or by considering nuclear or nucleus-independent chemical shielding values (NMR or NICS), have become useful theoretical tools to characterize the aromaticity of many molecules involving the main group elements. We have analyzed 21 metallabenzenes using variations of these techniques, which take account of the large core and metal orbital contributions which often lead to transition-metal-containing systems exhibiting anomalous shielding values. Analysis of individual orbital contributions to both the ring currents and chemical shielding values based upon the ipsocentric and CSGT (continuous set of gauge transformations) methods has shown that complexes such as the 18 electron Ir or Rh(C 5H 5)(PH 3) 2Cl 2 molecules should be classed as aromatic, whereas the 16 electron complexes such as Os or Ru(C 5H 5)(PH 3) 2Cl 2 should not, despite having the same occupancy of pi-MOs. The differences can be directly attributed to the HOMO/LUMO b 2 in-plane (d xy ) molecular orbital, which, when unoccupied, is available to disrupt the delocalized currents typical of aromatic systems. A range of Pd and Pt metallabenzenes with cyclopentadienyl and phosphine ligands is also discussed as having aromatic and nonaromatic character, respectively.

Ensembles generated from crystal structures of single distant homologues solve challenging molecular-replacement cases in AMPLE.

Molecular replacement (MR) is the predominant route to solution of the phase problem in macromolecular crystallography. Although routine in many cases, it becomes more effortful and often impossible when the available experimental structures typically used as search models are only distantly homologous to the target. Nevertheless, with current powerful MR software, relatively small core structures shared between the target and known structure, of 20-40% of the overall structure for example, can succeed as search models where they can be isolated. Manual sculpting of such small structural cores is rarely attempted and is dependent on the crystallographer's expertise and understanding of the protein family in question. Automated search-model editing has previously been performed on the basis of sequence alignment, in order to eliminate, for example, side chains or loops that are not present in the target, or on the basis of structural features (e.g. solvent accessibility) or crystallographic parameters (e.g. B factors). Here, based on recent work demonstrating a correlation between evolutionary conservation and protein rigidity/packing, novel automated ways to derive edited search models from a given distant homologue over a range of sizes are presented. A variety of structure-based metrics, many readily obtained from online webservers, can be fed to the MR pipeline AMPLE to produce search models that succeed with a set of test cases where expertly manually edited comparators, further processed in diverse ways with MrBUMP, fail. Further significant performance gains result when the structure-based distance geometry method CONCOORD is used to generate ensembles from the distant homologue. To our knowledge, this is the first such approach whereby a single structure is meaningfully transformed into an ensemble for the purposes of MR. Additional cases further demonstrate the advantages of the approach. CONCOORD is freely available and computationally inexpensive, so these novel methods offer readily available new routes to solve difficult MR cases.

Recent developments in MrBUMP: better search-model preparation, graphical interaction with search models, and solution improvement and assessment.

Increasing sophistication in molecular-replacement (MR) software and the rapid expansion of the PDB in recent years have allowed the technique to become the dominant method for determining the phases of a target structure in macromolecular X-ray crystallography. In addition, improvements in bioinformatic techniques for finding suitable homologous structures for use as MR search models, combined with developments in refinement and model-building techniques, have pushed the applicability of MR to lower sequence identities and made weak MR solutions more amenable to refinement and improvement. MrBUMP is a CCP4 pipeline which automates all stages of the MR procedure. Its scope covers everything from the sourcing and preparation of suitable search models right through to rebuilding of the positioned search model. Recent improvements to the pipeline include the adoption of more sensitive bioinformatic tools for sourcing search models, enhanced model-preparation techniques including better ensembling of homologues, and the use of phase improvement and model building on the resulting solution. The pipeline has also been deployed as an online service through CCP4 online, which allows its users to exploit large bioinformatic databases and coarse-grained parallelism to speed up the determination of a possible solution. Finally, the molecular-graphics application CCP4mg has been combined with MrBUMP to provide an interactive visual aid to the user during the process of selecting and manipulating search models for use in MR. Here, these developments in MrBUMP are described with a case study to explore how some of the enhancements to the pipeline and to CCP4mg can help to solve a difficult case.

SIMBAD: a sequence-independent molecular-replacement pipeline.

The conventional approach to finding structurally similar search models for use in molecular replacement (MR) is to use the sequence of the target to search against those of a set of known structures. Sequence similarity often correlates with structure similarity. Given sufficient similarity, a known structure correctly positioned in the target cell by the MR process can provide an approximation to the unknown phases of the target. An alternative approach to identifying homologous structures suitable for MR is to exploit the measured data directly, comparing the lattice parameters or the experimentally derived structure-factor amplitudes with those of known structures. Here, SIMBAD, a new sequence-independent MR pipeline which implements these approaches, is presented. SIMBAD can identify cases of contaminant crystallization and other mishaps such as mistaken identity (swapped crystallization trays), as well as solving unsequenced targets and providing a brute-force approach where sequence-dependent search-model identification may be nontrivial, for example because of conformational diversity among identifiable homologues. The program implements a three-step pipeline to efficiently identify a suitable search model in a database of known structures. The first step performs a lattice-parameter search against the entire Protein Data Bank (PDB), rapidly determining whether or not a homologue exists in the same crystal form. The second step is designed to screen the target data for the presence of a crystallized contaminant, a not uncommon occurrence in macromolecular crystallography. Solving structures with MR in such cases can remain problematic for many years, since the search models, which are assumed to be similar to the structure of interest, are not necessarily related to the structures that have actually crystallized. To cater for this eventuality, SIMBAD rapidly screens the data against a database of known contaminant structures. Where the first two steps fail to yield a solution, a final step in SIMBAD can be invoked to perform a brute-force search of a nonredundant PDB database provided by the MoRDa MR software. Through early-access usage of SIMBAD, this approach has solved novel cases that have otherwise proved difficult to solve.

Approaches to ab initio molecular replacement of α-helical transmembrane proteins.

α-Helical transmembrane proteins are a ubiquitous and important class of proteins, but present difficulties for crystallographic structure solution. Here, the effectiveness of the AMPLE molecular replacement pipeline in solving α-helical transmembrane-protein structures is assessed using a small library of eight ideal helices, as well as search models derived from ab initio models generated both with and without evolutionary contact information. The ideal helices prove to be surprisingly effective at solving higher resolution structures, but ab initio-derived search models are able to solve structures that could not be solved with the ideal helices. The addition of evolutionary contact information results in a marked improvement in the modelling and makes additional solutions possible.

Residue contacts predicted by evolutionary covariance extend the application of ab initio molecular replacement to larger and more challenging protein folds.

For many protein families, the deluge of new sequence information together with new statistical protocols now allow the accurate prediction of contacting residues from sequence information alone. This offers the possibility of more accurate ab initio (non-homology-based) structure prediction. Such models can be used in structure solution by molecular replacement (MR) where the target fold is novel or is only distantly related to known structures. Here, AMPLE, an MR pipeline that assembles search-model ensembles from ab initio structure predictions ('decoys'), is employed to assess the value of contact-assisted ab initio models to the crystallographer. It is demonstrated that evolutionary covariance-derived residue-residue contact predictions improve the quality of ab initio models and, consequently, the success rate of MR using search models derived from them. For targets containing ß-structure, decoy quality and MR performance were further improved by the use of a ß-strand contact-filtering protocol. Such contact-guided decoys achieved 14 structure solutions from 21 attempted protein targets, compared with nine for simple Rosetta decoys. Previously encountered limitations were superseded in two key respects. Firstly, much larger targets of up to 221 residues in length were solved, which is far larger than the previously benchmarked threshold of 120 residues. Secondly, contact-guided decoys significantly improved success with ß-sheet-rich proteins. Overall, the improved performance of contact-guided decoys suggests that MR is now applicable to a significantly wider range of protein targets than were previously tractable, and points to a direct benefit to structural biology from the recent remarkable advances in sequencing.

Routine phasing of coiled-coil protein crystal structures with AMPLE.

Coiled-coil protein folds are among the most abundant in nature. These folds consist of long wound α-helices and are architecturally simple, but paradoxically their crystallographic structures are notoriously difficult to solve with molecular-replacement techniques. The program AMPLE can solve crystal structures by molecular replacement using ab initio search models in the absence of an existent homologous protein structure. AMPLE has been benchmarked on a large and diverse test set of coiled-coil crystal structures and has been found to solve 80% of all cases. Successes included structures with chain lengths of up to 253 residues and resolutions down to 2.9 Å, considerably extending the limits on size and resolution that are typically tractable by ab initio methodologies. The structures of two macromolecular complexes, one including DNA, were also successfully solved using their coiled-coil components. It is demonstrated that both the ab initio modelling and the use of ensemble search models contribute to the success of AMPLE by comparison with phasing attempts using single structures or ideal polyalanine helices. These successes suggest that molecular replacement with AMPLE should be the method of choice for the crystallo-graphic elucidation of a coiled-coil structure. Furthermore, AMPLE may be able to exploit the presence of a coiled coil in a complex to provide a convenient route for phasing.