Search

The ABC efflux pump P-glycoprotein (P-gp) transports a wide variety of drugs and is inhibited by others. Some drugs stimulate ATP hydrolysis at the nucleotide binding domains (NBDs) and are transported, others uncouple ATP hydrolysis and transport, and others inhibit ATP hydrolysis. The molecular basis for the different behavior of these drugs is not well understood despite the availability of several structural models of P-gp complexes with ligands bound. Hypothetically, ligands differentially alter the conformational dynamics of peptide segments that mediate the coupling between the drug binding sites and the NBDs. Here, we explore by hydrogen-deuterium exchange mass spectrometry the dynamic consequences of a classic substrate and inhibitor, vinblastine and zosuquidar, binding to mouse P-gp (mdr1a) in lipid nanodiscs. The dynamics of P-gp in nucleotide-free, pre-hydrolysis, and post-hydrolysis states in the presence of each drug reveal distinct mechanisms of ATPase stimulation and implications for transport. For both drugs, there are common regions affected in a similar manner, suggesting that particular networks are the key to stimulating ATP hydrolysis. However, drug binding effects diverge in the post-hydrolysis state, particularly in the intracellular helices (ICHs 3 and 4) and neighboring transmembrane helices. The local dynamics and conformational equilibria in this region are critical for the coupling of drug binding and ATP hydrolysis and are differentially modulated in the catalytic cycle.

Antibody-based therapeutics are the fastest growing drug class on the market, used to treat aggressive forms of cancer, chronic autoimmune conditions, and numerous other disease states. While the specificity, affinity, and versatility of therapeutic antibodies can provide an advantage over traditional small molecule drugs, their development and optimization can be much more challenging and time-consuming. This is, in part, because the ideal formulation buffer systems used for in vitro characterization inadequately reflect the crowded biological environments (serum, endosomal lumen, etc.) that these drugs experience once administered to a patient. Such environments can perturb the binding of antibodies to their antigens and receptors, as well as homo- and hetero-aggregation, in ways that are incompletely understood, thereby altering therapeutic effect and disposition. While excluded volume effects are classically thought to favor binding, weak interactions with co-solutes in crowded conditions can inhibit binding. The second virial coefficient (B2) parameter quantifies such weak interactions and can be determined by a variety of techniques in dilute solution, but analogous methods in complex biological fluids are not well established. Here, we demonstrate that fluorescence correlation spectroscopy (FCS) is able to measure diffusive B2 values directly in undiluted serum. Apparent second virial coefficient (B2,app) measurements of antibodies in serum reveal that changes in the balance between attractive and repulsive interactions can dramatically impact global nonideality. Furthermore, our findings suggest that the common approach of isolating specific components and completing independent cross-term virial coefficient measurements is an incomplete representation of nonideality in serum. The approach presented here could enrich our understanding of the effects of biological environments on proteins in general, and advance the development of therapeutic antibodies and other protein-based therapeutics.STATEMENT OF SIGNIFICANCEWe present FCS as an orthogonal method to traditional methods for characterizing weak, nonspecific homo- and hetero-interactions through determination of self- and cross-term second virial coefficients, respectively. We also characterize weak interactions between therapeutic antibodies and serum components through determination of an apparent second virial coefficient (B2,app) directly in undiluted serum. Our results suggest that global nonideality effects are antibody-dependent, and that attractive and repulsive interactions with co-solutes are occurring simultaneously. This approach could advance our understanding of the impact of nonideality to the biophysical and pharmacological properties of therapeutic antibodies and other engineered proteins in relevant biological environments, and could accelerate the development and optimization of protein-based therapeutics.

Monoclonal antibodies (mAbs) have become an important class of therapeutics, particularly in the realm of anticancer immunotherapy. While the two antigen-binding fragments (Fabs) of an mAb allow for high-avidity binding to molecular targets, the crystallizable fragment (Fc) engages immune effector elements. mAbs of the IgG class are used for the treatment of autoimmune diseases and can elicit antitumor immune functions not only by several mechanisms including direct antigen engagement via their Fab arms but also by Fab binding to tumors combined with Fc engagement of complement component C1q and Fcγ receptors. Additionally, IgG binding to the neonatal Fc receptor (FcRn) allows for endosomal recycling and prolonged serum half-life. To augment the effector functions or half-life of an IgG1 mAb, we constructed a novel "2Fc" mAb containing two Fc domains in addition to the normal two Fab domains. Structural and functional characterization of this 2Fc mAb demonstrated that it exists in a tetrahedral-like geometry and retains binding capacity via the Fab domains. Furthermore, duplication of the Fc region significantly enhanced avidity for Fc receptors FcγRI, FcγRIIIa, and FcRn, which manifested as a decrease in complex dissociation rate that was more pronounced at higher densities of receptor. At intermediate receptor density, the dissociation rate for Fc receptors was decreased 6- to 130-fold, resulting in apparent affinity increases of 7- to 42-fold. Stoichiometric analysis confirmed that each 2Fc mAb may simultaneously bind two molecules of FcγRI or four molecules of FcRn, which is double the stoichiometry of a wild-type mAb. In summary, duplication of the IgG Fc region allows for increased avidity to Fc receptors that could translate into clinically relevant enhancement of effector functions or pharmacokinetics.

P-Glycoprotein (P-gp) is an ATP-dependent efflux pump that clears a wide variety of drugs and toxins from cells. P-gp undergoes large-scale structural changes and demonstrates conformational heterogeneity even within a single catalytic or drug-bound state, although the role of heterogeneity remains unclear. P-gp is found in a variety of cell types that vary in lipid composition, which modulates its activity. An understanding of structural or dynamic changes due to the lipid environment is lacking. We aimed to determine the effects of cholesterol in a membrane on the conformational behavior of P-gp in lipid nanodiscs. The presence of cholesterol stimulates ATP hydrolysis and alters lipid order and fluidity. Hydrogen/deuterium exchange mass spectrometry demonstrates that cholesterol in the membrane induces asymmetric, long-range changes in the distributions and exchange kinetics of conformations of the nucleotide-binding domains, correlating the effects of lipid composition on activity with specific changes in the P-gp conformational landscape.

Aromatase (CYP19A1) catalyzes the synthesis of estrogens from androgens and is an invaluable target of pharmacotherapy for estrogen-dependent cancers. CYP19A1 is also one of the most primordial human CYPs and, to the extent that its fundamental dynamics are conserved, is highly relevant to understanding those of the more recently evolved and promiscuous enzymes. A complementary approach employing molecular dynamics simulations and hydrogen-deuterium exchange mass spectrometry (HDX-MS) was employed to interrogate the changes in CYP19A1 dynamics coupled to binding androstenedione (ASD). Gaussian-accelerated molecular dynamics and HDX-MS agree that ASD globally suppresses CYP19A1 dynamics. Bimodal HDX patterns of the B'-C loop potentially arising from at least two conformations are present in free 19A1 only, supporting the possibility that conformational selection is operative. Random-acceleration molecular dynamics and adaptive biasing force simulations illuminate ASD's binding pathway, predicting ASD capture in the lipid headgroups and a pathway to the active site shielded from solvent. Intriguingly, the predicted access channel in 19A1 aligns well with the steroid binding sites of other human sterol-oxidizing CYPs.

Promiscuous and allosteric drug interactions with cytochrome P450 3A4 (CYP3A4) are ubiquitous but incompletely understood at the molecular level. A classic allosteric CYP3A4 drug interaction includes the benzodiazepine midazolam (MDZ). MDZ exhibits homotropic and heterotropic allostery when metabolized to 1'-hydroxy and 4-hydroxy metabolites in varying ratios. The combination of hydrogen-deuterium exchange mass spectrometry (HDX-MS) and Gaussian accelerated molecular dynamics (GaMD) simulations of CYP3A4 in lipid nanodiscs and in a lipid bilayer, respectively, reveals MDZ-dependent changes in dynamics in a membrane environment. The F-, G-, and intervening helices, as well as the loop preceding the β1-sheets, display the largest observed changes in HDX. The GaMD suggests a potential allosteric binding site for MDZ in the F'- and G'-regions, which undergo significant increases in HDX at near-saturating MDZ concentrations. The HDX-MS and GaMD results confirm that changes in dynamics are most significant near the developing consensus allosteric site, and these changes are distinct from those observed previously with the nonallosteric inhibitor ketoconazole. The results suggest that the allosteric MDZ remains mobile in its binding site at the Phe-cluster. The results further suggest that this binding site remains dynamic or changes the depth of insertion in the membrane.

Cytochrome P450 (CYP) monoxygenses utilize heme cofactors to catalyze oxidation reactions. They play a critical role in metabolism of many classes of drugs, are an attractive target for drug development, and mediate several prominent drug interactions. Many substrates and inhibitors alter the spin state of the ferric heme by displacing the heme's axial water ligand in the resting enzyme to yield a five-coordinate iron complex, or they replace the axial water to yield a nitrogen-ligated six-coordinate iron complex, which are traditionally assigned by UV-vis spectroscopy. However, crystal structures and recent pulsed electron paramagnetic resonance (EPR) studies find a few cases where molecules hydrogen bond to the axial water. The water-bridged drug-H2O-heme has UV-vis spectra similar to nitrogen-ligated, six-coordinate complexes, but are closer to "reverse type I" complexes described in older liteature. Here, pulsed and continuous wave (CW) EPR demonstrate that water-bridged complexes are remarkably common among a range of nitrogenous drugs or drug fragments that bind to CYP3A4 or CYP2C9. Principal component analysis reveals a distinct clustering of CW EPR spectral parameters for water-bridged complexes. CW EPR reveals heterogeneous mixtures of ligated states, including multiple directly-coordinated complexes and water-bridged complexes. These results suggest that water-bridged complexes are under-represented in CYP structural databases and can have energies similar to other ligation modes. The data indicates that water-bridged binding modes can be identified and distinguished from directly-coordinated binding by CW EPR.

Bispecific antibodies (bsAbs) combine the antigen specificities of two distinct Abs and demonstrate therapeutic promise based on novel mechanisms of action. Among the many platforms for creating bsAbs, controlled Fab-arm exchange (cFAE) has proven useful based on minimal changes to native Ab structure and the simplicity with which bsAbs can be formed from two parental Abs. Despite a published protocol for cFAE and its widespread use in the pharmaceutical industry, the reaction mechanism has not been determined. Knowledge of the mechanism could lead to improved yields of bsAb at faster rates as well as foster adoption of process control. In this work, a combination of Forster resonance energy transfer (FRET), nonreducing SDS-PAGE, and strategic mutation of the Ab hinge region was employed to identify and characterize the individual steps of cFAE. Fluorescence correlation spectroscopy (FCS) was used to determine the affinity of parental (homodimer) and bispecific (heterodimer) interactions within the C(H)3 domain, further clarifying the thermodynamic basis for bsAb formation. The result is a clear sequence of events with rate constants that vary with experimental conditions, where dissociation of the K409R parental Ab into half-Ab controls the rate of the reaction

Ligand-dependent changes in protein conformation are foundational to biology. Historical mechanistic models for substrate-specific proteins are induced fit (IF) and conformational selection (CS), which invoke a change in protein conformation after ligand binds or before ligand binds, respectively. These mechanisms have important, but rarely discussed, functional relevance because IF vs. CS can differentially affect a protein's substrate specificity or promiscuity, and its regulatory properties. The modern view of proteins as conformational ensembles in both ligand free and bound states, together with the realization that most proteins exhibit some substrate promiscuity, demands a deeper interpretation of the historical models and provides an opportunity to improve mechanistic analyses. Here we describe alternative analytical strategies for distinguishing the historical models, including the more complex expanded versions of IF and CS. Functional implications of the different models are described. We provide an alternative perspective based on protein ensembles interacting with ligand ensembles that clarifies how a single protein can ‘apparently’ exploit different mechanisms for different ligands. Mechanistic information about protein ensembles can be optimized when they are probed with multiple ligands.

The ATP binding cassette transporter P-glycoprotein (ABCB1 or P-gp) plays a major role in cellular resistance to drugs and drug interactions. Experimental studies support a mechanism with nucleotide-dependent fluctuation between inward-facing and outward-facing conformations, which are coupled to nucleotide hydrolysis. However, detailed insight into drug-dependent modulation of these conformational ensembles is lacking. Different drugs likely occupy partially overlapping but distinct sites and are therefore variably coupled to nucleotide binding and hydrolysis. Many fluorescent drug analogues are used in cell-based transport models; however, their specific interactions with P-gp have not been studied, and this limits interpretation of transport assays in terms of molecular models. Here we monitor binding of the fluorescent probe substrates BODIPY-verapamil, BODIPY-vinblastine, and Flutax-2 at low occupancy to murine P-gp in lipid nanodiscs via fluorescence correlation spectroscopy, in variable nucleotide-bound states. Changes in affinity for the different nucleotide-dependent conformations are probe-dependent. For BODIPY-verapamil and BODIPY-vinblastine, there are 2–10-fold increases in KD in the nucleotide-bound or vanadate-trapped state, compared to that in the nucleotide-free state. In contrast, the affinity of Flutax-2 is unaffected by nucleotide or vanadate trapping. In further contrast to BODIPY-verapamil and BODIPY-vinblastine, Flutax-2 does not cause stimulation of ATP hydrolysis despite the fact that it is transported in vesicle-based transport assays. Whereas the established substrates verapamil, paclitaxel, and vinblastine displace BODIPY-verapamil or BODIPY-vinblastine from their high-affinity sites, the transport substrate Flutax-2 is not displaced by any of these substrates. The results demonstrate a unique binding site for Flutax-2 that allows for transport without stimulation of ATP hydrolysis.

Submicrometer aggregates are frequently present at low levels in antibody-based therapeutics. Although intuition suggests that the fraction of the aggregate or the size of the aggregate present might correlate with deleterious clinical properties or formulation difficulties, it has been challenging to demonstrate which aggregate states, if any, trigger specific biological effects. One source of uncertainty about the putative linkage between aggregation and safety or efficacy lies in the likelihood that noncovalent aggregation differs in ideal buffers versus in serum and biological tissues; self-association or association with other proteins may vary widely with environment. Therefore, methods for monitoring aggregation and aggregate behavior in biologically relevant matrices could provide a tool for better predicting aggregate-dependent clinical outcomes and provide a basis for antibody engineering prior to clinical studies. Here, we generate models for soluble aggregates of THIOMABs and a bispecific antibody (bsAb) of defined size and exploit fluorescence correlation spectroscopy to monitor their diffusion properties in serum and viscosity-matched buffers. The monomers, dimers, and trimers of both THIOMABs and a bsAb reveal a modest increase in diffusion time in serum greater than expected for an increase in viscosity alone. A mixture of larger aggregates containing mostly bsAb pentamers exhibits a marked increase in diffusion time in serum and much greater intrasample variability, consistent with significant aggregation or interactions with serum components. The results indicate that small aggregates of several IgG platforms are not likely to aggregate with serum components, but nanometer-scale aggregates larger than trimers can interact with the serum in an Ab-dependent manner.

Lipid nanodiscs provide a native-like lipid environment for membrane proteins, and they have become a valuable platform for the study of membrane biophysics. A range of biophysical and biochemical analyses are enabled when membrane proteins are captured in lipid nanodiscs. Two parameters that can be controlled when capturing membrane proteins in lipid nanodiscs are the radius, and hence the surface area of the lipid surface, and the composition of the lipid bilayer. Despite their emergence as a versatile tool, most studies with lipid nanodiscs in the literature have focused on nanodiscs of a single radius with a single lipid. In light of the complexity of biological membranes, it is likely that nanodiscs with multiple membrane components would be more sophisticated models for membrane research. It is possible to prepare nanodiscs with more complex lipid mixtures to probe the effects of lipid composition on several aspects of membrane biochemistry. Detailed protocols are described here for the preparation of nanodiscs with mixtures of phospholipids, incorporation of cholesterol, and incorporation of a spectroscopic lipid probe. These protocols provide starting points for the construction of nanodiscs with more physiological membrane compositions or with useful biophysical probes. © 2019 by John Wiley & Sons, Inc. Basic Protocol 1: Assembly of mixed lipid nanodiscs Basic Protocol 2: Assembly of nanodiscs with cholesterol Basic Protocol 3: Incorporation of laurdan into nanodiscs for membrane fluidity measurements.

Detoxication, or 'drug-metabolizing', enzymes and drug transporters exhibit remarkable substrate promiscuity and catalytic promiscuity. In contrast to substrate-specific enzymes that participate in defined metabolic pathways, individual detoxication enzymes must cope with substrates of vast structural diversity, including previously unencountered environmental toxins. Presumably, evolution selects for a balance of 'adequate' kcat /KM values for a wide range of substrates, rather than optimizing kcat /KM for any individual substrate. However, the structural, energetic, and metabolic properties that achieve this balance, and hence optimize detoxication, are not well understood. Two features of detoxication enzymes that are frequently cited as contributions to promiscuity include the exploitation of highly reactive versatile cofactors, or cosubstrates, and a high degree of flexibility within the protein structure. This review examines these intuitive mechanisms in detail and clarifies the contributions of the classic ligand binding models 'induced fit' (IF) and 'conformational selection' (CS) to substrate promiscuity. The available literature data for drug metabolizing enzymes and transporters suggest that IF is exploited by these promiscuous detoxication enzymes, as it is with substrate-specific enzymes, but the detoxication enzymes uniquely exploit 'IFs' to retain a wide range of substrates at their active sites. In contrast, whereas CS provides no catalytic advantage to substrate-specific enzymes, promiscuous enzymes may uniquely exploit it to recruit a wide range of substrates. The combination of CS and IF, for recruitment and retention of substrates, can potentially optimize the promiscuity of drug metabolizing enzymes and drug transporters.

Antibody-based proteins have become an important class of biologic therapeutics, due in large part to the stability, specificity, and adaptability of the antibody framework. Indeed, antibodies not only have the inherent ability to bind both antigens and endogenous immune receptors; they have also proven extremely amenable to protein engineering. Thus, several derivatives of the monoclonal antibody format, including bispecific antibodies, antibody-drug conjugates, and antibody fragments, have demonstrated efficacy for treating human disease, particularly in the fields of immunology and oncology. Reviewed here are considerations for the design of antibody-based therapeutics, including immunological context, therapeutic mechanisms, and engineering strategies. First, characteristics of antibodies are introduced, with emphasis on structural domains, functionally important receptors, isotypic and allotypic differences, and modifications such as glycosylation. Next, aspects of therapeutic antibody design are discussed, including identification of antigen-specific variable regions, choice of expression system, utility of multispecific formats, and design of antibody derivatives based on fragmentation, oligomerization, or conjugation to other functional moieties. Finally, strategies to enhance antibody function through protein engineering are reviewed while highlighting the impact of fundamental biophysical properties on protein developability.

Drug-induced kidney injury, largely caused by proximal tubular intoxicants, limits development and clinical use of new and approved drugs. Assessing preclinical nephrotoxicity relies on animal models that are frequently insensitive; thus, potentially novel techniques — including human microphysiological systems, or “organs on chips” — are proposed to accelerate drug development and predict safety. Polymyxins are potent antibiotics against multidrug-resistant microorganisms; however, clinical use remains restricted because of high risk of nephrotoxicity and limited understanding of toxicological mechanisms. To mitigate risks, structural analogs of polymyxins (NAB739 and NAB741) are currently in clinical development. Using a microphysiological system to model human kidney proximal tubule, we exposed cells to polymyxin B (PMB) and observed significant increases of injury signals, including kidney injury molecule-1 KIM-1and a panel of injury-associated miRNAs (each P < 0.001). Surprisingly, transcriptional profiling identified cholesterol biosynthesis as the primary cellular pathway induced by PMB (P = 1.22 ×10–16), and effluent cholesterol concentrations were significantly increased after exposure (P < 0.01). Additionally, we observed no upregulation of the nuclear factor (erythroid derived-2)–like 2 pathway, despite this being a common pathway upregulated in response to proximal tubule toxicants. In contrast with PMB exposure, minimal changes in gene expression, injury biomarkers, and cholesterol concentrations were observed in response to NAB739 and NAB741. Our findings demonstrate the preclinical safety of NAB739 and NAB741 and reveal cholesterol biosynthesis as a potentially novel pathway for PMB-induced injury. To our knowledge, this is the first demonstration of a human-on-chip platform used for simultaneous safety testing of new chemical entities and defining unique toxicological pathway responses of an FDA-approved molecule.

The myeloablative agent busulfan (1,4-butanediol dimethanesulfonate) is an old drug that is used routinely to eliminate cancerous bone marrow prior to hematopoietic stem cell transplant. The myeloablative activity and systemic toxicity of busulfan have been ascribed to its ability to cross-link DNA. In contrast, here we demonstrate that incubation of busulfan with the thiol redox proteins glutaredoxin or thioredoxin at pH 7.4 and 37 °C results in the formation of putative S-tetrahydrothiophenium adducts at their catalytic Cys residues, followed by β-elimination to yield dehydroalanine. Both proteins contain a second Cys, in their catalytic C-X-X-C motif, which reacts with the dehydroalanine, the initial Cys adduct with busulfan, or the S-tetrahydrothiophenium, to form novel intramolecular cross-links. The reactivity of the dehydroalanine (DHA) formed is further demonstrated by adduction with glutathione to yield a lanthionine and by a novel reaction with the reducing agent tris(2-carboxyethyl)phosphine (TCEP), which yields a phosphine adduct via Michael addition to the DHA. Formation of a second quaternary organophosphonium salt via nucleophilic substitution with TCEP on the initial busulfan-protein adduct or on the THT(+)-Redoxin species is also observed. These results reveal a rich potential for reactions of busulfan with proteins in vitro, and likely in vivo. It is striking that several of the chemically altered protein products retain none of the atoms of busulfan, in contrast to typical drug-protein adducts or traditional protein modification reagents. In particular, the ability of a clinically used drug to convert Cys to dehydrolanine in intact proteins, and its subsequent reaction with biological thiols, is unprecedented.

The Mycobacterium tuberculosis (Mtb) genome encodes 20 different cytochrome P450 enzymes (CYPs), many of which serve essential biosynthetic roles. CYP51B1, the Mtb version of eukaryotic sterol demethylase, remains a potential therapeutic target. The binding of three drug fragments containing nitrogen heterocycles to CYP51B1 is studied here by continuous wave (CW) and pulsed electron paramagnetic resonance (EPR) techniques to determine how each drug fragment binds to the heme active-site. All three drug fragments form a mixture of complexes, some of which retain the axial water ligand from the resting state. Hyperfine sublevel correlation spectroscopy (HYSCORE) and electron-nuclear double resonance spectroscopy (ENDOR) observe protons of the axial water and on the drug fragments that reveal drug binding modes. Binding in CYP51B1 is complicated by the presence of multiple binding modes that coexist in the same solution. These results aid our understanding of CYP-inhibitor interactions and will help guide future inhibitor design.

Hydrogen deuterium exchange mass spectrometry (H/DX MS) provides a quantitative comparison of the relative rates of exchange of amide protons for solvent deuterons. In turn, the rate of amide exchange depends on a complex combination of the stability of local secondary structure, solvent accessibility, and dynamics. H/DX MS has, therefore, been widely used to probe structure and function of soluble proteins, but its application to membrane proteins was limited previously to detergent solubilized samples. The large excess of lipids from model membranes, or from membrane fractions derived from in vivo samples, presents challenges with mass spectrometry. The lipid nanodisc platform, consisting of apolipoprotein A-derived membrane scaffold proteins, provides a native like membrane environment in which to capture analyte membrane proteins with a well defined, and low, ratio of lipid to protein. Membrane proteins in lipid nanodiscs are amenable to H/DX MS, and this is expected to lead to a rapid increase in the number of membrane proteins subjected to this analysis. Here we review the few literature examples of the application of H/DX MS to membrane proteins in nanodiscs. The incremental improvements in the experimental workflow of the H/DX MS are described and potential applications of this approach to study membrane proteins are described.

The serum half-life and clearance of therapeutic monoclonal antibodies (mAbs) are critical factors that impact their efficacy and optimal dosing regimen. The pH-dependent binding of an mAb to the neonatal Fc receptor (FcRn) has long been recognized as an important determinant of its pharmacokinetics. However, FcRn affinity alone is not a reliable predictor of mAb half-life, suggesting that other biologic or biophysical mechanisms must be accounted for. mAb thermal stability, which reflects its unfolding and aggregation propensities, may also relate to its pharmacokinetic properties. However, no rigorous statistical regression methods have been used to identify combinations of physical parameters that best predict biologic properties. In this work, a panel of eight mAbs with published human pharmacokinetic data were selected for biophysical analyses of FcRn binding and thermal stability. Biolayer interferometry was used to characterize FcRn/mAb binding at acidic and neutral pH, while differential scanning calorimetry was used to determine thermodynamic unfolding parameters. Individual binding or stability parameters were generally weakly correlated with half-life and clearance values. Least absolute shrinkage and selection operator regression was used to identify the combination of two parameters with the best correlation to half-life and clearance as being the FcRn binding response at pH 7.0 and the change in heat capacity. Leave-one-out subsampling yielded a root mean square difference between observed and predicted half-life of just 2.7 days (16%). Thus, the incorporation of multiple biophysical parameters into a cohesive model may facilitate early-stage prediction of in vivo half-life and clearance based on simple in vitro experiments.

P-glycoprotein (P-gp) is a highly substrate-promiscuous efflux transporter that plays a critical role in drug disposition. P-gp utilizes ATP hydrolysis by nucleotide-binding domains (NBDs) to drive transitions between inward-facing (IF) conformations that bind drugs and outward-facing (OF) conformations that release them to the extracellular solution. However, the details of the protein dynamics within either macroscopic IF or OF conformation remain uncharacterized, and the functional role of local dynamics has not been determined. In this work we measured the local dynamics of the IF state of P-gp in lipid nanodiscs and in detergent solution by hydrogen-deuterium (H/D) exchange MS. We observed "EX1 exchange kinetics," or bimodal kinetics, for several peptides distributed in both NBDs, particularly for P-gp in the lipid nanodiscs. Remarkably, the EX1 kinetics occurred on several time scales, ranging from seconds to hours, suggesting highly complex, and correlated, motions. The results indicate at least three distinct conformational states in the ligand-free P-gp and suggest a rough conformational landscape. Addition of excess ATP and vanadate, to favor the OF conformations, caused a generalized, but modest, decrease in H/D exchange throughout the NBDs and slowed the EX1 kinetic transitions of several peptides. The functional implications of the results are consistent with the possibility that conformational selection provides a source of substrate promiscuity.

Cytosolic glutathione transferases (GSTs) comprise a large family of enzymes with canonical structures that diverge functionally and structurally among mammals, invertebrates and plants. Whereas mammalian GSTs have been characterized extensively with regard to their structure and function, invertebrate GSTs remain relatively unstudied. The invertebrate GSTs do, however, represent potentially important drug targets for infectious diseases and agricultural applications. In addition, it is essential to fully understand the structure and function of invertebrate GSTs, which play important roles in basic biological processes. Invertebrates harbor delta- and epsilon-class GSTs, which are not found in other organisms.Drosophila melanogasterGSTs (DmGSTs) are likely to contribute to detoxication or antioxidative stress during development, but they have not been fully characterized. Here, the structures of two epsilon-class GSTs fromDrosophila, DmGSTE6 and DmGSTE7, are reported at 2.1 and 1.5 Å resolution, respectively, and are compared with other GSTs to identify structural features that might correlate with their biological functions. The structures of DmGSTE6 and DmGSTE7 are remarkably similar; the structures do not reveal obvious sources of the minor functional differences that have been observed. The main structural difference between the epsilon- and delta-class GSTs is the longer helix (A8) at the C-termini of the epsilon-class enzymes.

P-glycoprotein (P-gp) is a member of the ABC transporter family that confers drug resistance to many tumors by catalyzing their efflux, and it is a major component of drug-drug interactions. P-gp couples drug efflux with ATP hydrolysis by coordinating conformational changes in the drug binding sites with the hydrolysis of ATP and release of ADP. To understand the relative rates of the chemical step for hydrolysis and the conformational changes that follow it, we exploited isotope exchange methods to determine the extent to which the ATP hydrolysis step is reversible. With γ(18)O4-labeled ATP, no positional isotope exchange is detectable at the bridging β-phosphorus-O-γ-phosphorus bond. Furthermore, the phosphate derived from hydrolysis includes a constant ratio of three (18)O/two (18)O/one (18)O that reflects the isotopic composition of the starting ATP in multiple experiments. Thus, H2O-exchange with HPO4(2-) (Pi) was negligible, suggesting that a [P-gp·ADP·Pi] is not long-lived. This further demonstrates that the hydrolysis is essentially irreversible in the active site. These mechanistic details of ATP hydrolysis are consistent with a very fast conformational change immediately following, or concomitant with, hydrolysis of the γ-phosphate linkage that ensures a high commitment to catalysis in both drug-free and drug-bound states.

The role of conformational ensembles in enzymatic reactions remains unclear. Discussion concerning "induced fit" versus "conformational selection" has, however, ignored detoxication enzymes, which exhibit catalytic promiscuity. These enzymes dominate drug metabolism and determine drug-drug interactions. The detoxication enzyme glutathione transferase A1-1 (GSTA1-1), exploits a molten globule-like active site to achieve remarkable catalytic promiscuity wherein the substrate-free conformational ensemble is broad with barrierless transitions between states. A quantitative index of catalytic promiscuity is used to compare engineered variants of GSTA1-1 and the catalytic promiscuity correlates strongly with characteristics of the thermodynamic partition function, for the substrate-free enzymes. Access to chemically disparate transition states is encoded by the substrate-free conformational ensemble. Pre-steady state catalytic data confirm an extension of the conformational selection model, wherein different substrates select different starting conformations. The kinetic liability of the conformational breadth is minimized by a smooth landscape. We propose that "local" molten globule behavior optimizes detoxication enzymes.

Cytochrome P4503A4 (CYP3A4) is a peripheral membrane protein that plays a major role in enzymatic detoxification of many drugs and toxins. CYP3A4 has an integral membrane N-terminal helix and a localized patch comprised of the G′ and F′ helix regions that are embedded in the membrane, but the effects of membrane composition on CYP3A4 function are unknown. Here, circular dichroism and differential scanning calorimetry were used to compare the stability of CYP3A4 in lipid bilayer nanodiscs with varying ratios of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine to 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). These lipids differ in the acyl-chain length and their degree of unsaturation. The thermal denaturation of CYP3A4 in nanodiscs occurs in a temperature range distinct from that of the nanodisc denaturation so it can be monitored calorimetrically. Melting temperatures (Tm), heat capacities (ΔCp), and calorimetric enthalpies (ΔHcal) for denaturation of CYP3A4 each increased with an increasing fraction of DMPC, with a maximum at 50% DMPC, before decreasing at 75% DMPC. Addition of the inhibitor ketoconazole results in increased thermal stability, and larger ΔCp and ΔHcal values, with different sensitivities to lipid composition. Effects of lipid composition on ligand binding dynamics were also studied. Equilibrium binding affinities of both ketoconazole (KTZ) and testosterone (TST) were minimally affected by lipid composition. However, stopped-flow analyses indicate that the rates of KTZ binding reach a maximum in membranes containing 50% DMPC, whereas the rate of TST binding decreases continuously with an increasing DMPC concentration. These results indicate that CYP3A4 is highly sensitive to the acyl-chain composition of the lipids and fluidity of the membrane in which it is embedded.

Membrane-bound cytochrome P4503A4 (CYP3A4) is the major source of enzymatic drug metabolism. Although several structural models of CYP3A4 in various ligand complexes are available, none includes a lipid bilayer. Details of the effects of the membrane on protein dynamics and solvation, and access channels for ligands, remain uncertain. H/D exchange mass spectrometry (H/DXMS) with ligand free CYP3A4 containing a deletion of residues 3-12, compared to that of the full length wild type, in lipid nanodiscs afforded 91% sequence coverage. Deuterium exchange was fast in the F- and G-helices, HI loop, and C-terminal loop. In contrast, there is very low exchange in the F'- and G'-helices. The results are consistent with the overall membrane orientation of CYP3A4 suggested by published MD simulations and spectroscopic results, and the solvent accessibility of the F/G loop suggests that it is not deeply membrane-embedded. Addition of ketoconazole results in only modest, but global, changes in solvent accessibility. Interestingly, with ketoconazole bound some peptides become less solvent accessible or dynamic, including the F- and G-helices, but several peptides demonstrate modestly increased accessibility. Differential scanning calorimetry (DSC) of CYP3A4-nanodiscs suggests membrane-induced stabilization compared to that of aggregated CYP3A4 in buffer, and this stabilization is enhanced upon addition of the ligand ketoconazole. This ligand-induced stabilization is accompanied by a very large increase in ΔH for CYP3A4 denaturation in nanodiscs, possibly due to increased CYP3A4-membrane interactions. Together, the results suggest a distinct orientation of CYP3A4 on the lipid membrane, and they highlight likely solvent access channels, which are consistent with several MD simulations.

In comparison to imidazole (IMZ) and 1,2,4-triazole (1,2,4-TRZ), the isosteric 1,2,3-triazole (1,2,3-TRZ) is unrepresented among cytochrome P450 (CYP) inhibitors. This is surprising because 1,2,3-TRZs are easily obtained via "click" chemistry. To understand this underrepresentation of 1,2,3-TRZs among CYP inhibitors, thermodynamic and density functional theory computational studies were performed with unsubstituted IMZ, 1,2,4-TRZ, and 1,2,3-TRZ. The results indicate that the lower affinity of 1,2,3-TRZ for the heme iron includes a large unfavorable entropy term likely originating in solvent-1,2,3-TRZ interactions; the difference is not solely due to differences in the enthalpy of heme-ligand interactions. In addition, the 1,2,3-TRZ fragment was incorporated into a well-established CYP3A4 substrate and mechanism-based inactivator, 17-α-ethynylestradiol (17EE), via click chemistry. This derivative, 17-click, yielded optical spectra consistent with low-spin ferric heme iron (type II) in contrast to 17EE, which yields a high-spin complex (type I). Furthermore, the rate of CYP3A4-mediated metabolism of 17-click was comparable to that of 17EE, with a different regioselectivity. Surprisingly, continuous-wave electron paramagnetic resonance (EPR) and HYSCORE EPR spectroscopy indicate that 17-click does not displace water from the sixth axial ligand position of CYP3A4 as expected for a type II ligand. We propose a binding model in which 17-click pendant 1,2,3-TRZ hydrogen bonds with the sixth axial water ligand. The results demonstrate the potential for 1,2,3-TRZ to form metabolically labile water-bridged low-spin heme complexes, consistent with recent evidence that nitrogenous type II ligands of CYPs can be efficiently metabolized. The specific case of [CYP3A4·17-click] highlights the risk of interpreting CYP-ligand complex structure on the basis of optical spectra.

Enzymological paradigms have shifted recently to acknowledge the biological importance of catalytic promiscuity. However, catalytic promiscuity is a poorly understood property, and no thermodynamic treatment has described the conformational landscape of promiscuous versus substrate-specific enzymes. Here, two structurally similar glutathione transferase (GST, glutathione S-transferase) isoforms with high specificity or high promiscuity are compared. Differential scanning calorimetry (DSC) indicates a reversible low temperature transition for the promiscuous GSTA1-1 that is not observed with substrate-specific GSTA4-4. This transition is assigned to rearrangement of the C terminus at the active site of GSTA1-1 based on the effects of ligands and mutations. Near-UV and far-UV circular dichroism indicate that this transition is due to repacking of tertiary contacts with the remainder of the subunit, rather than “unfolding” of the C terminus per se. Analysis of the DSC data using a modified Landau theory indicates that the local conformational landscape of the active site of GSTA1-1 is smooth, with barrierless transitions between states. The partition function of the C-terminal states is a broad unimodal distribution at all temperatures within this DSC transition. In contrast, the remainder of the GSTA1-1 subunit and the GSTA4-4 protein exhibit folded and unfolded macrostates with a significant energy barrier separating them. Their partition function includes a sharp unimodal distribution of states only at temperatures that yield either folded or unfolded macrostates. At intermediate temperatures the partition function includes a bimodal distribution. The barrierless rearrangement of the GSTA1-1 active site within a local smooth energy landscape suggests a thermodynamic basis for catalytic promiscuity.

The allosteric effect of fluconazole (effector) on the formation of 1'-hydroxymidazolam (1'-OH-MDZ) and 4-hydroxymidazolam (4-OH-MDZ) from midazolam (MDZ), a substrate of CYP3A4/5--members of the cytochrome P450 superfamily of enzymes--was examined in healthy volunteers. Following pretreatment with fluconazole, the ratio of the areas under the curve (AUCs) for 4-OH-MDZ and MDZ (AUC(4-OH)/AUC(MDZ)) increased by 35-62%, whereas the ratio AUC(1'-OH)/AUC(MDZ) decreased by 5-37%; the ratio AUC(1'-OH)/AUC(4-OH) decreased by 46-58% after fluconazole administration and had no association with the CYP3A5 genotype. The in vitro formation of 1'-OH-MDZ was more susceptible to inhibition by fluconazole than that of 4-OH-MDZ. Fluconazole decreased the intrinsic formation-clearance ratio of 1'-OH-MDZ/4-OH-MDZ to an extent that was quantitatively comparable to in vivo observations. The elimination clearance of MDZ metabolites appeared unaffected by fluconazole. This study demonstrated that fluconazole alters formation of MDZ metabolites, both in vivo and in vitro, in a manner consistent with an allosteric interaction. The 1'-OH-MDZ/4-OH-MDZ ratio may serve as a biomarker of such interactions among MDZ, CYP3A4/5, and other putative effectors.

The human ATP-binding cassette (ABC) transporter, P-glycoprotein (P-gp; ABCB1), mediates the ATP-dependent efflux of a variety of drugs. As a result, P-gp plays a critical role in tumor cell drug resistance and the pharmacokinetic properties of most drugs. P-gp exhibits extraordinary substrate and inhibitor promiscuity, resulting in a wide range of possible drug-drug interactions. Inhibitory antibodies have long been considered as a possible strategy to modulate P-gp-dependent cancer cell drug resistance, and it is widely suggested that the antibodies MRK16 and UIC2 inhibit P-gp by capturing a single isoform and preventing flux through the catalytic cycle. Although the crystal structures of many bacterial whole transporters, as well as isolated nucleotide-binding domains, have been solved, high resolution structural data for mammalian ABC transporters are currently lacking. It has been extremely difficult to determine the detailed mechanism of transport of P-gp, in part because it is difficult to obtain purified protein in well defined lipid systems. Here we exploit surface plasmon resonance (SPR) to probe conformational changes associated with these intermediate states for P-gp in lipid bilayer nanodiscs. The results indicate that P-gp in nanodiscs undergoes functionally relevant ligand-dependent conformational changes and that previously described inhibitory antibodies bind to multiple nucleotide-bound states but not the ADP-VO(4)-trapped state, which mimics the post-hydrolysis state. The results also suggest that the substrate drug vinblastine is released at stages that precede or follow the post-hydrolysis ADP-PO(4)·P-gp complex.

Cytochrome P450 3A4 (CYP3A4) is the dominant xenobiotic metabolizing CYP. Despite great interest in CYP enzymology, two in vitro aspects of CYP3A4 catalysis are still not well understood, namely, sequential metabolism and allosteric activation. We have therefore investigated such a system in which both phenomena are present. Here we report that the sequential metabolism of Nile Red (NR) is accelerated by the heterotropic allosteric effector α-naphthoflavone (ANF). ANF increases the rates of formation for NR metabolites M1 and M2 and also perturbs the metabolite ratio in favor of M2. Thus, ANF has as an allosteric effect on a kinetic branch point. Co-incubating deuterium-labeled NR and unlabeled M1, we show that ANF increases k(cat)/k(off) ~1.8-fold in favor of the k(cat) of M2 production. Steady-state metabolic experiments are analyzed using a kinetic model in which the enzyme and substrates are not in rapid equilibrium, and this distinction allows for the estimation of rates of catalysis for the formation of both the primary (M1) and secondary (M2) products, as well as the partitioning of enzyme between these states. These results are compared with those of earlier spectroscopic investigations of NR and ANF cooperativity, and a mechanism of ANF heteroactivation is presented that involves effects on substrate off rate and coupling efficiency.

It has been appreciated for a long time that enzymes exist as conformational ensembles throughout multiple stages of the reactions they catalyze, but there is renewed interest in the functional implications. The energy landscape that results from conformationlly diverse poteins is a complex surface with an energetic topography in multiple dimensions, even at the transition state(s) leading to product formation, and this represents a new paradigm. Nearly simultaneous with the renewed interest in conformational ensembles, a new paradigm concerning enzyme function has grown, wherein catalytic promiscuity has clear biological advantages in some cases. ‘Useful’, or biologically functional, promiscuity or the related behavior of ‘multifunctionality’, can be found in the immune system, enzymatic detoxification, signal transduction, and in the evolution of new function from an existing pool of folded protein scaffolds. Experimental evidence supports the widely held assumption that conformational heterogeneity promotes functional promiscuity. The common link between these co-evolving paradigms is the inherent structural plasticity and conformational dynamics of proteins that, on one hand, leads to complex but evolutionarily selected energy landscapes, and on the other hand promotes functional promiscuity. Here we consider a logical extension of the overlap between these two nascent paradigms: functionally promiscuous and multifunctional enzymes such as detoxification enzymes are expected to have an ensemble landscape with a greater number of states accessible on multiple time scales than substrate specific enzymes. Two attributes of detoxification enzymes become important in the context of conformational ensembles: these enzymes metabolize multiple substrates, often in substrate mixtures, and they can form multiple products from a single substrate. These properties, combined with complex conformational landscapes, lead to the possibility of interesting time-dependent, or emergent, properties. Here we demonstrate these properties with kinetic simulations of non-equilibrium steady state (NESS) behavior resulting from energy landscapes expected for detoxification enzymes. Analogous scenarios with other promiscuous enzymes may be worthy of consideration.

Electrophilic products of lipid peroxidation are important contributors to the progression of several pathological states. The prototypical α,β-unsaturated aldehyde, 4-hydroxynonenal (HNE), triggers cellular events associated with oxidative stress, which can be curtailed by the glutathione-dependent elimination of HNE. The glutathione transferases (GSTs) are a major determinate of the intracellular concentration of HNE and can influence susceptibility to toxic effects, particularly when HNE and GST levels are altered in disease states. In this article, we provide a brief summary of the cellular effects of HNE, followed by a review of its GST-catalyzed detoxification, with an emphasis on the structural attributes that play an important role in the interactions with alpha-class GSTs. Some of the key determining characteristics that impart high alkenal activity reside in the unique C-terminal interactions of the GSTA4-4 enzyme. Studies encompassing both kinetic and structural analyses of related isoforms will be highlighted, with additional attention to stereochemical aspects that demonstrate the capacity of GSTA4-4 to detoxify both enantiomers of the biologically relevant racemic mixture while generating a select set of diastereomeric products with subsequent implications. A summary of the literature that examines the interplay between GSTs and HNE in model systems relevant to oxidative stress will also be discussed to demonstrate the magnitude of importance of GSTs in the overall detoxification scheme.

Promiscuity at the molecular level plays a critical role in several cellular processes including signal transduction, immunity, and detoxification (1–9). More broadly, the evolution of new enzymatic activities from a pool of existing protein scaffolds, whether in nature or in vitro directed evolution, may proceed through ‘promiscuous’ intermediates that catalyze multiple reactions (10–19). A common description of this possible evolutionary process and relevant factors is schematized in Figure 1, wherein a specific enzyme undergoes mutation to yield a promiscuous variant that catalyzes a new reaction, without forfeiting its original activity (11, 20). At early stages along the evolutionary trajectory the mutation is neutral with respect to the original function. If this second enzymatic activity provides any survival advantage, then it is retained and, after gene duplication, one copy of the gene encoding the promiscuous enzyme can undergo further refinement to optimize the new function, without loss of the original function. It may be useful to distinguish between ‘promiscuous’ enzymes with catalytic activities unrelated to their presumed biological function, and ‘multifunctional’ enzymes with a clear functional advantage in accepting multiple substrates, as suggested (20). Here we use the term ‘promiscuity’ for both cases, because we aim to understand the physical and chemical properties of enzymes that process multiple substrates, regardless of the functional purpose. Recent studies have suggested that promiscuous intermediates easily evolve into efficient and specific enzymes (11, 20). While this model provides an extremely useful framework to conceptualize evolutionary structure-function relationships, ‘promiscuity’ is not a well-understood property, and it is difficult to quantify, although progress has been made (17, 19). Figure 1 A. Promiscuous Intermediate Model for Evolution of New Enzymatic Function. Enzymes with high catalytic specificity for one reaction (original catalytic activity) move along the indicated trajectory. At early stages they become promiscuous, gaining new ... For example, natural evolution of new function does not happen in isolation. Rather, each evolutionary intermediate must function in the context of the existing niche that includes other enzymes and a ‘background’ of potential ligands. Within a niche, a loss due to mutation in any useful trait is referred to as a ‘negative trade-off’(20). The evolutionary trajectory of an enzyme is unlikely to be a two dimensional process as depicted in the model in Figure 1a, but it will trace a multidimensional surface that includes the relationships between promiscuity and other properties that comprise the overall fitness of the enzyme. That is, there will be a third collective coordinate for negative trade-offs, as in Figure 1b. In this case, ‘promiscuity’ per se is not an obvious advantage in the evolutionary process unless it is achieved without significant negative trade-offs (20). Whereas a highly promiscuous enzyme would provide a very efficient template from which to evolve many new functions from a single protein sequence, such an intermediate would unlikely be neutral with respect to all of the original properties. If the magnitude of any negative trade-off is coupled exactly to the magnitude of promiscuity, then the mutational neutrality will not be maintained in all dimensions of fitness space and the path towards new function will have a different shape than if the putative negative trade-off is not correlated with the promiscuity in any dimension. Ultimately, the relationship between any potential negative trade-off and catalytic promiscuity will determine whether evolution will follow that path. These concepts underscore the need to define the relationship between catalytic promiscuity and other protein properties. Speculatively, for example, a new promiscuous enzyme would have increased susceptibility to inhibition by ligands that do not inhibit the ‘original’ enzyme and this would be a negative trade-off. This is based on the observation that structurally diverse drugs can bind to and inhibit a single enzyme, without yielding products (consider antibiotics, statins, etc.). In contrast catalysis requires proper alignment of multiple reactive functional groups. Each incremental step towards promiscuity, i.e. towards the diagonal of Figure 1a, will only be advantageous if the new catalytic activity is not accompanied by an equal or greater increase in susceptibility to ‘new’ inhibitory interactions that, at ambient ligand concentrations, limit the original enzyme activity. For the example here, inhibitory promiscuity is a potential negative trade-off that adds a new dimension to yield the scheme in Figure 1b, and can be substituted by any parameter hypothesized to be a negative trade-off. We emphasize that the current studies are not intended to establish specific evolutionary mechanisms, which may best be accomplished by successive rounds of random mutagenesis and relevant selection protocols. Rather we aim to understand more generally the relationships between catalytic promiscuity and other properties of enzymes that could limit or define evolutionary pathways. An understanding of these relationships is necessary to understand the role of catalytic promiscuity in evolution and to advance the incomplete model of Figure 1. The Cytochrome P450s provide an excellent test case to further understand these relationships because, as a family with highly conserved scaffold, they span a very wide range of catalytic promiscuity. Human hepatic CYPs dominate drug metabolism and have very broad substrate selectivities (19, 21). Here we use these multifunctional CYPs as a model for physico-chemical traits of enzymes near the limits of catalytic promiscuity. Other CYPs, including bacterial and human isoforms used here, contribute to specific biosynthetic pathways and are very substrate specific (22, 23). Collectively, the CYP superfamily provides insight into the relationship between catalytic and inhibitory promiscuity. The results indicate a surprising level of inhibitory promiscuity even for substrate-specific CYPs, and they further suggest that because inhibitory promiscuity is already present in highly specific enzymes, there is little cost associated with evolutionary steps toward catalytic promiscuity.

Cytochrome P450s (CYPs) are heme-containing monooxygenases that contribute to an enormous range of enzymatic function including biosynthetic and detoxification roles. This review summarizes recent studies concerning interactions of CYPs with ligands including substrates, inhibitors, and diatomic heme-ligating molecules. These studies highlight the complexity in the relationship between the heme spin state and active site occupancy, the roles of water in directing protein-ligand and ligand-heme interactions, and the details of interactions between heme and gaseous diatomic CYP ligands. Both kinetic and thermodynamic aspects of ligand binding are considered.

Inhibitors of anthrax lethal factor (LF) are currently being sought as effective therapeutics for the treatment of anthrax. Here we report a novel screening approach for inhibitors of LF, a yeast-hybrid-based assay system in which the expression of reporter genes from a Gal4 promoter is repressed by LF proteolytic activity. Yeast cells were co-transformed with LF and a chimeric transcription factor that contains an LF substrate sequence inserted between the DNA-binding and activation domains of Gal4. In the resulting yeast cells, LF cleaves the substrate, thus inactivating the chimeric Gal4 and resulting in lack of expression of reporter genes. Compounds that inhibit LF cleavage of its substrate are identified by changes in reporter gene activity. Relative to in vitro screens for inhibitors of LF proteolytic activity, this screen has the advantage of excluding compounds that are toxic or non-permeable to eukaryotic cells. Additionally, the screen has the advantage of being fast, easy and cheap because exogenous LF and substrate are not needed. An initial chemical library screen with this system has identified four candidate inhibitors which were confirmed to inhibit LF protease activity in an in vitro assay. Furthermore, FBS-00831, one of the compounds identified, protects Raw 264.7 macrophages from anthrax lethal toxin and the possible binding site on LF was also evaluated by molecular docking.

A previous study (i.e. Dickmann, L., et al. (2004) Mol. Pharmacol. 65, 842-850.) revealed some unusual properties of the R108H mutant of cytochrome P450 2C9 (CYP2C9), including elevated thermostability relative to CYP2C9, as well as a UV-visible absorbance spectrum that was indicative of nitrogenous ligation to the heme iron. In the present study, size-exclusion chromatography and UV-visible absorbance spectroscopy of CYP2C9 R108H monomers demonstrated that nitrogen ligation is indeed intramolecular. Pulsed electron paramagnetic resonance of CYP2C9 R108H monomers showed that a histidine is most likely bound to the heme as previously hypothesized. An energy-minimized model of the R108H mutant maintained a CYP fold, despite substantial movement of several loop regions of the mutant and, therefore, represents an extreme example of a closed conformation of the enzyme. Molecular dynamics (MD) simulations of CYP2C9 were performed in order to study the range of energetically accessible CYP2C9 conformations. These in silico studies showed that the B-C loop region of CYP2C9 moves away from the heme to a position resembling the putative open conformation described for rabbit CYP2B4. A model involving the movement of the B-C loop region and R108 between the open and closed conformations of CYP2C9 is presented, which helps to explain the enzyme’s ability to regio- and stereo-specifically metabolize some ligands, while allosterically activating others.

The cytosolic glutathione transferases (GSTs)1 comprise a family of multifunctional enzymes that share similar overall topologies, dimer interactions, and glutathione- (GSH) binding sites, yet retain the capacity to detoxify a broad spectrum of electrophiles (1-3). All members are functionally dimeric proteins that contain one active site per subunit that consists of a GSH- and substrate-binding site, referred to as the G-site and H-site, respectively. The conserved G-site is located in the N-terminal domain while the H-site is located between the bundle of α–helices in the C-terminal domain, variations of which allow GSTs to accommodate structurally diverse electrophiles. The alpha class GSTs, despite strong intra-class sequence similarity, are highly differentiated by their substrate specificities. While GSTA1-1 is regarded as a highly promiscuous enzyme with catalytic activity towards structurally diverse substrates, GSTA4-4 is recognized for its specificity toward alkenals such as 4-hydroxynonenal (HNE), an electrophile that has been described as one of the most toxic end products of oxidative stress (4-7). Both of these alpha class GSTs utilize three substrate recognition modules distributed throughout the primary structure to construct the H-site in the folded protein. The β1-α1 region, the end of the α4-helix, and the C-terminus all contribute to some extent to the H-site. However, the unique characteristics of the C-terminal α9-helix, a distinguishing feature of alpha class GSTs, have been the subject of numerous studies (5, 6, 8-20) and play a central role in the different specificities achieved within the constraints of a conserved scaffold (5, 6, 21, 22). While the position of the C-terminus of GSTA1-1 is highly heterogeneous, and converts to a more ordered helix that completes the H-site in a ligand-dependent manner, the well-defined α9-helix of GSTA4-4 creates a pre-existing H-site pocket shaped to accommodate alkenals (5). Aromatic-aromatic interactions can contribute to protein stabilization (23) and the human GSTA4-4 (hGSTA4-4) enzyme contains an edge-to-face aromatic-aromatic interaction between Phe111 in the α4-helix and Tyr217 in the α9-helix that is absent within the H-site region of hGSTA1-1. This interaction was previously hypothesized to contribute to a closer association between the α4-α5 helix-turn-helix “tower” and C-terminal regions of hGSTA4-4 that helps to pre-organize the C-terminus for alkenal binding (6). Incorporation of this interaction into hGSTA1-1 resulted in a tower mutant with a less dynamic C-terminus and increased selectivity for HNE, as well as stereoselectivity of product formation more reminiscent of hGSTA4-4 (24). This result also demonstrated that the conformational heterogeneity is highly coupled throughout the structure of the promiscuous hGSTA1-1 enzyme (6). Although the dynamics have yet to be explicitly characterized, a more extensive redesign of hGSTA1-1, made to mimic first-sphere interactions utilized by hGSTA4-4, was found to drastically alter the specificity of hGSTA1-1 from aromatic substitution reactions to Michael-like additions in the GIMFhelix mutant described by Nilsson et al. (22). While it is Tyr212 in the α9-helix that is proposed to interact with the aldehyde of alkenals, this area is occupied by the Cβ of Ala12 in hGSTA1-1 and the lack of a side chain at residue 12 is known to be required for high catalytic efficiency with alkenals (5, 22, 25). Residues lining the hydrophobic binding cavity in hGSTA4-4 were also selected for incorporation in the redesign of hGSTA1-1 (22). These include Ile107, Met108, and Phe111 in the α4-helix as well as Val216, Tyr217, and the aforementioned Tyr212. Since a number of the targeted residues comprise the C-terminus, which is known to be static and positioned to contribute the H-site in hGSTA4-4, the entire C-terminus (residues 208-222) was exchanged with that of hGSTA4-4. As a result, the final GIMFhelix mutant contains both residues necessary to potentially establish the tower region interaction, in addition to several other substitutions aimed to confer catalytic activity like that of hGSTA4-4. In fact the GIMFhelix mutant provides an excellent example of rational enzyme design and a model to further understand the structural and dynamic basis for the functional differences between hGSTA1-1 and hGSTA4-4. This redesign of hGSTA1-1, with only 6% of the residues mutated to key residues present in hGSTA4-4, resulted in 20- and > 300-fold increases in catalytic efficiency with HNE and nonenal, respectively (21, 22). An 11-fold decrease in the catalytic efficiency was observed for the aromatic substitution reaction involving 1-chloro-2,4-dinitrobenzene (CDNB) favored by hGSTA1-1, indicating the changes were also correlated with a general shift in specificity toward Michael-like additions similar to hGSTA4-4 and illustrating the importance of first-sphere interactions for altered specificity. Herein we describe the first crystal structures reported for GIMFhelix that offer a crystallographic look at the β1-α1 loop, α4-helix, and C-terminal region targeted to convey high alkenal activity to an essentially hGSTA1-1 protein scaffold. Comparisons of crystal structures also provide an opportunity to dissect interactions exploited by hGSTA1-1 and A4-4 that led to a unique situation inadvertently created by the Tyr212-Phe10-Phe220 network within the mutant's active site. Overall these structures help to elucidate the structural and conformational basis underlying the specificity and catalytic efficiency emerged in GIMFhelix.

Cytochrome P450 (P450) inhibition often occurs in a strongly substrate- and inhibitor-dependent manner, with a given inhibitor affecting the metabolism of different substrates to differing degrees and with a given substrate responding differently to different inhibitors. Traditionally, patterns of functional similarity and dissimilarity among substrates and inhibitors have been studied using clustering analysis of pair-wise correlation coefficients. Principal component analysis (PCA) is a widely used statistical technique that identifies the globally most significant independent trends in a set of data. Here, we show that PCA can be usefully applied to study the differential effects on a panel of P450 probe substrates by a panel of inhibitors, using published data on CYP3A4 (Kenworthy et al., 1999) and CYP2C9 (Kumar et al., 2006). PCA can detect functional similarities among substrates and inhibitors that are not readily apparent using pair-wise clustering analysis. PCA also allows identification of the functionally typical and atypical substrates that might be used in combination to fully explore the P450 functional landscape.

The contribution of ligand dynamics to CYP allosterism has not been considered in detail. On the basis of a previous study, we hypothesized that CYP2A6 and CYP2E1 accommodate multiple xylene ligands. As a result, the intramolecular ( k H/ k D) obs values observed for some xylene isomers are expected to be dependent on ligand concentration with contributions from [CYP.xylene] and [CYP.xylene.xylene], etc. To explore this possibility and the utility of kinetic isotope effects in characterizing allosteric CYP behavior, steady state kinetics, product ratios, and ( k H/ k D) obs values for CYP2E1 and CYP2A6 oxidation of m-xylene-alpha- (2)H 3 and p-xylene-alpha- (2)H 3 were determined. Evidence is presented that CYP2A6 accommodates multiple ligands and that intramolecular isotope effect experiments can provide insight into the mechanisms of multiple-ligand binding. CYP2A6 exhibited cooperative kinetics for m-xylene-alpha- (2)H 3 oxidation and a concentration-dependent decrease in the m-methylbenzylalcohol:2,4-dimethylphenol product ratio (9.8 +/- 0.1 and 4.8 +/- 0.3 at 2.5 microM and 1 mM, respectively). Heterotropic effects were observed as well, as incubations containing both 15 microM m-xylene-alpha- (2)H 3 and 200 microM p-xylene resulted in further reduction of the product ratio (2.4 +/- 0.2). When p-xylene (60 microM) was replaced with deuterium-labeled d 6- p-xylene (60 microM), an intermolecular competitive inverse isotope effect on 2,4-dimethylphenol formation [( k H/ k D) obs = 0.49] was observed, indicating that p-xylene exerts heterotropic effects by residing in the active site simultaneously with m-xylene. The data indicate that there is a concentration-dependent decrease in the reorientation rate of m-xylene, as no increase in ( k H/ k D) obs was observed in the presence of an increased level of metabolic switching. That is, the accommodation of a second xylene molecule in the active site leads to a decrease in substrate dynamics.

Cytochrome P450 3A4 (CYP3A4) plays a critical role in the metabolism of many drugs. CYP3A4 exhibits extraordinary substrate promiscuity and unusual allosteric kinetics. In addition, many CYPs catalyze sequential oxidations on a single substrate, but in most cases, mechanistic details of these processes are not well-established. As a result, in vivo clearance of many drugs and their metabolites is difficult to predict on the basis of the complex in vitro kinetics, and new in vitro probes are required to understand these behaviors. The near-IR fluorescent probe Nile Red, which has strong solvatochromic behavior, was investigated as a probe of allostery and sequential metabolism with CYP3A4. Nile Red binds with apparent Kd values of 0.05 and 2.3 muM, based on a sigmoidal dependence of heme spin state on Nile Red concentration, where the first equivalent of Nile Red increased the high-spin fraction by only 13% of the total change at saturation. Mass spectrometry analysis indicates that Nile Red is metabolized sequentially by CYP3A4 to the N-monoethyl and N-desethyl products, confirming that the immediate vicinity of the heme iron is one binding site. In the presence of CYP3A4, steady-state fluorescence emission and excitation spectra, as well as excited-state lifetimes at varying Nile Red concentrations, indicate a high-affinity site that modulates the fluorescent properties of Nile Red. The Nile Red binding site is competitively eliminated by itraconazole, which is a high-affinity ligand known to coordinate to the heme iron. Together, the data suggest that Nile Red binds to the active site with high affinity ( approximately 50 nM), where it is desolvated in a low-dielectric environment. In addition, Nile Red is sequentially oxidized at rates comparable to or faster than those of other in vitro probes, which emphasizes its utility in the further examination of this important kinetic phenomenon in vitro.

Catalytic promiscuity is a widespread, but poorly understood, phenomenon among enzymes with particular relevance to the evolution of new functions, drug metabolism, and in vitro biocatalyst engineering. However, there is at present no way to quantitatively measure or compare this important parameter of enzyme function. Here we define a quantitative index of promiscuity (I) that can be calculated from the catalytic efficiencies of an enzyme toward a defined set of substrates. A weighted promiscuity index (J) that accounts for patterns of similarity and dissimilarity among the substrates in the set is also defined. Promiscuity indices were calculated for three different enzyme classes: eight serine and cysteine proteases, two glutathione S-transferase (GST) isoforms, and three cytochrome P450 (CYP) isoforms. The proteases ranged from completely specific (granzyme B, J = 0.00) to highly promiscuous (cruzain, J = 0.83). The four drug-metabolizing enzymes studied (GST A1-1 and the CYP isoforms) were highly promiscuous, with J values between 0.72 and 0.92; GST A4-4, involved in the clearance of lipid peroxidation products, is moderately promiscuous (J = 0.37). Promiscuity indices also allowed for studies of correlation between substrate promiscuity and an enzyme's activity toward its most-favored substrate, for each of the three enzyme classes.

The membrane-bound protein cytochrome P450 3A4 (CYP3A4) is a major drug-metabolizing enzyme. Most studies of ligand binding by CYP3A4 are currently carried out in solution, in the absence of a model membrane. Therefore, there is little information concerning the membrane effects on CYP3A4 ligand binding behavior. Phospholipid bilayer Nanodiscs are a novel model membrane system derived from high density lipoprotein particles, whose stability, monodispersity, and consistency are ensured by their self-assembly. We explore the energetics of four ligands (6-(p-toluidino)-2-naphthalenesulfonic acid (TNS), alpha-naphthoflavone (ANF), miconazole, and bromocriptine) binding to CYP3A4 incorporated into Nanodiscs. Ligand binding to Nanodiscs was monitored by a combination of environment-sensitive ligand fluorescence and ligand-induced shifts in the fluorescence of tryptophan residues present in the scaffold proteins of Nanodiscs; binding to the CYP3A4 active site was monitored by ligand-induced shifts in the heme Soret band absorbance. The dissociation constants for binding to the active site in CYP3A4-Nanodiscs were 4.0 microm for TNS, 5.8 microm for ANF, 0.45 microm for miconazole, and 0.45 microm for bromocriptine. These values are for CYP3A4 incorporated into a lipid bilayer and are therefore presumably more biologically relevant that those measured using CYP3A4 in solution. In some cases, affinity measurements using CYP3A4 in Nanodiscs differ significantly from solution values. We also studied the equilibrium between ligand binding to CYP3A4 and to the membrane. TNS showed no marked preference for either environment; ANF preferentially bound to the membrane, and miconazole and bromocriptine preferentially bound to the CYP3A4 active site.

The DNA alkylating agent busulfan is used to 'precondition' patients with leukemia, lymphomas and other hematological disorders prior to hematopoietic stem cell transplants. Busulfan is metabolized via conjugation with glutathione (GSH) followed by intramolecular rearrangement to the GSH analog γ-glutamyl-dehydroalanyl -glycine (EdAG). EdAG contains the electrophilic dehydroalanine, which is expected to react with protein nucleophiles, particularly proteins with GSH binding sites such as glutaredoxins (Grx's). Incubation of EdAG with human Grx-1 or Grx-2 results in facile adduction of cys-23 and cys-77, respectively, as determined by ESI-MS/MS. The resulting modified proteins are catalytically inactive. In contrast, the glutathione transferase A1-1 includes a GSH binding site with a potentially reactive tyrosinate (Tyr-9) but it does not react with EdAG. Similarly, Cys-112 of GSTA1-1, which lies outside the active site and is known to form disulfides with GSH, does not react with EdAG. The results provide the first demonstration of the reactivity of any busulfan metabolites with intact proteins, and they suggest that GSH-binding sites containing thiolates are most susceptible. The adduction of Grx's by EdAG suggests the possible alteration of proteins that are normally regulated via Grx-dependent reversible glutathionylation or deglutathionylation. Dysregulation of Grx-dependent processes could contribute to cellular toxicity of busulfan.

Methods for the initial steps of surface plasmon resonance analysis of membrane proteins incorporated in lipid nanodiscs are described. Several types of Biacore sensor chips are available and require distinct strategies to immobilize proteonanodiscs on the chip surface. The procedures for immobilization on three of these chips (NTA, antibody coupled CM5, and L1) are described in this unit and results are demonstrated for a model system with cytochrome P4503A4 (CYP3A4) in nanodiscs binding to a polyclonal anti-CYP3A4 antibody. Advantages and disadvantages of each chip type are considered.

The human immunodeficiency virus (HIV) is the causative agent of acquired immune deficiency syndrome (AIDS) and is thus responsible for significant morbidity and mortality worldwide. Despite considerable effort, preparation of an effective vaccine for AIDS has been elusive and it has become clear that a fundamental understanding of the relevant antigenic targets on HIV is essential. The Env trimer spike is the only viral antigen present on the surface of the viral particle and it is the target of all broadly neutralizing antibodies isolated to date. Thus, a soluble, homogeneous, and well-defined preparation of Env trimers is an important first step towards biochemical and structural characterization of the antigenic spike. Phospholipid bilayer nanodiscs represent a relatively new technology that can serve as a platform for the assembly of membrane proteins into a native membrane-like environment. Here we describe the preparation and characterization of unprocessed full-length, natively glycoslyated gp160 Env proteins incorporated into nanodiscs (gp160-ND). The particles are soluble and well defined in the absence of detergent, and possess a morphology anticipated of Env incorporated into a lipid ND. Importantly, the gp160-NDs retain CD4 and Env antibody binding characteristics expected of a functional trimer spike and their incorporation into a lipid membrane allows interrogation of epitopes associated with the membrane-proximal ectodomain region of gp41. These studies provide the groundwork for the use of gp160-ND in more detailed biochemical and structural studies that may set the stage for their use in vaccine development.

Azoles and pyridines are commonly incorporated into small molecule inhibitor scaffolds that target cytochromes P450 (CYPs) as a strategy to increase drug binding affinity, impart isoform-dependent selectivity, and improve metabolic stability. Optical absorbance spectra of the CYP-inhibitor complex are widely used to infer whether these inhibitors are ligated directly to the heme iron as catalytically inert, low-spin (type II) complexes. Here, we show that the low-spin complex between a drug-metabolizing CYP2C9 variant and 4-(3-phenylpropyl)-1H-1,2,3-triazole (PPT) retains an axial water ligand despite exhibiting elements of "classic" type II optical behavior. Hydrogens of the axial water ligand are observed by pulsed electron paramagnetic resonance (EPR) spectroscopy for both inhibitor-free and inhibitor-bound species and show that inhibitor binding does not displace the axial water. A (15)N label incorporated into PPT is 0.444 nm from the heme iron, showing that PPT is also in the active site. The reverse type I inhibitor, LP10, of CYP125A1 from Mycobacterium tuberculosis, known from X-ray crystal structures to form a low-spin water-bridged complex, is found by EPR and by visible and near-infrared magnetic circular dichroism spectroscopy to retain the axial water ligand in the complex in solution.