Combined NMR and molecular dynamics conformational filter identifies unambiguously dynamic ensembles of Dengue protease NS2B/NS3pro.

The dengue protease NS2B/NS3pro has been reported to adopt either an 'open' or a 'closed' conformation. We have developed a conformational filter that combines NMR with MD simulations to identify conformational ensembles that dominate in solution. Experimental values derived from relaxation parameters for the backbone and methyl side chains were compared with the corresponding back-calculated relaxation parameters of different conformational ensembles obtained from free MD simulations. Our results demonstrate a high prevalence for the 'closed' conformational ensemble while the 'open' conformation is absent, indicating that the latter conformation is most probably due to crystal contacts. Conversely, conformational ensembles in which the positioning of the co-factor NS2B results in a 'partially' open conformation, previously described in both MD simulations and X-ray studies, were identified by our conformational filter. Altogether, we believe that our approach allows for unambiguous identification of true conformational ensembles, an essential step for reliable drug discovery.

Resolution enhancement of NMR by decoupling with the low-rank Hankel model.

Nuclear magnetic resonance (NMR) spectroscopy has become a formidable tool for biochemistry and medicine. Although J-coupling carries essential structural information it may also limit the spectral resolution. Homonuclear decoupling remains a challenging problem. In this work, we introduce a new approach that uses a specific coupling value as prior knowledge, and the Hankel property of the exponential NMR signal to achieve broadband heteronuclear decoupling using the low-rank method. Our results on synthetic and realistic HMQC spectra demonstrate that the proposed method not only effectively enhances resolution by decoupling, but also maintains sensitivity and suppresses spectral artefacts. The approach can be combined with non-uniform sampling, which means that the resolution can be further improved without any extra acquisition time.

NMR structure of emfourin, a novel protein metalloprotease inhibitor: Insights into the mechanism of action.

Emfourin (M4in) is a protein metalloprotease inhibitor recently discovered in the bacterium Serratia proteamaculans and the prototype of a new family of protein protease inhibitors with an unknown mechanism of action. Protealysin-like proteases (PLPs) of the thermolysin family are natural targets of emfourin-like inhibitors widespread in bacteria and known in archaea. The available data indicate the involvement of PLPs in interbacterial interaction as well as bacterial interaction with other organisms and likely in pathogenesis. Arguably, emfourin-like inhibitors participate in the regulation of bacterial pathogenesis by controlling PLP activity. Here, we determined the 3D structure of M4in using solution NMR spectroscopy. The obtained structure demonstrated no significant similarity to known protein structures. This structure was used to model the M4in-enzyme complex and the complex model was verified by small-angle X-ray scattering. Based on the model analysis, we propose a molecular mechanism for the inhibitor, which was confirmed by site-directed mutagenesis. We show that two spatially close flexible loop regions are critical for the inhibitor-protease interaction. One region includes aspartic acid forming a coordination bond with catalytic Zn2+ of the enzyme and the second region carries hydrophobic amino acids interacting with protease substrate binding sites. Such an active site structure corresponds to the noncanonical inhibition mechanism. This is the first demonstration of such a mechanism for protein inhibitors of thermolysin family metalloproteases, which puts forward M4in as a new basis for the development of antibacterial agents relying on selective inhibition of prominent factors of bacterial pathogenesis belonging to this family.

NMR spectrum reconstruction as a pattern recognition problem.

A new deep neural network based on the WaveNet architecture (WNN) is presented, which is designed to grasp specific patterns in the NMR spectra. When trained at a fixed non-uniform sampling (NUS) schedule, the WNN benefits from pattern recognition of the corresponding point spread function (PSF) pattern produced by each spectral peak resulting in the highest quality and robust reconstruction of the NUS spectra as demonstrated in simulations and exemplified in this work on 2D 1H-15N correlation spectra of three representative globular proteins with different sizes: Ubiquitin (8.6 kDa), Azurin (14 kDa), and Malt1 (44 kDa). The pattern recognition by WNN is also demonstrated for successful virtual homo-decoupling in a 2D methyl 1H-13C - HMQC spectrum of MALT1. We demonstrate using WNN that prior knowledge about the NUS schedule, which so far was not been fully exploited, can be used for designing new powerful NMR processing techniques that surpass the existing algorithmic methods.

Assignment of IVL-Methyl side chain of the ligand-free monomeric human MALT1 paracaspase-IgL3 domain in solution.

Mucosa-associated lymphoid tissue protein 1 (MALT1) plays a key role in adaptive immune responses by modulating specific intracellular signalling pathways that control the development and proliferation of both T and B cells. Dysfunction of these pathways is coupled to the progress of highly aggressive lymphoma as well as to potential development of an array of different immune disorders. In contrast to other signalling mediators, MALT1 is not only activated through the formation of the CBM complex together with the proteins CARMA1 and Bcl10, but also by acting as a protease that cleaves multiple substrates to promote lymphocyte proliferation and survival via the NF-κB signalling pathway. Herein, we present the partial 1H, 13C Ile/Val/Leu-Methyl resonance assignment of the monomeric apo form of the paracaspase-IgL3 domain of human MALT1. Our results provide a solid ground for future elucidation of both the three-dimensional structure and the dynamics of MALT1, key for adequate development of inhibitors, and a thorough molecular understanding of its function(s).

1H, 13C and 15N resonance assignment of backbone and IVL-methyl side chain of the S135A mutant NS3pro/NS2B protein of Dengue II virus reveals unique secondary structure features in solution.

The serotype II Dengue (DENV 2) virus is the most prevalent of all four known serotypes. Herein, we present nearly complete 1H, 15N, and 13C backbone and 1H, 13C isoleucine, valine, and leucine methyl resonance assignment of the apo S135A catalytically inactive variant of the DENV 2 protease enzyme folded as a tandem formed between the serine protease domain NS3pro and the cofactor NS2B, as well as the secondary structure prediction of this complex based on the assigned chemical shifts using the TALOS-N software. Our results provide a solid ground for future elucidation of the structure and dynamic of the apo NS3pro/NS2B complex, key for adequate development of inhibitors, and a thorough molecular understanding of their function(s).

NMR assignments and secondary structure distribution of emfourin, a novel proteinaceous protease inhibitor.

Emfourin (M4in) from Serratia proteamaculans is a new proteinaceous inhibitor of protealysin-like proteases (PLPs), a subgroup of the well-known and widely represented metallopeptidase M4 family. Although the biological role of PLPs is debatable, data published indicate their involvement in pathogenesis, including bacterial invasion into eukaryotic cells, suppression of immune defense of some animals, and destruction of plant cell walls. Gene colocalization into a bicistronic operon observed for some PLPs and their inhibitors (as in the case of M4in) implies a mutually consistent functioning of both entities. The originality of the amino acid sequence of M4in suggests it belongs to a previously unknown protein family and this encourages structural studies. In this work, we report a near-complete assignment of 1H, 13C, and 15N resonances of recombinant M4in and its structural-dynamic properties derived from the chemical shifts. According the NMR data analysis, the M4in molecule comprises 3-5 helical elements and 4-6 ß-strands, at least two of which are apparently antiparallel, ascribing this obviously globular protein to the α + ß structural class. Besides, two disordered regions also exist in the central loops between the regular secondary structural elements. The obtained data provide the basis for determining the high-resolution structure as well as functioning mechanism of M4in that can be used for development of new antibacterial therapeutic strategies.

Unambiguous Tracking of Protein Phosphorylation by Fast High-Resolution FOSY NMR*.

Dysregulation of post-translational modifications (PTMs) like phosphorylation is often involved in disease. NMR may elucidate exact loci and time courses of PTMs at atomic resolution and near-physiological conditions but requires signal assignment to individual atoms. Conventional NMR methods for this base on tedious global signal assignment that may often fail, as for large intrinsically disordered proteins (IDPs). We present a sensitive, robust alternative to rapidly obtain only the local assignment near affected signals, based on FOcused SpectroscopY (FOSY) experiments using selective polarisation transfer (SPT). We prove its efficiency by identifying two phosphorylation sites of glycogen synthase kinase 3 beta (GSK3ß) in human Tau40, an IDP of 441 residues, where the extreme spectral dispersion in FOSY revealed unprimed phosphorylation also of Ser409. FOSY may broadly benefit NMR studies of PTMs and other hotspots in IDPs, including sites involved in molecular interactions.

Backbone and side-chain chemical shift assignments for the ribosome-inactivating protein trichobakin (TBK).

Trichobakin (TBK) is a type-I ribosome-inactivating protein (RIP-I), acting as an extremely potent inhibitor of protein synthesis in the cell-free translation system of rabbit reticulocyte lysate (IC50: 3.5 pM). In this respect, TBK surpasses the well-studied highly homologous RIP-I trichosanthin (IC50: 20-27 pM), therefore creation of recombinant toxins based on it is of great interest. TBK needs to penetrate into cytosol through the cell membrane and specifically bind to α-sarcin/ricin loop of 28S ribosome RNA to perform the function of specific RNA depurination. At the moment, there is no detailed structural-dynamic information in solution about diverse states RIP-I can adopt at different stages on the way to protein synthesis inhibition. In this work, we report a near-complete assignment of 1H, 13C, and 15N TBK (27.3 kDa) resonances and analysis of the secondary structure based on the experimental chemical shifts data. This work will serve as a basis for further investigations of the structure, dynamics and interactions of the TBK with its molecular partners using NMR techniques.

Accurate measurement of dipole/dipole transverse cross-correlated relaxation [Formula: see text] in methylenes and primary amines of uniformly [Formula: see text]-labeled proteins.

Side chains possess a broader conformational space (compared to the backbone) and are directly affected by intra- and intermolecular interactions, hence their dynamics and the corresponding NMR relaxation data are more sensitive and informative. Nevertheless, transverse relaxation in [Formula: see text] ([Formula: see text] or [Formula: see text]) spin systems is predominantly non-measurable in uniformly [Formula: see text]-labeled proteins due to cross-correlation effects. In the present publication, we propose a number of pulse sequences for accurate and precise measurement of the dipole-dipole transverse cross-correlated relaxation rate [Formula: see text], which, similarly to [Formula: see text] measurements, provides information about the amplitudes of intramolecular dynamics. The suggested approach has allowed us to circumvent a number of obstacles that were limiting earlier applications of [Formula: see text]: (1) impossibility of transmission of the central component of the triplet of [Formula: see text] group to [Formula: see text]-acquisition via INEPT has been solved by transmission of the averaged signal of "inner" and "outer" components of the triplet; (2) direct recording of the entire triplets resulting in substantial overlap of side chain signals has been replaced by recording of individual singlets with the use of [Formula: see text]-modulated approach and constant-time evolution; (3) low sensitivity has been enhanced via proton acquisition which required special attention to a zero-quantum coherence evolution. The proposed method expands the set of "dynamics sensors" covering protein side chains and substantially improves the quality and the level of detail of experimental data describing dynamic processes in proteins and protein complexes.

Structural basis of the signal transduction via transmembrane domain of the human growth hormone receptor.

BACKGROUND: Prior studies of the human growth hormone receptor (GHR) revealed a distinct role of spatial rearrangements of its dimeric transmembrane domain in signal transduction across membrane. Detailed structural information obtained in the present study allowed elucidating the bases of such rearrangement and provided novel insights into receptor functioning. METHODS: We investigated the dimerization of recombinant TMD fragment GHR254-294 by means of high-resolution NMR in DPC micelles and molecular dynamics in explicit POPC membrane. RESULTS: We resolved two distinct dimeric structures of GHR TMD coexisting in membrane-mimicking micellar environment and providing left- and right-handed helix-helix association via different dimerization motifs. Based on the available mutagenesis data, the conformations correspond to the dormant and active receptor states and are distinguished by cis-trans isomerization of Phe-Pro266 bond in the transmembrane helix entry. Molecular dynamic relaxations of the structures in lipid bilayer revealed the role of the proline residue in functionally significant rearrangements of the adjacent juxtamembrane region supporting alternation between protein-protein and protein-lipid interactions of this region that can be triggered by ligand binding. Also, the importance of juxtamembrane SS bonding for signal persistency, and somewhat unusual aspects of transmembrane region interaction with water molecules were demonstrated. CONCLUSIONS: Two alternative dimeric structures of GHR TMD attributed to dormant and active receptor states interchange via allosteric rearrangements of transmembrane helices and extracellular juxtamembrane regions that support coordination between protein-protein and protein-lipid interactions. GENERAL SIGNIFICANCE: This study provides a holistic vision of GHR signal transduction across the membrane emphasizing the role of protein-lipid interactions.

Alternative packing of EGFR transmembrane domain suggests that protein-lipid interactions underlie signal conduction across membrane.

The human epidermal growth factor receptor (EGFR) of HER/ErbB receptor tyrosine kinase family mediates a broad spectrum of cellular responses transducing biochemical signals via lateral dimerization in plasma membrane, while inactive receptors can exist in both monomeric and dimeric forms. Recently, the dimeric conformation of the helical single-span transmembrane domains of HER/ErbB employing the relatively polar N-terminal motifs in a fashion permitting proper kinase activation was experimentally determined. Here we describe the EGFR transmembrane domain dimerization via an alternative weakly polar C-terminal motif A(661)xxxG(665) presumably corresponding to the inactive receptor state. During association, the EGFR transmembrane helices undergo a structural adjustment with adaptation of inter-molecular polar and hydrophobic interactions depending upon the surrounding membrane properties that directly affect the transmembrane helix packing. This might imply that signal transduction through membrane and allosteric regulation are inclusively mediated by coupled protein-protein and protein-lipid interactions, elucidating paradoxically loose linkage between ligand binding and kinase activation.

NMR-based approach to measure the free energy of transmembrane helix-helix interactions.

Knowledge of the energetic parameters of transmembrane helix-helix interactions is necessary for the establishment of a structure-energy relationship for α-helical membrane domains. A number of techniques have been developed to measure the free energies of dimerization and oligomerization of transmembrane α-helices, and all of these have their advantages and drawbacks. In this study we propose a methodology to determine the magnitudes of the free energy of interactions between transmembrane helices in detergent micelles. The suggested approach employs solution nuclear magnetic resonance (NMR) spectroscopy to determine the population of the oligomeric states of the transmembrane domains and introduces a new formalism to describe the oligomerization equilibrium, which is based on the assumption that both the dimerization of the transmembrane domains and the dissociation of the dimer can occur only upon the collision of detergent micelles. The technique has three major advantages compared with other existing approaches: it may be used to analyze both weak and relatively strong dimerization/oligomerization processes, it works well for the analysis of complex equilibria, e.g. when monomer, dimer and high-order oligomer populations are simultaneously present in the solution, and it can simultaneously yield both structural and energetic characteristics of the helix-helix interaction under study. The proposed methodology was applied to investigate the oligomerization process of transmembrane domains of fibroblast growth factor receptor 3 (FGFR3) and vascular endothelium growth factor receptor 2 (VEGFR2), and allowed the measurement of the free energy of dimerization of both of these objects. In addition the proposed method was able to describe the multi-state oligomerization process of the VEGFR2 transmembrane domain.

Structure of FGFR3 transmembrane domain dimer: implications for signaling and human pathologies.

Fibroblast growth factor receptor 3 (FGFR3) transduces biochemical signals via lateral dimerization in the plasma membrane, and plays an important role in human development and disease. Eight different pathogenic mutations, implicated in cancers and growth disorders, have been identified in the FGFR3 transmembrane segment. Here, we describe the dimerization of the FGFR3 transmembrane domain in membrane-mimicking DPC/SDS (9/1) micelles. In the solved NMR structure, the two transmembrane helices pack into a symmetric left-handed dimer, with intermolecular stacking interactions occurring in the dimer central region. Some pathogenic mutations fall within the helix-helix interface, whereas others are located within a putative alternative interface. This implies that although the observed dimer structure is important for FGFR3 signaling, the mechanism of FGFR3-mediated transduction across the membrane is complex. We propose an FGFR3 signaling mechanism that is based on the solved structure, available structures of isolated soluble FGFR domains, and published biochemical and biophysical data.

Specific membrane binding of neurotoxin II can facilitate its delivery to acetylcholine receptor.

The action of three-finger snake alpha-neurotoxins at their targets, nicotinic acetylcholine receptors (nAChR), is widely studied because of its biological and pharmacological relevance. Most such studies deal only with ligands and receptor models; however, for many ligand/receptor systems the membrane environment may affect ligand binding. In this work we focused on binding of short-chain alpha-neurotoxin II (NTII) from Naja oxiana to the native-like lipid bilayer, and the possible role played by the membrane in delivering the toxin to nAChR. Experimental (NMR and mutagenesis) and molecular modeling (molecular-dynamics simulation) studies revealed a specific interaction of the toxin molecule with the phosphatidylserine headgroup of lipids, resulting in the proper topology of NTII on lipid bilayers favoring the attack of nAChR. Analysis of short-chain alpha-neurotoxins showed that most of them possess a high positive charge and sequence homology in the lipid-binding motif of NTII, implying that interaction with the membrane surrounding nAChR may be common for the toxin family.

Modeling of 31P-NMR spectra of magnetically oriented phospholipid liposomes: A new analytical solution.

31P-NMR spectroscopy is widely used for studies of phospholipid liposomes, a commonly used model of a biological membrane. For the correct analysis of 31P-NMR spectra of the liposomes it is necessary to take into account that they are deformed by the magnetic field of the spectrometer. The liposomes become ellipsoidal and this affects the lineshape of the spectrum. In the present communication we suggest a new analytical formula for modeling of 31P-NMR spectra of the prolate phospholipid liposomes. The formula assumes a Lorentzian broadening function and exactly ellipsoidal shape of the liposomes. Based on the formula a program called P-FIT is designed for the practical analysis of the experimental multicomponent spectra of the prolate liposomes. The versatility of the program developed in a Mathematica environment is demonstrated by simulations of a number of 31P-NMR spectra with different complexity.

Interaction of three-finger toxins with phospholipid membranes: comparison of S- and P-type cytotoxins.

The CTs (cytotoxins) I and II are positively charged three-finger folded proteins from venom of Naja oxiana (the Central Asian cobra). They belong to S- and P-type respectively based on Ser-28 and Pro-30 residues within a putative phospholipid bilayer binding site. Previously, we investigated the interaction of CTII with multilamellar liposomes of dipalmitoylphosphatidylglycerol by wide-line (31)P-NMR spectroscopy. To compare interactions of these proteins with phospholipids, we investigated the interaction of CTI with the multilamellar liposomes of dipalmitoylphosphatidylglycerol analogously. The effect of CTI on the chemical shielding anisotropy and deformation of the liposomes in the magnetic field was determined at different temperatures and lipid/protein ratios. It was found that both the proteins do not affect lipid organization in the gel state. In the liquid crystalline state of the bilayer they disturb lipid packing. To get insight into the interactions of the toxins with membranes, Monte Carlo simulations of CTI and CTII in the presence of the bilayer membrane were performed. It was found that both the toxins penetrate into the bilayer with the tips of all the three loops. However, the free-energy gain on membrane insertion of CTI is smaller (by approximately 7 kcal/mol; 1 kcal identical with 4.184 kJ) when compared with CTII, because of the lower hydrophobicity of the membrane-binding site of CTI. These results clearly demonstrate that the P-type cytotoxins interact with membranes stronger than those of the S-type, although the mode of the membrane insertion is similar for both the types.

Interaction of the P-type cardiotoxin with phospholipid membranes.

The cardiotoxin (cytotoxin II, or CTII) isolated from cobra snake (Naja oxiana) venom is a 60-residue basic membrane-active protein featuring three-finger beta sheet fold. To assess possible modes of CTII/membrane interaction 31P- and 1H-NMR spectroscopy was used to study binding of the toxin and its effect onto multilamellar vesicles (MLV) composed of either zwitterionic or anionic phospholipid, dipalmitoylglycerophosphocholine (Pam2Gro-PCho) or dipalmitoylglycerophosphoglycerol (Pam2Gro-PGro), respectively. The analysis of 1H-NMR linewidths of the toxin and 31P-NMR spectral lineshapes of the phospholipid as a function of temperature, lipid-to-protein ratios, and pH values showed that at least three distinct modes of CTII interaction with membranes exist: (a) nonpenetrating mode; in the gel state of the negatively charged MLV the toxin is bound to the surface electrostatically; the binding to Pam2Gro-PCho membranes was not observed; (b) penetrating mode; hydrophobic interactions develop due to penetration of the toxin into Pam2Gro-PGro membranes in the liquid-crystalline state; it is presumed that in this mode CTII is located at the membrane/water interface deepening the side-chains of hydrophobic residues at the tips of the loops 1-3 down to the boundary between the glycerol and acyl regions of the bilayer; (c) the penetrating mode gives way to isotropic phase, stoichiometrically well-defined CTII/phospholipid complexes at CTII/lipid ratio exceeding a threshold value which was found to depend at physiological pH values upon ionization of the imidazole ring of His31. Biological implications of the observed modes of the toxin-membrane interactions are discussed.