The energetic cost of domain reorientation in maltose-binding protein as studied by NMR and fluorescence spectroscopy

Conformational Changes in Surface-Immobilized Proteins Measured Using Combined Atomic Force and Fluorescence Microscopy

Biological organisms rely on proteins to perform the majority of their functions. Most protein functions are based on their physical motions (conformational changes), which can be described as transitions between different conformational states in a multidimensional free-energy landscape. A comprehensive understanding of this free-energy landscape is therefore of paramount importance for understanding the biological functions of proteins. Protein dynamics includes both equilibrium and nonequilibrium motions, which typically exhibit a wide range of characteristic length and time scales. The relative probabilities of various conformational states in the energy landscape, the energy barriers between them, their dependence on external parameters such as force and temperature, and their connection to the protein function remain largely unknown for most proteins. In this paper, we present a multimolecule approach in which the proteins are immobilized at well-defined locations on Au substrates using an atomic force microscope (AFM)-based patterning method called nanografting. This method enables precise control over the protein location and orientation on the substrate, as well as the creation of biologically active protein ensembles that self-assemble into well-defined nanoscale regions (protein patches) on the gold substrate. We performed AFM-force compression and fluorescence experiments on these protein patches and measured the fundamental dynamical parameters such as protein stiffness, elastic modulus, and transition energies between distinct conformational states. Our results provide new insights into the processes that govern protein dynamics and its connection to protein function.

Protein dynamics: The future is bright and complicated!

Biological life depends on motion, and this manifests itself in proteins that display motion over a formidable range of time scales spanning from femtoseconds vibrations of atoms at enzymatic transition states, all the way to slow domain motions occurring on micro to milliseconds. An outstanding challenge in contemporary biophysics and structural biology is a quantitative understanding of the linkages among protein structure, dynamics, and function. These linkages are becoming increasingly explorable due to conceptual and methodological advances. In this Perspective article, we will point toward future directions of the field of protein dynamics with an emphasis on enzymes. Research questions in the field are becoming increasingly complex such as the mechanistic understanding of high-order interaction networks in allosteric signal propagation through a protein matrix, or the connection between local and collective motions. In analogy to the solution to the "protein folding problem," we argue that the way forward to understanding these and other important questions lies in the successful integration of experiment and computation, while utilizing the present rapid expansion of sequence and structure space. Looking forward, the future is bright, and we are in a period where we are on the doorstep to, at least in part, comprehend the importance of dynamics for biological function.

Huanglongbing Pandemic: Current Challenges and Emerging Management Strategies

Huanglongbing (HLB, aka citrus greening), one of the most devastating diseases of citrus, has wreaked havoc on the global citrus industry in recent decades. The culprit behind such a gloomy scenario is the phloem-limited bacteria "Candidatus Liberibacter asiaticus" (CLas), which are transmitted via psyllid. To date, there are no effective long-termcommercialized control measures for HLB, making it increasingly difficult to prevent the disease spread. To combat HLB effectively, introduction of multipronged management strategies towards controlling CLas population within the phloem system is deemed necessary. This article presents a comprehensive review of up-to-date scientific information about HLB, including currently available management practices and unprecedented challenges associated with the disease control. Additionally, a triangular disease management approach has been introduced targeting pathogen, host, and vector. Pathogen-targeting approaches include (i) inhibition of important proteins of CLas, (ii) use of the most efficient antimicrobial or immunity-inducing compounds to suppress the growth of CLas, and (iii) use of tools to suppress or kill the CLas. Approaches for targeting the host include (i) improvement of the host immune system, (ii) effective use of transgenic variety to build the host's resistance against CLas, and (iii) induction of systemic acquired resistance. Strategies for targeting the vector include (i) chemical and biological control and (ii) eradication of HLB-affected trees. Finally, a hypothetical model for integrated disease management has been discussed to mitigate the HLB pandemic.

Concurrent Identification and Characterization of Protein Structure and Continuous Internal Dynamics with REDCRAFT

Internal dynamics of proteins can play a critical role in the biological function of some proteins. Several well documented instances have been reported such as MBP, DHFR, hTS, DGCR8, and NSP1 of the SARS-CoV family of viruses. Despite the importance of internal dynamics of proteins, there currently are very few approaches that allow for meaningful separation of internal dynamics from structural aspects using experimental data. Here we present a computational approach named REDCRAFT that allows for concurrent characterization of protein structure and dynamics. Here, we have subjected DHFR (PDB-ID 1RX2), a 159-residue protein, to a fictitious, mixed mode model of internal dynamics. In this simulation, DHFR was segmented into 7 regions where 4 of the fragments were fixed with respect to each other, two regions underwent rigid-body dynamics, and one region experienced uncorrelated and melting event. The two dynamical and rigid-body segments experienced an average orientational modification of 7° and 12° respectively. Observable RDC data for backbone C′-N, N-HN, and C′-HN were generated from 102 uniformly sampled frames that described the molecular trajectory. The structure calculation of DHFR with REDCRAFT by using traditional Ramachandran restraint produced a structure with 29 Å of structural difference measured over the backbone atoms (bb-rmsd) over the entire length of the protein and an average bb-rmsd of more than 4.7 Å over each of the dynamical fragments. The same exercise repeated with context-specific dihedral restraints generated by PDBMine produced a structure with bb-rmsd of 21 Å over the entire length of the protein but with bb-rmsd of less than 3 Å over each of the fragments. Finally, utilization of the Dynamic Profile generated by REDCRAFT allowed for the identification of different dynamical regions of the protein and the recovery of individual fragments with bb-rmsd of less than 1 Å. Following the recovery of the fragments, our assembly procedure of domains (larger segments consisting of multiple fragments with a common dynamical profile) correctly assembled the four fragments that are rigid with respect to each other, categorized the two domains that underwent rigid-body dynamics, and identified one dynamical region for which no conserved structure could be defined. In conclusion, our approach was successful in identifying the dynamical domains, recovery of structure where it is meaningful, and relative assembly of the domains when possible.

Stability of Ligand-induced Protein Conformation Influences Affinity in Maltose-binding Protein

Our understanding of what determines ligand affinity of proteins is poor, even with high-resolution structures available. Both the non-covalent ligand-protein interactions and the relative free energies of available conformations contribute to the affinity of a protein for a ligand. Distant, non-binding site residues can influence the ligand affinity by altering the free energy difference between a ligand-free and ligand-bound conformation. Our hypothesis is that when different ligands induce distinct ligand-bound conformations, it should be possible to tweak their affinities by changing the free energies of the available conformations. We tested this idea for the maltose-binding protein (MBP) from Escherichia coli. We used single-molecule Förster resonance energy transfer (smFRET) to distinguish several unique ligand-bound conformations of MBP. We engineered mutations, distant from the binding site, to affect the stabilities of different ligand-bound conformations. We show that ligand affinity can indeed be altered in a conformation-dependent manner. Our studies provide a framework for the tuning of ligand affinity, apart from modifying binding site residues.

Exploring Large Domain Motions in Proteins Using Atomistic Molecular Dynamics with Enhanced Conformational Sampling

Conformational transitions in multidomain proteins are essential for biological functions. The Apo conformations are typically open and flexible, while the Holo states form more compact conformations stabilized by protein-ligand interactions. Unfortunately, the atomically detailed mechanisms for such open-closed conformational changes are difficult to be accessed experimentally as well as computationally. To simulate the transitions using atomistic molecular dynamics (MD) simulations, efficient conformational sampling algorithms are required. In this work, we propose a new approach based on generalized replica-exchange with solute tempering (gREST) for exploring the open-closed conformational changes in multidomain proteins. Wherein, selected surface charged residues in a target protein are defined as the solute region in gREST simulation and the solute temperatures are different in replicas and exchanged between them to enhance the domain motions. This approach is called gREST selected surface charged residues (gREST_SSCR) and is applied to the Apo and Holo states of ribose binding protein (RBP) in solution. The conformational spaces sampled with gREST_SSCR are much wider than those with the conventional MD, sampling open-closed conformational changes while maintaining RBP domains’ stability. The free-energy landscapes of RBP in the Apo and Holo states are drawn along with twist and hinge angles of the two moving domains. The inter-domain salt-bridges that are not observed in the experimental structures are also important in the intermediate states during the conformational changes.

Structure dictates the mechanism of ligand recognition in the histidine and maltose binding proteins

Two mechanisms, induced fit (IF) and conformational selection (CS), have been proposed to explain ligand recognition coupled conformational changes. The histidine binding protein (HisJ) adopts the CS mechanism, in which a pre-equilibrium is established between the open and the closed states with the ligand binding to the closed state. Despite being structurally similar to HisJ, the maltose binding protein (MBP) adopts the IF mechanism, in which the ligand binds the open state and induces a transition to the closed state. To understand the molecular determinants of this difference, we performed molecular dynamics (MD) simulations of coarse-grained dual structure based models. We find that intra-protein contacts unique to the closed state are sufficient to promote the conformational transition in HisJ, indicating a CS-like mechanism. In contrast, additional ligand-mimicking contacts are required to “induce” the conformational transition in MBP suggesting an IF-like mechanism. In agreement with experiments, destabilizing modifications to two structural features, the spine helix (SH) and the balancing interface (BI), present in MBP but absent in HisJ, reduce the need for ligand-mimicking contacts indicating that SH and BI act as structural restraints that keep MBP in the open state. We introduce an SH like element into HisJ and observe that this can impede the conformational transition increasing the importance of ligand-mimicking contacts. Similarly, simultaneous mutations to BI and SH in MBP reduce the barrier to conformational transitions significantly and promote a CS-like mechanism. Together, our results show that structural restraints present in the protein structure can determine the mechanism of conformational transitions and even simple models that correctly capture such structural features can predict their positions. MD simulations of such models can thus be used, in conjunction with mutational experiments, to regulate protein ligand interactions, and modulate ligand binding affinities.

•

MBP operates by induced fit, HisJ by the conformational selection mechanism.

•

Dual structure based models (dSBMs) encode two structures of a protein.

•

MD simulations of dSBMs can identify the mechanism of conformational transitions.

•

Locks, absent in HisJ, hold MBP open with ligand contacts required for closing.

•

Binding mechanisms can be modified by altering such structural locks.

Ligand-bound glutamine binding protein assumes multiple metastable binding sites with different binding affinities

Protein dynamics plays key roles in ligand binding. However, the microscopic description of conformational dynamics-coupled ligand binding remains a challenge. In this study, we integrate molecular dynamics simulations, Markov state model (MSM) analysis and experimental methods to characterize the conformational dynamics of ligand-bound glutamine binding protein (GlnBP). We show that ligand-bound GlnBP has high conformational flexibility and additional metastable binding sites, presenting a more complex energy landscape than the scenario in the absence of ligand. The diverse conformations of GlnBP demonstrate different binding affinities and entail complex transition kinetics, implicating a concerted ligand binding mechanism. Single molecule fluorescence resonance energy transfer measurements and mutagenesis experiments are performed to validate our MSM-derived structure ensemble as well as the binding mechanism. Collectively, our study provides deeper insights into the protein dynamics-coupled ligand binding, revealing an intricate regulatory network underlying the apparent binding affinity.

Zhang, Wu, Feng et al. show that ligand-bound glutamine binding protein assumes multiple metastable binding sites, presenting a more dynamic energy landscape than its ligand-free form. This study provides insights into the ligand-binding mechanisms coupled with protein dynamics that underly the apparent binding affinity.

Current Blockades of Proteins inside Nanopores for Real-Time Metabolome Analysis

Biological nanopores are emerging as powerful and low-cost sensors for real-time analysis of biological samples. Proteins can be incorporated inside the nanopore, and ligand binding to the protein adaptor yields changes in nanopore conductance. In order to understand the origin of these conductance changes and develop sensors for detecting metabolites, we tested the signal originating from 13 different protein adaptors. We found that the quality of the protein signal depended on both the size and charge of the protein. The engineering of a dipole within the surface of the adaptor reduced the current noise by slowing the protein dynamics within the nanopore. Further, the charge of the ligand and the induced conformational changes of the adaptor defined the conductance changes upon metabolite binding, suggesting that the protein resides in an electrokinetic minimum within the nanopore, the position of which is altered by the ligand. These results represent an important step toward understanding the dynamics of the electrophoretic trapping of proteins inside nanopores and will allow developing next-generation sensors for metabolome analysis.

Electrical Sensing of Phosphonates by Functional Coupling of Phosphonate Binding Protein PhnD to Solid-State Nanopores

Combining the stability of solid-state nanopores with the unique sensing properties of biological components in a miniaturized electrical hybrid nanopore device is a challenging approach to advance the sensitivity and selectivity of small-molecule detection in healthcare and environment analytics. Here, we demonstrate a simple method to design an electrical hybrid nanosensor comprising a bacterial binding protein tethered to a solid-state nanopore allowing high-affinity detection of phosphonates. The diverse family of bacterial substrate-binding proteins (SBPs) binds specifically and efficiently to various substances and has been implicated as an ideal biorecognition element for analyte detection in the design of hybrid bionanosensors. Here, we demonstrate that the coupling of the purified phosphonate binding protein PhnD via primary amines to the reactive NHS groups of P(DMAA-co-NMAS) polymers inside a single track-etched nanopore in poly(ethylene terephthalate) (PET) foils results in ligand-specific and concentration-dependent changes in the nanopore current. Application of the phosphonate 2-aminoethylphosphonate (2AEP) or ethylphosphonate (EP) induces a large conformational rearrangement in PnhD around the hinge in a venus flytrap mechanism resulting in a concentration depended on increase of the single pore current with binding affinities of 27 and 373 nM, respectively. Thus, the specificity and stability of this simple hybrid sensor concept combine the advantages of both, the diversity of ligand-specific substrate-binding proteins and solid-state nanopores encouraging further options to produce robust devices amenable to medical or environmental high-throughput-based applications in nanotechnology.

Both Ligands and Macromolecular Crowders Preferentially Bind to Closed Conformations of Maltose Binding Protein

In cellular environments, proteins not only interact with their specific partners but also encounter a high concentration of bystander macromolecules, or crowders. Nonspecific interactions with macromolecular crowders modulate the activities of proteins, but our knowledge about the rules of nonspecific interactions is still very limited. In previous work, we presented experimental evidence that macromolecular crowders acted competitively in inhibiting the binding of maltose binding protein (MBP) with its ligand maltose. Competition between a ligand and an inhibitor may result from binding to either the same site or different conformations of the protein. Maltose binds to the cleft between two lobes of MBP, and in a series of mutants, the affinities increased with an increase in the extent of lobe closure. Here we investigated whether macromolecular crowders also have a conformational or site preference when binding to MBP. The affinities of a polymer crowder, Ficoll70, measured by monitoring tryptophan fluorescence were 3-6-fold higher for closure mutants than for wild-type MBP. Competition between the ligand and crowder, as indicated by fitting of titration data and directly by nuclear magnetic resonance spectroscopy, and their similar preferences for closed MBP conformations further suggest the scenario in which the crowder, like maltose, preferentially binds to the interlobe cleft of MBP. Similar observations were made for bovine serum albumin as a protein crowder. Conformational and site preferences in MBP-crowder binding allude to the paradigm that nonspecific interactions can possess hallmarks of molecular recognition, which may be essential for intracellular organizations including colocalization of proteins and liquid-liquid phase separation.

Differential Substrate Recognition by Maltose Binding Proteins Influenced by Structure and Dynamics

The genome of the hyperthermophile Thermotoga maritima contains three isoforms of maltose binding protein (MBP) that are high-affinity receptors for di-, tri-, and tetrasaccharides. Two of these proteins (tmMBP1 and tmMBP2) share significant sequence identity, approximately 90%, while the third (tmMBP3) shares less than 40% identity. MBP from Escherichia coli (ecMBP) shares 35% sequence identity with the tmMBPs. This subset of MBP isoforms offers an interesting opportunity to investigate the mechanisms underlying the evolution of substrate specificity and affinity profiles in a genome where redundant MBP genes are present. In this study, the X-ray crystal structures of tmMBP1, tmMBP2, and tmMBP3 are reported in the absence and presence of oligosaccharides. tmMBP1 and tmMBP2 have binding pockets that are larger than that of tmMBP3, enabling them to bind to larger substrates, while tmMBP1 and tmMBP2 also undergo substrate-induced hinge bending motions (∼52°) that are larger than that of tmMBP3 (∼35°). Small-angle X-ray scattering was used to compare protein behavior in solution, and computer simulations provided insights into dynamics of these proteins. Comparing quantitative protein-substrate interactions and dynamical properties of tmMBPs with those of the promiscuous ecMBP and disaccharide selective Thermococcus litoralis MBP provides insights into the features that enable selective binding. Collectively, the results provide insights into how the structure and dynamics of tmMBP homologues enable them to differentiate between a myriad of chemical entities while maintaining their common fold.

Engineered synthetic antibodies as probes to quantify the energetic contributions of ligand binding to conformational changes in proteins

Conformational changes in proteins due to ligand binding are ubiquitous in biological processes and are integral to many biological systems. However, it is often challenging to link ligand-induced conformational changes to a resulting biological function because it is difficult to distinguish between the energetic components associated with ligand binding and those due to structural rearrangements. Here, we used a unique approach exploiting conformation-specific and regio-specific synthetic antibodies (sABs) to probe the energetic contributions of ligand binding to conformation changes. Using maltose-binding protein (MBP) as a model system, customized phage-display selections were performed to generate sABs that stabilize MBP in different conformational states, modulating ligand-binding affinity in competitive, allosteric, or peristeric manners. We determined that the binding of a closed conformation-specific sAB (sAB-11M) to MBP in the absence of maltose is entropically driven, providing new insight into designing antibody-stabilized protein interactions. Crystal structures of sABs bound to MBP, together with biophysical data, delineate the basis of free energy differences between different conformational states and confirm the use of the sABs as energy probes for dissecting enthalpic and entropic contributions to conformational transitions. Our work provides a foundation for investigating the energetic contributions of distinct conformational dynamics to specific biological outputs. We anticipate that our approach also may be valuable for analyzing the energy landscapes of regulatory proteins controlling biological responses to environmental changes.

Homodimerization enhances both sensitivity and dynamic range of the ligand-binding domain of type 1 metabotropic glutamate receptor

Cooperativity in ligand binding is a key emergent property of protein oligomers. Positive cooperativity (higher affinity for subsequent binding events than for initial binding) is frequent. However, the symmetrically homodimeric ligand-binding domain (LBD) of metabotropic glutamate receptor type 1 exhibits negative cooperativity. To investigate its origin and functional significance, we measured the response to glutamate in vitro of wild-type and C140S LBD as a function of the extent of dimerization. Our results indicate that homodimerization enhances the affinity of the first, but not the second, binding site, relative to the monomer, giving the dimeric receptor both greater sensitivity and a broader dynamic range.

Unraveling the Conformational Landscape of Ligand Binding to Glucose/Galactose-Binding Protein by Paramagnetic NMR and MD Simulations

Protein dynamics related to function can nowadays be structurally well characterized (i.e., instances obtained by high resolution structures), but they are still ill-defined energetically, and the energy landscapes are only accessible computationally. This is the case for glucose-galactose binding protein (GGBP), where the crystal structures of the apo and holo states provide structural information for the domain rearrangement upon ligand binding, while the time scale and the energetic determinants for such concerted dynamics have been so far elusive. Here, we use GGBP as a paradigm to define a functional conformational landscape, both structurally and energetically, by using an innovative combination of paramagnetic NMR experiments and MD simulations. Anisotropic NMR parameters induced by self-alignment of paramagnetic metal ions was used to characterize the ensemble of conformations adopted by the protein in solution while the rate of interconversion between conformations was elucidated by long molecular dynamics simulation on two states of GGBP, the closed-liganded (holo_cl) and open-unloaded (apo_op) states. Our results demonstrate that, in its apo state, the protein coexists between open-like (68%) and closed-like (32%) conformations, with an exchange rate around 25 ns. Despite such conformational heterogeneity, the presence of the ligand is the ultimate driving force to unbalance the equilibrium toward the holo_cl form, in a mechanism largely governed by a conformational selection mechanism.

Protein Conformational Changes Are Detected and Resolved Site Specifically by Second-Harmonic Generation

We present here a straightforward, broadly applicable technique for real-time detection and measurement of protein conformational changes in solution. This method is based on tethering proteins labeled with a second-harmonic generation (SHG) active dye to supported lipid bilayers. We demonstrate our method by measuring the conformational changes that occur upon ligand binding with three well-characterized proteins labeled at lysine residues: calmodulin (CaM), maltose-binding protein (MBP), and dihydrofolate reductase (DHFR). We also create a single-site cysteine mutant of DHFR engineered within the Met20 catalytic loop region and study the protein’s structural motion at this site. Using published x-ray crystal structures, we show that the changes in the SHG signals upon ligand binding are the result of structural motions that occur at the labeled sites between the apo and ligand-bound forms of the proteins, which are easily distinguished from each other. In addition, we demonstrate that different magnitudes of the SHG signal changes are due to different and specific ligand-induced conformational changes. Taken together, these data illustrate the potential of the SHG approach for detecting and measuring protein conformational changes for a wide range of biological applications.

Mapping the Free Energy Landscape of PKA Inhibition and Activation: A Double-Conformational Selection Model for the Tandem cAMP-Binding Domains of PKA RIα

Protein Kinase A (PKA) is the major receptor for the cyclic adenosine monophosphate (cAMP) secondary messenger in eukaryotes. cAMP binds to two tandem cAMP-binding domains (CBD-A and -B) within the regulatory subunit of PKA (R), unleashing the activity of the catalytic subunit (C). While CBD-A in RIα is required for PKA inhibition and activation, CBD-B functions as a “gatekeeper” domain that modulates the control exerted by CBD-A. Preliminary evidence suggests that CBD-B dynamics are critical for its gatekeeper function. To test this hypothesis, here we investigate by Nuclear Magnetic Resonance (NMR) the two-domain construct RIα (91–379) in its apo, cAMP₂, and C-bound forms. Our comparative NMR analyses lead to a double conformational selection model in which each apo CBD dynamically samples both active and inactive states independently of the adjacent CBD within a nearly degenerate free energy landscape. Such degeneracy is critical to explain the sensitivity of CBD-B to weak interactions with C and its high affinity for cAMP. Binding of cAMP eliminates this degeneracy, as it selectively stabilizes the active conformation within each CBD and inter-CBD contacts, which require both cAMP and W260. The latter is contributed by CBD-B and mediates capping of the cAMP bound to CBD-A. The inter-CBD interface is dispensable for intra-CBD conformational selection, but is indispensable for full activation of PKA as it occludes C-subunit recognition sites within CBD-A. In addition, the two structurally homologous cAMP-bound CBDs exhibit marked differences in their residual dynamics profiles, supporting the notion that conservation of structure does not necessarily imply conservation of dynamics.

Protein Kinase A (PKA) is the major receptor for the cAMP secondary messenger in eukaryotes. This study shows how PKA's regulatory subunit dynamically samples a degenerate free energy landscape that controls affinities for the catalytic subunit and cAMP; intra-domain conformational selection by cAMP controls inter-domain interactions and PKA activation.

Cyclic adenosine monophosphate (cAMP) is a messenger molecule produced within cells to control cellular metabolism in response to external stimuli. Protein Kinase A (PKA) is the major receptor for cAMP. cAMP binds to tandem cAMP-binding domains (CBD-A and -B) within the regulatory subunits of PKA (R), unleashing the activity of the catalytic subunit (C). While CBD-A is required for C-subunit inhibition and activation, in RIα CBD-B functions as a “gatekeeper” domain that modulates the control exerted by CBD-A. However, it is not currently clear how ligand binding and dynamics of CBD-B mediate its gatekeeper function. We comparatively analyzed by Nuclear Magnetic Resonance (NMR) a two-domain construct of the regulatory subunit RIα with no ligand, with cAMP₂ bound, and the C-bound form. These data show that both CBDs can exist in a system of uncorrelated conformational selection as both can independently sample activated and inactivated states (in what is known as a nearly degenerate free energy landscape). This explains why both RIα CBDs exhibit a higher cAMP-affinity than other cAMP receptors. Once cAMP has bound, the degeneracy is lost and dissociation of the kinase subunit is promoted through a combination of intra-domain conformational selection and changes in inter-CBD orientation. The proposed model—a double-conformational selection model—provides a general framework to interpret the effect of PKA mutations that have been reported in rare human disorders such as Carney complex and Acrodysostosis.

Development of a Reagentless Biosensor for Inorganic Phosphate, Applicable over a Wide Concentration Range

A fluorescent reagentless biosensor for inorganic phosphate (Pi), based on the E. coli PstS phosphate binding protein, was redesigned to allow measurements of higher Pi concentrations and at low, substoichiometric concentrations of biosensor. This was achieved by weakening Pi binding of the previous biosensor, and different approaches are described that could enable this change in properties. The readout, providing response to the Pi concentration, is delivered by tetramethylrhodamine fluorescence. In addition to two cysteine mutations for rhodamine labeling at positions 17 and 197, the final variant had an I76G mutation in the hinge region between the two lobes that make up the protein. Upon Pi binding, the lobes rotate on this hinge and the mutation on the hinge lowers affinity ∼200-fold, with a dissociation constant now in the tens to hundreds micromolar range, depending on solution conditions. The signal change on Pi binding was up to 9-fold, depending on pH. The suitability of the biosensor for steady-state ATPase assays was demonstrated with low biosensor usage and its advantage in ability to cope with Pi contamination.

Solution NMR studies of periplasmic binding proteins and their interaction partners

Periplasmic binding proteins (PBPs) are a crucial part of ATP-binding cassette import systems in Gram-negative bacteria. Central to their function is the ability to undergo a large-scale conformational rearrangement from open-unliganded to closed-liganded, which signals the presence of substrate and starts its translocation. Over the years, PBPs have been extensively studied not only owing to their essential role in nutrient uptake but also because they serve as excellent models for both practical applications (e.g., biosensor technology) and basic research (e.g., allosteric mechanisms). Although much of our knowledge at atomic level has been inferred from the detailed, static pictures afforded by crystallographic studies, nuclear magnetic resonance (NMR) has been able to fill certain gaps in such body of work, particularly with regard to dynamic processes. Here, we review NMR studies on PBPs, and their unique insights on conformation, dynamics, energetics, substrate binding, and interactions with related transport proteins. Based on the analysis of recent paramagnetic NMR results, as well as crystallographic and functional observations, we propose a mechanism that could explain the ability of certain PBPs to achieve a closed conformation in absence of ligand while others seem to remain open until ligand-mediated closure.

Binding proteins enhance specific uptake rate by increasing the substrate-transporter encounter rate

Microorganisms rely on binding-protein assisted, active transport systems to scavenge for scarce nutrients. Several advantages of using binding proteins in such uptake systems have been proposed. However, a systematic, rigorous and quantitative analysis of the function of binding proteins is lacking. By combining knowledge of selection pressure and physiochemical constraints, we derive kinetic, thermodynamic, and stoichiometric properties of binding-protein dependent transport systems that enable a maximal import activity per amount of transporter. Under the hypothesis that this maximal specific activity of the transport complex is the selection objective, binding protein concentrations should exceed the concentration of both the scarce nutrient and the transporter. This increases the encounter rate of transporter with loaded binding protein at low substrate concentrations, thereby enhancing the affinity and specific uptake rate. These predictions are experimentally testable, and a number of observations confirm them.

Rational design of a fusion protein to exhibit disulfide-mediated logic gate behavior

Synthetic cellular logic gates are primarily built from gene circuits owing to their inherent modularity. Single proteins can also possess logic gate functions and offer the potential to be simpler, quicker, and less dependent on cellular resources than gene circuits. However, the design of protein logic gates that are modular and integrate with other cellular components is a considerable challenge. As a step toward addressing this challenge, we describe the design, construction, and characterization of AND, ORN, and YES logic gates built by introducing disulfide bonds into RG13, a fusion of maltose binding protein and TEM-1 β-lactamase for which maltose is an allosteric activator of enzyme activity. We rationally designed these disulfide bonds to manipulate RG13’s allosteric regulation mechanism such that the gating had maltose and reducing agents as input signals, and the gates could be toggled between different gating functions using redox agents, although some gates performed suboptimally.

Role of water in ligand binding to maltose-binding protein: insight from a new docking protocol based on the 3D-RISM-KH molecular theory of solvation

Maltose-binding protein is a periplasmic binding protein responsible for transport of maltooligosaccarides through the periplasmic space of Gram-negative bacteria, as a part of the ABC transport system. The molecular mechanisms of the initial ligand binding and induced large scale motion of the protein's domains still remain elusive. In this study, we use a new docking protocol that combines a recently proposed explicit water placement algorithm based on the 3D-RISM-KH molecular theory of solvation and conventional docking software (AutoDock Vina) to explain the mechanisms of maltotriose binding to the apo-open state of a maltose-binding protein. We confirm the predictions of previous NMR spectroscopic experiments on binding modes of the ligand. We provide the molecular details on the binding mode that was not previously observed in the X-ray experiments. We show that this mode, which is defined by the fine balance between the protein-ligand direct interactions and solvation effects, can trigger the protein's domain motion resulting in the holo-closed structure of the maltose-binding protein with the maltotriose ligand in excellent agreement with the experimental data. We also discuss the role of water in blocking unfavorable binding sites and water-mediated interactions contributing to the stability of observable binding modes of maltotriose.

Both protein dynamics and ligand concentration can shift the binding mechanism between conformational selection and induced fit

This study aimed to shed light on the long debate over whether conformational selection (CS) or induced fit (IF) is the governing mechanism for protein-ligand binding. The main difference between the two scenarios is whether the conformational transition of the protein from the unbound form to the bound form occurs before or after encountering the ligand. Here we introduce the IF fraction (i.e., the fraction of binding events achieved via IF), to quantify the binding mechanism. Using simulations of a model protein-ligand system, we demonstrate that both the rate of the conformational transition and the concentration of ligand molecules can affect the IF fraction. CS dominates at slow conformational transition and low ligand concentration. An increase in either quantity results in a higher IF fraction. Despite the many-body nature of the system and the involvement of multiple, disparate types of dynamics (i.e., ligand diffusion, protein conformational transition, and binding reaction), the overall binding kinetics over wide ranges of parameters can be fit to a single exponential, with the apparent rate constant exhibiting a linear dependence on ligand concentration. The present study may guide future kinetics experiments and dynamics simulations in determining the IF fraction.

Protein conformational dynamics dictate the binding affinity for a ligand

Interactions between a protein and a ligand are essential to all biological processes. Binding and dissociation are the two fundamental steps of ligand-protein interactions, and determine the binding affinity. Intrinsic conformational dynamics of proteins have been suggested to play crucial roles in ligand binding and dissociation. Here, we demonstrate how protein dynamics dictate the binding and dissociation of a ligand through a single-molecule kinetic analysis for a series of maltose-binding protein mutants that have different intrinsic conformational dynamics and dissociation constants for maltose. Our results provide direct evidence that the ligand dissociation is determined by the intrinsic opening rate of the protein.

Subtle dynamics of holo glutamine binding protein revealed with a rigid paramagnetic probe

Bacterial periplasmic binding proteins (PBPs) are involved in the translocation of small molecules in the periplasm. To unload, the two domains of a PBP open up, allowing the ligand to exit. However, it is not clear whether there are dynamics near the binding site which can facilitate the rapid dissociation of a ligand. To visualize such dynamics, we utilized paramagnetic relaxation enhancement (PRE) NMR and introduced a rigid paramagnetic probe to a PBP, glutamine-binding protein (QBP) with its cognate ligand bound. A paramagnetic Cu(II) ion is sandwiched between an engineered di-histidine motif at a helix and an NTA capping molecule. The afforded paramagnetic probe is so rigid that PRE values calculated from a single structure of holo QBP largely agree with the observed values. The remaining PRE discrepancies, however, manifest dynamics of a loop in the opposite domain from the paramagnetic probe. This loop packs against the glutamine ligand in the holo QBP and undergoes fluctuations upon ligand dissociation, as assessed by steered molecular dynamics simulations. As such, the loop dynamics, occurring for a small population in nanosecond to microsecond time scale, may be related to the ligand dissociation process. The rigid paramagnetic probe described herein can be grafted to other protein systems for structure and dynamics studies.

Competitive interactions of ligands and macromolecular crowders with maltose binding protein

Cellular signaling involves a cascade of recognition events occurring in a complex environment with high concentrations of proteins, polysaccharides, and other macromolecules. The influence of macromolecular crowders on protein binding affinity through hard-core repulsion is well studied, and possible contributions of protein-crowder soft attraction have been implicated recently. Here we present direct evidence for weak association of maltose binding protein (MBP) with a polysaccharide crowder Ficoll, and that this association effectively competes with the binding of the natural ligand, maltose. Titration data over wide ranges of maltose and Ficoll concentrations fit well with a three-state competitive binding model. Broadening of MBP ¹H­¹⁵N TROSY spectra by the addition of Ficoll indicates weak protein-crowder association, and subsequent recovery of sharp NMR peaks upon addition of maltose indicates that the interactions of the crowder and the ligand with MBP are competitive. We hypothesize that, in the Escherichia coli periplasm, the competitive interactions of polysaccharides and maltose with MBP could allow MBP to shuttle between the peptidoglycan attached to the outer membrane and the ATP-binding cassette transporter in the inner membrane.

Molecular details of ligand selectivity determinants in a promiscuous β-glucan periplasmic binding protein

BACKGROUND

Members of the periplasmic binding protein (PBP) superfamily utilize a highly conserved inter-domain ligand binding site that adapts to specifically bind a chemically diverse range of ligands. This paradigm of PBP ligand binding specificity was recently altered when the structure of the Thermotoga maritima cellobiose-binding protein (tmCBP) was solved. The tmCBP binding site is bipartite, comprising a canonical solvent-excluded region (subsite one), adjacent to a solvent-filled cavity (subsite two) where specific and semi-specific ligand recognition occur, respectively.

RESULTS

A molecular level understanding of binding pocket adaptation mechanisms that simultaneously allow both ligand specificity at subsite one and promiscuity at subsite two has potentially important implications in ligand binding and drug design studies. We sought to investigate the determinants of ligand binding selectivity in tmCBP through biophysical characterization of tmCBP in the presence of varying β-glucan oligosaccharides. Crystal structures show that whilst the amino acids that comprise both the tmCBP subsite one and subsite two binding sites remain fixed in conformation regardless of which ligands are present, the rich hydrogen bonding potential of water molecules may facilitate the ordering and the plasticity of this unique PBP binding site.

CONCLUSIONS

The identification of the roles these water molecules play in ligand recognition suggests potential mechanisms that can be utilized to adapt a single ligand binding site to recognize multiple distinct ligands.

Interplay between conformational selection and induced fit in multidomain protein-ligand binding probed by paramagnetic relaxation enhancement

The binding of ligands and substrates to proteins has been extensively studied for many years and can be described, in its simplest form, by two limiting mechanisms: conformational selection and induced fit. Conformational selection involves the binding of ligand to a pre-existing sparsely-populated conformation of the free protein that is the same as that in the final protein-ligand complex. In the case of induced fit, the ligand binds to the major conformation of the free protein and only subsequent to binding undergoes a conformational change to the final protein-ligand complex. While these two mechanisms can be dissected and distinguished by transient kinetic measurements, direct direction, characterization and visualization of transient, sparsely-populated states of proteins are experimentally challenging. Unless trapped, sparsely-populated states are generally invisible to conventional structural and biophysical techniques, including crystallography and most NMR measurements. In this review we summarize some recent developments in the use of paramagnetic relaxation enhancement to directly study sparsely-populated states of proteins and illustrate the application of this approach to two proteins, maltose binding protein and calmodulin, both of which undergo large rigid body conformational rearrangements upon ligand binding from an open apo state to a closed ligand-bound holo state. We show that the apo state ensemble comprises a small population of partially-closed configurations that are similar but not identical to that of the holo state. These results highlight the complementarity and interplay of induced fit and conformational selection and suggest that the existence of partially-closed states in the absence of ligand facilitates the transition to the closed ligand-bound state.

Non-allosteric enzyme switches possess larger effector-induced changes in thermodynamic stability than their non-switch analogs

The ability to regulate cellular protein activity offers a broad range of biotechnological and biomedical applications. Such protein regulation can be achieved by modulating the specific protein activity or through processes that regulate the amount of protein in the cell. We have previously demonstrated that the nonhomologous recombination of the genes encoding maltose binding protein (MBP) and TEM1 β-lactamase (BLA) can result in genes that confer maltose-dependent resistance to β-lactam antibiotics even though the encoded proteins are not allosteric enzymes. We showed that these phenotypic switches-named based on their conferral of a switching phenotype to cells-resulted from a specific interaction with maltose in the cell that increased the switches cellular accumulation. Since phenotypic switches represent an important class of engineered proteins for basic science and biotechnological applications in vivo, we sought to elucidate the phenomena behind the increased accumulation and switching properties. Here, we demonstrate the key role for the linker region between the two proteins. Experimental evidence supports the hypothesis that in the absence of their effector, some phenotypic switches possess an increased rate of unfolding, decreased conformational stability, and increased protease susceptibility. These factors alone or in combination serve to decrease cellular accumulation. The effector functions to increase cellular accumulation by alleviating one or more of these defects. This perspective on the mechanism for phenotypic switching will aid the development of design rules for switch construction for applications and inform the study of the regulatory mechanisms of natural cellular proteins.

Allosteric signaling and dynamics of the clamshell-like NMDA receptor GluN1 N-terminal domain

N-methyl-D-aspartate receptors (NMDARs), neuronal glutamate-gated ion channels, are obligatory heterotetramers composed of GluN1 and GluN2 subunits. Each subunit contains two extracellular clamshell-like domains with an agonist-binding domain and a distal N-terminal domain (NTD). The GluN2 NTDs form mobile regulatory domains. In contrast, the dynamics of GluN1 NTD and its contribution to NMDAR function remain poorly understood. Here we show that GluN1 NTD is neither static nor functionally silent. Perturbing the conformation of GluN1 NTD affects both receptor gating and pharmacological properties. GluN1 NTD undergoes structural rearrangements that involve hinge bending and large twisting and untwisting motions, allowing for new intra- and intersubunit contacts. GluN1 NTD acts in trans with GluN2 NTD to influence binding of glutamate but, notably, not of GluN1 coagonist glycine. Our work uncovers a dynamic role of GluN1 NTD in controlling NMDAR function through new interdomain allosteric interactions.

Mapping the conformational stability of maltose binding protein at the residue scale using nuclear magnetic resonance hydrogen exchange experiments

Being able to differentiate local fluctuations from global folding-unfolding dynamics of a protein is of major interest for improving our understanding of structure-function determinants. The maltose binding protein (MBP), a protein that belongs to the maltose transport system, has a structure composed of two globular domains separated by a rigid-body "hinge bending". Here we determined, by using hydrogen exchange (HX) nuclear magnetic resonance experiments, the apparent stabilization free energies of 101 residues of MBP bound to β-cyclodextrin (MBP-βCD) under native conditions. We observed that the last helix of MBP (helix α14) has a lower protection factor than the rest of the protein. Further, HX experiments were performed using guanidine hydrochloride under subdenaturing conditions to discriminate between local fluctuations and global unfolding events and to determine the MBP-βCD energy landscape. The results show that helix α4 and a part of helices α5 and α6 are clearly grouped into a subdenaturing folding unit and represent a partially folded intermediate under native conditions. In addition, we observed that amide protons located in the hinge between the two globular domains share similar ΔG(gu)(app) and m values and should unfold simultaneously. These observations provide new points of view for improving our understanding of the thermodynamic stability and the mechanisms that drive folding-unfolding dynamics of proteins.

A master switch couples Mg²⁺-assisted catalysis to domain motion in B. stearothermophilus tryptophanyl-tRNA Synthetase

We demonstrate how tryptophanyl-tRNA synthetase uses conformation-dependent Mg(2+) activation to couple catalysis of tryptophan activation to specific, functional domain movements. Rate acceleration by Mg(2+) requires ∼-6.0 kcal/mol in protein⋅Mg(2+) interaction energy, none of which arises from the active site. A highly cooperative interaction between Mg(2+) and four residues from a remote, conserved motif that mediates the shear of domain movement (1) destabilizes the pretransition state conformation, thereby (2) inducing the Mg(2+) to stabilize the transition state for k(cat) by ∼-5.0 kcal/mol. Cooperative, long-range conformational effects on the metal therefore convert an inactive Mg(2+) coordination into one that can stabilize the transition state if, and only if, domain motion occurs. Transient, conformation-dependent Mg(2+) activation, analogous to the escapement in mechanical clocks, explains vectorial coupling.

Exploration of multi-state conformational dynamics and underlying global functional landscape of maltose binding protein

An increasing number of biological machines have been revealed to have more than two macroscopic states. Quantifying the underlying multiple-basin functional landscape is essential for understanding their functions. However, the present models seem to be insufficient to describe such multiple-state systems. To meet this challenge, we have developed a coarse grained triple-basin structure-based model with implicit ligand. Based on our model, the constructed functional landscape is sufficiently sampled by the brute-force molecular dynamics simulation. We explored maltose-binding protein (MBP) which undergoes large-scale domain motion between open, apo-closed (partially closed) and holo-closed (fully closed) states responding to ligand binding. We revealed an underlying mechanism whereby major induced fit and minor population shift pathways co-exist by quantitative flux analysis. We found that the hinge regions play an important role in the functional dynamics as well as that increases in its flexibility promote population shifts. This finding provides a theoretical explanation of the mechanistic discrepancies in PBP protein family. We also found a functional “backtracking” behavior that favors conformational change. We further explored the underlying folding landscape in response to ligand binding. Consistent with earlier experimental findings, the presence of ligand increases the cooperativity and stability of MBP. This work provides the first study to explore the folding dynamics and functional dynamics under the same theoretical framework using our triple-basin functional model.

A central goal of biology is to understand the function of the organism and its constituent parts at each of its scales of complexity. Function at the molecular level is often realized by changes in conformation. Unfortunately, experimental explorations of global motions critical for functional conformational changes are still challenging. In the present work, we developed a coarse grained triple-well structure-based model to explore the underlying functional landscape of maltose-binding protein (MBP). By quantitative flux analysis, we uncover the underlying mechanism by which the major induced fit and minor population shift pathways co-exist. Though we have previously lent credence to the assertion that dynamical equilibrium between open and minor closed conformations exist for all the free PBPs, the generality of this rule is still a matter of open debate. We found that the hinge flexibility is favorable to population shift mechanism. This finding provides a theoretical explanation of the mechanism discrepancies in PBP protein family. We also simulated the folding dynamics using this functional multi-basin model which successfully reproduced earlier protein melting experiment. This represents an exciting opportunity to characterize the interplay between folding and function, which is a long-standing question in the community. The theoretical approach employed in this study is general and can be applied to other systems.

Discovery of an auto-regulation mechanism for the maltose ABC transporter MalFGK2

The maltose transporter MalFGK₂, together with the substrate-binding protein MalE, is one of the best-characterized ABC transporters. In the conventional model, MalE captures maltose in the periplasm and delivers the sugar to the transporter. Here, using nanodiscs and proteoliposomes, we instead find that MalE is bound with high-affinity to MalFGK2 to facilitate the acquisition of the sugar. When the maltose concentration exceeds the transport capacity, MalE captures maltose and dissociates from the transporter. This mechanism explains why the transport rate is high when MalE has low affinity for maltose, and low when MalE has high affinity for maltose. Transporter-bound MalE facilitates the acquisition of the sugar at low concentrations, but also captures and dissociates from the transporter past a threshold maltose concentration. In vivo, this maltose-forced dissociation limits the rate of transport. Given the conservation of the substrate-binding proteins, this mode of allosteric regulation may be universal to ABC importers.

Structure-based model of allostery predicts coupling between distant sites

Allostery is a phenomenon that couples effector ligand binding at an allosteric site to a structural and/or dynamic change at a distant regulated site. To study an allosteric transition, we vary the size of the allosteric site and its interactions to construct a series of energy landscapes with pronounced minima corresponding to both the effector bound and unbound crystal structures. We use molecular dynamics to sample these landscapes. The degree of perturbation by the effector, modeled by the size of the allosteric site, provides an order parameter for allostery that allows us to determine how microscopic motions give rise to commonly discussed macroscopic mechanisms: (i) induced fit, (ii) population shift, and (iii) entropy driven. These mechanisms involve decreasing structural differences between the effector bound and unbound populations. A metric (ligand-induced cooperativity) can measure how cooperatively a given regulated site responds to effector binding and therefore what kind of allosteric mechanism is involved. We apply the model to three proteins with experimentally characterized transitions: (i) calmodulin-GFP Ca(2+) sensor protein, (ii) maltose binding protein, and (iii) CSL transcription factor. Remarkably, the model is able to reproduce allosteric motion and predict coupling in a manner consistent with experiment.

Structural, thermodynamic, and mechanistical studies in uroporphyrinogen III synthase: molecular basis of congenital erythropoietic porphyria

Congenital erythropoietic porphyria (CEP) is a rare autosomal disease ultimately related to deleterious mutations in uroporphyrinogen III synthase (UROIIIS), the fourth enzyme of the biosynthetic route of the heme group. UROIIIS catalyzes the cyclization of the linear tetrapyrrol hydroxymethylbilane (HMB), inverting the configuration in one of the aromatic rings. In the absence of the enzyme (or when ill-functioning), HMB spontaneously degrades to the by-product uroporphyrinogen I, which cannot lead to the heme group and accumulates in the body, producing some of the symptoms observed in CEP patients. In the present chapter, clinical, biochemical, and biophysical information has been compiled to provide an integrative view on the molecular basis of CEP. The high-resolution structure of UROIIIS sheds light on the enzyme reaction mechanism while thermodynamic analysis revealed that the protein is thermolabile. Pathogenic missense mutations are found throughout the primary sequence of the enzyme. All but one of these is rarely found in patients, whereas C73R is responsible for more than one-third of the reported cases. Most of the mutant proteins (C73R included) retain partial catalytic activity but the mutations often reduce the enzyme's stability. The stabilization of the protein in vivo is discussed in the context of a new line of intervention to complement existing treatments such as bone marrow transplantation and gene therapy.