High quality NMR structures: a new force field with implicit water and membrane solvation for Xplor-NIH

Strategies for elucidation of the structure and function of the large membrane protein complex, FoF1-ATP synthase, by nuclear magnetic resonance

Nuclear magnetic resonance (NMR) investigation of large membrane proteins requires well-focused questions and critical techniques. Here, research strategies for F_oF₁-ATP synthase, a membrane-embedded molecular motor, are reviewed, focusing on the β-subunit of F₁-ATPase and c-subunit ring of the enzyme. Segmental isotope-labeling provided 89% assignment of the main chain NMR signals of thermophilic Bacillus (T)F₁β-monomer. Upon nucleotide binding to Lys164, Asp252 was shown to switch its hydrogen-bonding partner from Lys164 to Thr165, inducing an open-to-closed bend motion of TF₁β-subunit. This drives the rotational catalysis. The c-ring structure determined by solid-state NMR showed that cGlu56 and cAsn23 of the active site took a hydrogen-bonded closed conformation in membranes. In 505 kDa TF_oF₁, the specifically isotope-labeled cGlu56 and cAsn23 provided well-resolved NMR signals, which revealed that 87% of the residue pairs took a deprotonated open conformation at the F_oa-c subunit interface, whereas they were in the closed conformation in the lipid-enclosed region.

Anionic lipids induce a fold-unfold transition in the membrane-translocating Engrailed homeodomain

Homeoprotein transcription factors have the property of interacting with membranes through their DNA-binding homeodomain, which is involved in unconventional internalization and secretion. Both processes depend on membrane-translocating events but their detailed molecular mechanisms are still poorly understood. We have previously characterized the conformational properties of Engrailed 2 homeodomain (EnHD) in aqueous solution and in micelles as membrane-mimetic environments. In the present study, we used small isotropic lipid bicelles as a more relevant membrane-mimetic model to characterize the membrane-bound state of EnHD. We show that lipid bicelles, in contrast to micelles, adequately reproduce the requirement of anionic lipids in the membrane binding and conformational transition of EnHD. The fold-unfold transition of EnHD induced by anionic lipids was characterized by NMR using ¹H, ¹³C, ¹⁵N chemical shifts, nuclear Overhauser effects, residual dipolar couplings, intramolecular and intermolecular paramagnetic relaxation enhancements induced by site-directed spin-label or paramagnetic lipid probe, respectively. A global unpacking of EnHD helices is observed leading to a loss of the native fold. However, near-native propensities of EnHD backbone conformation are maintained in membrane environment, including not only the three helices but also the turn connecting helices H2 and H3. NMR and coarse-grained molecular dynamics simulations reveal that the EnHD adopts a shallow insertion in the membrane, with the three helices oriented parallel to the membrane. EnHD explores extended conformations and closed U-shaped conformations, which are stabilized by anionic lipid recruitment.

Chemical Conformation of the Essential Glutamate Site of the c-Ring within Thermophilic Bacillus FoF1-ATP Synthase Determined by Solid-State NMR Based on its Isolated c-Ring Structure

Proton translocation through the membrane-embedded F_o component of F-type ATP synthase (F_oF₁) is facilitated by the rotation of the F_o c-subunit ring (c-ring), carrying protons at essential acidic amino acid residues. Cryo-electron microscopy (Cryo-EM) structures of F_oF₁ suggest a unique proton translocation mechanism. To elucidate it based on the chemical conformation of the essential acidic residues of the c-ring in F_oF₁, we determined the structure of the isolated thermophilic Bacillus F_o (tF_o) c-ring, consisting of 10 subunits, in membranes by solid-state NMR. This structure contains a distinct proton-locking conformation, wherein Asn23 (cN23) C^γO and Glu56 (cE56) C^δOH form a hydrogen bond in a closed form. We introduced stereo-array-isotope-labeled (SAIL) Glu and Asn into the tF_oc-ring to clarify the chemical conformation of these residues in tF_oF₁-ATP synthase (tF_oF₁). Two well-separated ¹³C signals could be detected for cN23 and cE56 in a 505 kDa membrane protein complex, respectively, thereby suggesting the presence of two distinct chemical conformations. Based on the signal intensity and structure of the tF_oc-ring and tF_oF₁, six pairs of cN23 and cE56 surrounded by membrane lipids take the closed form, whereas the other four in the a-c interface employ the deprotonated open form at a proportion of 87%. This indicates that the a-c interface is highly hydrophilic. The pK_a values of the four cE56 residues in the a-c interface were estimated from the cN23 signal intensity in the open and closed forms and distribution of polar residues around each cE56. The results favor a rotation of the c-ring for ATP synthesis.

Exploration of the Molecular Mechanism of Danzhi Xiaoyao Powder in Endometrial Cancer through Network Pharmacology

Endometrial cancer (EC) is a common malignant tumor of the female reproductive system. Current treatments such as surgery and long-term hormone therapy are ineffective and have side effects. Danzhi Xiaoyao powder (DXP) can inhibit the growth of EC cells and induce apoptosis, but the pharmacological and molecular mechanisms of anticancer effects are still unclear. In this study, active components and potential targets of DXP were obtained from public databases. Protein effects and regulatory pathways of common targets were analyzed by protein-protein interaction (PPI), GO and KEGG. The results of network pharmacology showed that there are 87 common targets between EC and DXP. GO enrichment analysis showed that these targets were associated with response to oxidative stress, response to nutrient levels, hormone receptor binding and nuclear hormone receptor binding, etc. The results of KEGG analysis indicated that IL-17, TNF, PI3K/AKT, and RAS/RAF/MEK/ERK (ERK) signaling pathway were enriched in the anti-EC of DXP. Additionally, we cultured HEC-1B and KLE cells for validate experiments. DXP showed an inhibition of proliferation, migration, and cell cycle of both cells. Moreover, the expression of RAS, p-RAF, p-MEK, ERK, and p-ERK related proteins were downregulated. In conclusion, DXP might inhibit the proliferation of EC cells via apoptosis. Furthermore, DXP-induced inhibition of EC development might involve RAS/RAF/MEK/ERK pathway.

A Cholesterol Dimer Stabilizes the Inactivated State of an Inward-Rectifier Potassium Channel

Cholesterol oligomers reside in multiple membrane protein X-ray crystal structures. Yet, there is no direct link between these oligomers and a biological function. Here we present the structural and functional details of a cholesterol dimer that stabilizes the inactivated state of an inward-rectifier potassium channel KirBac1.1. K⁺ efflux assays confirm that high cholesterol concentration reduces K⁺ conductance. We then determine the structure of the cholesterol-KirBac1.1 complex using Xplor-NIH simulated annealing calculations driven by solid-state NMR distance measurements. These calculations identified an α-α cholesterol dimer docked to a cleft formed by adjacent subunits of the homotetrameric protein. We compare these results to coarse grain molecular dynamics simulations. This is one of the first examples of a cholesterol oligomer performing a distinct biological function and structural characterization of a conserved promiscuous lipid binding region.

Water Accessibility Refinement of the Extended Structure of KirBac1.1 in the Closed State

NMR structures of membrane proteins are often hampered by poor chemical shift dispersion and internal dynamics which limit resolved distance restraints. However, the ordering and topology of these systems can be defined with site-specific water or lipid proximity. Membrane protein water accessibility surface area is often investigated as a topological function via solid-state NMR. Here we leverage water-edited solid-state NMR measurements in simulated annealing calculations to refine a membrane protein structure. This is demonstrated on the inward rectifier K+ channel KirBac1.1 found in Burkholderia pseudomallei. KirBac1.1 is homologous to human Kir channels, sharing a nearly identical fold. Like many existing Kir channel crystal structures, the 1p7b crystal structure is incomplete, missing 85 out of 333 residues, including the N-terminus and C-terminus. We measure solid-state NMR water proximity information and use this for refinement of KirBac1.1 using the Xplor-NIH structure determination program. Along with predicted dihedral angles and sparse intra- and inter-subunit distances, we refined the residues 1-300 to atomic resolution. All structural quality metrics indicate these restraints are a powerful way forward to solve high quality structures of membrane proteins using NMR.

Molecular basis of anticoagulant and anticomplement activity of the tick salivary protein Salp14 and its homologs

During feeding, a tick's mouthpart penetrates the host's skin and damages tissues and small blood vessels, triggering the extrinsic coagulation and lectin complement pathways. To elude these defense mechanisms, ticks secrete multiple anticoagulant proteins and complement system inhibitors in their saliva. Here, we characterized the inhibitory activities of the homologous tick salivary proteins tick salivary lectin pathway inhibitor, Salp14, and Salp9Pac from Ixodesscapularis in the coagulation cascade and the lectin complement pathway. All three proteins inhibited binding of mannan-binding lectin to the polysaccharide mannan, preventing the activation of the lectin complement pathway. In contrast, only Salp14 showed an appreciable effect on coagulation by prolonging the lag time of thrombin generation. We found that the anticoagulant properties of Salp14 are governed by its basic tail region, which resembles the C terminus of tissue factor pathway inhibitor alpha and blocks the assembly and/or activity of the prothrombinase complex in the same way. Moreover, the Salp14 protein tail contributes to the inhibition of the lectin complement pathway via interaction with mannan binding lectin-associated serine proteases. Furthermore, we identified BaSO₄-adsorbing protein 1 isolated from the tick Ornithodoros savignyi as a distant homolog of tick salivary lectin pathway inhibitor/Salp14 proteins and showed that it inhibits the lectin complement pathway but not coagulation. The structure of BaSO₄-adsorbing protein 1, solved here using NMR spectroscopy, indicated that this protein adopts a noncanonical epidermal growth factor domain-like structural fold, the first such report for tick salivary proteins. These data support a mechanism by which tick saliva proteins simultaneously inhibit both the host coagulation cascade and the lectin complement pathway.

The structure of a minimum amyloid fibril core formed by necroptosis-mediating RHIM of human RIPK3

Receptor-interacting protein kinases 3 (RIPK3), a central node in necroptosis, polymerizes in response to the upstream signals and then activates its downstream mediator to induce cell death. The active polymeric form of RIPK3 has been indicated as the form of amyloid fibrils assembled via its RIP homotypic interaction motif (RHIM). In this study, we combine cryogenic electron microscopy and solid-state NMR to determine the amyloid fibril structure of RIPK3 RHIM-containing C-terminal domain (CTD). The structure reveals a single protofilament composed of the RHIM domain. RHIM forms three β-strands (referred to as strands 1 through 3) folding into an S shape, a distinct fold from that in complex with RIPK1. The consensus tetrapeptide VQVG of RHIM forms strand 2, which zips up strands 1 and 3 via heterozipper-like interfaces. Notably, the RIPK3-CTD fibril, as a physiological fibril, exhibits distinctive assembly compared with pathological fibrils. It has an exceptionally small fibril core and twists in both handedness with the smallest pitch known so far. These traits may contribute to a favorable spatial arrangement of RIPK3 kinase domain for efficient phosphorylation.

The amyloid structure of mouse RIPK3 (receptor interacting protein kinase 3) in cell necroptosis

RIPK3 amyloid complex plays crucial roles during TNF-induced necroptosis and in response to immune defense in both human and mouse. Here, we have structurally characterized mouse RIPK3 homogeneous self-assembly using solid-state NMR, revealing a well-ordered N-shaped amyloid core structure featured with 3 parallel in-register β-sheets. This structure differs from previously published human RIPK1/RIPK3 hetero-amyloid complex structure, which adopted a serpentine fold. Functional studies indicate both RIPK1-RIPK3 binding and RIPK3 amyloid formation are essential but not sufficient for TNF-induced necroptosis. The structural integrity of RIPK3 fibril with three β-strands is necessary for signaling. Molecular dynamics simulations with a mouse RIPK1/RIPK3 model indicate that the hetero-amyloid is less stable when adopting the RIPK3 fibril conformation, suggesting a structural transformation of RIPK3 from RIPK1-RIPK3 binding to RIPK3 amyloid formation. This structural transformation would provide the missing link connecting RIPK1-RIPK3 binding to RIPK3 homo-oligomer formation in the signal transduction.

Receptor Interacting Protein Kinase 3 (RIPK3) has a key role in TNF-induced necroptosis. Here, the authors combine solid state NMR measurements, MD simulations and cell based assays to characterize mouse RIPK3 and they present the structure of the RIPK3 amyloid core.

Solution structures of the Shewanella woodyi H-NOX protein in the presence and absence of soluble guanylyl cyclase stimulator IWP-051

Heme-nitric oxide/oxygen binding (H-NOX) domains bind gaseous ligands for signal transduction in organisms spanning prokaryotic and eukaryotic kingdoms. In the bioluminescent marine bacterium Shewanella woodyi (Sw), H-NOX proteins regulate quorum sensing and biofilm formation. In higher animals, soluble guanylyl cyclase (sGC) binds nitric oxide with an H-NOX domain to induce cyclase activity and regulate vascular tone, wound healing and memory formation. sGC also binds stimulator compounds targeting cardiovascular disease. The molecular details of stimulator binding to sGC remain obscure but involve a binding pocket near an interface between H-NOX and coiled-coil domains. Here, we report the full NMR structure for CO-ligated Sw H-NOX in the presence and absence of stimulator compound IWP-051, and its backbone dynamics. Nonplanar heme geometry was retained using a semi-empirical quantum potential energy approach. Although IWP-051 binding is weak, a single binding conformation was found at the interface of the two H-NOX subdomains, near but not overlapping with sites identified in sGC. Binding leads to rotation of the subdomains and closure of the binding pocket. Backbone dynamics are similar across both domains except for two helix-connecting loops, which display increased dynamics that are further enhanced by compound binding. Structure-based sequence analyses indicate high sequence diversity in the binding pocket, but the pocket itself appears conserved among H-NOX proteins. The largest dynamical loop lies at the interface between Sw H-NOX and its binding partner as well as in the interface with the coiled coil in sGC, suggesting a critical role for the loop in signal transduction.

Increased usability, algorithmic improvements and incorporation of data mining for structure calculation of proteins with REDCRAFT software package

Background

Traditional approaches to elucidation of protein structures by Nuclear Magnetic Resonance spectroscopy (NMR) rely on distance restraints also known as Nuclear Overhauser effects (NOEs). The use of NOEs as the primary source of structure determination by NMR spectroscopy is time consuming and expensive. Residual Dipolar Couplings (RDCs) have become an alternate approach for structure calculation by NMR spectroscopy. In previous works, the software package REDCRAFT has been presented as a means of harnessing the information containing in RDCs for structure calculation of proteins. However, to meet its full potential, several improvements to REDCRAFT must be made.

Results

In this work, we present improvements to REDCRAFT that include increased usability, better interoperability, and a more robust core algorithm. We have demonstrated the impact of the improved core algorithm in the successful folding of the protein 1A1Z with as high as ±4 Hz of added error. The REDCRAFT computed structure from the highly corrupted data exhibited less than 1.0 Å with respect to the X-ray structure. We have also demonstrated the interoperability of REDCRAFT in a few instances including with PDBMine to reduce the amount of required data in successful folding of proteins to unprecedented levels. Here we have demonstrated the successful folding of the protein 1D3Z (to within 2.4 Å of the X-ray structure) using only N-H RDCs from one alignment medium.

Conclusions

The additional GUI features of REDCRAFT combined with the NEF compliance have significantly increased the flexibility and usability of this software package. The improvements of the core algorithm have substantially improved the robustness of REDCRAFT in utilizing less experimental data both in quality and quantity.

Membrane proteins in magnetically aligned phospholipid polymer discs for solid-state NMR spectroscopy

Well-hydrated phospholipid bilayers provide a near-native environment for membrane proteins. They enable the preparation of chemically-defined samples suitable for NMR and other spectroscopic experiments that reveal the structure, dynamics, and functional interactions of the proteins at atomic resolution. The synthetic polymer styrene maleic acid (SMA) can be used to prepare detergent-free samples that form macrodiscs with diameters greater than 30 nm at room temperature, and spontaneously align in the magnetic field of an NMR spectrometer at temperatures above 35 °C. Here we show that magnetically aligned macrodiscs are particularly well suited for solid-state NMR experiments of membrane proteins because the SMA-lipid assembly both immobilizes the embedded protein and provides uniaxial order for oriented sample (OS) solid-state NMR studies. We show that aligned macrodiscs incorporating four different membrane proteins with a wide range of sizes and topological complexity yield high-resolution OS solid-state NMR spectra. The work is dedicated to Michelle Auger who made key contributions to the field of membrane and membrane protein biophysics.

Defining the mobility range of a hinge-type connection using molecular dynamics and metadynamics

A designed disulfide-rich β-hairpin peptide that dimerizes spontaneously served as a hinge-type connection between proteins. Here, we analyze the range of dynamics of this hinge dimer with the aim of proposing new applications for the DNA-encodable peptide and establishing guidelines for the computational analysis of other disulfide hinges. A recent structural analysis based on nuclear magnetic resonance spectroscopy and ion mobility spectrometry revealed an averaged conformation in the hinge region which motivated us to investigate the dynamic behavior using a combination of molecular dynamics simulation, metadynamics and free energy surface analysis to characterize the conformational space available to the hinge. Principal component analysis uncovered two slow modes of the peptide, namely, the opening and closing motion and twisting of the two β-hairpins assembling the hinge. Applying a collective variable (CV) that mimics the first dominating mode, led to a major expansion of the conformational space. The description of the dynamics could be achieved by analysis of the opening angle and the twisting of the β-hairpins and, thus, offers a methodology that can also be transferred to other derivatives. It has been demonstrated that the hinge peptide’s lowest energy conformation consists of a large opening angle and strong twist but is separated by small energy barriers and can, thus, adopt a closed and untwisted structure. With the aim of proposing further applications for the hinge peptide, we simulated its behavior in the sterically congested environment of a four-helix bundle. Preliminary investigations show that one helix is pushed out and a three-helix bundle forms. The insights gained into the dynamics of the tetra-disulfide peptide and analytical guidelines developed in this study may contribute to the understanding of the structure and function of more complex hinge-type proteins, such as the IgG antibody family.

Comparing the Hinge-Type Mobility of Natural and Designed Intermolecular Bi-disulfide Domains

A pair of intermolecular disulfide bonds connecting two protein domains restricts their relative mobility in a systematic way. The bi-disulfide hinge cannot rotate like a single intermolecular disulfide bond yet is less restrained than three or more intermolecular disulfides which restrict the relative motion to a minimum. The intermediate mobility of bi-disulfide linked domains is characterized by their dominating opening and closing modes comparable to the mechanics of a door hinge on the macroscopic scale. Here we compare the central hinge region of Immunoglobulin G1 (IgG1) which is highly conserved among different species, with a recently designed hinge-type motif CHWECRGCRLVC from our lab, that was successfully used for the dimerization of the IgG1/κ-ab CL4 monocolonal antibody (mab). The minimal length of these synthetic hinges comprises only 12 amino acids, rendering them ideal models for computational studies. Well-tempered metadynamics was performed to adequately describe the available conformational space defined by the different hinges. In spite of the differences in amino acid composition and ring sizes, there are characteristic similarities of designed and natural hinges like the dependent mobility of the individual strands of each hinge domain.

Integrated NMR and cryo-EM atomic-resolution structure determination of a half-megadalton enzyme complex

Atomic-resolution structure determination is crucial for understanding protein function. Cryo-EM and NMR spectroscopy both provide structural information, but currently cryo-EM does not routinely give access to atomic-level structural data, and, generally, NMR structure determination is restricted to small (<30 kDa) proteins. We introduce an integrated structure determination approach that simultaneously uses NMR and EM data to overcome the limits of each of these methods. The approach enables structure determination of the 468 kDa large dodecameric aminopeptidase TET2 to a precision and accuracy below 1 Å by combining secondary-structure information obtained from near-complete magic-angle-spinning NMR assignments of the 39 kDa-large subunits, distance restraints from backbone amides and ILV methyl groups, and a 4.1 Å resolution EM map. The resulting structure exceeds current standards of NMR and EM structure determination in terms of molecular weight and precision. Importantly, the approach is successful even in cases where only medium-resolution cryo-EM data are available.

NMR structure determination is challenging for proteins with a molecular weight above 30 kDa and atomic-resolution structure determination from cryo-EM data is currently not the rule. Here the authors describe an integrated structure determination approach that simultaneously uses NMR and EM data and allows them to determine the structure of the 468 kDa dodecameric aminopeptidase TET2 complex.

Validation of protein backbone structures calculated from NMR angular restraints using Rosetta

Multidimensional solid-state NMR spectra of oriented membrane proteins can be used to infer the backbone torsion angles and hence the overall protein fold by measuring dipolar couplings and chemical shift anisotropies, which depend on the orientation of each peptide plane with respect to the external magnetic field. However, multiple peptide plane orientations can be consistent with a given set of angular restraints. This ambiguity is further exacerbated by experimental uncertainty in obtaining and interpreting such restraints. The previously developed algorithms for structure calculations using angular restraints typically involve a sequential walkthrough along the backbone to find the torsion angles between the consecutive peptide plane orientations that are consistent with the experimental data. This method is sensitive to experimental uncertainty in interpreting the peak positions of as low as ± 10 Hz, often yielding high structural RMSDs for the calculated structures. Here we present a significantly improved version of the algorithm which includes the fitting of several peptide planes at once in order to prevent propagation of error along the backbone. In addition, a protocol has been devised for filtering the structural solutions using Rosetta scoring functions in order to find the structures that both fit the spectrum and satisfy bioinformatics restraints. The robustness of the new algorithm has been tested using synthetic angular restraints generated from the known structures for two proteins: a soluble protein 2gb1 (56 residues), chosen for its diverse secondary structure elements, i.e. an alpha-helix and two beta-sheets, and a membrane protein 4a2n, from which the first two transmembrane helices (having a total of 64 residues) have been used. Extensive simulations have been performed by varying the number of fitted planes, experimental error, and the number of NMR dimensions. It has been found that simultaneously fitting two peptide planes always shifted the distribution of the calculated structures toward lower structural RMSD values as compared to fitting a single torsion-angle pair. For each protein, irrespective of the simulation parameters, Rosetta was able to distinguish the most plausible structures, often having structural RMSDs lower than 2 Å with respect to the original structure. This study establishes a framework for de-novo protein structure prediction using a combination of solid-state NMR angular restraints and bioinformatics.

Xplor-NIH for molecular structure determination from NMR and other data sources

Xplor-NIH is a popular software package for biomolecular structure determination from nuclear magnetic resonance (NMR) and other data sources. Here, some of Xplor-NIH's most useful data-associated energy terms are reviewed, including newer alternative options for using residual dipolar coupling data in structure calculations. Further, we discuss new developments in the implementation of strict symmetry for the calculation of symmetric homo-oligomers, and in the representation of the system as an ensemble of structures to account for motional effects. Finally, the different available force fields are presented, among other Xplor-NIH capabilities.

Structural Insights into the Yersinia pestis Outer Membrane Protein Ail in Lipid Bilayers

Yersinia pestis the causative agent of plague, is highly pathogenic and poses very high risk to public health. The outer membrane protein Ail (Adhesion invasion locus) is one of the most highly expressed proteins on the cell surface of Y. pestis, and a major target for the development of medical countermeasures. Ail is essential for microbial virulence and is critical for promoting the survival of Y. pestis in serum. Structures of Ail have been determined by X-ray diffraction and solution NMR spectroscopy, but the protein's activity is influenced by the detergents in these samples, underscoring the importance of the surrounding environment for structure-activity studies. Here we describe the backbone structure of Ail, determined in lipid bilayer nanodiscs, using solution NMR spectroscopy. We also present solid-state NMR data obtained for Ail in membranes containing lipopolysaccharide (LPS), a major component of the bacterial outer membranes. The protein in lipid bilayers, adopts the same eight-stranded β-barrel fold observed in the crystalline and micellar states. The membrane composition, however, appears to have a marked effect on protein dynamics, with LPS enhancing conformational order and slowing down the ¹⁵N transverse relaxation rate. The results provide information about the way in which an outer membrane protein inserts and functions in the bacterial membrane.

Applications of NMR to membrane proteins

Membrane proteins present a challenge for structural biology. In this article, we review some of the recent developments that advance the application of NMR to membrane proteins, with emphasis on structural studies in detergent-free, lipid bilayer samples that resemble the native environment. NMR spectroscopy is not only ideally suited for structure determination of membrane proteins in hydrated lipid bilayer membranes, but also highly complementary to the other principal techniques based on X-ray and electron diffraction. Recent advances in NMR instrumentation, spectroscopic methods, computational methods, and sample preparations are driving exciting new efforts in membrane protein structural biology.