Dynamic dysfunction in dihydrofolate reductase results from antifolate drug binding: modulation of dynamics within a structural state

The mechanistic insights into different aspects of promiscuity in metalloenzymes

Enzymes are nature's ultimate machinery to catalyze complex reactions. Though enzymes are evolved to catalyze specific reactions, they also show significant promiscuity in reactions and substrate selection. Metalloenzymes contain a metal ion or metal cofactor in their active site, which is crucial in their catalytic activity. Depending on the metal and its coordination environment, the metal ion or cofactor may function as a Lewis acid or base and a redox center and thus can catalyze a plethora of natural reactions. In fact, the versatility in the oxidation state of the metal ions provides metalloenzymes with a high level of catalytic adaptability and promiscuity. In this chapter, we discuss different aspects of promiscuity in metalloenzymes by using several recent experimental and theoretical works as case studies. We start our discussion by introducing the concept of promiscuity and then we delve into the mechanistic insight into promiscuity at the molecular level.

Allosteric regulatory control in dihydrofolate reductase is revealed by dynamic asymmetry

We investigated the relationship between mutations and dynamics in Escherichia coli dihydrofolate reductase (DHFR) using computational methods. Our study focused on the M20 and FG loops, which are known to be functionally important and affected by mutations distal to the loops. We used molecular dynamics simulations and developed position-specific metrics, including the dynamic flexibility index (DFI) and dynamic coupling index (DCI), to analyze the dynamics of wild-type DHFR and compared our results with existing deep mutational scanning data. Our analysis showed a statistically significant association between DFI and mutational tolerance of the DHFR positions, indicating that DFI can predict functionally beneficial or detrimental substitutions. We also applied an asymmetric version of our DCI metric (DCIasym ) to DHFR and found that certain distal residues control the dynamics of the M20 and FG loops, whereas others are controlled by them. Residues that are suggested to control the M20 and FG loops by our DCIasym metric are evolutionarily nonconserved; mutations at these sites can enhance enzyme activity. On the other hand, residues controlled by the loops are mostly deleterious to function when mutated and are also evolutionary conserved. Our results suggest that dynamics-based metrics can identify residues that explain the relationship between mutation and protein function or can be targeted to rationally engineer enzymes with enhanced activity.

DHFR Mutants Modulate Their Synchronized Dynamics with the Substrate by Shifting Hydrogen Bond Occupancies

Antibiotic resistance is a global health problem in which mutations occurring in functional proteins render drugs ineffective. The working mechanisms of the arising mutants are seldom apparent; a methodology to decipher these mechanisms systematically would render devising therapies to control the arising mutational pathways possible. Here we utilize C_α-C_β bond vector relaxations obtained from moderate length MD trajectories to determine conduits for functionality of the resistance conferring mutants of Escherichia coli dihydrofolate reductase. We find that the whole enzyme is synchronized to the motions of the substrate, irrespective of the mutation introducing gain-of-function or loss-of function. The total coordination of the motions suggests changes in the hydrogen bond dynamics with respect to the wild type as a possible route to determine and classify the mode-of-action of individual mutants. As a result, nine trimethoprim-resistant point mutations arising frequently in evolution experiments are categorized. One group of mutants that display the largest occurrence (L28R, W30G) work directly by modifying the dihydrofolate binding region. Conversely, W30R works indirectly by the formation of the E139-R30 salt bridge which releases energy resulting from tight binding by distorting the binding cavity. A third group (D27E, F153S, I94L) arising as single, resistance invoking mutants in evolution experiment trajectories allosterically and dynamically affects a hydrogen bonding motif formed at residues 59-69-71 which in turn modifies the binding site dynamics. The final group (I5F, A26T, R98P) consists of those mutants that have properties most similar to the wild type; these only appear after one of the other mutants is fixed on the protein structure and therefore display clear epistasis. Thus, we show that the binding event is governed by the entire enzyme dynamics while the binding site residues play gating roles. The adjustments made in the total enzyme in response to point mutations are what make quantifying and pinpointing their effect a hard problem. Here, we show that hydrogen bond dynamics recorded on sub-μs time scales provide the necessary fingerprints to decipher the various mechanisms at play.

Structure-guided functional studies of plasmid-encoded dihydrofolate reductases reveal a common mechanism of trimethoprim resistance in Gram-negative pathogens

Two plasmid-encoded dihydrofolate reductase (DHFR) isoforms, DfrA1 and DfrA5, that give rise to high levels of resistance in Gram-negative bacteria were structurally and biochemically characterized to reveal the mechanism of TMP resistance and to support phylogenic groupings for drug development against antibiotic resistant pathogens. Preliminary screening of novel antifolates revealed related chemotypes that showed high levels of inhibitory potency against Escherichia coli chromosomal DHFR (EcDHFR), DfrA1, and DfrA5. Kinetics and biophysical analysis, coupled with crystal structures of trimethoprim bound to EcDHFR, DfrA1 and DfrA5, and two propargyl-linked antifolates (PLA) complexed with EcDHFR, DfrA1 and DfrA5, were determined to define structural features of the substrate binding pocket and guide synthesis of pan-DHFR inhibitors.

Critical residue variations in two of the most clinically prevalent DHFR isoforms are identified as a common structural element in trimethoprim-resistant DHFR which impose changes on enzyme catalysis and ligand-cofactor cooperativity.

The chemotherapeutic drug methotrexate selects for antibiotic resistance

BACKGROUND

Understanding drivers of antibiotic resistance evolution is fundamental for designing optimal treatment strategies and interventions to reduce the spread of antibiotic resistance. Various cytotoxic drugs used in cancer chemotherapy have antibacterial properties, but how bacterial populations are affected by these selective pressures is unknown. Here we test the hypothesis that the widely used cytotoxic drug methotrexate affects the evolution and selection of antibiotic resistance.

METHODS

First, we determined methotrexate susceptibility (IC90) and selective abilities in a collection of Escherichia coli and Klebsiella pneumoniae strains with and without pre-existing trimethoprim resistance determinants. We constructed fluorescently labelled pairs of E. coli MG1655 differing only in trimethoprim resistance determinants and determined the minimum selective concentrations of methotrexate using flow-cytometry. We further used an experimental evolution approach to investigate the effects of methotrexate on de novo trimethoprim resistance evolution.

FINDINGS

We show that methotrexate can select for acquired trimethoprim resistance determinants located on the chromosome or a plasmid. Additionally, methotrexate co-selects for genetically linked resistance determinants when present together with trimethoprim resistance on a multi-drug resistance plasmid. These selective effects occur at concentrations 40- to >320-fold below the methotrexate minimal inhibitory concentration.

INTERPRETATION

Our results strongly suggest a selective role of methotrexate for virtually any antibiotic resistance determinant when present together with trimethoprim resistance on a multi-drug resistance plasmid. The presented results may have significant implications for patient groups strongly depending on effective antibiotic treatment.

FUNDING

PJJ was supported by UiT The Arctic University of Norway and the Northern Norway Regional Health Authority (SFP1292-16/HNF1586-21) and JPI-EC-AMR (Project 271,176/H10). DIA was supported by the Swedish Research Council (grant 2017-01,527). The publication charges for this article have been funded by a grant from the publication fund of UiT The Arctic University of Norway.

Distal Regions Regulate Dihydrofolate Reductase-Ligand Interactions

Protein motions play a fundamental role in enzyme catalysis and ligand binding. The relationship between protein motion and function has been extensively investigated in the model enzyme dihydrofolate reductase (DHFR). DHFR is an essential enzyme that catalyzes the reduction of dihydrofolate to tetrahydrofolate. Numerous experimental and computational methods have been used to probe the motions of DHFR through the catalytic cycle and to investigate the effect of distal mutations on DHFR motions and ligand binding. These experimental investigations have pushed forward the study of protein motions and their role in protein-ligand interactions. The introduction of mutations distal to the active site has been shown to have profound effects on ligand binding, hydride transfer rates and catalytic efficacy and these changes are captured by enzyme kinetics measurements. Distal mutations have been shown to exert their effects through a network of correlated amino acids and these effects have been investigated by NMR, protein dynamics, and analysis of coupled amino acids. The experimental methods and the findings that are reviewed here have broad implications for our understanding of enzyme mechanisms, ligand binding and for the future design and discovery of enzyme inhibitors.

Monoclonal antibody stability can be usefully monitored using the excitation-energy-dependent fluorescence edge-shift

Among the major challenges in the development of biopharmaceuticals are structural heterogeneity and aggregation. The development of a successful therapeutic monoclonal antibody (mAb) requires both a highly active and also stable molecule. Whilst a range of experimental (biophysical) approaches exist to track changes in stability of proteins, routine prediction of stability remains challenging. The fluorescence red edge excitation shift (REES) phenomenon is sensitive to a range of changes in protein structure. Based on recent work, we have found that quantifying the REES effect is extremely sensitive to changes in protein conformational state and dynamics. Given the extreme sensitivity, potentially this tool could provide a ‘fingerprint’ of the structure and stability of a protein. Such a tool would be useful in the discovery and development of biopharamceuticals and so we have explored our hypothesis with a panel of therapeutic mAbs. We demonstrate that the quantified REES data show remarkable sensitivity, being able to discern between structurally identical antibodies and showing sensitivity to unfolding and aggregation. The approach works across a broad concentration range (µg–mg/ml) and is highly consistent. We show that the approach can be applied alongside traditional characterisation testing within the context of a forced degradation study (FDS). Most importantly, we demonstrate the approach is able to predict the stability of mAbs both in the short (hours), medium (days) and long-term (months). The quantified REES data will find immediate use in the biopharmaceutical industry in quality assurance, formulation and development. The approach benefits from low technical complexity, is rapid and uses instrumentation which exists in most biochemistry laboratories without modification.

The role of NMR in leveraging dynamics and entropy in drug design

Nuclear magnetic resonance (NMR) spectroscopy has contributed to structure-based drug development (SBDD) in a unique way compared to the other biophysical methods. The potency of a ligand binding to a protein is dictated by the binding free energy, which is an intricate interplay between entropy and enthalpy. In addition to providing the atomic resolution structural information, NMR can help to identify protein-ligand interactions that potentially contribute to the enthalpic component of the free energy. NMR can also illuminate dynamic aspects of the interaction, which correspond to the entropic term of the free energy. The ability of NMR to access both terms in the free energy equation stems from the suite of experiments developed to shed light on various aspects that contribute to both entropy and enthalpy, deepening our understanding of the biological function of macromolecules and assisting to target them in physiological conditions. Here we provide a brief account of the contribution of NMR to SBDD, highlighting hallmark examples and discussing the challenges that demand further method development. In the era of integrated biology, the unique ability of NMR to directly ascertain structural and dynamical aspects of macromolecule and monitor changes in these properties upon engaging a ligand can be combined with computational and other structural and biophysical methods to provide a more complete picture of the energetics of drug engagement with the target. Such efforts can be used to engineer better drugs.

Transitions in DNA polymerase β μs-ms dynamics related to substrate binding and catalysis

DNA polymerase β (pol β) plays a central role in the DNA base excision repair pathway and also serves as an important model polymerase. Dynamic characterization of pol β from methyl-TROSY 13C-1H multiple quantum CPMG relaxation dispersion experiments of Ile and Met sidechains and previous backbone relaxation dispersion measurements, reveals transitions in μs-ms dynamics in response to highly variable substrates. Recognition of a 1-nt-gapped DNA substrate is accompanied by significant backbone and sidechain motion in the lyase domain and the DNA binding subdomain of the polymerase domain, that may help to facilitate binding of the apoenzyme to the segments of the DNA upstream and downstream from the gap. Backbone μs-ms motion largely disappears after formation of the pol β-DNA complex, giving rise to an increase in uncoupled μs-ms sidechain motion throughout the enzyme. Formation of an abortive ternary complex using a non-hydrolyzable dNTP results in sidechain motions that fit to a single exchange process localized to the catalytic subdomain, suggesting that this motion may play a role in catalysis.

Conformational Dynamics and Cooperativity Drive the Specificity of a Protein-Ligand Interaction

Molecular recognition is critical for the fidelity of signal transduction in biology. Conversely, the disruption of protein-protein interactions can lead to disease. Thus, comprehension of the molecular determinants of specificity is essential for understanding normal biological signaling processes and for the development of precise therapeutics. Although high-resolution structures have provided atomic details of molecular interactions, much less is known about the influence of cooperativity and conformational dynamics. Here, we used the Tiam2 PSD-95/Dlg/ZO-1 (PDZ) domain and a quadruple mutant (QM), engineered by swapping the identity of four residues important for specificity in the Tiam1 PDZ into the Tiam2 PDZ domain, as a model system to investigate the role of cooperativity and dynamics in PDZ ligand specificity. Surprisingly, equilibrium binding experiments found that the ligand specificity of the Tiam2 QM was switched to that of the Tiam1 PDZ. NMR-based studies indicated that Tiam2 QM PDZ, but not other mutants, had extensive microsecond to millisecond motions distributed throughout the entire domain suggesting structural cooperativity between the mutated residues. Thermodynamic analyses revealed energetic cooperativity between residues in distinct specificity subpockets that was dependent upon the identity of the ligand, indicating a context-dependent binding mechanism. Finally, isothermal titration calorimetry experiments showed distinct entropic signatures along the mutational trajectory from the Tiam2 wild-type to the QM PDZ domain. Collectively, our studies provide unique insights into how structure, conformational dynamics, and thermodynamics combine to modulate ligand-binding specificity and have implications for the evolution, regulation, and design of protein-ligand interactions.

The crystal structure of a tetrahydrofolate-bound dihydrofolate reductase reveals the origin of slow product release

Dihydrofolate reductase (DHFR) catalyzes the stereospecific reduction of 7,8-dihydrofolate (FH2) to (6s)-5,6,7,8-tetrahydrofolate (FH4) via hydride transfer from NADPH. The consensus Escherichia coli DHFR mechanism involves conformational changes between closed and occluded states occurring during the rate-limiting product release step. Although the Protein Data Bank (PDB) contains over 250 DHFR structures, the FH4 complex structure responsible for rate-limiting product release is unknown. We report to our knowledge the first crystal structure of an E. coli. DHFR:FH4 complex at 1.03 Å resolution showing distinct stabilizing interactions absent in FH2 or related (6R)-5,10-dideaza-FH4 complexes. We discover the time course of decay of the co-purified endogenous FH4 during crystal growth, with conversion from FH4 to FH2 occurring in 2–3 days. We also determine another occluded complex structure of E. coli DHFR with a slow-onset nanomolar inhibitor that contrasts with the methotrexate complex, suggesting a plausible strategy for designing DHFR antibiotics by targeting FH4 product conformations.

Hongnan Cao et al. present the X-ray crystal structure of E. coli dihydrofolate reductase (DHFR) in complex with its reduced substrate, (6s)-5,6,7,8-tetrahydrofolate (FH4). This structure provides the first glimpse of the rate-limiting product release step of the DHFR mechanism and suggests a strategy for designing DHFR-targeting antibiotics.

Increased substrate affinity in the Escherichia coli L28R dihydrofolate reductase mutant causes trimethoprim resistance

Dihydrofolate reductase (DHFR) is a ubiquitous enzyme with an essential role in cell metabolism. DHFR catalyzes the reduction of dihydrofolate to tetrahydrofolate, which is a precursor for purine and thymidylate synthesis. Several DHFR targeting antifolate drugs including trimethoprim, a competitive antibacterial inhibitor, have therefore been developed and are clinically used. Evolution of resistance against antifolates is a common public health problem rendering these drugs ineffective. To combat the resistance problem, it is important to understand resistance-conferring changes in the DHFR structure and accordingly develop alternative strategies. Here, we structurally and dynamically characterize Escherichia coli DHFR in its wild type (WT) and trimethoprim resistant L28R mutant forms in the presence of the substrate and its inhibitor trimethoprim. We use molecular dynamics simulations to determine the conformational space, loop dynamics and hydrogen bond distributions at the active site of DHFR for the WT and the L28R mutant. We also report their experimental kcat, Km, and Ki values, accompanied by isothermal titration calorimetry measurements of DHFR that distinguish enthalpic and entropic contributions to trimethoprim binding. Although mutations that confer resistance to competitive inhibitors typically make enzymes more promiscuous and decrease affinity to both the substrate and the inhibitor, strikingly, we find that the L28R mutant has a unique resistance mechanism. While the binding affinity differences between the WT and the mutant for the inhibitor and the substrate are small, the newly formed extra hydrogen bonds with the aminobenzoyl glutamate tail of DHF in the L28R mutant leads to increased barriers for the dissociation of the substrate and the product. Therefore, the L28R mutant indirectly gains resistance by enjoying prolonged binding times in the enzyme-substrate complex. While this also leads to slower product release and decreases the catalytic rate of the L28R mutant, the overall effect is the maintenance of a sufficient product formation rate. Finally, the experimental and computational analyses together reveal the changes that occur in the energetic landscape of DHFR upon the resistance-conferring L28R mutation. We show that the negative entropy associated with the binding of trimethoprim in WT DHFR is due to water organization at the binding interface. Our study lays the framework to study structural changes in other trimethoprim resistant DHFR mutants.

Distinct Roles for Conformational Dynamics in Protein-Ligand Interactions

Conformational dynamics has an established role in enzyme catalysis, but its contribution to ligand binding and specificity is largely unexplored. Here we used the Tiam1 PDZ domain and an engineered variant (QM PDZ) with broadened specificity to investigate the role of structure and conformational dynamics in molecular recognition. Crystal structures of the QM PDZ domain both free and bound to ligands showed structural features central to binding (enthalpy), while nuclear-magnetic-resonance-based methyl relaxation experiments and isothermal titration calorimetry revealed that conformational entropy contributes to affinity. In addition to motions relevant to thermodynamics, slower microsecond to millisecond switching was prevalent in the QM PDZ ligand-binding site consistent with a role in ligand specificity. Our data indicate that conformational dynamics plays distinct and fundamental roles in tuning the affinity (conformational entropy) and specificity (excited-state conformations) of molecular interactions. More broadly, our results have important implications for the evolution, regulation, and design of protein-ligand interactions.

Widespread Perturbation of Function, Structure, and Dynamics by a Conservative Single-Atom Substitution in Thymidylate Synthase

Thymidylate synthase (TSase) is responsible for synthesizing the sole de novo source of dTMP in all organisms. TSase is a drug target, and as such, it has been well studied in terms of both structure and reaction mechanism. Cysteine 146 in Escherichia coli TSase is universally conserved because it serves as the nucleophile in the enzyme mechanism. Here we use the C146S mutation to probe the role of the sulfur atom in early events in the catalytic cycle beyond serving as the nucleophile. Surprisingly, the single-atom substitution severely decreases substrate binding affinity, and the unfavorable ΔΔG°_bind is comprised of roughly equal enthalpic and entropic components at 25 °C. Chemical shifts in the free and dUMP-bound states show the mutation causes perturbations throughout TSase, including regions important for complex stability, in agreement with a less favorable enthalpy change. We measured the nuclear magnetic resonance methyl symmetry axis order parameter (S²_axis), a proxy for conformational entropy, for TSase at all vertices of the dUMP binding/C146S mutation thermodynamic cycle and found that the calculated TΔΔS°_conf is similar in sign and magnitude to the calorimetric TΔΔS°. Further, we ascribed minor resonances in wild-type-dUMP spectra to a state with a covalent bond between Sγ of C146 and C6 of dUMP and find S²_axis values are unaffected by covalent bond formation, indicating this reaction step is neutral with respect to ΔS°_conf. Lastly, the C146S mutation allowed us to measure cofactor analog binding by isothermal titration calorimetry without the confounding heat signature of covalent bond formation. Raltitrexed binds free and singly bound TSase with similar affinities, yet the two binding events have different enthalpy changes, providing further evidence of communication between the two active sites.

Protein motions and dynamic effects in enzyme catalysis

The role of protein motions in promoting the chemical step of enzyme catalysed reactions remains a subject of considerable debate. Here, a unified view of the role of protein dynamics in dihydrofolate reductase catalysis is described. Recently the role of such motions has been investigated by characterising the biophysical properties of isotopically substituted enzymes through a combination of experimental and computational analyses. Together with previous work, these results suggest that dynamic coupling to the chemical coordinate is detrimental to catalysis and may have been selected against during DHFR evolution. The full catalytic power of Nature's catalysts appears to depend on finely tuning protein motions in each step of the catalytic cycle.

Bacterial Thymidylate Synthase Binds Two Molecules of Substrate and Cofactor without Cooperativity

Thymidylate synthase (TSase) is a clinically important enzyme because it catalyzes synthesis of the sole de novo source of deoxy-thymidylate. Without this enzyme, cells die a "thymineless death" since they are starved of a crucial DNA synthesis precursor. As a drug target, TSase is well studied in terms of its structure and reaction mechanism. An interesting mechanistic feature of dimeric TSase is that it is "half-the-sites reactive", which is a form of negative cooperativity. Yet, the basis for this is not well-understood. Some experiments point to cooperativity at the binding steps of the reaction cycle as being responsible for the phenomenon, but the literature contains conflicting reports. Here we use ITC and NMR to resolve these inconsistencies. This first detailed thermodynamic dissection of multisite binding of dUMP to E. coli TSase shows the nucleotide binds to the free and singly bound forms of the enzyme with nearly equal affinity over a broad range of temperatures and in multiple buffers. While small but significant differences in ΔC°P for the two binding events show that the active sites are not formally equivalent, there is little-to-no allostery at the level of ΔG°bind. In addition NMR titration data reveal that there is minor intersubunit cooperativity in formation of a ternary complex with the mechanism based inhibitor, 5F-dUMP, and cofactor. Taken together, the data show that functional communication between subunits is minimal for both binding steps of the reaction coordinate.

Methotrexate-Loaded Four-Arm Star Amphiphilic Block Copolymer Elicits CD8+ T Cell Response against a Highly Aggressive and Metastatic Experimental Lymphoma

We have synthesized a well-defined four-arm star amphiphilic block copolymer [poly(DLLA)-b-poly(NVP)]4 [star-(PDLLA-b-PNVP)4] that consists of D,L-lactide (DLLA) and N-vinylpyrrolidone (NVP) via the combination of ring-opening polymerization (ROP) and xanthate-mediated reversible addition-fragmentation chain transfer (RAFT) polymerization. Synthesis of the polymer was verified by 1H NMR spectroscopy and gel permeation chromatography (GPC). The amphiphilic four-arm star block copolymer forms spherical micelles in water as demonstrated by transmission electron microscopy (TEM) and 1H NMR spectroscopy. Pyrene acts as a probe to ascertain the critical micellar concentration (cmc) by using fluorescence spectroscopy. Methotrexate (MTX)-loaded polymeric micelles of star-(PDLLA15-b-PNVP10)4 amphiphilic block copolymer were prepared and characterized by fluorescence and TEM studies. Star-(PDLLA15-b-PNVP10)4 copolymer was found to be significantly effective with respect to inhibition of proliferation and lysis of human and murine lymphoma cells. The amphiphilic block copolymer causes cell death in parental and MTX-resistant Dalton lymphoma (DL) and Raji cells. The formulation does not cause hemolysis in red blood cells and is tolerant to lymphocytes compared to free MTX. Therapy with MTX-loaded star-(PDLLA15-b-PNVP10)4 amphiphilic block copolymer micelles prolongs the life span of animals with neoplasia by reducing the tumor load, preventing metastasis and augmenting CD8+ T cell-mediated adaptive immune responses.

Dysfunctional conformational dynamics of protein kinase A induced by a lethal mutant of phospholamban hinder phosphorylation

The dynamic interplay between kinases and substrates is crucial for the formation of catalytically committed complexes that enable phosphoryl transfer. However, a clear understanding on how substrates modulate kinase structural dynamics to control catalytic efficiency is still missing. Here, we used solution NMR spectroscopy to study the conformational dynamics of two complexes of the catalytic subunit of the cAMP-dependent protein kinase A with WT and R14 deletion phospholamban, a lethal human mutant linked to familial dilated cardiomyopathy. Phospholamban is a central regulator of heart muscle contractility, and its phosphorylation by protein kinase A constitutes a primary response to β-adrenergic stimulation. We found that the single deletion of arginine in phospholamban's recognition sequence for the kinase reduces its binding affinity and dramatically reduces phosphorylation kinetics. Structurally, the mutant prevents the enzyme from adopting conformations and motions committed for catalysis, with concomitant reduction in catalytic efficiency. Overall, these results underscore the importance of a well-tuned structural and dynamic interplay between the kinase and its substrates to achieve physiological phosphorylation levels for proper Ca(2+) signaling and normal cardiac function.

Synchronous opening and closing motions are essential for cAMP-dependent protein kinase A signaling

Conformational fluctuations play a central role in enzymatic catalysis. However, it is not clear how the rates and the coordination of the motions affect the different catalytic steps. Here, we used NMR spectroscopy to analyze the conformational fluctuations of the catalytic subunit of the cAMP-dependent protein kinase (PKA-C), a ubiquitous enzyme involved in a myriad of cell signaling events. We found that the wild-type enzyme undergoes synchronous motions involving several structural elements located in the small lobe of the kinase, which is responsible for nucleotide binding and release. In contrast, a mutation (Y204A) located far from the active site desynchronizes the opening and closing of the active cleft without changing the enzyme's structure, rendering it catalytically inefficient. Since the opening and closing motions govern the rate-determining product release, we conclude that optimal and coherent conformational fluctuations are necessary for efficient turnover of protein kinases.

Current and emerging opportunities for molecular simulations in structure-based drug design

An overview of the current capabilities and limitations of molecular simulation of biomolecular complexes in the context of computer-aided drug design is provided. Steady improvements in computer hardware coupled with more refined representations of energetics are leading to a new appreciation of the driving forces of molecular recognition. Molecular simulations are poised to more frequently guide the interpretation of biophysical measurements of biomolecular complexes. Ligand design strategies emerge from detailed analyses of computed structural ensembles. The feasibility of routine applications to ligand optimization problems hinges upon successful extensive large scale validation studies and the development of protocols to intelligently automate computations.

Conformation-based restrictions and scaffold replacements in the design of hepatitis C virus polymerase inhibitors: discovery of deleobuvir (BI 207127)

Conformational restrictions of flexible torsion angles were used to guide the identification of new chemotypes of HCV NS5B inhibitors. Sites for rigidification were based on an acquired conformational understanding of compound binding requirements and the roles of substituents in the free and bound states. Chemical bioisosteres of amide bonds were explored to improve cell-based potency. Examples are shown, including the design concept that led to the discovery of the phase III clinical candidate deleobuvir (BI 207127). The structure-based strategies employed have general utility in drug design.

Ligand bioactive conformation plays a critical role in the design of drugs that target the hepatitis C virus NS3 protease

A ligand-focused strategy employed NMR, X-ray, modeling, and medicinal chemistry to expose the critical role that bioactive conformation played in the design of a variety of drugs that target the HCV protease. The bioactive conformation (bound states) were determined for key inhibitors identified along our drug discovery pathway from the hit to clinical compounds. All adopt similar bioactive conformations for the common core derived from the hit peptide DDIVPC. A carefully designed SAR analysis, based on the advanced inhibitor 1 in which the P1 to P3 side chains and the N-terminal Boc were sequentially truncated, revealed a correlation between affinity and the relative predominance of the bioactive conformation in the free state. Interestingly, synergistic conformation effects on potency were also noted. Comparisons with clinical and recently marketed drugs from the pharmaceutical industry showed that all have the same core and similar bioactive conformations. This suggested that the variety of appendages discovered for these compounds also properly satisfy the bioactive conformation requirements and allowed for a large variety of HCV protease drug candidates to be designed.

Allosteric activation of the Par-6 PDZ via a partial unfolding transition

Proteins exist in a delicate balance between the native and unfolded states, where thermodynamic stability may be sacrificed to attain the flexibility required for efficient catalysis, binding, or allosteric control. Partition-defective 6 (Par-6) regulates the Par polarity complex by transmitting a GTPase signal through the Cdc42/Rac interaction binding PSD-95/Dlg/ZO-1 (CRIB-PDZ) module that alters PDZ ligand binding. Allosteric activation of the PDZ is achieved by local rearrangement of the L164 and K165 side chains to stabilize the interdomain CRIB:PDZ interface and reposition a conserved element of the ligand binding pocket. However, microsecond to millisecond dynamics measurements revealed that L164/K165 exchange requires a larger rearrangement than expected. The margin of thermodynamic stability for the PDZ domain is modest (∼3 kcal/mol) and further reduced by transient interactions with the disordered CRIB domain. Measurements of local structural stability revealed that tertiary contacts within the PDZ are disrupted by a partial unfolding transition that enables interconversion of the L/K switch. The unexpected participation of partial PDZ unfolding in the allosteric mechanism of Par-6 suggests that native-state unfolding may be essential for the function of other marginally stable proteins.

Functional significance of evolving protein sequence in dihydrofolate reductase from bacteria to humans

With the rapidly growing wealth of genomic data, experimental inquiries on the functional significance of important divergence sites in protein evolution are becoming more accessible. Here we trace the evolution of dihydrofolate reductase (DHFR) and identify multiple key divergence sites among 233 species between humans and bacteria. We connect these sites, experimentally and computationally, to changes in the enzyme's binding properties and catalytic efficiency. One of the identified evolutionarily important sites is the N23PP modification (∼mid-Devonian, 415-385 Mya), which alters the conformational states of the active site loop in Escherichia coli dihydrofolate reductase and negatively impacts catalysis. This enzyme activity was restored with the inclusion of an evolutionarily significant lid domain (G51PEKN in E. coli enzyme; ∼2.4 Gya). Guided by this evolutionary genomic analysis, we generated a human-like E. coli dihydrofolate reductase variant through three simple mutations despite only 26% sequence identity between native human and E. coli DHFRs. Molecular dynamics simulations indicate that the overall conformational motions of the protein within a common scaffold are retained throughout evolution, although subtle changes to the equilibrium conformational sampling altered the free energy barrier of the enzymatic reaction in some cases. The data presented here provide a glimpse into the evolutionary trajectory of functional DHFR through its protein sequence space that lead to the diverged binding and catalytic properties of the E. coli and human enzymes.

A surprising role for conformational entropy in protein function

Formation of high-affinity complexes is critical for the majority of enzymatic reactions involving proteins. The creation of the family of Michaelis and other intermediate complexes during catalysis clearly involves a complicated manifold of interactions that are diverse and complex. Indeed, computing the energetics of interactions between proteins and small molecule ligands using molecular structure alone remains a great challenge. One of the most difficult contributions to the free energy of protein-ligand complexes to access experimentally is that due to changes in protein conformational entropy. Fortunately, recent advances in solution nuclear magnetic resonance (NMR) relaxation methods have enabled the use of measures-of-motion between conformational states of a protein as a proxy for conformational entropy. This review briefly summarizes the experimental approaches currently employed to characterize fast internal motion in proteins, how this information is used to gain insight into conformational entropy, what has been learned, and what the future may hold for this emerging view of protein function.

The dark energy of proteins comes to light: conformational entropy and its role in protein function revealed by NMR relaxation

Historically it has been virtually impossible to experimentally determine the contribution of residual protein entropy to fundamental protein activities such as the binding of ligands. Recent progress has illuminated the possibility of employing NMR relaxation methods to quantitatively determine the role of changes in conformational entropy in molecular recognition by proteins. The method rests on using fast internal protein dynamics as a proxy. Initial results reveal a large and variable role for conformational entropy in the binding of ligands by proteins. Such a role for conformational entropy in molecular recognition has significant implications for enzymology, signal transduction, allosteric regulation and the development of protein-directed pharmaceuticals.

Segmental motions, not a two-state concerted switch, underlie allostery in CheY

The switch between an inactive and active conformation is an important transition for signaling proteins, yet the mechanisms underlying such switches are not clearly understood. Escherichia coli CheY, a response regulator protein from the two-component signal transduction system that regulates bacterial chemotaxis, is an ideal protein for the study of allosteric mechanisms. By using 15N CPMG relaxation dispersion experiments, we monitored the inherent dynamic switching of unphosphorylated CheY. We show that CheY does not undergo a two-state concerted switch between the inactive and active conformations. Interestingly, partial saturation of Mg2+ enhances the intrinsic allosteric motions. Taken together with chemical shift perturbations, these data indicate that the μs-ms timescale motions underlying CheY allostery are segmental in nature. We propose an expanded allosteric network of residues, including W58, that undergo asynchronous, local switching between inactive and active-like conformations as the primary basis for the allosteric mechanism.

Ligand and receptor dynamics contribute to the mechanism of graded PPARγ agonism

Ligand binding to proteins is not a static process, but rather involves a number of complex dynamic transitions. A flexible ligand can change conformation upon binding its target. The conformation and dynamics of a protein can change to facilitate ligand binding. The conformation of the ligand, however, is generally presumed to have one primary binding mode, shifting the protein conformational ensemble from one state to another. We report solution nuclear magnetic resonance (NMR) studies that reveal peroxisome proliferator-activated receptor γ (PPARγ) modulators can sample multiple binding modes manifesting in multiple receptor conformations in slow conformational exchange. Our NMR, hydrogen/deuterium exchange and docking studies reveal that ligand-induced receptor stabilization and binding mode occupancy correlate with the graded agonist response of the ligand. Our results suggest that ligand and receptor dynamics affect the graded transcriptional output of PPARγ modulators.

Direct and allosteric inhibition of the FGF2/HSPGs/FGFR1 ternary complex formation by an antiangiogenic, thrombospondin-1-mimic small molecule

Fibroblast growth factors (FGFs) are recognized targets for the development of therapies against angiogenesis-driven diseases, including cancer. The formation of a ternary complex with the transmembrane tyrosine kinase receptors (FGFRs), and heparan sulphate proteoglycans (HSPGs) is required for FGF2 pro-angiogenic activity. Here by using a combination of techniques including Nuclear Magnetic Resonance, Molecular Dynamics, Surface Plasmon Resonance and cell-based binding assays we clarify the molecular mechanism of inhibition of an angiostatic small molecule, sm27, mimicking the endogenous inhibitor of angiogenesis, thrombospondin-1. NMR and MD data demonstrate that sm27 engages the heparin-binding site of FGF2 and induces long-range dynamics perturbations along FGF2/FGFR1 interface regions. The functional consequence of the inhibitor binding is an impaired FGF2 interaction with both its receptors, as demonstrated by SPR and cell-based binding assays. We propose that sm27 antiangiogenic activity is based on a twofold-direct and allosteric-mechanism, inhibiting FGF2 binding to both its receptors.

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 ¹⁵N and ²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.

A photochemical approach for controlled drug release in targeted drug delivery

Photochemistry provides a unique mechanism that enables the active control of drug release in cancer-targeting drug delivery. This study investigates the light-mediated release of methotrexate, an anticancer drug, using a photocleavable linker strategy based on o-nitrobenzyl protection. We evaluated two types of the o-nitrobenzyl-linked methotrexate for the drug release study and further extended the study to a fifth-generation poly(amidoamine) dendrimer carrier covalently conjugated with methotrexate via the o-nitrobenzyl linker. We performed the drug release studies by using a combination of three standard analytical methods that include UV/vis spectrometry, (1)H NMR spectroscopy, and anal. HPLC. This article reports that methotrexate is released by the photochemical mechanism in an actively controlled manner. The rate of the drug release varies in response to multiple control parameters, including linker design, light wavelength, exposure time, and the pH of the medium where the drug release occurs.

Evidence for dynamics in proteins as a mechanism for ligand dissociation

Signal transduction, regulatory processes, and pharmaceutical responses are highly dependent upon ligand residence times. Gaining insight into how physical factors influence residence times, or k_off, should enhance our ability to manipulate biological interactions. We report experiments that yield structural insight into k_off for a series of eight 2,4-diaminopyrimidine inhibitors of dihydrofolate reductase that vary by six orders of magnitude in binding affinity. NMR relaxation dispersion experiments revealed a common set of residues near the binding site that undergo a concerted, millisecond-timescale switching event to a previously unidentified conformation. The rate of switching from ground to excited conformations correlates exponentially with K_i and k_off, suggesting that protein dynamics serves as a mechanical initiator of ligand dissociation within this series and potentially for other macromolecule-ligand systems. Although k_conf,forward is faster than k_off, use of the ligand series allowed for connections to be drawn between kinetic events on different timescales.

Photochemical release of methotrexate from folate receptor-targeting PAMAM dendrimer nanoconjugate

Nanoparticle (NP)-based targeted drug delivery involves cell-specific targeting followed by a subsequent therapeutic action from the therapeutic carried by the NP system. NPs conjugated with methotrexate (MTX), a potent inhibitor of dihydrofolate reductase (DHFR) localized in cytosol, have been under investigation as a delivery system to target cancer cells to enhance the therapeutic index of methotrexate, which is otherwise non-selectively cytotoxic. Despite improved therapeutic activity from MTX-conjugated NPs in vitro and in vivo, the therapeutic action of these conjugates following cellular entry is poorly understood; in particular it is unclear whether the therapeutic activity requires release of the MTX. This study investigates whether MTX must be released from a nanoparticle in order to achieve the therapeutic activity. We report herein light-controlled release of methotrexate from a dendrimer-based conjugate and provide evidence suggesting that MTX still attached to the nanoconjugate system is fully able to inhibit the activity of its enzyme target and the growth of cancer cells.

Identification of endogenous ligands bound to bacterially expressed human and E. coli dihydrofolate reductase by 2D NMR

Dihydrofolate reductase (DHFR) is a well-studied drug target and a paradigm for understanding enzyme catalysis. Preparation of pure DHFR samples, in defined ligand-bound states, is a prerequisite for in vitro studies and drug discovery efforts. We use NMR spectroscopy to monitor ligand content of human and Escherichia coli DHFR (ecDHFR), which bind different co-purifying ligands during expression in bacteria. An alternate purification strategy yields highly pure DHFR complexes, containing only the desired ligands, in the quantities required for structural studies. Interestingly, ecDHFR is bound to endogenous THF while human DHFR is bound to NADP. Consistent with these findings, a designed "humanized" mutant of ecDHFR switches binding specificity in the cell.

Optimization of NMR spectroscopy of encapsulated proteins dissolved in low viscosity fluids

Comprehensive application of solution NMR spectroscopy to studies of macromolecules remains fundamentally limited by the molecular rotational correlation time. For proteins, molecules larger than 30 kDa require complex experimental methods, such as TROSY in conjunction with isotopic labeling schemes that are often expensive and generally reduce the potential information available. We have developed the reverse micelle encapsulation strategy as an alternative approach. Encapsulation of proteins within the protective nano-scale water pool of a reverse micelle dissolved in ultra-low viscosity nonpolar solvents overcomes the slow tumbling problem presented by large proteins. Here, we characterize the contributions from the various components of the protein-containing reverse micelle system to the rotational correlation time of the encapsulated protein. Importantly, we demonstrate that the protein encapsulated in the reverse micelle maintains a hydration shell comparable in size to that seen in bulk solution. Using moderate pressures, encapsulation in ultra-low viscosity propane or ethane can be used to magnify this advantage. We show that encapsulation in liquid ethane can be used to reduce the tumbling time of the 43 kDa maltose binding protein from ~23 to ~10 ns. These conditions enable, for example, acquisition of TOCSY-type data resolved on the adjacent amide NH for the 43 kDa encapsulated maltose binding protein dissolved in liquid ethane, which is typically impossible for proteins of such size without use of extensive deuteration or the TROSY effect.

NMR reveals novel mechanisms of protein activity regulation

NMR spectroscopy is one of the most powerful tools for the characterization of biomolecular systems. A unique aspect of NMR is its capacity to provide an integrated insight into both the structure and intrinsic dynamics of biomolecules. In addition, NMR can provide site-resolved information about the conformation entropy of binding, as well as about energetically excited conformational states. Recent advances have enabled the application of NMR for the characterization of supramolecular systems. A summary of mechanisms underpinning protein activity regulation revealed by the application of NMR spectroscopy in a number of biological systems studied in the lab is provided.

Dynamically committed, uncommitted, and quenched states encoded in protein kinase A revealed by NMR spectroscopy

Protein kinase A (PKA) is a ubiquitous phosphoryl transferase that mediates hundreds of cell signaling events. During turnover, its catalytic subunit (PKA-C) interconverts between three major conformational states (open, intermediate, and closed) that are dynamically and allosterically activated by nucleotide binding. We show that the structural transitions between these conformational states are minimal and allosteric dynamics encode the motions from one state to the next. NMR and molecular dynamics simulations define the energy landscape of PKA-C, with the substrate allowing the enzyme to adopt a broad distribution of conformations (dynamically committed state) and the inhibitors (high magnesium and pseudosubstrate) locking it into discrete minima (dynamically quenched state), thereby reducing the motions that allow turnover. These results unveil the role of internal dynamics in both kinase function and regulation.

Direct detection of structurally resolved dynamics in a multiconformation receptor-ligand complex

Structure-based drug design relies on static protein structures despite significant evidence for the need to include protein dynamics as a serious consideration. In practice, dynamic motions are neglected because they are not understood well enough to model, a situation resulting from a lack of explicit experimental examples of dynamic receptor-ligand complexes. Here, we report high-resolution details of pronounced ~1 ms time scale motions of a receptor-small molecule complex using a combination of NMR and X-ray crystallography. Large conformational dynamics in Escherichia coli dihydrofolate reductase are driven by internal switching motions of the drug-like, nanomolar-affinity inhibitor. Carr-Purcell-Meiboom-Gill relaxation dispersion experiments and NOEs revealed the crystal structure to contain critical elements of the high energy protein-ligand conformation. The availability of accurate, structurally resolved dynamics in a protein-ligand complex should serve as a valuable benchmark for modeling dynamics in other receptor-ligand complexes and prediction of binding affinities.

Importance of ligand bioactive conformation in the discovery of potent indole-diamide inhibitors of the hepatitis C virus NS5B

Significant advances have led to receptor induced-fit and conformational selection models for describing bimolecular recognition, but a more comprehensive view must evolve to also include ligand shape and conformational changes. Here, we describe an example where a ligand's "structural hinge" influences potency by inducing an "L-shape" bioactive conformation, and due to its solvent exposure in the complex, reasonable conformation-activity-relationships can be qualitatively attributed. From a ligand design perspective, this feature was exploited by successful linker hopping to an alternate "structural hinge" that led to a new and promising chemical series which matched the ligand bioactive conformation and the pocket bioactive space. Using a combination of X-ray crystallography, NMR and modeling with support from binding-site resistance mutant studies and photoaffinity labeling experiments, we were able to derive inhibitor-polymerase complexes for various chemical series.

Macromolecular NMR spectroscopy for the non-spectroscopist: beyond macromolecular solution structure determination

A strength of NMR spectroscopy is its ability to monitor, on an atomic level, molecular changes and interactions. In this review, which is intended for non-spectroscopist, we describe major uses of NMR in protein science beyond solution structure determination. After first touching on how NMR can be used to quickly determine whether a mutation induces structural perturbations in a protein, we describe the unparalleled ability of NMR to monitor binding interactions over a wide range of affinities, molecular masses and solution conditions. We discuss the use of NMR to measure the dynamics of proteins at the atomic level and over a wide range of timescales. Finally, we outline new and expanding areas such as macromolecular structure determination in multicomponent systems, as well as in the solid state and in vivo.

Mining electron density for functionally relevant protein polysterism in crystal structures

This review focuses on conceptual and methodological advances in our understanding and characterization of the conformational heterogeneity of proteins. Focusing on X-ray crystallography, we describe how polysterism, the interconversion of pre-existing conformational substates, has traditionally been analyzed by comparing independent crystal structures or multiple chains within a single crystal asymmetric unit. In contrast, recent studies have focused on mining electron density maps to reveal previously ‘hidden’ minor conformational substates. Functional tests of the importance of minor states suggest that evolutionary selection shapes the entire conformational landscape, including uniquely configured conformational substates, the relative distribution of these substates, and the speed at which the protein can interconvert between them. An increased focus on polysterism may shape the way protein structure and function is studied in the coming years.