Absence of the MSMEG-3762/63 pump enhances membrane potential
To evaluate the functional role of MSMEG-3762/63 in membrane potential (Δψ), cytofluorimetric analysis was performed on M. smegmatis (wt), M. smegmatis (ΔMSMEG_3763) and M. smegmatis (ΔMSMEG_3763 pBD04) complemented strain, using the DiOC2 fluorescent carbocyanine dye. DiOC2 exhibits green fluorescence in the monomeric state but shifts toward red emission as the dye molecules form aggregates at high cytosolic concentrations caused by higher membrane potentials; the green-red shift disappears in the presence of the depolarizing compound CCCP. Comparative analysis of the three strains showed that ΔMSMEG_3763 had a greater membrane potential than the wt strain (Fig. 1), and that the wt phenotype was restored in the complemented strain. The higher membrane potential in the ΔMSMEG_3763 null mutant suggests that the pump may be involved in both drug and ion transport, as described for other ABC pumps [21, 22].
MSMEG-3762/63 mediated drug tolerance is specific to rifampicin and ciprofloxacin
The increased efficacy of rifampicin and ciprofloxacin in null mutants of the MSMEG-3762/63 efflux pump was previously established [20]. To extend these findings to other first- and second-line anti-TB drugs, minimum inhibitory concentrations (MICs) were determined using the Resazurin Microtiter Assay (REMA). M. smegmatis mc2 155 (wt) and its isogenic deletion mutant (ΔMSMEG_3763) and complemented (ΔMSMEG_3763 pBD04) strain were grown to log phase and exposed to increasing concentrations of first-line (isoniazid, ethambutol, streptomycin) and second-line (levofloxacin, moxifloxacin, gatifloxacin) anti-TB drugs (Supplementary Table S2). In contrast to rifampicin and ciprofloxacin, no variations in MIC values of the null mutant compared to wt and complemented strains were detected, indicating that the MSMEG-3762/63 pump is not involved in the extrusion of these antibiotics. Therefore, we focused on the role of MSMEG-3762/63 in the efflux of rifampicin and ciprofloxacin.
The MSMEG_3762/63/65 operon is induced on exposure to sub-inhibitory antimicrobial drug concentrations
To investigate the possible role of rifampicin and ciprofloxacin as transcriptional inducers of the MSMEG_3762/63/65 operon, RT-qPCR analyses were performed on M. smegmatis wt cells exposed to sub-inhibitory concentrations of rifampicin and ciprofloxacin (1/3rd MIC). The expression of the operon increased 3.2-fold in the presence of rifampicin (Fig. 2a) and 2.8-fold in the presence of ciprofloxacin (Fig. 2b) compared to drug-free control bacilli. These data suggest that sub-inhibitory concentrations of these two antibiotics induce expression of this efflux pump system, we hypothesized by interfering with the DNA-binding activity of the MSMEG-3765 TetR transcriptional regulator that controls expression of the efflux pump.
The efflux pump operon MSMEG_3762/63/65 is overexpressed in macrophages
Cossu and coauthors [23] reported that the mRNA profiles of M. tuberculosis and M. smegmatis exposed to acid-nitrosative stress, mimicking the macrophage environment, showed a high number of differentially expressed genes. In particular, the MSMEG_3762/63/65 operon was upregulated, and this observation was confirmed by GFP promoter probe analysis [19]. To expand this finding, human THP-1 monocytes were differentiated into macrophages and exposed to M. smegmatis. Phagocytosis was permitted for 2 hours, after which the macrophages were lysed and RNA was extracted from the intracellular mycobacteria. The expression of the MSMEG_3762/63/65 operon was 2-fold higher compared to cells grown in RPMI alone (Fig. 2c). This result supports the hypothesis that the MSMEG-3762/63 efflux pump, and likely homologous Rv1687/86c in M. tuberculosis, are involved in the adaptation of mycobacteria to intracellular environments.
Rifampicin and ciprofloxacin dissociate the TetR-like regulator from the operator region
To further investigate the role of rifampicin and ciprofloxacin in the de-repression of the MSMEG_3762/63/65 operon, EMSA analyses were performed. The assay measured the binding of the MSMEG-3765 TetR recombinant protein to the cognate M. smegmatis 133 bp DNA fragment containing the operator site (36 bp palindromic motif) [19] in the presence and absence of these antimicrobial drugs. The presence of increasing concentrations of rifampicin or ciprofloxacin (Fig. 3a and b, respectively) in the DNA/protein binding reaction mix resulted in an increase of the free unbound DNA fragment in respect to the control mix without drug. Streptomycin, which the null mutant ΔMSMEG_3763 was not hypersensitive to (unlike rifampicin and ciprofloxacin, Supplementary Table S2), was used as a negative control. As expected, inhibition of the MSMEG-3765 binding to the TetR operator sequence with streptomycin was not observed (Fig. 3a, lanes 5 and 6). These data indicate that rifampicin or ciprofloxacin are able to directly inhibit the binding of the MSMEG-3765 TetR protein to its operator in the MSMEG_3762/63/65 operon, suggesting their role as specific ligands for this efflux pump repressor.
The dimer of MSMEG-3765 exhibits the canonical structural architecture of TetR family transcriptional regulators
To understand the biological activities of the MSMEG-3765 repressor, description of its structural and dynamic features is essential, because the shape and nature of the protein surface play crucial roles in protein function. Therefore, we computationally predicted its three-dimensional structure using I-TASSER software (see Mat and Met, and Supplementary methods). The Ramachandran plot in Supplementary Figure S1 reported over 97% residues in most favored and additional allowed regions demonstrating the good quality of the predicted 3D structure. The overall structural model of MSMEG-3765 is an Ω shaped-dimer (Fig. 4a) typical of a TetR family transcription regulator [24, 25]. Each polypeptide chain (residues 2-211) of the MSMEG-3765 homodimer is folded into nine α-helices (α1-α9 for chain A; α1´-α9´ for chain B) with connecting turns and loops. The tertiary structure of the MSMEG-3765 monomer is mainly stabilized by hydrophobic helix-to-helix interactions. As for other members of the TetR family, the global structure of the MSMEG-3765 homodimer can be divided into two DNA-binding domains at the N-terminal tail (NTDs) of each monomer, and a regulatory core at the C-terminal domain (CTD) involved in dimerization and ligand binding. Each NTD is composed of helices α1 to α3 (α1 = Ser21 - Arg37; α2 = Ile45 - Ala50; α3 = Pro56 - His60) (Fig. 4a and b) within which α2 and α3 constitute the classical helix-turn-helix (HTH) motif that is stabilized by the first α-helix (α1). Moreover, helix 1 is preceded by a positively charged region that, as observed for other TetR family members, may have a crucial role in the DNA recognition mechanism [24]. The regulatory domain is formed by helices α5 to α9 (residues α5 = Pro100 - Ser109, α6 = Lys113 - Ile122, α7 = Val128 - Asp148, α8 = Gly154 - Glu176, α9 = Val185 - Gln197) and their symmetric helices α5´ to α9´ and it is connected to the NTDs through helices α4 and α4´. Helices α5 to α7 and their symmetric counterparts α5´ to α7´ form a central triangle subdomain; whereas α8 and α9 and their symmetric helices α8´ and α9´ make up a four-helix bundle that constitute the dimerization interface (Fig. 4a and b). Analysis of the electrostatic surface potential indicated that the NTDs are rich of positively charged residues; on the contrary, the core of the regulatory domain is negatively charged with several hydrophobic pockets (Fig. 4c). We then analyzed the conservation of the MSMEG-3765 dimer structure with the Consurf server using as a query the chain A of the predicted homodimer. As expected, the residues located in the HTH motifs of the NTDs were highly conserved (Supplementary Fig. S2). In contrast, the residues situated at the outer surface surrounding the central triangle of the regulatory domain presented low conservation scores. Interestingly, the most conserved residues of the regulatory core domain belong mainly to helices forming the dimerization interface (α8, α9, α8´ and α9´), helices α5 to α7 and the dyad-related α5´ to α7´ that, as for other TetR regulators [24], may play important roles in ligand recognition and modulation of MSMEG-3765 homodimer binding.
The MSMEG-3765 N-terminal DNA domains are more dynamic than the C-terminal regulatory core domain
In order to understand the relationship between structure, dynamics and function of the MSMEG-3765 homodimer, we investigated the dynamic personality of the dimer using Molecular Dynamics (MD) simulation methods (Fig. 4d and e). These techniques detail the motions and conformational changes of proteins and nucleic acids [26]. We first generated a conformational ensemble describing protein intrinsic motions for the MSMEG-3765 dimer (Fig. 4d); we then analyzed the per-residue RMSF (Root Mean Square Fluctuation) values calculated from the simulated MD ensemble (Fig. 4e). The RMSFs indicate that in each MSMEG-3765 monomer the N-terminal DNA-binding domains (NTD) (residues 1–65) and the C-terminal regulatory domain (CTD) (residues 85–211) are characterized by a different conformational mobility, with the NTDs showing a higher degree of flexibility than CTD (RMSFavgNTD = 7.17 ± 3.27 Å; RMSFavgCTD = 3.25 ± 0.70 Å) (Fig. 4e and Supplementary Table S3). In particular, in both NTDs the region preceding the first α-helix (residues 1–20) presents higher RMSFs than the helices α1, α2 and α3 demonstrating that the N-terminal tail is more flexible than the rest of the domain. Helix α4 and its symmetric counterpart α4´, showing a per-residue averaged RMSF of 2.95 ± 0.39 Å and 3.06 ± 0.45 Å, respectively, are connected to the CTD by two ten-residue relatively flexible loops. Helices α5 to α7 and their symmetric counterparts α5´ to α7´ present higher RMSF values than the helices involved in the dimerization interface, demonstrating that the four-helix bundle is more rigid than the central triangle subdomain.
MSMEG-3765/DNA recognition is regulated by interactions between the HTH motifs and the major groove
To explore the molecular determinants driving the DNA recognition mechanism of MSMEG-3765 homodimer, we performed a series of molecular docking studies. At the N-terminal region of each MSMEG-3765 monomer the DNA-binding domain (NTD) is composed mainly of the first three helices. According to this observation, the structural model of the MSMEG-3765 dimer/DNA complex shows that the core of each NTD is composed of helices α1- α3, with helix α3 responsible for the majority of DNA contacts (Fig. 5a and b). Helix α3 and the symmetric counterpart α3´ recognize mainly major grooves. Interestingly, the majority of MSMEG-3765 dimer/DNA contacts are with the phosphate backbone. Indeed, because the operator sequence is partially palindromic the two NTDs (chain A = 1–65; chain B = 212–276) present a slightly different DNA binding interface producing a slight inclination (~ 15°) of the dimer with respect to the DNA long axis.
In particular, in monomer A, Asn38, Lys42 and Ser44 contact the phosphate backbone of C24, G25 and G26; in helix α2 Arg46 forms hydrogen bonds with the phosphate backbone of G27, C28, C8´and G9´ (symbol (´) indicates the anti-parallel operator strand) whereas Ala50 participates in DNA-binding by making hydrophobic contacts with C8´; Ala57, located within the recognition helix α3, is involved in hydrophobic base specific interactions fundamental for the DNA binding mechanism. In addition, Asp53, Val54 and Asp55 residues, located inside the turn of the HTH motif of the NTD monomer A, contribute to stabilization of the MSMEG-3765 dimer/DNA complex. For monomer B, in helix α3 Ala57´ and Leu58´ are involved in hydrophobic interactions with A17, C18 and T19 whereas His60´ is hydrogen-bonded with the phosphate backbone of A17´, A18´ and G19’; Ile45´ and Arg46´ of helix α3 form hydrophobic and hydrogen bond interactions with G16´of the anti-parallel DNA strand. Moreover, other interactions between the monomer B and the DNA are hydrogen bonds involving several residues located in the turn of the HTH motif and amino acids situated at the edge of helix α4. Notably, different to chain A, the DNA-binding of monomer B is further stabilized by a positively charged region preceding the first helix of the NTD. In particular, the side chain of Arg9 is deeply buried in the minor groove, forming hydrogen bonds with the phosphate backbone of T15, and G16, C24´ and A25´ (Fig. 5b).
Rifampicin and ciprofloxacin bind the MSMEG-3765 repressor at different binding sites
To uncover how rifampicin and ciprofloxacin might interact with the TetR-like regulator to abrogate DNA binding, we integrated these antibiotics into our molecular docking models of MSMEG-3765. The 3D structural model of MSMEG-3765 homodimer/rifampicin indicates that the ligand-binding surface of the homodimer is composed of residues from both monomers located within the C-terminal regulatory domain (CTD). Rifampicin binds into a hydrophobic groove of the protein with specific contacts to helices of the central triangle subdomain (α5 - α7) and the helix α8 (Fig. 6a). In detail, the binding of rifampicin is principally driven by an extensive network of hydrophobic interactions involving: Leu103 and Trp107 located in the middle of the α5; Gln136, Ile139 and Val140 encompassing the helix α7; Ala162, Leu165, Val166 and Val169 situated within the helix α8; Leu175 and Phe180 located in the helix α8´ of the other monomer (chain B). Additionally, the interaction of rifampicin with MSMEG-3765 is further stabilized by a series of hydrogen bonds between the ligand and different MSMEG-3765 residues: Arg120 at the C-terminal edge of the helix α6; Arg132 and Gln135 located inside the helix α7; and Ala162, Val166 of the helix α8 (Fig. 6b and c). In the case of ciprofloxacin, the structure of the MSMEG-3765 homodimer/ciprofloxacin complex shows that the ligand is accommodated into a shallow hydrophobic pocket on the protein surface formed mainly by Leu135 of the helix α7 and by two residues of the helix α8 (Ala162 and Val166).
.
Moreover, the guanidinium of Arg132 (α7) hydrogen binds with both carboxyl oxygens of ciprofloxacin (Fig. 6d, e and f). These latter interactions play an important role in the stabilization of the MSMEG-3765 homodimer/ciprofloxacin complex. Notably, comparison of the structural models obtained for rifampicin and ciprofloxacin clearly shows that the two antibiotics bind the MSMEG-3765 CTD at distinct binding sites suggesting that the two drugs likely present different structural mechanisms that mediate operon de-repression. To provide structural insight into the potential molecular mechanisms of efflux pump operon de-repression mediated by rifampicin and ciprofloxacin, we analyzed the 3D models of MSMEG-3765 in complex with either drug using the structural data available for other TetR family members [24]. This conformation analysis indicated that, upon ligand binding to CTD, de-repression of the operon occurs through conformational changes that abolish the ability of the homodimer to recognize the operator sequence by increasing the separation between NTDs of the two monomers.