SEARCH · Engineering Papers
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.
Protein-DNA Interaction Workflow (PDI Workflow), a pipeline that focuses on generating high-quality docking and molecular dynamics simulations for Protein-DNA complexes. This allows us to take DNA sequences with unknown tertiary structures, accurately predict their structure, dock them with the desired target protein, and then simulate their interactions using molecular dynamics simulations
Abstract The delicate interplay of covalent and non‐covalent interactions in proteins is inherently quantum mechanical and highly dynamic in nature. To directly interrogate the evolving nature of the electronic structure of proteins, we carry out 100‐ps‐scale ab initio molecular dynamics simulations of three representative small proteins with range‐separated hybrid density functional theory. We quantify the nature and length‐scale of the coupling of residue‐specific charge probability distributions in these proteins. While some nonpolar residues exhibit expectedly narrow charge distributions, most polar and charged residues exhibit broad, multimodal distributions. Even for nonpolar residues, we observe sequence‐specific deviations corresponding to charge accumulation or depletion that would be challenging to capture in a fixed charge force field. We quantify the effect of residue‐residue interactions on charge distributions first with linear cross‐correlations. We then show how additional insight can be gained from evaluating the mutual information of charge distributions. We show that a significant number of residues couple most strongly with residues that are distant in both sequence and space over a range of secondary structures including α‐helical, β‐sheet, disulfide bridging, and lasso motifs. The mutual information analysis is necessary to capture coupling between some polar and charged residues that would be otherwise missed.
Interactions between the major structural glycoprotein E2 of classical swine fever virus (CSFV) with host proteins have been identified as important factors affecting virus replication and virulence. Previously, using the yeast two-hybrid system, we identified swine host proteins specifically interacting with CSFV E2. In this report, we use a proximity ligation assay to demonstrate that swine host protein CCDC115 interacts with E2 in CSFV-infected swine cells. Using a randomly mutated E2 library in the context of a yeast two-hybrid methodology, specific amino acid mutations in the CSFV E2 protein responsible for disrupting the interaction with CCDC115 were identified. A recombinant CSFV mutant (E2ΔCCDC115v) harboring amino acid changes disrupting the E2 protein interaction with CCDC115 was produced and used as a tool to assess the role of the E2–CCDC115 interaction in viral replication and virulence in swine. CSFV E2ΔCCDC115v showed a slightly decreased ability to replicate in the SK6 swine cell line and a greater replication defect in primary swine macrophage cultures. A decreased E2–CCDC115 interaction detected by PLA is observed in cells infected with E2ΔCCDC115v. Importantly, animals intranasally infected with 105 TCID50 of E2ΔCCDC115v experienced a significantly longer survival period when compared with those infected with the parental Brescia strain. This result would indicate that the ability of CSFV E2 to bind host CCDC115 protein during infection plays an important role in virus replication in swine macrophages and in virus virulence during the infection in domestic swine.
SparcleQC is a Python package that, given a protein:ligand complex in the Protein Data Bank (PDB) file format, can create quantum mechanics/molecular mechanics (QM/MM)-like input files for the electronic structure theory packages PSI4, QChem, and NWChem. The resulting input files include quantum mechanical representations of the ligand and a small section of the protein, surrounded by point charges that represent the rest of the protein. Creation of these QM/MM input files includes cutting and capping the QM subregion, obtaining point charges for the protein, and adjusting charges at the QM/MM boundary; and each of these tasks are automated by the software. In this article, we describe the details of SparcleQC’s procedure, show examples of the Python API, and explain additional features that are helpful in protein:ligand interaction studies. Finally, we show that SparcleQC enables automated preparation of input files for QM/MM calculations, which can return can return accurate interaction energies in minutes, while a fully quantum mechanical computation on the protein:ligand complex could take days, if it is even possible.
Gram-positive bacteria use type VII secretion systems (T7SSs) to export effector proteins that manipulate the physiology of nearby prokaryotic and eukaryotic cells. Several mycobacterial T7SSs have established roles in virulence. By contrast, the genetically distinct T7SSb pathway found in Firmicutes bacteria more often functions to mediate bacterial competition. A lack of structural information on the T7SSb has limited the understanding of effector export by this protein secretion apparatus. Here, we present the 2.4 angstrom crystal structure of the extracellular region of the T7SSb subunit EsaA from Streptococcus gallolyticus. Our structure reveals that homodimeric EsaA is an elongated, arrow-shaped protein with a surface-accessible "tip", which in some species of bacteria serves as a receptor for lytic bacteriophages. Because it is the only T7SSb subunit large enough to traverse the peptidoglycan layer of Firmicutes, we propose that EsaA plays a critical role in transporting effectors across the entirety of the Gram-positive cell envelope.
Introduction VPS45 belongs to the Sec1/Munc18 family of proteins, which interact with and regulate Qa-SNARE function during membrane fusion. We have shown previously that Arabidopsis thaliana VPS45 interacts with the SYP61/SYP41/VTI12 SNARE complex, which locates on the trans -Golgi network (TGN). It is required for SYP41 stability, and it functions in cargo trafficking to the vacuole and in cell expansion. It is also required for correct auxin distribution during gravitropism and lateral root growth. Results As vps45 knockout mutation is lethal in Arabidopsis, we identified a mutant, vps45-3 , with a point mutation in the VPS45 gene causing a serine 284-to-phenylalanine substitution. The VPS45-3 protein is stable and maintains interaction with SYP61 and SYP41. However, vps45-3 plants display severe growth defects with significantly reduced organ and cell size, similar to vps45 RNAi transgenic lines that have reduced VPS45 protein levels. Root hair and pollen tube elongation, both processes of tip growth, are highly compromised in vps45-3 . Mutant root hairs are shorter and thicker than those of wild-type plants, and are wavy. These root hairs have vacuolar defects, containing many small vacuoles, compared with WT root hairs with a single large vacuole occupying much of the cell volume. Pollen tubes were also significantly shorter in vps45-3 compared to WT. Discussion We thus show that VPS45 is essential for proper tip growth and propose that the observed vacuolar defects lead to loss of the turgor pressure needed for tip growth.
Photoreceptor phosphodiesterase PDE6 is central for visual signal transduction. Maturation of PDE6 depends on a specialized chaperone complex of HSP90 with aryl hydrocarbon receptor-interacting protein-like 1 (AIPL1). Disruption of PDE6 maturation underlies a severe form of retina degeneration. Here, we report a 3.9 Å cryoelectron microscopy (cryo-EM) structure of the complex of HSP90 with AIPL1. This structure reveals a unique interaction of the FK506-binding protein (FKBP)-like domain of AIPL1 with HSP90 at its dimer interface. Unusually, the N terminus AIPL1 inserts into the HSP90 lumen in a manner that was observed previously for HSP90 clients. Deletion of the 7 N-terminal residues of AIPL1 decreased its ability to cochaperone PDE6. Multi-body refinement of the cryo-EM data indicated large swing-like movements of AIPL1-FKBP. Modeling the complex of HSP90 with AIPL1 using crosslinking constraints indicated proximity of the mobile tetratricopeptide repeat (TPR) domain with the C-terminal domain of HSP90. Our study establishes a framework for future structural studies of PDE6 maturation.
Adverse environmental conditions reduce crop productivity and often increase the load of unfolded or misfolded proteins in the endoplasmic reticulum (ER). This potentially lethal condition, known as ER stress, is buffered by the unfolded protein response (UPR), a set of signaling pathways designed to either recover ER functionality or ignite programmed cell death. Despite the biological significance of the UPR to the life of the organism, the regulatory transcriptional landscape underpinning ER stress management is largely unmapped, especially in crops. To fill this significant knowledge gap, we performed a large–scale systems–level analysis of the protein–DNA interaction (PDI) network in maize (Zea mays). Using 23 promoter fragments of six UPR marker genes in a high–throughput enhanced yeast one–hybrid assay, we identified a highly interconnected network of 262 transcription factors (TFs) associated with significant biological traits and 831 PDIs underlying the UPR. We established a temporal hierarchy of TF binding to gene promoters within the same family as well as across different families of TFs. Cistrome analysis revealed the dynamic activities of a variety of cis–regulatory elements (CREs) in ER stress–responsive gene promoters. By integrating the cistrome results into a TF network analysis, we mapped a subnetwork of TFs associated with a CRE that may contribute to UPR management. Lastly, we validated the role of a predicted network hub gene using the Arabidopsis system. The PDIs, TF networks, and CREs identified in our work are foundational resources for understanding transcription–regulatory mechanisms in the stress responses and crop improvement.
Prokaryotes and viruses have fought a long battle against each other. Prokaryotes use CRISPR–Cas-mediated adaptive immunity, while conversely, viruses evolve multiple anti-CRISPR (Acr) proteins to defeat these CRISPR–Cas systems. The type I-F CRISPR–Cas system in Pseudomonas aeruginosa requires the crRNA-guided surveillance complex (Csy complex) to recognize the invading DNA. Although some Acr proteins against the Csy complex have been reported, other relevant Acr proteins still need studies to understand their mechanisms. As such, here, we obtain three structures of previously unresolved Acr proteins (AcrF9, AcrF8, and AcrF6) bound to the Csy complex using electron cryo-microscopy (cryo-EM), with resolution at 2.57 Å, 3.42 Å, and 3.15 Å, respectively. The 2.57-Å structure reveals fine details for each molecular component within the Csy complex as well as the direct and water-mediated interactions between proteins and CRISPR RNA (crRNA). Our structures also show unambiguously how these Acr proteins bind differently to the Csy complex. AcrF9 binds to key DNA-binding sites on the Csy spiral backbone. AcrF6 binds at the junction between Cas7.6f and Cas8f, which is critical for DNA duplex splitting. AcrF8 binds to a distinct position on the Csy spiral backbone and forms interactions with crRNA, which has not been seen in other Acr proteins against the Csy complex. Our structure-guided mutagenesis and biochemistry experiments further support the anti-CRISPR mechanisms of these Acr proteins. Our findings support the convergent consequence of inhibiting degradation of invading DNA by these Acr proteins, albeit with different modes of interactions with the type I-F CRISPR–Cas system.
Researchers often measure changes in gene expression across conditions to better understand the shared functional roles and regulatory mechanisms of different genes. Analogous to this is comparing gene expression across species, which can improve our understanding of the evolutionary processes shaping the evolution of both individual genes and functional pathways. One area of interest is determining genes showing signals of coevolution, which can also indicate potential functional similarity, analogous to co-expression analysis often performed across conditions for a single species. However, as with any trait, comparing gene expression across species can be confounded by the non-independence of species due to shared ancestry, making standard hypothesis testing inappropriate. We compared RNA-Seq data across 18 fungal species using a multivariate Brownian Motion phylogenetic comparative method (PCM), which allowed us to quantify coevolution between protein pairs while directly accounting for the shared ancestry of the species. Our work indicates proteins which physically-interact show stronger signals of coevolution than randomly-generated pairs. Interactions with stronger empirical and computational evidence also showing stronger signals of coevolution. We examined the effects of number of protein interactions and gene expression levels on coevolution, finding both factors are overall poor predictors of the strength of coevolution between a protein pair. Simulations further demonstrate the potential issues of analyzing gene expression coevolution without accounting for shared ancestry in a standard hypothesis testing framework. Furthermore, our simulations indicate the use of a randomly-generated null distribution as a means of determining statistical significance for detecting coevolving genes with phylogenetically-uncorrected correlations, as has previously been done, is less accurate than PCMs, although is a significant improvement over standard hypothesis testing. These methods are further improved by using a phylogenetically-corrected correlation metric. Our work highlights potential benefits of using PCMs to detect gene expression coevolution from high-throughput omics scale data. This framework can be built upon to investigate other evolutionary hypotheses, such as changes in transcription regulatory mechanisms across species.
The serotonin transporter (SERT) is a member of the SLC6 neurotransmitter transporter family that mediates serotonin reuptake at presynaptic nerve terminals. SERT is the target of both therapeutic antidepressant drugs and psychostimulant substances such as cocaine and methamphetamines, which are small molecules that perturb normal serotonergic transmission by interfering with serotonin transport. Despite decades of studies, important functional aspects of SERT such as the oligomerization state of native SERT and its interactions with potential proteins remain unresolved. Here, we develop methods to isolate SERT from porcine brain (pSERT) using a mild, nonionic detergent, utilize fluorescence-detection size-exclusion chromatography to investigate its oligomerization state and interactions with other proteins, and employ single-particle cryo-electron microscopy to elucidate the structures of pSERT in complexes with methamphetamine or cocaine, providing structural insights into psychostimulant recognition and accompanying pSERT conformations. Methamphetamine and cocaine both bind to the central site, stabilizing the transporter in an outward open conformation. We also identify densities attributable to multiple cholesterol or cholesteryl hemisuccinate (CHS) molecules, as well as to a detergent molecule bound to the pSERT allosteric site. Under our conditions of isolation, we find that pSERT is best described as a monomeric entity, isolated without interacting proteins, and is ensconced by multiple cholesterol or CHS molecules.
Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.
Ovate family proteins (OFPs) are a family of plant growth regulators that play diverse roles in many aspects of physiological processes. OFPs have been characterized in various plant species including tomato, Arabidopsis , and rice. However, little is known about OFPs in woody species. Here, a total of 30 PtOFP genes were identified from the genome of Populus trichocarpa and were further grouped into four subfamilies based on their sequence similarities. Gene expression analysis indicated that some members of the PtOFP gene family displayed tissue/organ-specific patterns. Analysis of cis -acting elements in the promoter as well as gene expression by hormone treatment revealed putative involvement of PtOFPs in hormonal response. Furthermore, PtOFP1 (Potri.006G107700) was further experimentally demonstrated to act as a transcriptional repressor. Yeast two-hybrid assay showed physical interactions of PtOFP1 with other proteins, which suggests that they might function in various cellular processes by forming protein complexes. In addition, overexpression of PtOFP1 in Arabidopsis conferred enhanced tolerance to PEG-induced drought stress at seedling stage, as well as a higher survival rate than the wild type at mature stage. These results provide a systematic analysis of the Populus OFP gene family and lay a foundation for functional characterization of this gene family.
Highlights: • CCHCR1 and HAX1 proteins directly interact. • CCHCR1 is a psoriasis candidate gene located on PSORS1 locus. • HAX1 gene is located on PSORS4 locus and is overexpressed in psoriasis. • CCHCR1 and HAX1 co-localize with and affect specific mRNAs in vesicles. • Novel protein/mRNA networks involving CCHCR1 and HAX1 could play a role in psoriasis. CCHCR1 (Coiled-Coil alpha-Helical Rod 1), maps to chromosomal region 6p21.3, within the major psoriasis susceptibility locus PSORS1. CCHCR1 itself is a plausible psoriasis candidate gene, however its role in psoriasis pathogenesis remains unclear. We previously demonstrated that CCHCR1 protein acts as a cytoplasmic docking site for RNA polymerase II core subunit 3 (RPB3) in cycling cells, suggesting a role for CCHCR1 in vesicular trafficking between cellular compartments. Here, we report a novel interaction between CCHCR1 and the RNA binding protein HAX1. HAX1 maps to chromosomal region 1q21.3 within the PSORS4 locus and is over-expressed in psoriasis. Both CCHCR1 and HAX1 share subcellular co-localization with mitochondria, nuclei and cytoplasmic vesicles as P-bodies. By a series of ribonucleoprotein immunoprecipitation (RIP) assays, we isolated a pool of mRNAs complexed with HAX1 and/or CCHCR1 proteins. Among the mRNAs complexed with both CCHCR1 and HAX1 proteins, there are Vimentin mRNA, previously described to be bound by HAX1, and CAMP/LL37 mRNA, whose gene product is over-expressed in psoriasis.
Three-dimensional structure determination and prediction of proteins with intrinsically disordered regions, unstructured regions, conformational flexibility, and lacking homologous structures are challenging. We previously predicted and refined an in silico structure of a plant cellulose synthase from cotton (GhCESA1), and more recently, cryo-electron microscopy (cryo-EM) has resolved a majority of the lengths of two CESA structures from poplar (PttCESA8) and cotton (GhCESA7). However, 26–30% of these cryo-EM structures remain unresolved, including the N-terminal domain, half of the class-specific region, the gating loop region, and the C-terminal domain. Here, we describe the generation and evaluation of a full-length hybrid GhCESA1 model based on this cryo-EM PttCESA8 structure, with unresolved regions completed using this in silico refined GhCESA1 model. All-atom molecular dynamics simulations and subsequent energy minimizations were performed for the in silico and hybrid GhCESA1 models in a lipid bilayer-water-ion environment, and structural stability, dynamics, energetics, contacts, and quality were evaluated. The unresolved regions were found to be the most dynamic, in agreement with their poor electron density with cryo-EM. The hybrid model exhibited a higher total secondary structure content, more favorable intra-protein and protein-lipid interaction energies, and improved quality metrics. Moreover, hydrogen bonding was revealed to be a primary mechanism for intra-protein and protein-lipid contacts. These results demonstrate that in silico structure prediction and refinement may be useful to augment experimental structure determination, especially for disordered and unstructured regions. Furthermore, this hybrid model can serve as a steppingstone to derive full-length homology models of other CESAs found in more experimentally tractable organisms.
For billions of years, nature has optimized the photosynthetic machinery that converts light energy into chemical energy. Key primary reactions of photosynthesis occur in large membrane protein-cofactor complexes. The light-induced sequential electron transfer reactions occur through a chain of donor/acceptor cofactors embedded in the protein matrix resulting in a long-lived transmembrane charge-separated state. EPR is the method of choice to study electron transfer and the interaction of protein environment with redox-active cofactors. However, the spectra of organic cofactor radicals typically are not fully resolved and severely overlap at conventional X-band EPR. Even at Q-band EPR, this overlap is present and often a serious problem. As a result, there is a large variation of the reported EPR data and limited understanding of electronic structures of several redox-active cofactors. These serious problems can often be overcome by the excellent spectral resolution provided by high-frequency EPR (HF EPR). In this work, we study the electronic structure of the primary electron donor P 700 and the secondary electron acceptor A 1 of Photosystem I (PSI) using 130 GHz (D-band) EPR and Electron-Nuclear-Double-Resonance (ENDOR) spectroscopy. PSI was isotopically labeled with 15 N (I = $\frac{1/2}$) to avoid quadrupolar interactions in the most abundant nitrogen isotope 14 N (I = 1) and simplify the ENDOR spectra. ENDOR spectroscopy is central for determining the hyperfine coupling of nitrogen atoms of the two chlorophyll molecules comprising oxidized P 700 and the involvement of protein nitrogen atoms with reduced A 1 . While HF ENDOR of A 1 - allows identification of two nitrogen atoms, HF ENDOR of P 700 + still does not permit unique assignment of the recorded hyperfine couplings.
Renal interstitial fibrosis (RIF) is the common irreversible pathway by which chronic kidney disease (CKD) progresses to the end stage. The transforming growth factor-β (TGF-β)/signal transducer and activator of transcription 3 (STAT3) signaling pathway is a common factor leading to inflammation-mediated RIF, but its downstream regulatory mechanism is still unclear. Bioinformatics analysis predicted that serum amyloid A protein 1 (SAA1) was one of the target genes for transcriptional activation of STAT3 signaling. As an acute phase reaction protein, SAA1 plays an important role in many inflammatory reactions, and research has suggested that SAA1 is significantly elevated in the serum of patients with CKD. In this research, multiple experiments were performed to investigate the role of SAA1 in the process of RIF. SAA1 was abnormally highly expressed in kidney tissue from individuals who underwent unilateral ureteral obstruction (UUO) and TGF-β-induced HK2 cells, and the abnormal expression was directly related to the transcriptional activation of STAT3. Additionally, SAA1 can directly target and bind valosin-containing protein (VCP)-interacting membrane selenoprotein (VIMP) to inhibit the function of the Derlin-1/VCP/VIMP complex, preventing the transportation and degradation of the misfolded protein, resulting in endoplasmic reticulum (ER) stress characterized by an increase in glucose-regulated protein 78 (GRP78) levels and ultimately promoting the occurrence and development of RIF.
Reducing protein adhesion is a critical strategy in fouling-resistant material innovation, with broad applications spanning biomedical and healthcare devices, biosensors, industrial and environmental systems, and other important technological domains. Here, in this study, we elucidated protein adhesion behavior on polystyrene-based thin films by neutron reflectometry (NR) and quartz crystal microbalance with dissipation (QCM-D), using both lysozyme and bovine serum albumin (BSA) as model proteins. To this end, semifluorinated polystyrene thin films with gradient wettability and surface energy were fabricated through dry processing using plasma oxidation and gas-phase deposition. Although it is believed that a fully fluorinated alkyl chain offers extremely low surface energy, thus rejecting foulants, and has been used in many fouling-resistant surface designs, enhanced protein–surface interactions were observed consistently in NR and QCM-D results, due to the combined effects of surface morphology and chemistry. On the contrary, depositing shorter fluorinated silane onto a hydrophilic PS surface contributed to a more homogeneous nanoscale fluorine coating, resulting in less initial protein adsorption and improved surface recovery. Comparative analysis of proteins with different sizes on the nanopatterned semifluorinated surface revealed the influence of molecular characteristics on surface interactions. Lysozyme, being smaller and more compact, showed faster adsorption kinetics and higher surface coverage but largely reversible binding, whereas BSA, with its larger and more flexible structure, formed broader and more stable interfacial layers. This study fills the gap in understanding protein adhesion within the range of hydrophobicity (water contact angle ∼90°), as current strategies often associate with extreme hydrophilic and superhydrophobic surfaces due to hydration or low-surface-energy rejection mechanisms, respectively. It also provides in-depth insights into current combinatorial fouling-resistant surface design.