Characterization of Protein Flexibility Using Small-Angle X-Ray Scattering and Amplified Collective Motion Simulations

Large-scale flexibility within a multidomain protein often plays an important role in its biological function. Despite its inherent low resolution, small-angle x-ray scattering (SAXS) is well suited to investigate protein flexibility and determine, with the help of computational modeling, what kinds of protein conformations would coexist in solution. In this article, we develop a tool that combines SAXS data with a previously developed sampling technique called amplified collective motions (ACM) to elucidate structures of highly dynamic multidomain proteins in solution. We demonstrate the use of this tool in two proteins, bacteriophage T4 lysozyme and tandem WW domains of the formin-binding protein 21. The ACM simulations can sample the conformational space of proteins much more extensively than standard molecular dynamics (MD) simulations. Therefore, conformations generated by ACM are significantly better at reproducing the SAXS data than are those from MD simulations.     

Structural analysis of intrinsically disordered proteins by small-angle X-ray scattering

Small-angle scattering of X-rays (SAXS) is an established method to study the overall structure and structural transitions of biological macromolecules in solution. For folded proteins, the technique provides three-dimensional low resolution structures ab initio or it can be used to drive rigid-body modeling. SAXS is also a powerful tool for the quantitative analysis of flexible systems, including intrinsically disordered proteins (IDPs), and is highly complementary to the high resolution methods of X-ray crystallography and NMR. Here we present the basic principles of SAXS and review the main approaches to the characterization of IDPs and flexible multidomain proteins using SAXS. Together with the standard approaches based on the analysis of overall parameters, a recently developed Ensemble Optimization Method (EOM) is now available. The latter method allows for the co-existence of multiple protein conformations in solution compatible with the scattering data. Analysis of the selected ensembles provides quantitative information about flexibility and also offers insights into structural features. Examples of the use of SAXS and combined approaches with NMR, X-ray crystallography, and computational methods to characterize completely or partially disordered proteins are presented.     

Conformational dynamics of a multidomain protein by neutron scattering and computational analysis

The flexible conformations of a multidomain protein are responsible for its biological functions. Although MurD, a 47-kDa protein that consists of three domains, sequentially changes its domain conformation from an open form to a closed form through a semiclosed form in its enzymatic reaction, the domain dynamics in each conformation remains unclear. In this study, we verify the conformational dynamics of MurD in the corresponding three states (apo and ATP- and inhibitor-bound states) with a combination of small-angle x-ray and neutron scattering (SAXS and SANS), dynamic light scattering (DLS), neutron backscattering (NBS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations. Applying principal component analysis of the MD trajectories, twisting and open-closed domain modes are identified as the major collective coordinates. The deviations of the experimental SAXS profiles from the theoretical calculations based on the known crystal structures become smaller in the ATP-bound state than in the apo state, and a further decrease is evident upon inhibitor binding. These results suggest that domain motions of the protein are suppressed step by step of each ligand binding. The DLS and NBS data yield collective and self-translational diffusion constants, respectively, and we used them to extract collective domain motions in nanometer and nanosecond scales from the NSE data. In the apo state, MurD shows both twisting and open-closed domain modes, whereas an ATP binding suppresses twisting domain motions, and a further reduction of open-closed mode is seen in the inhibitor-binding state. These observations are consistent with the structure modifications measured by the small-angle scattering as well as the MD simulations. Such changes in the domain dynamics associated with the sequential enzymatic reactions should be related to the affinity and reaction efficiency with a ligand that binds specifically to each reaction state.     

SAXSDom: Modeling multidomain protein structures using small-angle X-ray scattering data

Many proteins are composed of several domains that pack together into a complex tertiary structure. Multidomain proteins can be challenging for protein structure modeling, particularly those for which templates can be found for individual domains but not for the entire sequence. In such cases, homology modeling can generate high quality models of the domains but not for the orientations between domains. Small-angle X-ray scattering (SAXS) reports the structural properties of entire proteins and has the potential for guiding homology modeling of multidomain proteins. In this article, we describe a novel multidomain protein assembly modeling method, SAXSDom that integrates experimental knowledge from SAXS with probabilistic Input-Output Hidden Markov model to assemble the structures of individual domains together. Four SAXS-based scoring functions were developed and tested, and the method was evaluated on multidomain proteins from two public datasets. Incorporation of SAXS information improved the accuracy of domain assembly for 40 out of 46 critical assessment of protein structure prediction multidomain protein targets and 45 out of 73 multidomain protein targets from the ab initio domain assembly dataset. The results demonstrate that SAXS data can provide useful information to improve the accuracy of domain-domain assembly. The source code and tool packages are available at https://github.com/jianlin-cheng/SAXSDom .     

Effect of interdomain dynamics on the structure determination of modular proteins by small-angle scattering

Multidomain proteins in which consecutive globular regions are connected by linkers are prevalent in nature (Levitt in Proc Natl Acad Sci USA 106:11079–11084, 2009). Some members of this family have largely resisted structural characterization as a result of challenges associated with their inherent flexibility. Small-angle scattering (SAS) is often the method of choice for their structural study. An extensive set of simulated data for both flexible and rigid multidomain systems was analyzed and modeled using standard protocols. This study clearly shows that SAXS profiles obtained from highly flexible proteins can be wrongly interpreted as arising from a rigid structure. In this context, it would be important to identify features from the SAXS data or from the derived structural models that indicate interdomain motions to differentiate between these two scenarios. Features of SAXS data that identify flexible proteins are: (1) general attenuation of fine structure in the scattering profiles, which becomes more dramatic in Kratky representations, and (2) a reduced number of interdomain correlation peaks in p(r) functions that also present large D _max values and a smooth decrease to 0. When modeling this dynamically averaged SAXS data, the structures obtained present characteristic trends: (1) ab initio models display a decrease in resolution, and (2) rigid-body models present highly extended conformations with a lack of interdomain contacts. The ensemble optimization method represents an excellent strategy to identify interdomain motions unambiguously. This study provides information that should help researchers to select the best modeling strategy for the structural interpretation of SAS experiments of multidomain proteins.     

Refining conformational ensembles of flexible proteins against small-angle x-ray scattering data

Intrinsically disordered proteins and flexible regions in multidomain proteins display substantial conformational heterogeneity. Characterizing the conformational ensembles of these proteins in solution typically requires combining one or more biophysical techniques with computational modeling or simulations. Experimental data can either be used to assess the accuracy of a computational model or to refine the computational model to get a better agreement with the experimental data. In both cases, one generally needs a so-called forward model (i.e., an algorithm to calculate experimental observables from individual conformations or ensembles). In many cases, this involves one or more parameters that need to be set, and it is not always trivial to determine the optimal values or to understand the impact on the choice of parameters. For example, in the case of small-angle x-ray scattering (SAXS) experiments, many forward models include parameters that describe the contribution of the hydration layer and displaced solvent to the background-subtracted experimental data. Often, one also needs to fit a scale factor and a constant background for the SAXS data but across the entire ensemble. Here, we present a protocol to dissect the effect of the free parameters on the calculated SAXS intensities and to identify a reliable set of values. We have implemented this procedure in our Bayesian/maximum entropy framework for ensemble refinement and demonstrate the results on four intrinsically disordered proteins and a protein with three domains connected by flexible linkers. Our results show that the resulting ensembles can depend on the parameters used for solvent effects and suggest that these should be chosen carefully. We also find a set of parameters that work robustly across all proteins.     

Structural Characterization of N-WASP Domain V Using MD Simulations with NMR and SAXS Data

Because of their large conformational heterogeneity, structural characterization of intrinsically disordered proteins (IDPs) is very challenging using classical experimental methods alone. In this study, we use NMR and small-angle x-ray scattering (SAXS) data with multiple molecular dynamics (MD) simulations to describe the conformational ensemble of the fully disordered verprolin homology domain of the neural Aldrich syndrome protein involved in the regulation of actin polymerization. First, we studied several back-calculation software of SAXS scattering intensity and optimized the adjustable parameters to accurately calculate the SAXS intensity from an atomic structure. We also identified the most appropriate force fields for MD simulations of this IDP. Then, we analyzed four conformational ensembles of neural Aldrich syndrome protein verprolin homology domain, two generated with the program flexible-meccano with or without NMR-derived information as input and two others generated by MD simulations with two different force fields. These four conformational ensembles were compared to available NMR and SAXS data for validation. We found that MD simulations with the AMBER-03w force field and the TIP4P/2005s water model are able to correctly describe the conformational ensemble of this 67-residue IDP at both local and global level.     

Validating Solution Ensembles from Molecular Dynamics Simulation by Wide-Angle X-ray Scattering Data

Wide-angle x-ray scattering (WAXS) experiments of biomolecules in solution have become increasingly popular because of technical advances in light sources and detectors. However, the structural interpretation of WAXS profiles is problematic, partly because accurate calculations of WAXS profiles from structural models have remained challenging. In this work, we present the calculation of WAXS profiles from explicit-solvent molecular dynamics (MD) simulations of five different proteins. Using only a single fitting parameter that accounts for experimental uncertainties because of the buffer subtraction and dark currents, we find excellent agreement to experimental profiles both at small and wide angles. Because explicit solvation eliminates free parameters associated with the solvation layer or the excluded solvent, which would require fitting to experimental data, we minimize the risk of overfitting. We further find that the influence from water models and protein force fields on calculated profiles are insignificant up to q ≈ 15 nm^?1. Using a series of simulations that allow increasing flexibility of the proteins, we show that incorporating thermal fluctuations into the calculations significantly improves agreement with experimental data, demonstrating the importance of protein dynamics in the interpretation of WAXS profiles. In addition, free MD simulations up to one microsecond suggest that the calculated profiles are highly sensitive with respect to minor conformational rearrangements of proteins, such as an increased flexibility of a loop or an increase of the radius of gyration by  <  1%. The present study suggests that quantitative comparison between MD simulations and experimental WAXS profiles emerges as an accurate tool to validate solution ensembles of biomolecules.     

Molecular insights on CALX-CBD12 interdomain dynamics from MD simulations,RDCs, and SAXS

Na⁺/Ca²⁺ exchangers (NCXs) are secondary active transporters that couple the translocation of Na⁺ with the transport of Ca²⁺ in the opposite direction. The exchanger is an essential Ca²⁺ extrusion mechanism in excitable cells. It consists of a transmembrane domain and a large intracellular loop that contains two Ca²⁺-binding domains, CBD1 and CBD2. The two CBDs are adjacent to each other and form a two-domain Ca²⁺ sensor called CBD12. Binding of intracellular Ca²⁺ to CBD12 activates the NCX but inhibits the NCX of Drosophila, CALX. NMR spectroscopy and SAXS studies showed that CALX and NCX CBD12 constructs display significant interdomain flexibility in the apo state but assume rigid interdomain arrangements in the Ca²⁺-bound state. However, detailed structure information on CBD12 in the apo state is missing. Structural characterization of proteins formed by two or more domains connected by flexible linkers is notoriously challenging and requires the combination of orthogonal information from multiple sources. As an attempt to characterize the conformational ensemble of CALX-CBD12 in the apo state, we applied molecular dynamics (MD) simulations, NMR (¹H-¹⁵N residual dipolar couplings), and small-angle x-ray scattering (SAXS) data in a combined strategy to select an ensemble of conformations in agreement with the experimental data. This joint approach demonstrated that CALX-CBD12 preferentially samples closed conformations, whereas the wide-open interdomain arrangement characteristic of the Ca²⁺-bound state is less frequently sampled. These results are consistent with the view that Ca²⁺ binding shifts the CBD12 conformational ensemble toward extended conformers, which could be a key step in the NCXs’ allosteric regulation mechanism. This strategy, combining MD with NMR and SAXS, provides a powerful approach to select ensembles of conformations that could be applied to other flexible multidomain systems.     

Entropic folding pathway of human epidermal growth factor explored by disulfide scrambling and amplified collective motion simulations

Epidermal growth factor (EGF) regulates cell proliferation and differentiation by binding to the EGF receptor (EGFR) extra-cellular domains. Human EGF is a small, single-chain protein comprising three distinct loops (A, B, and C), which are connected by three disulfide bridges (Cys6-Cys20, Cys14-Cys31, and Cys33-Cys42). These disulfide bridges are essential for structural stability and biological activity. EGF was extensively studied by disulfide scrambling, an experimental technique for the conformational entrapment of intermediate states, which allows us to study the folding pathway of proteins containing disulfide bonds. The experimental results showed that there is a major 2-disulfide intermediate (denoted EGF-II) and that the native disulfide bonding pattern is less prevalent in one of the mutants. In this article, we investigated for the first time the solution conformations of wild-type EGF, EGF-II, and the mutant S9C through extensive molecular dynamics (MD) simulations in water using both the standard MD technique and a recently developed amplified-collective-motion (ACM) sampling method. Compared to standard MD simulations, we achieved a much more enhanced sampling by the ACM simulations, and the structures were sufficiently relaxed to estimate configurational entropies. The simulation results suggest a predominantly entropic folding pathway governed by the disorder of three functional loop regions. Although EGF-II exhibits two native disulfide bonds (Cys14-Cys31 and Cys33- Cys42), its large configurational entropy inhibits a direct transition to the native structure in the folding process. When Ser9 is mutated into Cys, a non-native disulfide bridge Cys9- Cys20 is slightly more favorable than the native Cys6-Cys20 because a less constrained N-terminus affords larger entropy. Isomers that are functionally less active also exhibit a more localized dynamics of the functional loop regions, which may suggest a possible mechanism for the modulation of EGF activity.     

Structural analysis of flexible proteins in solution by small angle X-ray scattering combined with crystallography

In the last few years, SAXS of biological materials has been rapidly evolving and promises to move structural analysis to a new level. Recent innovations in SAXS data analysis allow ab initio shape predictions of proteins in solution. Furthermore, experimental scattering data can be compared to calculated scattering curves from the growing data base of solved structures and also identify aggregation and unfolded proteins. Combining SAXS results with atomic resolution structures enables detailed characterizations in solution of mass, radius, conformations, assembly, and shape changes associated with protein folding and functions. SAXS can efficiently reveal the spatial organization of protein domains, including domains missing from or disordered in known crystal structures, and establish cofactor or substrate-induced conformational changes. For flexible domains or unstructured regions that are not amenable for study by many other structural techniques, SAXS provides a unique technology. Here, we present SAXS shape predictions for PCNA that accurately predict a trimeric ring assembly and for a full-length DNA repair glycosylase with a large unstructured region. These new results in combination with illustrative published data show how SAXS combined with high resolution crystal structures efficiently establishes architectures, assemblies, conformations, and unstructured regions for proteins and protein complexes in solution.     

Dynamic Motion and Communication in the Streptococcal C1 Phage Lysin,PlyC

The growing problem of antibiotic resistance underlies the critical need to develop new treatments to prevent and control resistant bacterial infection. Exogenous application of bacteriophage lysins results in rapid and specific destruction of Gram-positive bacteria and therefore lysins represent novel antibacterial agents. The PlyC phage lysin is the most potent lysin characterized to date and can rapidly lyse Group A, C and E streptococci. Previously, we have determined the X-ray crystal structure of PlyC, revealing a complicated and unique arrangement of nine proteins. The scaffold features a multimeric cell-wall docking assembly bound to two catalytic domains that communicate and work synergistically. However, the crystal structure appeared to be auto-inhibited and raised important questions as to the mechanism underlying its extreme potency. Here we use small angle X-ray scattering (SAXS) and reveal that the conformational ensemble of PlyC in solution is different to that in the crystal structure. We also investigated the flexibility of the enzyme using both normal mode (NM) analysis and molecular dynamics (MD) simulations. Consistent with our SAXS data, MD simulations show rotational dynamics of both catalytic domains, and implicate inter-domain communication in achieving a substrate-ready conformation required for enzyme function. Our studies therefore provide insights into how the domains in the PlyC holoenzyme may act together to achieve its extraordinary potency.     

Backbone conformational flexibility of the lipid modified membrane anchor of the human N-Ras protein investigated by solid-state NMR and molecular dynamics simulation

The lipid modified human N-Ras protein, implicated in human cancer development, is of particular interest due to its membrane anchor that determines the activity and subcellular location of the protein. Previous solid-state NMR investigations indicated that this membrane anchor is highly dynamic, which may be indicative of backbone conformational flexibility. This article aims to address if a dynamic exchange between three structural models exist that had been determined previously. We applied a combination of solid-state nuclear magnetic resonance (NMR) methods and replica exchange molecular dynamics (MD) simulations using a Ras peptide that represents the terminal seven amino acids of the human N-Ras protein. Analysis of correlations between the conformations of individual amino acids revealed that Cys 181 and Met 182 undergo collective conformational exchange. Two major structures constituting about 60% of all conformations could be identified. The two conformations found in the simulation are in rapid exchange, which gives rise to low backbone order parameters and nuclear spin relaxation as measured by experimental NMR methods. These parameters were also determined from two 300 ns conventional MD simulations, providing very good agreement with the experimental data.     

Solvation Effects on the Conformational Behaviour of Gellan and Calcium Ion Binding to Gellan Double Helices

The contribution of the presence of solvent to the conformations adopted by disaccharide fragments within the repeat unit of gellan have been studied by molecular modelling techniques. Initial conformational energy searches, using a dielectric continuum to represent the solvent, provided starting geometries for a series of molecular dynamics simulations. The solution behaviour from these simulations was subsequently compared to fibre diffraction data of the potassium gellan salt. The present calculations indicate considerable flexibility of the glycosidic linkages, and this is discussed in relation to its effect on gel formation. One of the fragments was solvated with explicit water molecules. These calculations showed the same conformational behaviour as those simulations conducted in implicit solvent.Finally, a series of molecular dynamics (MD) simulations were performed to study the calcium binding to gellan. The results from this clearly showed a well defined binding site for this ion.Abbreviations MM molecular mechanics - MD molecular dynamics     

A review of multi-domain and flexible molecular chaperones studies by small-angle X-ray scattering

Intrinsic flexibility is closely related to protein function, and a plethora of important regulatory proteins have been found to be flexible, multi-domain or even intrinsically disordered. On the one hand, understanding such systems depends on how these proteins behave in solution. On the other, small-angle X-ray scattering (SAXS) is a technique that fulfills the requirements to study protein structure and dynamics relatively quickly with few experimental limitations. Molecular chaperones from Hsp70 and Hsp90 families are multi-domain proteins containing flexible and/or disordered regions that play central roles in cellular proteostasis. Here, we review the structure and function of these proteins by SAXS. Our general approach includes the use of SAXS data to determine size and shape parameters, as well as protein shape reconstruction and their validation by using accessory biophysical tools. Some remarkable examples are presented that exemplify the potential of the SAXS technique. Protein structure can be determined in solution even at limiting protein concentrations (for example, human mortalin, a mitochondrial Hsp70 chaperone). The protein organization, flexibility and function (for example, the J-protein co-chaperones), oligomeric status, domain organization, and flexibility (for the Hsp90 chaperone and the Hip and Hep1 co-chaperones) may also be determined. Lastly, the shape, structural conservation, and protein dynamics (for the Hsp90 chaperone and both p23 and Aha1 co-chaperones) may be studied by SAXS. We believe this review will enhance the application of the SAXS technique to the study of the molecular chaperones.     

Determination of conformational equilibrium of peptides in solution by NMR spectroscopy and theoretical conformational analysis: application to the calibration of mean-field solvation models.

Peptides occur in solution as ensembles of conformations rather than in a fixed conformation. The existing energy functions are usually inadequate to predict the conformational equilibrium in solution, because of failure to account properly for solvation, if the solvent is not considered explicitly (which is usually prohibitively expensive). NMR data are therefore widely incorporated into theoretical conformational analysis. Because of conformational flexibility, restrained molecular dynamics (with restraints derived from NMR data), which is usually applied to determine protein conformation is of limited use in the case of peptides. Instead, (a) the restraints are averaged within predefined time windows during molecular dynamics (MD) simulations (time averaging), (b) multiple-copy MD simulations are carried out and the restraints are averaged over the copies (ensemble averaging), or (c) a representative ensemble of sterically feasible conformations is generated and the weights of the conformations are then fitted so that the computed average observables match the experimental data (weight fitting). All these approaches are briefly discussed in this article. If an adequate force field is used, conformations with large statistical weights obtained from the weight-fitting procedure should also have low energies, which can be implemented in force field calibration. Such a procedure is particularly attractive regarding the parameterization of the solvation energy in nonaqueous solvents, e.g., dimethyl sulfoxide, for which thermodynamic solvation data are scarce. A method for calibration of solvation parameters in dimethyl sulfoxide, which is based on this principle was recently proposed by C. Baysal and H. Meirovitch (Journal of the American Chemical Society, 1998, Vol. 120, pp. 800--812), in which the energy gap between the conformations compatible with NMR data and the alternative conformations is maximized. In this work we propose an alternative method based on the principle that the best-fitting statistical weights of conformations should match the Boltzmann weights computed with the force field applied. Preliminary results obtained using three test peptides of varying conformational mobility: H-Ser(1)-Pro(2)-Lys(3)-Leu(4)-OH, Ac-Tyr(1)-D-Phe(2)-Ser(3)-Pro(4)-Lys(5)-Leu(6)-NH(2), and cyclo(Tyr(1)-D-Phe(2)-Ser(3)-Pro(4)-Lys(5)-Leu(6)) are presented.     

Rotamer Dynamics: Analysis of Rotamers in Molecular Dynamics Simulations of Proteins

Given by χ torsional angles, rotamers describe the side-chain conformations of amino acid residues in a protein based on the rotational isomers (hence the word rotamer). Constructed rotamer libraries, based on either protein crystal structures or dynamics studies, are the tools for classifying rotamers (torsional angles) in a way that reflect their frequency in nature. Rotamer libraries are routinely used in structure modeling and evaluation. In this perspective article, we would like to encourage researchers to apply rotamer analyses beyond their traditional use. Molecular dynamics (MD) of proteins highlight the in silico behavior of molecules in solution and thus can identify favorable side-chain conformations. In this article, we used simple computational tools to study rotamer dynamics (RD) in MD simulations. First, we isolated each frame in the MD trajectories in separate Protein Data Bank files via the cpptraj module in AMBER. Then, we extracted torsional angles via the Bio3D module in R language. The classification of torsional angles was also done in R according to the penultimate rotamer library. RD analysis is useful for various applications such as protein folding, study of rotamer-rotamer relationship in protein-protein interaction, real-time correlation between secondary structures and rotamers, study of flexibility of side chains in binding site for molecular docking preparations, use of RD as guide in functional analysis and study of structural changes caused by mutations, providing parameters for improving coarse-grained MD accuracy and speed, and many others. Major challenges facing RD to emerge as a new scientific field involve the validation of results via easy, inexpensive wet-lab methods. This realm is yet to be explored.     

Integration of small-angle X-ray scattering data into structural modeling of proteins and their assemblies

A major challenge in structural biology is to determine the configuration of domains and proteins in multidomain proteins and assemblies, respectively. All available data should be considered to maximize the accuracy and precision of these models. Small-angle X-ray scattering (SAXS) efficiently provides low-resolution experimental data about the shapes of proteins and their assemblies. Thus, we integrated SAXS profiles into our software for modeling proteins and their assemblies by satisfaction of spatial restraints. Specifically, we modeled the quaternary structures of multidomain proteins with structurally defined rigid domains as well as quaternary structures of binary complexes of structurally defined rigid proteins. In addition to SAXS profiles and the component structures, we used stereochemical restraints and an atomic distance-dependent statistical potential. The scoring function is optimized by a biased Monte Carlo protocol, including quasi-Newton and simulated annealing schemes. The final prediction corresponds to the best scoring solution in the largest cluster of many independently calculated solutions. To quantify how well the quaternary structures are determined based on their SAXS profiles, we used a benchmark of 12 simulated examples as well as an experimental SAXS profile of the homotetramer d-xylose isomerase. Optimization of the SAXS-dependent scoring function generally results in accurate models if sufficiently precise approximations for the constituent rigid bodies are available; otherwise, the best scoring models can have significant errors. Thus, SAXS profiles can play a useful role in the structural characterization of proteins and assemblies if they are combined with additional data and used judiciously. Our integration of a SAXS profile into modeling by satisfaction of spatial restraints will facilitate further integration of different kinds of data for structure determination of proteins and their assemblies.