Retention of Conformational Entropy upon Calmodulin Binding to Target Peptides Is Driven by Transient Salt Bridges

Calmodulin (CaM) is a highly flexible calcium-binding protein that mediates signal transduction through an ability to differentially bind to highly variable binding sequences in target proteins. To identify how binding affects CaM motions, and its relationship to conformational entropy and target peptide sequence, we have employed fully atomistic, explicit solvent molecular dynamics simulations of unbound CaM and CaM bound to five different target peptides. The calculated CaM conformational binding entropies correlate with experimentally derived conformational entropies with a correlation coefficient R² of 0.95. Selected side-chain interactions with target peptides restrain interhelical loop motions, acting to tune the conformational entropy of the bound complex via widely distributed CaM motions. In the complex with the most conformational entropy retention (CaM in complex with the neuronal nitric oxide synthase binding sequence), Lys-148 at the C-terminus of CaM forms transient salt bridges alternating between Glu side chains in the N-domain, the central linker, and the binding target. Additional analyses of CaM structures, fluctuations, and CaM-target interactions illuminate the interplay between electrostatic, side chain, and backbone properties in the ability of CaM to recognize and discriminate against targets by tuning its conformational entropy, and suggest a need to consider conformational dynamics in optimizing binding affinities.     

Photoactivation Reduces Side-Chain Dynamics of a LOV Photoreceptor

We used neutron-scattering experiments to probe the conformational dynamics of the light, oxygen, voltage (LOV) photoreceptor PpSB1-LOV from Pseudomonas putida in both the dark and light states. Global protein diffusion and internal macromolecular dynamics were measured using incoherent neutron time-of-flight and backscattering spectroscopy on the picosecond to nanosecond timescales. Global protein diffusion of PpSB1-LOV is not influenced by photoactivation. Observation-time-dependent global diffusion coefficients were found, which converge on the nanosecond timescale toward diffusion coefficients determined by dynamic light scattering. Mean-square displacements of localized internal motions and effective force constants, <k′>, describing the resilience of the proteins were determined on the respective timescales. Photoactivation significantly modifies the flexibility and the resilience of PpSB1-LOV. On the fast, picosecond timescale, small changes in the mean-square displacement and <k′> are observed, which are enhanced on the slower, nanosecond timescale. Photoactivation results in a slightly larger resilience of the photoreceptor on the fast, picosecond timescale, whereas in the nanosecond range, a significantly less resilient structure of the light-state protein is observed. For a residue-resolved interpretation of the experimental neutron-scattering data, we analyzed molecular dynamics simulations of the PpSB1-LOV X-ray structure. Based on these data, it is tempting to speculate that light-induced changes in the protein result in altered side-chain mobility mostly for residues on the protruding Jα helix and on the LOV-LOV dimer interface. Our results provide strong experimental evidence that side-chain dynamics play a crucial role in photoactivation and signaling of PpSB1-LOV via modulation of conformational entropy.     

Backbone dynamics of Tet repressor alpha8intersectionalpha9 loop

A set of single Trp mutants of class B Tet repressor (TetR), in which Trp residues are located from positions 159 to 167, has been engineered to investigate the dynamics of the loop joining the alpha-helices 8 and 9. The fluorescence anisotropy decay of most mutants can be described by the sum of three exponential components. The longest rotational correlation time, 30 ns at 10 degrees C, corresponds to the overall rotation of the protein. The shortest two components, on the subnanosecond and nanosecond time scale, are related to internal motions of the protein. The initial anisotropy, in the 0.16-0.22 range, indicates the existence of an additional ultrafast motion on the picosecond time scale. Examination of physical models for underlying motions indicates that librational motions of the Trp side chain within the rotameric chi(1) x chi(2) potential wells contribute to the picosecond depolarization process, whereas the subnanosecond and nanosecond depolarization processes are related to backbone dynamics. In the absence of inducer, the order parameters of these motions, about 0.90 and 0.80 for most positions, indicate limited flexibility of the loop backbone. Anhydrotetracycline binding to TetR induces an increased mobility of the loop on the nanosecond time scale. This suggests that entropic factors might play a role in the mechanism of allosteric transition.     

Concordance of residual dipolar couplings, backbone order parameters and crystallographic B-factors for a small alpha/beta protein: a unified picture of high probability, fast atomic motions in proteins

Using ensemble refinement of the third immunoglobulin binding domain (GB3) of streptococcal protein G (a small alpha/beta protein of 56 residues), we demonstrate that backbone (N-H, N-C', Calpha-Halpha, Calpha-C') residual dipolar coupling data in five independent alignment media, generalized order parameters from 15N relaxation data, and B-factors from a high-resolution (1.1A), room temperature crystal structure are entirely consistent with one another within experimental error. The optimal ensemble size representation is between four and eight, as assessed by complete cross-validation of the residual dipolar couplings. Thus, in the case of GB3, all three observables reflect the same low-amplitude anisotropic motions arising from fluctuations in backbone phi/psi torsion angles in the picosecond to nanosecond regime in both solution and crystalline environments, yielding a unified picture of fast, high-probability atomic motions in proteins. An understanding of these motions is crucial for understanding the impact of protein dynamics on protein function, since they provide part of the driving force for triggered conformational changes that occur, for example, upon ligand binding, signal transduction and enzyme catalysis.     

Side-chain flexibility in protein-ligand binding: the minimal rotation hypothesis

The goal of this work is to learn from nature about the magnitudes of side-chain motions that occur when proteins bind small organic molecules, and model these motions to improve the prediction of protein-ligand complexes. Following analysis of protein side-chain motions upon ligand binding in 63 complexes, we tested the ability of the docking tool SLIDE to model these motions without being restricted to rotameric transitions or deciding which side chains should be considered as flexible. The model tested is that side-chain conformational changes involving more atoms or larger rotations are likely to be more costly and less prevalent than small motions due to energy barriers between rotamers and the potential of large motions to cause new steric clashes. Accordingly, SLIDE adjusts the protein and ligand side groups as little as necessary to achieve steric complementarity. We tested the hypothesis that small motions are sufficient to achieve good dockings using 63 ligands and the apo structures of 20 different proteins and compared SLIDE side-chain rotations to those experimentally observed. None of these proteins undergoes major main-chain conformational change upon ligand binding, ensuring that side-chain flexibility modeling is not required to compensate for main-chain motions. Although more frugal in the number of side-chain rotations performed, this model substantially mimics the experimentally observed motions. Most side chains do not shift to a new rotamer, and small motions are both necessary and sufficient to predict the correct binding orientation and most protein-ligand interactions for the 20 proteins analyzed.     

Ester-linked hen egg white lysozyme shows a compact fold in a molecular dynamics simulation - possible causes and sensitivity of experimentally observable quantities to structural changes maintaining this compact fold

Prediction and understanding of the folding and stability of the 3D structure of proteins is still a challenge. The different atomic interactions, such as non polar contacts and hydrogen bonding, are known but their exact relative weights and roles when contributing to protein folding and stability are not identified. Initiated by a previous molecular dynamics simulation of fully ester-linked hen egg white lysozyme (HEWL), which showed a more compact fold of the ester-linked molecule compared to the native one, three variants of this protein are analyzed in the present study. These are 129-residue native HEWL, partly ester-linked HEWL, in which only 34 peptide linkages that are not involved in the helical or β-strand parts of native HEWL were replaced by ester linkages, and fully (126 residues) ester-linked HEWL. Native and partly ester-linked HEWL showed comparable behaviour, whereas fully ester-linked HEWL could not maintain the native secondary structure of HEWL in the simulation and adopted a more compact fold. The conformational changes were analyzed by comparing simulation averaged values of quantities that can be measured by NMR, such as (1)H-(15)N backbone order parameters, residual dipolar couplings, proton-proton NOE distances and (3)J-couplings with the corresponding values derived from experimental NMR data for native HEWL. The information content of the latter appeared to be insufficient to detect the local conformational rearrangements upon esterification of the loop regions of the protein. For fully ester-linked HEWL, a significantly reduced agreement was observed. Upon esterification, the backbone-side chain and side chain-side chain hydrogen-bonding pattern of HEWL changes to maintain its compactness and thus the structural stability of the ester-linked lysozymes.     

Main chain and side chain dynamics of the ubiquitin conjugating enzyme variant human Mms2 in the free and ubiquitin-bound States

Protein ubiquitination involves a cascade of enzymatic steps where ubiquitin (Ub) is sequentially transferred as a thiolester intermediate from an E1 enzyme to an E2 enzyme and finally to the protein target with the help of a Ub-protein ligase. Protein ubiquitination brought about by the Ubc13-Mms2 (E2-E2) complex has a unique role in the cell, unrelated to protein degradation. The Mms2-Ubc13 heterodimer links Ub molecules to one another through an isopeptide bond between its own C-terminus and Lys-63 on another Ub. The role of Mms2 is to orient a target-bound Ub molecule such that its Lys-63 is proximal to the C-terminus of the Ub molecule that is covalently linked to the active site of Ubc13. To gain insight into the influence of protein dynamics on the affinity of Ub for Mms2, we have determined pico- to nanosecond time scale fluctuations of the main chain and methyl side chains of human Mms2 in the free and Ub-bound states using solution state (15)N and (2)H nuclear magnetic resonance relaxation measurements. Analysis of the relaxation data allows for a semiquantitative estimation of the conformational entropy change (TDeltaS(NMR)) for the main chain and side chain methyl groups of Mms2 upon binding Ub. The value of TDeltaS(NMR) for the main chain and side chain methyl groups of Mms2 is -8 +/- 2 and -2 +/- 2 kcal mol(-)(1), respectively. The experimental DeltaG(binding) for the Mms2.Ub complex is -6 kcal mol(-)(1). Estimation of DeltaG(binding) using an empirical structure-based approach that does not account for changes in main chain entropy yields a value of -17 +/- 2 kcal mol(-)(1). However, inclusion of TDeltaS(NMR) for the main chain of Mms2 increases the estimated DeltaG(binding) to -9 +/- 3 kcal mol(-)(1). Assuming that changes in Ub main chain dynamics contribute to TDeltaS(NMR) to the same extent as Mms2, the estimated DeltaG(binding) is further reduced to -1 +/- 4 kcal mol(-)(1), a value close to the experimental DeltaG(binding).     

The role of slow and fast protein motions in allosteric interactions

Allostery is fundamentally thermodynamic in nature. Long-range communication in proteins may be mediated not only by changes in the mean conformation with enthalpic contribution but also by changes in dynamic fluctuations with entropic contribution. The important role of protein motions in mediating allosteric interactions has been established by NMR spectroscopy. By using CAP as a model system, we have shown how changes in protein structure and internal dynamics can allosterically regulate protein function and activity. The results indicate that changes in conformational entropy can give rise to binding enhancement, binding inhibition, or have no effect in the expected affinity, depending on the magnitude and sign of enthalpy–entropy compensation. Moreover, allosteric interactions can be regulated by the modulation a low-populated conformation states that serve as on-pathway intermediates for ligand binding. Taken together, the interplay between fast internal motions, which are intimately related to conformational entropy, and slow internal motions, which are related to poorly populated conformational states, can regulate protein activity in a way that cannot be predicted on the basis of the protein’s ground-state structure.     

Dynamic coupling and allosteric behavior in a nonallosteric protein

Long-range intraprotein interactions give rise to many important protein behaviors. Understanding how energy is transduced through protein structures to either transmit a signal or elicit conformational changes is therefore a current challenge in structural biology. In an effort to understand such linkages, multiple V --> A mutations were made in the small globular protein eglin c. The physical responses, as mapped by NMR spin relaxation, residual dipolar couplings (RDCs), and scalar couplings, illustrate that the interior of this nonallosteric protein forms a dynamic network and that local perturbations are transmitted as dynamic and structural changes to distal sites as far as 16 A away. Two basic types of propagation responses were observed: contiguous pathways of enhanced (attenuated) dynamics with no change in structure; and dispersed (noncontiguous) changes in methyl rotation rates that appear to result from subtle deformation of backbone structure. In addition, energy transmission is found to be unidirectional. In one mutant, an allosteric conformational change of a side chain is seen in the context of a pathway of propagated changes in picosecond to nanosecond dynamics. The observation of so many long-range interactions in a small, rigid system lends experimental weight to the idea that all well-folded proteins inherently possess allosteric features [Gunasekaran et al. (2004) Proteins 57, 433-443], and that dynamics are a rich source of information for mapping and gaining mechanistic insight into communication pathways in individual proteins.     

Conformational entropy changes upon lactose binding to the carbohydrate recognition domain of galectin-3

The conformational entropy of proteins can make significant contributions to the free energy of ligand binding. NMR spin relaxation enables site-specific investigation of conformational entropy, via order parameters that parameterize local reorientational fluctuations of rank-2 tensors. Here we have probed the conformational entropy of lactose binding to the carbohydrate recognition domain of galectin-3 (Gal3), a protein that plays an important role in cell growth, cell differentiation, cell cycle regulation, and apoptosis, making it a potential target for therapeutic intervention in inflammation and cancer. We used ¹⁵N spin relaxation experiments and molecular dynamics simulations to monitor the backbone amides and secondary amines of the tryptophan and arginine side chains in the ligand-free and lactose-bound states of Gal3. Overall, we observe good agreement between the experimental and computed order parameters of the ligand-free and lactose-bound states. Thus, the ¹⁵N spin relaxation data indicate that the molecular dynamics simulations provide reliable information on the conformational entropy of the binding process. The molecular dynamics simulations reveal a correlation between the simulated order parameters and residue-specific backbone entropy, re-emphasizing that order parameters provide useful estimates of local conformational entropy. The present results show that the protein backbone exhibits an increase in conformational entropy upon binding lactose, without any accompanying structural changes.     

Dynamics of a de novo designed three-helix bundle protein studied by 15N, 13C, and 2H NMR relaxation methods

Understanding how the amino acid sequence of a polypeptide chain specifies a unique, functional three-dimensional structure remains an important goal, especially in the context of the emerging discipline of de novo protein design. Alpha3D is a single chain protein of 73 amino acids resulting from a de novo design effort. Previous solution nuclear magnetic resonance studies of alpha3D confirm that the protein adopts the designed structure of a three-helix bundle. Furthermore, alpha3D has been previously shown to possess all of the major thermodynamic and structural characteristics of natural proteins, though it shares no sequence homology to any protein sequence in the database. In this work, the backbone and side-chain dynamics of alpha3D were investigated using 15N, 13C, and 2H nuclear magnetic resonance relaxation methods with the aim of assessing the character of the internal motions of this native-like protein of de novo design. At the backbone level, both 15N and 13C(alpha) relaxation studies indicate highly restrictive motion on the picosecond to nanosecond time scale in the alpha-helical regions of alpha3D, with increasing mobility at the ends of the alpha-helices and in the two loop regions. This is largely consistent with what is seen in proteins of natural origin. Overall, the view provided by both 2H and 13C methyl relaxation methods suggest that the side chains of alpha3D are more dynamic compared to natural proteins. Regions of relative flexibility bound clusters of rigid methyl-bearing side-chain groups that are interspersed with aromatic and beta-branched amino acids. The time scale of motions associated with methyl-bearing side chains of alpha3D are significantly longer than that seen in natural proteins. These results indicate that the strategies underlying the design of alpha3D have largely, but not completely, captured both the structural and dynamic character of natural proteins.     

Internal motions and exchange processes in human ileal bile acid binding protein as studied by backbone (15)N nuclear magnetic resonance spectroscopy

Human ileal bile acid binding protein (I-BABP), a member of the family of intracellular lipid binding proteins, is thought to play a role in the enterohepatic circulation of bile salts. Previously, we have shown by stopped-flow fluorescence analysis that positive binding cooperativity exhibited by I-BABP in its interactions with glycocholate (GCA) and glycochenodeoxycholate (GCDA), the two primary bile salts in humans, is related to a slow conformational change in the protein. In this study, we used backbone (15)N relaxation nuclear magnetic resonance (NMR) techniques to obtain residue-specific information about the internal dynamics of apo I-BABP and the doubly ligated I-BABP:GCA:GCDA complex on various time scales. According to our NMR data, bile salt binding is accompanied by a slight rigidification of the (15)N-(1)H bond vectors on the picosecond to nanosecond time scale, with most pronounced changes occurring in the C-D region. In contrast to the minor effects of ligation on fast motions, relaxation dispersion NMR experiments indicate a marked difference between the two protein states on the microsecond to millisecond time scale. In the apo form, an extensive network of conformational fluctuations is detected throughout segments of the EFGHIJ β-strands and the C-D loop, which cease upon complexation. Our NMR data are in agreement with a conformational selection model we proposed earlier for I-BABP and support the hypothesis of an allosteric mechanism of ligand binding. According to the NMR measurements, the helical cap region may have a less crucial role in mediating ligand entry and release than what has been indicated for fatty acid binding proteins.     

Dynamically Driven Protein Allostery Exhibits Disparate Responses for Fast and Slow Motions

There is considerable interest in the dynamic aspect of allosteric action, and in a growing list of proteins allostery has been characterized as being mediated predominantly by a change in dynamics, not a transition in conformation. For considering conformational dynamics, a protein molecule can be simplified into a number of relatively rigid microdomains connected by joints, corresponding to, e.g., communities and edges from a community network analysis. Binding of an allosteric activator strengthens intermicrodomain coupling, thereby quenching fast (e.g., picosecond to nanosecond) local motions but initiating slow (e.g., microsecond to millisecond), cross-microdomain correlated motions that are potentially of functional importance. This scenario explains allosteric effects observed in many unrelated proteins.     

Analysis of side chain mobility among protein G B1 domain mutants with widely varying stabilities

"Host-guest" studies of the B1 domain from Streptococcal protein G have been used previously to establish a thermodynamic scale for the beta-sheet-forming propensities of the 20 common amino acids. To investigate the contribution of side chain conformational entropy to the relative stabilities of B1 domain mutants, we have determined the dynamics of side chain methyl groups in 10 of the 20 mutants used in a previous study. Deuterium relaxation rates were measured using two-dimensional NMR techniques for 13CH2D groups. Analysis of the relaxation data using the Lipari-Szabo model-free formalism showed that mutations introduced at the guest position caused small but statistically significant changes in the methyl group dynamics. In addition, there was a low level of covariation of the Lipari-Szabo order parameters among the 10 mutants. The variations in conformational free energy estimated from the order parameters were comparable in magnitude to the variations in global stability of the 10 mutants but did not correlate with the global stability of the domain or with the structural properties of the guest amino acids. The data support the view that conformational entropy in the folded state is one of many factors that can influence the folding thermodynamics of proteins.     

Protein dynamics and allostery: an NMR view

Allostery, the process by which distant sites within a protein system are energetically coupled, is an efficient and ubiquitous mechanism for activity regulation. A purely mechanical view of allostery invoking only structural changes has developed over the decades as the classical view of the phenomenon. However, a fast growing list of examples illustrate the intimate link between internal motions over a wide range of time scales and function in protein-ligand interactions. Proteins respond to perturbations by redistributing their motions and they use fluctuating conformational states for binding and conformational entropy as a carrier of allosteric energy to modulate association with ligands. In several cases allosteric interactions proceed with minimal or no structural changes. We discuss emerging paradigms for the central role of protein dynamics in allostery.     

Conformational plasticity of the lipid transfer protein SCP2

The nonspecific lipid transfer protein sterol carrier protein 2 (SCP2) is involved in organellar fatty acid metabolism. A hydrophobic cavity in the structure of SCP2 accommodates a wide variety of apolar ligands such as cholesterol derivatives or fatty acyl-coenzyme A (CoA) conjugates. The properties of this nonspecific lipid binding pocket are explored using NMR chemical shift perturbations, paramagnetic relaxation enhancement, amide hydrogen exchange, and 15N relaxation measurements. A common binding cavity shared by different physiological ligands is identified. NMR relaxation measurements reveal that residues in the three C-terminal alpha-helices within the lipid binding region exhibit mobility at fast (picosecond to nanosecond) and slow (microsecond to millisecond) time scales. Ligand binding is associated with a considerable loss of peptide backbone mobility. The observed conformational dynamics in SCP2 may play a role for the access of hydrophobic ligands to an occluded binding pocket. The C-terminal peroxisomal targeting signal of SCP2 is specifically recognized by the Pex5p receptor protein, which conducts cargo proteins toward the peroxisomal organelle. Neither the C-terminal targeting signal nor the N-terminal precursor sequence interferes with lipid binding by SCP2. The alpha-helices involved in lipid binding also mediate a secondary interaction interface with the Pex5p receptor. Silencing of conformational dynamics of the peptide backbone in these helices upon either lipid or Pex5p binding might communicate the loading state of the cargo protein to the targeting receptor.     

Thermodynamic interpretation of protein dynamics from NMR relaxation measurements

Protein dynamics and thermodynamics can be characterized through measurements of relaxation rates of side chain (2)H and (13)C, and backbone (15)N nuclei using NMR spectroscopy. The rates reflect protein motions on timescales from picoseconds to milliseconds. Backbone and methyl side chain NMR relaxation measurements for several proteins are beginning to reveal the role of protein dynamics in protein stability and ligand binding.     

153 Wonderful roles of the entropy in protein dynamics,binding and folding

The entropy, which is central to the second law of thermodynamics, determines that the thermal energy always flows spontaneously from regions of higher temperature to regions of lower temperature. In the protein–solvent thermodynamic system, the entropy is defined as a measure of how evenly the thermal energy would distribute over the entire system (Liu et al., 2012). Such tendency to distribute energy as evenly as possible will reduce the state of order of the initial system, and hence, the entropy can be regarded as an expression of the disorder, or randomness of the system (Yang et al., 2012). For a protein–solvent system under a constant solvent condition, the origin of entropy is the thermal energy stored in atoms, which makes atoms jostle around and bump onto one another, thus leading to vibrations of the covalent bonds connecting two atoms (occurring on the fs timescale) and the rotational and translational motions of amino acid side chain groups (occurring on ps timescale) and water molecules. These motions break the noncovalent bonds around structural regions that are weakly constrained thereby triggering the competitive interactions among residues or between residues and water molecules leading ultimately to the loop motions (occurring on ns timescale) around the protein surface. The loop motions can further transmit either through the water network around the protein surface or via specific structural components (such as the hinge-bending regions) over the entire protein molecule leading to large concerted motions (occurring on μs to s timescales) that are most relevant to protein functions (Amadei, Linssen & Berendsen, 1993; Tao, Rao & Liu, 2010). Thus, the multiple hierarchies of the protein dynamics on distinct timescales (Henzler-Wildman & Kern, 2007) are a consequence of the cascade amplification of the microscopic motions of atoms and groups for which the entropy originating from atomic thermal energy is most fundamental. In the case of protein–ligand binding, the importance of the entropy is embodied in the following aspects. (i) The release of the water molecule kinetic energy (which is a process of the solvent entropy maximization) will cause Brownian motions of individual water molecules which result in strong Brownian bombardments to solute molecules causing molecule wanders/diffusions and subsequent accident contacts/collisions between proteins and ligands. (ii) Such collisions will inevitably cause water molecule displacement and, if the contact interfaces are properly complementary, the requirement to increase the solvent entropy would further displace the water network around the binding interfaces thus leading to the formation the initial protein-ligand complex. (iii) In the initial complex, the loose association of the two partners provide the opportunity for protein to increase conformational entropy, thus triggering the conformational adjustments through competitive interaction between protein residues and ligand, leading ultimately to the formation of tightly associated complex (Liu et al., 2012). In the protein folding process, the first stage, i.e. the rapid hydrophobic collapse (Agashe, Shastry & Udgaonkar, 1995; Dill, 1985), is in fact driven by the effect of the solvent entropy maximization. Specifically, the requirement to maintain as many as possible the dynamic hydrogen bonds among the water molecules will squeeze/sequestrate the hydrophobic amino acid side chains into the interior of the folding intermediates and expose the polar/charged side chains onto the intermediate surface. This will minimize the solvent accessible surface area of the folding intermediates and as thus maximize the entropy of the solvent. The resulting molten globule states (Ohgushi & Wada, 1983) may contain a few secondary structural components and native tertiary contacts, while many native contacts, or close residue–residue interactions present in the native state have not yet formed. However, the nature to increase the protein conformational entropy can trigger a further conformational adjustment process, i.e. the conformational entropy increase breaks the transient secondary or tertiary contacts and triggers the competitive interactions among protein residues and between residues and water. This process may repeat many rounds until the negative enthalpy change resulting from the noncovalent formations can overcompensate for protein conformational entropy loss. In summary, we consider that the tendency to maximize the entropy of the protein–solvent system, which originates from the atomic thermal energy, is the most fundamental driving factor for protein folding, binding, and dynamics, whereas the enthalpy reduction, an opposing factor that tends to make the system become ordered, can compensate for the effect of entropy loss to ultimately allow the system to reach equilibrium at the free energy minima, either global or local.     

Dynamics of the metallo-beta-lactamase from Bacteroides fragilis in the presence and absence of a tight-binding inhibitor

A significant determinant for the broad substrate specificity of the metallo-beta-lactamases from Bacteroides fragilis and other similar organisms is the presence of a plastic substrate binding site that is nevertheless capable of tight substrate binding in the Michaelis complex. To achieve these two competing ends, the molecule apparently employs a flexible flap that closes over the active site in the presence of substrate. These characteristics imply that dynamic changes are an important component of the mechanism of action of these enzymes. The backbone and tryptophan side chain dynamics of the metallo-beta-lactamase from B. fragilis have been examined using (15)N NMR relaxation measurements. Two states of the protein were examined, in the presence and absence of a tight-binding inhibitor. Relaxation measurements were analyzed by the model-free method. Overall, the metallo-beta-lactamase molecule is rigid and shows little flexibility except in loops. The flexibility of the loop that covers the active site is not unusually great as compared to the other loops of the protein. Local motion on a picosecond time scale was found to be very similar throughout the protein in the presence and absence of the inhibitor, but a significant difference was observed in the motions on a nanosecond time scale (tau(e)). Large-amplitude motions with a time constant of about 1.3 ns were observed for the flexible flap region (residues 45-55) in the absence of the inhibitor. These motions were completely damped out in the presence of the inhibitor. In addition, the motion of a tryptophan side chain at the tip of the beta-hairpin of the flap shows a very significant difference in motion on the ps time scale. These results indicate that the motions of the polypeptide chain in the flap region can be invoked to explain both the wide substrate specificity (the free form has considerable amplitude of motion in this region) and the catalytic efficiency of the metallo-beta-lactamase (the motions are damped out when the inhibitor and by implication a substrate binds in the active site).