Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

AI-Accelerated Design of Targeted Covalent Inhibitors for SARS-CoV-2

Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic’s evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (M pro ), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. Here, this work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic “warhead” to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of M pro with micromolar affinities (K I of 5.27 μM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

SAXS analysis of the intrinsic tenase complex bound to a lipid nanodisc highlights intermolecular contacts between factors VIIIa/IXa

Abstract The intrinsic tenase (Xase) complex, formed by factors (f) VIIIa and fIXa, forms on activated platelet surfaces and catalyzes the activation of factor X to Xa, stimulating thrombin production in the blood coagulation cascade. The structural organization of the membrane-bound Xase complex remains largely unknown, hindering our understanding of the structural underpinnings that guide Xase complex assembly. Here, we aimed to characterize the Xase complex bound to a lipid nanodisc with biolayer interferometry (BLI), Michaelis–Menten kinetics, and small-angle X-ray scattering (SAXS). Using immobilized lipid nanodiscs, we measured binding rates and nanomolar affinities for fVIIIa, fIXa, and the Xase complex. Enzyme kinetic measurements demonstrated the assembly of an active enzyme complex in the presence of lipid nanodiscs. An ab initio molecular envelope of the nanodisc-bound Xase complex allowed us to computationally model fVIIIa and fIXa docked onto a flexible lipid membrane and identify protein–protein interactions. Our results highlight multiple points of contact between fVIIIa and fIXa, including a novel interaction with fIXa at the fVIIIa A1–A3 domain interface. Lastly, we identified hemophilia A/B-related mutations with varying severities at the fVIIIa/fIXa interface that may regulate Xase complex assembly. Together, our results support the use of SAXS as an emergent tool to investigate the membrane-bound Xase complex and illustrate how mutations at the fVIIIa/fIXa dimer interface may disrupt or stabilize the activated enzyme complex.

Hematology↗
DOE OSTI2020

D3R grand challenge 4: blind prediction of protein–ligand poses, affinity rankings, and relative binding free energies

The Drug Design Data Resource (D3R) aims to identify best practice methods for computer aided drug design through blinded ligand pose prediction and affinity challenges. Herein, we report on the results of Grand Challenge 4 (GC4). GC4 focused on proteins beta secretase 1 and Cathepsin S, and was run in an analogous manner to prior challenges. In Stage 1, participant ability to predict the pose and affinity of BACE1 ligands were assessed. Following the completion of Stage 1, all BACE1 co-crystal structures were released, and Stage 2 tested affinity rankings with co-crystal structures. Finally, we provide an analysis of the results and discuss insights into determined best practice methods.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Revising the Role of Cortical Cytoskeleton during Secretion: Actin and Myosin XI Function in Vesicle Tethering

In plants, secretion of cell wall components and membrane proteins plays a fundamental role in growth and development as well as survival in diverse environments. Exocytosis, as the last step of the secretory trafficking pathway, is a highly ordered and precisely controlled process involving tethering, docking, and fusion of vesicles at the plasma membrane (PM) for cargo delivery. Although the exocytic process and machinery are well characterized in yeast and animal models, the molecular players and specific molecular events that underpin late stages of exocytosis in plant cells remain largely unknown. Here, by using the delivery of functional, fluorescent-tagged cellulose synthase (CESA) complexes (CSCs) to the PM as a model system for secretion, as well as single-particle tracking in living cells, we describe a quantitative approach for measuring the frequency of vesicle tethering events. Genetic and pharmacological inhibition of cytoskeletal function, reveal that the initial vesicle tethering step of exocytosis is dependent on actin and myosin XI. In contrast, treatments with the microtubule inhibitor, oryzalin, did not significantly affect vesicle tethering or fusion during CSC exocytosis but caused a minor increase in transient or aborted tethering events. With data from this new quantitative approach and improved spatiotemporal resolution of single particle events during secretion, we generate a revised model for the role of the cortical cytoskeleton in CSC trafficking.

Zhang, Weiwei↗
DOE OSTI2023

Structural insights into hydrolytic defluorination of difluoroacetate by microbial fluoroacetate dehalogenases

Fluorine forms the strongest single bond to carbon with the highest bond dissociation energy among natural products. However, fluoroacetate dehalogenases (FADs) have been shown to hydrolyze this bond in fluoroacetate under mild reaction conditions. Furthermore, two recent studies demonstrated that the FAD RPA1163 fromRhodopseudomonas palustriscan also accept bulkier substrates. In this study, we explored the substrate promiscuity of microbial FADs and their ability to defluorinate polyfluorinated organic acids. Enzymatic screening of eight purified dehalogenases with reported fluoroacetate defluorination activity revealed significant hydrolytic activity against difluoroacetate in three proteins. Product analysis using liquid chromatography‐mass spectrometry identified glyoxylic acid as the final product of enzymatic DFA defluorination. The crystal structures of DAR3835 fromDechloromonas aromaticaand NOS0089 fromNostocsp. were determined in the apo‐state along with the DAR3835 H274N glycolyl intermediate. Structure‐based site‐directed mutagenesis of DAR3835 demonstrated a key role for the catalytic triad and other active site residues in the defluorination of both fluoroacetate and difluoroacetate. Computational analysis of the dimer structures of DAR3835, NOS0089, and RPA1163 indicated the presence of one substrate access tunnel in each protomer. Moreover, protein‐ligand docking simulations suggested similar catalytic mechanisms for the defluorination of both fluoroacetate and difluoroacetate, with difluoroacetate being defluorinated via two consecutive defluorination reactions producing glyoxylate as the final product. Thus, our findings provide molecular insights into substrate promiscuity and catalytic mechanism of FADs, which are promising biocatalysts for applications in synthetic chemistry and bioremediation of fluorochemicals.

Biochemistry & Molecular Biology↗
DOE OSTI2022

Data Science-Driven Analysis of Substrate-Permissive Diketopiperazine Reverse Prenyltransferase NotF: Applications in Protein Engineering and Cascade Biocatalytic Synthesis of (-)-Eurotiumin A

Prenyltransfer is an early-stage carbon-hydrogen bond (C-H) functionalization prevalent in the biosynthesis of a diverse array of biologically active bacterial, fungal, plant, and metazoan diketopiperazine (DKP) alkaloids. Toward the development of a unified strategy for biocatalytic construction of prenylated DKP indole alkaloids, we sought to identify and characterize a substrate-permissive C2 reverse prenyltransferase (PT). As the first tailoring event within the biosynthesis of cytotoxic notoamide metabolites, PT NotF catalyzes C2 reverse prenyltransfer of brevianamide F. Solving a crystal structure of NotF (in complex with native substrate and prenyl donor mimic dimethylallyl S-thiolodiphosphate (DMSPP)) revealed a large, solvent-exposed active site, intimating NotF may possess a significantly broad substrate scope. To assess the substrate selectivity of NotF, we synthesized a panel of 30 sterically and electronically differentiated tryptophanyl DKPs, the majority of which were selectively prenylated by NotF in synthetically useful conversions (2 to > 99%). Quantitative representation of this substrate library and development of a descriptive statistical model provided insight into the molecular origins of NotF's substrate promiscuity. This approach enabled the identification of key substrate descriptors (electrophilicity, size, and flexibility) that govern the rate of NotF-catalyzed prenyltransfer, and the development of an "induced fit docking (IFD)-guided" engineering strategy for improved turnover of our largest substrates. We further demonstrated the utility of NotF in tandem with oxidative cyclization using flavin monooxygenase, BvnB. This one-pot, in vitro biocatalytic cascade enabled the first chemoenzymatic synthesis of the marine fungal natural product, (-)-eurotiumin A, in three steps and 60% overall yield.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

CoarsenConf: Equivariant Coarsening with Aggregated Attention for Molecular Conformer Generation

Molecular conformer generation (MCG) is an important task in cheminformatics and drug discovery. The ability to efficiently generate low-energy 3D structures can avoid expensive quantum mechanical simulations, leading to accelerated virtual screenings and enhanced structural exploration. Several generative models have been developed for MCG, but many struggle to consistently produce high-quality conformers for meaningful downstream applications. To address these issues, we introduce CoarsenConf, which coarse-grains molecular graphs based on torsional angles and integrates them into an SE(3)-equivariant hierarchical variational autoencoder. Through equivariant coarse-graining, we aggregate the fine-grained atomic coordinates of subgraphs connected via rotatable bonds, creating a variable-length coarse-grained latent representation. Our model uses a novel aggregated attention mechanism to restore fine-grained coordinates from the coarse-grained latent representation, enabling efficient generation of accurate conformers. Furthermore, we evaluate the chemical and biochemical quality of our generated conformers on multiple downstream applications, including property prediction and large-scale oracle-based protein docking. Overall, CoarsenConf generates more accurate conformer ensembles compared to prior generative models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Staphylococcus aureus vWF-binding protein triggers a strong interaction between clumping factor A and host vWF

The Staphylococcus aureus cell wall-anchored adhesin ClfA binds to the very large blood circulating protein, von Willebrand factor (vWF) via vWF-binding protein (vWbp), a secreted protein that does not bind the cell wall covalently. Here we perform force spectroscopy studies on living bacteria to unravel the molecular mechanism of this interaction. We discover that the presence of all three binding partners leads to very high binding forces (2000 pN), largely outperforming other known ternary complexes involving adhesins. Strikingly, our experiments indicate that a direct interaction involving features of the dock, lock and latch mechanism must occur between ClfA and vWF to sustain the extreme tensile strength of the ternary complex. Our results support a previously undescribed mechanism whereby vWbp activates a direct, ultra-strong interaction between ClfA and vWF. This intriguing interaction represents a potential target for therapeutic interventions, including synthetic peptides inhibiting the ultra-strong interactions between ClfA and its ligands.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

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↗
DOE OSTI2022

BIPSPI+: Mining Type-Specific Datasets of Protein Complexes to Improve Protein Binding Site Prediction

Computational approaches for predicting protein-protein interfaces are extremely useful for understanding and modelling the quaternary structure of protein assemblies. In particular, partner-specific binding site prediction methods allow delineating the specific residues that compose the interface of protein complexes. In recent years, new machine learning and other algorithmic approaches have been proposed to solve this problem. However, little effort has been made in finding better training datasets to improve the performance of these methods. With the aim of vindicating the importance of the training set compilation procedure, in this work we present BIPSPI+, a new version of our original server trained on carefully curated datasets that outperforms our original predictor. We show how prediction performance can be improved by selecting specific datasets that better describe particular types of protein interactions and interfaces (e.g. homo/hetero). In addition, our upgraded web server offers a new set of functionalities such as the sequence-structure prediction mode, hetero- or homo-complex specialization and the guided docking tool that allows to compute 3D quaternary structure poses using the predicted interfaces. BIPSPI+ is freely available at https://bipspi.cnb.csic.es.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Dirigent isoflavene-forming PsPTS2: 3D structure, stereochemical, and kinetic characterization comparison with pterocarpan-forming PsPTS1 homolog in pea

Pea (legume plant defense) phytoalexin pterocarpans (–)-maackiain and (+)-pisatin have opposite C6a and C11a configurations, but how this occurs biosynthetically is unknown. Herein, pea dirigent-protein (DP) PsPTS2, which generates 7,2'-dihydroxy-4',5'-methylenedioxyisoflav-3-ene (DMDIF), was investigated as regards stereoselectivity towards four possible 7,2'-dihydroxy-4',5'-methylenedioxyisoflavan-4-ol (DMDI) stereoisomers. Each DMDI stereoisomer configuration was determined using NMR spectroscopy, electronic circular dichroism, and molecular orbital analyses. PsPTS2 efficiently converted cis-(3R,4R)-DMDI into DMDIF, this being 20-fold faster than with the trans-(3R,4S)-isomer. The 4R-configured substrate’s near ß-axial OH orientation significantly enhanced its leaving group abilities in creating intermediate A-ring mono-quinone methide (QM), whereas the 4S-isomer’s 4-OH group α-equatorial orientation was a poorer leaving group. Substrate docking simulations also indicated that the 4R-configured ß-axial OH was closest to the Asp catalytic center, whereas the 4S-isomer’s α-equatorial OH was further away. Neither cis-(3S,4S)- nor trans-(3S, 4R)-DMDIs were substrates, even with the former having the same C3/C4 stereochemistry as (+)-pisatin. PsPTS2, assayed with cis-(3R,4R)-7,2'-dihydroxy-4'-methoxyisoflavan-4-ol [cis-(3R,4R)-DMI] and its C3/C4 stereoisomers, gave 2',7-dihydroxy-4'-methoxyisoflav-3-ene (DMIF), with the same substrate stereoselectivity. DP homologs may exist in licorice (Glycyrrhiza pallidiflora) and the tree legume Bolusanthus speciosus, as DMIF is found in both. PsPTS1 preferentially utilized cis-(3R,4R)-DMDI to give (–)-maackiain, this being 2200-fold more efficient than with cis-(3R,4R)-DMI. PsPTS1 also slowly converted trans-(3S,4R)-DMDI into (+)-maackiain, again reflecting the better 4R configured OH leaving group. Pea PsPTS2 and PsPTS1 provisionally provide a facile DP-engendered means to enable differing C6a and C11a configurations in (+)-pisatin and (–)-maackiain. Both DP transformations are considered to occur via DP-engendered generation of an identical mono-QM bound intermediate, which with PsPTS2 either re-aromatizes to give DMDIF or with PsPTS1 undergoes intramolecular cyclization to afford (–)-maackiain. Substrate docking simulations using PsPTS2 (2.5 Å resolution) and PsPTS1 (1.5 Å resolution) plausibly indicate that cis-(3R,4R)-DMDI binds in the anti-configuration in the PsPTS2 active site to afford DMDIF, and the syn-configuration in PsPTS1 to give maackiain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

CryoEM structures of the human CLC-2 voltage-gated chloride channel reveal a ball-and-chain gating mechanism

CLC-2 is a voltage-gated chloride channel that contributes to electrical excitability and ion homeostasis in many different tissues. Among the nine mammalian CLC homologs, CLC-2 is uniquely activated by hyperpolarization, rather than depolarization, of the plasma membrane. The molecular basis for the divergence in polarity of voltage gating among closely related homologs has been a long-standing mystery, in part because few CLC channel structures are available. Here, we report cryoEM structures of human CLC-2 at 2.46 – 2.76 Å, in the presence and absence of the selective inhibitor AK-42. AK-42 binds within the extracellular entryway of the Cl – -permeation pathway, occupying a pocket previously proposed through computational docking studies. In the apo structure, we observed two distinct conformations involving rotation of one of the cytoplasmic C-terminal domains (CTDs). In the absence of CTD rotation, an intracellular N-terminal 15-residue hairpin peptide nestles against the TM domain to physically occlude the Cl – -permeation pathway. This peptide is highly conserved among species variants of CLC-2 but is not present in other CLC homologs. Previous studies suggested that the N-terminal domain of CLC-2 influences channel properties via a “ball-and-chain” gating mechanism, but conflicting data cast doubt on such a mechanism, and thus the structure of the N-terminal domain and its interaction with the channel has been uncertain. Through electrophysiological studies of an N-terminal deletion mutant lacking the 15-residue hairpin peptide, we support a model in which the N-terminal hairpin of CLC-2 stabilizes a closed state of the channel by blocking the cytoplasmic Cl – -permeation pathway.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structural basis of Gabija anti-phage defence and viral immune evasion

Bacteria encode hundreds of diverse defence systems that protect them from viral infection and inhibit phage propagation. Gabija is one of the most prevalent anti-phage defence systems, occurring in more than 15% of all sequenced bacterial and archaeal genomes, but the molecular basis of how Gabija defends cells from viral infection remains poorly understood. Here we use X-ray crystallography and cryo-electron microscopy (cryo-EM) to define how Gabija proteins assemble into a supramolecular complex of around 500 kDa that degrades phage DNA. Gabija protein A (GajA) is a DNA endonuclease that tetramerizes to form the core of the anti-phage defence complex. Two sets of Gabija protein B (GajB) dimers dock at opposite sides of the complex and create a 4:4 GajA–GajB assembly (hereafter, GajAB) that is essential for phage resistance in vivo. We show that a phage-encoded protein, Gabija anti-defence 1 (Gad1), directly binds to the Gabija GajAB complex and inactivates defence. A cryo-EM structure of the virally inhibited state shows that Gad1 forms an octameric web that encases the GajAB complex and inhibits DNA recognition and cleavage. Our results reveal the structural basis of assembly of the Gabija anti-phage defence complex and define a unique mechanism of viral immune evasion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

30S subunit recognition and G1405 modification by the aminoglycoside-resistance 16S ribosomal RNA methyltransferase RmtC

Acquired ribosomal RNA (rRNA) methylation has emerged as a significant mechanism of aminoglycoside resistance in pathogenic bacterial infections. Modification of a single nucleotide in the ribosome decoding center by the aminoglycoside-resistance 16S rRNA (m 7 G1405) methyltransferases effectively blocks the action of all 4,6-deoxystreptamine ring-containing aminoglycosides, including the latest generation of drugs. To define the molecular basis of 30S subunit recognition and G1405 modification by these enzymes, we used a S-adenosyl-L-methionine analog to trap the complex in a postcatalytic state to enable determination of a global 3.0 Å cryo-electron microscopy structure of the m 7 G1405 methyltransferase RmtC bound to the mature Escherichia coli 30S ribosomal subunit. This structure, together with functional analyses of RmtC variants, identifies the RmtC N-terminal domain as critical for recognition and docking of the enzyme on a conserved 16S rRNA tertiary surface adjacent to G1405 in 16S rRNA helix 44 (h44). To access the G1405 N7 position for modification, a collection of residues across one surface of RmtC, including a loop that undergoes a disorder-to order transition upon 30S subunit binding, induces significant distortion of h44. This distortion flips G1405 into the enzyme active site where it is positioned for modification by two almost universally conserved RmtC residues. Further, these studies expand our understanding of ribosome recognition by rRNA modification enzymes and present a more complete structural basis for future development of strategies to inhibit m 7 G1405 modification to resensitize bacterial pathogens to aminoglycosides.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Molecular chaperone RAP interacts with LRP1 in a dynamic bivalent mode and enhances folding of ligand-binding regions of other LDLR family receptors

The low-density lipoprotein receptor (LDLR) family of receptors are cell-surface receptors that internalize numerous ligands and play crucial role in various processes, such as lipoprotein metabolism, hemostasis, fetal development, etc. Previously, receptor-associated protein (RAP) was described as a molecular chaperone for LDLR-related protein 1 (LRP1), a prominent member of the LDLR family. We aimed to verify this role of RAP for LRP1 and two other LDLR family receptors, LDLR and vLDLR, and to investigate the mechanisms of respective interactions using a cell culture model system, purified system, and in silico modelling. Upon coexpression of RAP with clusters of the ligand-binding complement repeats (CRs) of the receptors in secreted form in insect cells culture, the isolated proteins had increased yield, enhanced folding, and improved binding properties compared with proteins expressed without RAP, as determined by circular dichroism and surface plasmon resonance. Within LRP1 CR-clusters II and IV, we identified multiple sites comprised of adjacent CR doublets, which provide alternative bivalent binding combinations with specific pairs of lysines on RAP. Mutational analysis of these lysines within each of isolated RAP D1/D2 and D3 domains having high affinity to LRP1 and of conserved tryptophans on selected CR-doublets of LRP1, as well as in silico docking of a model LRP1 CR-triplet with RAP, indicated a universal role for these residues in interaction of RAP and LRP1. Consequently, we propose a new model of RAP interaction with LDLR family receptors based on switching of the bivalent contacts between molecules over time in a dynamic mode.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Molecular basis and biological relevance of bacterial and plant pinoresinol/lariciresinol reductase specificities

A bacterial pinoresinol/lariciresinol reductase (PLR) homolog named NrPinZ was obtained from a Novosphingobium rhizosphaerae sp. LY bacterial strain, with NrPinZ being part of its 5-step biochemical system catabolizing pinoresinol into coniferyl aldehyde and vanillin. Recombinant NrPinZ reduces racemic 8–8′ furanofuran lignans [(±)-pinoresinols, medioresinols, and syringaresinols] with similar overall catalytic efficiencies. In those reductions, only one of the two furan ring systems is reduced. Two other bacterial PLR homologs, NaPinZ and SlPinZ, from N. aromaticivorans F199 and Sphingobium lignivorans SYK-6, respectively, had comparable substrate versatilities and catalytic efficacies. Plant PLR homologs, by comparison, are either enantiospecific, enantioselective, or variants thereof, being able to reduce either one or both furan rings. For example, a recombinant enantioselective PLR (PLR_Tp2) from western red cedar (Thuja plicata) preferentially reduces both (+)-pinoresinol furan rings to afford (−)-secoisolariciresinol. BoltZ-2 modeling of NrPinZ and PLR_Tp2, together with substrate docking of (+)- and (−)-pinoresinols, medioresinols, and syringaresinols, was very instructive. The NrPinZ active site P1/P2 sub-pockets allow for both racemic forms to be catabolized. Conversely, the smaller P1 pocket in PLR_Tp2 preferentially positions (+)-pinoresinol for downstream metabolism into (−)-secoisolariciresinol, thereby providing a biochemical explanation for the different stereochemical outcomes. NrPinZ, NaPinZ, and SlPinZ, catalyzing substrate versatile catabolism of both racemic forms, may have important ramifications for gymnosperm and angiosperm lignin and lignan biodegradation, including its evolutionary significance and potential in enzyme engineering.

Boltz-2 molecular modeling↗
DOE OSTI2023

Fast and versatile sequence-independent protein docking for nanomaterials design using RPXDock

Computationally designed multi-subunit assemblies have shown considerable promise for a variety of applications, including a new generation of potent vaccines. One of the major routes to such materials is rigid body sequence-independent docking of cyclic oligomers into architectures with point group or lattice symmetries. Current methods for docking and designing such assemblies are tailored to specific classes of symmetry and are difficult to modify for novel applications. Here we describe RPXDock, a fast, flexible, and modular software package for sequence-independent rigid-body protein docking across a wide range of symmetric architectures that is easily customizable for further development. RPXDock uses an efficient hierarchical search and a residue-pair transform ( RPX ) scoring method to rapidly search through multidimensional docking space. We describe the structure of the software, provide practical guidelines for its use, and describe the available functionalities including a variety of score functions and filtering tools that can be used to guide and refine docking results towards desired configurations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Hallucination of closed repeat proteins containing central pockets

In pseudocyclic proteins, such as TIM barrels, β barrels, and some helical transmembrane channels, a single subunit is repeated in a cyclic pattern, giving rise to a central cavity that can serve as a pocket for ligand binding or enzymatic activity. Inspired by these proteins, we devised a deep-learning-based approach to broadly exploring the space of closed repeat proteins starting from only a specification of the repeat number and length. Biophysical data for 38 structurally diverse pseudocyclic designs produced in Escherichia coli are consistent with the design models, and the three crystal structures we were able to obtain are very close to the designed structures. Docking studies suggest the diversity of folds and central pockets provide effective starting points for designing small-molecule binders and enzymes.

59 BASIC BIOLOGICAL SCIENCES↗