1
|
Mackenzie HW, Hansen DF. A 13C-detected 15N double-quantum NMR experiment to probe arginine side-chain guanidinium 15N η chemical shifts. JOURNAL OF BIOMOLECULAR NMR 2017; 69:123-132. [PMID: 29127559 PMCID: PMC5711973 DOI: 10.1007/s10858-017-0137-2] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/23/2017] [Accepted: 09/25/2017] [Indexed: 05/25/2023]
Abstract
Arginine side-chains are often key for enzyme catalysis, protein-ligand and protein-protein interactions. The importance of arginine stems from the ability of the terminal guanidinium group to form many key interactions, such as hydrogen bonds and salt bridges, as well as its perpetual positive charge. We present here an arginine 13Cζ-detected NMR experiment in which a double-quantum coherence involving the two 15Nη nuclei is evolved during the indirect chemical shift evolution period. As the precession frequency of the double-quantum coherence is insensitive to exchange of the two 15Nη; this new approach is shown to eliminate the previously deleterious line broadenings of 15Nη resonances caused by the partially restricted rotation about the Cζ-Nε bond. Consequently, sharp and well-resolved 15Nη resonances can be observed. The utility of the presented method is demonstrated on the L99A mutant of the 19 kDa protein T4 lysozyme, where the measurement of small chemical shift perturbations, such as one-bond deuterium isotope shifts, of the arginine amine 15Nη nuclei becomes possible using the double-quantum experiment.
Collapse
Affiliation(s)
- Harold W Mackenzie
- Institute of Structural and Molecular Biology, Division of Biosciences, University College London, London, WC1E 6BT, UK
| | - D Flemming Hansen
- Institute of Structural and Molecular Biology, Division of Biosciences, University College London, London, WC1E 6BT, UK.
| |
Collapse
|
2
|
Structural basis for sulfation-dependent self-glycan recognition by the human immune-inhibitory receptor Siglec-8. Proc Natl Acad Sci U S A 2016; 113:E4170-9. [PMID: 27357658 DOI: 10.1073/pnas.1602214113] [Citation(s) in RCA: 66] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
Siglec-8 is a human immune-inhibitory receptor that, when engaged by specific self-glycans, triggers eosinophil apoptosis and inhibits mast cell degranulation, providing an endogenous mechanism to down-regulate immune responses of these central inflammatory effector cells. Here we used solution NMR spectroscopy to dissect the fine specificity of Siglec-8 toward different sialylated and sulfated carbohydrate ligands and determined the structure of the Siglec-8 lectin domain in complex with its prime glycan target 6'-sulfo sialyl Lewis(x) A canonical motif for sialic acid recognition, extended by a secondary motif formed by unique loop regions, recognizing 6-O-sulfated galactose dictates tight specificity distinct from other Siglec family members and any other endogenous glycan recognition receptors. Structure-guided mutagenesis revealed key contacts of both interfaces to be equally essential for binding. Our work provides critical structural and mechanistic insights into how Siglec-8 selectively recognizes its glycan target, rationalizes the functional impact of site-specific glycan sulfation in modulating this lectin-glycan interaction, and will enable the rational design of Siglec-8-targeted agonists to treat eosinophil- and mast cell-related allergic and inflammatory diseases, such as asthma.
Collapse
|
3
|
Biris N, Stavrakoudis A, Politou AS, Mikros E, Sakarellos-Daitsiotis M, Sakarellos C, Tsikaris V. The Ac-RGD-NH2 peptide as a probe of slow conformational exchange of short linear peptides in DMSO. Biopolymers 2003; 69:72-86. [PMID: 12717723 DOI: 10.1002/bip.10335] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
According to general belief, the conformational information on short linear peptides in solution derived at ambient temperature from NMR spectrometry represents a population-weighted average over all members of an ensemble of rapidly interconverting conformations. Usually the search for discrete conformations is concentrated at low temperatures especially when sharp NMR resonances are detected at room temperature. Using the peptide Ac-RGD-NH(2) (Ac-Arg-Gly-Asp-NH(2), Ac: acetyl) as a model system and following a new approach, we have been able to demonstrate that short linear peptides can adopt discrete conformational states in DMSO-d(6) (DMSO: dimethylsulfoxide) which vary in a way critically dependent on the reconstitution conditions used before their dissolution in DMSO-d(6). The conformers are stabilized by intramolecular hydrogen bonds, which persist at high temperatures and undergo a very slow exchange with their extended structures in the NMR chemical shift time scale. The reported findings provide clear evidence for the occurrence of solvent-induced conformational exchange and point to DMSO as a valuable medium for folding studies of short linear peptides.
Collapse
|
4
|
Birdsall B, Polshakov VI, Feeney J. NMR studies of ligand carboxylate group interactions with arginine residues in complexes of Lactobacillus casei dihydrofolate reductase with substrates and substrate analogues. Biochemistry 2000; 39:9819-25. [PMID: 10933799 DOI: 10.1021/bi000728s] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In a series of complexes of Lactobacillus casei dihydrofolate reductase (DHFR) formed with substrates and substrate analogues, the (1)H/(15)N NMR chemical shifts for the guanidino group of the conserved Arg 57 residue were found to be sensitive to the mode of binding of their H(eta) protons to the charged oxygen atoms in ligand carboxylate groups. In all cases, Arg 57 showed four nonequivalent H(eta) signals indicating hindered rotation about the N(epsilon)-C(zeta) and C(zeta)-N(eta) bonds. The H(eta)(12) and H(eta)(22) protons have large downfield shifts as expected for a symmetrical end-on interaction with the ligand carboxylate group. The chemical shifts are essentially the same in the complexes with folate and p-aminobenzoyl-L-glutamate (PABG) and similar to those found previously for the methotrexate complex reflecting the strong and similar hydrogen bonds formed with the carboxylate oxygens. Interestingly, the rates of rotation about the N(epsilon)-C(zeta) bond for the complexes containing the weakly binding PABG fragment are almost identical to those measured in the complex with methotrexate, which binds 10(7) times more tightly. In the methotrexate complex, this rotation depends on correlated rotations about the N(epsilon)-C(zeta) bond of Arg 57 and the C(alpha)-C' bond of the ligand glutamate alpha-carboxylate group. Thus, even in a fragment such as PABG, which has a much faster off-rate, the carboxylate group binds to the enzyme in a similar way to that in a parent molecule such as folate and methotrexate with the rotation about the N(epsilon)-C(zeta) bond of Arg 57 being essentially the same in all the different complexes.
Collapse
Affiliation(s)
- B Birdsall
- Molecular Structure Division, National Institute for Medical Research, The Ridgeway, Mill Hill, London NW7 1AA, U.K
| | | | | |
Collapse
|
5
|
|
6
|
Polshakov VI, Birdsall B, Feeney J. Characterization of rates of ring-flipping in trimethoprim in its ternary complexes with Lactobacillus casei dihydrofolate reductase and coenzyme analogues. Biochemistry 1999; 38:15962-9. [PMID: 10625463 DOI: 10.1021/bi9915263] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
NMR measurements have been used to investigate rates of ring-flipping and the activation parameters for the trimethoxybenzyl ring of the antibacterial drug trimethoprim (TMP) bound to Lactobacillus casei dihydrofolate reductase (DHFR) for a series of ternary complexes formed with analogues of the coenzyme NADPH. Rates were obtained at several temperatures from line shape analyses ((13)C-edited HSQC (1)H spectra) and transfer of magnetization measurements (zz-HSQC) on complexes containing 3'-O-[(13)C]trimethoprim. Examination of the structures of the complexes indicates that ring-flipping can only be achieved following major conformational changes and transient fluctuations of the protein and coenzyme structure around the trimethoxybenzyl ring. There is no simple correlation between rates of ring-flipping and binding constants. The presence of the coenzyme nicotinamide ring (in either its reduced or its oxidized forms) in the binding site close to the trimethoxybenzyl ring moiety is the major factor reducing the ring-flipping on coenzyme binding. Thus, the ternary complex with NADPH shows the largest reduction in the rate of ring-flipping (11 +/- 3 s(-)(1) at 298 K) as compared with the binary complex (793 +/- 80 s(-)(1) at 298 K). Complexes with NADPH analogues that either have no nicotinamide ring or are known to have their nicotinamide rings removed from the binding site show the smallest reductions. For the DHFR.TMP.NADP(+) complex where there are two conformations present, very different rates of ring-flipping were observed for the two forms. The activation parameters (DeltaH() and DeltaS()) for the ring-flipping in all the complexes are discussed in terms of the protein-ligand interactions and the possible constraints on the pathway through the transition state.
Collapse
Affiliation(s)
- V I Polshakov
- Division of Molecular Structure, National Institute for Medical Research, The Ridgeway, Mill Hill, London NW7 1AA, U.K., and Center for Drug Chemistry, Moscow 119815, Russia
| | | | | |
Collapse
|
7
|
Polshakov VI, Birdsall B, Frenkiel TA, Gargaro AR, Feeney J. Structure and dynamics in solution of the complex of Lactobacillus casei dihydrofolate reductase with the new lipophilic antifolate drug trimetrexate. Protein Sci 1999; 8:467-81. [PMID: 10091649 PMCID: PMC2144292 DOI: 10.1110/ps.8.3.467] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
We have determined the three-dimensional solution structure of the complex of Lactobacillus casei dihydrofolate reductase and the anticancer drug trimetrexate. Two thousand seventy distance, 345 dihedral angle, and 144 hydrogen bond restraints were obtained from analysis of multidimensional NMR spectra recorded for complexes containing 15N-labeled protein. Simulated annealing calculations produced a family of 22 structures fully consistent with the constraints. Several intermolecular protein-ligand NOEs were obtained by using a novel approach monitoring temperature effects of NOE signals resulting from dynamic processes in the bound ligand. At low temperature (5 degrees C) the trimethoxy ring of bound trimetrexate is flipping sufficiently slowly to give narrow signals in slow exchange, which give good NOE cross peaks. At higher temperature these broaden and their NOE cross peaks disappear thus allowing the signals in the lower-temperature spectrum to be identified as NOEs involving ligand protons. The binding site for trimetrexate is well defined and this was compared with the binding sites in related complexes formed with methotrexate and trimethoprim. No major conformational differences were detected between the different complexes. The 2,4-diaminopyrimidine-containing moieties in the three drugs bind essentially in the same binding pocket and the remaining parts of their molecules adapt their conformations such that they can make effective van der Waals interactions with essentially the same set of hydrophobic amino acids, the side-chain orientations and local conformations of which are not greatly changed in the different complexes (similar chi1 and chi2 values).
Collapse
Affiliation(s)
- V I Polshakov
- Division of Molecular Structure, National Institute for Medical Research, The Ridgeway, Mill Hill, London, United Kingdom
| | | | | | | | | |
Collapse
|
8
|
Morgan WD, Birdsall B, Nieto PM, Gargaro AR, Feeney J. 1H/15N HSQC NMR studies of ligand carboxylate group interactions with arginine residues in complexes of brodimoprim analogues and Lactobacillus casei dihydrofolate reductase. Biochemistry 1999; 38:2127-34. [PMID: 10026296 DOI: 10.1021/bi982359u] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
1H and 15N NMR studies have been undertaken on complexes of Lactobacillus casei dihydrofolate reductase (DHFR) formed with analogues of the antibacterial drug brodimoprim (2,4-diamino-5-(3', 5'-dimethoxy-4'-bromobenzyl)pyrimidine) in order to monitor interactions between carboxylate groups on the ligands and basic residues in the protein. These analogues had been designed by computer modeling with carboxylated alkyl chains introduced at the 3'-O position in order to improve their binding properties by making additional interactions with basic groups in the protein. Specific interactions between ligand carboxylate groups and the conserved Arg57 residue have been detected in studies of 1H/15N HSQC spectra of complexes of DHFR with both the 4-carboxylate and the 4, 6-dicarboxylate brodimoprim analogues. The spectra from both complexes showed four resolved signals for the four NHeta protons of the guanidino group of Arg57, and this is consistent with hindered rotation in the guanidino group resulting from interactions with the 4-carboxylate group in each analogue. In the spectra of each complex, one of the protons from each of the two NH2 groups and both nitrogens are considerably deshielded compared to the shielding values normally observed for such nuclei. This pattern of deshielding is that expected for a symmetrical end-on interaction of the carboxylate oxygens with the NHeta12 and NHeta22 guanidino protons. The differences in the degree of deshielding between the complexes of the two structurally similar brodimoprim analogues and the methotrexate indicates that the shielding is very sensitive to geometry, most probably to hydrogen bond lengths. The 1H/15N HSQC spectrum of the DHFR complex with the brodimoprim-6-carboxylate analogue does not feature any deshielded Arg NHeta protons and this argues against a similar interaction with the Arg57 in this case. It has not proved possible to determine whether the 6-carboxylate in this analogue is interacting directly with any residue in the protein. 1H/15N HSQC spectra have been fully assigned for the complexes with the three brodimoprim analogues and chemical shift mapping used to explore interactions in the binding site. The 1H signals of the bound ligands for all three brodimoprim analogues have been assigned. Their 1H chemical shifts were found to be fairly similar in the different complexes indicating that the 2, 4-diaminopyrimidine and the benzyl ring are binding in essentially the same binding sites and with the same overall conformation in the different complexes. The rotation rate about the NepsilonCzeta bond in the brodimoprim-4,6-dicarboxylate complex with DHFR has been determined from a zz-HSQC exchange experiment, and its value is quite similar to that observed in the DHFR.methotrexate complex (24 +/- 10 s-1 at 8 degrees C and 50 +/- 10 s-1 at 15 degrees C, respectively). The 1H and 15N chemical shift differences of selected amide and guanidino NH groups, measured between the DHFR complexes, provided further evidence about the interactions involving Arg57 with the 4-carboxylate and 4,6-dicarboxylate brodimoprim analogues.
Collapse
Affiliation(s)
- W D Morgan
- Molecular Structure Division, National Institute for Medical Research, London, U.K
| | | | | | | | | |
Collapse
|
9
|
Jansson M, Andersson G, Uhlén M, Nilsson B, Kördel J. The insulin-like growth factor (IGF)binding protein 1 binding epitope on IGF-I probed by heteronuclear NMR spectroscopy and mutational analysis. J Biol Chem 1998; 273:24701-7. [PMID: 9733769 DOI: 10.1074/jbc.273.38.24701] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
NMR spectroscopy studies and biosensor interaction analysis of native and site-directed mutants of insulin-like growth factor I (IGF-I) was applied to identify the involvement of individual residues in IGF-I binding to IGF-binding protein 1 (IGFBP-1). Backbone NMR chemical shifts were found to be affected by IGFBP-1 binding in the following residues: Pro2, Glu3, Cys6, Gly7, Gly19, Pro28-Gly30, Gly32, Arg36, Arg37, Gln40-Gly42, Pro63, Lys65, Pro66, and Lys68-Ala70. Three IGF-I arginine side chains were identified by NMR to participate in IGFBP-1 binding. All IGF-I arginine residues were replaced by alanines, using site-directed mutagenesis, in four single substituted variants, IGF-I(R21A), IGF-I(R50A), IGF-I(R55A), and IGF-I(R56A), and one double replacement mutant, IGF-I(R36A/R37A). Biosensor interaction analysis binding studies demonstrate the involvement of Arg36-Arg37 and Arg50 in IGFBP-1 binding, while experiments with the IGF-I receptor implicate Arg21, Arg36-Arg37, and Arg56 as part of the receptor binding epitope. These overlapping binding surfaces explain why IGF-I receptor and IGFBP-1 binding to IGF-I is competitive. The C terminus of free, but not IGFBP-1-bound, IGF-I is found to exist in two distinct, NMR-detectable conformations at 30 degreesC. One possible explanation for this structural heterogeneity could be cis-trans isomerization of the Cys6-Cys48 disulfide bond.
Collapse
Affiliation(s)
- M Jansson
- Department of Structural Chemistry, Pharmacia & Upjohn, SE-11287 Stockholm, Sweden
| | | | | | | | | |
Collapse
|
10
|
Dardel F, Ragusa S, Lazennec C, Blanquet S, Meinnel T. Solution structure of nickel-peptide deformylase. J Mol Biol 1998; 280:501-13. [PMID: 9665852 DOI: 10.1006/jmbi.1998.1882] [Citation(s) in RCA: 70] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
In the accompanying paper, we report that zinc is unlikely to be the co-factor supporting peptide deformylase activity in vivo. In contrast, nickel binding promotes full enzyme activity. The three-dimensional structure of the resulting nickel-containing peptide deformylase (catalytic domain, residues 1 to 147) was solved by NMR using a 13C-15N-doubly labelled protein sample. A set of 2261 restraints could be collected, with an average of 15.4 per amino acid. The resolution, which shows a good definition for the position of most side-chains, is greatly improved compared to that previously reported for the zinc-containing, inactive form. A comparison of the two stuctures indicates however that both share the same 3D organization. This shows that the nature of the bound metal is the primary determinant of the hydrolytic activity of this enzyme. Site-directed mutagenesis enabled us to determine the conserved residues of PDF involved in the structure of the active site. In particular, a buried arginine appears to be critical for the positioning of Cys90, one of the metal ligands. Furthermore, the 3D structure of peptide deformylase was compared to thermolysin and metzincins. Although the structural folds are very different, they all display a common structural motif involving an alpha-helix and a three-stranded beta-sheet. These conserved structural elements build a common scaffold which includes the active site, suggesting a common hydrolytic mechanism for these proteases. Finally, an invariant glycine shared by both PDF and metzincins enables us to extend the conserved motif from HEXXH to HEXXHXXG.
Collapse
Affiliation(s)
- F Dardel
- Unité Mixte de Recherche, Ecole Polytechnique, Palaiseau cedex, F-91128, France
| | | | | | | | | |
Collapse
|
11
|
Gargaro AR, Soteriou A, Frenkiel TA, Bauer CJ, Birdsall B, Polshakov VI, Barsukov IL, Roberts GC, Feeney J. The solution structure of the complex of Lactobacillus casei dihydrofolate reductase with methotrexate. J Mol Biol 1998; 277:119-34. [PMID: 9514736 DOI: 10.1006/jmbi.1997.1560] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
We have determined the three-dimensional solution structure of the complex of Lactobacillus casei dihydrofolate reductase (18.3 kDa, 162 amino acid residues) formed with the anticancer drug methotrexate using 2531 distance, 361 dihedral angle and 48 hydrogen bond restraints obtained from analysis of multidimensional NMR spectra. Simulated annealing calculations produced a family of 21 structures fully consistent with the constraints. The structure has four alpha-helices and eight beta-strands with two other regions, comprising residues 11 to 14 and 126 to 127, also interacting with each other in a beta-sheet manner. The methotrexate binding site is very well defined and the structure around its glutamate moiety was improved by including restraints reflecting the previously determined specific interactions between the glutamate alpha-carboxylate group with Arg57 and the gamma-carboxylate group with His28. The overall fold of the binary complex in solution is very similar to that observed in the X-ray studies of the ternary complex of L. casei dihydrofolate reductase formed with methotrexate and NADPH (the structures of the binary and ternary complexes have a root-mean-square difference over the backbone atoms of 0.97 A). Thus no major conformational change takes place when NADPH binds to the binary complex. In the binary complex, the loop comprising residues 9 to 23 which forms part of the active site has been shown to be in the "closed" conformation as defined by M. R. Sawaya & J. Kraut, who considered the corresponding loops in crystal structures of complexes of dihydrofolate reductases from several organisms. Thus the absence of the NADPH does not result in the "occluded" form of the loop as seen in crystal studies of some other dihydrofolate reductases in the absence of coenzyme. Some regions of the structure in the binary complex which form interaction sites for NADPH are less well defined than other regions. However, in general terms, the NADPH binding site appears to be essentially pre-formed in the binary complex. This may contribute to the tighter binding of coenzyme in the presence of methotrexate.
Collapse
Affiliation(s)
- A R Gargaro
- Division of Molecular Structure, National Institute for Medical Research, London, UK
| | | | | | | | | | | | | | | | | |
Collapse
|
12
|
Polshakov VI, Williams MA, Gargaro AR, Frenkiel TA, Westley BR, Chadwick MP, May FE, Feeney J. High-resolution solution structure of human pNR-2/pS2: a single trefoil motif protein. J Mol Biol 1997; 267:418-32. [PMID: 9096235 DOI: 10.1006/jmbi.1997.0896] [Citation(s) in RCA: 47] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
pNR-2/pS2 is a 60 residue extracellular protein, which was originally discovered in human breast cancer cells, and subsequently found in other tumours and normal gastric epithelial cells. We have determined the three-dimensional solution structure of a C58S mutant of human pNR-2/pS2 using 639 distance and 137 torsion angle constraints obtained from analysis of multidimensional NMR spectra. A series of simulated annealing calculations resulted in the unambiguous determination of the protein's disulphide bonding pattern and produced a family of 19 structures consistent with the constraints. The peptide contains a single "trefoil" sequence motif, a region of about 40 residues with a characteristic sequence pattern, which has been found, either singly or as a repeat, in about a dozen extracellular proteins. The trefoil domain contains three disulphide bonds, whose 1-5, 2-4 and 3-6 cysteine pairings form the structure into three closely packed loops with only a small amount of secondary structure, which consists of a short alpha-helix packed against a two-stranded antiparallel beta-sheet. The structure of the domain is very similar to those of the two trefoil domains that occur in porcine spasmolytic polypeptide (PSP), the only member of the trefoil family whose three-dimensional structure has been previously determined. Outside the trefoil domain, which forms the compact "head" of the molecule, the N and C-terminal strands are closely associated, forming an extended "tail", which has some beta-sheet character for part of its length and which becomes more disordered towards the termini as indicated by (15)N{(1)H} NOEs. We have considered the structural implications of the possible formation of a native C58-C58 disulphide-bonded homodimer. Comparison of the surface features of pNR-2/pS2 and PSP, and consideration of the sequences of the other human trefoil domains in the light of these structures, illuminates the possible role of specific residues in ligand/receptor binding.
Collapse
Affiliation(s)
- V I Polshakov
- Laboratory of Physical Methods, Chemical-Pharmaceutical Research Institute, Moscow, Russia
| | | | | | | | | | | | | | | |
Collapse
|
13
|
Verma CS, Caves LS, Hubbard RE, Roberts GC. Domain motions in dihydrofolate reductase: a molecular dynamics study. J Mol Biol 1997; 266:776-96. [PMID: 9102469 DOI: 10.1006/jmbi.1996.0818] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Molecular dynamics simulations have been carried out on the enzyme dihydrofolate reductase from Lactobacillus casei complexed with methotrexate, NADPH and 264 crystallographic water molecules. Analysis of correlations in atomic fluctuations reveal the presence of highly correlated motion (correlation coefficient > 0.6) in the region between residues 30 to 35 and 85 to 90 leading to the identification of two domains, an "adenosine-binding domain" and a "large domain", which rotate by 3 to 4 degrees with respect to each other. The strongest correlation (> 0.6) within the large domain involves a coupling between the motions of the "teen-loop", and the spatially contiguous loops linking beta 6-beta 7 and beta 7-beta 8. Moreover, there is a significant correlation (approximately 0.5) between the adenosine fragment of NADPH and the pteridine and p-aminobenzoyl fragments of methotrexate, which are separated by approximately 17 A, and is lost on removal of "rigid-body" motion from the original trajectory. This provides support for the idea that the relative motion of the two domains is a means by which the occupation of the binding site for the adenosine end of the coenzyme can affect methotrexate binding and vice versa. Quasiharmonic vibrational analysis of the trajectory reveals that the overall dynamics of the system are governed by domain motions whose contributions are dominant at low frequencies. In addition, different low-frequency modes are responsible for separately coupling the adenosine-binding site and parts of methotrexate.
Collapse
Affiliation(s)
- C S Verma
- Department of Chemistry, University of York, Heslington, UK
| | | | | | | |
Collapse
|
14
|
|