Mechanistic insights into ligand dissociation from the SARS-CoV-2 spike glycoprotein.

The COVID-19 pandemic, driven by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has spurred an urgent need for effective therapeutic interventions. The spike glycoprotein of the SARS-CoV-2 is crucial for infiltrating host cells, rendering it a key candidate for drug development. By interacting with the human angiotensin-converting enzyme 2 (ACE2) receptor, the spike initiates the infection of SARS-CoV-2. Linoleate is known to bind the spike glycoprotein, subsequently reducing its interaction with ACE2. However, the detailed mechanisms underlying the protein-ligand interaction remain unclear. In this study, we characterized the pathways of ligand dissociation and the conformational changes associated with the spike glycoprotein by using ligand Gaussian accelerated molecular dynamics (LiGaMD). Our simulations resulted in eight complete ligand dissociation trajectories, unveiling two distinct ligand unbinding pathways. The preference between these two pathways depends on the gate distance between two α-helices in the receptor binding domain (RBD) and the position of the N-linked glycan at N343. Our study also highlights the essential contributions of K417, N121 glycan, and N165 glycan in ligand unbinding, which are equally crucial in enhancing spike-ACE2 binding. We suggest that the presence of the ligand influences the motions of these residues and glycans, consequently reducing accessibility for spike-ACE2 binding. These findings enhance our understanding of ligand dissociation from the spike glycoprotein and offer significant implications for drug design strategies in the battle against COVID-19.

Exploring additives beyond phthalates: Release from plastic mulching films, biodegradation and occurrence in agricultural soils.

It is widely recognized that applications of plastic films result in plastic pollution in agroecosystems. However, there is limited knowledge on the release and occurrence of additives beyond phthalates in agricultural soil. In this study, the rates of release and biodegradation of various additives, including phthalates, bisphenols, organophosphate esters, phenolic antioxidants, and ultraviolet absorbents from mulching films in soil were quantified by laboratory incubation. The rates of release and biodegradation ranged from 0.069 d-1 to 5.893 d-1 and from 1.43 × 10-3 d-1 to 0.600 d-1, respectively. Both of these rates were affected by temperature, flooding, and the properties of additives, films, and soils. An estimated 4000 metric tons of these additives were released into soil annually in China exclusively. The total concentrations of these additives in 80 agricultural soils varied between 228 and 3455 µg kg-1, with phenolic antioxidants, phthalates, and bisphenols accounting for 54.1%, 25.2%, and 17.9% of the total concentrations, respectively. A preliminary risk assessment suggested that the current levels of these additives could potentially present moderate hazards to the soil ecosystem.

Hepatic SEL1L-HRD1 ER-associated degradation regulates systemic iron homeostasis via ceruloplasmin.

Iron homeostasis is critical for cellular and organismal function and is tightly regulated to prevent toxicity or anemia due to iron excess or deficiency, respectively. However, subcellular regulatory mechanisms of iron remain largely unexplored. Here, we report that SEL1L-HRD1 protein complex of endoplasmic reticulum (ER)-associated degradation (ERAD) in hepatocytes controls systemic iron homeostasis in a ceruloplasmin (CP)-dependent, and ER stress-independent, manner. Mice with hepatocyte-specific Sel1L deficiency exhibit altered basal iron homeostasis and are sensitized to iron deficiency while resistant to iron overload. Proteomics screening for a factor linking ERAD deficiency to altered iron homeostasis identifies CP, a key ferroxidase involved in systemic iron distribution by catalyzing iron oxidation and efflux from tissues. Indeed, CP is highly unstable and a bona fide substrate of SEL1L-HRD1 ERAD. In the absence of ERAD, CP protein accumulates in the ER and is shunted to refolding, leading to elevated secretion. Providing clinical relevance of these findings, SEL1L-HRD1 ERAD is responsible for the degradation of a subset of disease-causing CP mutants, thereby attenuating their pathogenicity. Together, this study uncovers the role of SEL1L-HRD1 ERAD in systemic iron homeostasis and provides insights into protein misfolding-associated proteotoxicity.

Molecular dynamics study reveals key disruptors of MEIG1-PACRG interaction.

Interactions between the meiosis-expressed gene 1 (MEIG1) and Parkin co-regulated gene (PACRG) protein are critical in the formation of mature sperm cells. Targeting either MEIG1 or PACRG protein could be a contraceptive strategy. The W50A and Y68A mutations on MEIG1 are known to interrupt the MEIG1-PACRG interactions resulting in defective sperm cells. However, the details about how the mutants disrupt the protein-protein binding are not clear. In this study, we reveal insights on MEIG1 and PACRG protein dynamics by applying Gaussian-accelerated molecular dynamics (GaMD) simulations and post-GaMD analysis. Our results show that the mutations destabilize the protein-protein interfacial interaction. The effect of the Y68A mutation is more significant than W50A as Y68 forms stronger polar interactions with PACRG. Because both human and mouse models demonstrate similar dynamic properties, the findings from mouse proteins can be applied to the human system. Moreover, we report a potential ligand binding pocket on the MEIG1 and PACRG interaction surface that could be a target for future drug design to inhibit the MEIG1-PACRG interaction. PACRG shows more qualified pockets along the protein-protein interface, implying that it is a better target than MEIG1. Our work provides a fundamental understanding of MEIG1 and PACRG protein dynamics, paving the way for drug discovery in male-based contraception.

Molecular Modeling of ABHD5 Structure and Ligand Recognition.

Alpha/beta hydrolase domain-containing 5 (ABHD5), also termed CGI-58, is the key upstream activator of adipose triglyceride lipase (ATGL), which plays an essential role in lipid metabolism and energy storage. Mutations in ABHD5 disrupt lipolysis and are known to cause the Chanarin-Dorfman syndrome. Despite its importance, the structure of ABHD5 remains unknown. In this work, we combine computational and experimental methods to build a 3D structure of ABHD5. Multiple comparative and machine learning-based homology modeling methods are used to obtain possible models of ABHD5. The results from Gaussian accelerated molecular dynamics and experimental data of the apo models and their mutants are used to select the most likely model. Moreover, ensemble docking is performed on representative conformations of ABHD5 to reveal the binding mechanism of ABHD5 and a series of synthetic ligands. Our study suggests that the ABHD5 models created by deep learning-based methods are the best candidate structures for the ABHD5 protein. The mutations of E41, R116, and G328 disturb the hydrogen bonding network with nearby residues and suppress membrane targeting or ATGL activation. The simulations also reveal that the hydrophobic interactions are responsible for binding sulfonyl piperazine ligands to ABHD5. Our work provides fundamental insight into the structure of ABHD5 and its ligand-binding mode, which can be further applied to develop ABHD5 as a therapeutic target for metabolic disease and cancer.

An Integrated Pharmacological, Structural, and Genetic Analysis of Extracellular Versus Intracellular ROS Production in Neutrophils.

The neutrophil NADPH oxidase produces both intracellular and extracellular reactive oxygen species (ROS). Although oxidase activity is essential for microbial killing, and ROS can act as signaling molecules in the inflammatory process, excessive extracellular ROS directly contributes to inflammatory tissue damage, as well as to cancer progression and immune dysregulation in the tumor microenvironment. How specific signaling pathways contribute to ROS localization is unclear. Here we used a systems pharmacology approach to identify the specific Class I PI3-K isoform p110ß, and PLD1, but not PLD2, as critical regulators of extracellular, but not intracellular ROS production in primary neutrophils. Combined crystallographic and molecular dynamics analysis of the PX domain of the oxidase component p47phox, which binds the lipid products of PI 3-K and PLD, was used to clarify the membrane-binding mechanism and guide the design of mutant mice whose p47phox is unable to bind 3-phosphorylated inositol phospholipids. Neutrophils from these K43A mutant animals were specifically deficient in extracellular, but not intracellular, ROS production, and showed increased dependency on signaling through the remaining PLD1 arm. These findings identify the PX domain of p47phox as a critical integrator of PLD1 and p110ß signaling for extracellular ROS production, and as a potential therapeutic target for modulating tissue damage and extracellular signaling during inflammation.

Gaussian accelerated molecular dynamics (GaMD): principles and applications.

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

Multiscale computational study of ligand binding pathways: Case of p38 MAP kinase and its inhibitors.

Protein kinases are one of the most important drug targets in the past 10 years. Understanding the inhibitor association processes will profoundly impact new binder designs with preferred binding kinetics. However, after more than a decade of effort, a complete atomistic-level study of kinase inhibitor binding pathways is still lacking. As all kinases share a similar scaffold, we used p38 kinase as a model system to investigate the conformational dynamics and free energy transition of inhibitor binding toward kinases. Two major kinase conformations, Asp-Phe-Gly (DFG)-in and DFG-out, and three types of inhibitors, type I, II, and III, were thoroughly investigated in this work. We performed Brownian dynamics simulations and up to 340 µs Gaussian-accelerated molecular dynamics simulations to capture the inhibitor binding paths and a series of conformational transitions of the p38 kinase from its apo to inhibitor-bound form. Eighteen successful binding trajectories, including all types of inhibitors, are reported herein. Our simulations suggest a mechanism of inhibitor recruitment, a faster ligand association step to a pre-existing DFG-in/DFG-out p38 protein, followed by a slower molecular rearrangement step to adjust the protein-ligand conformation followed by a shift in the energy landscape to reach the final bound state. The ligand association processes also reflect the energetic favor of type I and type II/III inhibitor binding through ATP and allosteric channels, respectively. These different binding routes are directly responsible for the fast (type I binders) and slow (type II/III binders) kinetics of different types of p38 inhibitors. Our findings also echo the recent study of p38 inhibitor dissociation, implying that ligand unbinding could undergo a reverse path of binding, and both processes share similar metastates. This study deepens the understanding of molecular and energetic features of kinase inhibitor-binding processes and will inspire future drug development from a kinetic point of view.

Enhanced Disrupting Effect of Benzophenone-1 Chlorination Byproducts to the Androgen Receptor: Cell-Based Assays and Gaussian Accelerated Molecular Dynamics Simulations.

Benzophenone-1 (BP-1), one of the commonly used ultraviolet filters, has caused increasing public concern due to frequently detected residues in environmental and recreational waters. Its susceptibility to residual chlorine and the potential to subsequently trigger endocrine disruption remain unknown. We herein investigated the chlorination of BP-1 in swimming pool water and evaluated the endocrine disruption toward the human androgen receptor (AR). The structures of monochlorinated (P1) and dichlorinated (P2) products were separated and characterized by mass spectrometry and 1H-1H NMR correlation spectroscopy. P1 and P2 exhibited significantly higher antiandrogenic activity in yeast two-hybrid assays (EC50, 6.13 µM and 9.30 µM) than did BP-1 (12.89 µM). Our 350 ns Gaussian accelerated molecular dynamics simulations showed the protein dynamics in a long-time scale equilibrium, and further energy calculations revealed that although increased hydrophobic interactions are primarily responsible for enhanced binding affinities between chlorinated products and the AR ligand binding domain, the second chloride in P2 still hinders the complex motion because of the solvation penalty. The mixture of BP-1-P1-P2 elicited additive antiandrogenic activity, well fitted by the concentration addition model. P1 and P2 at 1 µM consequently downregulated the mRNA expression of AR-regulated genes, NKX3.1 and KLK3, by 1.7-9.1-fold in androgen-activated LNCaP cells. Because chlorination of BP-1 occurs naturally by residual chlorine in aquatic environments, our results regarding enhanced antiandrogenic activity and disturbed AR signaling provided evidence linking the use of personal care products with potential health risks.

Elucidating degradation mechanisms of florfenicol in soil by stable-isotope assisted nontarget screening.

Antibiotics in soil environments are a growing concern. Identifying transformation products is key to elucidating degradation pathways and mechanisms of antibiotics and other organic micropollutants. The primary challenge of transformation product identification is the interference of matrices. In this study, stable-isotope assisted nontarget screening was used to identify biodegradation products of florfenicol in soil. A total of 74 candidates were prioritized from thousands of mass features observed by a tiered peak filtering approach. Moreover, with the support of in silico prediction tools, the structures of 12 transformation products were elucidated, and 9 of them were reported for the first time. A biodegradation map of florfenicol consisting of amide hydrolysis, dechlorination, dehydration, defluorination, and sulfone reduction was established based on these identified products. A total of 8 products were also found in 6 field soil samples with manure application. Because of the structural similarity to florfenicol, some transformation products might still keep antimicrobial activity toward a variety of bacterial species. The strategies demonstrated in this study provide a basis for efficient identification of transformation products of other organic micropollutants in a variety of environmental matrices. The results also shed light on the degradation mechanisms, risk assessments, and regulations of these compounds.

Docking simulation and antibiotic discovery targeting the MlaC protein in Gram-negative bacteria.

To maintain the lipid asymmetry of the cell envelope in Gram-negative bacteria, the MlaC protein serves as a lipid transfer factor and delivers phospholipids from the outer to the inner membrane. A strategy of antibiotic discovery is to design a proper compound that can tightly bind to the MlaC protein and inhibit the MlaC function. In this study, we performed virtual screening on multiple MlaC structures obtained from molecular dynamics simulations to identify potential MlaC binders. Our results suggested that clorobiocin is a compound that could bind to the MlaC protein. Through the comparison of the bound geometry between clorobiocin and novobiocin, we pointed out that the methyl-pyrrole group of the noviose sugar in clorobiocin forms hydrophobic interactions with amino acids in the phospholipid binding pocket, which allows the compound to bind deep in the active site. This also explains why clorobiocin shows a tighter binding affinity than novobiocin. Our study highlights a practical path of antibiotic development against Gram-negative bacteria.

Replica Exchange Gaussian Accelerated Molecular Dynamics: Improved Enhanced Sampling and Free Energy Calculation.

Through adding a harmonic boost potential to smooth the system potential energy surface, Gaussian accelerated molecular dynamics (GaMD) provides enhanced sampling and free energy calculation of biomolecules without the need of predefined reaction coordinates. This work continues to improve the acceleration power and energy reweighting of the GaMD by combining the GaMD with replica exchange algorithms. Two versions of replica exchange GaMD (rex-GaMD) are presented: force constant rex-GaMD and threshold energy rex-GaMD. During simulations of force constant rex-GaMD, the boost potential can be exchanged between replicas of different harmonic force constants with fixed threshold energy. However, the algorithm of threshold energy rex-GaMD tends to switch the threshold energy between lower and upper bounds for generating different levels of boost potential. Testing simulations on three model systems, including the alanine dipeptide, chignolin, and HIV protease, demonstrate that through continuous exchanges of the boost potential, the rex-GaMD simulations not only enhance the conformational transitions of the systems but also narrow down the distribution width of the applied boost potential for accurate energetic reweighting to recover biomolecular free energy profiles.

Ligand Binding Pathways and Conformational Transitions of the HIV Protease.

It is important to determine the binding pathways and mechanisms of ligand molecules to target proteins to effectively design therapeutic drugs. Molecular dynamics (MD) is a promising computational tool that allows us to simulate protein-drug binding at an atomistic level. However, the gap between the time scales of current simulations and those of many drug binding processes has limited the usage of conventional MD, which has been reflected in studies of the HIV protease. Here, we have applied a robust enhanced simulation method, Gaussian accelerated molecular dynamics (GaMD), to sample binding pathways of the XK263 ligand and associated protein conformational changes in the HIV protease. During two of 10 independent GaMD simulations performed over 500-2500 ns, the ligand was observed to successfully bind to the protein active site. Although GaMD-derived free energy profiles were not fully converged because of insufficient sampling of the complex system, the simulations still allowed us to identify relatively low-energy intermediate conformational states during binding of the ligand to the HIV protease. Relative to the X-ray crystal structure, the XK263 ligand reached a minimum root-mean-square deviation (RMSD) of 2.26 Å during 2.5 µs of GaMD simulation. In comparison, the ligand RMSD reached a minimum of only ∼5.73 Å during an earlier 14 µs conventional MD simulation. This work highlights the enhanced sampling power of the GaMD approach and demonstrates its wide applicability to studies of drug-receptor interactions for the HIV protease and by extension many other target proteins.

Brownian dynamic study of an enzyme metabolon in the TCA cycle: Substrate kinetics and channeling.

Malate dehydrogenase (MDH) and citrate synthase (CS) are two pacemaking enzymes involved in the tricarboxylic acid (TCA) cycle. Oxaloacetate (OAA) molecules are the intermediate substrates that are transferred from the MDH to CS to carry out sequential catalysis. It is known that, to achieve a high flux of intermediate transport and reduce the probability of substrate leaking, a MDH-CS metabolon forms to enhance the OAA substrate channeling. In this study, we aim to understand the OAA channeling within possible MDH-CS metabolons that have different structural orientations in their complexes. Three MDH-CS metabolons from native bovine, wild-type porcine, and recombinant sources, published in recent work, were selected to calculate OAA transfer efficiency by Brownian dynamics (BD) simulations and to study, through electrostatic potential calculations, a possible role of charges that drive the substrate channeling. Our results show that an electrostatic channel is formed in the metabolons of native bovine and recombinant porcine enzymes, which guides the oppositely charged OAA molecules passing through the channel and enhances the transfer efficiency. However, the channeling probability in a suggested wild-type porcine metabolon conformation is reduced due to an extended diffusion length between the MDH and CS active sites, implying that the corresponding arrangements of MDH and CS result in the decrease of electrostatic steering between substrates and protein surface and then reduce the substrate transfer efficiency from one active site to another.

Mechanism of the Association Pathways for a Pair of Fast and Slow Binding Ligands of HIV-1 Protease.

Equilibrium constants, together with kinetic rate constants of binding, are key factors in the efficacy and safety of drug compounds, informing drug design. However, the association pathways of protein-ligand binding, which contribute to their kinetic behaviors, are little understood. In this work, we used unbiased all-atom molecular dynamics (MD) simulations with an explicit solvent model to study the association processes of protein-ligand binding. Using the HIV protease (HIVp)-xk263 and HIVp-ritonavir protein-ligand systems as cases, we observed that ligand association is a multistep process involving diffusion, localization, and conformational rearrangements of the protein, ligand, and water molecules. Moreover, these two ligands preferred different routes of binding, which reflect two well-known binding mechanisms: induced-fit and conformation selection models. Our study shows that xk263 has a stronger capacity for desolvating surrounding water molecules, thereby inducing a semiopen conformation of the HIVp flaps (induced-fit model). In contrast, the slow dehydration characteristic of ritonavir allows for gradual association with the binding pocket of HIVp when the protein's flap conformation is fully open (conformation selection model). By studying the mechanism of ligand association and understanding the role of solvent molecules during the binding event, we can obtain a different perspective on the mechanism of macromolecule recognition, providing insights into drug discovery.

Investigation of Structural Dynamics of Enzymes and Protonation States of Substrates Using Computational Tools.

This review discusses the use of molecular modeling tools, together with existing experimental findings, to provide a complete atomic-level description of enzyme dynamics and function. We focus on functionally relevant conformational dynamics of enzymes and the protonation states of substrates. The conformational fluctuations of enzymes usually play a crucial role in substrate recognition and catalysis. Protein dynamics can be altered by a tiny change in a molecular system such as different protonation states of various intermediates or by a significant perturbation such as a ligand association. Here we review recent advances in applying atomistic molecular dynamics (MD) simulations to investigate allosteric and network regulation of tryptophan synthase (TRPS) and protonation states of its intermediates and catalysis. In addition, we review studies using quantum mechanics/molecular mechanics (QM/MM) methods to investigate the protonation states of catalytic residues of ß-Ketoacyl ACP synthase I (KasA). We also discuss modeling of large-scale protein motions for HIV-1 protease with coarse-grained Brownian dynamics (BD) simulations.

Characterization of Promiscuous Binding of Phosphor Ligands to Breast-Cancer-Gene 1 (BRCA1) C-Terminal (BRCT): Molecular Dynamics, Free Energy, Entropy and Inhibitor Design.

Inhibition of the protein-protein interaction (PPI) mediated by breast-cancer-gene 1 C-terminal (BRCT) is an attractive strategy to sensitize breast and ovarian cancers to chemotherapeutic agents that induce DNA damage. Such inhibitors could also be used for studies to understand the role of this PPI in DNA damage response. However, design of BRCT inhibitors is challenging because of the inherent flexibility associated with this domain. Several studies identified short phosphopeptides as tight BRCT binders. Here we investigated the thermodynamic properties of 18 phosphopeptides or peptide with phosphate mimic and three compounds with phosphate groups binding to BRCT to understand promiscuous molecular recognition and guide inhibitor design. We performed molecular dynamics (MD) simulations to investigate the interactions between inhibitors and BRCT and their dynamic behavior in the free and bound states. MD simulations revealed the key role of loops in altering the shape and size of the binding site to fit various ligands. The mining minima (M2) method was used for calculating binding free energy to explore the driving forces and the fine balance between configuration entropy loss and enthalpy gain. We designed a rigidified ligand, which showed unfavorable experimental binding affinity due to weakened enthalpy. This was because it lacked the ability to rearrange itself upon binding. Investigation of another phosphate group containing compound, C1, suggested that the entropy loss can be reduced by preventing significant narrowing of the energy well and introducing multiple new compound conformations in the bound states. From our computations, we designed an analog of C1 that introduced new intermolecular interactions to strengthen attractions while maintaining small entropic penalty. This study shows that flexible compounds do not always encounter larger entropy penalty, compared with other more rigid binders, and highlights a new strategy for inhibitor design.

Molecular dynamic study of MlaC protein in Gram-negative bacteria: conformational flexibility, solvent effect and protein-phospholipid binding.

The composition of the outer membrane in Gram-negative bacteria is asymmetric, with the lipopolysaccharides found in the outer leaflet and phospholipids in the inner leaflet. The MlaC protein transfers phospholipids from the outer to inner membrane to maintain such lipid asymmetry in the Mla pathway. In this work, we have performed molecular dynamics simulations on apo and phospholipid-bound systems to study the dynamical properties of MlaC. Our simulations show that the phospholipid forms hydrophobic interactions with the protein. Residues surrounding the entrance of the binding site exhibit correlated motions to control the site opening and closing. Lipid binding leads to increase of the binding pocket volume and precludes entry of the water molecules. However, in the absence of the phospholipid, water molecules can freely move in and out of the binding site when the pocket is open. Dehydration occurs when the pocket closes. This study provides dynamic information of the MlaC protein and may facilitate the design of antibiotics against the Mla pathway of Gram-negative bacteria.

Endosidin2 targets conserved exocyst complex subunit EXO70 to inhibit exocytosis.

The exocyst complex regulates the last steps of exocytosis, which is essential to organisms across kingdoms. In humans, its dysfunction is correlated with several significant diseases, such as diabetes and cancer progression. Investigation of the dynamic regulation of the evolutionarily conserved exocyst-related processes using mutants in genetically tractable organisms such as Arabidopsis thaliana is limited by the lethality or the severity of phenotypes. We discovered that the small molecule Endosidin2 (ES2) binds to the EXO70 (exocyst component of 70 kDa) subunit of the exocyst complex, resulting in inhibition of exocytosis and endosomal recycling in both plant and human cells and enhancement of plant vacuolar trafficking. An EXO70 protein with a C-terminal truncation results in dominant ES2 resistance, uncovering possible distinct regulatory roles for the N terminus of the protein. This study not only provides a valuable tool in studying exocytosis regulation but also offers a potentially new target for drugs aimed at addressing human disease.

Visualizing the tunnel in tryptophan synthase with crystallography: Insights into a selective filter for accommodating indole and rejecting water.

Four new X-ray structures of tryptophan synthase (TS) crystallized with varying numbers of the amphipathic N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F6) molecule are presented. These structures show one of the F6 ligands threaded into the tunnel from the ß-site and reveal a distinct hydrophobic region. Over this expanse, the interactions between F6 and the tunnel are primarily nonpolar, while the F6 phosphoryl group fits into a polar pocket of the ß-subunit active site. Further examination of TS structures reveals that one portion of the tunnel (T1) binds clusters of water molecules, whereas waters are not observed in the nonpolar F6 binding region of the tunnel (T2). MD simulation of another TS structure with an unobstructed tunnel also indicates the T2 region of the tunnel excludes water, consistent with a dewetted state that presents a significant barrier to the transfer of water into the closed ß-site. We conclude that hydrophobic molecules can freely diffuse between the α- and ß-sites via the tunnel, while water does not. We propose that exclusion of water serves to inhibit reaction of water with the α-aminoacrylate intermediate to form ammonium ion and pyruvate, a deleterious side reaction in the αß-catalytic cycle. Finally, while most TS structures show ßPhe280 partially blocking the tunnel between the α- and ß-sites, new structures show an open tunnel, suggesting the flexibility of the ßPhe280 side chain. Flexible docking studies and MD simulations confirm that the dynamic behavior of ßPhe280 allows unhindered transfer of indole through the tunnel, therefore excluding a gating role for this residue.