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) stalled at a DNA signal known as ops. 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 ops 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, most studies employ tertiary biases towards each native state, hampering the analysis of sequence-encoded interactions on fold-switching. 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 also by favoring β-barrel rupture into an intermediate from which fold-switching proceeds.

Download Full-text

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

PLoS Computational Biology ◽

10.1371/journal.pcbi.1008882 ◽

2021 ◽

Vol 17 (9) ◽

pp. e1008882

Author(s):

Pablo Galaz-Davison ◽

Ernesto A. Román ◽

César A. Ramírez-Sarmiento

Keyword(s):

Molecular Dynamics ◽

Virulence Factors ◽

Elongation Factor ◽

Active Role ◽

Structural Rearrangement ◽

Umbrella Sampling ◽

Protein Fold ◽

Terminal Domain ◽

Rnap Binding

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.

Download Full-text

Differential local stability governs the metamorphic fold-switch of bacterial virulence factor RfaH

10.1101/629477 ◽

2019 ◽

Author(s):

P. Galaz-Davison ◽

J.A. Molina ◽

S. Silletti ◽

E.A. Komives ◽

S.H. Knauer ◽

...

Keyword(s):

Molecular Dynamics ◽

Local Stability ◽

Large Scale ◽

Target Genes ◽

Structural Rearrangement ◽

Virulence Determinants ◽

Gram Negative ◽

Rapid Acquisition ◽

Terminal Domain ◽

Polymerase Binding

AbstractA regulatory factor RfaH, present in many Gram-negative bacterial pathogens, is required for transcription and translation of long operons encoding virulence determinants. Escherichia coli RfaH action is controlled by a unique large-scale structural rearrangement triggered by recruitment to transcription elongation complexes through a specific DNA sequence within these operons. Upon recruitment, the C-terminal domain of this two-domain protein refolds from an α-hairpin, which is bound to the RNA polymerase binding site within the N-terminal domain of RfaH, into an unbound β-barrel that interacts with the ribosome to enable translation. Although structures of the autoinhibited (α-hairpin) and active (β-barrel) states and plausible refolding pathways have been reported, how this reversible switch is encoded within RfaH sequence and structure is poorly understood. Here, we combined hydrogen-deuterium exchange measurements by mass spectrometry and nuclear magnetic resonance with molecular dynamics to evaluate the differential local stability between both RfaH folds. Deuteron incorporation reveals that the tip of the C-terminal hairpin (residues 125-145) is stably folded in the autoinhibited state (∼20% deuteron incorporation), while the rest of this domain is highly flexible (>40% deuteron incorporation) and its flexibility only decreases in the β-folded state. Computationally-predicted ΔGs agree with these results by displaying similar anisotropic stability within the tip of the α-hairpin and on neighboring N-terminal domain residues. Remarkably, the β-folded state shows comparable stability to non-metamorphic homologs. Our findings provide information critical for understanding the metamorphic behavior of RfaH and other chameleon proteins, and for devising targeted strategies to combat bacterial diseases.SignificanceInfections caused by Gram-negative bacteria are a worldwide health threat due to rapid acquisition of antibiotic resistance. RfaH, a protein essential for virulence in several Gram-negative pathogens, undergoes a large-scale structural rearrangement in which one RfaH domain completely refolds. Refolding transforms RfaH from an inactive state that restricts RfaH recruitment to a few target genes into an active state that binds to, and couples, transcription and translation machineries to elicit dramatic activation of gene expression. However, the molecular basis of this unique conformational change is poorly understood. Here, we combine molecular dynamics and structural biology to unveil the hotspots that differentially stabilize both states of RfaH. Our findings provide novel insights that will guide design of inhibitors blocking RfaH action.

Download Full-text

Surfactant desorption and scission free energies for cylindrical and spherical micelles from umbrella-sampling molecular dynamics simulations

Journal of Colloid and Interface Science ◽

10.1016/j.jcis.2021.04.138 ◽

2021 ◽

Author(s):

Boyao Wen ◽

Bofeng Bai ◽

Ronald G. Larson

Keyword(s):

Molecular Dynamics ◽

Molecular Dynamics Simulations ◽

Umbrella Sampling ◽

Free Energies ◽

Spherical Micelles ◽

Dynamics Simulations

Download Full-text

Study of Endogen Substrates, Drug Substrates and Inhibitors Binding Conformations on MRP4 and Its Variants by Molecular Docking and Molecular Dynamics

Molecules ◽

10.3390/molecules26041051 ◽

2021 ◽

Vol 26 (4) ◽

pp. 1051

Author(s):

Edgardo Becerra ◽

Giovanny Aguilera-Durán ◽

Laura Berumen ◽

Antonio Romo-Mancillas ◽

Guadalupe García-Alcocer

Keyword(s):

Molecular Dynamics ◽

Molecular Docking ◽

Binding Energy ◽

Binding Site ◽

Binding Sites ◽

Adenosine Monophosphate ◽

Competitive Inhibitor ◽

Umbrella Sampling ◽

Coarse Grained ◽

Nucleotide Polymorphisms

Multidrug resistance protein-4 (MRP4) belongs to the ABC transporter superfamily and promotes the transport of xenobiotics including drugs. A non-synonymous single nucleotide polymorphisms (nsSNPs) in the ABCC4 gene can promote changes in the structure and function of MRP4. In this work, the interaction of certain endogen substrates, drug substrates, and inhibitors with wild type-MRP4 (WT-MRP4) and its variants G187W and Y556C were studied to determine differences in the intermolecular interactions and affinity related to SNPs using protein threading modeling, molecular docking, all-atom, coarse grained, and umbrella sampling molecular dynamics simulations (AA-MDS and CG-MDS, respectively). The results showed that the three MRP4 structures had significantly different conformations at given sites, leading to differences in the docking scores (DS) and binding sites of three different groups of molecules. Folic acid (FA) had the highest variation in DS on G187W concerning WT-MRP4. WT-MRP4, G187W, Y556C, and FA had different conformations through 25 ns AA-MD. Umbrella sampling simulations indicated that the Y556C-FA complex was the most stable one with or without ATP. In Y556C, the cyclic adenosine monophosphate (cAMP) and ceefourin-1 binding sites are located out of the entrance of the inner cavity, which suggests that both cAMP and ceefourin-1 may not be transported. The binding site for cAMP and ceefourin-1 is quite similar and the affinity (binding energy) of ceefourin-1 to WT-MRP4, G187W, and Y556C is greater than the affinity of cAMP, which may suggest that ceefourin-1 works as a competitive inhibitor. In conclusion, the nsSNPs G187W and Y556C lead to changes in protein conformation, which modifies the ligand binding site, DS, and binding energy.

Download Full-text