Protein and Drug Design Algorithms Using Improved Biophysical Modeling

This thesis focuses on the development of algorithms that will allow protein design calculations to incorporate more realistic modeling assumptions. Protein design algorithms search large sequence spaces for protein sequences that are biologically and medically useful. Better modeling could improve the chance of success in designs and expand the range of problems to which these algorithms are applied. I have developed algorithms to improve modeling of backbone flexibility (DEEPer) and of more extensive continuous flexibility in general (EPIC and LUTE). I’ve also developed algorithms to perform multistate designs, which account for effects like specificity, with provable guarantees of accuracy (COMETS), and to accommodate a wider range of energy functions in design (EPIC and LUTE).

Anti-AIDS agents 81. Design, synthesis, and structure-activity relationship study of betulinic acid and moronic acid derivatives as potent HIV maturation inhibitors.

In our continuing study of triterpene derivatives as potent anti-HIV agents, different C-3 conformationally restricted betulinic acid (BA, 1) derivatives were designed and synthesized in order to explore the conformational space of the C-3 pharmacophore. 3-O-Monomethylsuccinyl-betulinic acid (MSB) analogues were also designed to better understand the contribution of the C-3' dimethyl group of bevirimat (2), the first-in-class HIV maturation inhibitor, which is currently in phase IIb clinical trials. In addition, another triterpene skeleton, moronic acid (MA, 3), was also employed to study the influence of the backbone and the C-3 modification toward the anti-HIV activity of this compound class. This study enabled us to better understand the structure-activity relationships (SAR) of triterpene-derived anti-HIV agents and led to the design and synthesis of compound 12 (EC(50): 0.0006 microM), which displayed slightly better activity than 2 as a HIV-1 maturation inhibitor.

Economics of new oncology drug development.

PURPOSE: Review existing studies and provide new results on the development, regulatory, and market aspects of new oncology drug development. METHODS: We utilized data from the US Food and Drug Administration (FDA), company surveys, and publicly available commercial business intelligence databases on new oncology drugs approved in the United States and on investigational oncology drugs to estimate average development and regulatory approval times, clinical approval success rates, first-in-class status, and global market diffusion. RESULTS: We found that approved new oncology drugs to have a disproportionately high share of FDA priority review ratings, of orphan drug designations at approval, and of drugs that were granted inclusion in at least one of the FDA's expedited access programs. US regulatory approval times were shorter, on average, for oncology drugs (0.5 years), but US clinical development times were longer on average (1.5 years). Clinical approval success rates were similar for oncology and other drugs, but proportionately more of the oncology failures reached expensive late-stage clinical testing before being abandoned. In relation to other drugs, new oncology drug approvals were more often first-in-class and diffused more widely across important international markets. CONCLUSION: The market success of oncology drugs has induced a substantial amount of investment in oncology drug development in the last decade or so. However, given the great need for further progress, the extent to which efforts to develop new oncology drugs will grow depends on future public-sector investment in basic research, developments in translational medicine, and regulatory reforms that advance drug-development science.

Cytocidal amino acid starvation of Saccharomyces cerevisiae and Candida albicans acetolactate synthase (ilv2{Delta}) mutants is influenced by the carbon source and rapamycin.

The isoleucine and valine biosynthetic enzyme acetolactate synthase (Ilv2p) is an attractive antifungal drug target, since the isoleucine and valine biosynthetic pathway is not present in mammals, Saccharomyces cerevisiae ilv2Delta mutants do not survive in vivo, Cryptococcus neoformans ilv2 mutants are avirulent, and both S. cerevisiae and Cr. neoformans ilv2 mutants die upon isoleucine and valine starvation. To further explore the potential of Ilv2p as an antifungal drug target, we disrupted Candida albicans ILV2, and demonstrated that Ca. albicans ilv2Delta mutants were significantly attenuated in virulence, and were also profoundly starvation-cidal, with a greater than 100-fold reduction in viability after only 4 h of isoleucine and valine starvation. As fungicidal starvation would be advantageous for drug design, we explored the basis of the starvation-cidal phenotype in both S. cerevisiae and Ca. albicans ilv2Delta mutants. Since the mutation of ILV1, required for the first step of isoleucine biosynthesis, did not suppress the ilv2Delta starvation-cidal defects in either species, the cidal phenotype was not due to alpha-ketobutyrate accumulation. We found that starvation for isoleucine alone was more deleterious in Ca. albicans than in S. cerevisiae, and starvation for valine was more deleterious than for isoleucine in both species. Interestingly, while the target of rapamycin (TOR) pathway inhibitor rapamycin further reduced S. cerevisiae ilv2Delta starvation viability, it increased Ca. albicans ilv1Delta and ilv2Delta viability. Furthermore, the recovery from starvation was dependent on the carbon source present during recovery for S. cerevisiae ilv2Delta mutants, reminiscent of isoleucine and valine starvation inducing a viable but non-culturable-like state in this species, while Ca. albicans ilv1Delta and ilv2 Delta viability was influenced by the carbon source present during starvation, supporting a role for glucose wasting in the Ca. albicans cidal phenotype.

Development of universal antidotes to control aptamer activity.

With an ever increasing number of people taking numerous medications, the need to safely administer drugs and limit unintended side effects has never been greater. Antidote control remains the most direct means to counteract acute side effects of drugs, but, unfortunately, it has been challenging and cost prohibitive to generate antidotes for most therapeutic agents. Here we describe the development of a set of antidote molecules that are capable of counteracting the effects of an entire class of therapeutic agents based upon aptamers. These universal antidotes exploit the fact that, when systemically administered, aptamers are the only free extracellular oligonucleotides found in circulation. We show that protein- and polymer-based molecules that capture oligonucleotides can reverse the activity of several aptamers in vitro and counteract aptamer activity in vivo. The availability of universal antidotes to control the activity of any aptamer suggests that aptamers may be a particularly safe class of therapeutics.

Conformational kinetics reveals affinities of protein conformational states.

Most biological reactions rely on interplay between binding and changes in both macromolecular structure and dynamics. Practical understanding of this interplay requires detection of critical intermediates and determination of their binding and conformational characteristics. However, many of these species are only transiently present and they have often been overlooked in mechanistic studies of reactions that couple binding to conformational change. We monitored the kinetics of ligand-induced conformational changes in a small protein using six different ligands. We analyzed the kinetic data to simultaneously determine both binding affinities for the conformational states and the rate constants of conformational change. The approach we used is sufficiently robust to determine the affinities of three conformational states and detect even modest differences in the protein's affinities for relatively similar ligands. Ligand binding favors higher-affinity conformational states by increasing forward conformational rate constants and/or decreasing reverse conformational rate constants. The amounts by which forward rate constants increase and reverse rate constants decrease are proportional to the ratio of affinities of the conformational states. We also show that both the affinity ratio and another parameter, which quantifies the changes in conformational rate constants upon ligand binding, are strong determinants of the mechanism (conformational selection and/or induced fit) of molecular recognition. Our results highlight the utility of analyzing the kinetics of conformational changes to determine affinities that cannot be determined from equilibrium experiments. Most importantly, they demonstrate an inextricable link between conformational dynamics and the binding affinities of conformational states.

Structural and Biochemical Dissection of the Trehalose Biosynthetic Complex in Pathogenic Fungi

Trehalose is a non-reducing disaccharide essential for pathogenic fungal survival and virulence. The biosynthesis of trehalose requires the trehalose-6-phosphate synthase, Tps1, and trehalose-6-phosphate phosphatase, Tps2. More importantly, the trehalose biosynthetic pathway is absent in mammals, conferring this pathway as an ideal target for antifungal drug design. However, lack of germane biochemical and structural information hinders antifungal drug design against these targets.

In this dissertation, macromolecular X-ray crystallography and biochemical assays were employed to understand the structures and functions of proteins involved in the trehalose biosynthetic pathway. I report here the first eukaryotic Tps1 structures from Candida albicans (C. albicans) and Aspergillus fumigatus (A. fumigatus) with substrates or substrate analogs. These structures reveal the key residues involved in substrate binding and catalysis. Subsequent enzymatic assays and cellular assays highlight the significance of these key Tps1 residues in enzyme function and fungal stress response. The Tps1 structure captured in its transition-state with a non-hydrolysable inhibitor demonstrates that Tps1 adopts an “internal return like” mechanism for catalysis. Furthermore, disruption of the trehalose biosynthetic complex formation through abolishing Tps1 dimerization reveals that complex formation has regulatory function in addition to trehalose production, providing additional targets for antifungal drug intervention.

I also present here the structure of the Tps2 N-terminal domain (Tps2NTD) from C. albicans, which may be involved in the proper formation of the trehalose biosynthetic complex. Deletion of the Tps2NTD results in a temperature sensitive phenotype. Further, I describe in this dissertation the structures of the Tps2 phosphatase domain (Tps2PD) from C. albicans, A. fumigatus and Cryptococcus neoformans (C. neoformans) in multiple conformational states. The structures of the C. albicans Tps2PD -BeF3-trehalose complex and C. neoformans Tps2PD(D24N)-T6P complex reveal extensive interactions between both glucose moieties of the trehalose involving all eight hydroxyl groups and multiple residues of both the cap and core domains of Tps2PD. These structures also reveal that steric hindrance is a key underlying factor for the exquisite substrate specificity of Tps2PD. In addition, the structures of Tps2PD in the open conformation provide direct visualization of the conformational changes of this domain that are effected by substrate binding and product release.

Last, I present the structure of the C. albicans trehalose synthase regulatory protein (Tps3) pseudo-phosphatase domain (Tps3PPD) structure. Tps3PPD adopts a haloacid dehydrogenase superfamily (HADSF) phosphatase fold with a core Rossmann-fold domain and a α/β fold cap domain. Despite lack of phosphatase activity, the cleft between the Tps3PPD core domain and cap domain presents a binding pocket for a yet uncharacterized ligand. Identification of this ligand could reveal the cellular function of Tps3 and any interconnection of the trehalose biosynthetic pathway with other cellular metabolic pathways.

Combined, these structures together with significant biochemical analyses advance our understanding of the proteins responsible for trehalose biosynthesis. These structures are ready to be exploited to rationally design or optimize inhibitors of the trehalose biosynthetic pathway enzymes. Hence, the work described in this thesis has laid the groundwork for the design of Tps1 and Tps2 specific inhibitors, which ultimately could lead to novel therapeutics to treat fungal infections.

Using electronic health record data for substance use Screening, Brief Intervention, and Referral to Treatment among adults with type 2 diabetes: Design of a National Drug Abuse Treatment Clinical Trials Network study.

BACKGROUND: The Affordable Care Act encourages healthcare systems to integrate behavioral and medical healthcare, as well as to employ electronic health records (EHRs) for health information exchange and quality improvement. Pragmatic research paradigms that employ EHRs in research are needed to produce clinical evidence in real-world medical settings for informing learning healthcare systems. Adults with comorbid diabetes and substance use disorders (SUDs) tend to use costly inpatient treatments; however, there is a lack of empirical data on implementing behavioral healthcare to reduce health risk in adults with high-risk diabetes. Given the complexity of high-risk patients' medical problems and the cost of conducting randomized trials, a feasibility project is warranted to guide practical study designs. METHODS: We describe the study design, which explores the feasibility of implementing substance use Screening, Brief Intervention, and Referral to Treatment (SBIRT) among adults with high-risk type 2 diabetes mellitus (T2DM) within a home-based primary care setting. Our study includes the development of an integrated EHR datamart to identify eligible patients and collect diabetes healthcare data, and the use of a geographic health information system to understand the social context in patients' communities. Analysis will examine recruitment, proportion of patients receiving brief intervention and/or referrals, substance use, SUD treatment use, diabetes outcomes, and retention. DISCUSSION: By capitalizing on an existing T2DM project that uses home-based primary care, our study results will provide timely clinical information to inform the designs and implementation of future SBIRT studies among adults with multiple medical conditions.

Design of Antibacterial Prochelators to Target Drug-Resistant Bacteria

Transition metals such as iron and copper are valued in biology for their redox activities because they are able to access various oxidation states. However, these transition metals are also implicated in a number of human disease states and play a role in bacterial infections. The ability to manipulate and monitor metal ions has vast implications on the fields of biology and human health. As such, the research described here covers two related goals: to manipulate metals in specific biological circumstances and to visualize this disturbance in cellular metal homeostasis.

Antibiotic resistance necessitates the development of drugs that exploit new mechanisms of action such as the disruption of metal homeostasis. In order to manipulate metals at the site of bacterial infection, two prochelators were developed around a β-lactam core such that the active chelator is released in the presence of bacteria that produce the resistance-causing β-lactamase enzyme. Both prochelators display enhanced activity toward resistant bacteria compared to clinical antibiotics.

Fluorescent sensors are a powerful tool for detecting small concentrations of biological analytes. Two analogs of a ratiometric fluorescent sensor were designed and synthesized to monitor cellular concentrations of copper and iron. These sensors were found to operate as designed in vitro; however the fluorescence intensity necessary for quantification of cellular metal pools has not yet been achieved.

Biologically Inspired Design of Protein-Silica Hybrid Nanoparticles for Drug Delivery Applications

The design and application of effective drug carriers is a fundamental concern in the delivery of therapeutics for the treatment of cancer and other vexing health problems. Traditionally utilized chemotherapeutics are limited in efficacy due to poor bioavailability as a result of their size and solubility as well as significant deleterious effects to healthy tissue through their inability to preferentially target pathological cells and tissues, especially in treatment of cancer. Thus, a major effort in the development of nanoscopic drug delivery vehicles for cancer treatment has focused on exploiting the inherent differences in tumor physiology and limiting the exposure of drugs to non-tumorous tissue, which is commonly achieved by encapsulation of chemotherapeutics within macromolecular or supramolecular carriers that incorporate targeting ligands and that enable controlled release. The overall aim of this work is to engineer a hybrid nanomaterial system comprised of protein and silica and to characterize its potential as an encapsulating drug carrier. The synthesis of silica, an attractive nanomaterial component because it is both biocompatible as well as structurally and chemically stable, within this system is catalyzed by self-assembled elastin-like polypeptide (ELP) micelles that incorporate of a class of biologically-inspired, silica-promoting peptides, silaffins. Furthermore, this methodology produces near-monodisperse, hybrid inorganic/micellar materials under mild reaction conditions such as temperature, pH and solvent. This work studies this material system along three avenues: 1) proof-of-concept silicification (i.e. the formation and deposition of silica upon organic materials) of ELP micellar templates, 2) encapsulation and pH-triggered release of small, hydrophobic chemotherapeutics, and 3) selective silicification of templates to potentiate retention of peptide targeting ability.

Mechanistic Models of Anti-HIV Microbicide Drug Delivery

A new modality for preventing HIV transmission is emerging in the form of topical microbicides. Some clinical trials have shown some promising results of these methods of protection while other trials have failed to show efficacy. Due to the relatively novel nature of microbicide drug transport, a rigorous, deterministic analysis of that transport can help improve the design of microbicide vehicles and understand results from clinical trials. This type of analysis can aid microbicide product design by helping understand and organize the determinants of drug transport and the potential efficacies of candidate microbicide products.

Microbicide drug transport is modeled as a diffusion process with convection and reaction effects in appropriate compartments. This is applied here to vaginal gels and rings and a rectal enema, all delivering the microbicide drug Tenofovir. Although the focus here is on Tenofovir, the methods established in this dissertation can readily be adapted to other drugs, given knowledge of their physical and chemical properties, such as the diffusion coefficient, partition coefficient, and reaction kinetics. Other dosage forms such as tablets and fiber meshes can also be modeled using the perspective and methods developed here.

The analyses here include convective details of intravaginal flows by both ambient fluid and spreading gels with different rheological properties and applied volumes. These are input to the overall conservation equations for drug mass transport in different compartments. The results are Tenofovir concentration distributions in time and space for a variety of microbicide products and conditions. The Tenofovir concentrations in the vaginal and rectal mucosal stroma are converted, via a coupled reaction equation, to concentrations of Tenofovir diphosphate, which is the active form of the drug that functions as a reverse transcriptase inhibitor against HIV. Key model outputs are related to concentrations measured in experimental pharmacokinetic (PK) studies, e.g. concentrations in biopsies and blood. A new measure of microbicide prophylactic functionality, the Percent Protected, is calculated. This is the time dependent volume of the entire stroma (and thus fraction of host cells therein) in which Tenofovir diphosphate concentrations equal or exceed a target prophylactic value, e.g. an EC50.

Results show the prophylactic potentials of the studied microbicide vehicles against HIV infections. Key design parameters for each are addressed in application of the models. For a vaginal gel, fast spreading at small volume is more effective than slower spreading at high volume. Vaginal rings are shown to be most effective if inserted and retained as close to the fornix as possible. Because of the long half-life of Tenofovir diphosphate, temporary removal of the vaginal ring (after achieving steady state) for up to 24h does not appreciably diminish Percent Protected. However, full steady state (for the entire stromal volume) is not achieved until several days after ring insertion. Delivery of Tenofovir to the rectal mucosa by an enema is dominated by surface area of coated mucosa and whether the interiors of rectal crypts are filled with the enema fluid. For the enema 100% Percent Protected is achieved much more rapidly than for vaginal products, primarily because of the much thinner epithelial layer of the mucosa. For example, 100% Percent Protected can be achieved with a one minute enema application, and 15 minute wait time.

Results of these models have good agreement with experimental pharmacokinetic data, in animals and clinical trials. They also improve upon traditional, empirical PK modeling, and this is illustrated here. Our deterministic approach can inform design of sampling in clinical trials by indicating time periods during which significant changes in drug concentrations occur in different compartments. More fundamentally, the work here helps delineate the determinants of microbicide drug delivery. This information can be the key to improved, rational design of microbicide products and their dosage regimens.