Chemically induced contraction and stretching of a linear rotaxane dimer

Copper(I)-induced assembly of two self-complementary identical units, which consist of a ring that incorporates a 1,10-phenanthroline group attached to a small filament containing a second 1,10-phenanthroline (phen) group, leads quantitatively to a doubly threaded complex. Each copper(I) center is four-coordinate and is located inside a ring and bound to a phen from the macrocyle. The two other coordination sites are occupied by a phen from the filament connected to the other ring. An X-ray structure of the dicopper(I) complex unambiguously demonstrates the doubly threaded nature of the system. The molecule has C(2) symmetry in the crystal. This is an extended form with a Cu small middle dot small middle dot small middle dotCu separation of 18.3 A and an overall length close to 40 A. Further synthetic work, which utilizes the two terminal phenolic functions of the previous dicopper(I) complex, gives rise to a more complex system in which both filaments have been prolonged in opposite directions by 2,2':6',2"-terpyridine (terpy) motifs and bulky stoppers. The organic backbone is that of a rotaxane dimer. Although redox cycling of Cu(I) to Cu(II) did not lead to intramolecular rearrangement, simple chemical reactions induced large conformational changes. The rotaxane dimer was set in motion as follows. The dicopper(I) complex, which is in an extended conformation, was demetallated by using KCN. From the free ligand, the dizinc complex was formed quantitatively at room temperature. (1)H NMR data show that a new conformation is obtained: each Zn(II) is five-coordinate (phen + terpy), and the molecule is in a contracted conformation. This process is reminiscent of biological muscles in the sense that the two filaments of this system can be moved along one another in a gliding motion that keeps the whole system together, but which converts a stretched compound (overall length approximately equal to 83 A) into a contracted species (overall length approximately equal to 65 A, according to CPK models). The motion is quantitatively reversed by the addition of an excess of copper(I) to the dizinc complex; this regenerates the extended starting form. Although the motivation of the present contribution was to illustrate that a muscle-like molecule may be stretched or contracted using electrochemistry and coordination chemistry, the main body of the work is organic synthesis. This is testified by the fact that the dicopper(I) rotaxane dimer was obtained in 23 steps from commercially available compounds.     

Light activates reduction of methotrexate by NADPH in the ternary complex with Escherichia coli dihydrofolate reductase

Methotrexate (MTX), a strong inhibitor of dihydrofolate reductase (DHFR), has been widely used for chemotherapy for many types of cancer as well as for juvenile rheumatoid arthritis. It mimics folate substrates and binds tightly to the active site of DHFR, perhaps in a conformation close to the transition state of the folate catalyzed reaction. Absorption, fluorescence and ultrasensitive Raman difference spectroscopies show that light-activated MTX reacts with NADPH in the enzyme active site, producing 5,8-dihydromethotrexate (5,8-dihydro-MTX) and NADP+. The reaction, which proceeds with a hydride transfer between C4 (pro-R side) of the nicotinamide ring and N5 of the pteridine ring, is similar to that between folate and NADPH except that the hydride is transferred to C6 in this case. Hence, MTX is catalytically competent in its excited state. Most experiments were performed on the Escherichia coli enzyme, but preliminary studies show that the reaction also occurs with human DHFR.     

How dihydrofolate reductase facilitates protonation of dihydrofolate

Dihydrofolate Reductase (DHFR) catalyzes the reduction of dihydrofolate (H2F) to tetrahydrofolate. On the basis of 10-12.5 ns molecular dynamics simulations of two conformations (closed and occluded) of the ternary DHFR/NADPH/H2F complex from Escherichia coli and a free energy perturbation approach, we have calculated the pKa value for the N5 atom in H2F. Our results suggest that the N5 atom in H2F is responsible for the pH dependency of the catalyzed reaction, meaning that DHFR facilitates protonation of H2F by approximately 4 pKa units. The mechanism behind this increase is due to favorable electrostatic interactions between the Asp27 residue and a proton at the N5 atom. The electrostatic interactions are enhanced by a hydrophobic active site, which to a large extent is made hydrophobic by the M20 loop in DHFR. Moreover, we find that the conformation imposed on H2F by DHFR to some extent also favors protonation of the N5 atom. Our results add support to previous findings and suggestions by Callender and co-workers [e.g., Deng, J.; Callender, R. J. Am. Chem. Soc. 1998, 120, 7730-7737] and explain why mutation of Asp27 may lead to severely reduced activity at neutral pH.     

Isotachophoretic analysis of the dihydrofolate reductase reaction in the presence of methotrexate and ascorbic acid

The antifolate methotrexate (MTX) is widely used in cancer chemotherapy. In this study, we show that MTX (MTX-Glu1) and MTX-polyglutamates (MTX-Glu2-5) strongly inhibited the growth of the leukemic cell line MOLT-4. This effect, however, was mitigated by ascorbic acid. We investigated whether ascorbic acid is able to reduce dihydrofolic acid (DHF) to tetrahydrofolic acid (THF) directly or by circumventing the MTX inhibition of dihydrofolate reductase (DHFR). The inhibition of this NADPH-dependent reduction of DHF by MTX-Glun in the absence or presence of ascorbate, was determined by analytical isotachophoresis. Using 0.01 M HCl/histidine, pH 6.0, as a leading electrolyte (L) and 0.005 M 2-(N-morpholino)ethanesulfonic acid (MES)/histidine, pH 6.0, as a terminating electrolyte (T), MTX-Glun derivatives including MTX-Glu1 could be easily separated, whereas the quantitative estimation of THF was not possible. A quantitative characterization of the DHFR reaction by measuring NADPH, NADP+ and ascorbate was achieved with another system (L: 0.01 M HCI/beta-alanine, pH 3.73; T: 0.01 M caproic acid, pH 3.27). Nanomolar concentrations of MTX-Glu1-5 inhibited consumption of NADPH and production of NADP+. Ascorbic acid was not able to reduce DHF, neither directly nor after inhibition of DHFR by MTX. However, ascorbic acid seemed to diminish the oxidation of THF and this may account for its capacity to reduce the inhibitory effect of MTX on MOLT-4 cells.     

高效毛细管电泳法测定二氢叶酸还原酶活力的研究

通过胶束电动毛细管电泳法研究分离二氢叶酸还原酶体系中二氢叶酸、四氢叶酸、 NADP、 NADPH和酶5种组分,在含0.002%Brij-35的pH 9.18 50 mmol/L 的硼砂缓冲溶液中,5种组分在18min内得到基线分离.通过对其产物四氢叶酸峰面积的定量测定,计算出二氢叶酸还原酶的米氏常数,建立了毛细管电泳法对二氢叶酸还原酶活力的测定方法.     

Freezing a single distal motion in dihydrofolate reductase

Constraining a single motion between distal residues separated by approximately 28 A in hybrid quantum/classical molecular dynamics simulations is found to increase the free energy barrier for hydride transfer in dihydrofolate reductase by approximately 3 kcal/mol. Our analysis indicates that a single distal constraint alters equilibrium motions throughout the enzyme on a wide range of time scales. This alteration of the conformational sampling of the entire system is sufficient to significantly increase the free energy barrier and decrease the rate of hydride transfer. Despite the changes in conformational sampling introduced by the constraint, the system assumes a similar transition state conformation with a donor-acceptor distance of approximately 2.72 A to enable the hydride transfer reaction. The modified thermal sampling leads to a substantial increase in the average donor-acceptor distance for the reactant state, however, thereby decreasing the probability of sampling the transition state conformations with the shorter distances required for hydride transfer. These simulations indicate that fast thermal fluctuations of the enzyme, substrate, and cofactor lead to conformational sampling of configurations that facilitate hydride transfer. The fast thermal motions are in equilibrium as the reaction progresses along the collective reaction coordinate, and the overall average equilibrium conformational changes occur on the slower time scale measured experimentally. Recent single molecule experiments suggest that at least some of these thermally averaged equilibrium conformational changes occur on the millisecond time scale of the hydride transfer reaction. Thus, introducing a constraint that modifies the conformational sampling of an enzyme could significantly impact its catalytic activity.     

The role of large-scale motions in catalysis by dihydrofolate reductase

Dihydrofolate reductase has long been used as a model system to study the coupling of protein motions to enzymatic hydride transfer. By studying environmental effects on hydride transfer in dihydrofolate reductase (DHFR) from the cold-adapted bacterium Moritella profunda (MpDHFR) and comparing the flexibility of this enzyme to that of DHFR from Escherichia coli (EcDHFR), we demonstrate that factors that affect large-scale (i.e., long-range, but not necessarily large amplitude) protein motions have no effect on the kinetic isotope effect on hydride transfer or its temperature dependence, although the rates of the catalyzed reaction are affected. Hydrogen/deuterium exchange studies by NMR-spectroscopy show that MpDHFR is a more flexible enzyme than EcDHFR. NMR experiments with EcDHFR in the presence of cosolvents suggest differences in the conformational ensemble of the enzyme. The fact that enzymes from different environmental niches and with different flexibilities display the same behavior of the kinetic isotope effect on hydride transfer strongly suggests that, while protein motions are important to generate the reaction ready conformation, an optimal conformation with the correct electrostatics and geometry for the reaction to occur, they do not influence the nature of the chemical step itself; large-scale motions do not couple directly to hydride transfer proper in DHFR.     

Radioenzymatic determination of subpicomole quantities of folic acid and dihydrofolate based on differential reduction by dihydrofolate reductase

A sensitive method for the determination of folic acid in the presence of dihydrofolate has been described. The analysis is based on the difference in reactivity of dihydrofolate reductase toward folic acid and dihydrofolate. The tetrahydrofolate formed was measured by stoichiometric entrapment into a covalent, ternary complex with thymidylate synthase and [³H]fluorodeoxyuridylate. Differentiation between folic acid and dihydrofolate was accomplished by the use of two levels of dihydrofolate reductase. Recovery of both folic acid and dihydrofolate was quantitative in the range of 0.1 to 2.0 pmol.     

Carbon-deuterium bonds as probes of dihydrofolate reductase

Much effort has been directed toward understanding the contributions of electrostatics and dynamics to protein function and especially to enzyme catalysis. Unfortunately, these studies have been limited by the absence of direct experimental probes. We have been developing the use of carbon-deuterium bonds as probes of proteins and now report the application of the technique to the enzyme dihydrofolate reductase, which catalyzes a hydride transfer and has served as a paradigm for biological catalysis. We observe that the stretching absorption frequency of (methyl- d 3) methionine carbon-deuterium bonds shows an approximately linear dependence on solvent dielectric. Solvent and computational studies support the empirical interpretation of the stretching frequency in terms of local polarity. To begin to explore the use of this technique to study enzyme function and mechanism, we report a preliminary analysis of (methyl- d 3) methionine residues within dihydrofolate reductase. Specifically, we characterize the IR absorptions at Met16 and Met20, within the catalytically important Met20 loop, and Met42, which is located within the hydrophobic core of the enzyme. The results confirm the sensitivity of the carbon-deuterium bonds to their local protein environment, demonstrate that dihydrofolate reductase is electrostatically and dynamically heterogeneous, and lay the foundation for the direct characterization protein electrostatics and dynamics and, potentially, their contribution to catalysis.     

A molecular dynamics study of the inhibition of chicken dihydrofolate reductase by a phenyl triazine

Molecular dynamics simulations have been used to investigate the ternary complex formed between chicken liver dihydrofolate reductase, a phenyl triazine inhibitor, and reduced nicotinamide adenine dinucleotide phosphate (NADPH). The solvent was represented by a sphere of water molecules encompassing the system. We report the results of quantum mechanical calculations of the rotational barrier in the pyrophosphate link and the barrier to inversion of the triazine ring. AMBER parameters for NADPH and the triazine are provided. Over the course of a 300-ps molecular dynamics simulation of the ternary complex in water, the triazine inhibitor maintains the same hydrogen bonding and hydrophobic interactions with the enzyme that are observed in the X-ray crystal structure. Despite the low calculated barrier to inversion of the triazine ring, a single puckered conformation is observed throughout the simulation. It is proposed that this is primarily due to interactions with Phe34, which maintains an approximately parallel orientation to the triazine ring. The nicotinamide portion of NADPH maintains the interactions observed in the crystal structure, but more conformational change is observed at the adenine end together with associated changes in the protein. Two conformations for the sidechain of Tyr31 are present in the X-ray structure. The main simulation reported here corresponds to the conformation characterized by (χ₁ = ? 161°, χ₂ = ? 103°). A separate simulation was also performed in which the sidechain of Tyr31 was initially set to the other conformation present in the crystal structure (χ₁ = 139°, χ₂ = ?99°). During this simulation, χ₁ of this sidechain gradually changed until it occupied the region characterized by χ₁ = ?160°, thereby suggesting that this is the preferred conformation for this residue. The simulation required 200 ps to reach structural equilibrium (as measured by the root mean square, rms, deviation from the initial crystal structure), thus reinforcing the view that simulations of at least several hundreds of picoseconds are desirable when studying such systems. © 1995 John Wiley & Sons, Inc.     

Identification of an anti-MRSA dihydrofolate reductase inhibitor from a diversity-oriented synthesis

The screening of a diversity-oriented synthesis library followed by structure-activity relationship investigations have led to the discovery of an anti-MRSA agent which operates as an inhibitor of Staphylococcus aureus dihydrofolate reductase.     

Synthesis of a neamine dimer targeting the dimerization initiation site of HIV-1 RNA

A neamine dimer designed to bind to a specific sequence of HIV-1 RNA has been synthesized. Starting from neomycin B (1), a five-step synthesis efficiently provided a key protected neamine monomer 6 (28%). From the latter, coupling reactions with activated diacids gave dimers. After deprotection, a neamine dimer was obtained as the hexachlorohydrate salt 15 with 13% overall yield over nine steps.     

Effect of mutation on enzyme motion in dihydrofolate reductase

Hybrid quantum-classical molecular dynamics simulations of a mutant Escherichia coli dihydrofolate reductase enzyme are presented. Although residue 121 is on the exterior of the enzyme, experimental studies have shown that the mutation of Gly-121 to valine reduces the rate of hydride transfer by a factor of 163. The simulations indicate that the decrease in the hydride transfer rate for the G121V mutant is due to an increase in the free energy barrier. The calculated free energy barrier is higher for the mutant than for the wild-type enzyme by an amount that is consistent with the experimentally observed rate reduction. The calculated transmission coefficients are comparable for the wild-type and mutant enzymes. The simulations suggest that this mutation may interrupt a network of coupled promoting motions proposed to play an important role in DHFR catalysis. This phenomenon has broad implications for protein engineering and drug design.     

Deformation of copolymer micelles induced by amphiphilic dimer particles

Combining self-consistent-field theory and density-functional theory,we systematically study the deformation of copolymer micelles induced by the presence of amphiphilic dimer particles.Due to the amphiphilic nature,dimer particles tend to accumulate onto the interface of the copolymer micelle.With increasing concentration of the symmetric dimer particles,which are made of two identical spherical particles,the micelle deforms from the initial sphere to ellipse,dumbbell, and finally separates into two micelles.Furthermore,asymmetric dimer particles,composed by two particles with different sizes,are considered to investigate the influence of geometry of dimer particles on the deformation of the micelle.It is found that the micelle inclines to deform into dumbbell due to the additional curvature originating in the gathering of asymmetric dimer particles onto the interface of the micelle.The present study on the deformation of micelles is useful to understand the possible shape variation in the course of cell division/fusion.     

A reverse turn structure induced by a D,L-alpha-aminoxy acid dimer

Our previous work revealed that two adjacent D-alpha-aminoxy acids could form two homochiral N-O turns, with the backbone folding into an extended helical structure (1.8(8)-helix). Here, we report the conformational studies of linear peptides 3-6, which contain a D,L-alpha-aminoxy acid dimer segment. The NMR and X-ray analysis of 3 showed that it folded into a loop conformation with two heterochiral N-O turns. This loop segment can be used to constrain tetrapeptides 4 and 6 to form a reverse turn structure. (1)H NMR dilution studies, DMSO-d6 addition studies, and 2D-NOESY data indicated that tetrapeptides 4 and 6 folded into reverse turn conformations featured by a head-to-tail 16-membered-ring intramolecular hydrogen bond. In contrast, tetrapeptide 5 with L-Ala instead of Gly or D-Ala as the N-terminal amino acid could not form the desired reverse turn structure for steric reasons. Quantum mechanics calculations showed that model pentamide 7, with the same substitution pattern of 4, adopted a novel reverse turn conformation featuring two heterochiral N-O turns (each of an 8-membered ring hydrogen bond), a cross-strand 16-membered ring hydrogen bond, and a 7-membered ring gamma-turn.     

Three-dimensional quantitative structure-activity and structure-selectivity relationships of dihydrofolate reductase inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) modelling using comparative molecular similarity indices analysis (CoMSIA) was applied to a series of 406 structurally diverse dihydrofolate reductase (DHFR) inhibitors from Pneumocystis carinii (pc) and rat liver (rl). X-ray crystal structures of three inhibitors bound to pcDHFR were used for defining the alignment rule. For pcDHFR, a QSAR model containing 6 components was selected using leave-10%-out cross-validation (n= 240, q2 = 0.65), while a 4-component model was selected for rlDHFR (n= 237, q2 = 0.63); both include steric, electrostatic and hydrophobic contributions. The models were validated using a large test set, designed to maximise its diversity and to verify the predictive accuracy of models for extrapolation. The pcDHFR model has r2 = 0.60 and mean absolute error (MAE) = 0.57 for the test set after removing 4 outliers, and the rlDHFR model has r2 = 0.60 and MAE = 0.69 after removing 4 test set outliers. In addition, classification models predicting selectivity for pcDHFR over rlDHFR were developed using soft independent modelling by class analogy (SIMCA), with a selectivity ratio of 2 (IC50,rlDHFR/ IC50,pcDHFR) used for delimiting classes. A 5-component model including steric and electrostatic contributions has cross-validated and test set classification rates of 0.67 and 0.68 for selective inhibitors, and 0.85 and 0.72 for unselective inhibitors. The predictive accuracy of models, together with the identification of important contributions in QSAR and classification models, offer the possibility of designing potent selective inhibitors and estimating their activity prior to synthesis.     

Inhibition of human dihydrofolate reductase by 2,4-diaminoquinazolines bearing simple substituents on the aromatic ring

A series of thirty eight 2,4-diaminoquinazolines having diverse substitution patterns on the aromatic ring was evaluated for inhibitory activity against dihydrofolate reductase (DHFR) obtained from a human lymphoblast cell line. Many of these compounds were also evaluated as inhibitors of rat liver DHFR under the same experimental conditions. In most instances the results obtained with each enzyme were comparable indicating that the rodent enzyme is a suitable model for the human DHFR as far as the determination of I₅₀ values is concerned. The results demonstrate that relatively simple 5-substituted- or 5,6-disubstituted-2,4-diaminoquinazolines can be potent DHFR inhibitors. The presence of a nonpolar substituent at position 7 or 8 was highly detrimental to inhibitory potency.