Search NASA⌕ Search

SEARCH · Search NASA

Results for “Docking”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76

Lucanthone and its derivative hycanthone inhibit apurinic endonuclease-1 (APE1) by direct protein binding

Lucanthone and hycanthone are thioxanthenone DNA intercalators used in the 1980s as antitumor agents. Lucanthone is in Phase I clinical trial, whereas hycanthone was pulled out of Phase II trials. Their potential mechanism of action includes DNA intercalation, inhibition of nucleic acid biosyntheses, and inhibition of enzymes like topoisomerases and the dual function base excision repair enzyme apurinic endonuclease 1 (APE1). Lucanthone inhibits the endonuclease activity of APE1, without affecting its redox activity. Our goal was to decipher the precise mechanism of APE1 inhibition as a prerequisite towards development of improved therapeutics that can counteract higher APE1 activity often seen in tumors. The IC{sub 50} values for inhibition of APE1 incision of depurinated plasmid DNA by lucanthone and hycanthone were 5 {mu}M and 80 nM, respectively. The K{sub D} values (affinity constants) for APE1, as determined by BIACORE binding studies, were 89 nM for lucanthone/10 nM for hycanthone. APE1 structures reveal a hydrophobic pocket where hydrophobic small molecules like thioxanthenones can bind, and our modeling studies confirmed such docking. Circular dichroism spectra uncovered change in the helical structure of APE1 in the presence of lucanthone/hycanthone, and notably, this effect was decreased (Phe266Ala or Phe266Cys or Trp280Leu) or abolished (Phe266Ala/Trp280Ala) when hydrophobic site mutants were employed. Reduced inhibition by lucanthone of the diminished endonuclease activity of hydrophobic mutant proteins (as compared to wild type APE1) supports that binding of lucanthone to the hydrophobic pocket dictates APE1 inhibition. The DNA binding capacity of APE1 was marginally inhibited by lucanthone, and not at all by hycanthone, supporting our hypothesis that thioxanthenones inhibit APE1, predominantly, by direct interaction. Finally, lucanthone-induced degradation was drastically reduced in the presence of short and long lived free radical scavengers, e.g., TRIS and DMSO, suggesting that the mechanism of APE1 breakdown may involve free radical-induced peptide bond cleavage.

59 BASIC BIOLOGICAL SCIENCES↗

Manipulating Fatty Acid Biosynthesis in Microalgae for Biofuel through Protein-Protein Interactions

Microalgae are a promising feedstock for renewable fuels, and algal metabolic engineering can lead to crop improvement, thus accelerating the development of commercially viable biodiesel production from algae biomass. We demonstrate that protein-protein interactions between the fatty acid acyl carrier protein (ACP) and thioesterase (TE) govern fatty acid hydrolysis within the algal chloroplast. Using green microalga Chlamydomonas reinhardtii (Cr) as a model, a structural simulation of docking CrACP to CrTE identifies a protein-protein recognition surface between the two domains. A virtual screen reveals plant TEs with similar in silico binding to CrACP. Employing an activity-based crosslinking probe designed to selectively trap transient protein-protein interactions between the TE and ACP, we demonstrate in vitro that CrTE must functionally interact with CrACP to release fatty acids, while TEs of vascular plants show no mechanistic crosslinking to CrACP. This is recapitulated in vivo, where overproduction of the endogenous CrTE increased levels of short-chain fatty acids and engineering plant TEs into the C. reinhardtii chloroplast did not alter the fatty acid profile. These findings highlight the critical role of protein-protein interactions in manipulating fatty acid biosynthesis for algae biofuel engineering as illuminated by activity-based probes.

59 BASIC BIOLOGICAL SCIENCES↗

Crystal Structure of BamB from Pseudomonas aeruginosa and Functional Evaluation of Its Conserved Structural Features

The assembly of β-barrel Outer Membrane Proteins (OMPs) in the outer membrane is essential for Gram-negative bacteria. The process requires the β-Barrel Assembly Machine (BAM), a multiprotein complex that, in E. coli, is composed of the OMP BamA and four lipoproteins BamB-E. Whereas BamA and BamD are essential, deletion of BamB, C or E produce membrane permeability defects. Here we present the high-resolution structure of BamB from Pseudomonas aeruginosa. This protein can complement the deletion of bamB in E. coli indicating that they are functionally equivalent. Conserved structural features include an eight-bladed β-propeller fold stabilized by tryptophan docking motifs with a central pore about 8 Å in diameter at the narrowest point. This pore distinguishes BamB from related b-propellers, such as quinoprotein dehydrogenases. However, a double mutation designed to block this pore was fully functional indicating that the opening is not essential. Two loops protruding from the bottom of the propeller are conserved and mediate binding to BamA. Conversely, an additional loop only present in E. coli BamB is not required for function. A cluster of highly conserved residues in a groove between blades 6 and 7 is crucial for proper BamB folding or biogenesis. It has been proposed that BamB may bind nascent OMPs by β-augmentation to its propeller outer strands, or recognize the aromatic residue signature at the C-terminus of OMPs. However, Isothermal Titration Calorimetry experiments and structural analysis do not support these proposals. The structural and mutagenesis analysis suggests that the main function of BamB is to bind and modulate BamA, rather than directly interact with nascent OMPs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-Ray Crystal Structure of the Ancestral 3-Ketosteroid Receptor-Progesterone-Mifepristone Complex Shows Mifepristone Bound at the Coactivator Binding Interface

Steroid receptors are a subfamily of nuclear receptors found throughout all metazoans. They are highly important in the regulation of development, inflammation, and reproduction and their misregulation has been implicated in hormone insensitivity syndromes and cancer. Steroid binding to SRs drives a conformational change in the ligand binding domain that promotes nuclear localization and subsequent interaction with coregulator proteins to affect gene regulation. SRs are important pharmaceutical targets, yet most SR-targeting drugs have off-target pharmacology leading to unwanted side effects. A better understanding of the structural mechanisms dictating ligand specificity and the evolution of the forces that created the SR-hormone pairs will enable the design of better pharmaceutical ligands. In order to investigate this relationship, we attempted to crystallize the ancestral 3-ketosteroid receptor (ancSR2) with mifepristone, a SR antagonist. Here, we present the x-ray crystal structure of the ancestral 3-keto steroid receptor (ancSR2)-progesterone complex at a resolution of 2.05 Å. This improves upon our previously reported structure of the ancSR2-progesterone complex, permitting unambiguous assignment of the ligand conformation within the binding pocket. Surprisingly, we find mifepristone, fortuitously docked at the protein surface, poised to interfere with coregulator binding. Recent attention has been given to generating pharmaceuticals that block the coregulator binding site in order to obstruct coregulator binding and achieve tissue-specific SR regulation independent of hormone binding. Mifepristone’s interaction with the coactivator cleft of this SR suggests that it may be a useful molecular scaffold for further coactivator binding inhibitor development.

Colucci, Jennifer K.↗

Understanding a Substrate’s Product Regioselectivity in a Family of Enzymes: A Case Study of Acetaminophen Binding in Cytochrome P450s

Product regioselectivity as influenced by molecular recognition is a key aspect of enzyme catalysis. We applied large-scale two-dimensional (2D) umbrella sampling (USP) simulations to characterize acetaminophen (APAP) binding in the active sites of the family of Cytochrome P450 (CYP) enzymes as a case study to show the different regioselectivity exhibited by a single substrate in comparative enzymes. Our results successfully explain the experimentally observed product regioselectivity for all five human CYPs included in this study, demonstrating that binding events play an important role in determining regioselectivity. In CYP2C9 and CYP3A4, weak interactions in an overall large active site cavity result in a fairly small binding free energy difference between APAP reactive binding states, consistent with experimental results that show little preference for resulting metabolites. In contrast, in CYP1A2 and CYP2E1, APAP is strongly restrained by a compact binding pocket, leading to a preferred binding conformation. The calculated binding equilibrium of APAP within the compact active site of CYP2A6 is able to predict the experimentally documented product ratios and is also applied to explain APAP regioselectivity in CYP1A2 and CYP2C9. APAP regioselectivity seems to be related to the selectivity for one binding conformation over another binding conformation as dictated by the size and shape of the active site. Additionally, unlike docking and molecular dynamics (MD), our free energy calculations successfully reproduced a unique APAP pose in CYP3A4 that had been reported experimentally, suggesting this approach is well suited to find the realistic binding pose and the lowest-energy starting structure for studying the chemical reaction step in the future.

59 BASIC BIOLOGICAL SCIENCES↗

Identification of a Novel Drug Lead That Inhibits HCV Infection and Cell-to-Cell Transmission by Targeting the HCV E2 Glycoprotein

We report that Hepatitis C Virus (HCV) infects 200 million individuals worldwide. Although several FDA approved drugs targeting the HCV serine protease and polymerase have shown promising results, there is a need for better drugs that are effective in treating a broader range of HCV genotypes and subtypes without being used in combination with interferon and/or ribavirin. Recently, two crystal structures of the core of the HCV E2 protein (E2c) have been determined, providing structural information that can now be used to target the E2 protein and develop drugs that disrupt the early stages of HCV infection by blocking E2’s interaction with different host factors. Using the E2c structure as a template, we have created a structural model of the E2 protein core (residues 421–645) that contains the three amino acid segments that are not present in either structure. Computational docking of a diverse library of 1,715 small molecules to this model led to the identification of a set of 34 ligands predicted to bind near conserved amino acid residues involved in the HCV E2: CD81 interaction. We used surface plasmon resonance detection to screen the ligand set for binding to recombinant E2 protein, and the best binders were subsequently tested to identify compounds that inhibit the infection of Huh-7 cells by HCV. One compound, 281816, blocked E2 binding to CD81 and inhibited HCV infection in a genotype-independent manner with IC50’s ranging from 2.2 µM to 4.6 µM. 281816 blocked the early and late steps of cell-free HCV entry and also abrogated the cell-to-cell transmission of HCV. Collectively the results obtained with this new structural model of E2c suggest the development of small molecule inhibitors such as 281816 that target E2 and disrupt its interaction with CD81 may provide a new paradigm for HCV treatment.

60 APPLIED LIFE SCIENCES↗

The interaction of RNA helicase DDX3 with HIV-1 Rev-CRM1-RanGTP complex during the HIV replication cycle

Molecular traffic between the nucleus and the cytoplasm is regulated by the nuclear pore complex (NPC), which acts as a highly selective channel perforating the nuclear envelope in eukaryotic cells. The human immunodeficiency virus (HIV) exploits the nucleocytoplasmic pathway to export its RNA transcripts across the NPC to the cytoplasm. Despite extensive study on the HIV life cycle and the many drugs developed to target this cycle, no current drugs have been successful in targeting the critical process of viral nuclear export, even though HIV’s reliance on a single host protein, CRM1, to export its unspliced and partially spliced RNA transcripts makes it a tempting target. Due to recent findings implicating a DEAD-box helicase, DDX3, in HIV replication and a member of the export complex, it has become an appealing target for anti-HIV drug inhibition. In the present research, we have applied a hybrid computational protocol to analyze protein-protein interactions in the HIV mRNA export cycle. This method is based on molecular docking followed by molecular dynamics simulation and accompanied by approximate free energy calculation (MM/GBSA), computational alanine scanning, clustering, and evolutionary analysis. We highlight here some of the most likely binding modes and interfacial residues between DDX3 and CRM1 both in the absence and presence of RanGTP. This work shows that although DDX3 can bind to free CRM1, addition of RanGTP leads to more concentrated distribution of binding modes and stronger binding between CRM1 and RanGTP.

59 BASIC BIOLOGICAL SCIENCES↗

The quinic acid derivative KZ-41 prevents glucose-induced caspase-3 activation in retinal endothelial cells through an IGF-1 receptor dependent mechanism

Retinal microaneurysms, an early disease manifestation of diabetic retinopathy, are associated with retinal endothelial cell (REC) death and macular edema. We previously demonstrated that a quinic acid (QA) analog, KZ-41, promoted REC survival by blunting stress-induced p38 MAPK activation. Herein, we sought to expand our understanding of the prosurvival signal transduction pathways actuated by KZ-41. Using human RECs exposed to high glucose (25 mM, 72 hours), we demonstrated that KZ-41 blocks caspase-3 activation by triggering phosphorylation of the PI3K regulatory subunit (p85; Tyr458) and its downstream target Akt (Ser473). Akt signal transduction was accompanied by autophosphorylation of the receptor tyrosine kinase, insulin growth factor-1 receptor (IGF-1R). IGF-1R knockdown using either the tyrosine kinase inhibitor AG1024 or silencing RNA abolished KZ-41’s pro-survival effect. Under high glucose stress, caspase-3 activation correlated with elevated ERK1/2 phosphorylation and decreased insulin receptor substrate-1 (IRS-1) levels. KZ-41 decreased ERK1/2 phosphorylation and reversed the glucose-dependent reduction in IRS-1. To gain insight into the mechanistic basis for IGF-1R activation by KZ-41, we used molecular modeling and docking simulations to explore a possible protein:ligand interaction between the IGF-1R kinase domain and KZ-41. Computational investigations suggest two possible KZ-41 binding sites within the kinase domain: a region with high homology to the insulin receptor contains one potential allosteric binding site, and another potential site on the other side of the kinase domain, near the hinge domain. These data, together with previous proof-of-concept efficacy studies demonstrating KZ-41 mitigates pathologic retinal neovascularization in the murine oxygen-induced retinopathy model, suggests that QA derivatives may offer therapeutic benefit in ischemic retinopathies.

59 BASIC BIOLOGICAL SCIENCES↗

Structure and Recognition of a Novel HIV-1 gp120-gp41 Interface Antibody that Caused MPER Exposure through Viral Escape

A comprehensive understanding of the regions on HIV-1 envelope trimers targeted by broadly neutralizing antibodies may contribute to rational design of an HIV-1 vaccine. We previously identified a participant in the CAPRISA cohort, CAP248, who developed trimer-specific antibodies capable of neutralizing 60% of heterologous viruses at three years post-infection. Here, we report the isolation by B cell culture of monoclonal antibody CAP248-2B, which targets a novel membrane proximal epitope including elements of gp120 and gp41. Despite low maximum inhibition plateaus, often below 50% inhibitory concentrations, the breadth of CAP248-2B significantly correlated with donor plasma. Site-directed mutagenesis, X-ray crystallography, and negative-stain electron microscopy 3D reconstructions revealed how CAP248-2B recognizes a cleavage-dependent epitope that includes the gp120 C terminus. While this epitope is distinct, it overlapped in parts of gp41 with the epitopes of broadly neutralizing antibodies PGT151, VRC34, 35O22, 3BC315, and 10E8. CAP248-2B has a conformationally variable paratope with an unusually long 19 amino acid light chain third complementarity determining region. Two phenylalanines at the loop apex were predicted by docking and mutagenesis data to interact with the viral membrane. Neutralization by CAP248-2B is not dependent on any single glycan proximal to its epitope, and low neutralization plateaus could not be completely explained by N- or O-linked glycosylation pathway inhibitors, furin co-transfection, or pre-incubation with soluble CD4. Viral escape from CAP248-2B involved a cluster of rare mutations in the gp120-gp41 cleavage sites. Simultaneous introduction of these mutations into heterologous viruses abrogated neutralization by CAP248-2B, but enhanced neutralization sensitivity to 35O22, 4E10, and 10E8 by 10-100-fold. Altogether, this study expands the region of the HIV-1 gp120-gp41 quaternary interface that is a target for broadly neutralizing antibodies and identifies a set of mutations in the gp120 C terminus that exposes the membrane-proximal external region of gp41, with potential utility in HIV vaccine design.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular basis of P[II] major human rotavirus VP8* domain recognition of histo-blood group antigens

Initial cell attachment of rotavirus (RV) to specific cell surface glycan receptors, which is the essential first step in RV infection, is mediated by the VP8* domain of the spike protein VP4. Recently, human histo-blood group antigens (HBGAs) have been identified as receptors or attachment factors for human RV strains. RV strains in the P[4] and P[8] genotypes of the P[II] genogroup share common recognition of the Lewis b (Le b ) and H type 1 antigens, however, the molecular basis of receptor recognition by the major human P[8] RVs remains unknown due to lack of experimental structural information. Here, we used nuclear magnetic resonance (NMR) spectroscopy-based titration experiments and NMR-derived high ambiguity driven docking (HADDOCK) methods to elucidate the molecular basis for P[8] VP8* recognition of the Leb (LNDFH I) and type 1 HBGAs. We also used X-ray crystallography to determine the molecular details underlying P[6] recognition of H type 1 HBGAs. Unlike P[6]/P[19] VP8*s that recognize H type 1 HBGAs in a binding surface composed of an α-helix and a β-sheet, referred as the “βα binding site”, the P[8] and P[4] VP8*s bind Le b HBGAs in a previously undescribed pocket formed by the edges of two β-sheets, referred to as the “ββ binding site”. Importantly, the P[8] and P[4] VP8*s retain binding capability to non-Le b type 1 HBGAs using the βα binding site. The presence of two distinct binding sites for Le b and non-Le b HBGA glycans in the P[8] and P[4] VP8* domains suggests host-pathogen co-evolution under structural and functional adaptation of RV pathogens to host glycan polymorphisms. Assessment and understanding of the precise impact of this co-evolutionary process in determining RV host ranges and cross-species RV transmission should facilitate improved RV vaccine development and prediction of future RV strain emergence and epidemics.

59 BASIC BIOLOGICAL SCIENCES↗

Viral Nuclease Inhibitors: Small molecule disruptors of the UL12 alkaline nuclease display broad anti-herpes virus activity

Herpes simplex virus 1 (HSV-1) UL12 encodes a highly conserved 5′ → 3′ alkaline exonuclease that is essential for the production of infectious virus. Together with the viral single-stranded DNA-binding/annealing protein ICP8, UL12 functions as a two-component recombinase that mediates recombination-dependent viral DNA replication. Here, we present the crystal structure of the catalytic domain of the HSV alkaline nuclease (UL12), which provides the first view of an α-herpesvirus alkaline nuclease. Using this structure, we optimized a series of small-molecule viral nuclease inhibitors (VNIs) that target the UL12 active site and potently inhibit UL12 exonuclease activity in vitro. We have thus established a robust platform for structure-based docking, SAR analysis and rational inhibitor design. Because UL12 orthologs are conserved across all human herpesviruses, we examined the activity of these compounds against the β- and γ-herpesvirus alkaline nucleases UL98 and SOX and found that they inhibit all three enzymes. The VNIs also exhibit antiviral activity against HSV-1 and HCMV in cell culture. EC 50 and IC 50 values were in the nanomolar to low micromolar range. Together, these findings establish herpesvirus alkaline nucleases as conserved, druggable antiviral targets and provide a foundation for the development of broad-spectrum anti-herpesvirus therapeutics, either as standalone agents or in combination with existing nucleoside analogs.

Sharma, Nidhi↗

Structural basis of divergent cyclin-dependent kinase activation by Spy1/RINGO proteins

Cyclin-dependent kinases (Cdks) are principal drivers of cell division and are an important therapeutic target to inhibit aberrant proliferation. Cdk enzymatic activity is tightly controlled through cyclin interactions, posttranslational modifications, and binding of inhibitors such as the p27 tumor suppressor protein. Spy1/RINGO (Spy1) proteins bind and activate Cdk but are resistant to canonical regulatory mechanisms that establish cell-cycle checkpoints. Cancer cells exploit Spy1 to stimulate proliferation through inappropriate activation of Cdks, yet the mechanism is unknown. We have determined crystal structures of the Cdk2-Spy1 and p27-Cdk2-Spy1 complexes that reveal how Spy1 activates Cdk. We find that Spy1 confers structural changes to Cdk2 that obviate the requirement of Cdk activation loop phosphorylation. Spy1 lacks the cyclin-binding site that mediates p27 and substrate affinity, explaining why Cdk-Spy1 is poorly inhibited by p27 and lacks specificity for substrates with cyclin-docking sites. We identify mutations in Spy1 that ablate its ability to activate Cdk2 and to proliferate cells. Our structural description of Spy1 provides important mechanistic insights that may be utilized for targeting upregulated Spy1 in cancer.

59 BASIC BIOLOGICAL SCIENCES↗

Structural basis for the specificity of PPM1H phosphatase for Rab GTPases

LRRK2 serine/threonine kinase is associated with inherited Parkinson’s disease. LRRK2 phosphorylates a subset of Rab GTPases within their switch 2 motif to control their interactions with effectors. Recent work has shown that the metal-dependent protein phosphatase PPM1H counteracts LRRK2 by dephosphorylating Rabs. PPM1H is highly selective for LRRK2 phosphorylated Rabs, and closely related PPM1J exhibits no activity towards substrates such as Rab8a phosphorylated at Thr72 (pThr72). Here, we have identified the molecular determinant of PPM1H specificity for Rabs. The crystal structure of PPM1H reveals a structurally conserved phosphatase fold that strikingly has evolved a 110-residue flap domain adjacent to the active site. The flap domain distantly resembles tudor domains that interact with histones in the context of epigenetics. Cellular assays, crosslinking and 3-D modelling suggest that the flap domain encodes the docking motif for phosphorylated Rabs. Consistent with this hypothesis, a PPM1J chimaera with the PPM1H flap domain dephosphorylates pThr72 of Rab8a both in vitro and in cellular assays. Therefore, PPM1H has acquired a Rab-specific interaction domain within a conserved phosphatase fold.

60 APPLIED LIFE SCIENCES↗

AmeriFlux US-Act Acton Lake Flux Tower Site

This is the AmeriFlux version of the carbon flux data for the site US-Act Acton Lake Flux Tower Site. Site Description - This tower is located at Acton Lake, a eutrophic reservoir immediately surrounded by forest but draining an agricultural watershed in Hueston Woods State Park outside of Oxford, OH. The EC instrumentation was originally deployed 25 m from shore on a dock piling. In April of 2018 a lake tower was installed >200 m from the closest shoreline. The EC instrumentation was deployed on this aquatic tower on 7 May 2018.

Waldo, Sarah↗

U.S. Freight Transload Facilities Dataset

The U.S. Freight Transload Facilities Dataset provides location information (latitude, longitude, zip, city, county, state)for more than 9,000 facilities across 50 U.S. States where freight may be transferred between waterways, railways, and roadways. The dataset lists the known modes and available direction(s) for freight transfers at each facility as of 2024. The U.S. Freight Transload Facilities dataset was built by mining and fusing several public sources, such as the USACE Master Docks Plus, the USDOT National Transportation Atlas Database (NTAD), files from the Intermodal Association of North America (IANA), and the industry publication Bulk Transloader. The dataset constitutes a key piece of a multimodal freight transportation network and routing algorithm developed by USACE-ERDC. The dataset is shared as a .csv file. The dataset is published for research purposes and should not be considered exhaustive or authoritative.

Peterson, Steven [ORNL] (ORCID:0000000287672998)↗

Renewable Energy Opportunities at the Kanto Installations, Japan

This document provides an overview of renewable resource development potential at the U.S. Army installations in the Kanto region in Japan, which includes Camp Zama, Yokohama North Dock, Sagamihara Family Housing Area (SFHA), Sagami General Depot, and Akasaka Press Center. This effort focuses on grid-connected generation of electricity from renewable energy sources and also on ground source heat pumps for heating and cooling buildings. The effort was funded by the Huntsville Army Corps of Engineers, and includes the development of a methodology for renewable resource assessment at Army installations located on foreign soil. The methodology is documented in Renewable Energy Assessment Methodology for Japanese OCONUS Army Installations. The site visit to the Kanto installations took place on April 5 and 6, 2010. At the current time, there are some renewable technologies that show economic potential. Because of siting restrictions and the small size of these installations, development of most renewable energy technologies will likely be limited to Camp Zama. Project feasibility is based on installation-specific resource availability and energy costs and projections based on accepted life-cycle cost methods. Development of any renewable energy project will be challenging, as it will require investigation into existing contractual obligations, new contracts that could be developed, the legality of certain partnerships, and available financing avenues, which involves the U.S. Forces Japan (USFJ), the Government of Japan (GOJ), and a number of other parties on both sides. The Army will not be able to implement a project without involvement and approval from the other services and multiple levels of Japanese government. However, implementation of renewable energy projects could be an attractive method for GOJ to reduce greenhouse gas emissions and lower annual utility payments to USFJ. This report recommends projects to pursue and offers approaches to use. The most promising opportunities include waste-to-energy and ground source heat pumps. Solar photovoltaics (PV) may also prove successful. Other resources were found to be insufficient on the Kanto installations.

15 GEOTHERMAL ENERGY↗

Quantum mechanical approaches to in silico enzyme characterization and drug design

The astonishing, exponentially increasing rates of genome sequencing has led to one of the most significant challenges for the biological and computational sciences in the 21 st century: assigning the likely functions of the encoded proteins. Enzymes represent a particular challenge, and a critical one, because the universe of enzymes is likely to contain many novel functions that may be useful for synthetic biology, or as drug targets. Current approaches to protein annotation are largely based on bioinformatics. At the simplest level, this annotation involves transferring the annotations of characterized enzymes to related sequences. In practice, however, there is no simple, sequence based criterion for transferring annotations, and bioinformatics alone cannot propose new enzymatic functions. Structure-based computational methods have the potential to address these limitations, by identifying potential substrates of enzymes, as we and others have shown. One successful approach has used in silico 'docking' methods, more commonly applied in structure-based drug design, to identify possible metabolite substrates. A major limitation of this approach is that it only considers substrate binding, and does not directly assess the potential of the enzyme to catalyze a particular reaction using a particular substrate. That is, substrate binding affinity is necessary but not sufficient to assign function. A reaction profile is ultimately what is needed for a more complete quantitative description of function. To address this rather fundamental limitation, they propose to use quantum mechanical methods to explicitly compute transition state barriers that govern the rates of catalysis. Although quantum mechanical, and mixed quantum/classical (QM/MM), methods have been used extensively to investigate enzymatic reactions, the focus has been primarily on elucidating complex reaction mechanisms. Here, the key catalytic steps are known, and they use these methods quantify substrate specificity. That is, we bring the power of quantum mechanics to bear on the problem of annotating enzyme function, which is a novel approach. Although it has been clear to us at the Jacobson group for some time that enzyme specificity may be encoded in transition states, rather than simply substrate recognition, the main limitation has always been computational expense. Using a hierarchy of different methods, they can reduce the list of plausible substrates of an enzyme to a small number in most cases, but even identifying the transition states for a dozen plausible substrates requires significant computational effort, beyond what is practical using standard QM/MM methods. For this project, they have chosen two enzyme superfamilies which they have used as 'model systems' for functional assignment. The enolase superfamily is a large group of {alpha}-{beta} barrel enzymes with highly diverse substrates and chemical transformations. Despite decades of work, over a third of the superfamily remains unassigned, which means that the remaining cases are by definition difficult to assign. They have focused on acid sugar dehydratases, and have considerable expertise on the matter. They are also interested in the isoprenoid synthase superfamily, which is of central interest to the synthetic biology community, because these enzymes are used by nature to create complex rare natural products of medicinal value. the most notable example of this is the artemisinin, an antimalarial compound that is found in trace amounts in the wormwod root. From the standpoint of enzyme function assignment, these enzymes are intriguing because they use a small number of chemically simple substrates to generate, potentially, tens of thousands of different products. Hence, substrate binding specificity is only a small part of the challenge; the key is determining how the enzyme directs the carbocation chemistry to specific products. These more complex modeling approaches clearly require quantum mechanical methods.

59 BASIC BIOLOGICAL SCIENCES↗