Atomic-Level Characterization of the Chain-Flipping Mechanism in Fatty-Acids Biosynthesis

Exploring Aspartic Protease Inhibitor Binding to Design-Selective Antimalarials

Selectivity is a major issue in the development of drugs targeting pathogen aspartic proteases. Here, we explore the selectivity-determining factors by studying specifically designed malaria aspartic protease (plasmepsin) open-flap inhibitors. Metadynamics simulations are used to uncover the complex binding/unbinding pathways of these inhibitors and describe the critical transition states in atomistic resolution. The simulation results are compared with experimentally determined enzymatic activities. Our findings demonstrate that plasmepsin inhibitor selectivity can be achieved by targeting the flap loop with hydrophobic substituents that enable ligand binding under the flap loop, as such a behavior is not observed for several other aspartic proteases. The ability to estimate the selectivity of compounds before they are synthesized is of considerable importance in drug design; therefore, we expect that our approach will be useful in selective inhibitor designs against not only aspartic proteases but also other enzyme classes.

Path to Actinorhodin: Regio- and Stereoselective Ketone Reduction by a Type II Polyketide Ketoreductase Revealed in Atomistic Detail

In type II polyketide synthases (PKSs), which typically biosynthesize several antibiotic and antitumor compounds, the substrate is a growing polyketide chain, shuttled between individual PKS enzymes, while covalently tethered to an acyl carrier protein (ACP): this requires the ACP interacting with a series of different enzymes in succession. During biosynthesis of the antibiotic actinorhodin, produced by Streptomyces coelicolor, one such key binding event is between an ACP carrying a 16-carbon octaketide chain (actACP) and a ketoreductase (actKR). Once the octaketide is bound inside actKR, it is likely cyclized between C7 and C12 and regioselective reduction of the ketone at C9 occurs: how these elegant chemical and conformational changes are controlled is not yet known. Here, we perform protein-protein docking, protein NMR, and extensive molecular dynamics simulations to reveal a probable mode of association between actACP and actKR; we obtain and analyze a detailed model of the C7-C12-cyclized octaketide within the actKR active site; and we confirm this model through multiscale (QM/MM) reaction simulations of the key ketoreduction step. Molecular dynamics simulations show that the most thermodynamically stable cyclized octaketide isomer (7R,12R) also gives rise to the most reaction competent conformations for ketoreduction. Subsequent reaction simulations show that ketoreduction is stereoselective as well as regioselective, resulting in an S-alcohol. Our simulations further indicate several conserved residues that may be involved in selectivity of C7-12 cyclization and C9 ketoreduction. Detailed insights obtained on ACP-based substrate presentation in type II PKSs can help design ACP-ketoreductase systems with altered regio- or stereoselectivity.

Control of Unsaturation in De Novo Fatty Acid Biosynthesis by FabA

Carrier protein-dependent biosynthesis provides a thiotemplated format for the production of natural products. Within these pathways, many reactions display exquisite substrate selectivity, a regulatory framework proposed to be controlled by protein-protein interactions (PPIs). In Escherichia coli, unsaturated fatty acids are generated within the de novo fatty acid synthase by a chain length-specific interaction between the acyl carrier protein AcpP and the isomerizing dehydratase FabA. To evaluate PPI-based control of reactivity, interactions of FabA with AcpP bearing multiple sequestered substrates were analyzed through NMR titration and guided high-resolution docking. Through a combination of quantitative binding constants, residue-specific perturbation analysis, and high-resolution docking, a model for substrate control via PPIs has been developed. The in silico results illuminate the mechanism of FabA substrate selectivity and provide a structural rationale with atomic detail. Helix III positioning in AcpP communicates sequestered chain length identity recognized by FabA, demonstrating a powerful strategy to regulate activity by allosteric control. These studies broadly illuminate carrier protein-dependent pathways and offer an important consideration for future inhibitor design and pathway engineering.

Probing allosteric regulations with coevolution-driven molecular simulations

Protein-mediated allosteric regulations are essential in biology, but their quantitative characterization continues to posit formidable challenges for both experiments and computations. Here, we combine coevolutionary information, multiscale molecular simulations, and free-energy methods to interrogate and quantify the allosteric regulation of functional changes in protein complexes. We apply this approach to investigate the regulation of adenylyl cyclase (AC) by stimulatory and inhibitory G proteins—a prototypical allosteric system that has long escaped from in-depth molecular characterization. We reveal a surprisingly simple ON/OFF regulation of AC functional dynamics through multiple pathways of information transfer. The binding of G proteins reshapes the free-energy landscape of AC following the classical population-shift paradigm. The model agrees with structural and biochemical data and reveals previously unknown experimentally consistent intermediates. Our approach showcases a general strategy to explore uncharted functional space in complex biomolecular regulations.

Probing the transport of Ni(II) ions through the internal tunnels of the Helicobacter pylori UreDFG multimeric protein complex

The survival of several pathogenic bacteria, such as Helicobacter pylori (Hp), relies on the activity of the nickel-dependent enzyme urease. Nickel insertion into urease is mediated by a multimeric chaperone complex (HpUreDFG) that is responsible for the transport of Ni(II) from a conserved metal binding motif located in the UreG dimer (CPH motif) to the catalytic site of the enzyme. The X-ray structure of HpUreDFG revealed the presence of water-filled tunnels that were proposed as a route for Ni(II) translocation. Here, we probe the transport of Ni(II) through the internal tunnels of HpUreDFG, from the CPH motif to the external surface of the complex, using microsecond-long enhanced molecular dynamics simulations. The results suggest a "bucket-brigade" mechanism whereby Ni(II) can be transported through a series of stations found along these internal pathways.

Decoding allosteric regulation by the acyl carrier protein

Enzymes in multistep metabolic pathways utilize an array of regulatory mechanisms to maintain a delicate homeostasis [K. Magnuson, S. Jackowski, C. O. Rock, J. E. Cronan, Jr, Microbiol. Rev. 57, 522-542 (1993)]. Carrier proteins in particular play an essential role in shuttling substrates between appropriate enzymes in metabolic pathways. Although hypothesized [E. Płoskoń et al., Chem. Biol. 17, 776-785 (2010)], allosteric regulation of substrate delivery has never before been demonstrated for any acyl carrier protein (ACP)-dependent pathway. Studying these mechanisms has remained challenging due to the transient and dynamic nature of protein-protein interactions, the vast diversity of substrates, and substrate instability [K. Finzel, D. J. Lee, M. D. Burkart, ChemBioChem 16, 528-547 (2015)]. Here we demonstrate a unique communication mechanism between the ACP and partner enzymes using solution NMR spectroscopy and molecular dynamics to elucidate allostery that is dependent on fatty acid chain length. We demonstrate that partner enzymes can allosterically distinguish between chain lengths via protein-protein interactions as structural features of substrate sequestration are translated from within the ACP four-helical bundle to the protein surface, without the need for stochastic chain flipping. These results illuminate details of cargo communication by the ACP that can serve as a foundation for engineering carrier protein-dependent pathways for specific, desired products.

Elucidation of transient protein-protein interactions within carrier protein-dependent biosynthesis

Fatty acid biosynthesis (FAB) is an essential and highly conserved metabolic pathway. In bacteria, this process is mediated by an elaborate network of protein•protein interactions (PPIs) involving a small, dynamic acyl carrier protein that interacts with dozens of other partner proteins (PPs). These PPIs have remained poorly characterized due to their dynamic and transient nature. Using a combination of solution-phase NMR spectroscopy and protein-protein docking simulations, we report a comprehensive residue-by-residue comparison of the PPIs formed during FAB in Escherichia coli. This technique describes and compares the molecular basis of six discrete binding events responsible for E. coli FAB and offers insights into a method to characterize these events and those in related carrier protein-dependent pathways.

Interfacial plasticity facilitates high reaction rate of E. coli FAS malonyl-CoA:ACP transacylase, FabD

Fatty acid synthases (FASs) and polyketide synthases (PKSs) iteratively elongate and often reduce two-carbon ketide units in de novo fatty acid and polyketide biosynthesis. Cycles of chain extensions in FAS and PKS are initiated by an acyltransferase (AT), which loads monomer units onto acyl carrier proteins (ACPs), small, flexible proteins that shuttle covalently linked intermediates between catalytic partners. Formation of productive ACP-AT interactions is required for catalysis and specificity within primary and secondary FAS and PKS pathways. Here, we use the Escherichia coli FAS AT, FabD, and its cognate ACP, AcpP, to interrogate type II FAS ACP-AT interactions. We utilize a covalent crosslinking probe to trap transient interactions between AcpP and FabD to elucidate the X-ray crystal structure of a type II ACP-AT complex. Our structural data are supported using a combination of mutational, crosslinking, and kinetic analyses, and long-timescale molecular dynamics (MD) simulations. Together, these complementary approaches reveal key catalytic features of FAS ACP-AT interactions. These mechanistic inferences suggest that AcpP adopts multiple, productive conformations at the AT binding interface, allowing the complex to sustain high transacylation rates. Furthermore, MD simulations support rigid body subdomain motions within the FabD structure that may play a key role in AT activity and substrate selectivity.

Activity Mapping the Acyl Carrier Protein: Elongating Ketosynthase Interaction in Fatty Acid Biosynthesis

Elongating ketosynthases (KSs) catalyze carbon-carbon bond-forming reactions during the committed step for each round of chain extension in both fatty acid synthases (FASs) and polyketide synthases (PKSs). A small α-helical acyl carrier protein (ACP) shuttles fatty acyl intermediates between enzyme active sites. To accomplish this task, the ACP relies on a series of dynamic interactions with multiple partner enzymes of FAS and associated FAS-dependent pathways. Recent structures of the Escherichia coli FAS ACP, AcpP, in covalent complexes with its two cognate elongating KSs, FabF and FabB, provide high-resolution details of these interfaces, but a systematic analysis of specific interfacial interactions responsible for stabilizing these complexes has not yet been undertaken. Here, we use site-directed mutagenesis with both in vitro and in vivo activity analyses to quantitatively evaluate these contacting surfaces between AcpP and FabF. We delineate the FabF interface into three interacting regions and demonstrate the effects of point mutants, double mutants, and region deletion variants. Results from these analyses reveal a robust and modular FabF interface capable of tolerating seemingly critical interface mutations with only the deletion of an entire region significantly compromising activity. Structure and sequence analyses of FabF orthologs from related type II FAS pathways indicate significant conservation of type II FAS KS interface residues and, overall, support its delineation into interaction regions. These findings strengthen our mechanistic understanding of molecular recognition events between ACPs and FAS enzymes and provide a blueprint for engineering ACP-dependent biosynthetic pathways.

Exploring Conformational Dynamics of the Extracellular Venus flytrap Domain of the GABAB Receptor: A Path-Metadynamics Study

γ-Aminobutyric acid (GABA) is the main inhibitory neurotransmitter in the central nervous system (CNS). Dysfunctional GABAergic neurotransmission is associated with numerous neurological and neuropsychiatric disorders. The GABAB receptor (GABAB-R) is a heterodimeric class C G protein-coupled receptor (GPCR) comprised of GABAB1a/b and GABAB2 subunits. The orthosteric binding site for GABA is located in the extracellular Venus flytrap (VFT) domain of the GABAB1a/b. Knowledge about molecular mechanisms and druggable receptor conformations associated with activation is highly important to understand the receptor function and for rational drug design. Currently, the conformational changes of the receptor upon activation are not well described. On the basis of other class C members, the VFT is proposed to fluctuate between an open/inactive and closed/active state and one of these conformations is stabilized upon ligand binding. In the present study, we investigated the dynamics of the GABAB1b-R VFT in the apo form by combining unbiased molecular dynamics with path-metadynamics. Our simulations confirmed the open/inactive and closed/active state as the main conformations adopted by the receptor. Sizeable energy barriers were found between stable minima, suggesting a relatively slow interconversion. Previously undisclosed metastable states were also identified, which might hold potential for future drug discovery efforts.

Structural Characterization of an ACP from Thermotoga maritima: Insights into Hyperthermal Adaptation

Thermotoga maritima, a deep-branching hyperthermophilic bacterium, expresses an extraordinarily stable Thermotoga maritima acyl carrier protein (Tm-ACP) that functions as a carrier in the fatty acid synthesis system at near-boiling aqueous environments. Here, to understand the hyperthermal adaptation of Tm-ACP, we investigated the structure and dynamics of Tm-ACP by nuclear magnetic resonance (NMR) spectroscopy. The melting temperature of Tm-ACP (101.4 °C) far exceeds that of other ACPs, owing to extensive ionic interactions and tight hydrophobic packing. The D59 residue, which replaces Pro/Ser of other ACPs, mediates ionic clustering between helices III and IV. This creates a wide pocket entrance to facilitate the accommodation of long acyl chains required for hyperthermal adaptation of the T. maritima cell membrane. Tm-ACP is revealed to be the first ACP that harbor an amide proton hyperprotected against hydrogen/deuterium exchange for I15. The hydrophobic interactions mediated by I15 appear to be the key driving forces of the global folding process of Tm-ACP. Our findings provide insights into the structural basis of the hyperthermal adaptation of ACP, which might have allowed T. maritima to survive in hot ancient oceans.

The Effect of Divalent Cations on the Thermostability of Type II Polyketide Synthase Acyl Carrier Proteins

The successful engineering of biosynthetic pathways hinges on understanding the factors that influence acyl carrier protein (ACP) stability and function. The stability and structure of ACPs can be influenced by the presence of divalent cations, but how this relates to primary sequence remains poorly understood. As part of a course-based undergraduate research experience, we investigated the thermostability of type II polyketide synthase (PKS) ACPs. We observed an approximate 40 °C range in the thermostability amongst the 14 ACPs studied, as well as an increase in stability (5 - 26 °C) of the ACPs in the presence of divalent cations. Distribution of charges in the helix II-loop-helix III region was found to impact the enthalpy of denaturation. Taken together, our results reveal clues as to how the sequence of type II PKS ACPs relates to their structural stability, information that can be used to study how ACP sequence relates to function.