The N-terminal domain of RfaH plays an active role in protein fold-switching

The bacterial elongation factor RfaH promotes the expression of virulence factors by specifically binding to RNA polymerases (RNAP) paused at a DNA signal. This behavior is unlike that of its paralog NusG, the major representative of the protein family to which RfaH belongs. Both proteins have an N-terminal domain (NTD) bearing an RNAP binding site, yet NusG C-terminal domain (CTD) is folded as a β-barrel while RfaH CTD is forming an α-hairpin blocking such site. Upon recognition of the specific DNA exposed by RNAP, RfaH is activated via interdomain dissociation and complete CTD structural rearrangement into a β-barrel structurally identical to NusG CTD. Although RfaH transformation has been extensively characterized computationally, little attention has been given to the role of the NTD in the fold-switching process, as its structure remains unchanged. Here, we used Associative Water-mediated Structure and Energy Model (AWSEM) molecular dynamics to characterize the transformation of RfaH, spotlighting the sequence-dependent effects of NTD on CTD fold stabilization. Umbrella sampling simulations guided by native contacts recapitulate the thermodynamic equilibrium experimentally observed for RfaH and its isolated CTD. Temperature refolding simulations of full-length RfaH show a high success towards α-folded CTD, whereas the NTD interferes with βCTD folding, becoming trapped in a β-barrel intermediate. Meanwhile, NusG CTD refolding is unaffected by the presence of RfaH NTD, showing that these NTD-CTD interactions are encoded in RfaH sequence. Altogether, these results suggest that the NTD of RfaH favors the α-folded RfaH by specifically orienting the αCTD upon interdomain binding and by favoring β-barrel rupture into an intermediate from which fold-switching proceeds.     

Interdomain contacts control folding of transcription factor RfaH

Escherichia coli RfaH activates gene expression by tethering the elongating RNA polymerase to the ribosome. This bridging action requires a complete refolding of the RfaH C-terminal domain (CTD) from an α-helical hairpin, which binds to the N-terminal domain (NTD) in the free protein, to a β-barrel, which interacts with the ribosomal protein S10 following RfaH recruitment to its target operons. The CTD forms a β-barrel when expressed alone or proteolytically separated from the NTD, indicating that the α-helical state is trapped by the NTD, perhaps co-translationally. Alternatively, the interdomain contacts may be sufficient to drive the formation of the α-helical form. Here, we use functional and NMR analyses to show that the denatured RfaH refolds into the native state and that RfaH in which the order of the domains is reversed is fully functional in vitro and in vivo. Our results indicate that all information necessary to determine its fold is encoded within RfaH itself, whereas accessory factors or sequential folding of NTD and CTD during translation are dispensable. These findings suggest that universally conserved RfaH homologs may change folds to accommodate diverse interaction partners and that context-dependent protein refolding may be widespread in nature.     

High-resolution structure of a modular hyperthermostable endo-β-1,4-mannanase from Thermotoga petrophila: The ancillary immunoglobulin-like module is a thermostabilizing domain

The endo-β-1,4-mannanase from the hyperthermostable bacterium Thermotoga petrophila (TpMan) is an enzyme that catalyzes the hydrolysis of mannan and heteromannan polysaccharides. Of the three domains that comprise TpMan, the N-terminal GH5 catalytic domain and the C-terminal carbohydrate-binding domain are connected through a central ancillary domain of unknown structure and function. In this study, we report the partial crystal structure of the TpMan at 1.45 Å resolution, so far, the first modular hyperthermostable endo-β-1,4-mannanase structure determined. The structure exhibits two domains, a (β/α)₈-barrel GH5 catalytic domain connected via a linker to the central domain with an immunoglobulin-like β-sandwich fold formed of seven β-strands. Functional analysis showed that whereas the immunoglobulin-like domain does not have the carbohydrate-binding function, it stacks on the GH5 catalytic domain acting as a thermostabilizing domain and allowing operation at hyperthermophilic conditions. The carbohydrate-binding domain is absent in the crystal structure most likely due to its high flexibility around the immunoglobulin-like domain which may act also as a pivot. These results represent new structural and functional information useful on biotechnological applications for biofuel and food industries.     

Flow-induced structural transition in the beta-switch region of glycoprotein Ib

The impact of fluid flow on structure and dynamics of biomolecules has recently gained much attention. In this article, we present a molecular-dynamics algorithm that serves to generate stable water flow under constant temperature, for the study of flow-induced protein behavior. Flow simulations were performed on the 16-residue β-switch region of platelet glycoprotein Ibα, for which crystal structures of its N-terminal domain alone and in complex with the A1 domain of von Willebrand factor have been solved. Comparison of the two structures reveals a conformational change in this region, which, upon complex formation, switches from an unstructured loop to a β-hairpin. Interaction between glycoprotein Ibα and von Willebrand factor initiates platelet adhesion to injured vessel walls, and the adhesion is enhanced by blood flow. It has been hypothesized that the loop to β-hairpin transition in glycoprotein Ibα is induced by flow before binding to von Willebrand factor. The simulations revealed clearly a flow-induced loop→β-hairpin transition. The transition is dominated by the entropy of the protein, and is seen to occur in two steps, namely a dihedral rotation step followed by a side-group packing step.     

A multidomain outer membrane protein from Pasteurella multocida: Modelling and simulation studies of PmOmpA

PmOmpA is a two-domain outer membrane protein from Pasteurella multocida. The N-terminal domain of PmOmpA is a homologue of the transmembrane β-barrel domain of OmpA from Escherichia coli, whilst the C-terminal domain of PmOmpA is a homologue of the extra-membrane Neisseria meningitidis RmpM C-terminal domain. This enables a model of a complete two domain PmOmpA to be constructed and its conformational dynamics explored via MD simulations of the protein embedded within two different phospholipid bilayers (DMPC and DMPE). The conformational stability of the transmembrane β-barrel is similar to that of a homology model of OprF from Pseudomonas aeruginosa in bilayer simulations. There is a degree of water penetration into the interior of the β-barrel, suggestive of a possible transmembrane pore. Although the PmOmpA model is stable over 20 ns simulations, retaining its secondary structure and fold integrity throughout, substantial flexibility is observed in a short linker region between the N- and the C-terminal domains. At low ionic strength, the C-terminal domain moves to interact electrostatically with the lipid bilayer headgroups. This study demonstrates that computational approaches may be applied to more complex, multi-domain outer membrane proteins, rather than just to transmembrane β-barrels, opening the possibility of in silico proteomics approaches to such proteins.     

Circular permutants in β-glucosidases (family 3) within a predicted double-domain topology that includes a (β/α)8-barrel

By predicting the general secondary structure for β-glucosidases (family 3), in conjunction with existing knowledge of the circular permutants present in B. fibrisolvens and R. albus, we were able to find the canonical elements of the secondary structure. The way these elements are linked suggests that there is a double-domain topology made up of a (β/α)₈-barrel domain and a “mainly all-β” domain. A number of already known conserved motifs are located within (or near) the C-terminal part of the putative parallel β-strands of the (β/α)₈-barrel, which is consistent with what is known about the location of catalytical sites for enzymes that have this domain topology. Within the circular permutants, two β/α units are located at the N-terminal part of the molecule, whereas the other six β/α units are located at the C-terminal end. In this way, the circular permutants can be seen to have a putative discontinuous double-domain topology. Proteins 31:214–223, 1998. © 1998 Wiley-Liss, Inc.     

Molecular dynamics simulation of mammalian 15S-lipoxygenase with AMBER force field

A molecular dynamics simulation study of mononuclear iron 15S-lipoxygenase (15S-LOX) from rabbit reticulocytes was performed to investigate its structure and dynamics; newly developed AMBER force field parameters were employed for the first coordination sphere of the catalytic iron (II). The results obtained from this study demonstrate that the structural features of the catalytic iron coordination site are in good agreement with available data obtained from experiments. The motional flexibility of the N-terminal β-barrel domain is greater than the C-terminal catalytic domain; flexibility was assessed in terms of B-factors and secondary structure calculations. The significant features obtained for the relative motional flexibility of these two domains of 15S-LOX in solution as well as the isolated C-terminal domain were analyzed in terms of radius of gyration and maximum diameter, which correlated well with the structural flexibility of 15-lipoxygenase-1 in solution as probed by small-angle X-ray scattering. The motional flexibility indicates interdomain motion between the N-terminal β-barrel and the C-terminal catalytic domain; this was further verified by the evaluation of central bending in the solvated LOX molecule, which identified an unstructured stretch of amino acids as the interdomain linker. The average bending angle confirmed significant central bending between these two domains, which was linked to the high degree of motional freedom of the N-terminal β-barrel domain in aqueous solutions. This can be considered to have biological relevance for membrane binding as well as for regulating the catalytic domain.     

Crystal structure of Arabidopsis cyclophilin38 reveals a previously uncharacterized immunophilin fold and a possible autoinhibitory mechanism

Cyclophilin38 (CYP38) is one of the highly divergent cyclophilins from Arabidopsis thaliana. Here, we report the crystal structure of the At-CYP38 protein (residues 83 to 437 of 437 amino acids) at 2.39-Å resolution. The structure reveals two distinct domains: an N-terminal helical bundle and a C-terminal cyclophilin β-barrel, connected by an acidic loop. Two N-terminal β-strands become part of the C-terminal cyclophilin β-barrel, thereby making a previously undiscovered domain organization. This study shows that CYP38 does not possess peptidyl-prolyl cis/trans isomerase activity and identifies a possible interaction of CYP38 with the E-loop of chlorophyll protein47 (CP47), a component of photosystem II. The interaction of CYP38 with the E-loop of CP47 is mediated through its cyclophilin domain. The N-terminal helical domain is closely packed together with the putative C-terminal cyclophilin domain and establishes a strong intramolecular interaction, thereby preventing the access of the cyclophilin domain to other proteins. This was further verified by protein–protein interaction assays using the yeast two-hybrid system. Furthermore, the non-Leucine zipper N-terminal helical bundle contains several new elements for protein–protein interaction that may be of functional significance. Together, this study provides the structure of a plant cyclophilin and explains a possible mechanism for autoinhibition of its function through an intramolecular interaction.     

Structural characterization of a rhamnose-binding glycoprotein (lectin) from Spanish mackerel (Scomberomorous niphonius) eggs

A rhamnose-binding glycoprotein (lectin), named SML, was isolated from the eggs of Spanish mackerel (Scomberomorous niphonius) by affinity and ion-exchange chromatographies. SML was composed of a non-covalently linked homodimer. The SML subunit was composed of 201 amino acid residues with two tandemly repeated domains, and contained 8 half-Cys residues in each domain, which is highly homologous to the N-terminal lectin domain of calcium-independent α-latrotoxin receptor in mammalian brains. Each domain has the same disulfide bonding pattern; Cys10–Cys40, Cys20–Cys99, Cys54–Cys86 and Cys67–Cys73 were located in the N-terminal domain, and Cys108–Cys138, Cys117–Cys195, Cys152–Cys182 and Cys163–Cys169 were in the C-terminal domain. SML was N-glycosylated at Asn168 in the C-terminal domain. The structure of the sugar chain was determined to be NeuAc-Galβ1-4GlcNAcβ1-2Manα1-6-(NeuAc-Galβ1-4GlcNAcβ1-2Manα1-3)Manβ1-4GlcNAcβ1-4GlcNAc-Asn.     

The C-terminal domain of Escherichia coli dihydrodipicolinate synthase (DHDPS) is essential for maintenance of quaternary structure and efficient catalysis

Dihydrodipicolinate synthase (DHDPS) catalyses the first committed step in the biosynthesis of (S)-lysine, an essential constituent of bacterial cell walls. Escherichia coli DHDPS is homotetrameric, and each monomer contains an N-terminal (β/α)₈-barrel, responsible for catalysis and regulation, and three C-terminal α-helices, the function of which is unknown. This study investigated the C-terminal domain of E. coli DHDPS by characterising a C-terminal truncated DHDPS (DHDPS-H225∗). DHDPS-H225∗ was unable to complement an (S)-lysine auxotroph, and showed significantly reduced solubility, stability, and maximum catalytic activity (k_cat = 1.20 ± 0.01 s⁻¹), which was only 1.6% of wild type E. coli DHDPS (DHDPS-WT). The affinity of DHDPS-H225∗ for substrates and the feedback inhibitor, (S)-lysine, remained comparable to DHDPS-WT. These changes were accompanied by disruption in the quaternary structure, which has previously been shown to be essential for efficient catalysis in this enzyme.     

Striped murrel S1 family serine protease: immune characterization,antibacterial property and enzyme activities

We report a molecular characterization of S1 family serine protease (SP-1) from snakehead murrel (or called striped murrel) Channa striatus (Cs). CsSP-1 polypeptide contained a catalytic core domain (otherwise known as serine protease trypsin domain) between H20 and I237 along with a catalytic triad at H⁶¹, D¹⁰⁴ and S¹⁹⁷. Phylogenetic analysis confirmed that CsSP-1 belongs to serine protease S1 family. The tertiary structure showed that CsSP-1 contains 14 β-sheets as 2 separate β-barrels (the first β-barrel consists of 8 β-sheets in the N-terminal region and the second β-barrel consists of 6 β-sheets in the C-terminal region) and 3 α-helical regions. Significantly (P < 0.05) the highest CsSP-1 mRNA expression was observed in intestine, liver and kidney, moderate expression was seen in spleen, head kidney, skin and blood, and the lowest one in brain, gill, muscle and heart. Further, the expression was induced in intestine with fungus Aphanomyces invadans and bacteria Aeromonas hydrophila. The recombinant CsSP-1 protein showed antibacterial activity against both gram-negative and gram-positive bacteria. The optimum CsSP-1 enzyme activity against the substrate casein was determined at 8 mM casein concentration. Moreover, the activity was highly influenced by 5 mM phenyl-methylsulfonyl fluoride followed by ethylenediaminetetraacetic acid, 4-(2-aminoethyl)benzenesulfonylfluoride hydrochloride and calpain inhibitor I. The CsSP-1 enzyme exhibited the highest activity at pH 7.5 and temperature 35°C. The overall results showed the potential involvement of CsSP-1 in the immune system of murrels. However, further research is necessary to study the mechanism of implicit trypsin association in the defence process.     

Crystal structure of BamD: an essential component of the β-Barrel assembly machinery of gram-negative bacteria

Folding and insertion of integral β-barrel proteins in the outer membrane (OM) is an essential process for Gram-negative bacteria that requires the β-barrel assembly machinery (BAM). Efficient OM protein (OMP) folding and insertion appears to require a consensus C-terminal signal in OMPs characterized by terminal F or W residues. The BAM complex is embedded in the OM and, in Escherichia coli, consists of the β-barrel BamA and four lipoproteins BamBCDE. BamA and BamD are broadly distributed across all species of Gram-negative bacteria, whereas the other components are present in only a subset of species. BamA and BamD are also essential for viability, suggesting that these two proteins constitute the functional core of the bacterial BAM complex. Here, we present the crystal structure of BamD from the thermophilic bacteria Rhodothermus marinus refined to 2.15 Å resolution. The protein contains five tetratricopeptide repeats (TPRs) organized into two offset tandems, each capped by a terminal helix. The N-terminal domain contains three TPRs and displays remarkable structural similarity with proteins that recognize targeting signals in extended conformations. The C-terminal domain harbors the remaining two TPRs and previously described mutations that impair binding to other BAM components map to this domain. Therefore, the structure suggests a model where the C-terminal domain provides a scaffold for interaction with BAM components, while the N-terminal domain participates in interaction with the substrates, either recognizing the C-terminal consensus sequence or binding unfolded OMP intermediates.     

Effects of Topology and Sequence in Protein Folding Linked via Conformational Fluctuations

Experiments have compared the folding of proteins with different amino acid sequences but the same basic structure, or fold. Results indicate that folding is robust to sequence variations for proteins with some nonlocal folds, such as all-β, whereas the folding of more local, all-α proteins typically exhibits a stronger sequence dependence. Here, we use a coarse-grained model to systematically study how variations in sequence perturb the folding energy landscapes of three model sequences with 3α, 4β + α, and β-barrel folds, respectively. These three proteins exhibit folding features in line with experiments, including expected rank order in the cooperativity of the folding transition and stability-dependent shifts in the location of the free-energy barrier to folding. Using a generalized-ensemble simulation approach, we determine the thermodynamics of around 2000 sequence variants representing all possible hydrophobic or polar single- and double-point mutations. From an analysis of the subset of stability-neutral mutations, we find that folding is perturbed in a topology-dependent manner, with the β-barrel protein being the most robust. Our analysis shows, in particular, that the magnitude of mutational perturbations of the transition state is controlled in part by the size or “width” of the underlying conformational ensemble. This result suggests that the mutational robustness of the folding of the β-barrel protein is underpinned by its conformationally restricted transition state ensemble, revealing a link between sequence and topological effects in protein folding.     

Experimental Evidence for the Existence of a Stable Half-Barrel Subdomain in the (β/α)8-Barrel Fold

The (β/α)₈-barrel is one of the most common folds functioning as enzymes. The emergence of two (β/α)₈-barrel enzymes involved in histidine biosynthesis, each of which has a twofold symmetric structure, has been proposed to be a consequence of tandem duplication and fusion of a (β/α)₄-half-barrel. However, little evidence has been found for the existence of an ancestral half-barrel in the evolution of other (β/α)₈-barrel proteins. In order to detect remnants of an ancestral half-barrel in the (β/α)₈-barrel structure of Escherichia coli N-(5′-phosphoribosyl)anthranilate isomerase, we engineered three potential half-barrel units, (β/α)_1-4, (β/α)_3-6, and (β/α)_5-8. Among these three arrangements, only (β/α)_3-6 is stable; it exists in equilibrium between monomeric and dimeric forms. Thus, the central segment of N-(5′-phosphoribosyl)anthranilate isomerase from E. coli can serve as a half-barrel precursor. A tandem duplication of (β/α)_3-6 yielded predominantly monomeric structures that were quite stable. This result exemplified that the structural characteristics of noncovalently assembled half-barrels could be improved by duplication and fusion. Moreover, our results may provide information regarding the local structural units that encompass interactions important for the early folding events of this ubiquitous protein conformation.     

Alpha-Synuclein Lipid-Dependent Membrane Binding and Translocation through the α-Hemolysin Channel

Gauging the interactions of a natively unfolded Parkinson disease-related protein, alpha-synuclein (α-syn) with membranes and its pathways between and within cells is important for understanding its pathogenesis. Here, to address these questions, we use a robust β-barrel channel, α-hemolysin, reconstituted into planar lipid bilayers. Transient, ∼95% blockage of the channel current by α-syn was observed when 1), α-syn was added from the membrane side where the shorter (stem) part of the channel is exposed; and 2), the applied potential was lower on the side of α-syn addition. While the on-rate of α-syn binding to the channel strongly increased with the applied field, the off-rate displayed a turnover behavior. Statistical analysis suggests that at voltages >50 mV, a significant fraction of the α-syn molecules bound to the channel undergoes subsequent translocation. The observed on-rate varied by >100 times depending on the bilayer lipid composition. Removal of the last 25 amino acids from the highly negatively charged C-terminal of α-syn resulted in a significant decrease in the binding rates. Taken together, these results demonstrate that β-barrel channels may serve as sensitive probes of α-syn interactions with membranes as well as model systems for studies of channel-assisted protein transport.     

Structure of the BamC two-domain protein obtained by Rosetta with a limited NMR data set

The CS-RDC-NOE Rosetta program was used to generate the solution structure of a 27-kDa fragment of the Escherichia coli BamC protein from a limited set of NMR data. The BamC protein is a component of the essential five-protein β-barrel assembly machine in E. coli. The first 100 residues in BamC were disordered in solution. The Rosetta calculations showed that BamC_101-344 forms two well-defined domains connected by an ∼ 18-residue linker, where the relative orientation of the domains was not defined. Both domains adopt a helix-grip fold previously observed in the Bet v 1 superfamily. ¹⁵N relaxation data indicated a high degree of conformational flexibility for the linker connecting the N-terminal domain and the C-terminal domain in BamC. The results here show that CS-RDC-NOE Rosetta is robust and has a high tolerance for misassigned nuclear Overhauser effect restraints, greatly simplifying NMR structure determinations.     

Crystal structure of α-galactosidase from Lactobacillus acidophilus NCFM: insight into tetramer formation and substrate binding

Lactobacillus acidophilus NCFM is a probiotic bacterium known for its beneficial effects on human health. The importance of α-galactosidases (α-Gals) for growth of probiotic organisms on oligosaccharides of the raffinose family present in many foods is increasingly recognized. Here, the crystal structure of α-Gal from L. acidophilus NCFM (LaMel36A) of glycoside hydrolase (GH) family 36 (GH36) is determined by single-wavelength anomalous dispersion. In addition, a 1.58-Å-resolution crystallographic complex with α-d-galactose at substrate binding subsite − 1 was determined. LaMel36A has a large N-terminal twisted β-sandwich domain, connected by a long α-helix to the catalytic (β/α)₈-barrel domain, and a C-terminal β-sheet domain. Four identical monomers form a tightly packed tetramer where three monomers contribute to the structural integrity of the active site in each monomer. Structural comparison of LaMel36A with the monomeric Thermotoga maritima α-Gal (TmGal36A) reveals that O2 of α-d-galactose in LaMel36A interacts with a backbone nitrogen in a glycine-rich loop of the catalytic domain, whereas the corresponding atom in TmGal36A is from a tryptophan side chain belonging to the N-terminal domain. Thus, two distinctly different structural motifs participate in substrate recognition. The tetrameric LaMel36A furthermore has a much deeper active site than the monomeric TmGal36A, which possibly modulates substrate specificity. Sequence analysis of GH36, inspired by the observed structural differences, results in four distinct subgroups having clearly different active-site sequence motifs. This novel subdivision incorporates functional and architectural features and may aid further biochemical and structural analyses within GH36.