A hybrid NMR/SAXS-based approach for discriminating oligomeric protein interfaces using Rosetta

Recent Advances in NMR Protein Structure Prediction with ROSETTA

Nuclear magnetic resonance (NMR) spectroscopy is a powerful method for studying the structure and dynamics of proteins in their native state. For high-resolution NMR structure determination, the collection of a rich restraint dataset is necessary. This can be difficult to achieve for proteins with high molecular weight or a complex architecture. Computational modeling techniques can complement sparse NMR datasets (<1 restraint per residue) with additional structural information to elucidate protein structures in these difficult cases. The Rosetta software for protein structure modeling and design is used by structural biologists for structure determination tasks in which limited experimental data is available. This review gives an overview of the computational protocols available in the Rosetta framework for modeling protein structures from NMR data. We explain the computational algorithms used for the integration of different NMR data types in Rosetta. We also highlight new developments, including modeling tools for data from paramagnetic NMR and hydrogen-deuterium exchange, as well as chemical shifts in CS-Rosetta. Furthermore, strategies are discussed to complement and improve structure predictions made by the current state-of-the-art AlphaFold2 program using NMR-guided Rosetta modeling.

Lack of evidence that Pseudomonas aeruginosa AmpDh3-PA0808 constitute a type VI secretion system effector-immunity pair

Type VI secretion systems (T6SSs) are cell envelope-spanning protein complexes that Gram-negative bacteria use to inject a diverse arsenal of antibacterial toxins into competitor cells. Recently, Wang et al. reported that the H2-T6SS of Pseudomonas aeruginosa delivers the peptidoglycan recycling amidase, AmpDh3, into the periplasm of recipient cells where it is proposed to act as a peptidoglycan degrading toxin. They further reported that PA0808, the open reading frame downstream of AmpDh3, encodes an immunity protein that localizes to the periplasm where it binds to and inactivates intercellularly delivered AmpDh3, thus protecting against its toxic activity. Given that AmpDh3 has an established role in cell wall homeostasis and that no precedent exists for cytosolic enzymes moonlighting as T6SS effectors, we attempted to replicate these findings. We found that cells lacking PA0808 are not susceptible to bacterial killing by AmpDh3 and that PA0808 and AmpDh3 do not physically interact in vitro or in vivo. Additionally, we found no evidence that AmpDh3 is exported from cells, including by strains with a constitutively active H2-T6SS. Finally, subcellular fractionation experiments and a 1.97 Å crystal structure reveal that PA0808 does not contain a canonical signal peptide or localize to the correct cellular compartment to confer protection against a cell wall targeting toxin. Taken together, these results cast doubt on the assertion that AmpDh3-PA0808 constitutes an H2-T6SS effector-immunity pair.

Conjuring up a ghost: structural and functional characterization of FhuF, a ferric siderophore reductase from E. coli

Iron is a fundamental element for virtually all forms of life. Despite its abundance, its bioavailability is limited, and thus, microbes developed siderophores, small molecules, which are synthesized inside the cell and then released outside for iron scavenging. Once inside the cell, iron removal does not occur spontaneously, instead this process is mediated by siderophore-interacting proteins (SIP) and/or by ferric-siderophore reductases (FSR). In the past two decades, representatives of the SIP subfamily have been structurally and biochemically characterized; however, the same was not achieved for the FSR subfamily. Here, we initiate the structural and functional characterization of FhuF, the first and only FSR ever isolated. FhuF is a globular monomeric protein mainly composed by α-helices sheltering internal cavities in a fold resembling the "palm" domain found in siderophore biosynthetic enzymes. Paramagnetic NMR spectroscopy revealed that the core of the cluster has electronic properties in line with those of previously characterized 2Fe-2S ferredoxins and differences appear to be confined to the coordination of Fe(III) in the reduced protein. In particular, the two cysteines coordinating this iron appear to have substantially different bond strengths. In similarity with the proteins from the SIP subfamily, FhuF binds both the iron-loaded and the apo forms of ferrichrome in the micromolar range and cyclic voltammetry reveals the presence of redox-Bohr effect, which broadens the range of ferric-siderophore substrates that can be thermodynamically accessible for reduction. This study suggests that despite the structural differences between FSR and SIP proteins, mechanistic similarities exist between the two classes of proteins.

Macromolecular modeling and design in Rosetta: recent methods and frameworks

The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at http://www.rosettacommons.org.

Integrative Protein Modeling in RosettaNMR from Sparse Paramagnetic Restraints

Computational methods to predict protein structure from nuclear magnetic resonance (NMR) restraints that only require assignment of backbone signals, hold great potential to study larger proteins. Ideally, computational methods designed to work with sparse data need to add atomic detail that is missing in the experimental restraints. We introduce a comprehensive framework into the Rosetta suite that uses NMR restraints derived from paramagnetic labeling. Specifically, RosettaNMR incorporates pseudocontact shifts, residual dipolar couplings, and paramagnetic relaxation enhancements. It continues to use backbone chemical shifts and nuclear Overhauser effect distance restraints. We assess RosettaNMR for protein structure prediction by folding 28 monomeric proteins and 8 homo-oligomeric proteins. Furthermore, the general applicability of RosettaNMR is demonstrated on two protein-protein and three protein-ligand docking examples. Paramagnetic restraints generated more accurate models for 85% of the benchmark proteins and, when combined with chemical shifts, sampled high-accuracy models (≤2Å) in 50% of the cases.

Fast NMR method to probe solvent accessibility and disordered regions in proteins

Understanding protein structure and dynamics, which govern key cellular processes, is crucial for basic and applied research. Intrinsically disordered protein (IDP) regions display multifunctionality via alternative transient conformations, being key players in disease mechanisms. IDP regions are abundant, namely in small viruses, allowing a large number of functions out of a small proteome. The relation between protein function and structure is thus now seen from a different perspective: as IDP regions enable transient structural arrangements, each conformer can play different roles within the cell. However, as IDP regions are hard and time-consuming to study via classical techniques (optimized for globular proteins with unique conformations), new methods are required. Here, employing the dengue virus (DENV) capsid (C) protein and the immunoglobulin-binding domain of streptococcal protein G, we describe a straightforward NMR method to differentiate the solvent accessibility of single amino acid N-H groups in structured and IDP regions. We also gain insights into DENV C flexible fold region biological activity. The method, based on minimal pH changes, uses the well-established ¹H-¹⁵N HSQC pulse sequence and is easily implementable in current protein NMR routines. The data generated are simple to interpret, with this rapid approach being an useful first-choice IDPs characterization method.

The role of small-angle scattering in structure-based screening applications

In many biomolecular interactions, changes in the assembly states and structural conformations of participants can act as a complementary reporter of binding to functional and thermodynamic assays. This structural information is captured by a number of structural biology and biophysical techniques that are viable either as primary screens in small-scale applications or as secondary screens to complement higher throughput methods. In particular, small-angle X-ray scattering (SAXS) reports the average distance distribution between all atoms after orientational averaging. Such information is important when for example investigating conformational changes involved in inhibitory and regulatory mechanisms where binding events do not necessarily cause functional changes. Thus, we summarise here the current and prospective capabilities of SAXS-based screening in the context of other methods that yield structural information. Broad guidelines are also provided to assist readers in preparing screening protocols that are tailored to available X-ray sources.

Rosetta Protein Structure Prediction from Hydroxyl Radical Protein Footprinting Mass Spectrometry Data

In recent years mass spectrometry-based covalent labeling techniques such as hydroxyl radical footprinting (HRF) have emerged as valuable structural biology techniques, yielding information on protein tertiary structure. These data, however, are not sufficient to predict protein structure unambiguously, as they provide information only on the relative solvent exposure of certain residues. Despite some recent advances, no software currently exists that can utilize covalent labeling mass spectrometry data to predict protein tertiary structure. We have developed the first such tool, which incorporates mass spectrometry derived protection factors from HRF labeling as a new centroid score term for the Rosetta scoring function to improve the prediction of protein tertiary structures. We tested our method on a set of four soluble benchmark proteins with known crystal structures and either published HRF experimental results or internally acquired data. Using the HRF labeling data, we rescored large decoy sets of structures predicted with Rosetta for each of the four benchmark proteins. As a result, the model quality improved for all benchmark proteins as compared to when scored with Rosetta alone. For two of the four proteins we were even able to identify atomic resolution models with the addition of HRF data.

Solution structure of CXCL13 and heparan sulfate binding show that GAG binding site and cellular signalling rely on distinct domains

Chemokines promote directional cell migration through binding to G-protein-coupled receptors, and as such are involved in a large array of developmental, homeostatic and pathological processes. They also interact with heparan sulfate (HS), the functional consequences of which depend on the respective location of the receptor- and the HS-binding sites, a detail that remains elusive for most chemokines. Here, to set up a biochemical framework to investigate how HS can regulate CXCL13 activity, we solved the solution structure of CXCL13. We showed that it comprises an unusually long and disordered C-terminal domain, appended to a classical chemokine-like structure. Using three independent experimental approaches, we found that it displays a unique association mode to HS, involving two clusters located in the α-helix and the C-terminal domain. Computational approaches were used to analyse the HS sequences preferentially recognized by the protein and gain atomic-level understanding of the CXCL13 dimerization induced upon HS binding. Starting with four sets of 254 HS tetrasaccharides, we identified 25 sequences that bind to CXCL13 monomer, among which a single one bound to CXCL13 dimer with high consistency. Importantly, we found that CXCL13 can be functionally presented to its receptor in a HS-bound form, suggesting that it can promote adhesion-dependent cell migration. Consistently, we designed CXCL13 mutations that preclude interaction with HS without affecting CXCR5-dependent cell signalling, opening the possibility to unambiguously demonstrate the role of HS in the biological function of this chemokine.

Chemical shift-based methods in NMR structure determination

Chemical shifts are highly sensitive probes harnessed by NMR spectroscopists and structural biologists as conformational parameters to characterize a range of biological molecules. Traditionally, assignment of chemical shifts has been a labor-intensive process requiring numerous samples and a suite of multidimensional experiments. Over the past two decades, the development of complementary computational approaches has bolstered the analysis, interpretation and utilization of chemical shifts for elucidation of high resolution protein and nucleic acid structures. Here, we review the development and application of chemical shift-based methods for structure determination with a focus on ab initio fragment assembly, comparative modeling, oligomeric systems, and automated assignment methods. Throughout our discussion, we point out practical uses, as well as advantages and caveats, of using chemical shifts in structure modeling. We additionally highlight (i) hybrid methods that employ chemical shifts with other types of NMR restraints (residual dipolar couplings, paramagnetic relaxation enhancements and pseudocontact shifts) that allow for improved accuracy and resolution of generated 3D structures, (ii) the utilization of chemical shifts to model the structures of sparsely populated excited states, and (iii) modeling of sidechain conformations. Finally, we briefly discuss the advantages of contemporary methods that employ sparse NMR data recorded using site-specific isotope labeling schemes for chemical shift-driven structure determination of larger molecules. With this review, we aim to emphasize the accessibility and versatility of chemical shifts for structure determination of challenging biological systems, and to point out emerging areas of development that lead us towards the next generation of tools.

Enhancing Protein Conformational Space Sampling Using Distance Profile-Guided Differential Evolution

De novo protein structure prediction aims to search for low-energy conformations as it follows the thermodynamics hypothesis that places native conformations at the global minimum of the protein energy surface. However, the native conformation is not necessarily located in the lowest-energy regions owing to the inaccuracies of the energy model. This study presents a differential evolution algorithm using distance profile-based selection strategy to sample conformations with reasonable structure effectively. In the proposed algorithm, besides energy, the residue-residue distance is considered another measure of the conformation. The average distance errors of decoys between the distance of each residue pair and the corresponding distance in the distance profiles are first calculated when the trial conformation yields a larger energy value than that of the target. Then, the distance acceptance probability of the trial conformation is designed based on distance profiles if the trial conformation obtains a lower average distance error compared with that of the target conformation. The trial conformation is accepted to the next generation in accordance with its distance acceptance probability. By using the dual constraints of energy and distance in guiding sampling, the algorithm can sample conformations with lower energies and more reasonable structures. Experimental results of 28 benchmark proteins show that the proposed algorithm can effectively predict near-native protein structures.

Combining NMR and small angle X-ray scattering for the study of biomolecular structure and dynamics

Small-angle X-ray scattering (SAXS) and Nuclear Magnetic Resonance (NMR) are established methods to analyze the structure and structural transitions of biological macromolecules in solution. Both methods are directly applicable to near-native macromolecular solutions and allow one to study structural responses to physical and chemical changes or ligand additions. Whereas SAXS is applied to elucidate overall structure, interactions and flexibility over a wide range of particle sizes, NMR yields atomic resolution detail for moderately sized macromolecules. NMR is arguably the most powerful technique for the experimental analysis of dynamics. The joint application of these two highly complementary techniques provides an extremely useful approach that facilitates comprehensive characterization of biomacromolecular solutions.

•

SAXS and NMR are effective and highly complementary techniques in structural biology.

•

Constraints from SAXS can be readily incorporated in NMR structure calculations.

•

High resolution NMR models of domains can serve as building blocks for SAXS-based rigid body modeling.

•

Flexible systems can be well described using ensemble approaches combining SAXS and NMR.

•

Dynamics studies can be enhanced by combining SAXS and NMR.

Hybrid Applications of Solution Scattering to Aid Structural Biology

Biomolecular applications of solution X-ray and neutron scattering (SAXS and SANS, respectively) started in late 1960s - early 1970s but were relatively limited in their ability to provide a detailed structural picture and lagged behind what became the two primary methods of experimental structural biology - X-ray crystallography and NMR. However, improvements in both data analysis and instrumentation led to an explosive growth in the number of studies that used small-angle scattering (SAS) for investigation of macromolecular structure, often in combination with other biophysical techniques. Such hybrid applications are nowadays quickly becoming a norm whenever scattering data are used for two reasons. First, it is generally accepted that SAS data on their own cannot lead to a uniquely defined high-resolution structural model, creating a need for supplementing them with information from complementary techniques. Second, solution scattering data are frequently applied in situations when a method such NMR or X-ray crystallography cannot provide a satisfactory structural picture, which makes these additional restraints highly desirable. Maturation of the hybrid bio-SAS approaches brings to light new questions including completeness of the conformational space sampling, model validation, and data compatibility.

15N and 13C- SOFAST-HMQC editing enhances 3D-NOESY sensitivity in highly deuterated, selectively [1H,13C]-labeled proteins

The ongoing NMR method development effort strives for high quality multidimensional data with reduced collection time. Here, we apply 'SOFAST-HMQC' to frequency editing in 3D NOESY experiments and demonstrate the sensitivity benefits using highly deuterated and ¹⁵N, methyl labeled samples in H₂O. The experiments benefit from a combination of selective T ₁ relaxation (or L-optimized effect), from Ernst angle optimization and, in certain types of experiments, from using the mixing time for both NOE buildup and magnetization recovery. This effect enhances sensitivity by up to 2.4× at fast pulsing versus reference HMQC sequences of same overall length and water suppression characteristics. Representative experiments designed to address interesting protein NMR challenges are detailed. Editing capabilities are exploited with heteronuclear ¹⁵N,¹³C-edited, or with diagonal-free ¹³C aromatic/methyl-resolved 3D-SOFAST-HMQC-NOESY-HMQC. The latter experiment is used here to elucidate the methyl-aromatic NOE network in the hydrophobic core of the 19 kDa FliT-FliJ flagellar protein complex. Incorporation of fast pulsing to reference experiments such as 3D-NOESY-HMQC boosts digital resolution, simplifies the process of NOE assignment and helps to automate protein structure determination.

Protocols for Molecular Modeling with Rosetta3 and RosettaScripts

Previously, we published an article providing an overview of the Rosetta suite of biomacromolecular modeling software and a series of step-by-step tutorials [Kaufmann, K. W., et al. (2010) Biochemistry 49, 2987–2998]. The overwhelming positive response to this publication we received motivates us to here share the next iteration of these tutorials that feature de novo folding, comparative modeling, loop construction, protein docking, small molecule docking, and protein design. This updated and expanded set of tutorials is needed, as since 2010 Rosetta has been fully redesigned into an object-oriented protein modeling program Rosetta3. Notable improvements include a substantially improved energy function, an XML-like language termed “RosettaScripts” for flexibly specifying modeling task, new analysis tools, the addition of the TopologyBroker to control conformational sampling, and support for multiple templates in comparative modeling. Rosetta’s ability to model systems with symmetric proteins, membrane proteins, noncanonical amino acids, and RNA has also been greatly expanded and improved.

Recognition and targeting mechanisms by chaperones in flagellum assembly and operation

The flagellum is a complex bacterial nanomachine that requires the proper assembly of several different proteins for its function. Dedicated chaperones are central in preventing aggregation or undesired interactions of flagellar proteins, including their targeting to the export gate. FliT is a key flagellar chaperone that binds to several flagellar proteins in the cytoplasm, including its cognate filament-capping protein FliD. We have determined the solution structure of the FliT chaperone in the free state and in complex with FliD and the flagellar ATPase FliI. FliT adopts a four-helix bundle and uses a hydrophobic surface formed by the first three helices to recognize its substrate proteins. We show that the fourth helix constitutes the binding site for FlhA, a membrane protein at the export gate. In the absence of a substrate protein FliT adopts an autoinhibited structure wherein both the binding sites for substrates and FlhA are occluded. Substrate binding to FliT activates the complex for FlhA binding and thus targeting of the chaperone-substrate complex to the export gate. The activation and targeting mechanisms reported for FliT appear to be shared among the other flagellar chaperones.

Hybrid Approaches to Structural Characterization of Conformational Ensembles of Complex Macromolecular Systems Combining NMR Residual Dipolar Couplings and Solution X-ray Scattering

Solving structures or structural ensembles of large macromolecular systems in solution poses a challenging problem. While NMR provides structural information at atomic resolution, increased spectral complexity, chemical shift overlap, and short transverse relaxation times (associated with slow tumbling) render application of the usual techniques that have been so successful for medium sized systems (<50 kDa) difficult. Solution X-ray scattering, on the other hand, is not limited by molecular weight but only provides low resolution structural information related to the overall shape and size of the system under investigation. Here we review how combining atomic resolution structures of smaller domains with sparse experimental data afforded by NMR residual dipolar couplings (which yield both orientational and shape information) and solution X-ray scattering data in rigid-body simulated annealing calculations provides a powerful approach for investigating the structural aspects of conformational dynamics in large multidomain proteins. The application of this hybrid methodology is illustrated for the 128 kDa dimer of bacterial Enzyme I which exists in a variety of open and closed states that are sampled at various points in the catalytic cycles, and for the capsid protein of the human immunodeficiency virus.

Protein-protein interactions: a supra-structural phenomenon demanding trans-disciplinary biophysical approaches

Responsive formation of protein:protein interaction (PPI) upon diverse stimuli is a fundament of cellular function. As a consequence, PPIs are complex, adaptive entities, and exist in structurally heterogeneous interplays defined by the energetic states of the free and complexed protomers. The biophysical and structural investigations of PPIs consequently demand hybrid approaches, implementing orthogonal methods and strategies for global data analysis. Currently, impressive developments in hardware and software within several methodologies define a new era for the biostructural community. Data can be obtained at increasing resolution, at relevant time-scales and under increasingly relevant experimental conditions, intricate data are interpreted reliably, and the questions posed and answered grow in complexity. With this review, highlights from the study of PPIs using a multitude of biophysical methods, are reported. The aim is to depict how the elucidation of the interplay of structures requires the interplay of methods.

A Novel MHC-I Surface Targeted for Binding by the MCMV m06 Immunoevasin Revealed by Solution NMR

As part of its strategy to evade detection by the host immune system, murine cytomegalovirus (MCMV) encodes three proteins that modulate cell surface expression of major histocompatibility complex class I (MHC-I) molecules: the MHC-I homolog m152/gp40 as well as the m02-m16 family members m04/gp34 and m06/gp48. Previous studies of the m04 protein revealed a divergent Ig-like fold that is unique to immunoevasins of the m02-m16 family. Here, we engineer and characterize recombinant m06 and investigate its interactions with full-length and truncated forms of the MHC-I molecule H2-L(d) by several techniques. Furthermore, we employ solution NMR to map the interaction footprint of the m06 protein on MHC-I, taking advantage of a truncated H2-L(d), "mini-H2-L(d)," consisting of only the α1α2 platform domain. Mini-H2-L(d) refolded in vitro with a high affinity peptide yields a molecule that shows outstanding NMR spectral features, permitting complete backbone assignments. These NMR-based studies reveal that m06 binds tightly to a discrete site located under the peptide-binding platform that partially overlaps with the β2-microglobulin interface on the MHC-I heavy chain, consistent with in vitro binding experiments showing significantly reduced complex formation between m06 and β2-microglobulin-associated MHC-I. Moreover, we carry out NMR relaxation experiments to characterize the picosecond-nanosecond dynamics of the free mini-H2-L(d) MHC-I molecule, revealing that the site of interaction is highly ordered. This study provides insight into the mechanism of the interaction of m06 with MHC-I, suggesting a structural manipulation of the target MHC-I molecule at an early stage of the peptide-loading pathway.

A community resource of experimental data for NMR / X-ray crystal structure pairs

We have developed an online NMR / X-ray Structure Pair Data Repository. The NIGMS Protein Structure Initiative (PSI) has provided many valuable reagents, 3D structures, and technologies for structural biology. The Northeast Structural Genomics Consortium was one of several PSI centers. NESG used both X-ray crystallography and NMR spectroscopy for protein structure determination. A key goal of the PSI was to provide experimental structures for at least one representative of each of hundreds of targeted protein domain families. In some cases, structures for identical (or nearly identical) constructs were determined by both NMR and X-ray crystallography. NMR spectroscopy and X-ray diffraction data for 41 of these "NMR / X-ray" structure pairs determined using conventional triple-resonance NMR methods with extensive sidechain resonance assignments have been organized in an online NMR / X-ray Structure Pair Data Repository. In addition, several NMR data sets for perdeuterated, methyl-protonated protein samples are included in this repository. As an example of the utility of this repository, these data were used to revisit questions about the precision and accuracy of protein NMR structures first outlined by Levy and coworkers several years ago (Andrec et al., Proteins 2007;69:449-465). These results demonstrate that the agreement between NMR and X-ray crystal structures is improved using modern methods of protein NMR spectroscopy. The NMR / X-ray Structure Pair Data Repository will provide a valuable resource for new computational NMR methods development.