Atomistic modelling of scattering data in the Collaborative Computational Project for Small Angle Scattering (CCP-SAS)1
aDepartment of Structural and Molecular Biology, University College London, Darwin Building, Gower Street, London WC1E 6BT, UK, bCenter for Neutron Research, National Institute of Standards and Technology, Gaithersburg, MD 20899-8562, USA, cDepartment of Biochemistry, University of Texas Health Science Center at San Antonio, San Antonio, TX 78229-3900, USA, dDepartment of Biochemistry and Molecular Biophysics, Kansas State University, Manhattan, KS 66506, USA, eDepartment of Biology, Illinois Institute of Technology, 3101 S. Dearborn, Chicago, IL 60616, USA, fPharmacy Department, Franklin-Wilkins Building, King's College London, 150 Stamford Street, London SE1 9NH, UK, gDepartment of Chemistry, University of Bath, Claverton Down, Bath BA2 7AY, UK, hSchool of Biosciences, University of Nottingham, Sutton Bonington Campus, Leicestershire LE12 5RD, UK, iResearch Complex at Harwell, STFC Rutherford Appleton Laboratory, Harwell Campus, Didcot, Oxfordshire OX11 0FA, UK, jISIS Facility, STFC Rutherford Appleton Laboratory, Harwell Campus, Didcot, Oxfordshire OX11 0QX, UK, kDiamond Light Source Ltd, Diamond House, Harwell Science and Innovation Campus, Chilton, Didcot, Oxfordshire OX11 0DE, UK, and lDepartment of Chemistry, The University of Tennessee, Knoxville, TN 37996-1600, USA
*Correspondence e-mail: email@example.com, firstname.lastname@example.org
The capabilities of current computer simulations provide a unique opportunity to model small-angle scattering (SAS) data at the atomistic level, and to include other structural constraints ranging from molecular and atomistic energetics to crystallography, electron microscopy and NMR. This extends the capabilities of solution scattering and provides deeper insights into the physics and chemistry of the systems studied. Realizing this potential, however, requires integrating the experimental data with a new generation of modelling software. To achieve this, the CCP-SAS collaboration (http://www.ccpsas.org/) is developing open-source, high-throughput and user-friendly software for the atomistic and coarse-grained molecular modelling of scattering data. Robust state-of-the-art molecular simulation engines and molecular dynamics and Monte Carlo force fields provide constraints to the solution structure inferred from the small-angle scattering data, which incorporates the known physical chemistry of the system. The implementation of this software suite involves a tiered approach in which GenApp provides the deployment infrastructure for running applications on both standard and high-performance computing hardware, and SASSIE provides a workflow framework into which modules can be plugged to prepare structures, carry out simulations, calculate theoretical scattering data and compare results with experimental data. GenApp produces the accessible web-based front end termed SASSIE-web, and GenApp and SASSIE also make community SAS codes available. Applications are illustrated by case studies: (i) inter-domain flexibility in two- to six-domain proteins as exemplified by HIV-1 Gag, MASP and ubiquitin; (ii) the hinge conformation in human IgG2 and IgA1 antibodies; (iii) the complex formed between a hexameric protein Hfq and mRNA; and (iv) synthetic `bottlebrush' polymers.
Small-angle X-ray scattering (SAXS) and neutron scattering (SANS) are diffraction techniques for investigating a broad range of science. Here, we are particularly interested in their use in investigations of the structural properties of biomaterials, including proteins, nucleic acids and polysaccharides, and soft condensed matter systems, including synthetic polymers, micelles and liquid crystals (Perkins et al., 2008, 2011; Svergun et al., 2013; Higgins & Benoit, 1996; Gabrys, 2000). Like macromolecular crystallography, SAXS and SANS experiments have benefited hugely from major infrastructural investment in high-brilliance multiuser X-ray synchrotrons and neutron sources during the past two decades. Instrumental improvements that include higher beam intensities, the automation of data collection and more sensitive detector technologies have vastly increased both the throughput of samples and the quality of SAXS and SANS data (Hura et al., 2009; Pernot et al., 2013; Round et al., 2015; Heenan et al., 2011; Dewhurst et al., 2016).
As these are low-resolution techniques, detailed analysis of small-angle scattering (SAS) data requires the use of as much a priori knowledge about the system as possible. Thus SAS analysis often begins with semi-quantitative techniques such as radius of gyration RG fits or pair density distribution P(r) functions (Feigin & Svergun, 1987). The traditional and most straightforward approach to quantitative structural analysis of SAS data makes assumptions about the class of structure under investigation. The SAS from that object (the form factors and/or structure factors) is then calculated analytically using the parameters of the class (e.g. radius and density/composition for a sphere) as the unknown pieces of information to extract from the data. A large number of form factors have been derived and new ones continue to be published, but such analytical solutions are restricted to shape classes with sufficiently simple symmetry. In order to address the need to fit data from objects with little or no symmetry (ubiquitous in biomaterials such as proteins), non-atomistic real space and ab initio approaches were developed, most notably in the ATSAS suite of programs from EMBL Hamburg in which shape envelopes from spherical harmonics or assemblies of small spheres are used to fit experimental SAS data (Petoukhov et al., 2012; https://www.embl-hamburg.de/biosaxs/software.html). However, other than the minimal assumptions built into the shape optimization algorithm, no a priori information is used other than that gained from long experience and expert knowledge, often limiting the level of detail accessible. Increases in computer power and sophistication have led a number of groups to use known atomic coordinates (obtained for example from a crystal structure) to create atomistic models of macromolecules to fit SAXS and SANS data, including the use of automated fit procedures; other groups add further a priori information by including some type of energetic constraint (reviewed by Perkins et al., 2008, 2009, 2011; Lipfert & Doniach, 2007; Putnam et al., 2007; Rambo & Tainer, 2013). Indeed, despite the inherently low resolution of SAS and the loss of information due to orientational averaging, atomistic level structural models can be derived, with the bonding and energetic constraints they impose, when advanced molecular simulation is used to derive physically plausible configurations. These models are consistent not only with the scattering data but with at least some of the known physical chemistry of the systems, and provide functional insights for tests by further experimentation. An extensive but not exhaustive list of SAS data fitting packages, organized into categories, is available for download on the canSAS-maintained SAS portal at http://smallangle.org/content/Software.
The use of atomistic modelling to inform SAS analysis dates back several decades. For example, the first atomistic modelling, based on crystal structures but constrained by scattering curves, of antibody structures for the three Fab and Fc fragments (Fig. 1) showed that only a limited number of fragment conformations would fit the X-ray or neutron scattering curves (Perkins et al., 1991; Mayans et al., 1995). The first automation of atomistic scattering modelling was reported for the X-ray and neutron scattering fits for human IgA1 antibody using the SCT and SCTPL software tools with a commercial molecular dynamics (MD) package (Boehm et al., 1999). The hinge (or linkers) between crystal structures for the Fab and Fc regions were moved to explore relevant configuration space using MD to generate 12 000 trial full structures, of which only 102 gave good X-ray and neutron fits. This resulted in the first atomistic solution scattering structure to be deposited in the Protein Data Bank (PDB code 1iga). Since that time, SCT and SCTPL modelling has resulted in 77 such structures (Wright & Perkins, 2015), including the antibody classes of adaptive immunity (Fig. 1), the complement proteins of innate immunity with as many as 30 small domains, and linear anionic oligosaccharides containing up to 36 carbohydrate rings (Perkins et al., 2008). When comparisons were made, these structures compared well with those from other methods, such as protein crystallography (Perkins et al., 2008, 2011). During the past decade, several other groups have also pursued such modelling approaches (Whitten et al., 2008; Pelikan et al., 2009; Yang et al., 2009, 2010; Schneidman-Duhovny et al., 2010, 2016; Poitevin et al., 2011; Różycki et al., 2011; Evrard et al., 2011; Ihms & Foster, 2015; Chen & Hub, 2015; Jimenez-Garcia et al., 2015; Knight & Hub, 2015).
Despite these efforts, there has still been no significant shift toward atomistic or coarse grained modelling for SAXS and SANS. Further, there have been very few efforts to keep pace with rapidly evolving simulation methods through increases in computer power, coupled with new realistic force fields and robust sophisticated simulation engines. To begin addressing this issue, a modelling framework, termed SASSIE, emerged in 2004 at the NIST Center for Neutron Research (Datta et al., 2007; https://sassie-web.chem.utk.edu/sassie2/). SASSIE was developed to provide a general modular framework that enabled modern simulation methods to be applied to model scattering data using physical constraints (Curtis et al., 2012). It is important to understand that SASSIE is a framework built on a plugin architecture. It is meant to be agnostic to the particular MD engine, force field, type of material or approach to solving the molecular structure. Any limitations in that regard come strictly from what modules are available. In fact, many if not all of the atomistic efforts developed so far could in principle be wrapped into modules that take advantage of the integrated SASSIE workflow. The other new element in SASSIE was the incorporation of Monte Carlo (MC) simulation methods to create ensembles of biomolecular structures by sampling user-selected backbone dihedral angles to model experimental X-ray and neutron data. Without this advance, the generation of atomistic structures using modern force-field-based simulations could take months or even be inaccessible. Because robust, force-field-based structures are used throughout the SASSIE workflow, most modern simulation packages can be used to model SAS data in detail. Through the generation of robust and complete physics models using best-practice simulation methods, the resulting ensemble of structures can be fitted to SAXS and SANS data within the integrated workflow.
Despite the advances through the development of SASSIE, four key issues remain to make various SAS modelling problems more tractable. Firstly, regarding accessibility, because the SASSIE framework supports an array of plugin modules, the end user installation of the software could be frustrating and require a level of technical knowledge usually lacking in a typical non-expert SAS end user. New versions can quickly overwhelm the users and the small development team. Secondly, it became quickly evident that the largest barrier for the non-simulation end users was the need for expert knowledge to prepare the correct starting trial structural models and the associated protein structure files. The general difficulty for biological non-experts is creating a well constructed starting structure, while that for soft matter non-experts is that the appropriate force fields may not be available. Thirdly, the increasing complexity of systems of interest is leading to an increasing need to run these simulations on high-performance computing (HPC) resources, something outside the skill set of most SAS users. Fourthly, it has become clear that developing and maintaining the SASSIE framework, while also wrapping or developing the possible and desired packages and tools as new modules plugged into the SASSIE framework, is beyond the ability of a single small group to manage. This requires a larger community effort.
The 17 Collaborative Computational Projects (CCPs) in the UK (as of August 2016; http://www.ccp.ac.uk/) provide a software infrastructure to build individual research projects and to maintain and distribute code libraries. In order to reveal how atomic level molecular structures in biological or soft matter systems account for experimental scattering data, the Collaborative Computational Project for Small Angle Scattering (CCP-SAS), jointly funded by the EPSRC research council in the UK and the National Science Foundation in the USA, was created in 2012 to address these issues of access and long-term sustainability. The specific initial goals of the consortium were to (i) significantly lower the barrier for bench scientists to access the power of high-end state-of-the-art molecular modelling and computational chemistry tools; (ii) provide a user-friendly software environment that integrates SAS data with those tools for purposes of structural refinement, further informed by data from complementary techniques such as analytical ultracentrifugation, electron microscopy or NMR; and (iii) build a long-term development and maintenance support structure through community development and engagement with large-scale SAS user facilities as well as other CCPs. Here, we provide an overview of the CCP-SAS project, focusing heavily on its current core activities. These comprise the development of a new GenApp infrastructure for deployment of computational code, the ongoing development of the SASSIE framework and its implementation as SASSIE-web powered by the new GenApp package to provide a web front end and HPC back end, and the ongoing development of the workflow of modules required to address molecular simulations, scattering calculators and the analyses of their output. For some soft matter systems, the extension of SASSIE to coarse-grained and hybrid methods (mixing shapes with atomistic structures) will be important. To illustrate some representative atomistic modelling workflows, we summarize applications of SASSIE-web to a broad range of systems in biology and soft matter (Fig. 1) (Datta et al., 2007; Nan et al., 2017, unpublished work; Castañeda, Chaturvedi et al., 2016; Castañeda, Dixon et al., 2016; Clark et al., 2013; Hui et al., 2015; Peng et al., 2014; Zhang et al., 2014).
The initial goal of CCP-SAS is to provide an open-source cloud-based software environment that not only makes clear how the modelling fit analyses were performed, and permits experimental teams to understand complex chemical interactions and structural organizations, but is flexible enough to incorporate additional different experimental constraints into the modelling workflow. The CCP-SAS project also aims to provide documentation and training, and ultimately to foster a sustainable community of users. This user base includes experimental research groups, software developers and instrument scientists at multiuser scattering facilities.
Current CCP-SAS activities include these nine tasks: development of the GenApp infrastructure and SASSIE framework; deployment of SASSIE as SASSIE-web to the community; wrapping existing code and developing new code as new modules for SASSIE; developing new methods for eventual incorporation as new modules in SASSIE; working with members of the SAS community to implement their relevant methods and codes into the SASSIE framework; providing help and guidance to members of the SAS community to wrap their standalone codes using GenApp for separate web deployment outside of SASSIE-web; where feasible and reasonable, hosting such separate web applications on the CCP-SAS cluster for the benefit of the community; running tutorials and workshops; and working to engage various community stakeholders.
Two core principles of CCP-SAS are to use both existing and open-source software as much as possible. If a critical non-open-source component is needed, it can be incorporated, but then an alternative open-source solution is identified to replace this as quickly as possible. This policy accommodates the drive for open-source software for proper validation and transparency increasingly requested by funding bodies and helps engage community support. Thus all CCP-SAS software, including SASSIE and GenApp, is freely available and open source. While one closed-source package currently remains (August 2016), this will be removed as soon as the alternative modules are validated.
The GenApp infrastructure was developed to simplify the deployment of CCP-SAS software (Brookes et al., 2015). Common issues addressed by GenApp include easing the deployment of a workflow of modules, support for legacy codes and the reduction of dependencies on dedicated software teams. This is achieved by enabling the generation of web-based and standalone graphical user interface (GUI) applications over the same underlying executable software while providing transparent access to back end computational resources and connections to high-performance computing gateways (Fig. 2a). Long-term sustainability questions are addressed by decoupling the GUI and back end interfaces from the core computational codes, such as the SASSIE suite being developed. GenApp is thus the core technology to address the accessibility issues as well as the long-term sustainability issues.
In GenApp, an application is defined as a collection of executable modules which are presented through a common user interface (Fig. 2a). This provides a powerful paradigm to combine both existing and new codes in order to perform novel workflows or develop different types of modelling applications. The addition of a module in GenApp is simple, and only requires the writing of a short JSON wrapper (a module) to detail the input and output, and the editing of two JSON files, one to specify where the module should appear in the applications menu system, and the other to specify how the application itself is to be presented. The modules themselves can be written in any supported language, independent of the choice of the target GUI implementation. Separating the scientific code from the GUI not only facilitates the linking of component modules into larger workflows and applications, but also reduces the burden in supporting legacy codes. GenApp also facilitates the creation of applications as web servers or gateways. This includes remote file management and the execution and management of lengthy non-interactive jobs. The latter capability, provided through integration with Apache Airavata (https://airavata.apache.org/), allows GenApp applications to harness a range of high-performance computing resources including local clusters, supercomputers, national grids, and academic and commercial clouds. We anticipate that GenApp will be useful to generate a wide range of scientific applications beyond the scope covered by CCP-SAS.
GenApp was designed to be generic, and thus its power is available to any developers seeking to take advantage of the ease of deployment and transparent access to high-end computing resources it offers. GenApp modules can be part of a module developer's standalone application or hosted in CCP-SAS computer resources as a public web-based science gateway. GenApp web applications can in principle be deployed on any cloud resource, and instances have been tested on XSEDE (https://www.xsede.org/) and AWS (https://aws.amazon.com/). The developed GenApp module can be added to our open repository. Currently the project is working with SAS developers including WillItFit (Pedersen et al., 2013) and QuaFit (Spinozzi & Beltramini, 2012). Both packages are deployed for alpha testing as web applications hosted on our CCP-SAS resources. Interested parties are invited to send an email message to email@example.com.
The aim of SASSIE-web is to allow experimentalists (including novice users) to construct their own modelling workflows from a set of simulation and analysis modules, then run them transparently on centrally maintained back end resources for scattering curve comparisons, from nothing more than a web browser (https://sassie-web.chem.utk.edu/sassie2/) (Fig. 3). The provision of a web interface avoids the need for users to install and maintain large complex software on their own machines, and facilitates the provision of a high-performance computing back end to accelerate the computationally expensive steps of the modelling process. The SASSIE-web menu organizes the workflow in terms of six sets of modules (Fig. 3a): (i) Tools, which includes utilities to predict scattering length densities, interpolate experimental scattering data files when required, and extract or merge macromolecular structures; (ii) Build, which includes utilities to check PDB-formatted coordinate files; (iii) Interact, which provides a molecular viewer to present an interactive display of a specified structure using JSMol (Hanson et al., 2013); (iv) Simulate, which provides the modules that create the representative ensemble of trial structures for test against the data (Fig. 3b); (v) Calculate, which provides a range of scattering curve calculators; and (vi) Analyze, which determines the goodness of fit between the simulated and experimental scattering curves in order to identify the best-fit scattering structure (Fig. 3c) and provide visualizations to display the trial structures and the best-fit subset of these as envelopes.
The modular design of SASSIE not only gives the user the freedom to employ any combination of existing modules but also allows them to plug in new modules, and import coordinate models generated with other packages at any stage of the workflow. This modular nature of SASSIE, combined with the ease of deployment and end user accessibility, makes SASSIE-web an attractive option for SAS computational groups wishing to contribute their codes. For example, the Capriqorn software to calculate scattering curves from molecular simulations with explicit water models is being integrated into the SASSIE framework (Köfinger & Hummer, 2013). Interested parties are invited to email firstname.lastname@example.org.
The SASSIE workflow is summarized in Fig. 2(b). The modelling of scattering data is crucially dependent on correct starting atomistic models. Even though over 121 000 structures (August 2016) are available in the PDB, it does not follow that these are ready for molecular simulation and scattering curve calculations. Common problems are gaps in the protein structure caused through disorder (especially at surface loops), errors in the amino acid sequence, missing structures such as incomplete glycan chains or N-terminal or C-terminal sequences, and missing or misnamed atoms (e.g. hydrogen atoms, carbon atoms and disulphide bridges). If the structure of interest is not available in the PDB, standard homology (comparative) modelling techniques which are necessarily outside SASSIE (such as MODELLER; Šali & Blundell, 1993) can be used to generate the input structure from the most closely related protein structure in the PDB. In this process, the amino acid sequence will need to be replaced by the sequence of interest. PDB-scan assesses whether a PDB file is ready for a scattering curve simulation, and where possible provides files enabling CHARMM force-field parameterization (MacKerell et al., 1998; Best et al., 2012). Scans provide information on missing atoms and residues and those not covered as standard by the CHARMM force field. PDB-scan also reports on whether symmetry information present in the PDB header can be used to create a dimer or higher-order oligomer that is the actual biological unit to be modelled. Suitable coordinate files can be derived from the output that are ready to be used by a wide range of simulation and modelling software packages, including those focused on soft matter systems. To complement the capabilities of PDB-scan, a new module termed PDB-Rx is in preparation to correct mistakes discovered by PDB-scan (Wright et al., 2016).
In SASSIE (Fig. 2b), the key stage of modelling SAXS and SANS data is the generation of ensembles of atomistic structures that sample the configuration space of physically realistic models. Early approaches used various MD or MC methods to vary the appropriate segment in the system of interest (Boehm et al., 1999; Datta et al., 2007; Khan et al., 2010). For biological work, to generate structural models of protein or nucleic acids rapidly, SASSIE-web offers dihedral angle MC simulations through the Markov sampling of backbone torsion angles in user-specified regions of the input model (Curtis et al., 2012). MC simulations can be performed on any PDB structure which contains all atoms in the model. However, in order to make use of the full range of simulation and analysis options in SASSIE, it is recommended that the input PDB file is prepared for MD simulation using the CHARMM force field. It is not necessary to obtain the CHARMM simulation package in order to perform this process. One common approach is to use the structure building tool PSFGEN which is distributed openly as a plugin with the VMD visualization program or the NAMD simulation package (Humphrey et al., 1996; Phillips et al., 2005). As an alternative, access to CHARMM force-field parameterization is provided by CHARMM-GUI (Jo et al., 2008). The starting input structure must be a complete structure without missing residues (see above) and atom and residue names must be compatible with those defined in the CHARMM force field (MacKerell et al., 1998).
The modelling strategy completely depends on the system of interest. During a typical simulation workflow for a multidomain protein with linkers to be varied between the domains, about three to six linker regions in the starting structure are sampled in the simulation. Depending on the system of interest, around 10 000 to 50 000 structures might be required to sample adequate configuration space for most problems; see Table 1 for examples. Since steric clashes can easily occur during the simulation, the avoidance of atomic overlap is achieved by specifying an overlap distance cutoff (typically 0.3 nm) and the atom name(s) to which this applies. Other options include the selection of simulated structures to remain within a fixed range of RG values and/or satisfy intra- and intermolecular distance constraints. Collections of output structures are stored in the DCD file format used by the CHARMM, NAMD and X-PLOR MD packages (this binary format stores multiple structures much more efficiently than text-based PDB files; Brunger, 1992). These files can also be visualized in many molecular viewers, such as VMD or Chimera (Humphrey et al., 1996; Pettersen et al., 2004). Presently, two interfaces to MC simulations are provided in the Simulate module, namely Monomer MC and Complex MC. As their names suggest, the former provides a simplified interface focusing on single-chain biosystems, while the latter facilitates the simulation of more complex topologies. A tutorial using the Monomer MC module, based on its original use case of the HIV-1 Gag protein (Fig. 1) (Datta et al., 2007), can be found at https://sassie-web.chem.utk.edu/sassie2/docs/sassie-web-quick-start/quick-start.html. While this example covered a workflow for a protein, the other simulation engines such as NAMD and CHARMM within SASSIE enable any molecular system to be simulated, including in particular soft matter systems, something that is the focus of ongoing work.
The outcome of the MC simulations is available to another module that uses energy minimization and MD to sample degrees of freedom not sampled in the MC trajectories from biomolecular models as parameterized in the CHARMM force field (Fig. 2b). NAMD (version 2.9) is used as the simulation engine (Phillips et al., 2005). A reference PDB file name is input, together with the matching starting structure in either PDB or DCD format, and the CHARMM topology (PSF) file. The four optional modes of operation are as follows: (i) minimization alone; (ii) minimization followed by MD; (iii) minimization followed by MD leading to a second round of minimization; and (iv) molecular simulation (energy minimization and/or MD) with a user supplied input file. Both the minimization and MD are performed using the generalized Born implicit solvent model. If a DCD file is selected as the input file then the simulations are run on each frame.
Structural models generated by the MC simulations can also be sent to the Torsional Angle MD (TAMD) module for refinement (Zhang et al., 2016, unpublished work). TAMD samples molecular configurations in torsion angle space, and allows the convenient specification of rigid domains and flexible degrees of freedom consistent with the MC sampling stage (Chen et al., 2005). For this, the ensemble generated by MC simulations is first sub-sampled to select representative configurations that provide a thorough coverage. Each selected configuration is then used to initiate TAMD simulations, which allows refinement of the local structural features and provides improved sampling of conformational degrees of freedom that are not included in the MC moves. Atomistic implicit solvent force fields available in CHARMM are used to provide a balance between computational efficiency and efficacy (Chen et al., 2008). By default, the module currently uses an efficient solvent-accessible surface area (SASA) implicit solvent model (Ferrara et al., 2002) that can handle proteins, nucleic acids and carbohydrates.
Theoretical scattering curves for modelled structures are computed from atomistic positions, such as via the Debye equation. As this requires the calculation of distances between every pair of atoms or scattering centres, the computing effort increases with the square of the number of atoms or centres, making this hugely time consuming. These computations become even more computer intensive if the pair distances are convoluted with the scattering length densities of each pair of scattering centres for neutron contrast variation work or for X-ray work with high and low electron densities. Scattering curve calculation can be accelerated by the use of coarse graining of the original atomic structures (resulting in fewer scattering centres), and the use of binning algorithms to reduce the number of distances to be processed. An alternative strategy is to use high-performance computing and graphics processing unit technology to accelerate the computations, both of which have been pursued within CCP-SAS. Several calculators are available in the SASSIE-web framework, such as CRYSOL and EM_to_SANS (Svergun et al., 2012; Curtis et al., 2012). Two of the most commonly used are described briefly here.
The SasCalc module (Watson & Curtis, 2013) calculates neutron and X-ray scattering profiles from a user-supplied structure file using an exact all-atom expression for the scattering intensity in which the orientations of the Q vectors are taken from a quasi-uniform spherical grid generated by the golden ratio. This `golden vector' method is currently configured to handle atomic trajectory input files (DCD or PDB). Our implementation of the `golden vector' method within SasCalc includes corrections for contrast for both X-ray and neutron scattering and harnesses graphical processing units for massively parallel calculations. The use of an atomistic scattering calculator is a vital step towards supporting the full range of soft matter systems beyond biological systems.
The SCT module (Wright & Perkins, 2015) first converts the atomistic structure into a coarse grained sphere model, where each sphere represents about 4–5 atoms, and then employs the Debye equation adapted to small spheres of diameter below the structural resolution of the scattering curves (about 1 nm diameter) to calculate the scattering curve. The simulations utilize single-density spheres. X-ray scattering simulations involve the addition of a hydration monolayer of water molecules in the sphere model creation step, and assume pinhole geometries with X-ray data that are already slit desmeared. Neutron scattering simulations for proteins in heavy water do not require a hydration shell (Perkins, 2001), but require a smearing correction for wavelength spread and beam divergence in the final scattering curve, as well as a linear buffer background correction for residual incoherent scattering from protons. As part of CCP-SAS, SCT was made open source and publicly available and is downloadable from https://github.com/dww100/sct (Wright & Perkins, 2015). The new SCT release in the SASSIE-web Calculate module (Fig. 2b) features an improved user interface (including the acceptance of DCD coordinate files) and modules which facilitate its use in modelling workflows (e.g. comparing the theoretical and experimental curves).
The final stage of the modelling (Fig. 2b) is of course the comparison of the simulated scattering profiles from the structural coordinates with the experimental SAXS and/or SANS scattering curves in order to identify the best-fit structures. In the Analyze module of SASSIE, the χ2 filter module offers one approach. This module compares the theoretical scattering profiles with the interpolated experimental data. The user supplies an input experimental data file containing three columns – Q, I(Q) and error of I(Q) at each Q value – together with the I(0) value from the Guinier RG analysis or the P(r) analysis. In the output, three mathematical options are provided to evaluate the quality of the comparison, namely the reduced χ2, χ2 and the Pearson χ2. To process the simulated scattering curves produced by SCT, the module that was originally termed SCTPL (Perkins et al., 2008, 2011) is now renamed as SCT Analyze to clarify its purpose. As with the χ2 filter module the user is given a choice of comparison metric (R factor or reduced χ2). Typical R factors for best-fit SAXS or SANS analyses are between 2 and 8%, and typical best-fit χ2 values are around 1–2 (Table 1).
The Analyze module also helps deduce the biological significance of the final best-fit structures (Fig. 2b). Starting from the DCD frames file or the PDB coordinate files, the density plot module generates files (using the Gaussian cube file format) with volumetric data. Often, this is used to visualize sub-ensembles of structures that give the best agreement with experimental data, in addition to views of all the trial structures generated. The resulting envelopes can be visualized in molecular viewing packages such as VMD. Given the atomistic nature of the best-fit structures, further analyses of these best-fit structures become possible; e.g. these may include the calculation of electrostatic surface charge effects.
The CCP-SAS web site with background, reports and publications is at http://www.ccpsas.org/. At its inception, CCP-SAS benefitted from resources from the NIST Center for Neutron Research in the USA, and the Diamond Light Source and the ISIS Pulsed Neutron Source in the UK. But, being a CCP, CCP-SAS aims to create a community of users and provide a training infrastructure, in addition to developing software that suits the need of experimentalists. Its outreach strategy involves regular meetings, tutorials and workshops with users at scattering facilities, and engaging wider audiences at conferences through targeted training activities. In terms of maintenance, most users of CCP-SAS software will not have direct access to support staff. Consequently, detailed yet comprehensible documentation is required. Each completed package or module developed within CCP-SAS has its own online documentation. The documentation outlines the elements of each module and its interface and how to use it. In addition, worked examples are provided with accompanying files based on previous successful projects. This documentation is provided online at https://sassie-web.chem.utk.edu/sassie2/docs/ via a prominent Docs tab on the home page (Fig. 3a). No matter how well documented and tested the software is, there will always be new user issues. A CCP-SAS Google Group provides user support (https://groups.google.com/forum/#!forum/madscatt). This is linked to the SASSIE-web interface so that users can directly report issues in specific modules, raise new features to be added and propose the best ways to tackle new projects. CCP-SAS is also actively engaged in supporting the larger SAS software developer community.
The first six examples of SASSIE-driven atomistic scattering modelling below illustrate the breadth of its applications in structural biology and how biologically useful information is obtained. Both single-chain proteins and multimers have been analysed in a range of SASSIE workflows (Table 1), but all driven from the same user interface. Common to all six biology projects (Fig. 1 and Table 1) is the definition of a correctly formatted initial structure (including the addition of hydrogen atoms), which is energy minimized with NAMD or CHARMM using the CHARMM27 or CHARMM36 force field (MacKerell et al., 1998). This initial structure is then subjected to MC simulation to generate an ensemble of physically relevant structures. Theoretical scattering profiles of each structure in the ensemble are compared with the experimental SAXS and/or SANS data to define the best-fit structures, from which new biological insight is obtained, as exemplified in Figs. 4 and 5 (§3.5). The resulting ensembles can then form the basis of further studies. The seventh example illustrates an application of SASSIE to a synthetic polymer system.
The first system to be studied by SASSIE was HIV-1 Gag, a long polypeptide chain which consists of four domains (MA, CA, p2 and NC) connected to four long flexible linkers. The human immunodeficiency virus type 1 (HIV-1) Gag polyprotein leads to the efficient assembly of virus-like particles in mammalian cells after this is cleaved by a viral protease. Gag is composed of three well defined immature proteins, matrix (MA), capsid (CA) and nucleocapsid (NC), alongside p2 whose structure is not well known (Fig. 1). These form virus-like particles when exposed to nucleic acids and their assembly will be governed by the solution structure of Gag.
Because molecular structures for the MA, CA and NC domains have been determined by crystallography and NMR, the solution structure of monomeric Gag could be modelled from SANS data and SASSIE modelling (Datta et al., 2007). The primary unknown was the long flexible linkers that join these three domains, including flexible linkers within both the CA and NC domains. The Monomer MC module was used to generate eight groups of 600 structures (totalling 4800) by varying specified main chain φ and ψ angles in five peptide linkers (marked with arrows in Fig. 1). This large ensemble of trial structures was output into a DCD file, and then these were energy minimized using CHARMM. The resulting scattering analysis of the 4800 models gave χ2 values that ranged between 1.2 for four groups of best-fit structures to as high as 30 for the other 16 groups of structures. While no single correct structure was identified, the key result from the best-fit ensembles showed that unbound Gag was folded over into a compact shape with the N-terminal MA and C-terminal NC domains close to each other, i.e. such a structure undergoes a conformational change when this is assembled into a virus. This modelling workflow is described in more detail in the online tutorial for SASSIE.
In a recent similar study, the examination of the bacterial single-stranded DNA binding protein from SAXS and SANS data showed that the protein's long disordered C-terminal tails were relatively collapsed around the well defined N-terminal core protein structure that binds single-stranded DNA, and compacted further upon binding single-stranded DNA, which was at variance with the previously hypothesized model (Green et al., 2016). To visualize this outcome by atomistic modelling, SASSIE was first used to generate the most compact atomistic structure possible for the core and tails. Next, 10 000 structures for the full-length protein were generated from MD, in which the tails were allowed to adopt all stereochemically permitted disordered conformations. Careful allowance for hydration was required for the SAXS fits; however, the SANS data could be directly compared with the unhydrated models because the surface hydration shell of bound water molecules does not contribute to the SANS data. Finally, the resulting curve fits showed that the most collapsed ensemble of tail structures best represented the experimental scattering curves.
Polyubiquitination is a post-translational modification of an intracellular protein by an ubiquitin dimer that signals major cellular events. Different signalling pathways result, depending on the isopeptide linkages formed between the C-terminus of the so-termed `distal' ubiquitin (Ub) and the ∊-amine of any one of seven lysine residues on the other `proximal' Ub. For example, Ub dimers linked by Lys48 or Lys63 mediate proteosomal degradation and DNA repair, respectively. The Ub dimer formed through Lys27, termed K27-Ub2 (Fig. 1), has unique biochemical properties (Castañeda, Chaturvedi et al., 2016; Castañeda, Dixon et al., 2016).
The ubiquitin dimer is joined by an isopeptide bond between two monomers. Of interest is that the structure and dynamics of K27-Ub2 were examined jointly by NMR structural constraints, SANS data and ensemble modelling by SASSIE (Castañeda, Dixon et al., 2016). The sparse ensemble selection method was used to determine representative conformational ensembles for K27-Ub2 for the NMR analyses. The ensembles of 23 000 sterically allowed structures were generated using the monomer MC routine in SASSIE to vary residues 72–76 of the distal Ub monomer that were connected to Lys27 in the proximal Ub monomer (arrow in Fig. 1). Residues 72–76 were considered to be flexible. Importantly, the fit for the residual dipolar couplings from NMR was significantly improved if two K27-Ub2 conformers and not one were considered. Independent SASSIE fits of the experimental SANS data using these 23 000 structures showed that a one-conformer ensemble gave good agreement, and two-conformer ensembles slightly improved the agreement. Similarity with the structure of the ligand-bound state of the Ub dimer linked via Lys48 suggested a possible receptor for K27-Ub2, which was then confirmed experimentally. The biological importance of this dimer was revealed by studying the interaction between K27-Ub2 and its receptor by molecular docking based on NMR signal perturbation and paramagnetic spin labelling. In this receptor complex, surface-exposed hydrophobic patches on each of the two Ub proteins formed a V-shaped groove in the dimer that interacted specifically with the receptor inside its V, thus rationalizing the distinct biochemistry of this particular Ub dimer.
The lectin pathway of the complement system in plasma is activated by complexes on pathogen surfaces that comprise a recognition component (MBL: mannose-binding lectin) that binds multivalently to mannose residues on the pathogen. MBL forms complexes with an MBL-associated serine protease (MASP) that leads to MASP activation in a manner that is unclear. MASP exists as a homodimer, with two six-domain monomers tightly bound in an antiparallel arrangement by their N-terminal domains (Fig. 1).
The dimeric MASP proteins are composed of two copies of six protein domains joined by short linkers; the dimer is formed by the tight noncovalent pairing of the N-terminal CUB–EGF domain pair (Fig. 1). The solution structure of the full MASP dimer with 12 domains was not known, and unravelling this structure was critical to understanding MASP activation (Nan et al., unpublished work). The SASSIE modelling of SAXS data for MASP was performed alongside crystallographic investigation of the same proteins and was achievable despite the large size of this protein. An initial linear six-domain monomer model was built using a combination of existing domain structures and new crystal structures for the three N-terminal domains that form the dimer. When full-length dimeric MASP was studied by analytical ultracentrifugation, the experimental sedimentation coefficients s20,w of 5.4–5.9 S were notably larger than those of 5.3–5.5 S calculated from the initial linear model. From the X-ray RG analyses, the experimental values of 7.5–7.9 nm for full-length MASP were likewise notably smaller than those of 8.2–9.0 nm calculated from the initial linear model. Both differences indicated that the solution structures of MASP were more bent than the initial linear model.
In the SASSIE workflow, the Monomer MC module was used to sample MASP configurations by varying one or four inter-domain linkers in the three- and six-domain proteins (arrows in Fig. 1). As many as 30 910 trial conformations (Table 1) were generated using maximum rotation steps for the peptide φ and ψ angles of up to 80°. The SCT scattering curves were calculated from coarse-grained sphere models using a cutoff of four atoms per sphere in a grid with a cube side of 0.530 nm. A hydration sphere shell corresponding to 0.3 g of water per gram was added to each unhydrated model. Three N-linked glycan chains are present in MASP. These could not be considered during the SASSIE modelling, which was performed without explicit glycans added. Once the best-fit MASP structures had been identified, glycans in extended conformations were added to these, whereupon the R factors were improved to reduced final values of 4.6–5.2% (Table 1). The SASSIE modelling showed that much improved curve fits resulted from bent full-length atomistic MASP structures, compared to the extended initial structure. This key result revealed that MASP existed as flexible structures in which the two SP domains at the tips of the MASP dimer were able to move towards each other. Although this hypothesis is not proven, the modelling suggests that the MASP domains are flexible and that the two SP domains at the tips of the MASP dimer may come sufficiently close to explain how MASP auto-activation may take place. The MASP example, being constrained by crystal structures, showed that as many as eight variable linkers can be analysed using SASSIE. The incorporation of greater numbers of variable linkers requires other approaches, such as that for intrinsically disordered proteins discussed above (Green et al., 2016).
IgG antibodies are central to the adaptive immune response against pathogens. As therapeutics, over 300 IgG monoclonal antibodies have been approved for clinical use. The four human IgG subclasses IgG1–IgG4 in serum differ primarily in their hinges, where their lengths are 15, 12, 62 and 12 amino acids, respectively. Atomistic antibody modelling by SASSIE is an ideal method to investigate how the Fab regions are connected to the Fc regions through two long flexible linkers (or hinges) (Fig. 1). Because the two Fab and Fc regions are largely independent of each other, antibody modelling is distinct from the examples of the linear two- to six-domain proteins above.
SASSIE was used to study the structure of a monoclonal human IgG2 antibody that was characterized by SANS (Clark et al., 2013). The atomistic modelling was initiated by a homology model for human IgG2 that was generated from a crystal structure for mouse IgG2a. This homology model was subjected to MC simulation by sampling three residues in the upper hinge in random rotational steps of up to 10°. Each of the resulting 56 511 conformations was subjected to energy minimization, followed by a generalized Born-implicit solvent MD simulation, and another round of energy minimization. From comparison between the experimental and calculated SANS curves, the standard plot of χ2 versus RG values showed a U-shaped distribution in which the best-fit structures appeared at the minimum (e.g. Fig. 4a below). Visual inspection showed that the conformational space about the hinge had been well sampled, but only a small set of conformations were in agreement with experiment with χ2 < 2. An asymmetric arrangement of the two Fab regions compared to the Fc region was identified. This ensemble of structures was consistent with the scattering data; however the configurations may or may not be energetically plausible. To complete this study, energetic information from simulations was used to refine the ensemble of best-fit structures for IgG2. The widely used Adaptive Poisson Boltzmann Solver implemented in SASSIE-web was used to calculate solvation free energies from the ensemble models (Baker et al., 2001; Dolinsky et al., 2004, 2007). These solvation energies acted as a further filter on acceptable models and produced a reduced subset of structures exhibiting lower free energies. The use of free-energy constraints meant that the final models corresponded to more physically reasonable structures. This software technology was able to identify specific interactions known to affect function and/or chemical stability, and illustrated a new approach made possible because of the modules in SASSIE.
In related studies that used the older SCT and SCTPL approach, atomistic analyses for monoclonal human IgG1 and IgG4 were based on known crystal structures for the Fab and Fc regions. MD simulations were used to vary the hinge conformations in order to interpret the SAXS and SANS data (Rayner et al., 2014, 2015). The outcome also revealed asymmetric IgG structures. The resulting atomistic structures explained why human IgG1 binds to its receptors and complement more readily than human IgG4.
IgA1 and IgA2 antibodies are important in mucosal immunity. IgA nephropathy is a leading cause of chronic kidney disease, in which the deposition of IgA1-containing immune complexes in the kidney often leads to renal failure. The structure of IgA1 is unusual in possessing two long 23-residue hinges between the Fab and Fc regions. These hinges are O-glycosylated with GalNAc.Gal.NeuNAc moieties, and these O-glycans are often found at reduced levels in patients with IgA nephropathy.
The SASSIE modelling of SAXS and SANS data (Fig. 4) investigated the impact of glycosylation on the IgA1 solution structure. To clarify whether variations in these O-glycans affect IgA1 function and disease, human IgA1 was studied with four different O-glycosylation levels (Hui et al., 2015). Analytical ultracentrifugation showed that all four IgA1 samples were monomeric with similar sedimentation coefficients s020,w. SAXS and SANS data in light and heavy water, respectively, for the four IgA1 samples revealed no conformational changes between the four IgA1 samples. Interestingly, the SANS data acquired in heavy water suggested that a reduction in O-glycan content was correlated with an increase in the propensity for IgA1 to aggregate, i.e. this may be related to the onset of IgA nephropathy. The SASSIE modelling workflow for IgA1 proceeded in two stages. First, a truncated IgA1 structure was modelled from crystal structures for each of the human IgA1 Fab and Fc regions, with the hinge and C-terminal regions modelled de novo. Two N-glycan chains at Asn263 in the Fc region in the crystal structure were also varied. This IgA1 structural model was energy minimized (Fig. 5b). MC simulations of the hinge conformations resulted in 172 833 trial models whose calculated scattering curves were compared to the SANS data to identify the best-fit truncated IgA1 models. Second, these best-fit truncated IgA1 structures were completed by adding structures for the two N-glycosylated C-terminal tailpieces obtained from MD simulations to give 146 484 trial models for full-length IgA1 (Fig. 4a). Principal component analysis identified four major IgA1 conformations. One of these conformations gave very good SAXS and SANS curve fits (Figs. 4b and 4c). Whilst no structural variation was found with differing glycosylation levels, in agreement with experiment, the addition of six O-glycans to the hinges improved the SAXS fit and resulted in final R factors of 4.8–6.2%. The final ensemble of 113 best-fit models showed that the solution structures of full-length IgA1 possessed extended hinges and asymmetrically positioned Fab and Fc regions. Ample space in IgA1 was revealed for the functionally important binding of two FcαR receptors to its Fc region (Fig. 5a).
The modelling of a large protein–RNA complex was achieved by the combination of SASSIE structural modelling with evidence from chemical footprinting. The hexameric RNA-binding protein Hfq from Escherichia coli (Fig. 1) enables the regulation of mRNA by bacterial small noncoding RNAs (sRNAs) in response to stress and other environmental signals. In order to determine how the Hfq hexamer brings sRNA and mRNA together in the proper orientation for this regulatory function, SASSIE was used to model SAXS data for unbound Hfq (6 × 102 residues) and mRNA (284 nucleotides) and their complex (Peng et al., 2014). The SASSIE modelling was based on the crystal structure of the core hexameric Hfq complexed with two small RNA heptamers (Fig. 1; Wang et al., 2013).
Kratky plots of the SAXS data showed that free mRNA has an extended structure, while Hfq has a compact globular structure. When Hfq was added to mRNA, the change in the appearance of the Kratky plot showed that the extended mRNA structure had compacted and wrapped itself around the Hfq protein (Fig. 1). To verify these scattering results, initial atomistic models were built. That for the full-length E. coli Hfq hexamer was created by appending its disordered N- and C-terminal residues (residues 1–5 and 66–102, respectively) to the crystal structure of the HfQ core. Those for the mRNA models (273 nucleotides) were generated from three-dimensional structures for six RNA fragments using the MC-Sym web server (Parisien & Major, 2008). These six RNA models were merged to give an L-shaped mRNA structure with a flexible pivot between nucleotides 128 and 129. This was energy minimized. The full Hfq–mRNA complex was formed by the superimposition of the initial structures for Hfq and mRNA with the crystal structure of the Hfq core complexed with the two small RNA heptamers. The scattering modelling of the complex was based on variations of both the Hfq and mRNA models. Thus MC simulations were performed in which the terminal residues 1–5 and 66–102 in Hfq were allowed to move, while the Hfq core was held fixed. The simulations held the mRNA structure fixed except for the pivot between nucleotides 128 and 129. From this, 917 models from the 19 132 generated for the complex were accepted to give good scattering fits after comparison with the scattering data. The key result from the 917 models showed that the full-sized mRNA structure could bend around both sides of the Hfq hexamer (Fig. 1). This outcome from the SASSIE modelling provided evidence for how Hfq–mRNA binding could specifically distort mRNA such that sRNA could bind to exposed regions of mRNA, thus explaining the translational control achieved by the sRNAs.
Bottlebrush polymers are a technologically interesting class of macromolecules. As the name suggests, multiple side chain grafts radiate from the polymer backbone, impacting chain flexibility, interactions, self-assembly and dynamics (Zhang et al., 2014). Unlike the biological examples above, for which the starting coordinates were generated from crystal structures, this study used the AMBER MD package (Case et al., 2016) to create a norbornenyl end-functionalized poly(lactide) macromonomer (NB-PLA) with a poly NB (PNB) backbone and PLA side chains. This monomer was replicated 25 times. The resulting polymer was then solvated in tetrahydrofuran, energy minimized in a periodic box and brought to equilibrium. The largest simulated system comprised PNB25–g-PLA19 and 34 298 tetrahydrofuran molecules. Trajectories were output every 100 ps and used to compute the scattering curve, which could then be compared with experimental SANS data. SASSIE was used to automate the processing of many simulation frames and to filter for those conformations whose statistically averaged structures showed better agreement. This analysis demonstrated that structures with an RG value of ∼3.7 nm gave better fits to the SANS data. Moreover, the best-fit trajectories also suggested that the scattering form factor could be well approximated by a short rigid cylinder or ellipsoid of revolution.
SAXS and SANS experiments are powerful experimental methods for elucidating the solution structures of biological macromolecules at low structural resolution. The future development of SAXS and SANS will require the inclusion of advanced atomistic modelling to analyse scattering data properly. There are three main advances offered by CCP-SAS: (i) providing an open-source software environment for developers that (ii) facilitates the easy uptake by users of advanced molecular modelling for the interpretation of SAXS and SANS experiments, and (iii) is seamlessly linked to high-performance computing resources offered through the SASSIE-web front end. SASSIE-web provides a unified workflow framework to molecular simulation engines for both MC and MD, scattering calculators (SasCalc and SCT), and other structure building, job management and general analysis modules (Fig. 2b). The foundation of this platform is the GenApp deployment infrastructure, developed within the CCP-SAS project, which enables the generation of web and standalone GUI applications from the underlying code and provides interfaces to high performance computing resources (Fig. 2a). The SASSIE applications summarized above illustrate how advanced computational modelling will assist scattering projects in addition to the traditional experimental SAXS and SANS data analyses. While daunting at first sight, it is an important challenge for SASSIE to make the atomistic modelling as easy as possible.
The uptake of CCP-SAS software is increasing, with over 200 users registered on the SASSIE-web server to date. At the SAS2015 Conference in Berlin, about 24 protein, nine protein–lipid, eight protein–DNA and two chemical physics projects were reported to be under way. Meanwhile, the number of third-party web applications such as WillItFit is growing, with five groups in various stages of using GenApp and CCP-SAS compute resources to deploy their code. Further enhancements for SASSIE will include additional modules to make the modelling workflow easier for the user. For example a new module, termed PDB-Rx, is in preparation for the Build set of SASSIE modules to help rectify errors identified in PDB-scan or in the user-supplied PDB coordinate file, or to complete any omissions (Wright et al., 2016). The goal of PDB-Rx will be to automate not only the `tidying' and completion of PDB files, but also the preparation of structures using the CHARMM force field for use in the SASSIE simulations. The analysis of new macromolecular systems with intrinsic disorder is becoming increasingly important, following the recognition that many human proteins show disorder. For this, it becomes necessary to develop models that represent ensembles of disordered structures, which is what SASSIE does well (Datta et al., 2007; Green et al., 2016). As illustrated by the MASP and IgA1 examples above, new modules will also be needed to incorporate glycan chains more easily into trial structures, rather than adding the glycans after the best-fit structure is determined. It is a limitation of standard molecular modelling software that the non-protein or non-nucleic acid elements of many systems are not included in many biologically focused packages. And as solution-derived structures become even more commonplace there will also be a need to revisit the deposition of best-fit atomistic structures in their own right in public databases (Wright & Perkins, 2015), together with the experimental data used to derive these structures.
Modules that are specific for soft condensed matter systems important in physical chemistry (Higgins & Benoit, 1996; Gabrys, 2000) have not been described in this article either. However, the development of atomistic models for polymers, surfactant micelles and lipid nanodiscs appropriate for SAS modelling is in progress within CCP-SAS. Besides the usual difficulty in generating a representative starting structure, these systems suffer from a lack of good appropriate force fields. For soft matter systems, it will also be necessary to account for concentration-dependent inter-particle effects, unlike the case of biological systems where scattering data are often extrapolated to infinite dilution. Models of soft matter systems will need to be large enough to allow several micelles to form, and to allow for models showing a realistic degree of polydispersity (not required in biological systems) to be generated. Coarse graining will be essential to achieve molecular models of such systems, particularly when large-scale movements of molecules (rather than just torsional or bending motions within one molecule) are required to generate potential structures for comparison with scattering patterns. This work will also be driven by the desire to model more complex mixed systems, such as surfactant micelles with polymers or colloidal particles, which are the focus of typical standard soft matter small-angle scattering studies.
An exciting prospect going forward is the development of ever more robust, easy to use tools that will eventually enable the SAS user community to routinely take full advantage of combining rapid SAXS and SANS atomistic modelling with data from complementary disciplines such as analytical ultracentrifugation, X-ray crystallography and NMR spectroscopy (§§2.2–2.4), as well as electron microscopy. Indeed the combination of different experimental methods provides new insights not available from one method alone, as demonstrated by the SASSIE applications (Table 1). In addition, the availability of experimentally founded atomistic models allows us to make use of the many programs available in the molecular modelling community for statistical analyses, the evaluation of energetics and the calculation of parameters relevant for other structural techniques such as NMR.
One benefit of being part of the UK's CCP system, which provides an effective means of focusing computational resources for selected communities, is the opportunity to interact regularly with other CCP groups, many of which have overlapping and intersecting interests. The best-known of these is CCP4 (crystallography) (Winn et al., 2011). As well as CCP4, CCP-SAS overlaps with CCP-EM (electron cryomicroscopy) (Wood et al., 2015), CCPN (macromolecular NMR spectroscopy), CCP5 (simulation of soft condensed phases) and CCPBioSim (biomolecular simulation) (Lonsdale et al., 2014). CCP-SAS provides an ideal path forward, not only to tackle the advancement of co-refinement of data from various techniques, but also to advance the soft matter agenda. Although CCP-SAS was initially funded as a joint US/UK venture, the CCP-SAS consortium viewed the project from the outset as more global. As such the project is actively seeking collaborations with and engagement by the global SAS community and welcomes inquiries into creating joint efforts for the benefit of that community.
1This article will form part of a virtual special issue of the journal, presenting some highlights of the 16th International Conference on Small-Angle Scattering (SAS2015).
The authors are supported by the CCP-SAS project, a joint EPSRC (EP/K039121/1) and NSF (CHE-1265821) grant. SJP is supported by the MRC (MR/K011715/1). EHB is supported by NSF (CHE-1265817) and NIH (GM090154). JC is supported by NSF (CHE 1265850). We are most thankful to the many colleagues from the scattering community who have provided bug reports and feedback on SASSIE to date. Certain commercial equipment, instruments, materials, suppliers or software are identified in this paper to foster understanding. Such identification does not imply recommendation or endorsement by the National Institute of Standards and Technology, nor does it imply that the materials or equipment identified are necessarily the best available for the purpose.
Baker, N. A., Sept, D., Joseph, S., Holst, M. J. & McCammon, J. A. (2001). Proc. Natl Acad. Sci. USA, 98, 10037–10041. Web of Science CrossRef PubMed CAS
Best, R. B., Zhu, X., Shim, J., Lopes, P. E. M., Mittal, J., Feig, M. & MacKerell, A. D. Jr (2012). J. Chem. Theory Comput. 8, 3257–3273. Web of Science CrossRef CAS PubMed
Boehm, M. K., Woof, J. M., Kerr, M. A. & Perkins, S. J. (1999). J. Mol. Biol. 286, 1421–1447. Web of Science CrossRef PubMed CAS
Brookes, E. H., Anjum, N., Curtis, J. E., Marru, S., Singh, R. & Pierce, M. (2015). Concurrency Comput. Pract. Exper. 27, 4292–4303. Web of Science CrossRef
Brunger, A. T. (1992). X-PLOR, Version 3.1. A System for X-ray Crystallography and NMR. New Haven: Yale University Press.
Case, D. A. et al. (2016). AMBER 2016. University of California, San Francisco, California, USA.
Castañeda, C. A., Chaturvedi, A., Camara, C. M., Curtis, J. E., Krueger, S. & Fushman, D. (2016). Phys. Chem. Chem. Phys. 18, 5771–5788. Web of Science PubMed
Castañeda, C. A., Dixon, E. K., Walker, O., Chaturvedi, A., Nakasone, M. A., Curtis, J. E., Reed, M. R., Krueger, S., Cropp, T. A. & Fushman, D. (2016). Structure, 24, 1–14. Web of Science PubMed
Chen, J., Brooks, C. L. & Khandogin, J. (2008). Curr. Opin. Struct. Biol. 18, 140–148. Web of Science CrossRef PubMed CAS
Chen, J., Im, W. & Brooks, C. L. (2005). J. Comput. Chem. 26, 1565–1578. Web of Science CrossRef PubMed CAS
Chen, P. C. & Hub, J. S. (2015). Biophys. J. 108, 2573–2584. Web of Science CrossRef CAS PubMed
Clark, N. J., Zhang, H., Krueger, S., Lee, H. J., Ketchem, R. R., Kerwin, B., Kanapuram, S. R., Treuheit, M. J., McAuley, A. & Curtis, J. E. (2013). J. Phys. Chem. B, 117, 14029–14038. Web of Science CrossRef CAS PubMed
Curtis, J. E., Raghunandan, S., Nanda, H. & Krueger, S. (2012). Comput. Phys. Commun. 183, 382–389. Web of Science CrossRef CAS
Datta, S. A. K., Curtis, J. E., Ratcliff, W., Clark, P. K., Crist, R. M., Lebowitz, J., Krueger, S. & Rein, A. (2007). J. Mol. Biol. 365, 812–824. Web of Science CrossRef PubMed CAS
Dewhurst, C. D., Grillo, I., Honecker, D., Bonnaud, M., Jacques, M., Amrouni, C., Perillo-Marcone, A., Manzin, G. & Cubitt, R. (2016). J. Appl. Cryst. 49, 1–14. Web of Science CrossRef CAS IUCr Journals
Dolinsky, T. J., Czodrowski, P., Li, H., Nielsen, J. E., Jensen, J. H., Klebe, G. & Baker, N. A. (2007). Mol. Simul. Nucl. Acids Res. 35, W522–W525. Web of Science CrossRef
Dolinsky, T. J., Nielsen, J. E., McCammon, J. A. & Baker, N. A. (2004). Nucl. Acids Res. 32, W665–W667. Web of Science CrossRef PubMed CAS
Evrard, G., Mareuil, F., Bontems, F., Sizun, C. & Perez, J. (2011). J. Appl. Cryst. 44, 1264–1271. Web of Science CrossRef CAS IUCr Journals
Feigin, L. A. & Svergun, D. I. (1987). Structure Analysis by Small-Angle X-ray and Neutron Scattering. New York: Plenum Press.
Ferrara, P., Apostolakis, J. & Caflisch, A. (2002). Proteins, 46, 24–33. Web of Science CrossRef PubMed CAS
Gabrys, B. J. (2000). Editor. Applications of Neutron Scattering to Soft Condensed Matter. Amsterdam: Gordon and Breach.
Green, M., Hatter, L., Brookes, E., Soultanas, P. & Scott, D. J. (2016). J. Mol. Biol. 428, 357–364. Web of Science CrossRef CAS PubMed
Hanson, R. M., Prilusky, J., Renjian, Z., Nakane, T. & Sussman, J. L. (2013). Isr. J. Chem. 53, 207–216. Web of Science CrossRef CAS
Heenan, R. K., Rogers, S. E., Turner, D., Terry, A. E., Treadgold, J. & King, S. M. (2011). Neutron News, 22(2), 19–21. CrossRef
Herr, A. B., Ballister, E. R. & Bjorkman, P. J. (2003). Nature, 423, 614–620. CrossRef PubMed CAS
Higgins, J. S. & Benoit, H. C. (1996). Polymers and Neutron Scattering. Oxford Series on Neutron Scattering in Condensed Matter, 8. Oxford: Clarendon Press.
Hui, G. K., Wright, D. W., Vennard, O. L., Rayner, L. E., Pang, M., Yeo, S. C., Gor, J., Molyneux, K., Barratt, J. & Perkins, S. J. (2015). Biochem. J. 471, 255–265. Web of Science CrossRef PubMed
Humphrey, W., Dalke, A. & Schulten, K. (1996). J. Mol. Graph. 14, 33–38. Web of Science CrossRef CAS PubMed
Hura, G. L., Menon, A. L., Hammel, M., Rambo, R. P., Poole, F. L., Tsutakawa, S. E., Jenney, F. E. Jr, Classen, S., Frankel, K. A., Hopkins, R. C., Yang, S., Scott, J. W., Dillard, B. D., Adams, M. W. W. & Tainer, J. A. (2009). Nat. Methods, 6, 606–612. Web of Science CrossRef PubMed CAS
Ihms, E. C. & Foster, M. P. (2015). Bioinformatics, 231, 1951–1958. Web of Science CrossRef
Jiménez-García, B., Pons, C., Svergun, D. I., Bernadó, P. & Fernández-Recio, J. (2015). Nucleic Acids Res. 43, W356–W361. Web of Science PubMed
Jo, S., Kim, T., Iyer, V. G. & Im, W. (2008). J. Comput. Chem. 29, 1859–1865. Web of Science CrossRef PubMed CAS
Khan, S., Gor, J., Mulloy, B. & Perkins, S. J. (2010). J. Mol. Biol. 395, 504–521. Web of Science CrossRef PubMed CAS
Knight, C. J. & Hub, J. S. (2015). Nucl. Acids Res. 43, W225–W230. Web of Science CrossRef PubMed
Köfinger, J. & Hummer, G. (2013). Phys. Rev. E, 87, 052712.
Lipfert, J. & Doniach, S. (2007). Annu. Rev. Biophys. Biomol. Struct. 36, 307–327. Web of Science CrossRef PubMed CAS
Lonsdale, R., Rouse, S. L., Sansom, M. S. P. & Mulholland, A. J. (2014). PLoS Comput. Biol. 10, e1003714. Web of Science CrossRef PubMed
MacKerell, A. D. et al. (1998). J. Phys. Chem. B, 102, 3586–3616. CrossRef CAS PubMed
Mayans, M. O., Coadwell, W. J., Beale, D., Symons, D. B. A. & Perkins, S. J. (1995). Biochem. J. 311, 283–291. CrossRef CAS PubMed
Nan, R., Furze, C. M., Wright, D. W., Gor, J., Wallis, R. & Perkins, S. J. (2017). Structure. In the press.
Parisien, D. & Major, F. (2008). Nature, 452, 51–55. Web of Science CrossRef PubMed CAS
Pedersen, M. C., Arleth, L. & Mortensen, K. (2013). J. Appl. Cryst. 46, 1894–1898. Web of Science CrossRef CAS IUCr Journals
Pelikan, M., Hura, G. L. & Hammel, M. (2009). Gen. Physiol. Biophys. 28, 174–189. Web of Science CrossRef PubMed CAS
Peng, Y., Curtis, J. E., Fang, X. & Woodson, S. A. (2014). Proc. Natl Acad. Sci. USA, 111, 17134–17139. Web of Science CrossRef CAS PubMed
Perkins, S. J. (2001). Biophys. Chem. 93, 129–139. Web of Science CrossRef PubMed CAS
Perkins, S. J., Nan, R., Li, K., Khan, S. & Abe, Y. (2011). Methods, 54, 181–199. Web of Science CrossRef CAS PubMed
Perkins, S. J., Nealis, A. S., Sutton, B. J. & Feinstein, A. (1991). J. Mol. Biol. 221, 1345–1366. CrossRef PubMed CAS
Perkins, S. J., Okemefuna, A. I., Fernando, A. N., Bonner, A., Gilbert, H. E. & Furtado, P. B. (2008). Methods Cell Biol. 84, 375–423. Web of Science CrossRef PubMed CAS
Perkins, S. J., Okemefuna, A. I., Nan, R., Li, K. & Bonner, A. (2009). J. R. Soc. Interface, 6, S679–S696. Web of Science CrossRef PubMed CAS
Pernot, P. et al. (2013). J. Synchrotron Rad. 20, 660–664. Web of Science CrossRef CAS IUCr Journals
Petoukhov, M. V., Franke, D., Shkumatov, A. V., Tria, G., Kikhney, A. G., Gajda, M., Gorba, C., Mertens, H. D. T., Konarev, P. V. & Svergun, D. I. (2012). J. Appl. Cryst. 45, 342–350. Web of Science CrossRef CAS IUCr Journals
Pettersen, E. F., Goddard, T. D., Huang, C. C., Couch, G. S., Greenblatt, D. M., Meng, E. C. & Ferrin, T. E. (2004). J. Comput. Chem. 25, 1605–1612. Web of Science CrossRef PubMed CAS
Phillips, J. C., Braun, R., Wang, W., Gumbart, J., Tajkhorshid, E., Villa, E., Chipot, C., Skeel, R. D., Kalé, L. & Schulten, K. (2005). J. Comput. Chem. 26, 1781–1802. Web of Science CrossRef PubMed CAS
Poitevin, F., Orland, H., Doniach, S., Koehl, P. & Delarue, M. (2011). Nucleic Acids Res. 39, W184–W189. Web of Science CrossRef CAS PubMed
Putnam, C. D., Hammel, M., Hura, G. L. & Tainer, J. A. (2007). Q. Rev. Biophys. 40, 191–285. Web of Science CrossRef PubMed CAS
Rambo, R. P. & Tainer, J. A. (2013). Annu. Rev. Biophys. 42, 415–441. Web of Science CrossRef CAS PubMed
Rayner, L. E., Hui, G. K., Gor, J., Heenan, R. K., Dalby, P. A. & Perkins, S. J. (2014). J. Biol. Chem. 289, 20740–20756. Web of Science CrossRef CAS PubMed
Rayner, L. E., Hui, G. K., Gor, J., Heenan, R. K., Dalby, P. A. & Perkins, S. J. (2015). J. Biol. Chem. 290, 8420–8438. Web of Science CrossRef CAS PubMed
Round, A., Felisaz, F., Fodinger, L., Gobbo, A., Huet, J., Villard, C., Blanchet, C. E., Pernot, P., McSweeney, S., Roessle, M., Svergun, D. I. & Cipriani, F. (2015). Acta Cryst. D71, 67–75. Web of Science CrossRef IUCr Journals
Różycki, B., Kim, Y. C. & Hummer, G. (2011). Structure, 19, 109–116. Web of Science PubMed
Šali, A. & Blundell, T. L. (1993). J. Mol. Biol. 234, 779–815. PubMed Web of Science
Schneidman-Duhovny, D., Hammel, M. & Sali, A. (2010). Nucleic Acids Res. 38, W540–W544. Web of Science CAS PubMed
Schneidman-Duhovny, D., Hammel, M., Tainer, J. A. & Sali, A. (2016). Nucleic Acids Res. 44, W424–W429. Web of Science PubMed
Spinozzi, F. & Beltramini, M. (2012). Biophys. J. 103, 511–521. Web of Science CrossRef CAS PubMed
Svergun, D., Barberato, C. & Koch, M. H. J. (1995). J. Appl. Cryst. 28, 768–773. CrossRef CAS Web of Science IUCr Journals
Svergun, D. I., Koch, M. H. J., Timmins, P. A. & May, R. P. (2013). Small Angle X-ray and Neutron Scattering from Solutions of Biological Macromolecules. IUCr Texts on Crystallography, 19. Oxford University Press.
Wang, W., Wang, L., Wu, J., Gong, Q. & Shi, Y. (2013). Nucleic Acids Res. 41, 5938–5948. Web of Science CrossRef CAS PubMed
Watson, M. C. & Curtis, J. E. (2013). J. Appl. Cryst. 46, 1171–1177. Web of Science CrossRef CAS IUCr Journals
Whitten, A. E., Cai, S. & Trewhella, J. (2008). J. Appl. Cryst. 41, 222–226. Web of Science CrossRef CAS IUCr Journals
Winn, M. D. et al. (2011). Acta Cryst. D67, 235–242. Web of Science CrossRef CAS IUCr Journals
Wood, C., Burnley, T., Patwardhan, A., Scheres, S., Topf, M., Roseman, A. & Winn, M. (2015). Acta Cryst. D71, 123–126. Web of Science CrossRef IUCr Journals
Wright, D. W. & Perkins, S. J. (2015). J. Appl. Cryst. 48, 953–961. Web of Science CrossRef CAS IUCr Journals
Wright, D. W., Zhang, H., Perkins, S. J. & Curtis, J. E. (2016). In preparation.
Yang, S., Blachowicz, L., Makowski, L. & Roux, B. (2010). Proc. Natl Acad. Sci. USA, 107, 15757–15762. Web of Science CrossRef CAS PubMed
Yang, S., Park, S., Makowski, L. & Roux, B. (2009). Biophys. J. 96, 4449–4463. Web of Science CrossRef PubMed CAS
Zhang, W., Chen, J., Howell, S. C., Heindel, A., Wright, D. W., Perkins, S. J. & Curtis, J. E. (2016). In preparation.
Zhang, Z., Carrillo, J.-M. Y., Ahn, S., Wu, B., Hong, K., Smith, G. S. & Do, C. (2014). Macromolecules, 47, 5808–5814. Web of Science CrossRef CAS
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.