Conformational dynamics of L-lysine, L-arginine, L-ornithine binding protein reveals ligand-dependent plasticity

The role of C-terminal helix in the conformational transition of an arginine binding protein

•

Probe the role of C-ter. helix (CTH) in conformational transition of TmArgBP.

•

Presence of CTH almost doubles the barrier to access the closed-state.

•

In the absence of CTH, the protein can fluctuate between the two conformations.

•

CTH not only constraints the open-state conformation but also guides in accessing it.

The thermotoga maritima arginine binding protein (TmArgBP) is a periplasmic binding protein that has a short helix at the C-terminal end (CTH), which is swapped between the two chains. We apply a coarse-grained structure-based model (SBM) and all-atom MD simulation on this protein to understand the mechanism and the role of CTH in the conformational transition. When the results of SBM simulations of TmArgBP in the presence and absence of CTH are compared, we find that CTH is strategically located at the back of the binding pocket restraining the open-state conformation thereby disengaging access to the closed-state. We also ran all-atom MD simulations of open-state TmArgBP with and without CTH and discovered that in the absence of CTH the protein could reach the closed-state within 250 ns, while in its presence, the protein remained predominantly in its open-state conformation. In the simulation started from unliganded closed-state conformation without CTH, the protein exhibited multiple transitions between the two states, suggesting CTH as an essential structural element to stabilize the open-state conformation. In another simulation that began with an unliganded closed-state conformation with CTH, the protein was able to access the open-state. In this simulation the CTH was observed to reorient itself to interact with the protein emphasizing its role in assisting the conformational change. Based on our findings, we believe that CTH not only acts as a structural element that constraints the protein in its open-state but it may also guide the protein back to its open-state conformation upon ligand unbinding.

Evolution of cyclohexadienyl dehydratase from an ancestral solute-binding protein

The emergence of enzymes through the neofunctionalization of noncatalytic proteins is ultimately responsible for the extraordinary range of biological catalysts observed in nature. Although the evolution of some enzymes from binding proteins can be inferred by homology, we have a limited understanding of the nature of the biochemical and biophysical adaptations along these evolutionary trajectories and the sequence in which they occurred. Here we reconstructed and characterized evolutionary intermediate states linking an ancestral solute-binding protein to the extant enzyme cyclohexadienyl dehydratase. We show how the intrinsic reactivity of a desolvated general acid was harnessed by a series of mutations radiating from the active site, which optimized enzyme-substrate complementarity and transition-state stabilization and minimized sampling of noncatalytic conformations. Our work reveals the molecular evolutionary processes that underlie the emergence of enzymes de novo, which are notably mirrored by recent examples of computational enzyme design and directed evolution.

Role of cis-trans proline isomerization in the function of pathogenic enterobacterial Periplasmic Binding Proteins

Periplasmic Binding Proteins (PBPs) trap nutrients for their internalization into bacteria by ABC transporters. Ligand binding triggers PBP closure by bringing its two domains together like a Venus flytrap. The atomic determinants that control PBP opening and closure for nutrient capture and release are not known, although it is proposed that opening and ligand release occur while in contact with the ABC transporter for concurrent substrate translocation. In this paper we evaluated the effect of the isomerization of a conserved proline, located near the binding site, on the propensity of PBPs to open and close. ArgT/LAO from Salmonella typhimurium and HisJ from Escherichia coli were studied through molecular mechanics at two different temperatures: 300 and 323 K. Eight microseconds were simulated per protein to analyze protein opening and closure in the absence of the ABC transporter. We show that when the studied proline is in trans, closed empty LAO and HisJ can open. In contrast, with the proline in cis, opening transitions were much less frequent and characterized by smaller changes. The proline in trans also renders the open trap prone to close over a ligand. Our data suggest that the isomerization of this conserved proline modulates the PBP mechanism: the proline in trans allows the exploration of conformational space to produce trap opening and closure, while in cis it restricts PBP movement and could limit ligand release until in productive contact with the ABC transporter. This is the first time that a proline isomerization has been related to the control of a large conformational change like the PBP flytrap mechanism.

Designing, construction and characterization of genetically encoded FRET-based nanosensor for real time monitoring of lysine flux in living cells

BACKGROUND

Engineering microorganisms in order to improve the metabolite flux needs a detailed knowledge of the concentrations and flux rates of metabolites and metabolic intermediates in vivo. Fluorescence resonance energy transfer (FRET) based genetically encoded nanosensors represent a promising tool for measuring the metabolite levels and corresponding rate changes in live cells. Here, we report the development of a series of FRET based genetically encoded nanosensor for real time measurement of lysine at cellular level, as the improvement of microbial strains for the production of L-lysine is of major interest in industrial biotechnology.

RESULTS

The lysine binding periplasmic protein (LAO) from Salmonella enterica serovar typhimurium LT2 strain was used as the reporter element for the sensor. The LAO was sandwiched between GFP variants i.e. cyan fluorescent protein (CFP) and yellow fluorescent protein (YFP). Affinity, pH stability, specificity and metal ions effects was scrutinized for the in vitro characterization of this nanosensor, named as FLIPK. The FLIPK is specific to lysine and found to be stable with the pH within the physiological range. The calculated affinity (K d ) of FLIPK was 97 µM. For physiological applications, mutants with different binding affinities were also generated and investigated in vitro. The developed nanosensor efficiently monitored the intracellular level of lysine in bacterial as well as yeast cell.

CONCLUSION

The developed novel lysine fluorescence resonance energy transfer sensors can be used for in vivo monitoring of lysine levels in prokaryotes as well as eukaryotes. The potential of these sensors is that they can be used as reporter tools in the development of metabolically engineered microbial strains or for real-time monitoring of intracellular lysine during fermentation.

Ancestral Protein Reconstruction Yields Insights into Adaptive Evolution of Binding Specificity in Solute-Binding Proteins

The promiscuous functions of proteins are an important reservoir of functional novelty in protein evolution, but the molecular basis for binding promiscuity remains elusive. We used ancestral protein reconstruction to experimentally characterize evolutionary intermediates in the functional expansion of the polar amino acid-binding protein family, which has evolved to bind a variety of amino acids with high affinity and specificity. High-resolution crystal structures of an ancestral arginine-binding protein in complex with l-arginine and l-glutamine show that the promiscuous binding of l-glutamine is enabled by multi-scale conformational plasticity, water-mediated interactions, and selection of an alternative conformational substate productive for l-glutamine binding. Evolution of specialized glutamine-binding proteins from this ancestral protein was achieved by displacement of water molecules from the protein-ligand interface, reducing the entropic penalty associated with the promiscuous interaction. These results provide a structural and thermodynamic basis for the co-option of a promiscuous interaction in the evolution of binding specificity.

A consistent force field parameter set for zwitterionic amino acid residues

Isolated amino acids play an important role in biochemistry and are therefore an interesting object of study. Atomistic molecular dynamics (MD) simulations can provide a high-resolution picture of the dynamic features of these species, especially in their biological environment. Unfortunately, most standard force field packages lack libraries for isolated amino acids in their zwitterionic form. Although several studies have used ad-hoc parameterizations for single amino acids, a consistent force-field parameter set for these molecules is still missing. Here, we present such a parameter library derived from the widely used parm99SB set from the AMBER program package. The parameter derivation for all 20 proteinogenic amino acids transparently followed established procedures with histidine treated in three different protonation states. All amino acids were subjected to MD simulations in four different forms for comparison: zwitterionic, N-teminally capped with acetyl, C-terminally capped with N-methyl, and capped at both termini. Simulation results show similarities between the different forms. Five zwitterionic amino acids-arginine, glutamate, glycine, phenylalanine, leucine-were simulated in a protein environment. Proteins and ligands generally retained their initial structure. The new parameter set will thus facilitate future atomistic simulations of these species.

Quantitatively characterizing the ligand binding mechanisms of choline binding protein using Markov state model analysis

Protein-ligand recognition plays key roles in many biological processes. One of the most fascinating questions about protein-ligand recognition is to understand its underlying mechanism, which often results from a combination of induced fit and conformational selection. In this study, we have developed a three-pronged approach of Markov State Models, Molecular Dynamics simulations, and flux analysis to determine the contribution of each model. Using this approach, we have quantified the recognition mechanism of the choline binding protein (ChoX) to be ∼90% conformational selection dominant under experimental conditions. This is achieved by recovering all the necessary parameters for the flux analysis in combination with available experimental data. Our results also suggest that ChoX has several metastable conformational states, of which an apo-closed state is dominant, consistent with previous experimental findings. Our methodology holds great potential to be widely applied to understand recognition mechanisms underlining many fundamental biological processes.

Role of the two structural domains from the periplasmic Escherichia coli histidine-binding protein HisJ

Escherichia coli HisJ is a type II periplasmic binding protein that functions to reversibly capture histidine and transfer it to its cognate inner membrane ABC permease. Here, we used NMR spectroscopy to determine the structure of apo-HisJ (26.5 kDa) in solution. HisJ is a bilobal protein in which domain 1 (D1) is made up of two noncontiguous subdomains, and domain 2 (D2) is expressed as the inner domain. To better understand the roles of D1 and D2, we have isolated and characterized each domain separately. Structurally, D1 closely resembles its homologous domain in apo- and holo-HisJ, whereas D2 is more similar to the holo-form. NMR relaxation experiments reveal that HisJ becomes more ordered upon ligand binding, whereas isolated D2 experiences a significant reduction in slower (millisecond to microsecond) motions compared with the homologous domain in apo-HisJ. NMR titrations reveal that D1 is able to bind histidine in a similar manner as full-length HisJ, albeit with lower affinity. Unexpectedly, isolated D1 and D2 do not interact with each other in the presence or absence of histidine, which indicates the importance of intact interdomain-connecting elements (i.e. hinge regions) for HisJ functioning. Our results shed light on the binding mechanism of type II periplasmic binding proteins where ligand is initially bound by D1, and D2 plays a supporting role in this dynamic process.

Behavior of solvent-exposed hydrophobic groove in the anti-apoptotic Bcl-XL protein: clues for its ability to bind diverse BH3 ligands from MD simulations

Bcl-XL is a member of Bcl-2 family of proteins involved in the regulation of intrinsic pathway of apoptosis. Its overexpression in many human cancers makes it an important target for anti-cancer drugs. Bcl-XL interacts with the BH3 domain of several pro-apoptotic Bcl-2 partners. This helical bundle protein has a pronounced hydrophobic groove which acts as a binding region for the BH3 domains. Eight independent molecular dynamics simulations of the apo/holo forms of Bcl-XL were carried out to investigate the behavior of solvent-exposed hydrophobic groove. The simulations used either a twin-range cut-off or particle mesh Ewald (PME) scheme to treat long-range interactions. Destabilization of the BH3 domain-containing helix H2 was observed in all four twin-range cut-off simulations. Most of the other major helices remained stable. The unwinding of H2 can be related to the ability of Bcl-XL to bind diverse BH3 ligands. The loss of helical character can also be linked to the formation of homo- or hetero-dimers in Bcl-2 proteins. Several experimental studies have suggested that exposure of BH3 domain is a crucial event before they form dimers. Thus unwinding of H2 seems to be functionally very important. The four PME simulations, however, revealed a stable helix H2. It is possible that the H2 unfolding might occur in PME simulations at longer time scales. Hydrophobic residues in the hydrophobic groove are involved in stable interactions among themselves. The solvent accessible surface areas of bulky hydrophobic residues in the groove are significantly buried by the loop LB connecting the helix H2 and subsequent helix. These observations help to understand how the hydrophobic patch in Bcl-XL remains stable in the solvent-exposed state. We suggest that both the destabilization of helix H2 and the conformational heterogeneity of loop LB are important factors for binding of diverse ligands in the hydrophobic groove of Bcl-XL.

A role for both conformational selection and induced fit in ligand binding by the LAO protein

Molecular recognition is determined by the structure and dynamics of both a protein and its ligand, but it is difficult to directly assess the role of each of these players. In this study, we use Markov State Models (MSMs) built from atomistic simulations to elucidate the mechanism by which the Lysine-, Arginine-, Ornithine-binding (LAO) protein binds to its ligand. We show that our model can predict the bound state, binding free energy, and association rate with reasonable accuracy and then use the model to dissect the binding mechanism. In the past, this binding event has often been assumed to occur via an induced fit mechanism because the protein's binding site is completely closed in the bound state, making it impossible for the ligand to enter the binding site after the protein has adopted the closed conformation. More complex mechanisms have also been hypothesized, but these have remained controversial. Here, we are able to directly observe roles for both the conformational selection and induced fit mechanisms in LAO binding. First, the LAO protein tends to form a partially closed encounter complex via conformational selection (that is, the apo protein can sample this state), though the induced fit mechanism can also play a role here. Then, interactions with the ligand can induce a transition to the bound state. Based on these results, we propose that MSMs built from atomistic simulations may be a powerful way of dissecting ligand-binding mechanisms and may eventually facilitate a deeper understanding of allostery as well as the prediction of new protein-ligand interactions, an important step in drug discovery.