Pivotal role of Gly 121 in dihydrofolate reductase from Escherichia coli: the altered structure of a mutant enzyme may form the basis of its diminished catalytic performance

The structure and folding of dihydrofolate reductase (DHFR) from Escherichia coli and the mutant G121V-DHFR, in which glycine 121 in the exterior FG loop was replaced with valine, were studied by molecular dynamics simulations and CD and fluorescence spectroscopy. The importance of residue 121 for the chemical step during DHFR catalysis had been demonstrated previously. High-temperature MD simulations indicated that while DHFR and G121V-DHFR followed similar unfolding pathways, the strong contacts between the M20 loop and the FG loop in DHFR were less stable in the mutant. These contacts have been proposed to be involved in a coupled network of interactions that influence the protein dynamics and promote catalysis [Benkovic, S. J., and Hammes-Schiffer, S. (2003) Science 301, 1196-1202]. CD spectroscopy of DHFR and G121V-DHFR indicated that the two proteins existed in different conformations at room temperature. While the thermally induced unfolding of DHFR was highly cooperative with a midpoint at 51.6 +/- 0.7 degrees C, G121V-DHFR exhibited a gradual decrease in its level of secondary structure without a clear melting temperature. Temperature-induced unfolding and renaturation from the urea-denatured state revealed that both proteins folded via highly fluorescent intermediates. The formation of these intermediates occurred with relaxation times of 149 +/- 4.5 and 256 +/- 13 ms for DHFR and G121V-DHFR, respectively. The fluorescence intensity for the intermediates formed during refolding of G121V-DHFR was approximately twice that of the wild-type. While the fluorescence intensity then slowly decayed for DHFR toward a state representing the native protein, G121V-DHFR appeared to be trapped in a highly fluorescent state. These results suggest that the reduced catalytic activity of G121V-DHFR is the consequence of nonlocal structural effects that may result in a perturbation of the network of promoting motions.     

Modeling of Hidden Structures Using Sparse Chemical Shift Data from NMR Relaxation Dispersion

NMR relaxation dispersion measurements report on conformational changes occurring on the μs-ms timescale. Chemical shift information derived from relaxation dispersion can be used to generate structural models of weakly populated alternative conformational states. Current methods to obtain such models rely on determining the signs of chemical shift changes between the conformational states, which are difficult to obtain in many situations. Here, we use a “sample and select” method to generate relevant structural models of alternative conformations of the C-terminal-associated region of Escherichia coli dihydrofolate reductase (DHFR), using only unsigned chemical shift changes for backbone amides and carbonyls (¹H, ¹⁵N, and ¹³C′). We find that CS-Rosetta sampling with unsigned chemical shift changes generates a diversity of structures that are sufficient to characterize a minor conformational state of the C-terminal region of DHFR. The excited state differs from the ground state by a change in secondary structure, consistent with previous predictions from chemical shift hypersurfaces and validated by the x-ray structure of a partially humanized mutant of E. coli DHFR (N23PP/G51PEKN). The results demonstrate that the combination of fragment modeling with sparse chemical shift data can determine the structure of an alternative conformation of DHFR sampled on the μs-ms timescale. Such methods will be useful for characterizing alternative states, which can potentially be used for in silico drug screening, as well as contributing to understanding the role of minor states in biology and molecular evolution.     

Coupling interactions of distal residues enhance dihydrofolate reductase catalysis: mutational effects on hydride transfer rates

Recently, the participation in catalysis of residues spatially removed from the enzyme's active site has received considerable attention. The influence of the distal Gly-121 residue on the chemical step of hydride transfer in dihydrofolate reductase (DHFR) catalysis had been demonstrated previously [Cameron, C. E., and Benkovic, S. J. (1997) Biochemistry 36, 15792-15800]. In our continuing effort to identify catalytically important residues that are distal from the active site, we used sequence conservation information, kinetic data on site-directed mutants, dynamic motion information from NMR methods, and correlated motions from MD simulations to identify a subset of residues. Among them, the region spanning positions 41-45 is distal to the active site and was chosen as the focus for the mutagenesis and kinetic studies reported here. Specifically, the highly conserved Met-42 was selected for site-directed mutagenesis. While the reaction kinetics for the M42F mutant enzyme did not deviate from wild-type behavior, a 41-fold reduction in the forward hydride transfer rate was found for the M42W mutant. Given the established role of Gly-121 in the hydride transfer process, double mutant enzymes involving positions 42 and 121 were constructed and characterized. These double mutant enzymes generally showed little changes in substrate and cofactor binding but synergistic decreases in forward hydride transfer rates, while the decreases in reverse rates were additive. Along with supporting information from mixed quantum/classical MD simulations [Agarwal, P. K., Billeter, S. R., Rajagopalan, P. T., Benkovic, S. J., and Hammes-Schiffer, S. (2002) Proc. Natl. Acad. Sci. U.S.A. 99, 2794-2799], the data suggest that a coupled dynamic motion of these distal residues enhances crossing of the chemical reaction barrier and imply an expanded nonstatic role for the protein fold in catalysis.     

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.     

Site-directed deletion mutants of a carboxyl-terminal region of human dihydrofolate reductase.

Substrate and inhibitor binding to dihydrofolate reductase (DHFR) primarily involves residues in the amino-terminal half of the enzyme; however, antibody binding studies performed in this laboratory suggested that the loop region located in the carboxyl terminus of human DHFR (hDHFR; residues 140-186) is involved in conformational changes that occur upon ligand binding and affect enzyme function (Ratnam, M., Tan, X., Prendergast, N.J., Smith, P.L. & Freisheim, J.H. (1988) Biochemistry 27, 4800-4804). To investigate this observation further, site-directed mutagenesis was used to construct deletion mutants of hDHFR missing 1 (del-1), 2 (del-2), 4 (del-4), and 6 (del-6) residues from loops in the carboxyl terminus of the enzyme. The del-1 mutant enzyme has a two-amino acid substitution in addition to the one-amino acid deletion. Deletion of only one amino acid resulted in a 35% decrease in the specific activity of the enzyme. The del-6 mutant enzyme was inactive. Surprisingly, the del-4 mutant enzyme retained a specific activity almost 33% that of the wild type. The specific activity of the del-2 mutant enzyme was slightly higher (38% wild-type activity) than that of the del-4 mutant. All three active deletion mutants were much less stable than the wild-type enzyme, and all three showed at least a 10-fold increase in Km values for both substrates. The del-1 and del-2 mutants exhibited a similar increase in KD values for both substrate and cofactor. The three active deletion mutants lost activity at concentrations of activating agents such as KCl, urea, and p-hydroxymercuribenzoate that continued to stimulate the wild-type enzyme. Antibody binding studies revealed conformational differences between the wild-type and mutant enzymes both in the absence and presence of bound folate. Thus, although the loops near the carboxyl terminus are far removed from the active site, small deletions of this region significantly affect DHFR function, indicating that the loop structure in mammalian DHFR plays an important functional role in its conformation and catalysis.     

Mapping Protein-Protein Proximity in the Purinosome

The enzymes in the human de novo purine synthesis pathway were found to form a cellular complex, the purinosome, upon culturing cells in purine-depleted medium (An, S., Kumar R., Sheets, E. D., and Benkovic, S. J. (2008) Science 320, 103–106). Purinosome formation and dissociation were found to be modulated by several factors, including the microtubule network and cell signaling involving protein phosphorylation. To determine whether the pathway enzymes are in physical contact, we probed for the protein-protein interactions (PPIs) within the purinosome with a novel application of the Tango PPI reporter system (Barnea, G., Strapps, W., Herrada, G., Berman, Y., Ong, J., Kloss, B., Axel, R., and Lee, K. J. (2008) Proc. Natl. Acad. Sci. U.S.A. 105, 64–69). We found PPIs among all six enzymes within the pathway and evidence for a core involving the first three enzymes. We also captured purinosomes under both purine-rich and purine-depleted conditions. The results provide additional insights into the transient nature and topography of the purinosome.     

Kinetic mechanism of DNA polymerase I (Klenow fragment): identification of a second conformational change and evaluation of the internal equilibrium constant.

In a previously determined minimal kinetic scheme for DNA polymerization catalyzed by the Klenow fragment (KF) of Escherichia coli DNA polymerase I, a nonchemical step that interconverted the KF'.DNAn+1.PPi and KF.DNAn+1PPi complexes was not observed in correct incorporation [Kuchta, R. D., Mizrahi, V., Benkovic, P.A., Johnson, K.A., & Benkovic, S.J. (1987) Biochemistry 26, 8410-8417] but was detected in misincorporation [Kuchta, R. D., Benkovic, P.A., & Benkovic, S.J. (1988) Biochemistry 27, 6716-6725]. In a pulse-chase experiment in this study, a burst amplitude of 100% of the enzyme concentration is observed; under pulse-quench conditions, the burst amplitude is 80%, indicative of the accumulation of the KF'.DNA.dNTP species owing to a slow step subsequent to chemical bond formation. This latter step was unequivocally identified by single-turnover pyrophosphorolysis and pyrophosphate-exchange experiments as one interconverting KF'.DNAn+1.PPi and KF.DNAn+1.PPi. The rate constants for this step in both directions were established through the rate constants for processive synthesis and pyrophosphorolysis. Pyrophosphorolysis of a 3'-phosphorothioate DNA duplex confirmed that the large elemental effect observed previously [Mizrahi, V., Henrie, R. N., Marlier, J.F., Johnson, K.A., & Benkovic, S.J. (1985) Biochemistry 24, 4010-4018] in this direction but not in polymerization is due to a marked decrease in the affinity of KF for the phosphorothioate-substituted duplex and not to the chemical step. The combination of the experimentally measured equilibrium constant for the bound KF.DNA species with the collective kinetic measurements further extends previous insights into the dynamics of the polymerization process catalyzed by KF.     

Structural Model of Weak Binding Actomyosin in the Prepowerstroke State

We present the first in silico model of the weak binding actomyosin in the initial powerstroke state, representing the actin binding-induced major structural changes in myosin. First, we docked an actin trimer to prepowerstroke myosin then relaxed the complex by a 100-ns long unrestrained molecular dynamics. In the first few nanoseconds, actin binding induced an extra primed myosin state, i.e. the further priming of the myosin lever by 18° coupled to a further closure of switch 2 loop. We demonstrated that actin induces the extra primed state of myosin specifically through the actin N terminus-activation loop interaction. The applied in silico methodology was validated by forming rigor structures that perfectly fitted into an experimentally determined EM map of the rigor actomyosin. Our results unveiled the role of actin in the powerstroke by presenting that actin moves the myosin lever to the extra primed state that leads to the effective lever swing.     

Mapping the antigenic epitopes of human dihydrofolate reductase by systematic synthesis of peptides on solid supports.

All of the 181 possible overlapping hexapeptides as well as 179 octapeptides covering the amino acid sequence of human dihydrofolate reductase (hDHFR) were synthesized on polyethylene supports. The synthetic procedure of Geysen et al. (Geysen, H. M., Rodda, S. J., Mason, T. J., Tribbick, G., and Schoofs, P. G. (1987) J. Immunol. Methods 102, 259-274) was modified to obtain up to 100 nmol of peptide on each pin. Peptides constituting antigenic epitopes on hDHFR were identified by examining the binding of antibodies raised against both native and denatured hDHFR to these peptides by enzyme-linked immunosorbent assay. The peptides bound in a similar pattern to polyclonal antibodies against both native and denatured dihydrofolate reductase (DHFR). Six major epitopes were located corresponding to residues 27-33, 45-51, 67-74, 133-139, 153-158, and 176-181 using both hexapeptides and octapeptides. An additional epitope, constituting residues 14-21, was found by the use of octapeptides. Most of the epitopes are hydrophilic and reside largely in "loop" regions at the boundaries of secondary structural elements of hDHFR. This observation is consistent with our previous results which suggested that ligand binding at the active site of the enzyme can cause a dramatic reduction in antibody binding to DHFR due to conformational constraints in flexible loop regions in various parts of the molecule. The similarity of the immunogenic profiles of native versus denatured hDHFR indicates that the two forms of the antigen share the same amino acid sequence-specific epitopes. Competitive enzyme-linked immunosorbent assay showed that the binding of anti-hDHFR antiserum to both native and denatured hDHFR was inhibited by approximately 30% by the seven antigenic peptides, indicating that a significant proportion of the antibodies elicited by this enzyme is specific for short peptides. Besides revealing the antigenic structure of DHFR our results provide a rational basis for the design of mutant DHFRs to study the importance of loop residues in the conformational dynamics of the enzyme.     

Structural Elucidation of the Cyclization Mechanism of α-1,6-Glucan by Bacillus circulans T-3040 Cycloisomaltooligosaccharide Glucanotransferase

Bacillus circulans T-3040 cycloisomaltooligosaccharide glucanotransferase belongs to the glycoside hydrolase family 66 and catalyzes an intramolecular transglucosylation reaction that produces cycloisomaltooligosaccharides from dextran. The crystal structure of the core fragment from Ser-39 to Met-738 of B. circulans T-3040 cycloisomaltooligosaccharide glucanotransferase, devoid of its N-terminal signal peptide and C-terminal nonconserved regions, was determined. The structural model contained one catalytic (β/α)₈-barrel domain and three β-domains. Domain N with an immunoglobulin-like β-sandwich fold was attached to the N terminus; domain C with a Greek key β-sandwich fold was located at the C terminus, and a carbohydrate-binding module family 35 (CBM35) β-jellyroll domain B was inserted between the 7th β-strand and the 7th α-helix of the catalytic domain A. The structures of the inactive catalytic nucleophile mutant enzyme complexed with isomaltohexaose, isomaltoheptaose, isomaltooctaose, and cycloisomaltooctaose revealed that the ligands bound in the catalytic cleft and the sugar-binding site of CBM35. Of these, isomaltooctaose bound in the catalytic site extended to the second sugar-binding site of CBM35, which acted as subsite −8, representing the enzyme·substrate complex when the enzyme produces cycloisomaltooctaose. The isomaltoheptaose and cycloisomaltooctaose bound in the catalytic cleft with a circular structure around Met-310, representing the enzyme·product complex. These structures collectively indicated that CBM35 functions in determining the size of the product, causing the predominant production of cycloisomaltooctaose by the enzyme. The canonical sugar-binding site of CBM35 bound the mid-part of isomaltooligosaccharides, indicating that the original function involved substrate binding required for efficient catalysis.     

Fundamental and biotechnological applications of neutron scattering measurements for macromolecular dynamics

To explore macromolecular dynamics on the picosecond timescale, we used neutron spectroscopy. First, molecular dynamics were analyzed for the hyperthermophile malate dehydrogenase from Methanococcus jannaschii and a mesophilic homologue, the lactate dehydrogenase from Oryctolagus cunniculus muscle. Hyperthermophiles have elaborate molecular mechanisms of adaptation to extremely high temperature. Using a novel elastic neutron scattering approach that provides independent measurements of the global flexibility and of the structural resilience (rigidity), we have demonstrated that macromolecular dynamics represents one of these molecular mechanisms of thermoadaptation. The flexibilities were found to be similar for both enzymes at their optimal activity temperature and the resilience is higher for the hyperthermophilic protein. Secondly, macromolecular motions were examined in a native and immobilized dihydrofolate reductase (DHFR) from Escherichia coli. The immobilized mesophilic enzyme has increased stability and decreased activity, so that its properties are changed to resemble those of the thermophilic enzyme. Are these changes reflected in dynamical behavior? For this study, we performed quasielastic neutron scattering measurements to probe the protein motions. The residence time is 7.95 ps for the native DHFR and 20.36 ps for the immobilized DHFR. The average height of the potential barrier to local motions is therefore increased in the immobilized DHFR, with a difference in activation energy equal to 0.54 kcal/mol, which is, using the theoretical rate equation, of the same order than expected from calculation.     

An antibody probe to determine the native species of glycinamide ribonucleotide transformylase in chicken liver

Antibody probes of Western blots [Renart, J., Reiser, J., & Stark, G. (1979) Proc. Natl. Acad. Sci. U.S.A. 76, 3116] of chicken liver homogenates under various conditions revealed that glycinamide ribonucleotide transformylase can be rapidly proteolyzed in such homogenates. These findings, along with molecular weight measurements by ultracentrifugation, identify the true form of glycinamide ribonucleotide transformylase as a monomeric protein of 117000 daltons. This protein has been purified 400-fold in 44% yield from chicken liver in one step on an affinity column of 10-formyl-5,8-dideazafolate-Sepharose. Native glycinamide ribonucleotide transformylase retains full activity after proteolytic cleavage to a form (Mr 55000) similar to fragments seen in the Western blot of the homogenates. This phenomenon may be responsible for the previous identification of glycinamide ribonucleotide (GAR) transformylase as a dimer of 55000-dalton subunits. Similar analyses using antibodies to 5-aminoimidazole-4-carboxamide ribonucleotide (AICAR) transformylase [Mueller, W. T., & Benkovic, S. J. (1981) Biochemistry 20, 337] and trifunctional enzyme [Smith, G. K., Mueller, W. T., Wasserman, G. F., Taylor, W. D., & Benkovic, S. J. (1980) Biochemistry 19, 4313] confirm that these two proteins were isolated in their native forms.     

Reversible alkylation of an active site methionine residue in dehydroquinase

Iodoacetic acid inactivates dehydroquinase by simultaneously alkylating 2 methionine residues (Met-23 and Met-205), presumed to be active site residues (described in Kleanthous, C., Campbell, D. G., and Coggins, J. R. (1990) J. Biol. Chem. 265, 10929-10934). Although both sites are carboxymethylated to the same degree in the inactivated enzyme, the modification of Met-205 may be reversed by treatment with mercaptoethanol at alkaline pH, as shown by the stoichiometric loss of label from this site. This, in turn, leads to partial reactivation of the inactive enzyme. Alkylation of Met-23 is not reversible under these conditions. The chemistry of the cleavage reaction at Met-205 was investigated by isolating the cleavage product which was identified by mass spectrometry as the ammonium salt of 2-hydroxyethyl thioacetate. This result is consistent with nucleophilic attack by the thiolate anion of mercaptoethanol on the alpha-carbon of the carboxymethyl moiety, which restores the side chain of the methionine residue (Met-205) and liberates 2-hydroxyethyl thioacetate. The differential reactivity of the 2 carboxymethylated methionine residues toward mercaptoethanol is likely to be a reflection of their different microenvironments in the folded protein. This assertion is borne out by unfolding experiments which indicate that neither of the carboxymethylated methionine residues in dicarboxymethylated dehydroquinase is susceptible to mercaptoethanol cleavage if the protein is first denatured by either guanidine hydrochloride or urea. Furthermore, this denatured material refolds after removal of denaturant to yield protein with reactivation properties similar to untreated, dicarboxymethylated enzyme.     

Insight into the molecular mechanism about lowered dihydrofolate binding affinity to dihydrofolate reductase-like 1 (DHFRL1)

Human dihydrofolate reductase-like 1 (DHFRL1) has been identified as a second human dihydrofolate reductase (DHFR) enzyme. Although DHFRL1 have high sequence homology with human DHFR, dihydrofolate (DHF) exhibits a lowered binding affinity to DHFRL1 and the corresponding molecular mechanism is still unknown. To address this question, we studied the binding of DHF to DHFRL1 and DHFR by using molecular dynamics simulation. Moreover, to investigate the role the 24th residue of DHFR/DHFRL1 plays in DHF binding, R24W DHFRL1 mutant was also studied. The van der Waals interaction are more crucial for the total DHF binding energies, while the difference between the DHF binding energies of human DHFR and DHFRL1 can be attributed to the electrostatic interaction and the polar desolvation free energy. More specifically, lower DHF affinity to DHFRL1 can be mainly attributed to the reduction of net electrostatic interactions of residues Arg32 and Gln35 of DHFRL1 with DHF as being affected by Arg24. The side chain of Arg24 in DHFRL1 can extend deeply into the binding sites of DHF and NADPH, and disturb the DHF binding by steric effect, which rarely happens in human DHFR and R24W DHFRL1 mutant. Additionally, the conformation of loop I in DHFRL1 was also studied in this work. Interestingly, the loop conformation resemble to normal closed state of Escherichia coli DHFR other than the closed state of human DHFR. We hope this work will be useful to understand the general characteristics of DHFRL1.     

Numerous Classes of General Anesthetics Inhibit Etomidate Binding to ��-Aminobutyric Acid Type A (GABAA) Receptors

Enhancement of γ-aminobutyric acid type A receptor (GABA_AR)-mediated inhibition is a property of most general anesthetics and a candidate for a molecular mechanism of anesthesia. Intravenous anesthetics, including etomidate, propofol, barbiturates, and neuroactive steroids, as well as volatile anesthetics and long-chain alcohols, all enhance GABA_AR function at anesthetic concentrations. The implied existence of a receptor site for anesthetics on the GABA_AR protein was supported by identification, using photoaffinity labeling, of a binding site for etomidate within the GABA_AR transmembrane domain at the β-α subunit interface; the etomidate analog [³H]azietomidate photolabeled in a pharmacologically specific manner two amino acids, α1Met-236 in the M1 helix and βMet-286 in the M3 helix (Li, G. D., Chiara, D. C., Sawyer, G. W., Husain, S. S., Olsen, R. W., and Cohen, J. B. (2006) J. Neurosci. 26, 11599–11605). Here, we use [³H]azietomidate photolabeling of bovine brain GABA_ARs to determine whether other structural classes of anesthetics interact with the etomidate binding site. Photolabeling was inhibited by anesthetic concentrations of propofol, barbiturates, and the volatile agent isoflurane, at low millimolar concentrations, but not by octanol or ethanol. Inhibition by barbiturates, which was pharmacologically specific and stereospecific, and by propofol was only partial, consistent with allosteric interactions, whereas isoflurane inhibition was nearly complete, apparently competitive. Protein sequencing showed that propofol inhibited to the same extent the photolabeling of α1Met-236 and βMet-286. These results indicate that several classes of general anesthetics modulate etomidate binding to the GABA_AR: isoflurane binds directly to the site with millimolar affinity, whereas propofol and barbiturates inhibit binding but do not bind in a mutually exclusive manner with etomidate.     

A Triple Mutant in the Ω-loop of TEM-1 β-Lactamase Changes the Substrate Profile via a Large Conformational Change and an Altered General Base for Catalysis

β-Lactamases are bacterial enzymes that hydrolyze β-lactam antibiotics. TEM-1 is a prevalent plasmid-encoded β-lactamase in Gram-negative bacteria that efficiently catalyzes the hydrolysis of penicillins and early cephalosporins but not oxyimino-cephalosporins. A previous random mutagenesis study identified a W165Y/E166Y/P167G triple mutant that displays greatly altered substrate specificity with increased activity for the oxyimino-cephalosporin, ceftazidime, and decreased activity toward all other β-lactams tested. Surprisingly, this mutant lacks the conserved Glu-166 residue critical for enzyme function. Ceftazidime contains a large, bulky side chain that does not fit optimally in the wild-type TEM-1 active site. Therefore, it was hypothesized that the substitutions in the mutant expand the binding site in the enzyme. To investigate structural changes and address whether there is an enlargement in the active site, the crystal structure of the triple mutant was solved to 1.44 Å. The structure reveals a large conformational change of the active site Ω-loop structure to create additional space for the ceftazidime side chain. The position of the hydroxyl group of Tyr-166 and an observed shift in the pH profile of the triple mutant suggests that Tyr-166 participates in the hydrolytic mechanism of the enzyme. These findings indicate that the highly conserved Glu-166 residue can be substituted in the mechanism of serine β-lactamases. The results reveal that the robustness of the overall β-lactamase fold coupled with the plasticity of an active site loop facilitates the evolution of enzyme specificity and mechanism.     

Theoretical Investigation on the Binding Specificity of Sialyldisaccharides with Hemagglutinins of Influenza A Virus by Molecular Dynamics Simulations

Recognition of cell-surface sialyldisaccharides by influenza A hemagglutinin (HA) triggers the infection process of influenza. The changes in glycosidic torsional linkage and the receptor conformations may alter the binding specificity of HAs to the sialylglycans. In this study, 10-ns molecular dynamics simulations were carried out to examine the structural and dynamic behavior of the HAs bound with sialyldisaccharides Neu5Acα(2–3)Gal (N23G) and Neu5Acα(2–6)Gal (N26G). The analysis of the glycosidic torsional angles and the pair interaction energy between the receptor and the interacting residues of the binding site reveal that N23G has two binding modes for H1 and H5 and a single binding mode for H3 and H9. For N26G, H1 and H3 has two binding modes, and H5 and H9 has a single binding mode. The direct and water-mediated hydrogen bonding interactions between the receptors and HAs play dominant roles in the structural stabilization of the complexes. It is concluded from pair interaction energy and Molecular Mechanic-Poisson-Boltzmann Surface Area calculations that N26G is a better receptor for H1 when compared with N23G. N23G is a better receptor for H5 when compared with N26G. However, H3 and H9 can recognize N23G and N26G in equal binding specificity due to the marginal energy difference (≈2.5 kcal/mol). The order of binding specificity of N23G is H3 > H5 > H9 > H1 and N26G is H1 > H3 > H5 > H9, respectively. The proposed conformational models will be helpful in designing inhibitors for influenza virus.     

Unique Structural Characteristics of the Rabbit Prion Protein

Rabbits are one of the few mammalian species that appear to be resistant to transmissible spongiform encephalopathies due to the structural characteristics of the rabbit prion protein (RaPrP^C) itself. Here, we determined the solution structures of the recombinant protein RaPrP^C-(91–228) and its S173N variant and detected the backbone dynamics of their structured C-terminal domains-(121–228). In contrast to many other mammalian PrP^Cs, loop 165–172, which connects β-sheet-2 and α-helix-2, is well-defined in RaPrP^C. For the first time, order parameters S² are obtained for residues in this loop region, indicating that loop 165–172 of RaPrP^C is highly ordered. Compared with the wild-type RaPrP^C, less hydrogen bonds form in the S173N variant. The NMR dynamics analysis reveals a distinct increase in the structural flexibility of loop 165–172 and helix-3 after the S173N substitution, implying that the S173N substitution disturbs the long range interaction of loop 165–172 with helix-3, which further leads to a marked decrease in the global conformational stability. Significantly, RaPrP^C possesses a unique charge distribution, carrying a continuous area of positive charges on the surface, which is distinguished from other PrP^Cs. The S173N substitution causes visible changes of the charge distribution around the recognition sites for the hypothetical protein X. Our results suggest that the ordered loop 165–172 and its interaction with helix-3, together with the unique distribution of surface electrostatic potential, significantly contribute to the unique structural characteristics of RaPrP^C.     

Multi-timescale dynamics study of FKBP12 along the rapamycin-mTOR binding coordinate

Drugs can affect function in proteins by modulating their flexibility. Despite this possibility, there are very few studies on how drug binding affects the dynamics of target macromolecules. FKBP12 (FK506 binding protein 12) is a prolyl cis-trans isomerase and a drug target. The immunosuppressant drug rapamycin exerts its therapeutic effect by serving as an adaptor molecule between FKBP12 and the cell proliferation regulator mTOR (mammalian target of rapamycin). To understand the role of dynamics in rapamycin-based immunosuppression and to gain insight into the role of dynamics in the assembly of supramolecular complexes, we used ¹⁵N, ¹³C, and ²H NMR spin relaxation to characterize FKBP12 along the binding coordinate that leads to cell cycle arrest. We show that sequential addition of rapamycin and mTOR leads to incremental rigidification of the FKBP12 backbone on the picosecond-nanosecond timescale. Both binding events lead to perturbation of main-chain and side-chain dynamics at sites distal to the binding interfaces, suggesting tight coupling interactions dispersed throughout the FKBP12-rapamycin interface. Binding of the first molecule, rapamycin, quenches microsecond-millisecond motions of the FKBP12 80's loop. This loop provides much of the surface buried at the protein-protein interface of the ternary complex, leading us to assert that preorganization upon rapamycin binding facilitates binding of the second molecule, mTOR. Widespread microsecond-millisecond motions of the backbone persist in the drug-bound enzyme, and we provide evidence that these slow motions represent coupled dynamics of the enzyme and isomerization of the bound drug. Finally, the pattern of microsecond-millisecond dynamics reported here in the rapamycin complex is dramatically different from the pattern in the complex with the structurally related drug FK506. This raises the important question of how two complexes that are highly isomorphic based on high-resolution static models have such different flexibilities in solution.