NMR characterization of the metallo-beta-lactamase from Bacteroides fragilis and its interaction with a tight-binding inhibitor: role of an active-site loop.

Understanding the structure and dynamics of the enzymes that mediate antibiotic resistance of pathogenic bacteria will allow us to take steps to combat this increasingly serious public health hazard. Complete backbone NMR resonance assignments have been made for the broad-specificity metallo-beta-lactamase CcrA from Bacteroides fragilis in the presence and absence of a tight-binding inhibitor. Chemical shift indices show that the secondary structure of the CcrA molecule in solution is very similar to that in published crystal structures. A loop adjacent to the two-zinc catalytic site exhibits significant structural variation in the published structures, but appears from the NMR experiments to be a regular beta-hairpin. Backbone heteronuclear NOE measurements indicate that this region has slightly greater flexibility on a picosecond to nanosecond time scale than the molecule as a whole. The loop appears to have an important role in the binding of substrates and inhibitors. Binding of the inhibitor 3-[2'-(S)-benzyl-3'-mercaptopropanoyl]-4-(S)-carboxy-5, 5-dimethylthiazolidine causes a marked increase in the stability of the protein toward unfolding and aggregation, and causes changes in the NMR resonance frequencies of residues close to the active (zinc-binding) site, including the beta-hairpin loop. There is a small but significant increase in the heteronuclear NOE for this loop upon inhibitor binding, indicative of a decrease in flexibility. In particular, the NOE of the indole ring of tryptophan 49, at the tip of the beta-hairpin loop, changes from a low value characteristic of a random coil chain to a significantly higher value, close to that observed for the backbone amides in this region of the protein. These results strongly suggest that the hairpin loop participates in the binding of substrate and in the shielding of the zinc sites from solvent. The broad specificity of the CcrA metallo-beta-lactamase may in fact reside in the plasticity of this part of the protein, which allows it to accommodate and bind tightly to substrates of a variety of shapes and sizes.     

Effect of cofactor binding and loop conformation on side chain methyl dynamics in dihydrofolate reductase

Dihydrofolate reductase (DHFR) has several flexible active site loops that facilitate ligand binding and catalysis. Previous studies of backbone dynamics in several complexes of DHFR indicate that the time scale and amplitude of motion depend on the conformation of the active site loops. In this study, information on dynamics is extended to methyl-containing side chains. To understand the role of side chain dynamics in ligand binding and loop conformation, methyl deuterium relaxation rates of Escherichia coli DHFR in binary folate and ternary folate:NADP+ complexes have been measured, together with chi(1) rotamer populations for threonine, isoleucine, and valine residues, determined from measurements of 3J(CgammaCO) and 3J(CgammaN) coupling constants. The results indicate that, in addition to backbone motional restriction in the adenosine-binding site, side chain flexibility in the active site and the surrounding active site loops is diminished upon binding NADP+. Resonances for several methyls in the active site and the surrounding active site loops were severely broadened in the folate:NADP+ ternary complex, suggesting the presence of motion on the chemical shift time scale. The side chains of Ile14 and Ile94, which pack against the nicotinamide and pterin rings of the cofactor and substrate, respectively, exhibit rotamer disorder in the ternary folate:NADP+ complex. Conformational fluctuations of these side chains may play a role in transition state stabilization; the observed line broadening for Ile14 suggests motions on a microsecond/millisecond time scale.     

Role of a solvent-exposed tryptophan in the recognition and binding of antibiotic substrates for a metallo-beta-lactamase

Numerous X-ray crystal structures of the metallo-beta-lactamase from Bacteroides fragilis and related organisms show a beta-hairpin loop immediately adjacent to the active-site zinc atom(s). Both crystallographic and NMR information show that the end of this beta-hairpin loop, which contains a solvent exposed tryptophan residue, Trp49, is highly flexible in the absence of substrates or other ligands, giving rise in some of the X-ray structures to a lack of observable electron density in this region. We report an investigation of the role of this mobile, solvent-exposed tryptophan using site-directed mutagenesis, steady state kinetics measurements and characterization by NMR. Trp49 appears to have a role both in substrate binding and in promotion of catalysis. Substitution of this residue with a number of different side chains indicates that the binding interaction depends on the bulky hydrophobic and aromatic nature of the indole ring, which can provide relatively non-specific interactions with a variety of antibiotic substrates. In this way, the tryptophan at this position provides a large degree of the breadth of substrate specificity for the metallo-beta-lactamase. Previous studies established that the antibiotic binding site was sufficiently plastic that the derivatization of existing antibiotics is unlikely to result in the successful treatment of bacterial infections incorporating this resistance element. Rather, a more productive approach may be to design therapeutics directed towards this solvent-exposed tryptophan residue.     

Gating of the active site of triose phosphate isomerase: Brownian dynamics simulations of flexible peptide loops in the enzyme.

The enzyme triose phosphate isomerase has flexible peptide loops at its active sites. The loops close over these sites upon substrate binding, suggesting that the dynamics of the loops could be of mechanistic and kinetic importance. To investigate these issues, the loop motions in the dimeric enzyme were simulated by Brownian dynamics. The two loops, one on each monomer, were represented by linear chains of appropriately parameterized spheres, each sphere corresponding to an amino acid residue. The loops moved in the electrostatic field of the rest of the enzyme, which was held rigid in its crystallographically observed conformation. In the absence of substrate, the loops exhibited gating of the active site with a period of about 1 ns and occupied "closed" conformations for about half of the time. As the period of gating is much shorter than the enzyme-substrate relaxation time, the motion of the loops does not reduce the rate constant for the approach of substrate from its simple diffusion-controlled value. This suggests that the flexible loops may have evolved to create the appropriate environment for catalysis while, at the same time, minimizing the kinetic penalty for gating the active site.     

Backbone dynamics in dihydrofolate reductase complexes: role of loop flexibility in the catalytic mechanism

To elucidate the influence of local motion of the polypeptide chain on the catalytic mechanism of an enzyme, we have measured (15)N relaxation data for Escherichia coli dihydrofolate reductase in three different complexes, representing different stages in the catalytic cycle of the enzyme. NMR relaxation data were analyzed by the model-free approach, corrected for rotational anisotropy, to provide insights into the backbone dynamics. There are significant differences in the backbone dynamics in the different complexes. Complexes in which the cofactor binding site is occluded by the Met20 loop display large amplitude motions on the picosecond/nanosecond time scale for residues in the Met20 loop, the adjacent betaF-betaG loop and for residues 67-69 in the adenosine binding loop. Formation of the closed Met20 loop conformation in the ternary complex with folate and NADP(+), results in attenuation of the motions in the Met20 loop and the betaF-betaG loop but leads to increased flexibility in the adenosine binding loop. New fluctuations on a microsecond/millisecond time scale are observed in the closed E:folate:NADP(+) complex in regions that form hydrogen bonds between the Met20 and the betaF-betaG loops. The data provide insights into the changes in backbone dynamics during the catalytic cycle and point to an important role of the Met20 and betaF-betaG loops in controlling access to the active site. The high flexibility of these loops in the occluded conformation is expected to promote tetrahydrofolate-assisted product release and facilitate binding of the nicotinamide ring to form the Michaelis complex. The backbone fluctuations in the Met20 loop become attenuated once it closes over the active site, thereby stabilizing the nicotinamide ring in a geometry conducive to hydride transfer. Finally, the relaxation data provide evidence for long-range motional coupling between the adenosine binding loop and distant regions of the protein.     

Interactions of a novel inhibitor from an extremophilic Bacillus sp. with HIV-1 protease: implications for the mechanism of inactivation

The active site cleft of the HIV-1 protease (PR) is bound by two identical conformationally mobile loops known as flaps, which are important for substrate binding and catalysis. The present article reports, for the first time, an HIV-1 PR inhibitor, ATBI, from an extremophilic Bacillus sp. The inhibitor is found to be a hydrophilic peptide with Mr of 1147, and an amino acid sequence of Ala-Gly-Lys-Lys-Asp-Asp-Asp-Asp-Pro-Pro-Glu. Sequence homology exhibited no similarity with the reported peptidic inhibitors of HIV-1 PR. Investigation of the kinetics of the enzyme-inhibitor interactions revealed that ATBI is a noncompetitive and tight binding inhibitor with the IC(50) and K(i) values 18.0 and 17.8 nm, respectively. The binding of the inhibitor with the enzyme and the subsequent induction of the localized conformational changes in the flap region of the HIV-1 PR were monitored by exploiting the intrinsic fluorescence of the surface exposed Trp-42 residues, which are present at the proximity of the flaps. We have demonstrated by fluorescence and circular dichroism studies that ATBI binds in the active site of the HIV-1 PR and thereby leads to the inactivation of the enzyme. Based on our results, we propose that the inactivation is due to the reorganization of the flaps impairing its flexibility leading toward inaccessibility of the substrate to the active site of the enzyme.     

Nanosecond dynamics of acetylcholinesterase near the active center gorge

To delineate the role of peptide backbone flexibility and rapid molecular motion in acetylcholinesterase catalysis and inhibitor association, we investigated the decay of fluorescence anisotropy at three sites of fluorescein conjugation to cysteine-substitution mutants of the enzyme. One cysteine was placed in a loop at the peripheral site near the rim of the active center gorge (H287C); a second was in a helical region outside of the active center gorge (T249C); a third was at the tip of a small, flexible omega loop well separated from the gorge (A262C). Mutation and fluorophore conjugation did not appreciably alter catalytic or inhibitor binding parameters of the enzyme. The results show that each site examined was associated with a high degree of segmental motion; however, the A262C and H287C sites were significantly more flexible than the T249C site. Association of the active center inhibitor, tacrine, and the peripheral site peptide inhibitor, fasciculin, had no effect on the anisotropy decay of fluorophores at positions 249 and 262. Fasciculin, but not tacrine, on the other hand, dramatically altered the decay profile of the fluorophore at the 287 position, in a manner consistent with fasciculin reducing the segmental motion of the peptide chain in this local region. The results suggest that the motions of residues near the active center gorge and across from the Cys(69)-Cys(96) omega loop are uncoupled and that ligand binding at the active center or the peripheral site does not influence acetylcholinesterase conformational dynamics globally, but induces primarily domain localized decreases in flexibility proximal to the bound ligand.     

Dynamics of the substrate binding pocket in the presence of an inhibitor covalently attached to a fungal lipase

To gain insight into the mobility of the occupied ligand-binding pocket of the Rhizomucor miehei lipase we have conducted a rigorous molecular dynamics analysis. The covalently attached inhibitor, ethylhexylphosphonate, was employed as a mimic of the putative tetrahedral intermediate in the esterolytic reaction. Our results show that in this lipase, ligand recognition is influenced by the flexibility of the binding pocket, a feature that is common to many other enzymes. Several regions around the active site were found to move significantly to adapt to the inhibitor. These motions are correlated to the flexibility of the inhibitor. In particular, the hexyl chain of the inhibitor shows considerable mobility, and adjacent residues in the binding cleft accommodate to this flexibility. Pronounced fluctuations in the binding pocket induced by the flexibility of the inhibitor are observed in the hinge region F79-S82, the active site loop region W88-V95 and the protein regions P209-F215/H257-Y260. The flexibility in the regions F79-S82 and H257-Y260, where the shorter ethyl chain is located, indicates that additional space in this binding cleft region is available for accommodating a larger moiety. Fluctuations in the region W88-V95 and P209-F215 are due to the relatively short flexible hexyl carbon chain. This part of the binding pocket could be stiffened by the presence of a longer carbon chain. Though the inhibitor is covalently attached through the phosphonate moiety, interaction of the remainder of the molecule and the enzyme are determined by hydrophobic interactions, where the Van der Waals energies are approximately 25% lower than the electrostatic contributions.     

Molecular dynamics study of the connection between flap closing and binding of fullerene-based inhibitors of the HIV-1 protease

The complementary spatial relationship between fullerene C(60) and the hydrophobic cavity region of the human immunodeficiency virus (HIV) protease, which houses the active site of the enzyme, has led to the suggestion that fullerene-based derivatives could have potential use as effective HIV protease inhibitors. The ability of such compounds to desolvate the cavity region leads to a strong hydrophobic interaction between the C(60) moiety and residues in the cavity region. In this study, the connection between the motion of the so-called flexible flaps of the cavity and favorable binding of a fullerene-based protease inhibitor is explored using multiple-time scale molecular dynamics simulations and free energy techniques. In addition, the effect of the interaction between the C(60) moiety and the residues in the cavity region on the water content of the cavity is also investigated. Conformational free energy profiles along a suitably chosen flap opening coordinate show a considerable barrier to flap opening in the presence of the inhibitor, while no such barrier exists for the protease alone. This result is interpreted in terms of a strong hydrophobic interaction between the C(60) moiety and the flexible flaps, which cause the latter to close tightly around the inhibitor, thereby expelling water from the cavity and leading to a favorable binding interaction. This interpretation is rationalized by direct analysis of the water content in the cavity in the presence and absence of the inhibitor.     

Changes in structure and dynamics of the Fv fragment of a catalytic antibody upon binding of inhibitor

Binding of the product inhibitor p-nitrophenol to the monoclonal esterolytic antibody NPN43C9 has been investigated by performing NMR spectroscopy of the heterodimeric variable-domain fragment (Fv) of the antibody in the presence and absence of inhibitor. Structural information from changes in chemical shift upon binding has been related to the changes in local dynamics in the active site of the catalytic antibody using NMR relaxation measurements. Significant changes in the chemical shifts of the backbone resonances upon binding extend beyond the immediate vicinity of the antigen binding site into the interface between the two associated polypeptides that form the Fv heterodimer, a possible indication that the binding of ligand causes a change in the relative orientations of the component light (V(L)) and heavy (V(H)) chain polypeptides. Significant differences in backbone dynamics were observed between the free Fv and the complex with p-nitrophenol. A number of resonances, including almost all of the third hypervariable loop of the light chain (L3), were greatly broadened in the free form of the protein. Other residues in the antigen-binding site showed less broadening of resonances, but still required exchange terms (R(ex)) in the model-free dynamics analysis, consistent with motion on a slow timescale in the active site region of the free Fv. Binding of p-nitrophenol caused these resonances to sharpen, but some R(ex) terms are still required in the analysis of the backbone dynamics. We conclude that the slow timescale motions in the antigen-binding site are very different in the bound and free forms of the Fv, presumably due to the damping of large-amplitude motions by the bound inhibitor.     

Coupling between Catalytic Loop Motions and Enzyme Global Dynamics

Catalytic loop motions facilitate substrate recognition and binding in many enzymes. While these motions appear to be highly flexible, their functional significance suggests that structure-encoded preferences may play a role in selecting particular mechanisms of motions. We performed an extensive study on a set of enzymes to assess whether the collective/global dynamics, as predicted by elastic network models (ENMs), facilitates or even defines the local motions undergone by functional loops. Our dataset includes a total of 117 crystal structures for ten enzymes of different sizes and oligomerization states. Each enzyme contains a specific functional/catalytic loop (10–21 residues long) that closes over the active site during catalysis. Principal component analysis (PCA) of the available crystal structures (including apo and ligand-bound forms) for each enzyme revealed the dominant conformational changes taking place in these loops upon substrate binding. These experimentally observed loop reconfigurations are shown to be predominantly driven by energetically favored modes of motion intrinsically accessible to the enzyme in the absence of its substrate. The analysis suggests that robust global modes cooperatively defined by the overall enzyme architecture also entail local components that assist in suitable opening/closure of the catalytic loop over the active site.     

Mapping backbone dynamics in solution with site-directed spin labeling: GCN4-58 bZip free and bound to DNA

In site-directed spin labeling, a nitroxide-containing side chain is introduced at selected sites in a protein. The EPR spectrum of the labeled protein encodes information about the motion of the nitroxide on the nanosecond time scale, which has contributions from the rotary diffusion of the protein, from internal motions in the side chain, and from backbone fluctuations. In the simplest model for the motion of noninteracting (surface) side chains, the contribution from the internal motion is sequence independent, as is that from protein rotary diffusion. Hence, differences in backbone motions should be revealed by comparing the sequence-dependent motions of nitroxides at structurally homologous sites. To examine this model, nitroxide side chains were introduced, one at a time, along the GCN4-58 bZip sequence, for which NMR (15)N relaxation experiments have identified a striking gradient of backbone mobility along the DNA-binding region [Bracken et al. (1999) J. Mol. Biol. 285, 2133]. Spectral simulation techniques and a simple line width measure were used to extract dynamical parameters from the EPR spectra, and the results reveal a mobility gradient similar to that observed in NMR relaxation, indicating that side chain motions mirror backbone motions. In addition, the sequence-dependent side chain dynamics were analyzed in the DNA/protein complex, which has not been previously investigated by NMR relaxation methods. As anticipated, the backbone motions are damped in the DNA-bound state, although a gradient of motion persists with residues at the DNA-binding site being the most highly ordered, similar to those of helices on globular proteins.     

Toward an understanding of the role of dynamics on enzymatic catalysis in lactate dehydrogenase

The motions of key residues at the substrate binding site of lactate dehydrogenase (LDH) were probed on the 10 ns to 10 ms time scale using laser-induced temperature-jump relaxation spectroscopy employing both UV fluorescence and isotope-edited IR absorption spectroscopy as structural probes. The dynamics of the mobile loop, which closes over the active site and is important for catalysis and binding, were characterized by studies of the inhibitor oxamate binding to the LDH/NADH binary complex monitoring the changes in emission of bound NADH. The bound NAD-pyruvate adduct, whose pyruvate moiety likely interacts with the same residues that interact with pyruvate in its ternary complex with LDH, served as a probe for any relative motions of active site residues against the substrate. The frequencies of its C=O stretch and -COO(-) antisymmetric stretch shift substantially should any relative motion of the polar moieties at the active site (His-195, Asp-168, Arg-109, and Arg-171) occur. The dynamics associated with loop closure are observed to involve several steps with motions from 1 to 300 microms. Apart from the "melting" of a few residues on the protein's surface, no kinetics were observed on any time scale in experiments of the bound NAD-pyr adduct although the measurements were made with a high degree of accuracy, even for final temperatures close to the unfolding transition of the protein. This is contrary to simple physical considerations and models. These results show that, once a productive protein/substrate complex is formed, the binding pocket is very rigid with very little, if any, motion apart from the mobile loop. The results also show that loop opening involves concomitant movement of the substrate out of the binding pocket.     

Structure and dynamics of the G121V dihydrofolate reductase mutant: lessons from a transition-state inhibitor complex

It is well known that enzyme flexibility is critical for function. This is due to the observation that the rates of intramolecular enzyme motions are often matched to the rates of intermolecular events such as substrate binding and product release. Beyond this role in progression through the reaction cycle, it has been suggested that enzyme dynamics may also promote the chemical step itself. Dihydrofolate reductase (DHFR) is a model enzyme for which dynamics have been proposed to aid in both substrate flux and catalysis. The G121V mutant of DHFR is a well studied form that exhibits a severe reduction in the rate of hydride transfer yet there remains dispute as to whether this defect is caused by altered structure, dynamics, or both. Here we address this by presenting an NMR study of the G121V mutant bound to reduced cofactor and the transition state inhibitor, methotrexate. NMR chemical shift markers demonstrate that this form predominantly adopts the closed conformation thereby allowing us to provide the first glimpse into the dynamics of a catalytically relevant complex. Based on (15)N and (2)H NMR spin relaxation, we find that the mutant complex has modest changes in ps-ns flexibility with most affected residues residing in the distal adenosine binding domain rather than the active site. Thus, aberrant ps-ns dynamics are likely not the main contributor to the decreased catalytic rate. The most dramatic effect of the mutation involves changes in μs-ms dynamics of the F-G and Met20 loops. Whereas loop motion is quenched in the wild type transition state inhibitor complex, the F-G and Met20 loops undergo excursions from the closed conformation in the mutant complex. These excursions serve to decrease the population of conformers having the correct active site configuration, thus providing an explanation for the G121V catalytic defect.     

On the mechanism of the metallo-beta-lactamase from Bacteroides fragilis.

The catalytic mechanism of metallo-beta-lactamase from Bacteroides fragilis, a dinuclear Zn(II)-containing enzyme responsible for multiple antibiotic resistance, has been investigated by using nitrocefin as a substrate. Rapid-scanning and single-wavelength stopped-flow studies revealed the accumulation during turnover of an enzyme-bound intermediate with intense absorbance at 665 nm (epsilon = 30 000 M(-1) cm(-1)). The proposed minimum kinetic mechanism for the B. fragilis metallo-beta-lactamase-catalyzed nitrocefin hydrolysis [Wang, Z., and Benkovic, S. J. (1998) J. Biol. Chem. 273, 22402-22408] was confirmed, and more accurate kinetic parameters were obtained from computer simulations and fitting. The intermediate was shown to be a novel anionic species bound to the enzyme through a Zn-acyl linkage and contains a negatively charged nitrogen leaving group. This is the first time such an intermediate was observed in the catalytic cycle of a Zn(II)-containing hydrolase and is evidence for a unique beta-lactam hydrolysis mechanism in which the amine can leave as an anion; prior protonation is not required. The electrostatic interaction between the negatively charged intermediate and the positively charged dinuclear Zn(II) center of the enzyme is important for stabilization of the intermediate. The catalytic reaction was accelerated in the presence of exogenous nucleophiles or anions, and neither the product nor the enzyme was modified during turnover, indicating that a Zn-bound hydroxide (rather than Asp-103) is the active site nucleophile. On the basis of all the information on hand, a catalytic mechanism of the B. fragilis metallo-beta-lactamase is proposed.     

An essential histidine residue for the activity of UDPglucose 4-epimerase from Kluyveromyces fragilis.

UDPglucose 4-epimerase from Kluyveromyces fragilis was completely inactivated by diethylpyrocarbonate following pseudo-first order reaction kinetics. The pH profile of diethylpyrocarbonate inhibition and reversal of inhibition by hydroxylamine suggested specific modification of histidyl residues. Statistical analysis of the residual enzyme activity and the extent of modification indicated modification of 1 essential histidine residue to be responsible for loss in catalytic activity of yeast epimerase. No major structural change in the quarternary structure was observed in the modified enzyme as shown by the identical elution pattern on a calibrated Sephacryl 200 column and association of coenzyme NAD to the apoenzyme. Failure of the substrates to afford any protection against diethylpyrocarbonate inactivation indicated the absence of the essential histidyl residue at the substrate binding region of the active site. Unlike the case of native enzyme, sodium borohydride failed to reduce the pyridine moiety of the coenzyme in the diethylpyrocarbonate-modified enzyme. This indicated the presence of the essential histidyl residue in close proximity to the coenzyme binding region of the active site. The abolition of energy transfer phenomenon between the tryptophan and coenzyme fluorophore on complete inactivation by diethylpyrocarbonate without any loss of protein or coenzyme fluorescence are also added evidences in this direction.     

Coarse-grained dynamics of the receiver domain of NtrC: fluctuations, correlations and implications for allosteric cooperativity

Receiver domains are key molecular switches in bacterial signaling. Structural studies have shown that the receiver domain of the nitrogen regulatory protein C (NtrC) exists in a conformational equilibrium encompassing both inactive and active states, with phosphorylation of Asp54 allosterically shifting the equilibrium towards the active state. To analyze dynamical fluctuations and correlations in NtrC as it undergoes activation, we have applied a coarse-grained dynamics algorithm using elastic network models. Normal mode analysis reveals possible dynamical pathways for the transition of NtrC from the inactive state to the active state. The diagonalized correlation between the inactive and the active (phosphorylated) state shows that most correlated motions occur around the active site of Asp54 and in the region Thr82 to Tyr101. This indicates a coupled correlation of dynamics in the "Thr82-Tyr101" motion. With phosphorylation inducing significant flexibility changes around the active site and alpha3 and alpha4 helices, we find that this activation makes the active-site region and the loops of alpha3/beta4 and alpha4/beta5 more stable. This means that phosphorylation entropically favors the receiver domain in its active state, and the induced conformational changes occur in an allosteric manner. Analyses of the local flexibility and long-range correlated motion also suggest a dynamics criterion for determining the allosteric cooperativity of NtrC, and may be applicable to other proteins.     

Structural basis for ligand binding and processivity in cellobiohydrolase Cel6A from Humicola insolens

The enzymatic digestion of cellulose entails intimate involvement of cellobiohydrolases, whose characteristic active-center tunnel contributes to a processive degradation of the polysaccharide. The cellobiohydrolase Cel6A displays an active site within a tunnel formed by two extended loops, which are known to open and close in response to ligand binding. Here we present five structures of wild-type and mutant forms of Cel6A from Humicola insolens in complex with nonhydrolyzable thio-oligosaccharides, at resolutions from 1.7-1.1 A, dissecting the structural accommodation of a processing substrate chain through the active center during hydrolysis. Movement of ligand is facilitated by extensive solvent-mediated interactions and through flexibility in the hydrophobic surfaces provided by a sheath of tryptophan residues.     

Dynamics of the active site loops in catalyzing aminoacylation reaction in seryl and histidyl tRNA synthetases

Aminoacylation reaction is the first step of protein biosynthesis. The catalytic reorganization at the active site of aminoacyl tRNA synthetases (aaRSs) is driven by the loop motions. There remain lacunae of understanding concerning the catalytic loop dynamics in aaRSs. We analyzed the functional loop dynamics in seryl tRNA synthetase from Methanopyrus kandleri (^mkSerRS) and histidyl tRNA synthetases from Thermus thermophilus (^ttHisRS), respectively, using molecular dynamics. Results confirm that the motif 2 loop and other active site loops are flexible spots within the catalytic domain. Catalytic residues of the loops form a network of interaction with the substrates to form a reactive state. The loops undergo transitions between closed state and open state and the relaxation of the constituent residues occurs in femtosecond to nanosecond time scale. Order parameters are higher for constituent catalytic residues which form a specific network of interaction with the substrates to form a reactive state compared to the Gly residues within the loop. The development of interaction is supported from mutation studies where the catalytic domain with mutated loop exhibits unfavorable binding energy with the substrates. During the open-close motion of the loops, the catalytic residues make relaxation by ultrafast librational motion as well as fast diffusive motion and subsequently relax rather slowly via slower diffusive motion. The Gly residues act as a hinge to facilitate the loop closing and opening by their faster relaxation behavior. The role of bound water is analyzed by comparing implicit solvent-based and explicit solvent-based simulations. Loops fail to form catalytically competent geometry in absence of water. The present result, for the first time reveals the nature of the active site loop dynamics in aaRS and their influence on catalysis.     

Molecular dynamics studies of alanine racemase: a structural model for drug design

Alanine racemase (AlaR) is a bacterial enzyme that catalyzes the interconversion of L- and D-alanine, which is an essential constituent of the peptidoglycan layer of the bacterial cell wall and requires pyridoxal 5'-phosphate (PLP) as a cofactor. The enzyme is universal to bacteria, including mycobacteria, making it an attractive target for drug design. To investigate the effects of flexibility on the binding modes of the substrate and an inhibitor and to analyze how the active site is affected by the presence of the substrate versus inhibitor, a molecular dynamics simulation on the full AlaR dimer from Bacillus stearothermophilus (pdb code: 1SFT) with a D-alanine molecule in one active site and the noncovalent inhibitor, propionate, in the second site has been carried out. Within the time scale of the simulation, we show that the active site becomes more stabilized in the presence of substrate versus inhibitor. The results of this simulation are in agreement with the proposed mechanism of alanine racemase reaction in which the substrate carboxyl group directly participates in the catalysis by acting cooperatively with Tyr 265' and Lys 39. A structural water molecule in contact with both substrate and inhibitor (i.e., in both active sites) and bridging residues in both active sites was identified. It shows a remarkably low mobility and does not exchange with bulk water. This water molecule can be taken into account for the design of specific AlaR inhibitors by either utilizing it as a bridging group or displacing it with an inhibitor atom. The results presented here provide insights into the dynamics of the alanine racemase in the presence of substrate/inhibitor, which will be used for the rational design of novel inhibitors.