Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ligand Conformation”

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 19 records

Probing ligand conformation and net dimensionality in a series of tetraphenylethene-based metal–organic frameworks

Tetraphenylethene-based ligands with lowered symmetry are promising building blocks for the construction of novel luminescent metal–organic frameworks (MOFs). However, few examples have been reported, and predicting the ligand conformation and the dimensionality of the resulting MOF remains challenging. In order to uncover how synthetic conditions and accessible ligand conformations may affect the resulting MOF structure, four new MOF structures were synthesized under solvothermal conditions using the meta-coordinated tetraphenylethene-based ligand m-ETTC and paddlewheel SBUs composed of Co(II), Cu(II), and Zn(II). WSU-10 (WSU = Washington State University) is formed with either Zn or Cu comprising stacked psuedo-2D layers. The dimensionality of WSU-10 can be intentionally increased through the addition of pyrazine as a pillar ligand into the synthesis, forming the 3D structure WSU-11. The third structure, WSU-20, is formed by the combination of Zn or Co with m-ETTC and is intrinsically 3D without the use of a pillar ligand; interestingly, this is the result of a distortion in the paddlewheel SBU. Finally, Cu was also found to form a new structure (WSU-12), which displays an m-ETTC conformation unique from that found in the other isolated MOFs. Structural features are compared across the series and a mechanistic relationship between WSU-10 and -20 is proposed, providing insight into the factors that can encourage the generation of frameworks with increased dimensionality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Role of Ligand–Ligand Interactions in Multimodal Ligand Conformational Equilibria and Surface Pattern Formation

Multimodal chromatography uses multiple modes of interaction such as charge, hydrophobic, or hydrogen bonding to separate proteins. Recently, we used molecular dynamics (MD) simulations to show that ligands immobilized on surfaces can interact and associate with neighboring ligands to form hydrophobic and charge patches, which may have important implications for the nature of protein–surface interactions. Here, we study interfacial systems of increasing complexity—from a single immobilized multimodal ligand to high density surfaces—to better understand how ligand behavior is affected by the presence of a surface and the presence of other ligands in the vicinity, and how this behavior scales to larger systems. Furthermore, we find that tethering a ligand to a surface restricts its conformations to a subset of those observed in free solution, yet the ligand maintains flexibility in the plane of the surface and can form contacts with neighboring ligands. We find that although the formation of a contact between two neighboring ligands is slightly unfavorable, three neighboring ligands exhibit a preference for the formation of a fully connected cluster. To explore how these trends in ligand association extend to a larger surface with high density of ligands, we performed coarse-grained Monte Carlo (MC) simulations of a 132-ligand surface using ligand interactions parametrized based on free energies obtained from the three-ligand MD simulations. Despite their simplicity, the coarse-grained simulations qualitatively capture the cluster size distribution of ligands observed in detailed MD simulations. Quantitative differences between the two suggest opportunities for improvements in the coarse-grained energy function for efficient predictions of cluster and pattern formations. Our approach presents a promising route to the engineering of multimodal patterns for future chromatographic resin design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conformational Equilibria of Multimodal Chromatography Ligands in Water and Bound to Protein Surfaces

Multimodal chromatography uses small ligands with multiple modes of interaction, e.g., charged, hydrophobic or hydrogen bonding, to separate proteins from complex mixtures. The mechanism by which multimodal ligands interact with proteins is expected to be affected by ligand conformations, among other factors. In this report, we determine conformational equilibria of two commercially used multimodal cation exchange ligands, Capto MMC and Nuvia cPrime, in a range of solvents, a Lennard-Jones (LJ) liquid, ethanol, and water, using molecular dynamics (MD) simulations. By mapping ligand conformations onto two key torsion angles, ω and φ, in these solvents and in low and high dielectric media, we quantify the relative importance of intramolecular and solvent-mediated interactions. In a high dielectric medium, Capto MMC preferentially samples three conformations, which are stabilized by a combination of an intramolecular torsion potential (on ω) and LJ interactions. In an LJ liquid, solvent molecules compete with intramolecular interactions while simultaneously providing an osmotic force, stabilizing both closer and farther distances between ligand sites. This has the overall result of “flattening out” the conformational landscape. Interestingly, in ethanol and water, hydrogen bonding among the amide hydrogen and solvent molecules stabilizes two additional conformations of Capto MMC in which ω takes on less favorable cis-like configurations. MD simulations of ligands in free solution with three therapeutic antibody fragments show that ligand conformational equilibria remain effectively unchanged upon binding to proteins. Although, there is 20–30% dehydration of the overall ligand upon binding, the hydrogen-bonding sites are dehydrated to a much smaller extent, particularly in cis-like configurations. Conformational preferences of Nuvia cPrime are similar to that of Capto MMC, except for the effect of symmetry arising from the absence of an alkyl thiol tail. Characterizing the conformational equilibria of these two ligands in free solution and bound to a protein provides a foundation for developing a mechanistic understanding of protein–multimodal ligand interactions.

59 BASIC BIOLOGICAL SCIENCES↗

Measuring Relative Energies of Ligand Binding Conformations on Nanocluster Surfaces with Temperature-Dependent FTIR Spectroscopy

We present a method to measure relative energies between binding conformations of carboxylate ligands on InP magic-sized clusters in solution. Using Markov chain Monte Carlo global fitting analysis on temperature-dependent vibrational spectra of cluster-bound ligands, we observe significantly different relative energies between various bidentate and monodentate binding motifs. Relative to the monodentate motif, the chelating conformation is 0.7 ± 0.3 kcal/mol more stable and the syn–syn bridging conformation is 1.1 ± 0.5 kcal/mol more stable, but the syn–anti bridging conformation exhibits no significant difference. Our results demonstrate that the relative energy between monodentate-bound carboxylates and unbound carboxylic acids is 4.52 ± 0.05 kcal/mol, or 1582 ± 19 cm –1 , nearly identical to the carboxylate asymmetric stretching frequency. Here, we suggest that the ligand vibrational energy may play a key role in ligand dissociation by compensating for energy differences between bound and dissociated ligand states. This approach gives important experimental insights into ligand binding and can inform future nanocrystal surface engineering.

Conformation↗

Structural basis of the UDP-diacylglucosamine pyrophosphohydrolase LpxH inhibition by sulfonyl piperazine antibiotics

The UDP-2,3-diacylglucosamine pyrophosphate hydrolase LpxH is an essential lipid A biosynthetic enzyme that is conserved in the majority of gram-negative bacteria. It has emerged as an attractive novel antibiotic target due to the recent discovery of an LpxH-targeting sulfonyl piperazine compound (referred to as AZ1) by AstraZeneca. However, the molecular details of AZ1 inhibition have remained unresolved, stymieing further development of this class of antibiotics. In this paper we report the crystal structure of Klebsiella pneumoniae LpxH in complex with AZ1. We show that AZ1 fits snugly into the L -shaped acyl chain-binding chamber of LpxH with its indoline ring situating adjacent to the active site, its sulfonyl group adopting a sharp kink, and its N -CF 3 –phenyl substituted piperazine group reaching out to the far side of the LpxH acyl chain-binding chamber. Intriguingly, despite the observation of a single AZ1 conformation in the crystal structure, our solution NMR investigation has revealed the presence of a second ligand conformation invisible in the crystalline state. Together, these distinct ligand conformations delineate a cryptic inhibitor envelope that expands the observed footprint of AZ1 in the LpxH-bound crystal structure and enables the design of AZ1 analogs with enhanced potency in enzymatic assays. These designed compounds display striking improvement in antibiotic activity over AZ1 against wild-type K. pneumoniae , and coadministration with outer membrane permeability enhancers profoundly sensitizes Escherichia coli to designed LpxH inhibitors. Remarkably, none of the sulfonyl piperazine compounds occupies the active site of LpxH, foretelling a straightforward path for rapid optimization of this class of antibiotics.

59 BASIC BIOLOGICAL SCIENCES↗

X-ray crystal structure of a designed rigidified imaging scaffold in the ligand-free conformation

Imaging scaffolds composed of designed protein cages fused to designed ankyrin repeat proteins (DARPins) have enabled the structure determination of small proteins by cryogenic electron microscopy (cryo-EM). One particularly well characterized scaffold type is a symmetric tetrahedral assembly composed of 24 subunits, 12 A and 12 B, which has three cargo-binding DARPins positioned on each vertex. Here, the X-ray crystal structure of a representative tetrahedral scaffold in the apo state is reported at 3.8 Å resolution. The X-ray crystal structure complements recent cryo-EM findings on a closely related scaffold, while also suggesting potential utility for crystallographic investigations. As observed in this crystal structure, one of the three DARPins, which serve as modular adaptors for binding diverse `cargo' proteins, present on each of the vertices is oriented towards a large solvent channel. The crystal lattice is unusually porous, suggesting that it may be possible to soak crystals of the scaffold with small (≤30 kDa) protein cargo ligands and subsequently determine cage–cargo structures via X-ray crystallography. The results suggest the possibility that cryo-EM scaffolds may be repurposed for structure determination by X-ray crystallography, thus extending the utility of electron-microscopy scaffold designs for alternative structural biology applications.

59 BASIC BIOLOGICAL SCIENCES↗

Structural Characterization and Ligand-Induced Conformational Changes of SenB, a Se-Glycosyltransferase Involved in Selenoneine Biosynthesis

Selenium (Se) is an essential micronutrient that is found naturally in proteins, nucleic acids, and natural products. Unlike selenoproteins and selenonucleic acids, little is known about the structures of biosynthetic enzymes that incorporate Se into small molecules. Here, we report the X-ray crystal structure of SenB, the first known Se-glycosyltransferase that was recently found to be involved in the biosynthesis of the Se-containing metabolite selenoneine. SenB catalyzes C–Se bond formation using selenophosphate and an activated uridine diphosphate sugar as a Se and glycosyl donor, respectively, making it the first known selenosugar synthase and one of only four bona fide C–Se bond-forming enzymes discovered to date. Our crystal structure, determined to 2.25 Å resolution, reveals that SenB is a type B glycosyltransferase, displaying the prototypical fold with two globular Rossmann-like domains and a catalytic interdomain cleft. By employing complementary structural biology techniques, we find that SenB undergoes both local and global substrate-induced conformational changes, demonstrating a significant increase in α-helicity and a transition to a more compact conformation. Our results provide the first structure of SenB and set the stage for further biochemical characterization in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

qFit-ligand Reveals Widespread Conformational Heterogeneity of Drug-Like Molecules in X-Ray Electron Density Maps

Here, proteins and ligands sample a conformational ensemble that governs molecular recognition, activity, and dissociation. In structure-based drug design, access to this conformational ensemble is critical to understand the balance between entropy and enthalpy in lead optimization. However, ligand conformational heterogeneity is currently severely underreported in crystal structures in the Protein Data Bank, owing in part to a lack of automated and unbiased procedures to model an ensemble of protein–ligand states into X-ray data. Here, we designed a computational method, qFit-ligand , to automatically resolve conformationally averaged ligand heterogeneity in crystal structures, and applied it to a large set of protein receptor–ligand complexes. In an analysis of the cancer related BRD4 domain, we found that up to 29% of protein crystal structures bound with drug-like molecules present evidence of unmodeled, averaged, relatively isoenergetic conformations in ligand–receptor interactions. In many retrospective cases, these alternate conformations were adventitiously exploited to guide compound design, resulting in improved potency or selectivity. Combining qFit-ligand with high-throughput screening or multitemperature crystallography could therefore augment the structure-based drug design toolbox.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Engineered synthetic antibodies as probes to quantify the energetic contributions of ligand binding to conformational changes in proteins

Conformational changes in proteins due to ligand binding are ubiquitous in biological processes and are integral to many biological systems. However, it is often challenging to link ligand-induced conformational changes to a resulting biological function because it is difficult to distinguish between the energetic components associated with ligand binding and those due to structural rearrangements. Consequently, here, we used a unique approach exploiting conformation-specific and regio-specific synthetic antibodies (sABs) to probe the energetic contributions of ligand binding to conformation changes. Using maltose-binding protein (MBP) as a model system, customized phage-display selections were performed to generate sABs that stabilize MBP in different conformational states, modulating ligand-binding affinity in competitive, allosteric, or peristeric manners. We determined that the binding of a closed conformation–specific sAB (sAB-11M) to MBP in the absence of maltose is entropically driven, providing new insight into designing antibody-stabilized protein interactions. Crystal structures of sABs bound to MBP, together with biophysical data, delineate the basis of free energy differences between different conformational states and confirm the use of the sABs as energy probes for dissecting enthalpic and entropic contributions to conformational transitions. Our work provides a foundation for investigating the energetic contributions of distinct conformational dynamics to specific biological outputs. We anticipate that our approach also may be valuable for analyzing the energy landscapes of regulatory proteins controlling biological responses to environmental changes.

59 BASIC BIOLOGICAL SCIENCES↗

Effects of Catalytic Action and Ligand Binding on Conformational Ensembles of Adenylate Kinase

Crystal structures of adenylate kinase (AdK) from Escherichia coli capture two states: an “open” conformation (apo) obtained in the absence of ligands and a “closed” conformation in which ligands are bound. Other AdK crystal structures suggest intermediate conformations that may lie on the transition pathway between these two states. To characterize the transition from open to closed states in solution, X-ray solution scattering data were collected from AdK in the apo form and with progressively increasing concentrations of five different ligands. Scattering data from apo AdK are consistent with scattering predicted from the crystal structure of AdK in the open conformation. In contrast, data from AdK samples saturated with Ap5A do not agree with that calculated from AdK in the closed conformation. Using cluster analysis of available structures, we selected representative structures in five conformational states: open, partially open, intermediate, partially closed, and closed. We used these structures to estimate the relative abundances of these states for each experimental condition. X-ray solution scattering data obtained from AdK with AMP are dominated by scattering from AdK in the open conformation. For AdK in the presence of high concentrations of ATP and ADP, the conformational ensemble shifts to a mixture of partially open and closed states. Even when AdK is saturated with Ap5A, a significant proportion of AdK remains in a partially open conformation. Here, these results are consistent with an induced-fit model in which the transition of AdK from an open state to a closed state is initiated by ATP binding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Speed vs Accuracy: Effect on Ligand Pose Accuracy of Varying Box Size and Exhaustiveness in AutoDock Vina

Abstract Structure‐based virtual high‐throughput screening involves docking chemical libraries to targets of interest. A parameter pertinent to the accuracy of the resulting pose is the root mean square deviation (RMSD) from a known crystallographic structure, i. e., the ‘docking power’. Here, using a popular algorithm, Autodock Vina, as a model program, we evaluate the effects of varying two common docking parameters: the box size (the size of docking search space) and the exhaustiveness of the global search (the number of independent runs starting from random ligand conformations) on the RMSD from the PDBbind v2017 refined dataset of experimental protein‐ligand complexes. Although it is clear that exhaustiveness is an important parameter, there is wide variation in the values used, with variation between 1 and >100. We, therefore, evaluated a combination of cubic boxes of different sizes and five exhaustiveness values (1, 8, 25, 50, 75, 100) within the range of those commonly adopted. The results show that the default exhaustiveness value of 8 performs well overall for most box sizes. In contrast, for all box sizes, but particularly for large boxes, an exhaustiveness value of 1 led to significantly higher median RMSD (mRMSD) values. The docking power was slightly improved with an exhaustiveness of 25, but the mRMSD changes little with values higher than 25. Therefore, although low exhaustiveness is computationally faster, the results are more likely to be far from reality, and, conversely, values >25 led to little improvement at the expense of computational resources. Overall, we recommend users to use at least the default exhaustiveness value of 8 for virtual screening calculations.

Agarwal, Rupesh↗

tinyIFD: A High-Throughput Binding Pose Refinement Workflow Through Induced-Fit Ligand Docking

A critical step in structure-based drug discovery is predicting whether and how a candidate molecule binds to a model of a therapeutic target. However, substantial protein side chain movements prevent current screening methods, such as docking, from accurately predicting the ligand conformations and require expensive refinements to produce viable candidates. Here, we present the development of a high-throughput and flexible ligand pose refinement workflow, called “tinyIFD”. The main features of the workflow include the use of specialized high-throughput, small-system MD simulation code mdgx.cuda and an actively learning model zoo approach. We show the application of this workflow on a large test set of diverse protein targets, achieving 66% and 76% success rates for finding a crystal-like pose within the top-2 and top-5 poses, respectively. We also applied this workflow to the SARS-CoV-2 main protease (M pro ) inhibitors, where we demonstrate the benefit of the active learning aspect in this workflow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uranyl Naphthylsalophen and Pyrasal Complexes: Oxo Ligands Acting as Hydrogen Bond Acceptors in the Solid State

Uranium is most stable when it is exposed to oxygen or water in its +6 oxidation state as the uranyl (UO 2 2+ ) ion. This ion is subsequently particularly stable and very resistant to functionalization due to the inverse trans effect. Uranyl oxo ligands are typically not considered good hydrogen bond acceptors due to their weak Lewis basicity; however, the ligands bound in the equatorial plane greatly affect the strength of the oxo ligands’ hydrogen bonding. Here, in this work, new naphthylsalophen and pyrasal complexes of uranium were synthesized and crystallized for characterization in the solid state. The bond lengths and angles of the uranyl ion and the ligand conformation are compared. In the solid state, one of the pyrasal complexes showed a hydrogen bond directly from a water molecule to the uranyl oxo ligand, which resulted in an asymmetric lengthening of the U–O yl bonds from 1.789 to 1.862 Å and 1.784 to 1.844 Å.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ligand-induced transmembrane conformational coupling in monomeric $\mathrm{EGFR}$

Single pass cell surface receptors regulate cellular processes by transmitting ligand-encoded signals across the plasma membrane via changes to their extracellular and intracellular conformations. This transmembrane signaling is generally initiated by ligand binding to the receptors in their monomeric form. While subsequent receptor-receptor interactions are established as key aspects of transmembrane signaling, the contribution of monomeric receptors has been challenging to isolate due to the complexity and ligand-dependence of these interactions. By combining membrane nanodiscs produced with cell-free expression, single-molecule Förster Resonance Energy Transfer measurements, and molecular dynamics simulations, we report that ligand binding induces intracellular conformational changes within monomeric, full-length epidermal growth factor receptor (EGFR). Our observations establish the existence of extracellular/intracellular conformational coupling within a single receptor molecule. We implicate a series of electrostatic interactions in the conformational coupling and find the coupling is inhibited by targeted therapeutics and mutations that also inhibit phosphorylation in cells. Collectively, these results introduce a facile mechanism to link the extracellular and intracellular regions through the single transmembrane helix of monomeric EGFR, and raise the possibility that intramolecular transmembrane conformational changes upon ligand binding are common to single-pass membrane proteins.

59 BASIC BIOLOGICAL SCIENCES↗

SARS-CoV2 billion-compound docking

Abstract This dataset contains ligand conformations and docking scores for 1.4 billion molecules docked against 6 structural targets from SARS-CoV2, representing 5 unique proteins: MPro, NSP15, PLPro, RDRP, and the Spike protein. Docking was carried out using the AutoDock-GPU platform on the Summit supercomputer and Google Cloud. The docking procedure employed the Solis Wets search method to generate 20 independent ligand binding poses per compound. Each compound geometry was scored using the AutoDock free energy estimate, and rescored using RFScore v3 and DUD-E machine-learned rescoring models. Input protein structures are included, suitable for use by AutoDock-GPU and other docking programs. As the result of an exceptionally large docking campaign, this dataset represents a valuable resource for discovering trends across small molecule and protein binding sites, training AI models, and comparing to inhibitor compounds targeting SARS-CoV-2. The work also gives an example of how to organize and process data from ultra-large docking screens.

60 APPLIED LIFE SCIENCES↗

Structural characterization of ligand binding and pH-specific enzymatic activity of mouse Acidic Mammalian Chitinase

Chitin is an abundant biopolymer and pathogen-associated molecular pattern that stimulates a host innate immune response. Mammals express chitin-binding and chitin-degrading proteins to remove chitin from the body. One of these proteins, Acidic Mammalian Chitinase (AMCase), is an enzyme known for its ability to function under acidic conditions in the stomach but is also active in tissues with more neutral pHs, such as the lung. Here, we used a combination of biochemical, structural, and computational modeling approaches to examine how the mouse homolog (mAMCase) can act in both acidic and neutral environments. We measured kinetic properties of mAMCase activity across a broad pH range, quantifying its unusual dual activity optima at pH 2 and 7. We also solved high-resolution crystal structures of mAMCase in complex with oligomeric GlcNAcn, the building block of chitin, where we identified extensive conformational ligand heterogeneity. Leveraging these data, we conducted molecular dynamics simulations that suggest how a key catalytic residue could be protonated via distinct mechanisms in each of the two environmental pH ranges. These results integrate structural, biochemical, and computational approaches to deliver a more complete understanding of the catalytic mechanism governing mAMCase activity at different pH. Engineering proteins with tunable pH optima may provide new opportunities to develop improved enzyme variants, including AMCase, for therapeutic purposes in chitin degradation.

59 BASIC BIOLOGICAL SCIENCES↗