Validation of macromolecular flexibility in solution by small-angle X-ray scattering (SAXS)

Isoleucine binding and regulation of Escherichia coli and Staphylococcus aureus threonine dehydratase (IlvA)

In Staphylococcus aureus, the branched-chain amino acid biosynthetic pathway provides essential intermediates for membrane biosynthesis. Threonine deaminase (IlvA) is the first enzyme in the pathway, and isoleucine feedback-regulates the enzyme in Escherichia coli. These studies on E. coli IlvA (EcIlvA) introduced the concept of allosteric regulation. To investigate the regulation of S. aureus IlvA (SaIlvA), we first conducted additional studies on EcIlvA. The previously determined crystal structure of EcIlvA revealed a tetrameric assembly of protomers, each with catalytic and regulatory domains, but the structural basis of isoleucine regulation was not characterized. Here, we present the crystal structure of the EcIlvA regulatory domain bound to isoleucine, which reveals the isoleucine binding site and conformational changes that initiate at Phe352 and propagate 23 Angstrom across the domain. This suggests an allosteric pathway that extends to the active site of the adjacent protomer, mediating regulation across the protomer-protomer interface. The EcIlvA(F352A) mutant binds isoleucine but is feedback-resistant due to the absence of the initiating Phe352. In contrast, SaIlvA is not feedback-regulated by isoleucine and does not bind it. The structure of the SaIlvA regulatory domain reveals a different organization that lacks the isoleucine binding site. Other potential allosteric inhibitors of SaIlvA, including phospholipid intermediates, do not affect enzyme activity. We propose that the absence of feedback inhibition in SaIlvA is due to its role in membrane biosynthesis. These findings enhance our understanding of IlvA's allosteric regulation and offer opportunities for engineering feedback-resistant IlvA variants for biotechnological use.

Predicting RNA structure and dynamics with deep learning and solution scattering

Advanced deep learning and statistical methods can predict structural models for RNA molecules. However, RNAs are flexible, and it remains difficult to describe their macromolecular conformations in solutions where varying conditions can induce conformational changes. Small-angle x-ray scattering (SAXS) in solution is an efficient technique to validate structural predictions by comparing the experimental SAXS profile with those calculated from predicted structures. There are two main challenges in comparing SAXS profiles to RNA structures: the absence of cations essential for stability and charge neutralization in predicted structures and the inadequacy of a single structure to represent RNA's conformational plasticity. We introduce a solution conformation predictor for RNA (SCOPER) to address these challenges. This pipeline integrates kinematics-based conformational sampling with the innovative deep learning model, IonNet, designed for predicting Mg2+ ion binding sites. Validated through benchmarking against 14 experimental data sets, SCOPER significantly improved the quality of SAXS profile fits by including Mg2+ ions and sampling of conformational plasticity. We observe that an increased content of monovalent and bivalent ions leads to decreased RNA plasticity. Therefore, carefully adjusting the plasticity and ion density is crucial to avoid overfitting experimental SAXS data. SCOPER is an efficient tool for accurately validating the solution state of RNAs given an initial, sufficiently accurate structure and provides the corrected atomistic model, including ions.

Predicting RNA Structure and Dynamics with Deep Learning and Solution Scattering

Advanced deep learning and statistical methods can predict structural models for RNA molecules. However, RNAs are flexible, and it remains difficult to describe their macromolecular conformations in solutions where varying conditions can induce conformational changes. Small-angle X-ray scattering (SAXS) in solution is an efficient technique to validate structural predictions by comparing the experimental SAXS profile with those calculated from predicted structures. There are two main challenges in comparing SAXS profiles to RNA structures: the absence of cations essential for stability and charge neutralization in predicted structures and the inadequacy of a single structure to represent RNA's conformational plasticity. We introduce Solution Conformation Predictor for RNA (SCOPER) to address these challenges. This pipeline integrates kinematics-based conformational sampling with the innovative deep-learning model, IonNet, designed for predicting Mg2+ ion binding sites. Validated through benchmarking against fourteen experimental datasets, SCOPER significantly improved the quality of SAXS profile fits by including Mg2+ ions and sampling of conformational plasticity. We observe that an increased content of monovalent and bivalent ions leads to decreased RNA plasticity. Therefore, carefully adjusting the plasticity and ion density is crucial to avoid overfitting experimental SAXS data. SCOPER is an efficient tool for accurately validating the solution state of RNAs given an initial, sufficiently accurate structure and provides the corrected atomistic model, including ions.

Structural and Functional Biology of Mammalian ALOX Isoforms with Particular Emphasis on Enzyme Dimerization and Their Allosteric Properties

The human genome involves six functional arachidonic acid (AA) lipoxygenase (ALOX) genes, and the corresponding enzymes (ALOX15, ALOX15B, ALOX12, ALOX12B, ALOXE3, ALOX5) have been implicated in cell differentiations and in the pathogenesis of inflammatory, hyperproliferative, metabolic, and neurological disorders. Humans express two different AA 15-lipoxygenating ALOX isoforms, and these enzymes are called ALOX15 (15-LOX1) and ALOX15B (15-LOX2). Chromosomal localization, sequence alignments, and comparison of the enzyme properties suggest that pig and mouse ALOX15 orthologs (leukocyte-type 12-LOX) on the one hand and rabbit and human ALOX15 orthologs on the other (reticulocyte-type 15-LOX1) belong to the same enzyme family despite their different reaction specificities with AA as a substrate. In contrast, human ALOX12 (platelet-type 12-LOX), as well as pig and mouse ALOX15 (leukocyte-type 12-LOX), belong to different enzyme families, although they exhibit a similar reaction specificity with AA as a substrate. The complex multiplicity of mammalian ALOX isoforms and the controversial enzyme nomenclatures are highly confusing and prompted us to summarize the current knowledge on the biological functions, enzymatic properties, and allosteric regulation mechanisms of mammalian ALOX15, ALOX15B, and ALOX12 orthologs that belong to three different enzyme sub-families.

Perspectives on solution-based small angle X-ray scattering for protein and biological macromolecule structural biology

Small-angle X-ray scattering (SAXS) is used to extract structural information from a wide variety of non-crystalline samples in different fields (e.g., materials science, physics, chemistry, and biology). This review provides an overview of SAXS as applied to structural biology, specifically for proteins and other biomacromolecules in solution with an emphasis on extracting key structural parameters and the interpretation of SAXS data using a diverse array of techniques. These techniques cover aspects of building and assessing models to describe data measured from monodispersed and ideal dilute samples through to more complicated structurally polydisperse systems. Ab initio modelling, rigid body modelling as well as normal-mode analysis, molecular dynamics, mixed component and structural ensemble modelling are discussed. Dealing with polydispersity both physically in terms of component separation as well as approaching the analysis and modelling of data of mixtures and evolving systems are described, including methods for data decomposition such as single value decomposition/principle component analysis and evolving factor analysis. This review aims to highlight that solution SAXS, with the cohort of developments in data analysis and modelling, is well positioned to build upon the traditional 'single particle view' foundation of structural biology to take the field into new areas for interpreting the structures of proteins and biomacromolecules as population-states and dynamic structural systems.

Effects of N-glycans on the structure of human IgA2

The transition of IgA antibodies into clinical development is crucial because they have the potential to create a new class of therapeutics with superior pathogen neutralization, cancer cell killing, and immunomodulation capacity compared to IgG. However, the biological role of IgA glycans in these processes needs to be better understood. This study provides a detailed biochemical, biophysical, and structural characterization of recombinant monomeric human IgA2, which varies in the amount/locations of attached glycans. Monomeric IgA2 antibodies were produced by removing the N-linked glycans in the CH1 and CH2 domains. The impact of glycans on oligomer formation, thermal stability, and receptor binding was evaluated. In addition, we performed a structural analysis of recombinant IgA2 in solution using Small Angle X-Ray Scattering (SAXS) to examine the effect of glycans on protein structure and flexibility. Our results indicate that the absence of glycans in the Fc tail region leads to higher-order aggregates. SAXS, combined with atomistic modeling, showed that the lack of glycans in the CH2 domain results in increased flexibility between the Fab and Fc domains and a different distribution of open and closed conformations in solution. When binding with the Fcα-receptor, the dissociation constant remains unaltered in the absence of glycans in the CH1 or CH2 domain, compared to the fully glycosylated protein. These results provide insights into N-glycans' function on IgA2, which could have important implications for developing more effective IgA-based therapeutics in the future.

Resolving domain positions of cellobiose dehydrogenase by small angle X-ray scattering

The interdomain electron transfer (IET) between the catalytic flavodehydrogenase domain and the electron-transferring cytochrome domain of cellobiose dehydrogenase (CDH) plays an essential role in biocatalysis, biosensors and biofuel cells, as well as in its natural function as an auxiliary enzyme of lytic polysaccharide monooxygenase. We investigated the mobility of the cytochrome and dehydrogenase domains of CDH, which is hypothesised to limit IET in solution by small angle X-ray scattering (SAXS). CDH from Myriococcum thermophilum (syn. Crassicarpon hotsonii, syn. Thermothelomyces myriococcoides) was probed by SAXS to study the CDH mobility at different pH and in the presence of divalent cations. By comparison of the experimental SAXS data, using pair-distance distribution functions and Kratky plots, we show an increase in CDH mobility at higher pH, indicating alterations of domain mobility. To further visualise CDH movement in solution, we performed SAXS-based multistate modelling. Glycan structures present on CDH partially masked the resulting SAXS shapes, we diminished these effects by deglycosylation and studied the effect of glycoforms by modelling. The modelling shows that with increasing pH, the cytochrome domain adopts a more flexible state with significant separation from the dehydrogenase domain. On the contrary, the presence of calcium ions decreases the mobility of the cytochrome domain. Experimental SAXS data, multistate modelling and previously reported kinetic data show how pH and divalent ions impact the closed state necessary for the IET governed by the movement of the CDH cytochrome domain.

The solution structure of the unbound IgG Fc receptor CD64 resembles its crystal structure: Implications for function

FcγRI (CD64) is the only high-affinity Fcγ receptor found on monocytes, macrophages, eosinophils, neutrophils and dendritic cells. It binds immunoglobulin G (IgG) antibody-antigen complexes at its Fc region to trigger key immune responses. CD64 contains three immunoglobulin-fold extracellular domains (D1, D2 and D3) and a membrane-spanning region. Despite the importance of CD64, no solution structure for this is known to date. To investigate this, we used analytical ultracentrifugation, small-angle X-ray scattering, and atomistic modelling. Analytical ultracentrifugation revealed that CD64 was monomeric with a sedimentation coefficient s020,w of 2.53 S, together with some dimer. Small-angle X-ray scattering showed that its radius of gyration RG was 3.3-3.4 nm and increased at higher concentrations to indicate low dimerization. Monte Carlo modelling implemented in the SASSIE-web package generated 279,162 physically-realistic trial CD64 structures. From these, the scattering best-fit models at the lowest measured concentrations that minimised dimers revealed that the D1, D2 and D3 domains were structurally similar to those seen in three CD64 crystal structures, but showed previously unreported flexibility between D1, D2 and D3. Despite the limitations of the scattering data, the superimposition of the CD64 solution structures onto crystal structures of the IgG Fc-CD64 complex showed that the CD64 domains do not sterically clash with the IgG Fc region, i.e. the solution structure of CD64 was sufficiently compact to allow IgG to bind to its high-affinity Fcγ receptor. This improved understanding may result in novel approaches to inhibit CD64 function, and opens the way for the solution study of the full-length CD64-IgG complex.

Conformational multiplicity of bacterial ferric binding protein revealed by small angle x-ray scattering and molecular dynamics calculations

This study combines molecular dynamics (MD) simulations with small angle x-ray scattering (SAXS) measurements to investigate the range of conformations that can be adopted by a pH/ionic strength (IS) sensitive protein and to quantify its distinct populations in solution. To explore how the conformational distribution of proteins may be modified in the environmental niches of biological media, we focus on the periplasmic ferric binding protein A (FbpA) from Haemophilus influenzae involved in the mechanism by which bacteria capture iron from higher organisms. We examine iron-binding/release mechanisms of FbpA in varying conditions simulating its biological environment. While we show that these changes fall within the detectable range for SAXS as evidenced by differences observed in the theoretical scattering patterns calculated from the crystal structure models of apo and holo forms, detection of conformational changes due to the point mutation D52A and changes in ionic strength (IS) from SAXS scattering profiles have been challenging. Here, to reach conclusions, statistical analyses with SAXS profiles and results from different techniques were combined in a complementary fashion. The SAXS data complemented by size exclusion chromatography point to multiple and/or alternative conformations at physiological IS, whereas they are well-explained by single crystallographic structures in low IS buffers. By fitting the SAXS data with unique conformations sampled by a series of MD simulations under conditions mimicking the buffers, we quantify the populations of the occupied substates. We also find that the D52A mutant that we predicted by coarse-grained computational modeling to allosterically control the iron binding site in FbpA, responds to the environmental changes in our experiments with conformational selection scenarios that differ from those of the wild type.

Evolution of Rev7 interactions in eukaryotic TLS DNA polymerase Polζ

Translesion synthesis (TLS) DNA polymerase Polζ is crucial for the bypass replication over sites of DNA damage. The Rev7 subunit of Polζ is a HORMA (Hop1, Rev7, Mad2) protein that facilitates recruitment of Polζ to the replication fork via interactions with the catalytic subunit Rev3 and the translesion synthesis scaffold protein Rev1. Human Rev7 (hRev7) interacts with two Rev7-binding motifs (RBMs) of hRev3 by a mechanism conserved among HORMA proteins whereby the safety-belt loop of hRev7 closes on the top of the ligand. The two copies of hRev7 tethered by the two hRev3-RBMs form a symmetric head-to-head dimer through the canonical HORMA dimerization interface. Recent cryo-EM structures reveal that Saccharomyces cerevisiae Polζ (scPolζ) also includes two copies of scRev7 bound to distinct regions of scRev3. Surprisingly, the HORMA dimerization interface is not conserved in scRev7, with the two scRev7 protomers forming an asymmetric head-to-tail dimer with a much smaller interface than the hRev7 dimer. Here, we validated the two adjacent RBM motifs in scRev3, which bind scRev7 with affinities that differ by two orders of magnitude and confirmed the 2:1 stoichiometry of the scRev7:Rev3 complex in solution. However, our biophysical studies reveal that scRev7 does not form dimers in solution either on its own accord or when tethered by the two RBMs in scRev3. These findings imply that the scRev7 dimer observed in the cryo-EM structures is induced by scRev7 interactions with other Polζ subunits and that Rev7 homodimerization via the HORMA interface is a mechanism that emerged later in evolution.

Interdomain Linkers Regulate Histidine Kinase Activity by Controlling Subunit Interactions

Bacterial chemoreceptors regulate the cytosolic multidomain histidine kinase CheA through largely unknown mechanisms. Residue substitutions in the peptide linkers that connect the P4 kinase domain to the P3 dimerization and P5 regulatory domain affect CheA basal activity and activation. To understand the role that these linkers play in CheA activity, the P3-to-P4 linker (L3) and P4-to-P5 linker (L4) were extended and altered in variants of Thermotoga maritima (Tm) CheA. Flexible extensions of the L3 and L4 linkers in CheA-LV1 (linker variant 1) allowed for a well-folded kinase domain that retained wild-type (WT)-like binding affinities for nucleotide and normal interactions with the receptor-coupling protein CheW. However, CheA-LV1 autophosphorylation activity registered ∼50-fold lower compared to WT. Neither a WT nor LV1 dimer containing a single P4 domain could autophosphorylate the P1 substrate domain. Autophosphorylation activity was rescued in variants with extended L3 and L4 linkers that favor helical structure and heptad spacing. Autophosphorylation depended on linker spacing and flexibility and not on sequence. Pulse-dipolar electron-spin resonance (ESR) measurements with spin-labeled adenosine 5'-triphosphate (ATP) analogues indicated that CheA autophosphorylation activity inversely correlated with the proximity of the P4 domains within the dimers of the variants. Despite their separation in primary sequence and space, the L3 and L4 linkers also influence the mobility of the P1 substrate domains. In all, interactions of the P4 domains, as modulated by the L3 and L4 linkers, affect domain dynamics and autophosphorylation of CheA, thereby providing potential mechanisms for receptors to regulate the kinase.

Structural interpretations of a flexible cold-active AMS8 lipase by combining small-angle X-ray scattering and molecular dynamics simulation (SAXS-MD)

Determining structure of highly flexible protein with multiple conformations can be challenging. This paper aims to combine molecular dynamics (MD) and small angle X-ray diffraction (SAX) techniques as a solution to overcome issues related to protein conformation in hardly crystallized protein. Based on prior studies, a cold-active lipase AMS8 was simulated in solvents showing stability in its N-terminal and high flexibility in its C-terminal. However, MD in its own algorithm could not explain the basis of macromolecule conformational transitions or changes related to protein through folding. Hence, by combining SAXS with MD, it is possible to understand the structure of flexible AMS8 lipase in natural space. Based on the findings, SAXS ab-initio model of AMS8 lipase was identified as a monomeric protein in which the optimized model of cold-active lipase AMS8 derived from SAXS data was found to be aligned with AMS8 homology model under series of MD timeframe.

SAXS analysis of the intrinsic tenase complex bound to a lipid nanodisc highlights intermolecular contacts between factors VIIIa/IXa

The intrinsic tenase (Xase) complex, formed by factors (f) VIIIa and fIXa, forms on activated platelet surfaces and catalyzes the activation of factor X to Xa, stimulating thrombin production in the blood coagulation cascade. The structural organization of the membrane-bound Xase complex remains largely unknown, hindering our understanding of the structural underpinnings that guide Xase complex assembly. Here, we aimed to characterize the Xase complex bound to a lipid nanodisc with biolayer interferometry (BLI), Michaelis-Menten kinetics, and small-angle X-ray scattering (SAXS). Using immobilized lipid nanodiscs, we measured binding rates and nanomolar affinities for fVIIIa, fIXa, and the Xase complex. Enzyme kinetic measurements demonstrated the assembly of an active enzyme complex in the presence of lipid nanodiscs. An ab initio molecular envelope of the nanodisc-bound Xase complex allowed us to computationally model fVIIIa and fIXa docked onto a flexible lipid membrane and identify protein-protein interactions. Our results highlight multiple points of contact between fVIIIa and fIXa, including a novel interaction with fIXa at the fVIIIa A1-A3 domain interface. Lastly, we identified hemophilia A/B-related mutations with varying severities at the fVIIIa/fIXa interface that may regulate Xase complex assembly. Together, our results support the use of SAXS as an emergent tool to investigate the membrane-bound Xase complex and illustrate how mutations at the fVIIIa/fIXa dimer interface may disrupt or stabilize the activated enzyme complex.

Untangling the Conformational Plasticity of V66M Human proBDNF Polymorphism as a Modifier of Psychiatric Disorder Susceptibility

The human genetic variant BDNF (V66M) represents the first example of neurotrophin family member that has been linked to psychiatric disorders. In order to elucidate structural differences that account for the effects in cognitive function, this hproBDNF polymorph was expressed, refolded, purified, and compared directly to the WT variant for the first time for differences in their 3D structures by DSF, limited proteolysis, FT-IR, and SAXS measurements in solution. Our complementary studies revealed a deep impact of V66M polymorphism on hproBDNF conformations in solution. Although the mean conformation in solution appears to be more compact in the V66M variant, overall, we demonstrated a large increase in flexibility in solution upon V66M mutation. Thus, considering that plasticity in IDR is crucial for protein function, the observed alterations may be related to the functional alterations in hproBDNF binding to its receptors p75NTR, sortilin, HAP1, and SorCS2. These effects can provoke altered intracellular neuronal trafficking and/or affect proBDNF physiological functions, leading to many brain-associated diseases and conditions such as cognitive impairment and anxiety. The structural alterations highlighted in the present study may pave the way to the development of drug discovery strategies to provide greater therapeutic responses and of novel pharmacologic strategy in human populations with this common polymorphism, ultimately guiding personalized medicine for neuropsychiatric disorders.

Domain architecture and catalysis of the Staphylococcus aureus fatty acid kinase

Fatty acid kinase (Fak) is a two-component enzyme that generates acyl-phosphate for phospholipid synthesis. Fak consists of a kinase domain protein (FakA) that phosphorylates a fatty acid enveloped by a fatty acid binding protein (FakB). The structural basis for FakB function has been established, but little is known about FakA. Here, we used limited proteolysis to define three separate FakA domains: the amino terminal FakA_N, the central FakA_L, and the carboxy terminal FakA_C. The isolated domains lack kinase activity, but activity is restored when FakA_N and FakA_L are present individually or connected as FakA_NL. The X-ray structure of the monomeric FakA_N captures the product complex with ADP and two Mg2+ ions bound at the nucleotide site. The FakA_L domain encodes the dimerization interface along with conserved catalytic residues Cys240, His282, and His284. AlphaFold analysis of FakA_L predicts the catalytic residues are spatially clustered and pointing away from the dimerization surface. Furthermore, the X-ray structure of FakA_C shows that it consists of two subdomains that are structurally related to FakB. Analytical ultracentrifugation demonstrates that FakA_C binds FakB, and site-directed mutagenesis confirms that a positively charged wedge on FakB meshes with a negatively charged groove on FakA_C. Finally, small angle X-ray scattering analysis is consistent with freely rotating FakA_N and FakA_C domains tethered by flexible linkers to FakA_L. These data reveal specific roles for the three independently folded FakA protein domains in substrate binding and catalysis.

NaV1.2 EFL domain allosterically enhances Ca2+ binding to sites I and II of WT and pathogenic calmodulin mutants bound to the channel CTD

Neuronal voltage-gated sodium channel NaV1.2 C-terminal domain (CTD) binds calmodulin (CaM) constitutively at its IQ motif. A solution structure (6BUT) and other NMR evidence showed that the CaM N domain (CaMN) is structurally independent of the C-domain (CaMC) whether CaM is bound to the NaV1.2IQp (1,901-1,927) or NaV1.2CTD (1,777-1,937) with or without calcium. However, in the CaM + NaV1.2CTD complex, the Ca2+ affinity of CaMN was more favorable than in free CaM, while Ca2+ affinity for CaMC was weaker than in the CaM + NaV1.2IQp complex. The CTD EF-like (EFL) domain allosterically widened the energetic gap between CaM domains. Cardiomyopathy-associated CaM mutants (N53I(N54I), D95V(D96V), A102V(A103V), E104A(E105A), D129G(D130G), and F141L(F142L)) all bound the NaV1.2 IQ motif favorably under resting (apo) conditions and bound calcium normally at CaMN sites. However, only N53I and A102V bound calcium at CaMC sites at [Ca2+] < 100 μM. Thus, they are expected to respond like wild-type CaM to Ca2+ spikes in excitable cells.

The role of SAXS and molecular simulations in 3D structure elucidation of a DNA aptamer against lung cancer

Aptamers are short, single-stranded DNA or RNA oligonucleotide molecules that function as synthetic analogs of antibodies and bind to a target molecule with high specificity. Aptamer affinity entirely depends on its tertiary structure and charge distribution. Therefore, length and structure optimization are essential for increasing aptamer specificity and affinity. Here, we present a general optimization procedure for finding the most populated atomistic structures of DNA aptamers. Based on the existed aptamer LC-18 for lung adenocarcinoma, a new truncated LC-18 (LC-18t) aptamer LC-18t was developed. A three-dimensional (3D) shape of LC-18t was reported based on small-angle X-ray scattering (SAXS) experiments and molecular modeling by fragment molecular orbital or molecular dynamic methods. Molecular simulations revealed an ensemble of possible aptamer conformations in solution that were in close agreement with measured SAXS data. The aptamer LC-18t had stronger binding to cancerous cells in lung tumor tissues and shared the binding site with the original larger aptamer. The suggested approach reveals 3D shapes of aptamers and helps in designing better affinity probes.

Albumin in patients with liver disease shows an altered conformation

Human serum albumin (HSA) constitutes the primary transporter of fatty acids, bilirubin, and other plasma compounds. The binding, transport, and release of its cargos strongly depend on albumin conformation, which is affected by bound ligands induced by physiological and pathological conditions. HSA is both highly oxidized and heavily loaded with fatty acids and bilirubin in chronic liver disease. By employing small-angle X-ray scattering we show that HSA from the plasma of chronic liver disease patients undergoes a distinct opening compared to healthy donors. The extent of HSA opening correlates with clinically relevant variables, such as the model of end-stage liver disease score, bilirubin, and fatty acid levels. Although the mild oxidation of HSA in vitro does not alter overall structure, the alteration of patients’ HSA correlates with its redox state. This study connects clinical data with structural visualization of albumin dynamicity in solution and underlines the functional importance of albumin’s inherent flexibility.

Paar et al. propose a SAXS-based approach to study conformations of human serum albumin (HSA) from patients with liver disease and a structural understanding of HSA dynamicity and its correlation with clinical variables are provided. Using it on real clinical samples, this study has concrete practical implications too.

X-ray scattering reveals disordered linkers and dynamic interfaces in complexes and mechanisms for DNA double-strand break repair impacting cell and cancer biology

Evolutionary selection ensures specificity and efficiency in dynamic metastable macromolecular machines that repair DNA damage without releasing toxic and mutagenic intermediates. Here we examine non‐homologous end joining (NHEJ) as the primary conserved DNA double‐strand break (DSB) repair process in human cells. NHEJ has exemplary key roles in networks determining the development, outcome of cancer treatments by DSB‐inducing agents, generation of antibody and T‐cell receptor diversity, and innate immune response for RNA viruses. We determine mechanistic insights into NHEJ structural biochemistry focusing upon advanced small angle X‐ray scattering (SAXS) results combined with X‐ray crystallography (MX) and cryo‐electron microscopy (cryo‐EM). SAXS coupled to atomic structures enables integrated structural biology for objective quantitative assessment of conformational ensembles and assemblies in solution, intra‐molecular distances, structural similarity, functional disorder, conformational switching, and flexibility. Importantly, NHEJ complexes in solution undergo larger allosteric transitions than seen in their cryo‐EM or MX structures. In the long‐range synaptic complex, X‐ray repair cross‐complementing 4 (XRCC4) plus XRCC4‐like‐factor (XLF) form a flexible bridge and linchpin for DNA ends bound to KU heterodimer (Ku70/80) and DNA‐PKcs (DNA‐dependent protein kinase catalytic subunit). Upon binding two DNA ends, auto‐phosphorylation opens DNA‐PKcs dimer licensing NHEJ via concerted conformational transformations of XLF‐XRCC4, XLF–Ku80, and LigIV^BRCT–Ku70 interfaces. Integrated structures reveal multifunctional roles for disordered linkers and modular dynamic interfaces promoting DSB end processing and alignment into the short‐range complex for ligation by LigIV. Integrated findings define dynamic assemblies fundamental to designing separation‐of‐function mutants and allosteric inhibitors targeting conformational transitions in multifunctional complexes.

Molecular model of a sensor of two-component signaling system

Two-component systems (TCS) are widespread signaling systems present in all domains of life. TCS typically consist of a signal receptor/transducer and a response regulator. The receptors (histidine kinases, chemoreceptors and photoreceptors) are often embedded in the membrane and have a similar modular structure. Chemoreceptors were shown to function in highly ordered arrays, with trimers of dimers being the smallest functional unit. However, much less is known about photoreceptors. Here, we use small-angle scattering (SAS) to show that detergent-solubilized sensory rhodopsin II in complex with its cognate transducer forms dimers at low salt concentration, which associate into trimers of dimers at higher buffer molarities. We then fit an atomistic model of the whole complex into the SAS data. The obtained results suggest that the trimer of dimers is "tripod"-shaped and that the contacts between the dimers occur only through their cytoplasmic regions, whereas the transmembrane regions remain unconnected.

Assessment of AMBER Force Fields for Simulations of ssDNA

Single-stranded DNA (ssDNA) plays an important role in biological processes and is used in DNA nanotechnology and other novel applications. Many important research questions can be addressed with molecular simulations of ssDNA molecules; however, no dedicated force field for ssDNA has been developed, and there is limited experimental information about ssDNA structures. This study assesses the accuracy and applicability of existing Amber force fields for all-atom simulations of ssDNA, such as ff99, bsc0, bsc1, and OL15, in implicit and explicit solvents via comparison to available experimental data, such as Forster resonance energy transfer and small angle X-ray scattering. We observed that some force fields agree better with experiments than others mainly due to the difference in parameterization of the propensity for hydrogen bonding and base stacking. Overall, the Amber ff99 force field in the IGB5 or IGB8 implicit solvent and the bsc1 force field in the explicit TIP3P solvent had the best agreement with experiment.

Rigid monoclonal antibodies improve detection of SARS-CoV-2 nucleocapsid protein

Monoclonal antibodies (mAbs) are the basis of treatments and diagnostics for pathogens and other biological phenomena. We conducted a structural characterization of mAbs against the N-terminal domain of nucleocapsid protein (NPNTD) from SARS-CoV-2 using small-angle X-ray scattering and transmission electron microscopy. Our solution-based results distinguished the mAbs' flexibility and how this flexibility affects the assembly of multiple mAbs on an antigen. By pairing two mAbs that bind different epitopes on the NPNTD, we show that flexible mAbs form a closed sandwich-like complex. With rigid mAbs, a juxtaposition of the antigen-binding fragments is prevented, enforcing a linear arrangement of the mAb pair, which facilitates further mAb polymerization. In a modified sandwich enzyme-linked immunosorbent assay, we show that rigid mAb-pairings with linear polymerization led to increased NPNTD detection sensitivity. These enhancements can expedite the development of more sensitive and selective antigen-detecting point-of-care lateral flow devices, which are critical for early diagnosis and epidemiological studies of SARS-CoV-2 and other pathogens.

Visualizing functional dynamicity in the DNA-dependent protein kinase holoenzyme DNA-PK complex by integrating SAXS with cryo-EM

Assembly of KU and DNA-dependent protein kinase catalytic subunit (DNA-PKcs) at DNA double strand breaks (DSBs) forms DNA-PK holoenzyme as a critical initiating step for non-homologous end joining (NHEJ) repair of DSBs produced by radiation and chemotherapies. Advanced cryo-electron microscopy (cryo-EM) imaging together with breakthrough macromolecular X-ray crystal (MX) structures of KU and DNA-PKcs recently enabled visualization of the ∼600 kDa DNA-PK assembly at near atomic resolution. These important static structures provide the foundation for definition and interpretation of functional movements crucial to mechanistic understanding that can be tested through solution state structure analysis. We herein therefore leverage Cryo-EM and MX structures for the interpretation of synchrotron small-angle X-ray scattering (SAXS) data on DNA-PK conformations in solution to inform the structural mechanism for NHEJ initiation. SAXS, which measures thermodynamic solution-state conformational states and assemblies outside of cryo- and solid-state conditions, unveils the inherent flexibility of KU, DNA-PKcs and DNA-PK. The combined structural measurements reveal mobility of KU80 C-terminal region (KU80CTR), motion/plasticity of HEAT (DNA-PKcs Huntingtin, Elongation Factor 3, PP2 A, and TOR1) regions, allosteric switching upon DNA-PKcs autophosphorylation, and dimeric arrangements of DNA-PK assembly. Importantly, the results uncover displacement of the N-terminal HEAT domain during autophosphorylation as suitable for a regulated release mechanism of DNA-PKcs from DNA-PK to control unproductive access to toxic and mutagenic DNA repair intermediates. These integrated analyses show that the marriage of SAXS with cryo-EM leverages the strengths of both techniques to enable assessment of functional conformations and flexibility defining atomic-resolution molecular mechanisms for DSB repair.

Nucleoid remodeling during environmental adaptation is regulated by HU-dependent DNA bundling

Bacterial nucleoid remodeling dependent on conserved histone-like protein, HU is one of the determining factors in global gene regulation. By imaging of near-native, unlabeled E. coli cells by soft X-ray tomography, we show that HU remodels nucleoids by promoting the formation of a dense condensed core surrounded by less condensed isolated domains. Nucleoid remodeling during cell growth and environmental adaptation correlate with pH and ionic strength controlled molecular switch that regulated HUαα dependent intermolecular DNA bundling. Through crystallographic and solution-based studies we show that these effects mechanistically rely on HUαα promiscuity in forming multiple electrostatically driven multimerization interfaces. Changes in DNA bundling consequently affects gene expression globally, likely by constrained DNA supercoiling. Taken together our findings unveil a critical function of HU–DNA interaction in nucleoid remodeling that may serve as a general microbial mechanism for transcriptional regulation to synchronize genetic responses during the cell cycle and adapt to changing environments.

HU is among the most conserved and abundant nucleoid-associated proteins in eubacteria. Here the authors investigate the role of histone-like proteins (HU) in the 3D organization of the bacteria DNA and show via soft X-ray tomography the process of nucleoid remodeling.

Crystal and solution structures of fragments of the human leucocyte common antigen-related protein

Leucocyte common antigen-related protein (LAR) is a post-synaptic type I transmembrane receptor protein that is important for neuronal functionality and is genetically coupled to neuronal disorders such as attention deficit hyperactivity disorder (ADHD). To understand the molecular function of LAR, structural and biochemical studies of protein fragments derived from the ectodomain of human LAR have been performed. The crystal structure of a fragment encompassing the first four FNIII domains (LAR^FN1-4) showed a characteristic L shape. SAXS data suggested limited flexibility within LAR^FN1-4, while rigid-body refinement of the SAXS data using the X-ray-derived atomic model showed a smaller angle between the domains defining the L shape compared with the crystal structure. The capabilities of the individual LAR fragments to interact with heparin was examined using microscale thermophoresis and heparin-affinity chromatography. The results showed that the three N-terminal immunoglobulin domains (LAR^Ig1-3) and the four C-terminal FNIII domains (LAR^FN5-8) both bound heparin, while LAR^FN1-4 did not. The low-molecular-weight heparin drug Innohep induced a shift in hydrodynamic volume as assessed by size-exclusion chromatography of LAR^Ig1-3 and LAR^FN5-8, while the chemically defined pentameric heparin drug Arixtra did not. Together, the presented results suggest the presence of an additional heparin-binding site in human LAR.

Generation of the configurational ensemble of an intrinsically disordered protein from unbiased molecular dynamics simulation

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes, play a major role in cell signaling, and are associated with human diseases. To understand IDP function it is critical to determine their configurational ensemble, i.e., the collection of 3-dimensional structures they adopt, and this remains an immense challenge in structural biology. Attempts to determine this ensemble computationally have been hitherto hampered by the necessity of reweighting molecular dynamics (MD) results or biasing simulation in order to match ensemble-averaged experimental observables, operations that reduce the precision of the generated model because different structural ensembles may yield the same experimental observable. Here, by employing enhanced sampling MD we reproduce the experimental small-angle neutron and X-ray scattering profiles and the NMR chemical shifts of the disordered N terminal (SH4UD) of c-Src kinase without reweighting or constraining the simulations. The unbiased simulation results reveal a weakly funneled and rugged free energy landscape of SH4UD, which gives rise to a heterogeneous ensemble of structures that cannot be described by simple polymer theory. SH4UD adopts transient helices, which are found away from known phosphorylation sites and could play a key role in the stabilization of structural regions necessary for phosphorylation. Our findings indicate that adequately sampled molecular simulations can be performed to provide accurate physical models of flexible biosystems, thus rationalizing their biological function.

Four steps for revealing and adjusting the 3D structure of aptamers in solution by small-angle X-ray scattering and computer simulation

Nucleic acid (NA) aptamers bind to their targets with high affinity and selectivity. The three-dimensional (3D) structures of aptamers play a major role in these non-covalent interactions. Here, we use a four-step approach to determine a true 3D structure of aptamers in solution using small-angle X-ray scattering (SAXS) and molecular structure restoration (MSR). The approach consists of (i) acquiring SAXS experimental data of an aptamer in solution, (ii) building a spatial distribution of the molecule's electron density using SAXS results, (iii) constructing a 3D model of the aptamer from its nucleotide primary sequence and secondary structure, and (iv) comparing and refining the modeled 3D structures with the experimental SAXS model. In the proof-of-principle we analyzed the 3D structure of RE31 aptamer to thrombin in a native free state at different temperatures and validated it by circular dichroism (CD). The resulting 3D structure of RE31 has the most energetically favorable conformation and the same elements such as a B-form duplex, non-complementary region, and two G-quartets which were previously reported by X-ray diffraction (XRD) from a single crystal. More broadly, this study demonstrates the complementary approach for constructing and adjusting the 3D structures of aptamers, DNAzymes, and ribozymes in solution, and could supply new opportunities for developing functional nucleic acids. Graphical abstract.

A super-Gaussian Poisson-Boltzmann model for electrostatic free energy calculation: smooth dielectric distribution for protein cavities and in both water and vacuum states

Calculations of electrostatic potential and solvation free energy of macromolecules are essential for understanding the mechanism of many biological processes. In the classical implicit solvent Poisson-Boltzmann (PB) model, the macromolecule and water are modeled as two-dielectric media with a sharp border. However, the dielectric property of interior cavities and ion-channels is difficult to model realistically in a two-dielectric setting. In fact, the detection of water molecules in a protein cavity remains to be an experimental challenge. This introduces an uncertainty, which affects the subsequent solvation free energy calculation. In order to compensate this uncertainty, a novel super-Gaussian dielectric PB model is introduced in this work, which devices an inhomogeneous dielectric distribution to represent the compactness of atoms and characterizes empty cavities via a gap dielectric value. Moreover, the minimal molecular surface level set function is adopted so that the dielectric profile remains to be smooth when the protein is transferred from water phase to vacuum. An important feature of this new model is that as the order of super-Gaussian function approaches the infinity, the dielectric distribution reduces to a piecewise constant of the two-dielectric model. Mathematically, an effective dielectric constant analysis is introduced in this work to benchmark the dielectric model and select optimal parameter values. Computationally, a pseudo-time alternative direction implicit (ADI) algorithm is utilized for solving the super-Gaussian PB equation, which is found to be unconditionally stable in a smooth dielectric setting. Solvation free energy calculation of a Kirkwood sphere and various proteins is carried out to validate the super-Gaussian model and ADI algorithm. One macromolecule with both water filled and empty cavities is employed to demonstrate how the cavity uncertainty in protein structure can be bypassed through dielectric modeling in biomolecular electrostatic analysis.

Functional Relevance of Interleukin-1 Receptor Inter-domain Flexibility for Cytokine Binding and Signaling

The interleukin 1 (IL-1) receptor family, whose members contain three immunoglobulin-like domains (D1-D3) in the extracellular region, is responsible for transmitting pleiotropic signals of IL-1 cytokines. The inter-domain flexibility of IL-1 receptors and its functional roles have not been fully elucidated. In this study, we used small-angle X-ray scattering to show that ligand-binding primary receptors and co-receptors in the family all have inherent inter-domain flexibility due to the D2/D3 linker. Variants of the IL-1RAcP and IL-18Rβ co-receptors with mutated D2/D3 linkers cannot form a cytokine-receptor complex and mediate signaling. Our analysis further revealed that these mutated co-receptors exhibited a changed conformational ensemble, suggesting that loss of function is due to the alteration of receptor dynamics. Taken together, our results demonstrate that the D2/D3 linker is a critical functional determinant of IL-1 receptor and underscore the important roles of the inter-domain flexibility in cytokine/receptor binding and signaling.

CH2 domain orientation of human immunoglobulin G in solution: Structural comparison of glycosylated and aglycosylated Fc regions using small-angle X-ray scattering

The N-linked glycan in immunoglobulin G is critical for the stability and function of the crystallizable fragment (Fc) region. Alteration of these protein properties upon the removal of the N-linked glycan has often been explained by the alteration of the CH2 domain orientation in the Fc region. To confirm this hypothesis, we examined the small-angle X-ray scattering (SAXS) profile of the glycosylated Fc region (gFc) and aglycosylated Fc region (aFc) in solution. Conformational characteristics of the CH2 domain orientation were validated by comparison with SAXS profiles theoretically calculated from multiple crystal structures of the Fc region with different CH2 domain orientations. The reduced chi-square values from the fitting analyses of gFc and aFc associated with the degree of openness or closure of each crystal structure, as determined from the first principal component that partially governed the variation of the CH2 domain orientation extracted by a singular value decomposition analysis. For both gFc and aFc, the best-fitted SAXS profiles corresponded to ones calculated based on the crystal structure of gFc that formed a "semi-closed" CH2 domain orientation. Collectively, the data indicated that the removal of the N-linked glycan only negligibly affected the CH2 domain orientation in solution. These findings will guide the development of methodology for the production of highly refined functional Fc variants.

Domain cross-talk within a bifunctional enzyme provides catalytic and allosteric functionality in the biosynthesis of aromatic amino acids

Because of their special organization, multifunctional enzymes play crucial roles in improving the performance of metabolic pathways. For example, the bacterium Prevotella nigrescens contains a distinctive bifunctional protein comprising a 3-deoxy-d-arabino heptulosonate-7-phosphate synthase (DAH7PS), catalyzing the first reaction of the biosynthetic pathway of aromatic amino acids, and a chorismate mutase (CM), functioning at a branch of this pathway leading to the synthesis of tyrosine and phenylalanine. In this study, we characterized this P. nigrescens enzyme and found that its two catalytic activities exhibit substantial hetero-interdependence and that the separation of its two distinct catalytic domains results in a dramatic loss of both DAH7PS and CM activities. The protein displayed a unique dimeric assembly, with dimerization solely via the CM domain. Small angle X-ray scattering (SAXS)-based structural analysis of this protein indicated a DAH7PS-CM hetero-interaction between the DAH7PS and CM domains, unlike the homo-association between DAH7PS domains normally observed for other DAH7PS proteins. This hetero-interaction provides a structural basis for the functional interdependence between the two domains observed here. Moreover, we observed that DAH7PS is allosterically inhibited by prephenate, the product of the CM-catalyzed reaction. This allostery was accompanied by a striking conformational change as observed by SAXS, implying that altering the hetero-domain interaction underpins the allosteric inhibition. We conclude that for this C-terminal CM-linked DAH7PS, catalytic function and allosteric regulation appear to be delivered by a common mechanism, revealing a distinct and efficient evolutionary strategy to utilize the functional advantages of a bifunctional enzyme.

Biochemical characterization and low-resolution SAXS structure of two-domain endoglucanase BlCel9 from Bacillus licheniformis

Lignocellulose feedstock constitutes the most abundant carbon source in the biosphere; however, its recalcitrance remains a challenge for microbial conversion into biofuel and bioproducts. Bacillus licheniformis is a microbial mesophilic bacterium capable of secreting a large number of glycoside hydrolase (GH) enzymes, including a glycoside hydrolase from GH family 9 (BlCel9). Here, we conducted biochemical and biophysical studies of recombinant BlCel9, and its low-resolution molecular shape was retrieved from small angle X-ray scattering (SAXS) data. BlCel9 is an endoglucanase exhibiting maximum catalytic efficiency at pH 7.0 and 60 °C. Furthermore, it retains 80% of catalytic activity within a broad range of pH values (5.5-8.5) and temperatures (up to 50 °C) for extended periods of time (over 48 h). It exhibits the highest hydrolytic activity against phosphoric acid swollen cellulose (PASC), followed by bacterial cellulose (BC), filter paper (FP), and to a lesser extent carboxymethylcellulose (CMC). The HPAEC-PAD analysis of the hydrolytic products demonstrated that the end product of the enzymatic hydrolysis is primarily cellobiose, and also small amounts of glucose, cellotriose, and cellotetraose are produced. SAXS data analysis revealed that the enzyme adopts a monomeric state in solution and has a molecular mass of 65.8 kDa as estimated from SAXS data. The BlCel9 has an elongated shape composed of an N-terminal family 3 carbohydrate-binding module (CBM3c) and a C-terminal GH9 catalytic domain joined together by 20 amino acid residue long linker peptides. The domains are closely juxtaposed in an extended conformation and form a relatively rigid structure in solution, indicating that the interactions between the CBM3c and GH9 catalytic domains might play a key role in cooperative cellulose biomass recognition and hydrolysis.

Bayesian inference of protein conformational ensembles from limited structural data

Many proteins consist of folded domains connected by regions with higher flexibility. The details of the resulting conformational ensemble play a central role in controlling interactions between domains and with binding partners. Small-Angle Scattering (SAS) is well-suited to study the conformational states adopted by proteins in solution. However, analysis is complicated by the limited information content in SAS data and care must be taken to avoid constructing overly complex ensemble models and fitting to noise in the experimental data. To address these challenges, we developed a method based on Bayesian statistics that infers conformational ensembles from a structural library generated by all-atom Monte Carlo simulations. The first stage of the method involves a fast model selection based on variational Bayesian inference that maximizes the model evidence of the selected ensemble. This is followed by a complete Bayesian inference of population weights in the selected ensemble. Experiments with simulated ensembles demonstrate that model evidence is capable of identifying the correct ensemble and that correct number of ensemble members can be recovered up to high level of noise. Using experimental data, we demonstrate how the method can be extended to include data from Nuclear Magnetic Resonance (NMR) and structural energies of conformers extracted from the all-atom energy functions. We show that the data from SAXS, NMR chemical shifts and energies calculated from conformers can work synergistically to improve the definition of the conformational ensemble.

Structure of an H1-Bound 6-Nucleosome Array Reveals an Untwisted Two-Start Chromatin Fiber Conformation

Chromatin adopts a diversity of regular and irregular fiber structures in vitro and in vivo. However, how an array of nucleosomes folds into and switches between different fiber conformations is poorly understood. We report the 9.7 Å resolution crystal structure of a 6-nucleosome array bound to linker histone H1 determined under ionic conditions that favor incomplete chromatin condensation. The structure reveals a flat two-start helix with uniform nucleosomal stacking interfaces and a nucleosome packing density that is only half that of a twisted 30-nm fiber. Hydroxyl radical footprinting indicates that H1 binds the array in an on-dyad configuration resembling that observed for mononucleosomes. Biophysical, cryo-EM, and crosslinking data validate the crystal structure and reveal that a minor change in ionic environment shifts the conformational landscape to a more compact, twisted form. These findings provide insights into the structural plasticity of chromatin and suggest a possible assembly pathway for a 30-nm fiber.

Conformational Flexibility of Ubiquitin-Modified and SUMO-Modified PCNA Shown by Full-Ensemble Hybrid Methods

Ubiquitin-modified proliferating cell nuclear antigen (PCNA) and small ubiquitin-like modifier (SUMO)-modified PCNA regulate DNA damage tolerance pathways. X-ray crystal structures of these proteins suggested that they do not have much conformational flexibility because the modifiers have preferred binding sites on the surface of PCNA. By contrast, small-angle X-ray scattering analyses of these proteins suggested that they have different degrees of conformational flexibility, with SUMO-modified PCNA being more flexible. These conclusions were based on minimal-ensemble hybrid approaches, which produce unrealistic models by representing flexible proteins with only a few static structures. To overcome the limitations of minimal-ensemble hybrid approaches and to determine the degree of conformational flexibility of ubiquitin-modified PCNA and SUMO-modified PCNA, we utilized a novel full-ensemble hybrid approach. We carried out molecular simulations and small-angle X-ray scattering analyses of both proteins and obtained outstanding agreement between the full ensembles generated by the simulations and the experimental data. We found that both proteins have a high degree of conformational flexibility. The modifiers occupy many positions around the back and side of the PCNA ring. Moreover, we found no preferred ubiquitin-binding or SUMO-binding sites on PCNA. This conformational flexibility likely facilitates the recognition of downstream effector proteins and the formation of PCNA tool belts.

A methodology to calculate small-angle scattering profiles of macromolecular solutions from molecular simulations in the grand-canonical ensemble

The theoretical framework to evaluate small-angle scattering (SAS) profiles for multi-component macromolecular solutions is re-examined from the standpoint of molecular simulations in the grand-canonical ensemble, where the chemical potentials of all species in solution are fixed. This statistical mechanical ensemble resembles more closely scattering experiments, capturing concentration fluctuations that arise from the exchange of molecules between the scattering volume and the bulk solution. The resulting grand-canonical expression relates scattering intensities to the different intra- and intermolecular pair distribution functions, as well as to the distribution of molecular concentrations on the scattering volume. This formulation represents a generalized expression that encompasses most of the existing methods to evaluate SAS profiles from molecular simulations. The grand-canonical SAS methodology is probed for a series of different implicit-solvent, homogeneous systems at conditions ranging from dilute to concentrated. These systems consist of spherical colloids, dumbbell particles, and highly flexible polymer chains. Comparison of the resulting SAS curves against classical methodologies based on either theoretical approaches or canonical simulations (i.e., at a fixed number of molecules) shows equivalence between the different scattering intensities so long as interactions between molecules are net repulsive or weakly attractive. On the other hand, for strongly attractive interactions, grand-canonical SAS profiles deviate in the low- and intermediate-q range from those calculated in a canonical ensemble. Such differences are due to the distribution of molecules becoming asymmetric, which yields a higher contribution from configurations with molecular concentrations larger than the nominal value. Additionally, for flexible systems, explicit discrimination between intra- and inter-molecular SAS contributions permits the implementation of model-free, structural analysis such as Guinier's plots at high molecular concentrations, beyond what the traditional limits are for such analysis.

Conformational Plasticity of the Immunoglobulin Fc Domain in Solution

Fragment crystallizable (Fc) region of immunoglobulin G (IgG) antibody binds to specific Fc receptors (FcγRs) to control antibody effector functions. Currently, engineered specific Fc-FcγR interactions are validated with a static conformation derived from the crystal structure. However, computational evidence suggests that the conformational variability of Fcs plays an important role in receptor recognition. Here we elucidate Fc flexibility of IgG1, IgG2, and IgG1 Fc with mutations (M255Y/S257T/T259E) in solution by small-angle X-ray scattering (SAXS). Measured SAXS profiles and experimental parameters show variations in flexibility between Fc isotypes. We develop and apply a modeling tool for an accurate conformational sampling of Fcs followed by SAXS fitting. Revealed conformational variability of the CH2 domain as low as 10 Å in displacement, illustrates the power of the atomistic modeling combined with SAXS. This inexpensive SAXS-based approach offers to improve the engineering of antibodies for tailoring Fc receptor interactions through altering and measuring Fc flexibility.

A Calmodulin C-Lobe Ca2+-Dependent Switch Governs Kv7 Channel Function

Kv7 (KCNQ) voltage-gated potassium channels control excitability in the brain, heart, and ear. Calmodulin (CaM) is crucial for Kv7 function, but how this calcium sensor affects activity has remained unclear. Here, we present X-ray crystallographic analysis of CaM:Kv7.4 and CaM:Kv7.5 AB domain complexes that reveal an Apo/CaM clamp conformation and calcium binding preferences. These structures, combined with small-angle X-ray scattering, biochemical, and functional studies, establish a regulatory mechanism for Kv7 CaM modulation based on a common architecture in which a CaM C-lobe calcium-dependent switch releases a shared Apo/CaM clamp conformation. This C-lobe switch inhibits voltage-dependent activation of Kv7.4 and Kv7.5 but facilitates Kv7.1, demonstrating that mechanism is shared by Kv7 isoforms despite the different directions of CaM modulation. Our findings provide a unified framework for understanding how CaM controls different Kv7 isoforms and highlight the role of membrane proximal domains for controlling voltage-gated channel function. VIDEO ABSTRACT.

Modeling Structure and Dynamics of Protein Complexes with SAXS Profiles

Small-angle X-ray scattering (SAXS) is an increasingly common and useful technique for structural characterization of molecules in solution. A SAXS experiment determines the scattering intensity of a molecule as a function of spatial frequency, termed SAXS profile. SAXS profiles can be utilized in a variety of molecular modeling applications, such as comparing solution and crystal structures, structural characterization of flexible proteins, assembly of multi-protein complexes, and modeling of missing regions in the high-resolution structure. Here, we describe protocols for modeling atomic structures based on SAXS profiles. The first protocol is for comparing solution and crystal structures including modeling of missing regions and determination of the oligomeric state. The second protocol performs multi-state modeling by finding a set of conformations and their weights that fit the SAXS profile starting from a single-input structure. The third protocol is for protein-protein docking based on the SAXS profile of the complex. We describe the underlying software, followed by demonstrating their application on interleukin 33 (IL33) with its primary receptor ST2 and DNA ligase IV-XRCC4 complex.

Hybrid Methods for Modeling Protein Structures Using Molecular Dynamics Simulations and Small-Angle X-Ray Scattering Data

Small-angle X-ray scattering (SAXS) is an efficient experimental tool to measure the overall shape of macromolecular structures in solution. However, due to the low resolution of SAXS data, high-resolution data obtained from X-ray crystallography or NMR and computational methods such as molecular dynamics (MD) simulations are complementary to SAXS data for understanding protein functions based on their structures at atomic resolution. Because MD simulations provide a physicochemically proper structural ensemble for flexible proteins in solution and a precise description of solvent effects, the hybrid analysis of SAXS and MD simulations is a promising method to estimate reasonable solution structures and structural ensembles in solution. Here, we review typical and useful in silico methods for modeling three dimensional protein structures, calculating theoretical SAXS profiles, and analyzing ensemble structures consistent with experimental SAXS profiles. We also review two examples of the hybrid analysis, termed MD-SAXS method in which MD simulations are carried out without any knowledge of experimental SAXS data, and the experimental SAXS data are used only to assess the consistency of the solution model from MD simulations with those observed in experiments. One example is an investigation of the intrinsic dynamics of EcoO109I using the computational method to obtain a theoretical profile from the trajectory of an MD simulation. The other example is a structural investigation of the vitamin D receptor ligand-binding domain using snapshots generated by MD simulations and assessment of the snapshots by experimental SAXS data.

Small Angle Scattering and Structural Biology: Data Quality and Model Validation

This chapter provides a brief review of the current state-of-the-art in small-angle scattering (SAS) from biomolecules in solution in regard to: (1) sample preparation and instrumentation, (2) data reduction and analysis, and (3) three-dimensional structural modelling and validation. In this context, areas of ongoing research in regard to the interpretation of SAS data will be discussed with a particular focus on structural modelling using computational methods and data from different experimental techniques, including SAS (hybrid methods). Finally, progress made in establishing community accepted publication guidelines and a standard reporting framework that includes SAS data deposition in a public data bank will be described. Importantly, SAS data with associated meta-data can now be held in a format that supports exchange between data archives and seamless interoperability with the world-wide Protein Data Bank (wwPDB). Biomolecular SAS is thus well positioned to contribute to an envisioned federation of data archives in support of hybrid structural biology.

Orchestrated Domain Movement in Catalysis by Cytochrome P450 Reductase

NADPH-cytochrome P450 reductase is a multi-domain redox enzyme which is a key component of the P450 mono-oxygenase drug-metabolizing system. We report studies of the conformational equilibrium of this enzyme using small-angle neutron scattering, under conditions where we are able to control the redox state of the enzyme precisely. Different redox states have a profound effect on domain orientation in the enzyme and we analyse the data in terms of a two-state equilibrium between compact and extended conformations. The effects of ionic strength show that the presence of a greater proportion of the extended form leads to an enhanced ability to transfer electrons to cytochrome c. Domain motion is intrinsically linked to the functionality of the enzyme, and we can define the position of the conformational equilibrium for individual steps in the catalytic cycle.

Structural Mechanism for Modulation of Synaptic Neuroligin-Neurexin Signaling by MDGA Proteins

Neuroligin-neurexin (NL-NRX) complexes are fundamental synaptic organizers in the central nervous system. An accurate spatial and temporal control of NL-NRX signaling is crucial to balance excitatory and inhibitory neurotransmission, and perturbations are linked with neurodevelopmental and psychiatric disorders. MDGA proteins bind NLs and control their function and interaction with NRXs via unknown mechanisms. Here, we report crystal structures of MDGA1, the NL1-MDGA1 complex, and a spliced NL1 isoform. Two large, multi-domain MDGA molecules fold into rigid triangular structures, cradling a dimeric NL to prevent NRX binding. Structural analyses guided the discovery of a broad, splicing-modulated interaction network between MDGA and NL family members and helped rationalize the impact of autism-linked mutations. We demonstrate that expression levels largely determine whether MDGAs act selectively or suppress the synapse organizing function of multiple NLs. These results illustrate a potentially brain-wide regulatory mechanism for NL-NRX signaling modulation.

The structural basis for the functional comparability of factor VIII and the long-acting variant recombinant factor VIII Fc fusion protein

Essentials Recombinant factor VIII (rFVIII) Fc fusion protein has a 1.5-fold longer half-life than rFVIII. Five orthogonal methods were used to characterize the structure of rFVIIIFc compared to rFVIII. The C-terminal Fc fusion does not perturb the structure of FVIII in rFVIIIFc. The FVIII and Fc components of rFVIIIFc are flexibly tethered and functionally independent.

SUMMARY

Background Fusion of the human IgG1 Fc domain to the C-terminal C2 domain of B-domain-deleted (BDD) factor VIII (FVIII) results in the recombinant FVIII Fc (rFVIIIFc) fusion protein, which has a 1.5-fold longer half-life in humans. Objective To assess the structural properties of rFVIIIFc by comparing its constituent FVIII and Fc elements with their respective isolated components, and evaluating their structural independence within rFVIIIFc. Methods rFVIIIFc and its isolated FVIII and Fc components were compared by the use of hydrogen-deuterium exchange mass spectrometry (HDX-MS). The structure of rFVIIIFc was also evaluated by the use of X-ray crystallography, small-angle X-ray scattering (SAXS), and electron microscopy (EM). The degree of steric interference by the appended Fc domain was assessed by EM and surface plasmon resonance (SPR). Results HDX-MS analysis of rFVIIIFc revealed that fusion caused no structural perturbations in FVIII or Fc. The rFVIIIFc crystal structure showed that the FVIII component is indistinguishable from published BDD FVIII structures. The Fc domain was not observed, indicating high mobility. SAXS analysis was consistent with an ensemble of rigid-body models in which the Fc domain exists in a largely extended orientation relative to FVIII. Binding of Fab fragments of anti-C2 domain antibodies to BDD FVIII was visualized by EM, and the affinities of the corresponding intact antibodies for BDD FVIII and rFVIIIFc were comparable by SPR analysis. Conclusions The FVIII and Fc components of rFVIIIFc are structurally indistinguishable from their isolated constituents, and show a high degree of structural independence, consistent with the functional comparability of rFVIIIFc and unmodified FVIII.

Analysis of Functional Dynamics of Modular Multidomain Proteins by SAXS and NMR

Multiprotein machines drive virtually all primary cellular processes. Modular multidomain proteins are widely distributed within these dynamic complexes because they provide the flexibility needed to remodel structure as well as rapidly assemble and disassemble components of the machinery. Understanding the functional dynamics of modular multidomain proteins is a major challenge confronting structural biology today because their structure is not fixed in time. Small-angle X-ray scattering (SAXS) and nuclear magnetic resonance (NMR) spectroscopy have proven particularly useful for the analysis of the structural dynamics of modular multidomain proteins because they provide highly complementary information for characterizing the architectural landscape accessible to these proteins. SAXS provides a global snapshot of all architectural space sampled by a molecule in solution. Furthermore, SAXS is sensitive to conformational changes, organization and oligomeric states of protein assemblies, and the existence of flexibility between globular domains in multiprotein complexes. The power of NMR to characterize dynamics provides uniquely complementary information to the global snapshot of the architectural ensemble provided by SAXS because it can directly measure domain motion. In particular, NMR parameters can be used to define the diffusion of domains within modular multidomain proteins, connecting the amplitude of interdomain motion to the architectural ensemble derived from SAXS. Our laboratory has been studying the roles of modular multidomain proteins involved in human DNA replication using SAXS and NMR. Here, we present the procedure for acquiring and analyzing SAXS and NMR data, using DNA primase and replication protein A as examples.

Biochemical Characterization and Low-Resolution SAXS Molecular Envelope of GH1 β-Glycosidase from Saccharophagus degradans

The marine bacteria Saccharophagus degradans (also known as Microbulbifer degradans), are rod-shaped and gram-negative motile γ-proteobacteria, capable of both degrading a variety of complex polysaccharides and fermenting monosaccharides into ethanol. In order to obtain insights into structure-function relationships of the enzymes, involved in these biochemical processes, we characterized a S. degradans β-glycosidase from glycoside hydrolase family 1 (SdBgl1B). SdBgl1B has the optimum pH of 6.0 and a melting temperature T _m of approximately 50 °C. The enzyme has high specificity toward short D-glucose saccharides with β-linkages with the following preferences β-1,3 > β-1,4 ≫ β-1,6. The enzyme kinetic parameters, obtained using artificial substrates p-β-NPGlu and p-β-NPFuc and also the disaccharides cellobiose, gentiobiose and laminaribiose, revealed SdBgl1B preference for p-β-NPGlu and laminaribiose, which indicates its affinity for glucose and also preference for β-1,3 linkages. To better understand structural basis of the enzyme activity its 3D model was built and analysed. The 3D model fits well into the experimentally retrieved low-resolution SAXS-based envelope of the enzyme, confirming monomeric state of SdBgl1B in solution.