Kinetic Barrier to Enzyme Inhibition Is Manipulated by Dynamical Local Interactions in E. coli DHFR.

Dihydrofolate reductase (DHFR) is an important drug target and a highly studied model protein for understanding enzyme dynamics. DHFR's crucial role in folate synthesis renders it an ideal candidate to understand protein function and protein evolution mechanisms. In this study, to understand how a newly proposed DHFR inhibitor, 4'-deoxy methyl trimethoprim (4'-DTMP), alters evolutionary trajectories, we studied interactions that lead to its superior performance over that of trimethoprim (TMP). To elucidate the inhibition mechanism of 4'-DTMP, we first confirmed, both computationally and experimentally, that the relative binding free energy cost for the mutation of TMP and 4'-DTMP is the same, pointing the origin of the characteristic differences to be kinetic rather than thermodynamic. We then employed an interaction-based analysis by focusing first on the active site and then on the whole enzyme. We confirmed that the polar modification in 4'-DTMP induces additional local interactions with the enzyme, particularly, the M20 loop. These changes are propagated to the whole enzyme as shifts in the hydrogen bond networks. To shed light on the allosteric interactions, we support our analysis with network-based community analysis and show that segmentation of the loop domain of inhibitor-bound DHFR must be avoided by a successful inhibitor.

Conformational multiplicity of bacterial ferric binding protein revealed by small angle x-ray scattering and molecular dynamics calculations.

This study combines molecular dynamics (MD) simulations with small angle x-ray scattering (SAXS) measurements to investigate the range of conformations that can be adopted by a pH/ionic strength (IS) sensitive protein and to quantify its distinct populations in solution. To explore how the conformational distribution of proteins may be modified in the environmental niches of biological media, we focus on the periplasmic ferric binding protein A (FbpA) from Haemophilus influenzae involved in the mechanism by which bacteria capture iron from higher organisms. We examine iron-binding/release mechanisms of FbpA in varying conditions simulating its biological environment. While we show that these changes fall within the detectable range for SAXS as evidenced by differences observed in the theoretical scattering patterns calculated from the crystal structure models of apo and holo forms, detection of conformational changes due to the point mutation D52A and changes in ionic strength (IS) from SAXS scattering profiles have been challenging. Here, to reach conclusions, statistical analyses with SAXS profiles and results from different techniques were combined in a complementary fashion. The SAXS data complemented by size exclusion chromatography point to multiple and/or alternative conformations at physiological IS, whereas they are well-explained by single crystallographic structures in low IS buffers. By fitting the SAXS data with unique conformations sampled by a series of MD simulations under conditions mimicking the buffers, we quantify the populations of the occupied substates. We also find that the D52A mutant that we predicted by coarse-grained computational modeling to allosterically control the iron binding site in FbpA, responds to the environmental changes in our experiments with conformational selection scenarios that differ from those of the wild type.

DHFR Mutants Modulate Their Synchronized Dynamics with the Substrate by Shifting Hydrogen Bond Occupancies.

Antibiotic resistance is a global health problem in which mutations occurring in functional proteins render drugs ineffective. The working mechanisms of the arising mutants are seldom apparent; a methodology to decipher these mechanisms systematically would render devising therapies to control the arising mutational pathways possible. Here we utilize Cα-Cß bond vector relaxations obtained from moderate length MD trajectories to determine conduits for functionality of the resistance conferring mutants of Escherichia coli dihydrofolate reductase. We find that the whole enzyme is synchronized to the motions of the substrate, irrespective of the mutation introducing gain-of-function or loss-of function. The total coordination of the motions suggests changes in the hydrogen bond dynamics with respect to the wild type as a possible route to determine and classify the mode-of-action of individual mutants. As a result, nine trimethoprim-resistant point mutations arising frequently in evolution experiments are categorized. One group of mutants that display the largest occurrence (L28R, W30G) work directly by modifying the dihydrofolate binding region. Conversely, W30R works indirectly by the formation of the E139-R30 salt bridge which releases energy resulting from tight binding by distorting the binding cavity. A third group (D27E, F153S, I94L) arising as single, resistance invoking mutants in evolution experiment trajectories allosterically and dynamically affects a hydrogen bonding motif formed at residues 59-69-71 which in turn modifies the binding site dynamics. The final group (I5F, A26T, R98P) consists of those mutants that have properties most similar to the wild type; these only appear after one of the other mutants is fixed on the protein structure and therefore display clear epistasis. Thus, we show that the binding event is governed by the entire enzyme dynamics while the binding site residues play gating roles. The adjustments made in the total enzyme in response to point mutations are what make quantifying and pinpointing their effect a hard problem. Here, we show that hydrogen bond dynamics recorded on sub-µs time scales provide the necessary fingerprints to decipher the various mechanisms at play.

N-Terminus of the Third PDZ Domain of PSD-95 Orchestrates Allosteric Communication for Selective Ligand Binding.

PDZ domains constitute common models to study single-domain allostery without significant structural changes. The third PDZ domain of PSD-95 (PDZ3) is known to have selective structural features that confer unique modulatory roles to this unit. In this model system, two residues, H372 directly connected to the binding site and G330 holding an off-binding-site position, were designated to assess the effect of mutations on binding selectivity. It has been observed that the H372A and G330T-H372A mutations change ligand preferences from class I (T/S amino acid at position -2 of the ligand) to class II (hydrophobic amino acid at the same position). Alternatively, the G330T single mutation leads to the recognition of both ligand classes. We have performed a series of molecular dynamics (MD) simulations for wild-type, H372A, and G330T single mutants and a double mutant of PDZ3 in the absence and presence of both types of ligands. With the combination of free-energy difference calculations and a detailed analysis of MD trajectories, "class switching" and "class bridging" behavior of PDZ3 mutants, as well as their effects on ligand selection and binding affinities are explained. We show that the dynamics of the charged N-terminus plays a fundamental role in determining the binding preferences in PDZ3 by altering the electrostatic energy. These findings are corroborated by simulations on N-terminus-truncated versions of these systems. The dynamical allostery orchestrated by the N-terminus offers a fresh perspective to the study of communication pathways in proteins.

High-Order Epistasis in Catalytic Power of Dihydrofolate Reductase Gives Rise to a Rugged Fitness Landscape in the Presence of Trimethoprim Selection.

Evolutionary fitness landscapes of several antibiotic target proteins have been comprehensively mapped showing strong high-order epistasis between mutations, but understanding these effects at the biochemical and structural levels remained open. Here, we carried out an extensive experimental and computational study to quantitatively understand the evolutionary dynamics of Escherichia coli dihydrofolate reductase (DHFR) enzyme in the presence of trimethoprim-induced selection. To facilitate this, we developed a new in vitro assay for rapidly characterizing DHFR steady-state kinetics. Biochemical and structural characterization of resistance-conferring mutations targeting a total of ten residues spanning the substrate binding pocket of DHFR revealed distinct changes in the catalytic efficiencies of mutated DHFR enzymes. Next, we measured biochemical parameters (Km, Ki, and kcat) for a mutant library carrying all possible combinations of six resistance-conferring DHFR mutations and quantified epistatic interactions between them. We found that the high-order epistasis in catalytic power of DHFR (kcat and Km) creates a rugged fitness landscape under trimethoprim selection. Taken together, our data provide a concrete illustration of how epistatic coupling at the level of biochemical parameters can give rise to complex fitness landscapes, and suggest new strategies for developing mutant specific inhibitors.

MODE-TASK: large-scale protein motion tools.

Summary: MODE-TASK, a novel and versatile software suite, comprises Principal Component Analysis, Multidimensional Scaling, and t-Distributed Stochastic Neighbor Embedding techniques using Molecular Dynamics trajectories. MODE-TASK also includes a Normal Mode Analysis tool based on Anisotropic Network Model so as to provide a variety of ways to analyse and compare large-scale motions of protein complexes for which long MD simulations are prohibitive. Beside the command line function, a GUI has been developed as a PyMOL plugin. Availability and implementation: MODE-TASK is open source, and available for download from https://github.com/RUBi-ZA/MODE-TASK. It is implemented in Python and C++. It is compatible with Python 2.x and Python 3.x and can be installed by Conda. Supplementary information: Supplementary data are available at Bioinformatics online.

Unraveling the Motions behind Enterovirus 71 Uncoating.

Enterovirus 71 can be a severe pathogen in small children and immunocompromised adults. Virus uncoating is a critical step in the infection of the host cell; however, the mechanisms that control this process remain poorly understood. We applied normal mode analysis and perturbation response scanning to several complexes of the virus capsid and present a coarse-graining approach to analyze the full capsid. We show that our method offers an alternative to expressing the system as a set of rigid blocks and accounts for the interconnection between nodes within each subunit and protein interfaces across the capsid. In our coarse-grained approach, the modes associated with capsid expansion are captured in the first three nondegenerate modes and correspond to the changes observed in structural studies of the virus. We show that the resolution of the analysis may be modified without losing information on the global motions leading to uncoating. Perturbation response scanning revealed that a protomer cannot serve as a functional unit to explain deformations of the capsid. Instead, we define a pentamer as the minimum functional unit to investigate changes within the capsid. From the modal analysis and perturbation response scanning, we locate a hotspot region surrounding the fivefold axis. The range of the effect of these single, hotspot residues extend to 140 Å. The perturbation of internal capsid residues in this region displayed greatest propensity to capsid expansion, thus indicating the significant role that the RNA genome may play in triggering uncoating.

MD-TASK: a software suite for analyzing molecular dynamics trajectories.

SUMMARY: Molecular dynamics (MD) determines the physical motions of atoms of a biological macromolecule in a cell-like environment and is an important method in structural bioinformatics. Traditionally, measurements such as root mean square deviation, root mean square fluctuation, radius of gyration, and various energy measures have been used to analyze MD simulations. Here, we present MD-TASK, a novel software suite that employs graph theory techniques, perturbation response scanning, and dynamic cross-correlation to provide unique ways for analyzing MD trajectories. AVAILABILITY AND IMPLEMENTATION: MD-TASK has been open-sourced and is available for download from https://github.com/RUBi-ZA/MD-TASK , implemented in Python and supported on Linux/Unix. CONTACT: o.tastanbishop@ru.ac.za.

FbpA iron storage and release are governed by periplasmic microenvironments.

Ferric binding protein (FbpA) is part of an elaborate iron piracy mechanism evolved in Gram-negative bacteria, shuttling iron in the periplasmic space, from the outer to the cytoplasmic membrane side. We address how the dissociation process of iron is facilitated, since the binding constant of iron is on the order of 1018 M-1 at 6.5 pH and 200 mM ionic strength (IS). We monitor the conformational preferences of FbpA by extensive molecular dynamics (MD) simulations under conditions where IS, charge states of iron coordinating tyrosines and pH are varied, as well as when a mutation is introduced at an allosteric site. Steered MD is utilized to predict the binding affinity of iron. After triggering lobe opening by changing the charge states of tyrosines, the conformations adopted and the iron binding affinity still depend on pH, IS and allosteric interactions. To relate the observed conformational changes to the environmental conditions that might be encountered in the periplasmic space, we offer a plausible model that couples electrostatic potential distribution to the mechanical motions invoked. Although low pH/IS and allosteric perturbations decrease the affinity of iron, it remains high for spontaneous dissociation. However, the conformational changes modulated by the environmental conditions expose iron for chelation. Our study provides a quantitative dimension and molecular details to interpret the contribution of possible environmental conditions present in the periplasmic space to iron dissociation from FbpA, opening up the opportunity of modulating function via allosteric mutations or altering environmental conditions, thus offering a new route to developing strategies towards antibiotic resistance by targeting nutritional requirements.

Increased substrate affinity in the Escherichia coli L28R dihydrofolate reductase mutant causes trimethoprim resistance.

Dihydrofolate reductase (DHFR) is a ubiquitous enzyme with an essential role in cell metabolism. DHFR catalyzes the reduction of dihydrofolate to tetrahydrofolate, which is a precursor for purine and thymidylate synthesis. Several DHFR targeting antifolate drugs including trimethoprim, a competitive antibacterial inhibitor, have therefore been developed and are clinically used. Evolution of resistance against antifolates is a common public health problem rendering these drugs ineffective. To combat the resistance problem, it is important to understand resistance-conferring changes in the DHFR structure and accordingly develop alternative strategies. Here, we structurally and dynamically characterize Escherichia coli DHFR in its wild type (WT) and trimethoprim resistant L28R mutant forms in the presence of the substrate and its inhibitor trimethoprim. We use molecular dynamics simulations to determine the conformational space, loop dynamics and hydrogen bond distributions at the active site of DHFR for the WT and the L28R mutant. We also report their experimental kcat, Km, and Ki values, accompanied by isothermal titration calorimetry measurements of DHFR that distinguish enthalpic and entropic contributions to trimethoprim binding. Although mutations that confer resistance to competitive inhibitors typically make enzymes more promiscuous and decrease affinity to both the substrate and the inhibitor, strikingly, we find that the L28R mutant has a unique resistance mechanism. While the binding affinity differences between the WT and the mutant for the inhibitor and the substrate are small, the newly formed extra hydrogen bonds with the aminobenzoyl glutamate tail of DHF in the L28R mutant leads to increased barriers for the dissociation of the substrate and the product. Therefore, the L28R mutant indirectly gains resistance by enjoying prolonged binding times in the enzyme-substrate complex. While this also leads to slower product release and decreases the catalytic rate of the L28R mutant, the overall effect is the maintenance of a sufficient product formation rate. Finally, the experimental and computational analyses together reveal the changes that occur in the energetic landscape of DHFR upon the resistance-conferring L28R mutation. We show that the negative entropy associated with the binding of trimethoprim in WT DHFR is due to water organization at the binding interface. Our study lays the framework to study structural changes in other trimethoprim resistant DHFR mutants.

Detailed molecular dynamics simulations of human transferrin provide insights into iron release dynamics at serum and endosomal pH.

Human serum transferrin (hTf) transports ferric ions in the blood stream and delivers them to cells via receptor-mediated endocytosis. hTf is folded into two homologous lobes; we utilize three of the available crystal structures delineating large conformational changes involved in iron binding/dissociation. We address the problems of whether the release process follows the same trend at serum (~7.4) and endosomal (~5.6) pH, and if there is communication between the lobes. In the absence of the transferrin receptor, we study the dynamics of the full structure as well as the separate lobes in different closed, partially open, and open conformations under neutral and endosomal pH conditions. Results corroborate those experimental observations underscoring the distinguishing effect of pH on the dynamics of hTf. Furthermore, in a total of 2 µs molecular dynamics simulations, residue fluctuations elucidate the cross talk between the lobes correlated by the peptide linker bridging them at serum pH, while their correlations are lost under endosomal conditions. At serum pH, interplay between relative mobility of the lobes is correlated with iron release rates, rendering the initial conformational change an important contributor to the dynamics under these conditions. Interestingly, C-lobe opening lags behind that of the N-lobe as long as there is at least one iron bound, making the more stable C-lobe an attractive target for recognition by receptors. At endosomal pH, both lobes readily open, making irons available for delivery.

Perturbation response scanning specifies key regions in subtilisin serine protease for both function and stability.

Can one infer the amino acids of the enzymes that are responsible for the stability or the level of the catalytic activity by computationally experimenting on the inhibited enzyme in the enzyme-inhibitor complex? In this article, we answer this question positively both by designing molecular dynamics simulations and by devising coarse-grained methodologies on the subtilisin serine protease. Both methodologies are based on the cross-correlations of the fluctuations of the residues, obtained either by monitoring the trajectories from the simulation or by constructing the inverse Laplacian of the elastic network model, of the complex. A perturbation scanning is applied to the complex using these correlations. The results indicate that the two methods almost point out the same regions on the flexible of the enzyme. These regions are: (i) 50-61, (ii) 155-164 and (iii) 192-194, all of which are designated to be important by experimental studies in the literature.

Designing molecular dynamics simulations to shift populations of the conformational states of calmodulin.

We elucidate the mechanisms that lead to population shifts in the conformational states of calcium-loaded calmodulin (Ca(2+)-CaM). We design extensive molecular dynamics simulations to classify the effects that are responsible for adopting occupied conformations available in the ensemble of NMR structures. Electrostatic interactions amongst the different regions of the protein and with its vicinal water are herein mediated by lowering the ionic strength or the pH. Amino acid E31, which is one of the few charged residues whose ionization state is highly sensitive to pH differences in the physiological range, proves to be distinctive in its control of population shifts. E31A mutation at low ionic strength results in a distinct change from an extended to a compact Ca(2+)-CaM conformation within tens of nanoseconds, that otherwise occur on the time scales of microseconds. The kinked linker found in this particular compact form is observed in many of the target-bound forms of Ca(2+)-CaM, increasing the binding affinity. This mutation is unique in controlling C-lobe dynamics by affecting the fluctuations between the EF-hand motif helices. We also monitor the effect of the ionic strength on the conformational multiplicity of Ca(2+)-CaM. By lowering the ionic strength, the tendency of nonspecific anions in water to accumulate near the protein surface increases, especially in the vicinity of the linker. The change in the distribution of ions in the vicinal layer of water allows N- and C- lobes to span a wide variety of relative orientations that are otherwise not observed at physiological ionic strength. E31 protonation restores the conformations associated with physiological environmental conditions even at low ionic strength.

A stochastic programming approach to the antibiotics time machine problem.

Antibiotics Time Machine is an important problem to understand antibiotic resistance and how it can be reversed. Mathematically, it can be modeled as follows: Consider a set of genotypes, each of which contain a set of mutated and unmutated genes. Suppose that a set of growth rate measurements of each genotype under a set of antibiotics is given. The transition probabilities of a 'realization' of a Markov chain associated with each arc under each antibiotic are computable via a predefined function given the growth rate realizations. The aim is to maximize the expected probability of reaching to the genotype with all unmutated genes given the initial genotype in a predetermined number of transitions, considering the following two sources of uncertainties: (i) the randomness in growth rates, (ii) the randomness in transition probabilities, which are functions of growth rates. We develop stochastic mixed-integer linear programming and dynamic programming approaches to solve static and dynamic versions of the Antibiotics Time Machine Problem under the aforementioned uncertainties. We adapt a Sample Average Approximation approach that exploits the special structure of the problem and provide accurate solutions that perform very well in an out-of-sample analysis.

Computational strategies for protein conformational ensemble detection.

Protein function is constrained by the three-dimensional structure but is delineated by its dynamics. This framework must satisfy specificity of function along with adaptability to changing environments and evolvability under external constraints. The accessibility of the available regions of the energy landscape for a set of conditions and shifts in the populations upon their modulation have effects propagating across scales, from biomolecular interactions, to organisms, to populations. Developing the ability to detect and juggle protein conformations supplemented by a physics-based understanding has implications for not only in vivo problems but also for resistance impeding drug discovery and bionano-sensor design.

Subtle pH differences trigger single residue motions for moderating conformations of calmodulin.

This study reveals the essence of ligand recognition mechanisms by which calmodulin (CaM) controls a variety of Ca(2+) signaling processes. We study eight forms of calcium-loaded CaM each with distinct conformational states. Reducing the structure to two degrees of freedom conveniently describes main features of the conformational changes of CaM via simultaneous twist-bend motions of the two lobes. We utilize perturbation-response scanning (PRS) technique, coupled with molecular dynamics simulations. PRS is based on linear response theory, comprising sequential application of directed forces on selected residues followed by recording the resulting protein coordinates. We analyze directional preferences of the perturbations and resulting conformational changes. Manipulation of a single residue reproduces the structural change more effectively than that of single/pairs/triplets of collective modes of motion. Our findings also give information on how the flexible linker acts as a transducer of binding information to distant parts of the protein. Furthermore, by perturbing residue E31 located in one of the EF hand motifs in a specific direction, it is possible to induce conformational change relevant to five target structures. Independently, using four different pK(a) calculation strategies, we find this particular residue to be the charged residue (out of a total of 52), whose ionization state is most sensitive to subtle pH variations in the physiological range. It is plausible that at relatively low pH, CaM structure is less flexible. By gaining charged states at specific sites at a pH value around 7, such as E31 found in the present study, local conformational changes in the protein will lead to shifts in the energy landscape, paving the way to other conformational states. These findings are in accordance with Fluorescence Resonance Energy Transfer (FRET) measured shifts in conformational distributions towards more compact forms with decreased pH. They also corroborate mutational studies and proteolysis results which point to the significant role of E31 in CaM dynamics.

Dynamic Community Composition Unravels Allosteric Communication in PDZ3.

The third domain of PSD-95 (PDZ3) is a model for investigating allosteric communication in protein and ligand interactions. While motifs contributing to its binding specificity have been scrutinized, a conformational dynamical basis is yet to be established. Despite the miniscule structural changes due to point mutants, the observed significant binding affinity differences have previously been assessed with a focus on two α-helices located at the binding groove (α2) and the C-terminus (α3). Here, we employ a new computational approach to develop a generalized view on the molecular basis of PDZ3 binding selectivity and interaction communication for a set of point mutants of the protein (G330T, H372A, G330T-H372A) and its ligand (CRIPT, named L1, and its T-2F variant, L2) along with the wild type (WT). To analyze the dynamical aspects hidden in the conformations that are produced by molecular dynamics simulations, we utilize variations in community composition calculated based on the betweenness centrality measure from graph theory. We find that the highly charged N-terminus, which is located far from the ligand, has the propensity to share the same community with the ligand in the biologically functional complexes, indicating a distal segment might mediate the binding dynamics. N- and C-termini of PDZ3 share communities, and α3 acts as a hub for the whole protein by sustaining the communication with all structural segments, albeit being a trait not unique to the functional complexes. Moreover, α2 which lines the binding cavity frequently parts communities with the ligand and is not a controller of the binding but is rather a slave to the overall dynamics coordinated by the N-terminus. Thus, ligand binding fate in PDZ3 is traced to the population of community compositions extracted from dynamics despite the lack of significant conformational changes.

A trimethoprim derivative impedes antibiotic resistance evolution.

The antibiotic trimethoprim (TMP) is used to treat a variety of Escherichia coli infections, but its efficacy is limited by the rapid emergence of TMP-resistant bacteria. Previous laboratory evolution experiments have identified resistance-conferring mutations in the gene encoding the TMP target, bacterial dihydrofolate reductase (DHFR), in particular mutation L28R. Here, we show that 4'-desmethyltrimethoprim (4'-DTMP) inhibits both DHFR and its L28R variant, and selects against the emergence of TMP-resistant bacteria that carry the L28R mutation in laboratory experiments. Furthermore, antibiotic-sensitive E. coli populations acquire antibiotic resistance at a substantially slower rate when grown in the presence of 4'-DTMP than in the presence of TMP. We find that 4'-DTMP impedes evolution of resistance by selecting against resistant genotypes with the L28R mutation and diverting genetic trajectories to other resistance-conferring DHFR mutations with catalytic deficiencies. Our results demonstrate how a detailed characterization of resistance-conferring mutations in a target enzyme can help identify potential drugs against antibiotic-resistant bacteria, which may ultimately increase long-term efficacy of antimicrobial therapies by modulating evolutionary trajectories that lead to resistance.

How orientational order governs collectivity of folded proteins.

The past decade has witnessed the development and success of coarse-grained network models of proteins for predicting many equilibrium properties related to collective modes of motion. Curiously, the results are usually robust toward the different cutoff distances used for constructing the residue networks from the knowledge of the experimental coordinates. In this study, we present a systematical study of network construction and their effect on the predicted properties. Probing bond orientational order around each residue, we propose a natural partitioning of the interactions into an essential and a residual set. In this picture, the robustness originates from the way with which new contacts are added, so that an unusual local orientational order builds up. These residual interactions have a vanishingly small effect on the force vectors on each residue. The stability of the overall force balance then translates into the Hessian as small shifts in the slow modes of motion and an invariance of the corresponding eigenvectors. We introduce a rescaled version of the Hessian matrix and point out a link between the matrix Frobenius norm based on spectral stability arguments and orientational local order. A recipe for the optimal choice of partitioning the interactions into essential and residual components is prescribed. Implications for the study of biologically relevant properties of proteins are discussed with specific examples.

Perturbation-response scanning reveals ligand entry-exit mechanisms of ferric binding protein.

We study apo and holo forms of the bacterial ferric binding protein (FBP) which exhibits the so-called ferric transport dilemma: it uptakes iron from the host with remarkable affinity, yet releases it with ease in the cytoplasm for subsequent use. The observations fit the "conformational selection" model whereby the existence of a weakly populated, higher energy conformation that is stabilized in the presence of the ligand is proposed. We introduce a new tool that we term perturbation-response scanning (PRS) for the analysis of remote control strategies utilized. The approach relies on the systematic use of computational perturbation/response techniques based on linear response theory, by sequentially applying directed forces on single-residues along the chain and recording the resulting relative changes in the residue coordinates. We further obtain closed-form expressions for the magnitude and the directionality of the response. Using PRS, we study the ligand release mechanisms of FBP and support the findings by molecular dynamics simulations. We find that the residue-by-residue displacements between the apo and the holo forms, as determined from the X-ray structures, are faithfully reproduced by perturbations applied on the majority of the residues of the apo form. However, once the stabilizing ligand (Fe) is integrated to the system in holo FBP, perturbing only a few select residues successfully reproduces the experimental displacements. Thus, iron uptake by FBP is a favored process in the fluctuating environment of the protein, whereas iron release is controlled by mechanisms including chelation and allostery. The directional analysis that we implement in the PRS methodology implicates the latter mechanism by leading to a few distant, charged, and exposed loop residues. Upon perturbing these, irrespective of the direction of the operating forces, we find that the cap residues involved in iron release are made to operate coherently, facilitating release of the ion.