Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Protein Structure Prediction”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Novel sterol binding domains in bacteria

Sterol lipids are widely present in eukaryotes and play essential roles in signaling and modulating membrane fluidity. Although rare, some bacteria also produce sterols, but their function in bacteria is not known. Moreover, many more species, including pathogens and commensal microbes, acquire or modify sterols from eukaryotic hosts through poorly understood molecular mechanisms. The aerobic methanotroph Methylococcus capsulatus was the first bacterium shown to synthesize sterols, producing a mixture of C-4 methylated sterols that are distinct from those observed in eukaryotes. C-4 methylated sterols are synthesized in the cytosol and localized to the outer membrane, suggesting that a bacterial sterol transport machinery exists. Until now, the identity of such machinery remained a mystery. In this study, we identified three novel proteins that may be the first examples of transporters for bacterial sterol lipids. The proteins, which all belong to well-studied families of bacterial metabolite transporters, are predicted to reside in the inner membrane, periplasm, and outer membrane of M. capsulatus, and may work as a conduit to move modified sterols to the outer membrane. Quantitative analysis of ligand binding revealed their remarkable specificity for 4-methylsterols, and crystallographic structures coupled with docking and molecular dynamics simulations revealed the structural bases for substrate binding by two of the putative transporters. Their striking structural divergence from eukaryotic sterol transporters signals that they form a distinct sterol transport system within the bacterial domain. Finally, bioinformatics revealed the widespread presence of similar transporters in bacterial genomes, including in some pathogens that use host sterol lipids to construct their cell envelopes. The unique folds of these bacterial sterol binding proteins should now guide the discovery of other proteins that handle this essential metabolite.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Structure based virtual screening identifies small molecule effectors for the sialoglycan binding protein Hsa

Infective endocarditis (IE) is a cardiovascular disease often caused by bacteria of the viridans group of streptococci, which includes Streptococcus gordonii and Streptococcus sanguinis. Previous research has found that serine-rich repeat (SRR) proteins on the S. gordonii bacterial surface play a critical role in pathogenesis by facilitating bacterial attachment to sialylated glycans displayed on human platelets. Despite their important role in disease progression, there are currently no anti-adhesive drugs available on the market. Here, we performed structure-based virtual screening using an ensemble docking approach followed by consensus scoring to identify novel small molecule effectors against the sialoglycan binding domain of the SRR adhesin protein Hsa from the S. gordonii strain DL1. The screening successfully predicted nine compounds which were able to displace the native ligand (sialyl-T antigen) in an in vitro assay and bind competitively to Hsa. Furthermore, hierarchical clustering based on the MACCS fingerprints showed that eight of these small molecules do not share a common scaffold with the native ligand. This study indicates that SRR family of adhesin proteins can be inhibited by diverse small molecules and thus prevent the interaction of the protein with the sialoglycans. This opens new avenues for discovering potential drugs against IE.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Removal of Non-Specifically Bound Proteins Using Rayleigh Waves Generated on ST-Quartz Substrates

Label-free biosensors are plagued by the issue of non-specific protein binding which negatively affects sensing parameters such as sensitivity, selectivity, and limit-of-detection. In the current work, we explore the possibility of using the Rayleigh waves in ST-Quartz devices to efficiently remove non-specifically bound proteins via acoustic streaming. A coupled-field finite element (FE) fluid structure interaction (FSI) model of a surface acoustic wave (SAW) device based on ST-Quartz substrate in contact with a liquid loading was first used to predict trends in forces related to SAW-induced acoustic streaming. Based on model predictions, it is found that the computed SAW body force is sufficient to overcome adhesive forces between particles and a surface while lift and drag forces prevent reattachment for a range of SAW frequencies. We further performed experiments to validate the model predictions and observe that the excitation of Rayleigh SAWs removed non-specifically bound (NSB) antigens and antibodies from sensing and non-sensing regions, while rinsing and blocking agents were ineffective. An amplified RF signal applied to the device input disrupted the specific interactions between antigens and their capture antibody as well. ST-quartz allows propagation of Rayleigh and leaky SH-SAW waves in orthogonal directions. Thus, the results reported here could allow integration of three important biosensor functions on a single chip, i.e., removal of non-specific binding, mixing, and sensing in the liquid phase.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2021

Predicting the Functional State of Protein Kinases Using Interpretable Graph Neural Networks

Kinases are a family of proteins that function as molecular switches, regulating several essential cellular activities such as cell proliferation. Dysfunctional kinases are implicated in several types of cancers and hence they are actively pursued as drug targets. Given the vast number of complex kinase structures that are available in the protein data bank (PDB), there is a necessity to develop methodologies that can identify structurally important moieties of the kinases in an automated fashion, for such techniques can be instrumental in identifying novel drug targets. In this work, we develop a graph neural network (GNN) based deep learning framework for classifying the functionally active and inactive states of a large set of eukaryotic protein kinases, making use of their 3D structure from the PDB. We show that GNN based machine learning models can classify protein states with an accuracy greater than 97%. We further use the GNN models to automatically identify regions of the kinases that are important for its function. For this purpose, Gradient-weighted Class Activation Mapping (Grad-CAM) was implemented on the protein graphs. Remarkably, Grad-CAM consistently identifies the highly conserved DFG motif as the most important part of the protein across the entire kinome, without any prior input. Other regions of the hydrophobic core such as the HRD motif were also identified by the interpretable GNN framework, consistent with the literature. We discuss the significance of each of these regions in detail.

Ashwin Ravichandran↗
DOE OSTI2020

A Clustering-based biased Monte Carlo Approach to Protein Titration Curve Prediction

We develop and implement a novel approach to computing the ensemble averages in systems characterized by pair-wise interactions between the entities. Methods involving full enumeration of the configuration space result in exponential complexity. Sampling methods such as Markov Chain Monte Carlo (MCMC) algorithms have been proposed to tackle the exponential complexity of these problems. In certain scenarios where significant energetic coupling exists between the entities, the accuracy of the such algorithms can be diminished. We propose a strategy to improve the accuracy of the MCMC runs by taking advantage of the cluster structure in the interaction energy matrix. We propose two different schemes for performing the biased MCMC runs on the partitioned systems and show that they are valid MCMC schemes. We then apply these algorithms to the problem of computing the protonation fractions and hence the titration curves of titratable protein residues that constitute a given protein. We leverage both synthesized and real-world systems and show the improved performance of our biased MCMC methods when compared to the regular MCMC method.

Visweswara Sathanur, Arun↗
DOE OSTI2022

Design principles for site-selective hydroxylation by a Rieske oxygenase

Rieske oxygenases exploit the reactivity of iron to perform chemically challenging C–H bond functionalization reactions. Thus far, only a handful of Rieske oxygenases have been structurally characterized and remarkably little information exists regarding how these enzymes use a common architecture and set of metallocenters to facilitate a diverse range of reactions. Herein, we detail how two Rieske oxygenases SxtT and GxtA use different protein regions to influence the site-selectivity of their catalyzed monohydroxylation reactions. We present high resolution crystal structures of SxtT and GxtA with the native β-saxitoxinol and saxitoxin substrates bound in addition to a Xenon-pressurized structure of GxtA that reveals the location of a substrate access tunnel to the active site. Ultimately, this structural information allowed for the identification of six residues distributed between three regions of SxtT that together control the selectivity of the C–H hydroxylation event. Substitution of these residues produces a SxtT variant that is fully adapted to exhibit the non-native site-selectivity and substrate scope of GxtA. Importantly, we also found that these selectivity regions are conserved in other structurally characterized Rieske oxygenases, providing a framework for predictively repurposing and manipulating Rieske oxygenases as biocatalysts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Structure of HIV-1 RRE stem-loop II identifies two conformational states of the high-affinity Rev binding site

During HIV infection, specific RNA-protein interaction between the Rev response element (RRE) and viral Rev protein is required for nuclear export of intron-containing viral mRNA transcripts. Rev initially binds the high-affinity site in stem-loop II, which promotes oligomerization of additional Rev proteins on RRE. Here, we present the crystal structure of RRE stem-loop II in distinct closed and open conformations. The high-affinity Rev-binding site is located within the three-way junction rather than the predicted stem IIB. The closed and open conformers differ in their non-canonical interactions within the three-way junction, and only the open conformation has the widened major groove conducive to initial Rev interaction. Rev binding assays show that RRE stem-loop II has high- and low-affinity binding sites, each of which binds a Rev dimer. We propose a binding model, wherein Rev-binding sites on RRE are sequentially created through structural rearrangements induced by Rev-RRE interactions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Analysis of the molecular determinants for furin cleavage of the spike protein S1/S2 site in defined strains of the prototype coronavirus murine hepatitis virus (MHV)

We analyzed the spike protein S1/S2 cleavage of selected strains of a prototype coronavirus, mouse hepatitis virus (MHV) by the cellular protease furin, in order to understand the structural requirements underlying the sequence selectivity of the scissile segment. The probability of cleavage of selected MHV strains was first evaluated from furin cleavage scores predicted by the ProP computer software, and then cleavage was measured experimentally with a fluorogenic peptide cleavage assay consisting of S1/S2 peptide mimics and purified furin. We found that in vitro cleavability varied across MHV strains in line with predicted results—but with the notable exception of MHV-A59, which was not cleaved despite a high score predicted for its sequence. Using the known X-Ray structure of furin in complex with a substrate-like inhibitor as an initial structural reference, we carried out molecular dynamics (MD) simulations to learn the modes of binding of the peptides in the furin active site, and the suitability of the complex for initiation of the enzymatic cleavage. We identified the 3D structural requirements of the furin active site configuration that enable bound peptides to undergo cleavage, and the way in which the various strains tested experimentally are fulfilling these requirements. We find that despite some flexibility in the organization of the peptide bound to the active site of the enzyme, the presence of a histidine at P2 of MHV-A59 fails to properly orient the sidechain of His194 of the furin catalytic triad and therefore produces a distortion that renders the peptide/complex structural configuration in the active site incompatible with requirements for cleavage initiation. The Ser/Thr in P1 of MHV-2 and MHV-S has a similar effect of distorting the conformation of the furin active site residues produced by the elimination of the canonical salt-bridge formed by arginine in P1 position. This work informs a study of coronavirus infection and pathogenesis with respect to the function of the viral spike protein, and suggests an important process of viral adaptation and evolution within the spike S1/S2 structural loop.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Extending quantum-mechanical benchmark accuracy to biological ligand-pocket interactions

Predicting the binding affinity of ligands to protein pockets is key in the drug design pipeline. The flexibility of ligand-pocket motifs arises from a range of attractive and repulsive electronic interactions during binding. Accurately accounting for all interactions requires robust quantum-mechanical (QM) benchmarks, which are scarce for ligand-pocket systems. Additionally, disagreement between “gold standard” Coupled Cluster (CC) and Quantum Monte Carlo (QMC) methods casts doubt on many benchmarks for larger non-covalent systems. We introduce the “QUantum Interacting Dimer” (QUID) benchmark framework containing 170 non-covalent (non-)equilibrium systems modeling chemically and structurally diverse ligand-pocket motifs. Symmetry-adapted perturbation theory shows that QUID broadly covers non-covalent binding motifs and energetic contributions. Robust binding energies are obtained using complementary CC and QMC methods, achieving agreement of 0.5 kcal/mol. The benchmark data analysis reveals that several dispersion-inclusive density functional approximations provide accurate energy predictions, though their atomic van der Waals forces differ in magnitude and orientation. Contrarily, semiempirical methods and empirical force fields require improvements in capturing non-covalent interactions (NCIs) for out-of-equilibrium geometries. The wide span of NCIs, highly accurate interaction energies, and analysis of molecular properties take QUID beyond the “gold standard” for QM benchmarks of ligand-protein systems.

Puleva, Mirela [University of Luxembourg, Luxembou↗
DOE OSTI2024

Ultraviolet Superradiance from Mega-Networks of Tryptophan in Biological Architectures

Networks of tryptophan (Trp)–an aromatic amino acid with strong fluorescence response–are ubiquitous in biological systems, forming diverse architectures in transmembrane proteins, cytoskeletal filaments, subneuronal elements, photoreceptor complexes, virion capsids, and other cellular structures. We analyze the cooperative effects induced by ultraviolet (UV) excitation of several biologically relevant Trp mega-networks, thus giving insights into novel mechanisms for cellular signaling and control. Our theoretical analysis in the single-excitation manifold predicts the formation of strongly superradiant states due to collective interactions among organized arrangements of up to >10 5 Trp UV-excited transition dipoles in microtubule architectures, which leads to an enhancement of the fluorescence quantum yield (QY) that is confirmed by our experiments. We demonstrate the observed consequences of this superradiant behavior in the fluorescence QY for hierarchically organized tubulin structures, which increases in different geometric regimes at thermal equilibrium before saturation, highlighting the effect’s persistence in the presence of disorder. Our work thus showcases the many orders of magnitude across which the brightest (hundreds of femtoseconds) and darkest (tens of seconds) states can coexist in these Trp lattices.

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

Efficient production of oxidized terpenoids via engineering fusion proteins of terpene synthase and cytochrome P450

The functionalization of terpenes using cytochrome P450 enzymes is a versatile route to the production of useful derivatives that can be further converted to value-added products. Many terpenes are hydrophobic and volatile making their availability as a substrate for P450 enzymes significantly limited during microbial production. In this study, we developed a strategy to improve the accessibility of terpene molecules for the P450 reaction by linking terpene synthase and P450 together. As a model system, fusion proteins of 1,8-cineole synthase (CS) and P450 cin were investigated and it showed an improved hydroxylation of the monoterpenoid 1,8-cineole up to 5.4-fold. Structural analysis of the CS-P450 cin fusion proteins by SEC-SAXS indicated a dimer formation with preferred orientations of the active sites of the two domains. We also applied the enzyme fusion strategy to the oxidation of a sesquiterpene epi-isozizaene and the fusion enzymes significantly improved albaflavenol production in engineered E. coli. From the analysis of positive and negative examples of the fusion strategy, we proposed key factors in structure-based prediction and evaluation of fusion enzymes. Developing fusion enzymes for terpene synthase and P450 presents an efficient strategy toward oxidation of hydrophobic terpene compounds. This strategy could be widely applicable to improve the biosynthetic titer of the functionalized products from hydrophobic terpene intermediates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

In silico λ-dynamics predicts protein binding specificities to modified RNAs

Abstract RNA modifications shape gene expression through a variety of chemical changes to canonical RNA bases. Although numbering in the hundreds, only a few RNA modifications are well characterized, in part due to the absence of methods to identify modification sites. Antibodies remain a common tool to identify modified RNA and infer modification sites through straightforward applications. However, specificity issues can result in off-target binding and confound conclusions. This work utilizes in silico λ-dynamics to efficiently estimate binding free energy differences of modification-targeting antibodies between a variety of naturally occurring RNA modifications. Crystal structures of inosine and N6-methyladenosine (m6A) targeting antibodies bound to their modified ribonucleosides were determined and served as structural starting points. λ-Dynamics was utilized to predict RNA modifications that permit or inhibit binding to these antibodies. In vitro RNA-antibody binding assays supported the accuracy of these in silico results. High agreement between experimental and computed binding propensities demonstrated that λ-dynamics can serve as a predictive screen for antibody specificity against libraries of RNA modifications. More importantly, this strategy is an innovative way to elucidate how hundreds of known RNA modifications interact with biological molecules without the limitations imposed by in vitro or in vivo methodologies.

Biochemistry & Molecular Biology↗
DOE OSTI2026

Structural Heterogeneity and Hydrodynamics of an Intrinsically Disordered Protein Condensate

Biology demonstrates precise control over the free-energy landscape through the selective partitioning of biomacromolecules into membraneless organelles, enabling essential functions such as biochemical transformations, signaling cascades, and mechanical reinforcement. Although the function of these condensates depends on their underlying structure and hydrodynamics, molecular-scale information on these systems remains sparse. Here, in this study, neutron scattering is used to probe the organization and dynamics of the intrinsically disordered N-terminal domain of Galectin-3, an extracellular lectin responsible for facilitating liquid–liquid phase separation on the cellular surface, in both dilute and condensed phases. Dilute solutions contain isolated protein chains in equilibrium with mesoscopic clusters, whereas the condensed phase adopts a bicontinuous, microemulsion-like morphology. The dilute phase behavior is quantitatively described by coarse-grained polymer models from soft-matter physics, demonstrating their predictive power for complex biological proteins. At elevated concentrations, the proteins self-assemble akin to block copolymers, microphase separating through the aggregation of hydrophobic domains along the protein contour. The resulting condensate remains fluid-like despite a 25-fold increase in concentration; its internal hydrodynamics slow by only a factor of 3 relative to dilute protein chains. These results provide a molecular-level framework for how disordered proteins achieve both the structural complexity and dynamic fluidity of biomolecular condensates.

Carrick, Brian R. [Massachusetts Inst. of Technolo↗
DOE OSTI2024

Feeling the Strain: Quantifying Ligand Deformation in Photosynthesis

Structural distortion of protein-bound ligands can play a critical role in enzyme function by tuning the electronic and chemical properties of the ligand molecule. However, quantifying these effects is difficult due to the limited resolution of protein structures and the difficulty of generating accurate structural restrains for non-protein ligands. Here, we seek to quantify these effects through a statistical analysis of ligand distortion in Chlorophyll (Chl) proteins (CP), where ring deformation is thought to play a role in energy and electron transfer. To assess the accuracy of ring-deformation estimates from available structural data, we take advantage of the C 2 symmetry of Photosystem II (PSII), comparing ring-deformation estimates for equivalent sites both within and between 113 distinct X-ray and Cryogenic electron microscopy (CryoEM) PSII structures. Significantly, we find that several deformation modes exhibit considerable variability in predictions, even for equivalent monomers, down to 2 Å resolution, to an extent that probably prevents their utilization in optical calculations. We further find that refinement restrains play a critical role in determining deformation values to resolution as low as 2 Å. However, for those modes that are well-resolved in the structural data, ring deformation in PSII is strongly conserved across all species tested, from cyanobacteria to algae. Furthermore, these results highlight both the opportunities and limitations inherent in the structure-based analysis of the bioenergetic and optical properties of CPs and other protein-ligand complexes.

14 SOLAR ENERGY↗
DOE OSTI2024

Annual Report for Structure-Aware Unsupervised, Transformational Machine Learning for Drug Discovery

The major goal of this project is to develop machine learning (ML) methods to enable improved predictive power on real drug discovery for novel targets. More specifically, we plan to demonstrate the capability and effectiveness of ML tools utilizing unlabeled large-volume protein-ligand datasets. We also plan to demonstrate the capability and effectiveness of the developed methods by testing on a realistic drug discovery task to identify pan-coronavirus protease inhibitors such as SARS-CoV-2. While the overall goals and milestones remain consistent with the original proposal, certain technical details have been modified, which we will describe in this report.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

De Novo Design of Molecular Recognition for Sequence Defined Polymers

The proposed research aims at developing a computational framework that allows the de novo design of sequence defined polymers for molecular recognition. This framework will include the prediction of the 1) accessible backbone conformational space in solution environments for small to medium polymer length and 2) how side chains of different size and functionality can be used to design convergent or divergent receptor configurations. A structural database inspired by the protein database will be build that can be used to construct larger scaffolds based on smaller sequences

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

Evidence of the Involvement of a Cyclase Gene in the Biosynthesis of Ochratoxin A in Aspergillus carbonarius

Ochratoxin A (OTA) is a well-known mycotoxin with wide distribution in food and feed. Fungal genome sequencing has great utility for identifying secondary metabolites gene clusters for known and novel compounds. A comparative analysis of the OTA-biosynthetic cluster in A. steynii, A. westerdijkiae, A. niger, A. carbonarius, and P. nordicum has revealed a high synteny in OTA cluster organization in five structural genes (otaA, otaB, ota, otaR1, and otaD). Moreover, a recent detailed comparative genome analysis of Aspergilli OTA producers led to the identification of a cyclase gene, otaY, located in the OTA cluster between the otaA and otaB genes, encoding for a predicted protein with high similarity to SnoaLs domain. These proteins have been shown to catalyze ring closure steps in the biosynthesis of polyketide antibiotics produced in Streptomyces. In the present study, we demonstrated an upregulation of the cyclase gene in A. carbonarius under OTA permissive conditions, consistent with the expression trends of the other OTA cluster genes and their role in OTA biosynthesis by complete gene deletion. Our results pointed out the involvement of a cyclase gene in OTA biosynthetic pathway for the first time. They represent a step forward in the understanding of the molecular basis of OTA biosynthesis in A. carbonarius.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The role of filamentation in activation and DNA sequence specificity of the sequence-specific endonuclease SgrAI

Filament formation by metabolic, biosynthetic, and other enzymes has recently come into focus as a mechanism to fine-tune enzyme activity in the cell. Filamentation is key to the function of SgrAI, a sequence-specific DNA endonuclease that has served as a model system to provide some of the deepest insights into the biophysical characteristics of filamentation and its functional consequences. Structure-function analyses reveal that, in the filamentous state, SgrAI stabilizes an activated enzyme conformation that leads to accelerated DNA cleavage activity and expanded DNA sequence specificity. The latter is thought to be mediated by sequence-specific DNA structure, protein–DNA interactions, and a disorder-to-order transition in the protein, which collectively affect the relative stabilities of the inactive, non-filamentous conformation and the active, filamentous conformation of SgrAI bound to DNA. Full global kinetic modeling of the DNA cleavage pathway reveals a slow, rate-limiting, second-order association rate constant for filament assembly, and simulations of in vivo activity predict that filamentation is superior to non-filamenting mechanisms in ensuring rapid activation and sequestration of SgrAI's DNA cleavage activity on phage DNA and away from the host chromosome. In vivo studies demonstrate the critical requirement for accelerated DNA cleavage by SgrAI in its biological role to safeguard the bacterial host. Collectively, these data have advanced our understanding of how filamentation can regulate enzyme structure and function, while the experimental strategies used for SgrAI can be applied to other enzymatic systems to identify novel functional roles for filamentation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗