Integrating NMR, SAXS, and Atomistic Simulations: Structure and Dynamics of a Two-Domain Protein

Structural Basis of Non-Latent Signaling by the Anti-Müllerian Hormone Procomplex

Most TGFβ family ligands exist as procomplexes consisting of a prodomain noncovalently bound to a growth factor (GF); Whereas some prodomains confer latency, the Anti-Müllerian Hormone (AMH) prodomain maintains a remarkably high affinity for the GF yet remains active. Using single particle EM methods, we show the AMH prodomain consists of two subdomains: a vestigial TGFβ prodomain-like fold and a novel, helical bundle GF-binding domain, the result of an exon insertion 450 million years ago, that engages both receptor epitopes. When associated with the prodomain, the AMH GF is distorted into a strained, open conformation whose closure upon bivalent binding of AMHR2 displaces the prodomain through a conformational shift mechanism to allow for signaling.

Small angle X-ray scattering analysis of thermophilic cytochrome P450 CYP119 and the effects of the N-terminal histidine tag

Combining size exclusion chromatography-small angle X-ray scattering (SEC-SAXS) and molecular dynamics (MD) analysis is a promising approach to investigate protein behavior in solution, particularly for understanding conformational changes due to substrate binding in cytochrome P450s (CYPs). This study investigates conformational changes in CYP119, a thermophilic CYP from Sulfolobus acidocaldarius that exhibits structural flexibility similar to mammalian CYPs. Although the crystal structure of ligand-free (open state) and ligand-bound (closed state) forms of CYP119 is known, the overall structure of the enzyme in solution has not been explored until now. It was found that theoretical scattering profiles from the crystal structures of CYP119 did not align with the SAXS data, but conformers from MD simulations, particularly starting from the open state (46 % of all frames), agreed well. Interestingly, a small percentage of closed-state conformers also fit the data (9 %), suggesting ligand-free CYP119 samples ligand-bound conformations. Ab initio SAXS models for N-His tagged CYP119 revealed a tail-like unfolded structure impacting protein flexibility, which was confirmed by in silico modeling. SEC-SAXS analysis of N-His CYP119 indicated pentameric structures in addition to monomers in solution, affecting the stability and activity of the enzyme. This study adds insights into the conformational dynamics of CYP119 in solution.

Principles Governing the Phase Separation of Multidomain Proteins

A variety of membraneless organelles, often termed "biological condensates", play an important role in the regulation of cellular processes such as gene transcription, translation, and protein quality control. On the basis of experimental and theoretical investigations, liquid-liquid phase separation (LLPS) has been proposed as a possible mechanism for the origin of biological condensates. LLPS requires multivalent macromolecules that template the formation of long-range, intermolecular interaction networks and results in the formation of condensates with defined composition and material properties. Multivalent interactions driving LLPS exhibit a wide range of modes from highly stereospecific to nonspecific and involve both folded and disordered regions. Multidomain proteins serve as suitable macromolecules for promoting phase separation and achieving disparate functions due to their potential for multivalent interactions and regulation. Here, we aim to highlight the influence of the domain architecture and interdomain interactions on the phase separation of multidomain protein condensates. First, the general principles underlying these interactions are illustrated on the basis of examples of multidomain proteins that are predominantly associated with nucleic acid binding and protein quality control and contain both folded and disordered regions. Next, the examples showcase how LLPS properties of folded and disordered regions can be leveraged to engineer multidomain constructs that form condensates with the desired assembly and functional properties. Finally, we highlight the need for improvements in coarse-grained computational models that can provide molecular-level insights into multidomain protein condensates in conjunction with experimental efforts.

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.

Small-angle X-ray scattering for the proteomics community: current overview and future potential

Introduction: Proteins are biological nanoparticles. For structural proteomics and hybrid structural biology, complementary methods are required that allow both high throughput and accurate automated data analysis. Small-angle X-ray scattering (SAXS) is a method for observing the size and shape of particles, such as proteins and complexes, in solution. SAXS data can be used to model both the structure, oligomeric state, conformational changes, and flexibility of biomolecular samples.Areas covered: The key principles of SAXS, its sample requirements, and its current and future applications for structural proteomics are briefly reviewed. Recent technical developments in SAXS experiments are discussed, and future potential of the method in structural proteomics is evaluated.Expert opinion: SAXS is a method suitable for several aspects of integrative structural proteomics, with current technical developments allowing for higher throughput and time-resolved studies, as well as the analysis of complex samples, such as membrane proteins. Increasing automation and streamlined data analysis are expected to equip SAXS for structure-based screening workflows. Originally, structural genomics had a heavy focus on folded, crystallizable proteins and complexes - SAXS is a method allowing an expansion of this focus to flexible and disordered systems.

NMR-Based Methods for Protein Analysis

Nuclear magnetic resonance (NMR) spectroscopy is a well-established method for analyzing protein structure, interaction, and dynamics at atomic resolution and in various sample states including solution state, solid state, and membranous environment. Thanks to rapid NMR methodology development, the past decade has witnessed a growing number of protein NMR studies in complex systems ranging from membrane mimetics to living cells, which pushes the research frontier further toward physiological environments and offers unique insights in elucidating protein functional mechanisms. In particular, in-cell NMR has become a method of choice for bridging the huge gap between structural biology and cell biology. Herein, we review the recent developments and applications of NMR methods for protein analysis in close-to-physiological environments, with special emphasis on in-cell protein structural determination and the analysis of protein dynamics, both difficult to be accessed by traditional methods.

Quantitative comparison between sub-millisecond time resolution single-molecule FRET measurements and 10-second molecular simulations of a biosensor protein

Molecular Dynamics (MD) simulations seek to provide atomic-level insights into conformationally dynamic biological systems at experimentally relevant time resolutions, such as those afforded by single-molecule fluorescence measurements. However, limitations in the time scales of MD simulations and the time resolution of single-molecule measurements have challenged efforts to obtain overlapping temporal regimes required for close quantitative comparisons. Achieving such overlap has the potential to provide novel theories, hypotheses, and interpretations that can inform idealized experimental designs that maximize the detection of the desired reaction coordinate. Here, we report MD simulations at time scales overlapping with in vitro single-molecule Förster (fluorescence) resonance energy transfer (smFRET) measurements of the amino acid binding protein LIV-BP^SS at sub-millisecond resolution. Computationally efficient all-atom structure-based simulations, calibrated against explicit solvent simulations, were employed for sampling multiple cycles of LIV-BP^SS clamshell-like conformational changes on the time scale of seconds, examining the relationship between these events and those observed by smFRET. The MD simulations agree with the smFRET measurements and provide valuable information on local dynamics of fluorophores at their sites of attachment on LIV-BP^SS and the correlations between fluorophore motions and large-scale conformational changes between LIV-BP^SS domains. We further utilize the MD simulations to inform the interpretation of smFRET data, including Förster radius (R₀) and fluorophore orientation factor (κ²) determinations. The approach we describe can be readily extended to distinct biochemical systems, allowing for the interpretation of any FRET system conjugated to protein or ribonucleoprotein complexes, including those with more conformational processes, as well as those implementing multi-color smFRET.

Förster (fluorescence) resonance energy transfer (FRET) has been used extensively by biophysicists as a molecular-scale ruler that yields fundamental structural and kinetic insights into transient processes including complex formation and conformational rearrangements required for biological function. FRET techniques require the identification of informative fluorophore labeling sites, spaced at defined distances to inform on a reaction coordinate of interest and consideration of noise sources that have the potential to obscure quantitative interpretations. Here, we describe an approach to leverage advancements in computationally efficient all-atom structure-based molecular dynamics simulations in which structural dynamics observed via FRET can be interpreted in full atomistic detail on commensurate time scales. We demonstrate the potential of this approach using a model FRET system, the amino acid binding protein LIV-BP^SS conjugated to self-healing organic fluorophores. LIV-BP^SS exhibits large scale, sub-millisecond clamshell-like conformational changes between open and closed conformations associated with ligand unbinding and binding, respectively. Our findings inform on the molecular basis of the dynamics observed by smFRET and on strategies to optimize fluorophore labeling sites, the manner of fluorophore attachment, and fluorophore composition.

Evolving SAXS versatility: solution X-ray scattering for macromolecular architecture, functional landscapes, and integrative structural biology

Small-angle X-ray scattering (SAXS) has emerged as an enabling integrative technique for comprehensive analyses of macromolecular structures and interactions in solution. Over the past two decades, SAXS has become a mainstay of the structural biologist's toolbox, supplying multiplexed measurements of molecular shape and dynamics that unveil biological function. Here, we discuss evolving SAXS theory, methods, and applications that extend the field of small-angle scattering beyond simple shape characterization. SAXS, coupled with size-exclusion chromatography (SEC-SAXS) and time-resolved (TR-SAXS) methods, is now providing high-resolution insight into macromolecular flexibility and ensembles, delineating biophysical landscapes, and facilitating high-throughput library screening to assess macromolecular properties and to create opportunities for drug discovery. Looking forward, we consider SAXS in the integrative era of hybrid structural biology methods, its potential for illuminating cellular supramolecular and mesoscale structures, and its capacity to complement high-throughput bioinformatics sequencing data. As advances in the field continue, we look forward to proliferating uses of SAXS based upon its abilities to robustly produce mechanistic insights for biology and medicine.

Coarse-grain simulations on NMR conformational ensembles highlight functional residues in proteins

Dynamics are a key feature of protein function, and this is especially true of gating residues, which occupy cavity or tunnel lining positions in the protein structure, and will reversibly switch between open and closed conformations in order to control the diffusion of small molecules within a protein's internal matrix. Earlier work on globins and hydrogenases have shown that these gating residues can be detected using a multiscale scheme combining all-atom classic molecular dynamics simulations and coarse-grain calculations of the resulting conformational ensemble mechanical properties. Here, we show that the structural variations observed in the conformational ensembles produced by NMR spectroscopy experiments are sufficient to induce noticeable mechanical changes in a protein, which in turn can be used to identify residues important for function and forming a mechanical nucleus in the protein core. This new approach, which combines experimental data and rapid coarse-grain calculations and no longer needs to resort to time-consuming all-atom simulations, was successfully applied to five different protein families.

Modified Potential Functions Result in Enhanced Predictions of a Protein Complex by All-Atom Molecular Dynamics Simulations, Confirming a Stepwise Association Process for Native Protein-Protein Interactions

The relative prevalence of native protein-protein interactions (PPIs) are the cornerstone for understanding the structure, dynamics and mechanisms of function of protein complexes. In this study, we develop a scheme for scaling the protein-water interaction in the CHARMM36 force field, in order to better fit the solvation free energy of amino acids side-chain analogues. We find that the molecular dynamics simulation with the scaled force field, CHARMM36s, as well as a recently released version, CHARMM36m, effectively improve on the overly sticky association of proteins, such as ubiquitin. We investigate the formation of a heterodimer protein complex between the SAM domains of the EphA2 receptor and the SHIP2 enzyme by performing a combined total of 48 μs simulations with the different potential functions. While the native SAM heterodimer is only predicted at a low rate of 6.7% with the original CHARMM36 force field, the yield is increased to 16.7% with CHARMM36s, and to 18.3% with CHARMM36m. By analyzing the 25 native SAM complexes formed in the simulations, we find that their formation involves a preorientation guided by Coulomb interactions, consistent with an electrostatic steering mechanism. In 12 cases, the complex could directly transform to the native protein interaction surfaces with only small adjustments in domain orientation. In the other 13 cases, orientational and/or translational adjustments are needed to reach the native complex. Although the tendency for non-native complexes to dissociate has nearly doubled with the modified potential functions, a dissociation followed by a reassociation to the correct complex structure is still rare. Instead, the remaining non-native complexes undergo configurational changes/surface searching, which, however, rarely leads to native structures on a time scale of 250 ns. These observations provide a rich picture of the mechanisms of protein-protein complex formation and suggest that computational predictions of native complex PPIs could be improved further.

Accuracy and precision of protein structures determined by magic angle spinning NMR spectroscopy: for some 'with a little help from a friend'

We present a systematic investigation into the attainable accuracy and precision of protein structures determined by heteronuclear magic angle spinning solid-state NMR for a set of four proteins of varied size and secondary structure content. Structures were calculated using synthetically generated random sets of C-C distances up to 7 Å at different degrees of completeness. For single-domain proteins, 9-15 restraints per residue are sufficient to derive an accurate model structure, while maximum accuracy and precision are reached with over 15 restraints per residue. For multi-domain proteins and protein assemblies, additional information on domain orientations, quaternary structure and/or protein shape is needed. As demonstrated for the HIV-1 capsid protein assembly, this can be accomplished by integrating MAS NMR with cryoEM data. In all cases, inclusion of TALOS-derived backbone torsion angles improves the accuracy for small number of restraints, while no further increases are noted for restraint completeness above 40%. In contrast, inclusion of TALOS-derived torsion angle restraints consistently increases the precision of the structural ensemble at all degrees of distance restraint completeness.

Rotational Dynamics of Proteins from Spin Relaxation Times and Molecular Dynamics Simulations

Conformational fluctuations and rotational tumbling of proteins can be experimentally accessed with nuclear spin relaxation experiments. However, interpretation of molecular dynamics from the experimental data is often complicated, especially for molecules with anisotropic shape. Here, we apply classical molecular dynamics simulations to interpret the conformational fluctuations and rotational tumbling of proteins with arbitrarily anisotropic shape. The direct calculation of spin relaxation times from simulation data did not reproduce the experimental data. This was successfully corrected by scaling the overall rotational diffusion coefficients around the protein inertia axes with a constant factor. The achieved good agreement with experiments allowed the interpretation of the internal and overall dynamics of proteins with significantly anisotropic shape. The overall rotational diffusion was found to be Brownian, having only a short subdiffusive region below 0.12 ns. The presented methodology can be applied to interpret rotational dynamics and conformation fluctuations of proteins with arbitrary anisotropic shape. However, a water model with more realistic dynamical properties is probably required for intrinsically disordered proteins.