Diagnostic chemical shift markers for loop conformation and substrate and cofactor binding in dihydrofolate reductase complexes

Heteronuclear NMR methods have been used to probe the conformation of four complexes of Escherichia coli dihydrofolate reductase (DHFR) in solution. (1)H(N), (15)N, and (13)C(alpha) resonance assignments have been made for the ternary complex with folate and oxidized NADP(+) cofactor and the ternary complex with folate and a reduced cofactor analog, 5,6-dihydroNADPH. The backbone chemical shifts have been compared with those of the binary complex of DHFR with the substrate analog folate and the binary complex with NADPH (the holoenzyme). Analysis of (1)H(N) and (15)N chemical shifts has led to the identification of marker resonances that report on the active site conformation of the enzyme. Other backbone amide resonances report on the presence of ligands in the pterin binding pocket and in the adenosine and nicotinamide-ribose binding sites of the NADPH cofactor. The chemical shift data indicate that the enzyme populates two dominant structural states in solution, with the active site loops in either the closed or occluded conformations defined by X-ray crystallography; there is no evidence that the open conformation observed in some X-ray structures of E. coli DHFR are populated in solution.     

Comparison of the structures of operator DNA free and in complex with lambda repressor

The structures of operator DNA unbound and in complex with lambda repressor protein are compared. The conformation of the left 10 base pairs of a lambda right regulatory operator DNA sequence has been previously determined in solution using nuclear magnetic resonance techniques and the structure of a homologous left regulatory operator DNA bound to lambda repressor N-terminal domain had been previously solved using X-ray crystallography. The DNA adopts an overall linear B-form DNA both in the absence and presence of lambda repressor. Superimpositioning of the DNA structures reveals small differences between them that are due to the binding of protein and not to the different techniques used for their determination.     

Conformational changes in the active site loops of dihydrofolate reductase during the catalytic cycle

Escherichia coli dihydrofolate reductase (DHFR) has several flexible loops surrounding the active site that play a functional role in substrate and cofactor binding and in catalysis. We have used heteronuclear NMR methods to probe the loop conformations in solution in complexes of DHFR formed during the catalytic cycle. To facilitate the NMR analysis, the enzyme was labeled selectively with [(15)N]alanine. The 13 alanine resonances provide a fingerprint of the protein structure and report on the active site loop conformations and binding of substrate, product, and cofactor. Spectra were recorded for binary and ternary complexes of wild-type DHFR bound to the substrate dihydrofolate (DHF), the product tetrahydrofolate (THF), the pseudosubstrate folate, reduced and oxidized NADPH cofactor, and the inactive cofactor analogue 5,6-dihydroNADPH. The data show that DHFR exists in solution in two dominant conformational states, with the active site loops adopting conformations that closely approximate the occluded or closed conformations identified in earlier X-ray crystallographic analyses. A minor population of a third conformer of unknown structure was observed for the apoenzyme and for the disordered binary complex with 5,6-dihydroNADPH. The reactive Michaelis complex, with both DHF and NADPH bound to the enzyme, could not be studied directly but was modeled by the ternary folate:NADP(+) and dihydrofolate:NADP(+) complexes. From the NMR data, we are able to characterize the active site loop conformation and the occupancy of the substrate and cofactor binding sites in all intermediates formed in the extended catalytic cycle. In the dominant kinetic pathway under steady-state conditions, only the holoenzyme (the binary NADPH complex) and the Michaelis complex adopt the closed loop conformation, and all product complexes are occluded. The catalytic cycle thus involves obligatory conformational transitions between the closed and occluded states. Parallel studies on the catalytically impaired G121V mutant DHFR show that formation of the closed state, in which the nicotinamide ring of the cofactor is inserted into the active site, is energetically disfavored. The G121V mutation, at a position distant from the active site, interferes with coupled loop movements and appears to impair catalysis by destabilizing the closed Michaelis complex and introducing an extra step into the kinetic pathway.     

Nuclear magnetic resonance detection of bound water molecules in the active site of Lactobacillus casei dihydrofolate reductase in aqueous solution.

Proton nuclear magnetic resonance spectroscopy has been used to detect two water molecules bound to residues in the active site of the Lactobacillus casei dihydrofolate reductase (DHFR). Their presence was detected by measuring nuclear Overhauser effects between NH protons in protein residues and protons in the individual bound water molecules in two-dimensional nuclear Overhauser effect spectroscopy (NOESY), in nuclear Overhauser effect spectroscopy in the rotating frame (ROESY) and three-dimensional 1H-15N ROESY-heteronuclear multiple quantum coherence spectra recorded on samples containing appropriately 15N-labelled DHFR. For the DHFR-methotrexate-NADPH complex, two bound molecules were found, one close to the Trp5 amide NH proton and the other near to the Trp21 indole HE1 proton: these correspond to two of the water molecules (Wat201 and Wat253) detected in the crystal structure studies described by Bolin and co-workers. However, the nuclear magnetic resonance experiments did not detect any of the other bound water molecules observed in the X-ray studies. The nuclear magnetic resonance results indicate that the two bound water molecules that were detected have lifetimes in the solution state that are longer than approximately two nanoseconds. This is of considerable interest, since one of these water molecules (Wat253) has been implicated as the likely proton donor in the catalytic reduction of dihydrofolate to tetrahydrofolate.     

Solution nuclear magnetic resonance structure of a protein disulfide oxidoreductase from Methanococcus jannaschii

The solution structure of the protein disulfide oxidoreductase Mj0307 in the reduced form has been solved by nuclear magnetic resonance. The secondary and tertiary structure of this protein from the archaebacterium Methanococcus jannaschii is similar to the structures that have been solved for the glutaredoxin proteins from Escherichia coli, although Mj0307 also shows features that are characteristic of thioredoxin proteins. Some aspects of Mj0307's unique behavior can be explained by comparing structure-based sequence alignments with mesophilic bacterial and eukaryotic glutaredoxin and thioredoxin proteins. It is proposed that Mj0307, and similar archaebacterial proteins, may be most closely related to the mesophilic bacterial NrdH proteins. Together these proteins may form a unique subgroup within the family of protein disulfide oxidoreductases.     

Recent developments in structural proteomics for protein structure determination

The major challenges in structural proteomics include identifying all the proteins on the genome-wide scale, determining their structure-function relationships, and outlining the precise three-dimensional structures of the proteins. Protein structures are typically determined by experimental approaches such as X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. However, the knowledge of three-dimensional space by these techniques is still limited. Thus, computational methods such as comparative and de novo approaches and molecular dynamic simulations are intensively used as alternative tools to predict the three-dimensional structures and dynamic behavior of proteins. This review summarizes recent developments in structural proteomics for protein structure determination; including instrumental methods such as X-ray crystallography and NMR spectroscopy, and computational methods such as comparative and de novo structure prediction and molecular dynamics simulations.     

Nuclear magnetic resonance solution and X-ray structures of squash trypsin inhibitor exhibit the same conformation of the proteinase binding loop

A comparison of the solution nuclear magnetic resonance (n.m.r.) structures of squash trypsin inhibitor from seeds of the squash Cucurbita maxima with the X-ray structure of a trypsin complex of the inhibitor shows that the n.m.r. and X-ray structures are similar in terms of the global folding and secondary structure. The average atomic root-mean-square difference between the 36 n.m.r. structures on the one hand and the X-ray structure is 0.96 A for the backbone atoms and 1.95 A for all heavy atoms. The n.m.r. and X-ray structures exhibit extremely similar conformations of the primary proteinase binding loop. Despite the overall similarity, there are small differences between the mean computed structure and the X-ray structure. The n.m.r. structures have slightly different positions of the segments from residues 16 to 18, and 24 and 25. The n.m.r. results show that the inclusion of stereospecific assignments and precise distance constraints results in a significant improvement in the definition of the n.m.r. structure, making possible a detailed analysis of the local conformations in the protein.     

Nuclear magnetic resonance determination of metal-proton distances in a synthetic calcium binding site of rabbit skeletal troponin C

The binding of gadolinium to a synthetic peptide of 13 amino acid residues representing the calcium binding loop of site 3 of rabbit skeletal troponin C [AcSTnC(103-115)amide] has been studied by using proton nuclear magnetic resonance (1H NMR) spectroscopy. In particular, the proton line broadening and enhanced spin-lattice relaxation have been used to determine proton-metal ion distances for several assigned nuclei in the peptide-metal ion complex. These distances have been used in conjunction with other constraints and a distance algorithm procedure to demonstrate that the structure of the peptide-metal complex as shown by 1H NMR is consistent with the structure of the EF calcium binding loop in the X-ray structure of parvalbumin but that the available 1H NMR distances do not uniquely define the solution structure.     

Ligand-induced conformational changes in the crystal structures of Pneumocystis carinii dihydrofolate reductase complexes with folate and NADP+

Structural data from two independent crystal forms (P212121 and P21) of the folate (FA) binary complex and from the ternary complex with the oxidized coenzyme, NADP+, and recombinant Pneumocystis carinii dihydrofolate reductase (pcDHFR) refined to an average of 2.15 A resolution, show the first evidence of ligand-induced conformational changes in the structure of pcDHFR. These data are also compared with the crystal structure of the ternary complex of methotrexate (MTX) with NADPH and pcDHFR in the monoclinic lattice with data to 2.5 A resolution. Comparison of the data for the FA binary complex of pcDHFR with those for the ternary structures reveals significant differences, with a >7 A movement of the loop region near residue 23 that results in a new "flap-open" position for the binary complex, and a "closed" position in the ternary complexes, similar to that reported for Escherichia coli (ec) DHFR complexes. In the orthorhombic lattice for the binary FA pcDHFR complex, there is also an unwinding of a short helical region near residue 47 that places hydrophobic residues Phe-46 and Phe-49 toward the outer surface, a conformation that is stabilized by intermolecular packing contacts. The pyrophosphate moiety of NADP+ in the ternary folate pcDHFR complexes shows significant differences in conformation compared with that observed in the MTX-NADPH-pcDHFR ternary complex. Additionally, comparison of the conformations among these four pcDHFR structures reveals evidence for subdomain movement that correlates with cofactor binding states. The larger binding site access in the new "flap-open" loop 23 conformation of the binary FA complex is consistent with the rapid release of cofactor from the product complex during catalysis as well as the more rapid release of substrate product from the binary complex as a result of the weaker contacts of the closed loop 23 conformation, compared to ecDHFR.     

NMR structures of apo L. casei dihydrofolate reductase and its complexes with trimethoprim and NADPH: contributions to positive cooperative binding from ligand-induced refolding, conformational changes, and interligand hydrophobic interactions

In order to examine the origins of the large positive cooperativity (ΔG(0)(coop) = -2.9 kcal mol(-1)) of trimethoprim (TMP) binding to a bacterial dihydrofolate reductase (DHFR) in the presence of NADPH, we have determined and compared NMR solution structures of L. casei apo DHFR and its binary and ternary complexes with TMP and NADPH and made complementary thermodynamic measurements. The DHFR structures are generally very similar except for the A-B loop region and part of helix B (residues 15-31) which could not be directly detected for L. casei apo DHFR because of line broadening from exchange between folded and unfolded forms. Thermodynamic and NMR measurements suggested that a significant contribution to the cooperativity comes from refolding of apo DHFR on binding the first ligand (up to -0.95 kcals mol(-1) if 80% of A-B loop requires refolding). Comparisons of Cα-Cα distance differences and domain rotation angles between apo DHFR and its complexes indicated that generally similar conformational changes involving domain movements accompany formation of the binary complexes with either TMP or NADPH and that the binary structures are approaching that of the ternary complex as would be expected for positive cooperativity. These favorable ligand-induced structural changes upon binding the first ligand will also contribute significantly to the cooperative binding. A further substantial contribution to cooperative binding results from the proximity of the bound ligands in the ternary complex: this reduces the solvent accessible area of the ligand and provides a favorable entropic hydrophobic contribution (up to -1.4 kcal mol(-1)).     

Relative stability of protein structures determined by X-ray crystallography or NMR spectroscopy: a molecular dynamics simulation study

The relative stability of protein structures determined by either X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy has been investigated by using molecular dynamics simulation techniques. Published structures of 34 proteins containing between 50 and 100 residues have been evaluated. The proteins selected represent a mixture of secondary structure types including all alpha, all beta, and alpha/beta. The proteins selected do not contain cysteine-cysteine bridges. In addition, any crystallographic waters, metal ions, cofactors, or bound ligands were removed before the systems were simulated. The stability of the structures was evaluated by simulating, under identical conditions, each of the proteins for at least 5 ns in explicit solvent. It is found that not only do NMR-derived structures have, on average, higher internal strain than structures determined by X-ray crystallography but that a significant proportion of the structures are unstable and rapidly diverge in simulations.     

RNA structure: experimental analysis

Among all of the biological macromolecules, the functional versatility of RNAs is unique including encoding or transferring genetic information and performing catalysis. These biological functions are highly dependent upon RNA folding and structure. Since the discovery of catalytic RNAs in the early 1980s, a recent breakthrough came from the identification of a wealth of micro RNAs, small interfering RNAs and regulatory RNAs, all involved in modulation of gene expression. The structure of these novel RNAs, either free or in complex with specific ligands, can be analyzed using various experimental strategies, including X-ray crystallography, cryo-electron microscopy, nuclear magnetic resonance spectroscopy, structure-specific probes, with some that can be used in living cells, RNA engineering, thermal denaturation and mass spectrometry. Among these, X-ray crystallography has recently enabled determination of the structures of several large and complex RNAs, as well as of ribonucleoprotein complexes. The database of RNA structure has grown tremendously since the recent crystal structure analyses of the prokaryotic ribosome and its subunits. These methods are now widely applied to a variety of biologically relevant RNAs.     

Protein thermal stability: the role of protein structure and aqueous environment

A comprehensive bioinformatic analysis was performed on all protein homologous pairs from mesophilic and thermophilic microorganisms present in the RCSB Protein Data Bank in order to yield a clue on the role of protein structure and aqueous environment. Subsequently self-assembly and LB studies were carried out at increasing temperature by nanogravimetry with thermostable thioredoxin (Trx) from Alicyclobacillus acidocaldarius (BacTrx) versus the mesophilic Escherichia coli counterpart (EcTrx). The comparison with earlier 3D atomic structure determined on the same proteins by X-ray crystallographic diffraction and nuclear magnetic resonance confirm the role inner bound water in determining protein thermostability, as suggested by the bioinformatic and nanogravimetric analysis. The above comparative characterizations in protein solution, thin film and crystal allow to draw a possible coherent explanation for the origin and the molecular mechanisms of both heat stability and radiation resistance in proteins.     

Prospects for robust biocatalysis: engineering of novel specificity in a halophilic amino acid dehydrogenase

Heat- and solvent-tolerant enzymes from halophiles, potentially important industrially, offer a robust framework for protein engineering, but few solved halophilic structures exist to guide this. Homology modelling has guided mutations in glutamate dehydrogenase (GDH) from Halobacterium salinarum to emulate conversion of a mesophilic GDH to a methionine dehydrogenase. Replacement of K89, A163 and S367 by leucine, glycine and alanine converted halophilic GDH into a dehydrogenase accepting l-methionine, l-norleucine and l-norvaline as substrates. Over-expression in the halophilic expression host Haloferax volcanii and three-step purification gave ~98 % pure protein exhibiting maximum activity at pH 10. This enzyme also showed enhanced thermostability and organic solvent tolerance even at 70 °C, offering a biocatalyst resistant to harsh industrial environments. To our knowledge, this is the first reported amino acid specificity change engineered in a halophilic enzyme, encouraging use of mesophilic models to guide engineering of novel halophilic biocatalysts for industrial application. Calibrated gel filtration experiments show that both the mutant and the wild-type enzyme are stable hexamers.     

Mechanistic analysis of allosteric and non-allosteric effects arising from nanobody binding to two epitopes of the dihydrofolate reductase of Escherichia coli

Although allosteric effector antibodies are used widely as modulators of receptors and enzymes, experimental analysis of their mechanism remains highly challenging. Here, we investigate the molecular mechanisms of allosteric and non-allosteric effector antibodies in an experimentally tractable system, consisting of single-domain antibodies (nanobodies) that target the model enzyme dihydrofolate reductase (DHFR) from Escherichia coli. A panel of thirty-five nanobodies was isolated using several strategies to increase nanobody diversity. The nanobodies exhibit a variety of effector properties, including partial inhibition, strong inhibition and stimulation of DHFR activity. Despite these diverse effector properties, chemical shift perturbation NMR epitope mapping identified only two epitope regions: epitope α is a new allosteric site that is over 10 Å from the active site, while epitope β is located in the region of the Met20 loop. The structural basis for DHFR allosteric inhibition or activation upon nanobody binding to the α epitope was examined by solving the crystal structures of DHFR in complex with Nb113 (an allosteric inhibitor) and Nb179 (an allosteric activator). The structures suggest roles for conformational constraint and altered protein dynamics, but not epitope distortion, in the observed allosteric effects. The crystal structure of a β epitope region binder (ca1698) in complex with DHFR is also reported. Although CDR3 of ca1698 occupies the substrate binding site, ca1698 displays linear mixed inhibition kinetics instead of simple competitive inhibition kinetics. Two mechanisms are proposed to account for this apparent anomaly. Evidence for structural convergence of ca1698 and Nb216 during affinity maturation is also presented.     

Solution structure of the Grb2 SH2 domain complexed with a high-affinity inhibitor

The solution structure of the growth factor receptor-bound protein 2 (Grb2) SH2 domain complexed with a high-affinity inhibitor containing a non-phosphorus phosphate mimetic within a macrocyclic platform was determined by nuclear magnetic resonance (NMR) spectroscopy. Unambiguous assignments of the bound inhibitor and intermolecular NOEs between the Grb2 SH2 domain and the inhibitor was accomplished using perdeuterated Grb2 SH2 protein. The well-defined solution structure of the complex was obtained and compared to those by X-ray crystallography. Since the crystal structure of the Grb2 SH2 domain formed a domain-swapped dimer and several inhibitors were bound to a hinge region, there were appreciable differences between the solution and crystal structures. Based on the binding interactions between the inhibitor and the Grb2 SH2 domain in solution, we proposed a design of second-generation inhibitors that could be expected to have higher affinity.     

Tryptophan fluorescence monitors multiple conformational changes required for glutamine phosphoribosylpyrophosphate amidotransferase interdomain signaling and catalysis.

Single tryptophan residues were incorporated into each of three peptide segments that play key roles in the structural transition of ligand-free, inactive glutamine phosphoribosylpyrophosphate (PRPP) amidotransferase to the active enzyme-substrate complex. Intrinsic tryptophan fluorescence and fluorescence quenching were used to monitor changes in a phosphoribosyltransferase (PRTase) "flexible loop", a "glutamine loop", and a C-terminal helix. Steady state fluorescence changes resulting from substrate binding were used to calculate binding constants and to detect the structural rearrangements that coordinate reactions at active sites for glutamine hydrolysis and PRTase catalysis. Pre-steady state kinetics of enzyme.PRPP and enzyme.PRPP.glutamine complex formation were determined from stopped-flow fluorescence measurements. The kinetics of the formation of the enzyme.PRPP complex were consistent with a model with two or more steps in which rapid equilibrium binding of PRPP is followed by a slow enzyme isomerization. This isomerization is ascribed to the closing of the PRTase flexible loop and is likely the rate-limiting step in the reaction of PRPP with NH(3). The pre-steady state kinetics for binding glutamine to the binary enzyme. PRPP complex could also be fit to a model involving rapid equilibrium binding of glutamine followed by an enzyme isomerization step. The changes monitored by fluorescence account for the interconversions between "end state" structures determined previously by X-ray crystallography and define an intermediate enzyme.PRPP conformer.     

Comparison of multiple crystal structures with NMR data for engrailed homeodomain

Two methods are currently available to solve high resolution protein structures—X-ray crystallography and nuclear magnetic resonance (NMR). Both methods usually produce highly similar structures, but small differences between both solutions are always observed. Here the raw NMR data as well as the solved NMR structure were compared to the multiple crystal structures solved for the WT 60 residue three helix bundle engrailed homeodomain (EnHD) and single point mutants. There was excellent agreement between TALOS-predicted and crystal structure-observed dihedral angles and a good agreement for the ³ J(H ^N H ^α ) couplings for the multiple crystal structures. Around 1% of NOEs were violated for any crystal structure, but no NOE was inconsistent with all of the crystal structures. Violations usually occurred for surface residues or for residues for which multiple discreet conformations were observed between the crystal structures. Comparison of the disorder shown in the multiple crystal structures shows little correlation with dynamics under native conditions for this protein.