Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Biosci., 02 May 2022
Sec. Biological Modeling and Simulation
This article is part of the Research Topic Understanding Biomolecular Folding and Assembly: The Energy Landscape Perspective View all 5 articles

Comparative Structural Dynamics of Isoforms of Helicobacter pylori Adhesin BabA Bound to Lewis b Hexasaccharide via Multiple Replica Molecular Dynamics Simulations

  • Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, India

BabA of Helicobacter pylori is the ABO blood group antigen-binding adhesin. Despite considerable diversity in the BabA sequence, it shows an extraordinary adaptation in attachment to mucosal layers. In the current study, multiple replica molecular dynamics simulations were conducted in a neutral aqueous solution to elucidate the conformational landscape of isoforms of BabA bound to Lewis b (Leb) hexasaccharide. In addition, we also investigated the underlying molecular mechanism of the BabA-glycan complexation using the MM/GBSA scheme. The conformational dynamics of Leb in the free and protein-bound states were also studied. The carbohydrate-binding site across the four isoforms was examined, and the conformational variability of several vital loops was observed. The cysteine–cysteine loops and the two diversity loops (DL1 and DL2) were identified to play an essential role in recognizing the glycan molecule. The flexible crown region of BabA was stabilized after association with Leb. The outward movement of the DL2 loop vanished upon ligand binding for the Spanish specialist strain (S381). Our study revealed that the S831 strain shows a stronger affinity to Leb than other strains due to an increased favorable intermolecular electrostatic contribution. Furthermore, we showed that the α1-2-linked fucose contributed most to the binding by forming several hydrogen bonds with key amino acids. Finally, we studied the effect of the acidic environment on the BabA-glycan complexation via constant pH MD simulations, which showed a reduction in the binding free energy in the acidic environment. Overall, our study provides a detailed understanding of the molecular mechanism of Leb recognition by four isoforms of H. pylori that may help the development of therapeutics targeted at inhibiting H. pylori adherence to the gastric mucosa.

Introduction

The mucus layer of the stomach in humans represents the first barrier for gastrointestinal bacteria. It consists of mucin MUC5AC and MUC6, secreted from the surface mucosa and gland mucosa, respectively (Ho et al., 1995). These secreting mucins are heavily O-glycosylated glycoproteins (Strous and Dekker, 1992; Nordman et al., 1995). One of the most common gastric pathogens colonizing the stomach mucosa is Helicobacter pylori. It affects more than half of the global population (Cave, 1996). It is a Gram-negative, microaerophilic, spiral bacteria that has been identified as a major etiological factor in gastritis, gastric and duodenal ulcers, gastric adenocarcinoma, and gastric-mucosa-associated lymphoid tissue (MALT) lymphoma (Schreiber et al., 2004; Polk and Peek, 2010). It can thrive in harsh conditions of the stomach by the adaptive mechanism of using a variety of outer membrane proteins to adhere with glycan moieties on the gastric epithelium. This mechanism helps H. pylori to sustain local colonization and adherence onto gastric epithelium (Ilver, 1998).

Adhesin of H. pylori to mucin and host cells is highly relevant for the progression of the disease, colonization, and mucin-mediated response (Sheu, 2003; Yamaoka, 2008; Ishijima et al., 2011). Several molecules have been implicated as receptors for H. pylori adhesins (Evans and Evans, 2000). H. pylori adheres to the glycans present at high density in glycolipids, glycoproteins, and mucins of the gastrointestinal tract (GI). Primarily, it interacts via the blood group antigen-binding adhesin (BabA), which binds to mono (ABO) and di-fucosylated (Lewis b; Leb) derivatives of type 1 O-glycan structures. On the other hand, the sialic acid-binding adhesin (SabA) binds to sialyl-Lewis x (SLex) and sialyl-Lewis a (SLea) (Lindén et al., 2002, 2004). The most studied interaction is between BabA and Leb (Ilver, 1998). This di-fucosylated oligosaccharide is abundantly expressed in gastric mucosa, including the O phenotype population. They are most susceptible to peptic ulcer disease (Boren et al., 1993; Lindén et al., 2002). BabA belongs to the Helicobacter outer membrane porins (Hop), adhering to glycan moieties on the gastric epithelium (Alm et al., 2000; Hage et al., 2015).

The H. pylori strain variants determine the longevity of the established infection and disease outcomes (Kang and Blaser, 2006; Suerbaum and Josenhans, 2007). Furthermore, the increased genetic diversity of H. pylori, particularly the high adaptive evolution in the BabA adhesin, affects binding preference to the blood group strains such as the ABO versus O, depicting as “generalist” versus “specialist” preference, respectively (Kang and Blaser, 2006). It is also evident that BabA’s adaptive traits, such as the specialist strains, are more likely to bind to Leb and blood group O dominant populations. In contrast, the generalist strains can bind Leb and the GalNAc- and Gal-forms of blood group A (ALeb) and B (BLeb) derivatives. In particular, H. pylori strains have displayed a 1,000-fold range of affinities against Leb glycan receptors (Ka ∼108 to 1011 M−1) (Aspholm-Hurtig, 2004).

Structurally, BabA contains an N-terminal extracellular host-binding domain and a C-terminal outer membrane domain, forming a β-barrel structure (Alm et al., 2000). As depicted in Figure 1A, BabA also contains two α-helical regions-a handle and a head region. On top of the head region, there is a β-sheet known as the crown. The handle region containing the α and β unit of the extracellular domain, the highly conserved β-barrel transmembrane domain, succeeds the α-C2 helix of the handle region. The core of the head region is primarily comprised of four helices, which are at a perpendicular angle to the handle region, making a markedly kinked tertiary structure. The crown is composed of four antiparallel β-sheets at the highest tip of the protein, which consists of the important Leb binding sites (Doohan et al., 2021).

FIGURE 1
www.frontiersin.org

FIGURE 1. (A) Ribbon representation of BabA protein with the bound Leb (ball and stick representation). Three different regions of the protein, namely, crown, head, and handle, are shown in distinct colors. (B) Percentage probability of the sequence alignment between all four isoforms. (C) SNFG representation of the Lewis sugar b (Leb).

In contrast to BabA, there is a complete absence of a four-strand antiparallel β-sheet of the crown region in SabA, which recognizes sialyl-Lewisx (SLex) found on the cancerous and infected gastric tissues (Pang et al., 2014). The multiple sequence analysis of BabA isoforms shows ∼93% sequence diversity (Figure 1B and Supplementary Figure S1). It has been observed that Leb interacts with the crown of BabA (Hage et al., 2015). Leb is comprised of two fucose residues (Fuc1 and Fuc4), two galactose residues (Gal2 and Gal5), an N-acetyl glucosamine residue (GlcNAc3), and a glucose residue (Glc6) (see Figure 1C). However, each isoform’s structure to molecular determinants and specificity toward Leb is unknown.

This study focuses on understanding the BabA adhesin sequence diversity among different strains (strain 17875; Isoform 1, strain P436; Isoform 2, strain S831; Isoform 3, and strain A730; Isoform 4) toward Leb-type glycan receptor binding and elucidating the biophysical basis of their molecular recognition and conformational dynamics. The atomistic description provided in the study helps to understand structural and binding differences in BabA’s carbohydrate-binding domain (CBD) for each isoform. We investigated the conformational dynamics of various functional loops, such as Cysteine bound loops (CL) and Diversification loops (DL), which exhibit high polymorphism in H. pylori isolates toward Leb binding strength. In addition, we studied the conformational dynamics of Leb and elucidated the glycan-binding mechanism via molecular mechanics generalized Born surface area (MM/GBSA).

The H. pylori infection occurs in the gastric epithelium, where its mucus layer contains pH ranges from neutral to acidic. However, using chemotaxis motions, H. pylori moves from an acidic to a neutral environment. The binding phenomenon is acid sensitive, and the binding affinity reduces with a lowering pH value. To investigate the effect of an acidic environment on the binding, we conducted constant pH MD simulations for Iso1 as a test case. Altogether, this study provides molecular insights into the recognition and binding preference of Leb across H. pylori BabA variants.

Materials and Methods

MD Simulations of BabA–Leb Complexes

The following four isoforms of BabA complexed with Leb were used in the current study: 1) Iso1 (PDB ID 5F7M), 2) Iso2 (PDB ID 5F9A), 3) Iso3 (PDB ID 5F8R), and 4) Iso4 (PDB ID 5F93) (Moonens et al., 2016). We also simulated apo-BabA (without Leb) and Leb in an aqueous solution. The AMBER ff14SB force field (Maier et al., 2015) was used for the protein modeling, and glycam06j-1 (Kirschner et al., 2008) was employed to describe Leb. The protonation state of each residue was estimated using the Propka 3.1 webserver (Olsson et al., 2011), corresponding to pH seven for all four different Isoforms of BabA. Each system was solvated in a truncated octahedral box using the TIP3P (Price and Brooks, 2004) water molecules. The distance between the wall and solute was kept as 10 Å. Each solvated system was neutralized by adding an appropriate number of counterions. All bonds, including hydrogen atoms, were restrained using the SHAKE algorithm (Kräutler et al., 2001). The long-range electrostatic interactions were treated using the particle mesh-Ewald (PME) method (Darden et al., 1993). A cut-off distance of 10 Å was set for non-bonded interactions. The temperature was maintained using the Langevin thermostat with a collision frequency of 2 ps−1 and the pressure using the Berendsen barostat (Berendsen et al., 1984). Two stages of minimization were performed using the steepest descent followed by the conjugate gradient algorithm for 500 steps with a harmonic restraint of 2 kcal mol−1 Å−2. It was then subjected to gradual heating from 0 to 300 K with a restraint force of 2 kcal mol−1 Å−2. Before the production run, we conducted a 1 ns simulation at the NPT ensemble without restraint on the solute for equilibration. Finally, the production simulations for complexes and apo-BabA were performed in the NPT ensemble with a timestep of 2 fs for 1 µs using the pmemd.cuda module of AMBER18(Case et al., 2018). We recorded the coordinates every 10 ps, resulting in 100,000 configurations. To ensure better sampling, we conducted two more replica simulations for apo and complex systems. We also conducted 2 × 3 µs long simulations for the free glycan in the aqueous solution using a similar protocol as discussed in our previous work (Roy et al., 2022).

In addition, to investigate the role of an acidic environment on the complexation of BabA/Leb, we conducted constant pH molecular dynamics (CpHMD) simulation at pH 4.5 for 2 × 250 ns for Iso1 by titrating acidic amino acids in the crown region (Donnini et al., 2011). In the selected region, six ionizable residues were found (E80, E158, D161, D172, D199, and D213), which we attempted to titrate during the simulation. After every 180 fs of explicit MD simulation, the protonation change event was attempted for all mentioned residues in the GB implicit model (Mongan et al., 2004).

Analysis

All the trajectories were analyzed using the cpptraj module of AmberTools19 (Roe and Cheatham, 2013). The first 200 ns data from each run were discarded for the trajectory analysis to minimize initial noises. After that, the combined trajectories for three replicas were prepared and subjected to various calculations related to the dynamics of the systems. The initial stability and flexibility were estimated by calculating the root mean squared deviations (RMSD) and the root mean squared fluctuations (RMSF) concerning the well-equilibrated conformations. Conformational space of the sugar molecules in the free and bound state was also estimated with the help of glycosidic dihedral angle, Cremer–Pople (CP) (Cremer and Pople, 1975) parameters of ring puckering. The glycosidic dihedral angles for the α/β-n linkages are defined as φ=O5C1OCn and ψ=C1OCnC(n1). The puckering state of any sugar as per the Cremer–Pople convention was estimated using the spherical coordinate set, Q, θ, and φ. The spherical space was further subdivided into 38 IUPAC recognized regions, where each space is dedicated to a specific ring conformation. To estimate the global and local conformations of several loops related to binding, the Cartesian principal component analysis (PCA) and dihedral principal component analysis (dPCA) were performed, which were discussed in detail in our previous study (Roy et al., 2021). The coordinate of Cα atom of each amino acid and dihedral angle of peptide bonds of the selected loop region were used for PCA and dPCA, respectively. Furthermore, the free energy surface (FES) was constructed using the Boltzmann equation: ΔG = −kBT ln(ρ) - k, where kB is the Boltzmann constant, T is the temperature, and ρ is the probability density of the geometric coordinate x. The parameter (k) was chosen such that the global minimum was located at 0 kcal/mol. In addition, a distance cut-off of ≤3.5 Å and an angle cut-off of ≥120° were used to evaluate the hydrogen bonding (Jeffrey, 1997).

Binding Free Energy Calculations Using the MM/GBSA Scheme

The molecular mechanics/generalized Born surface area (MM/GBSA) scheme (Onufriev et al., 2004) was employed to estimate binding free energy. For protein-glycan systems, it has been reported that MM/GBSA outperforms the computationally expensive MM/PBSA scheme (Mishra and Koča, 2018; Sood et al., 2018; Roy et al., 2020a). It offers a good trade-off between speed and accuracy compared to other computationally expensive methods, such as thermodynamics integration (TI) and free energy perturbation (FEP). The methodological details of the MM/GBSA protocol were discussed in our previous works (Kar et al., 2013; Roy et al., 2020b; Jonniya et al., 2021). The binding free energy (ΔGbinding) was determined from the free energies of the receptor/BabA (protein), the ligand (carbohydrate), and the complex (complex) according to the following equation (Kollman et al., 2000):

ΔGbinding= Gcomplex(Gcarbohydrate+ Gprotein)(1)

The free energy of each species (protein, carbohydrate, and complex) was estimated using the following formula:

G=EvdW+Eele+GGB+GSA(2)

where EvdW and Eele represent the van der Waals and electrostatic interaction energy, respectively. The polar solvation free energy (GGB) was estimated using the generalized Born equation, while the non-polar solvation free energy (GSA)  was determined using the following equation:

GSA= γA+b(3)

with γ=0.00542 kcal.mol1.Å2 and b=0.92 kcal.mol1. The symbol A denotes the solvent-accessible surface area (SASA). The GB model (igb = 2) developed by Onufriev et al. was used in this study (Onufriev et al., 2004). The configurational entropy was not considered in the current study due to the high computational cost. We considered 25000 frames from the last 500 ns of each production simulation for the free energy calculation. The overall binding free energy was estimated by averaging the binding free energies obtained from the three independent runs. In addition, we performed a per-residue decomposition of the binding free energy by the MM/GBSA scheme. On the other hand, 7500 frames from the last 150 ns of the trajectory were employed for the binding free energy calculation at pH 4.5. The final values were determined by averaging both replica runs.

Results and Discussions

Structural Stability and Flexibility Analyses in Neutral Aqueous Solution

First, we investigated the structural stability of each system by analyzing the root mean squared deviations (RMSD) with respect to the corresponding crystal structure. The time evolution of RMSD for multiple runs is presented in Supplementary Figures S2-S3. It demonstrates that simulations are stable for each BabA isoform in its apo (see Supplementary Figure S2) and glycan-bound forms (Supplementary Figure S3). The handle region of BabA was more flexible as it was not fully crystallized, which led to an increase in the overall RMSD. Next, we determined the RMSD of each BabA–Leb complex isoform by excluding the handle region and computed the corresponding probability distribution by combining all three replica runs (see Figure 2). In addition, each complex was compared with the respective free form. The average RMSD values for multiple runs were computed and listed in Table 1. It is evident from Table 1 that the average value varies between 2.96 Å and 3.60 Å for all apo isoforms, while the RMSD value decreases upon complexation, with the values ranging from 2.27 to 3.16 Å for all four complexes. It is evident from Figure 2A-D and Table 1 that compared to apo, BabA–Leb complexes got relatively stabilized. However, for Iso4, a slightly higher RMSD was observed. These results indicate that the interaction between Leb and BabA leads to overall structural stability. All complexes except Iso3 showed similar distribution for apo and bound states. In Iso3, the apo form showed a stable conformation across all three independent runs, which got relatively flexible after the binding of Leb. In the apo case, the distribution was narrow with a sharp peak, while a broader distribution was observed when a complex was formed. However, the peak shifted to a lower RMSD value compared to apo Iso3.

FIGURE 2
www.frontiersin.org

FIGURE 2. Distribution of root mean squared deviation (RMSD) of the crown and head region. (A) Iso1 (Apo and Complex), (B) Iso2 (Apo and Complex), (C) Iso3 (Apo and Complex), (D) Iso4 (Apo and Complex), and (E) free and bound Leb.

TABLE 1
www.frontiersin.org

TABLE 1. The average root mean squared deviation (RMSD), radius of gyration (RoG), and solvent accessible surface area (SASA) over three replica runs using the block average method. The standard error of the mean is given in parentheses.

Next, we calculated the RMSD of the ring atoms of Leb in the free and protein-bound states concerning the initial coordinates. The corresponding RMSD distribution is displayed in Figure 2E. A broad distribution was observed for the free glycan with a mean value of ∼0.9 Å. In the case of the Iso3 complex, a very sharp unimodal distribution with a peak at ∼0.5 Å was observed, indicating a firm binding. A similar distribution was observed for Iso2/Leb. On the other hand, moderate flexibility was observed for the glycan when bound to Iso1 and Iso4. In the case of Iso1, the most probable peak is located at 1.0 Å, with a secondary peak at ∼2.6 Å. In the case of Iso4, the RMSD distribution is characterized by two peaks, situated at ∼0.5 Å and ∼1.1 Å. The degree of flexibility of the glycan in the protein-bound form may reflect its binding affinity, as discussed in the later section.

Next, the average values of the radius of gyration (RoG) were computed to check the compactness of the protein and are listed in Table 1. The time evolution of RoG of each isoform in free and ligand-bound states is shown in Supplementary Figures S4, S5. It is evident from Supplementary Figure S5 that in the case of Iso4, run2 showed distinct time series compared to other runs. It is worth mentioning that the deviation in RoG was mostly because of the flexible handle region, which is located far from the binding region. Overall, the calculated average RoG values indicated no significant changes for each isoform. Next, we computed the SASA and found that the SASA value of apo variants was more significant than their respective complex form, as seen in Table 1 and Supplementary Figures S6, S7. The distribution of RoG and SASA for apo and complexes are shown in Supplementary Figure S8. These observations agree with the RMSD profile, suggesting that non-polar amino acid residues form interactions during complexation with the glycan. In addition, it suggests that the overall stability of BabA is enhanced upon complexation with the glycan, indicating that the hydrophobic core of BabA is more protected from the external environment upon binding with Leb (see Supplementary Figure S8).

Finally, to provide further insights into the flexibility of different regions of BabA, the root means square fluctuations (RMSF) of Cα atoms were computed. The difference RMSF (∆RMSF) plot of the complex from the respective apo conformation was computed and shown in Figure 3. The N-terminal domain of BabA was relatively more flexible, and it corresponds to the handle region. However, this region is located far from the binding site and can be ignored for our current study. Furthermore, it is evident from Figure 3 that among different significant loop regions of BabA, DL2 of apo Iso3 exhibited large fluctuations compared to the corresponding complex. This indicates that the DL2 loop underwent structural changes upon interaction with Leb and could play a crucial role in the molecular specificity of Leb across isoforms of the H. pylori adhesin BabA.

FIGURE 3
www.frontiersin.org

FIGURE 3. (A) Difference in root mean squared fluctuation (RMSF) of four isoforms with their respective Apo conformation. Different loop regions are highlighted in the plot. (B) Location of all key loop regions in BabA.

Overall, Iso3 showed differences in dynamics and flexibility upon complexation, especially near the binding site, compared to the other three isoforms. This structural change may influence the Leb recognition process.

Conformational Dynamics of Free and Protein-Bound Leb in Neutral Aqueous Solution

The flexible glycosidic linkage of Leb yields several conformations, which may play a critical role in protein-glycan recognition. The flexibility was investigated by calculating the free energy surface (FES) using φ/ψ angles of each glycosidic linkage and shown in Figure 4. The flexibility of each linkage stays similar for free and bound states. Although the conformational sampling across the dihedral space remained almost unchanged after Leb binding, a minor variation for specific linkages, such as [Gal (5)β1-4Glc (6)] and [GlcNAc(3)β1-3Gal (5)], was observed. Both linkages displayed higher flexibility compared to others. Furthermore, linkages such as Fuc (4)α1-4GlcNAc(3), Fuc (1)α1-2Gal (2), and Gal (2)β1-3GlcNAc(3) stayed at a single minimum along with a very narrow conformational basin irrespective of free or bound states. The GlcNAc(3)β1-3Gal (5) linkage showed two equiprobable minima around (−80, −150) and (−80, −80), which got diminished after ligand binding. In the bound state, the flexibility of Leb remained almost similar across all isoforms, indicating a stable conformation around the binding pocket. However, in the case of Iso4, the [Gal (5)β1-4Glc (6)] linkage showed two closely separated minima regions, indicating a slightly different conformation.

FIGURE 4
www.frontiersin.org

FIGURE 4. 2D free energy surface constructed using dihedral angles φ/ψ of glycosidic angle in the free and bound form of Leb.

Another important parameter for characterizing the carbohydrate conformation is the ring geometry, for which we evaluated the puckering profile for each monosaccharide. We used the Cremer–Pople parameter (θ) to construct the conformational free energy to show the transition from 4C1 to 1C4. In our study, no monosaccharide was able to complete the chair (4C1)-to-chair (1C4) transition, except for glucose (6) and galactose (5) in the free system (see Figure 5). However, other galactose, such as Gal (2), stayed around the 4C1 conformer, indicating a similar rigidity as other studies (Alibay et al., 2018). In the case of β-glucose, the complete transition was observed in the free simulation where the 1C4 state is energetically much higher than the other chair form 4C1 (∼2 kcal/mol). On the contrary, in the case of Iso4, both the pucker states were sampled in an equiprobable way. The flexible nature of this ring conformation was also reflected in the dihedral space, as shown in the FES of the [Gal (5)β1-4Glc (6)] linkage (Figure 4). An earlier study also reported a change in the puckering state for glucose in a microsecond timescale (Alibay and Bryce, 2019). Fucose residues stayed only in the 1C4 state, as observed in another protein–carbohydrate recognition study (Roy et al., 2020a). Overall, Leb displayed similar flexibility across all four isoforms.

FIGURE 5
www.frontiersin.org

FIGURE 5. 1-D free energy surface generated using Cremer–Pople θ puckering profiles of each monosaccharide of Leb.

Conformational Flexibility of the Crown Region

A thorough structural investigation was conducted to gain insights into the dynamicity of the crown region of BabA, which mainly interacts with the glycan molecule. This region is comprised of several key loops, such as CL1-3 and DL1-2, which are the key functional regions for the glycan binding. The principal component analysis (PCA) of BabA was conducted, and the free energy landscape was generated using the first two principal components (see Figure 6). The apo structures of all four BabA isoforms showed several co-existing minima connected via shallow energy barriers, suggesting a very dynamical system in solution. However, upon Leb binding, these regions shrunk to a certain extent, indicating the stability of the protein–carbohydrate complexes, especially for Iso1 and Iso2. In addition, for Iso4, both in apo and bound states, several minima regions were separated by a high energy barrier (more than 4 kcal/mol). However, the Cartesian coordinate-based PCA fails to provide a correct separation of internal and overall motions. The projection on the low-dimensional space offers only limited information about the multidimensional energy landscape (Krivov and Karplus, 2004).

FIGURE 6
www.frontiersin.org

FIGURE 6. Free energy analysis of complete BabA with respect to Cα principal components, where PC1 and PC2 represent the first and second principal components of motion. The first row represents the Apo cases, (A) Iso1, (B) Iso2, (C) Iso3, and (D) Iso4, whereas the second row represents the complex simulations, (E) Iso1, (F) Iso2, (G) Iso3, and (H) Iso4.

We also generated the porcupine plot corresponding to PC1 displaying the directions of movements (see Figure 7). In Figure 7, the arrow corresponds to the direction of the movement, while the length of the arrow represents the movement’s strength. Except for the terminal region, most regions showed stability or minor movements. However, the crown region or the carbohydrate-binding region displayed the least movement for apo and complexes except for Iso1. In the case of apo Iso1, DL2 showed a movement toward the binding region. However, in the case of the bound state, the direction of the movement was completely reversed. A similar trend was observed for apo Iso3. Apart from the DL2 loop, some parts of the CL3 loop showed upward movements or toward the crown region. However, this cysteine–cysteine loop region is far from the binding region. So, the flexible movement of this region may not affect the glycan binding.

FIGURE 7
www.frontiersin.org

FIGURE 7. Porcupine plot corresponding to PC1 obtained by performing principal component analysis on MD trajectories, (A) Iso1 (Apo), (B) Iso2 (Apo), (C) Iso3 (Apo), (D) Iso4 (Apo), (E) Iso1 (Complex), (F) Iso2 (Complex), (G) Iso3 (Complex), and (H) Iso4 (Complex).

Further, the structural stability of the crown region was investigated by estimating the number of intra-residue hydrogen bonds (H-bonds). We evaluated the occupancy of H-bonds of the crown region using the Bridge2 software (Siemers and Bondar, 2021). The H-bonds heatmap corresponding to the crown region is shown in Figure 8, where only H-bonds with more than 70% occupancy are listed. For all cases, two to three stable H-bonds were observed, which may determine the structural integrity of the region. The key residues obtained from this analysis are primarily located in the DL2 loop and its surrounding region.

FIGURE 8
www.frontiersin.org

FIGURE 8. Heatmap representation of the possible H-bond paths of the crown region of BabA using Bridge2 (Siemers and Bondar, 2021). Here, the residue pairs with more than 70% occupancy were listed. (A) Iso1 (Complex), (B) Iso2 (Complex), (C) Iso3 (Complex), and (D) Iso4 (Complex).

Furthermore, we also estimated the joint occupancy of H-bonds, which reports true when all the estimated H-bonds are present and false otherwise (Siemers and Bondar, 2021). The joint occupancy for Iso2 was 83.9%, the highest among all isoforms (Iso1 ∼61.5%, Iso3 ∼56.3%, and Iso4 ∼57.3%). Despite the high occupancy of different hydrogen bond pairs, the DL2 region was more flexible across several isoforms of BabA.

Furthermore, to provide a deeper understanding of the dynamicity of significant loop regions, an inspection of loop regions was conducted by estimating the RMSD. The time evolution of all three cysteine–cysteine loops (CL1, CL2, and CL3) RMSD was estimated over all three MD runs and shown in Supplementary Figures S9–S11. In all three replicas, no major conformational changes were observed. This behavior was also true for the diversity loop 1 (DL1), as shown in Supplementary Figure S12. We performed the backbone dihedral angles principal component analysis (dPCA) to extract more information about the conformational changes. CL1 and CL3 loops, away from the binding region, showed a moderate conformational basin across all four isoforms (see Supplementary Figures S13 and 14). For apo structures, the principal minima of CL1 were separated from the rest of the conformational spaces by a high energy barrier. However, after the complex formation, these regions were connected by the shallow free energy barrier. On the other hand, for CL3, the conformation variability remains unchanged before and after the complex formation. As these regions are far apart from the carbohydrate-binding site, the dynamics of these regions may not interfere with Leb binding.

Among all the cysteine-cysteine and diversity loops, the major loops are CL2, DL1, and DL2, which play a crucial role in ligand binding (Moonens et al., 2016). For CL2, the free energy surface was constructed using the first two dPCs, as shown in Figure 9. For all the cases, the free energy surface was characterized by two equiprobable minima, separated by a high energy barrier (Figure 9). However, the minima regions shrunk in the bound state compared to the apo form. The key change was observed in Iso3, where the apo BabA was sampled over four different areas separated by a high energy barrier. However, this intrinsic flexibility diminished after binding with Leb, indicating a better ligand binding phenomenon. Furthermore, to investigate the conformational dynamics of DL1 and DL2 in detail, dPCA was conducted, and the resultant FES is shown in Supplementary Figures S15 and S26. The DL1 loop showed a very stable profile throughout the replica runs, as is evident from the time series of RMSD (Supplementary Figure S12). However, flexibility was observed for all apo and complexes. Among all isoforms, Iso3 showed a more stable conformational space, and the corresponding free energy surface was characterized by two principal minima in the apo and complex forms. It should be noted here that the location of the minima region shifted after the complexation.

FIGURE 9
www.frontiersin.org

FIGURE 9. 2D free energy surface (FES) generated using the first two dihedral PCA (dPCA) components of the CL2 region. (A) Iso1 (Apo), (B) Iso2 (Apo), (C) Iso3 (Apo), (D) Iso4 (Apo), (E) Iso1 (Complex), (F) Iso2 (Complex), (G) Iso3 (Complex), and (H) Iso4 (Complex).

In DL2, the flexibility was observed only in cases of apo forms of Iso1 and Iso3, which decreased after the binding of Leb, as evident from Supplementary Figure S16A and C (for apo cases) and Supplementary Figure S16E and G (for complexes). Similar behavior was observed from the RMSD plot of the DL2 loop in apo and complex forms (see Figure 10A). In apo cases, all three replica runs showed flexibility over the entire 1.0 µs for Iso1 and Iso3. The RMSF analysis (Figure 3) also supports the same observation.

FIGURE 10
www.frontiersin.org

FIGURE 10. (A) Time evolution of root mean square deviation (RMSD) of DL2 loop over three replicas run. Apo and complex cases are shown in the first and second column, respectively. (B) Time evolution of secondary structure content (Beta sheet and turn) of DL2 segment in Iso1 (Apo and complex). Structural superimposition of BabA from the mentioned regions for (C) Iso1 and (D) Iso3.

To investigate the reason behind the shift in RMSD, we conducted the DSSP analysis using the cpptraj module of AmberTools19 (Roe and Cheatham, 2013), where specifically a change in β-sheet and turn was observed. The time evolution of the secondary structure in Iso1 for three replica runs is shown in Figure 10B. The flexible behavior of the DL2 loop in apo cases resonates with the percentage change in the turn of this segment. However, DL2 did not form any β-sheet in the apo simulation of BabA isoforms. In contrast, after the complexation with Leb, a change in coordination of sheet and turn was observed, agreeing with the RMSD profile. Specifically, when the RMSD was stable, a steady percentage of turn was observed, which diminished with the fluctuation. Contrarily, we observed the appearance of β-sheet with the increase in RMSD. Therefore the occurrence of β-sheet was inversely proportional to the decrease in the percentage of turn in the loop structure. The location of the β-sheet could be seen in the representative structure from the respective segment of RMSD (Figure 10C, Iso1). Similarly, the 3D conformation was also shown for Iso3 to indicate the change in the structure of the DL2 loop region (Figure 10D).

Overall, our study revealed that the intrinsic flexible nature of Iso1 persists after the glycan binding. In the case of apo Iso3, the flexibility of the crown region was relatively higher than the other three isoforms, which got drastically reduced after the complex formation. It indicates a stable Leb binding for Iso3 compared to the other isoforms.

The Binding Energetics Behind the Molecular Recognition of Leb Toward Variants of BabA

To elucidate the binding mechanism of Leb against variants of BabA, we calculated the binding free energy for each complex and its different components (see Table 2). We calculated the average binding free energy over the three replica runs and found −37.70, −34.96, −40.26, and −38.28 kcal/mol for Iso1, Iso2, Iso3, and Iso4, respectively. Overall, it suggests that Leb has a stronger affinity toward Iso3 of BabA than the other three isoforms. In contrast, Iso2 showed the lowest affinity for Leb (see Figure 11A). It is to be noted that the electrostatic contribution (ΔEelec) was more favorable than the van der Waals (ΔEvdW) interactions for all complexes. ΔEelec varied between −46.72 kcal/mol and −76.25 kcal/mol, while ΔEvdW varied between −41.77 kcal/mol and −44.71 kcal/mol. The total favorable interaction energy (ΔEMM), which comprises ΔEvdW and ΔEelec, was found to be more favorable for Iso3 (ΔEMM = −120.96 kcal/mol) than other complexes. The lowest value was obtained for Iso2 (ΔEMM = -89.65 kcal/mol). In the case of Iso3, both ΔEvdW (−44.71 kcal/mol) and ΔEelec (−76.25 kcal/mol) were more favorable than other complexes. Although ΔEelec was more favorable to the complexation compared to ΔEvdW, the net electrostatic component, i.e., the sum of intermolecular electrostatic and polar solvation free energy (ΔEelec+ ΔGpol), was unfavorable to Leb binding for all complexes. It suggests that the intermolecular van der Waals interactions mainly governed the binding of Leb. Overall, for Iso3, both components (ΔEelec and ΔEvdW) favor the BabA–Leb complexation the most compared to other isoforms, resulting in a better affinity for Leb.

TABLE 2
www.frontiersin.org

TABLE 2. The binding free energies of all four complexes along with their different energy components from the MM/GBSA scheme in kcal/mol. Standard deviation is given in parentheses.

FIGURE 11
www.frontiersin.org

FIGURE 11. (A) Energetic components of binding free energy for all four complexes using MM/GBSA scheme. (B) Percentage contribution of each monosaccharide for all four isoform complexes.

It is worth noting here that the binding free energy across the triplicate runs varied to a certain extent for all except Iso2. To investigate the origin of variation, we calculated the RMSF of the bound Leb (Supplementary Figure S17). It is evident from Supplementary Figure S17 that in the case of Iso2, Leb showed the least flexibility compared to the other three isoforms. Consequently, Iso2 showed almost similar binding free energy in all three replicas. However, Iso2 showed the lowest binding free energy among all four cases. It may indicate the fluctuation of a few monosaccharides can influence the protein-glycan binding free energy. However, the key residue, such as fucose, displayed the least flexibility across all four isoforms that govern the protein-glycan recognition.

The percentage contribution of each monosaccharide to the total binding was plotted and shown in Figure 11B. Fuc (1) and GlcNAc(3) contributed more favorably to binding than other monosaccharides for all complexes. Furthermore, our study revealed that Fuc (1) contributed ∼40–53% to the total binding free energy among monosaccharides in each isoform, and GlcNAc(3) contributed ∼15–27%. In Iso2, Fuc (1) contributed the highest among other isoforms, accounting for 53%, while the contribution of GlcNAc(3) was the lowest (∼15%) compared to other isoforms. It agrees with a previous study where it was reported that Fuc and GlcNAc contributed more toward the total binding free energy that favors the cholera toxin (CT)/LewisY complexation (Roy et al., 2020a).

Per-Residue Decomposition of ΔGbind

Next, we computed the contribution of each amino acid of BabA to the total binding free energy to find out the hotspot residues. All amino acids contributing more than −1.5 kcal/mol of energy are listed in Supplementary Table S1 and graphically shown in Figure 12. We noted that the highest contributing residues toward Leb binding (>−3.5 kcal/mol) come from Iso3 and Iso4 in agreement with the binding energetics given previously, especially the residue S440 from Iso3 and S247 from Iso4. Similarly, S244 of Iso1 contributed high in the binding with a -3.10 kcal/mol value. In contrast, in the case of Iso2, the highest contributing residue, A240, showed a contribution of -2.52 kcal/mol, which also agrees with the lesser binding affinity than others toward Leb. The other contributing residues involved in the binding included T246, G191, V243, N194 (Iso1), A240, Q235, A241, T243, G191 (Iso2), V243, G191, T246, N194 (Iso3), and G192, T249, Q246, N195 (Iso4). Overall, these results suggest that primary interacting residues for recognizing Leb come from CL2 (G191 and N194) and DL2 (V243, Q235, and S244) loops, which agrees with the above mentioned RMSF and free energy surface analysis. Moreover, these results are consistent with the experimental findings in explaining the significant role of CL2 and DL2 loops in recognizing glycans for the attachment of BabA of H. pylori (Moonens et al., 2016).

FIGURE 12
www.frontiersin.org

FIGURE 12. Per residue-wise contribution to the binding free energy of BabA and Leb complexes: (A) Iso1, (B) Iso2, (C) Iso3, and (D) Iso4. Interaction profile of glycan and protein in the binding site location where glycans and proteins are shown in ball and stick model and ribbon representation, respectively, (E) Iso1, (F) Iso2, (G) Iso3, and (H) Iso4.

In addition, we calculated the contributions of each monosaccharide to ΔGbind for all complexes and listed them in Supplementary Table S2. Fuc (1) contributed more favorably to the protein-glycan binding than other monosaccharides for all complexes. The contribution of Fuc (1) ranged from -7.65 to -8.33 kcal/mol and was more favorable for Iso3 and Iso4 than Iso1 and Iso2. In addition, GlcNAc also contributed more favorably to the binding. The contribution of GlcNAc was found to be the highest for Iso3 (-5.33 kcal/mol), followed by Iso4 (-5.11 kcal/mol), Iso1 (−4.77 kcal/mol), and Iso2 (−2.25 kcal/mol). However, in Iso2, Gal (1) contributed more significantly (−3.19 kcal/mol) than other monosaccharides. Therefore, it suggests that GlcNAc is the next significant monosaccharide favoring the association between Leb and BabA. Furthermore, it is revealed from Supplementary Table S2 that the electrostatic contribution (Tele) of Fuc (1) to the total binding was more favorable than TvdW for each complex ranging from -12.04 kcal/mol to -19.83 kcal/mol. However, Tele was overcompensated by the disfavourable polar solvation of energy (TGB) in each complex. Therefore, the key contributing component of Fuc (1) was vdW (TvdW) which favored the binding.

Subsequently, we computed the H-bond occupancy formed between Leb and BabA for all isoforms, as listed in Supplementary Table S3. Supplementary Figure S18 shows that isoforms 1, 3, and 4 showed high H-bond frequency compared to Iso2. The critical BabA–Leb H-bond interactions are noted in Supplementary Table S3. It suggests that Fuc (1) formed multiple H-bonds with BabA residues (G191, T246, N194, and C189), with the occupancy varying between ∼67 and ∼90% for Iso1, while it ranged between ∼62 and ∼88% for Iso3. Similarly, T249, C190, G192, and N195 from Iso4 interacted with Fuc (1) via H-bond, with the occupancy varying between ∼56 and ∼86%. In the case of Iso2, Fuc (1) formed H-bond with C156, N161, G158, and T210, with an occupancy ranging from ∼40 to ∼87%. Compared to other isoforms, the number of residues forming H-bonds with Leb was less in Iso2. It is evident from Supplementary Figure S18A that the average number of H-bonds in isoforms 1, 2, and 3 was estimated as 8, while on average, 6 H-bonds were observed for Iso2.

Further, we computed the hydrophobic contacts between BabA–Leb as shown in Supplementary Figure S18B. The average value of the hydrophobic contact was the highest for Iso2 (∼12), while for other isoforms, on average, five to six contacts were formed. Altogether, it suggests that the electrostatic interactions favored the association between Leb and BabA for isoforms 1, 3, and 4, while the hydrophobic interaction favors the complexation in Iso2.

The abovementioned analysis was complemented by plotting the interaction diagram of BabA/Leb for all four complexes (see Supplementary Figure S19). As seen from Supplementary Figure S19C, in Iso3, Fuc (1) formed H-bond interactions with G191, C189, T246, and N194. Similarly, GluNAc(3) formed H-bonds with D233, Gal (5) with S244, and Fuc (4) with N208. The spike representation was used to depict the hydrophobic interactions. Overall, Supplementary Figure S19A–D indicate fewer interactions in Iso2 than other isoforms, which agrees with the aforementioned binding affinity. Overall, the interaction profile suggests a strong hydrogen bonding and hydrophobic interactions between BabA and Leb for Iso3 compared to other isoforms of the BabA family, which could unravel the molecular specificity of hexasaccharide Leb against isoforms of H. pylori adhesin BabA.

Influence of Acidic Environment in the Binding Mechanism

Finally, we estimated the binding affinity of Iso1 (17875 strain) at acidic pH (4.5) to elucidate the effect of pH on glycan binding. The RMSD of BabA for both runs was estimated and shown in Supplementary Figure S20. Both runs attained equilibrium after 100 ns, and the last 150 ns were used to calculate the binding free energy using the MM/GBSA protocol. The binding affinity of Iso1 against Leb in the acidic condition was estimated as −34.20 kcal/mol, which is ∼4 kcal/mol lower than what was evaluated in neutral solution (−37.70 kcal/mol). This finding was also supported by the recent experimental results where BabA senses the change in pH in the environment (Bugaytsova et al., 2017). In acidic conditions, the binding affinity of BabA against Leb was dramatically reduced for the strain 17875, agreeing with our MM/GBSA calculation.

Details of the binding free energy and its components are listed in Supplementary Table S4. An increase in the polar solvation energy and decrease in the van der Waals component mainly lowered the binding free energy. However, Iso1 lacks the key salt bridge interaction between D198 and R207, which plays a significant role in acid sensitivity; still, it senses the change in environment as expected by the experimental findings (Hatakeyama, 2017). The short α helix, or the key-coil (residue 199-202) region located in the DL1 loop, is present in Iso1. However, the central residue in the key-coil region for acid sensitivity and binding affinity, i.e., P199, was not present in the 17875 strain, which may reduce the sensitivity against the pH change (Bugaytsova et al., 2017). The critical residues involved in the binding are shown in Supplementary Figure S21, indicating a similar pattern relative to the neutral condition. Irrespective of pH, the key residues are the same, with a minor variation determining the strength of the binding. In acidic conditions, direct contributions from the amino acids of the key-coil region were not observed.

Conclusion

Herein, we performed all-atom multiple replicas molecular dynamics simulations to elucidate the recognition mechanism of the O blood group antigen (Leb) by BabA of H. pylori. This study involved the O specialist (Iso3) and ABO generalist strains of H. pylori. We investigated the conformational preferences of key amino acids present in the crown region of BabA upon Leb binding. Our study revealed no significant change in the flexibility of Leb in the free and BabA-bound states, which was evident from the ring puckering profile. In addition, we showed that the difference in the recognition process between Leb and BabA was due to the distinct movement of loops (CL2, DL1, and DL2) near the glycan-binding sites. Among all isoforms, Iso1 (strain 17875) showed the maximum motions near the DL2 loop, which persisted after the glycan binding. Interestingly, despite the flexibility of the DL2 loop, glycan stays close to the binding pocket due to the changes in secondary conformations. For Iso3 (strain S831), the outward movement of the DL2 loop vanished upon ligand binding, stipulating a better binding partner. Furthermore, structural changes in CL2 were observed for the S831 strain (Iso3) upon glycan binding. This implies that the recognition process must rely on the receptor’s conformational dynamics, which were investigated by characterizing several key loops, such as CL2, DL1, and DL2. Estimating the conformational flexibility of the crown region indicates the stable 2–3 H-bonds located mainly in the DL2 region, which determines the binding region’s structural integrity. In our study, the flexibility of the CL2 loop was drastically reduced for Iso3, which indicates an increased binding affinity. The higher flexibility of the DL2 loop in the bound state for Iso1 results from the interplay between beta-sheet and turn conformations.

To quantify the conformational changes with the binding affinity of Leb, we estimated the binding free energy using the MM/GBSA method. It ranks Iso3 as a strong binding partner, whereas Iso2 shows the least binding affinity among all four isoforms. In the case of Iso3, the stabilization of the critical loops helps achieve a better binding affinity than other isoforms. As this strain is O specialist in nature, better binding supports the experimental characteristics of this strain. The Fuc (1) residue, linked by the α1-2 linkage, provides a significant contribution to the binding. The role of this specific fucose was also blended with the crystallographic information of the BabA–Leb interaction (Moonens et al., 2016). In addition, critical residues involved in the binding, i.e., G191, N194, V243, Q235, and S244, mainly lie in the CL2 and DL2 loop regions implying the importance of these regions in the recognition process. Finally, we investigated the effect of pH on the BabA–Leb association by performing constant pH MD simulations, which revealed a reduction in the binding free energy in an acidic environment. Overall, our study provides a detailed understanding of the interaction of O blood group Lewis hexasaccharide across different isoforms of BabA.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

PK conceived and supervised the project. RR, NAJ, and MFS conducted molecular dynamic simulations and data analyses. RR and NAJ wrote the manuscript. PK edited the manuscript.

Funding

The Department of Science and Technology (DST), Govt. of India (grant number DST/NSM/R&D_HPC_Applications/2021/03.18).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.852895/full#supplementary-material

References

Alibay, I., and Bryce, R. A. (2019). Ring Puckering Landscapes of Glycosaminoglycan-Related Monosaccharides from Molecular Dynamics Simulations. J. Chem. Inf. Model. 59, 4729–4741. doi:10.1021/acs.jcim.9b00529

PubMed Abstract | CrossRef Full Text | Google Scholar

Alibay, I., Burusco, K. K., Bruce, N. J., and Bryce, R. A. (2018). Identification of Rare Lewis Oligosaccharide Conformers in Aqueous Solution Using Enhanced Sampling Molecular Dynamics. J. Phys. Chem. B 122, 2462–2474. doi:10.1021/acs.jpcb.7b09841

PubMed Abstract | CrossRef Full Text | Google Scholar

Alm, R. A., Bina, J., Andrews, B. M., Doig, P., Hancock, R. E. W., and Trust, T. J. (2000). Comparative Genomics of Helicobacter pylori : Analysis of the Outer Membrane Protein Families. Infect. Immun. 68, 4155–4168. doi:10.1128/IAI.68.7.4155-4168.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Aspholm-Hurtig, M., Dailide, G., Lahmann, M., Kalia, A., Ilver, D., Roche, N., et al. (2004). Functional Adaptation of BabA, the H. pylori ABO Blood Group Antigen Binding Adhesin. Science 305, 519–522. doi:10.1126/science.1098801

PubMed Abstract | CrossRef Full Text | Google Scholar

Berendsen, H. J. C., Postma, J. P. M., van Gunsteren, W. F., DiNola, A., and Haak, J. R. (1984). Molecular Dynamics with Coupling to an External bath. J. Chem. Phys. 81, 3684–3690. doi:10.1063/1.448118

CrossRef Full Text | Google Scholar

Borén, T., Falk, P., Roth, K. A., Larson, G., and Normark, S. (1993). Attachment of Helicobacter pylori to Human Gastric Epithelium Mediated by Blood Group Antigens. Science 262, 1892–1895. doi:10.1126/science.8018146

PubMed Abstract | CrossRef Full Text | Google Scholar

Bugaytsova, J. A., Björnham, O., Chernov, Y. A., Gideonsson, P., Henriksson, S., Mendez, M., et al. (2017). Helicobacter pylori Adapts to Chronic Infection and Gastric Disease via pH-Responsive BabA-Mediated Adherence. Cell Host & Microbe 21, 376–389. doi:10.1016/j.chom.2017.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Case, D. A., Ben-Shalom, I. Y., Brozell, S. R., Cerutti, D. S., Cheatham, T. E., Cruzeiro, V. W. D., et al. (2018). AMBER 2018. San Francisco: University of California.

Google Scholar

Cave, D. R. (1996). Transmission and Epidemiology of Helicobacter pylori. Am. J. Med. 100, 12S–18S. doi:10.1016/S0002-9343(96)80224-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Cremer, D., and Pople, J. A. (1975). General Definition of Ring Puckering Coordinates. J. Am. Chem. Soc. 97, 1354–1358. doi:10.1021/ja00839a011

CrossRef Full Text | Google Scholar

Darden, T., York, D., and Pedersen, L. (1993). Particle Mesh Ewald: AnN⋅Log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 98, 10089–10092. doi:10.1063/1.464397

CrossRef Full Text | Google Scholar

Donnini, S., Tegeler, F., Groenhof, G., and Grubmüller, H. (2011). Constant pH Molecular Dynamics in Explicit Solvent with λ-Dynamics. J. Chem. Theor. Comput. 7, 1962–1978. doi:10.1021/ct200061r

CrossRef Full Text | Google Scholar

Doohan, D., Rezkitha, Y. A. A., Waskito, L. A., Yamaoka, Y., and Miftahussurur, M. (2021). Helicobacter pylori BabA-SabA Key Roles in the Adherence Phase: The Synergic Mechanism for Successful Colonization and Disease Development. Toxins 13, 485. doi:10.3390/toxins13070485

PubMed Abstract | CrossRef Full Text | Google Scholar

Evans, D. J., and Evans, D. G. (2000). Helicobacter pylori Adhesins: Review and Perspectives. Helicobacter 5, 183–195. doi:10.1046/j.1523-5378.2000.00029.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hage, N., Howard, T., Phillips, C., Brassington, C., Overman, R., Debreczeni, J., et al. (2015). Structural Basis of Lewis B Antigen Binding by the Helicobacter pylori Adhesin BabA. Sci. Adv. 1, e1500315. doi:10.1126/sciadv.1500315

PubMed Abstract | CrossRef Full Text | Google Scholar

Hatakeyama, M. (2017). A Sour Relationship between BabA and Lewis B. Cell Host & Microbe 21, 318–320. doi:10.1016/j.chom.2017.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Ho, S. B., Roberton, A. M., Shekels, L. L., Lyftogt, C. T., Niehans, G. A., and Toribara, N. W. (1995). Expression Cloning of Gastric Mucin Complementary DNA and Localization of Mucin Gene Expression. Gastroenterology 109, 735–747. doi:10.1016/0016-5085(95)90380-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Ilver, D., Arnqvist, A., Ögren, J., Frick, I.-M., Kersulyte, D., Incecik, E. T., et al. (1998). Helicobacter pylori Adhesin Binding Fucosylated Histo-Blood Group Antigens Revealed by Retagging. Science 279, 373–377. doi:10.1126/science.279.5349.373

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishijima, N., Suzuki, M., Ashida, H., Ichikawa, Y., Kanegae, Y., Saito, I., et al. (2011). BabA-mediated Adherence Is a Potentiator of the Helicobacter pylori Type IV Secretion System Activity. J. Biol. Chem. 286, 25256–25264. doi:10.1074/jbc.M111.233601

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeffrey, G. A. (1997). An Introduction to Hydrogen Bonding. New York Oxford: Oxford University Press.

Google Scholar

Jonniya, N. A., Zhang, J., and Kar, P. (2021). Molecular Mechanism of Inhibiting WNK Binding to OSR1 by Targeting the Allosteric Pocket of the OSR1-CCT Domain with Potential Antihypertensive Inhibitors: An In Silico Study. J. Phys. Chem. B 125, 9115–9129. doi:10.1021/acs.jpcb.1c04672

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, J., and Blaser, M. J. (2006). Bacterial Populations as Perfect Gases: Genomic Integrity and Diversification Tensions in Helicobacter pylori. Nat. Rev. Microbiol. 4, 826–836. doi:10.1038/nrmicro1528

PubMed Abstract | CrossRef Full Text | Google Scholar

Kar, P., Lipowsky, R., and Knecht, V. (2013). Importance of Polar Solvation and Configurational Entropy for Design of Antiretroviral Drugs Targeting HIV-1 Protease. J. Phys. Chem. B 117, 5793–5805. doi:10.1021/jp3085292

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirschner, K. N., Yongye, A. B., Tschampel, S. M., González-Outeiriño, J., Daniels, C. R., Foley, B. L., et al. (2008). GLYCAM06: A Generalizable Biomolecular Force Field. Carbohydrates: GLYCAM06. J. Comput. Chem. 29, 622–655. doi:10.1002/jcc.20820

PubMed Abstract | CrossRef Full Text | Google Scholar

Kollman, P. A., Massova, I., Reyes, C., Kuhn, B., Huo, S., Chong, L., et al. (2000). Calculating Structures and Free Energies of Complex Molecules: Combining Molecular Mechanics and Continuum Models. Acc. Chem. Res. 33, 889–897. doi:10.1021/ar000033j

PubMed Abstract | CrossRef Full Text | Google Scholar

Kräutler, V., Van Gunsteren, W. F., and Hünenberger, P. H. (2001). A Fast SHAKE Algorithm to Solve Distance Constraint Equations for Small Molecules in Molecular Dynamics Simulations. J. Comput. Chem. 22, 501–508. doi:10.1002/1096-987x(20010415)22:5<501:aid-jcc1021>3.0.co;2-v

CrossRef Full Text | Google Scholar

Krivov, S. V., and Karplus, M. (2004). Hidden Complexity of Free Energy Surfaces for Peptide (Protein) Folding. Proc. Natl. Acad. Sci. U.S.A. 101, 14766–14770. doi:10.1073/pnas.0406234101

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindén, S., Nordman, H., Hedenbro, J., Hurtig, M., Borén, T., and Carlstedt, I. (2002). Strain- and Blood Group-dependent Binding of Helicobacter pylori to Human Gastric MUC5AC Glycoforms. Gastroenterology 123, 1923–1930. doi:10.1053/gast.2002.37076

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindén, S., Borén, T., Dubois, A., and Carlstedt, I. (2004). Rhesus Monkey Gastric Mucins: Oligomeric Structure, Glycoforms and Helicobacter pylori Binding. Biochem. J. 379, 765–775. doi:10.1042/bj20031557

PubMed Abstract | CrossRef Full Text | Google Scholar

Maier, J. A., Martinez, C., Kasavajhala, K., Wickstrom, L., Hauser, K. E., and Simmerling, C. (2015). ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theor. Comput. 11, 3696–3713. doi:10.1021/acs.jctc.5b00255

PubMed Abstract | CrossRef Full Text | Google Scholar

Mishra, S. K., and Koča, J. (2018). Assessing the Performance of MM/PBSA, MM/GBSA, and QM-MM/GBSA Approaches on Protein/Carbohydrate Complexes: Effect of Implicit Solvent Models, QM Methods, and Entropic Contributions. J. Phys. Chem. B 122, 8113–8121. doi:10.1021/acs.jpcb.8b03655

PubMed Abstract | CrossRef Full Text | Google Scholar

Mongan, J., Case, D. A., and McCammon, J. A. (2004). Constant pH Molecular Dynamics in Generalized Born Implicit Solvent. J. Comput. Chem. 25, 2038–2048. doi:10.1002/jcc.20139

PubMed Abstract | CrossRef Full Text | Google Scholar

Moonens, K., Gideonsson, P., Subedi, S., Bugaytsova, J., Romaõ, E., Mendez, M., et al. (2016). Structural Insights into Polymorphic ABO Glycan Binding by Helicobacter pylori. Cell Host Microbe 19, 55–66. doi:10.1016/j.chom.2015.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Nordman, H., Davies, J. R., Lindell, G., and Carlstedt, I. (1995). Human Gastric Mucins - a Major Population Identified as MUC5. Biochem. Soc. Trans. 23, 533S. doi:10.1042/bst023533s

PubMed Abstract | CrossRef Full Text | Google Scholar

Olsson, M. H. M., Søndergaard, C. R., Rostkowski, M., and Jensen, J. H. (2011). PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. J. Chem. Theor. Comput. 7, 525–537. doi:10.1021/ct100578z

PubMed Abstract | CrossRef Full Text | Google Scholar

Onufriev, A., Bashford, D., and Case, D. A. (2004). Exploring Protein Native States and Large-Scale Conformational Changes with a Modified Generalized Born Model. Proteins 55, 383–394. doi:10.1002/prot.20033

PubMed Abstract | CrossRef Full Text | Google Scholar

Pang, S. S., Nguyen, S. T. S., Perry, A. J., Day, C. J., Panjikar, S., Tiralongo, J., et al. (2014). The Three-Dimensional Structure of the Extracellular Adhesion Domain of the Sialic Acid-Binding Adhesin SabA from Helicobacter pylori. J. Biol. Chem. 289, 6332–6340. doi:10.1074/jbc.M113.513135

PubMed Abstract | CrossRef Full Text | Google Scholar

Polk, D. B., and Peek, R. M. (2010). Helicobacter pylori: Gastric Cancer and beyond. Nat. Rev. Cancer 10, 403–414. doi:10.1038/nrc2857

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, D. J., and Brooks, C. L. (2004). A Modified TIP3P Water Potential for Simulation With Ewald Summation. J. Chem. Phys. 121, 10096–10103. doi:10.1063/1.1808117

PubMed Abstract | CrossRef Full Text | Google Scholar

Roe, D. R., and Cheatham, T. E. (2013). PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J. Chem. Theor. Comput. 9, 3084–3095. doi:10.1021/ct400341p

PubMed Abstract | CrossRef Full Text | Google Scholar

Roy, R., Ghosh, B., and Kar, P. (2020a). Investigating Conformational Dynamics of Lewis Y Oligosaccharides and Elucidating Blood Group Dependency of Cholera Using Molecular Dynamics. ACS Omega 5, 3932–3942. doi:10.1021/acsomega.9b03398

PubMed Abstract | CrossRef Full Text | Google Scholar

Roy, R., Mishra, A., Poddar, S., Nayak, D., and Kar, P. (2020b). Investigating the mechanism of recognition and structural dynamics of nucleoprotein-RNA complex from Peste des petits ruminants virus via Gaussian accelerated molecular dynamics simulations. J. Biomol. Struct. Dyn. 40, 2302–2315. doi:10.1080/07391102.2020.1838327

PubMed Abstract | CrossRef Full Text | Google Scholar

Roy, R., Jonniya, N. A., Poddar, S., Sk, M. F., and Kar, P. (2021). Unraveling the Molecular Mechanism of Recognition of Human Interferon-Stimulated Gene Product 15 by Coronavirus Papain-Like Proteases: A Multiscale Simulation Study. J. Chem. Inf. Model. 61, 6038–6052. doi:10.1021/acs.jcim.1c00918

PubMed Abstract | CrossRef Full Text | Google Scholar

Roy, R., Poddar, S., Sk, M. F., and Kar, P. (2022). Conformational Preferences of Triantennary and Tetraantennary Hybrid N-Glycans in Aqueous Solution: Insights from 20 μs Long Atomistic Molecular Dynamic Simulations. J. Biomol. Struct. Dyn., 1–16. doi:10.1080/07391102.2022.2047109

PubMed Abstract | CrossRef Full Text | Google Scholar

Schreiber, S., Konradt, M., Groll, C., Scheid, P., Hanauer, G., Werling, H.-O., et al. (2004). The Spatial Orientation of Helicobacter pylori in the Gastric Mucus. Proc. Natl. Acad. Sci. U.S.A. 101, 5024–5029. doi:10.1073/pnas.0308386101

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheu, B.-S. (2003). Host Gastric Lewis Expression Determines the Bacterial Density of Helicobacter pylori in babA2 Genopositive Infection. Gut 52, 927–932. doi:10.1136/gut.52.7.927

PubMed Abstract | CrossRef Full Text | Google Scholar

Siemers, M., and Bondar, A.-N. (2021). Interactive Interface for Graph-Based Analyses of Dynamic H-Bond Networks: Application to Spike Protein S. J. Chem. Inf. Model. 61, 2998–3014. doi:10.1021/acs.jcim.1c00306

CrossRef Full Text | Google Scholar

Sood, A., Gerlits, O. O., Ji, Y., Bovin, N. V., Coates, L., and Woods, R. J. (2018). Defining the Specificity of Carbohydrate-Protein Interactions by Quantifying Functional Group Contributions. J. Chem. Inf. Model. 58, 1889–1901. doi:10.1021/acs.jcim.8b00120

PubMed Abstract | CrossRef Full Text | Google Scholar

Strous, G. J., and Dekker, J. (1992). Mucin-Type Glycoproteins. Crit. Rev. Biochem. Mol. Biol. 27, 57–92. doi:10.3109/10409239209082559

PubMed Abstract | CrossRef Full Text | Google Scholar

Suerbaum, S., and Josenhans, C. (2007). Helicobacter pylori Evolution and Phenotypic Diversification in a Changing Host. Nat. Rev. Microbiol. 5, 441–452. doi:10.1038/nrmicro1658

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamaoka, Y. (2008). Roles of Helicobacter pylori BabA in Gastroduodenal Pathogenesis. WJG 14, 4265. doi:10.3748/wjg.14.4265

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: conformational dynamics, multiple replica molecular dynamics simulations, BabA, free energy landscape, principal component analysis, phi/psi maps, loop dynamics

Citation: Roy R, Jonniya NA, Sk MF and Kar P (2022) Comparative Structural Dynamics of Isoforms of Helicobacter pylori Adhesin BabA Bound to Lewis b Hexasaccharide via Multiple Replica Molecular Dynamics Simulations. Front. Mol. Biosci. 9:852895. doi: 10.3389/fmolb.2022.852895

Received: 11 January 2022; Accepted: 30 March 2022;
Published: 02 May 2022.

Edited by:

Neelanjana Sengupta, Indian Institute of Science Education and Research Kolkata, India

Reviewed by:

Suman Chakrabarty, S.N. Bose National Centre for Basic Sciences, India
Manuel Hitzenberger, Technical University of Munich, Germany

Copyright © 2022 Roy, Jonniya, Sk and Kar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Parimal Kar, parimal@iiti.ac.in

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.