SEARCH · Engineering Papers
Results for “homology modeling”
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.
Structural Analysis of Diheme Cytochrome c by Hydrogen–Deuterium Exchange Mass Spectrometry and Homology Modeling
Explore the source record for details and available documents.
Reduction of a Heme Cofactor Initiates N -Nitroglycine Degradation by NnlA
Linear nitramines are potentially carcinogenic environmental contaminants. The NnlA enzyme from Variovorax sp. strain JS1663 degrades the nitramine N-nitroglycine (NNG)—a natural product produced by some bacteria—to glyoxylate and nitrite (NO 2 –). Ammonium (NH 4 +) was predicted as the third product of this reaction. A source of nonheme Fe II was shown to be required for initiation of NnlA activity. However, the role of this Fe II for NnlA activity was unclear. This work reveals that NnlA contains a b-type heme cofactor. Reduction of this heme—either by a nonheme iron source or dithionite—is required to initiate NnlA activity. Therefore, Fe II is not an essential substrate for holoenzyme activity. Our data show that reduced NnlA (Fe II -NnlA) catalyzes at least 100 turnovers and does not require O 2 . Finally, NH 4 + was verified as the third product, accounting for the complete nitrogen mass balance. Size exclusion chromatography showed that NnlA is a dimer in solution. Additionally, Fe II -NnlA is oxidized by O 2 and NO 2 – and stably binds carbon monoxide (CO) and nitric oxide (NO). These are characteristics shared with heme-binding PAS domains. Furthermore, a structural homology model of NnlA was generated using the PAS domain from Pseudomonas aeruginosa Aer2 as a template. The structural homology model suggested His 73 is the axial ligand of the NnlA heme. Site-directed mutagenesis of His 73 to alanine decreased the heme occupancy of NnlA and eliminated NNG activity, validating the homology model. We conclude that NnlA forms a homodimeric heme-binding PAS domain protein that requires reduction for initiation of the activity.
Centrifugal Pump Model for System Codes for Advanced NPP Designs
In this document, a homologous model and one-dimensional line model for centrifugal pump are described to provide the theoretical background of a pump model to be developed in system-level codes. Input and output parameters are proposed, along with a summary of the governing equations for each model. Furthermore, hydraulic performance degradation of centrifugal pump is also modeled in the homologous pump model. This document suggests that system codes provide a homologous and one-dimensional line model for centrifugal pumps, allowing the users to choose the model that best fits their purposes.
X-ray crystallographic and hydrogen deuterium exchange studies confirm alternate kinetic models for homolog insulin monomers
Despite the crucial role of various insulin analogs in achieving satisfactory glycemic control, a comprehensive understanding of their in-solution dynamic mechanisms still holds the potential to further optimize rapid insulin analogs, thus significantly improving the well-being of individuals with Type 1 Diabetes. Here, we employed hydrogen-deuterium exchange mass spectrometry to decipher the molecular dynamics of newly modified and functional insulin analog. A comparative analysis of H/D dynamics demonstrated that the modified insulin exchanges deuterium atoms faster and more extensively than the intact insulin aspart. Additionally, we present new insights derived from our 2.5 Å resolution X-ray crystal structure of modified hexamer insulin analog at ambient temperature. Furthermore, we obtained a distinctive side-chain conformation of the Asn3 residue on the B chain (AsnB3) by operating a comparative analysis with a previously available cryogenic rapid-acting insulin structure (PDB_ID: 4GBN). The experimental conclusions have demonstrated compatibility with modified insulin’s distinct cellular activity, comparably to aspart. Additionally, the hybrid structural approach combined with computational analysis employed in this study provides novel insight into the structural dynamics of newly modified and functional insulin vs insulin aspart monomeric entities. It allows further molecular understanding of intermolecular interrelations driving dissociation kinetics and, therefore, a fast action mechanism.
Modeling the Influenza A NP-vRNA-Polymerase Complex in Atomic Detail
Seasonal flu is an acute respiratory disease that exacts a massive toll on human populations, healthcare systems and economies. The disease is caused by an enveloped Influenza virus containing eight ribonucleoprotein (RNP) complexes. Each RNP incorporates multiple copies of nucleoprotein (NP), a fragment of the viral genome (vRNA), and a viral RNA-dependent RNA polymerase (POL), and is responsible for packaging the viral genome and performing critical functions including replication and transcription. A complete model of an Influenza RNP in atomic detail can elucidate the structural basis for viral genome functions, and identify potential targets for viral therapeutics. In this work we construct a model of a complete Influenza A RNP complex in atomic detail using multiple sources of structural and sequence information and a series of homology-modeling techniques, including a motif-matching fragment assembly method. Our final model provides a rationale for experimentally-observed changes to viral polymerase activity in numerous mutational assays. Further, our model reveals specific interactions between the three primary structural components of the RNP, including potential targets for blocking POL-binding to the NP-vRNA complex. The methods developed in this work open the possibility of elucidating other functionally-relevant atomic-scale interactions in additional RNP structures and other biomolecular complexes.
Structural basis for a highly conserved RNA-mediated enteroviral genome replication
Abstract Enteroviruses contain conserved RNA structures at the extreme 5′ end of their genomes that recruit essential proteins 3CD and PCBP2 to promote genome replication. However, the high-resolution structures and mechanisms of these replication-linked RNAs (REPLRs) are limited. Here, we determined the crystal structures of the coxsackievirus B3 and rhinoviruses B14 and C15 REPLRs at 1.54, 2.2 and 2.54 Å resolution, revealing a highly conserved H-type four-way junction fold with co-axially stacked sA-sD and sB-sC helices that are stabilized by a long-range A•C•U base-triple. Such conserved features observed in the crystal structures also allowed us to predict the models of several other enteroviral REPLRs using homology modeling, which generated models almost identical to the experimentally determined structures. Moreover, our structure-guided binding studies with recombinantly purified full-length human PCBP2 showed that two previously proposed binding sites, the sB-loop and 3′ spacer, reside proximally and bind a single PCBP2. Additionally, the DNA oligos complementary to the 3′ spacer, the high-affinity PCBP2 binding site, abrogated its interactions with enteroviral REPLRs, suggesting the critical roles of this single-stranded region in recruiting PCBP2 for enteroviral genome replication and illuminating the promising prospects of developing therapeutics against enteroviral infections targeting this replication platform.
Graph-Directed Approach for Downselecting Toxins for Experimental Structure Determination
Conotoxins are short, cysteine-rich peptides of great interest as novel therapeutic leads and of great concern as lethal biological agents due to their high affinity and specificity for various receptors involved in neuromuscular transmission. Currently, of the approximately 6000 known conotoxin sequences, only about 3% have associated structural characterization, which leads to a bottleneck in rapid high-throughput screening (HTS) for identification of potential leads or threats. In this work, we combine a graph-based approach with homology modeling to expand the library of conotoxin structures and to identify those conotoxin sequences that are of the greatest value for experimental structural characterization. The latter would allow for the rapid expansion of the known structural space for generating high quality template-based models. Our approach generalizes to other evolutionarily-related, short, cysteine-rich venoms of interest. Overall, we present and validate an approach for venom structure modeling and experimental guidance and employ it to produce a 290%-larger library of approximate conotoxin structures for HTS. We also provide a set of ranked conotoxin sequences for experimental structure determination to further expand this library.
ppdx : Automated modeling of protein–protein interaction descriptors for use with machine learning
This paper describes ppdx, a python workflow tool that combines protein sequence alignment, homology modeling, and structural refinement, to compute a broad array of descriptors for characterizing protein–protein interactions. The descriptors can be used to predict various properties of interest, such as protein–protein binding affinities, or inhibitory concentrations (IC 50 ), using approaches that range from simple regression to more complex machine learning models. The software is highly modular. It supports different protocols for generating structures, and 95 descriptors can be currently computed. More protocols and descriptors can be easily added. The implementation is highly parallel and can fully exploit the available cores in a single workstation, or multiple nodes on a supercomputer, allowing many systems to be analyzed simultaneously. As an illustrative application, ppdx is used to parametrize a model that predicts the IC 50 of a set of antigens and a class of antibodies directed to the influenza hemagglutinin stalk.
Structure of the RECK CC domain, an evolutionary anomaly
Explore the source record for details and available documents.
Structural characterization and computational analysis of PDZ domains in Monosiga brevicollis
Abstract Identification of the molecular networks that facilitated the evolution of multicellular animals from their unicellular ancestors is a fundamental problem in evolutionary cellular biology. Choanoflagellates are recognized as the closest extant nonmetazoan ancestors to animals. These unicellular eukaryotes can adopt a multicellular‐like “rosette” state. Therefore, they are compelling models for the study of early multicellularity. Comparative studies revealed that a number of putative human orthologs are present in choanoflagellate genomes, suggesting that a subset of these genes were necessary for the emergence of multicellularity. However, previous work is largely based on sequence alignments alone, which does not confirm structural nor functional similarity. Here, we focus on the PDZ domain, a peptide‐binding domain which plays critical roles in myriad cellular signaling networks and which underwent a gene family expansion in metazoan lineages. Using a customized sequence similarity search algorithm, we identified 178 PDZ domains in the Monosiga brevicollis proteome. This includes 11 previously unidentified sequences, which we analyzed using Rosetta and homology modeling. To assess conservation of protein structure, we solved high‐resolution crystal structures of representative M. brevicollis PDZ domains that are homologous to human Dlg1 PDZ2, Dlg1 PDZ3, GIPC, and SHANK1 PDZ domains. To assess functional conservation, we calculated binding affinities for mbGIPC, mbSHANK1, mbSNX27, and mbDLG‐3 PDZ domains from M. brevicollis . Overall, we find that peptide selectivity is generally conserved between these two disparate organisms, with one possible exception, mbDLG‐3. Overall, our results provide novel insight into signaling pathways in a choanoflagellate model of primitive multicellularity.
Mercury methylation by metabolically versatile and cosmopolitan marine bacteria
Microbes transform aqueous mercury (Hg) into methylmercury (MeHg), a potent neurotoxin that accumulates in terrestrial and marine food webs, with potential impacts on human health. This process requires the gene pair hgcAB, which encodes for proteins that actuate Hg methylation, and has been well described for anoxic environments. However, recent studies report potential MeHg formation in suboxic seawater, although the microorganisms involved remain poorly understood. In this study, we conducted large-scale multi-omic analyses to search for putative microbial Hg methylators along defined redox gradients in Saanich Inlet, British Columbia, a model natural ecosystem with previously measured Hg and MeHg concentration profiles. Analysis of gene expression profiles along the redoxcline identified several putative Hg methylating microbial groups, including Calditrichaeota, SAR324 and Marinimicrobia, with the last the most active based on hgc transcription levels. Marinimicrobia hgc genes were identified from multiple publicly available marine metagenomes, consistent with a potential key role in marine Hg methylation. Computational homology modelling predicts that Marinimicrobia HgcAB proteins contain the highly conserved amino acid sites and folding structures required for functional Hg methylation. Furthermore, a number of terminal oxidases from aerobic respiratory chains were associated with several putative novel Hg methylators. Our findings thus reveal potential novel marine Hg-methylating microorganisms with a greater oxygen tolerance and broader habitat range than previously recognized.
Short-wavelength-sensitive 2 (Sws2) visual photopigment models combined with atomistic molecular simulations to predict spectral peaks of absorbance
For many species, vision is one of the most important sensory modalities for mediating essential tasks that include navigation, predation and foraging, predator avoidance, and numerous social behaviors. The vertebrate visual process begins when photons of the light interact with rod and cone photoreceptors that are present in the neural retina. Vertebrate visual photopigments are housed within these photoreceptor cells and are sensitive to a wide range of wavelengths that peak within the light spectrum, the latter of which is a function of the type of chromophore used and how it interacts with specific amino acid residues found within the opsin protein sequence. Minor differences in the amino acid sequences of the opsins are known to lead to large differences in the spectral peak of absorbance (i.e. the λ max value). In our prior studies, we developed a new approach that combined homology modeling and molecular dynamics simulations to gather structural information associated with chromophore conformation, then used it to generate statistical models for the accurate prediction of λ max values for photopigments derived from Rh1 and Rh2 amino acid sequences. In the present study, we test our novel approach to predict the λ max of phylogenetically distant Sws2 cone opsins. To build a model that can predict the λ max using our approach presented in our prior studies, we selected a spectrally-diverse set of 11 teleost Sws2 photopigments for which both amino acid sequence information and experimentally measured λ max values are known. The final first-order regression model, consisting of three terms associated with chromophore conformation, was sufficient to predict the λ max of Sws2 photopigments with high accuracy. This study further highlights the breadth of our approach in reliably predicting λ max values of Sws2 cone photopigments, evolutionary-more distant from template bovine RH1, and provided mechanistic insights into the role of known spectral tuning sites.
Mechanisms of Polyethylene Terephthalate Pellet Fragmentation into Nanoplastics and Assimilable Carbons by Wastewater Comamonas
Comamonadaceae bacteria are enriched on poly(ethylene terephthalate) (PET) microplastics in wastewaters and urban rivers, but the PET-degrading mechanisms remain unclear. Here, we investigated these mechanisms with Comamonas testosteroniKF-1, a wastewater isolate, by combining microscopy, spectroscopy, proteomics, protein modeling, and genetic engineering. Compared to minor dents on PET films, scanning electron microscopy revealed significant fragmentation of PET pellets, resulting in a 3.5-fold increase in the abundance of small nanoparticles (<100 nm) during 30-day cultivation. Infrared spectroscopy captured primarily hydrolytic cleavage in the fragmented pellet particles. Solution analysis further demonstrated double hydrolysis of a PET oligomer, bis(2-hydroxyethyl) terephthalate, to the bioavailable monomer terephthalate. Supplementation with acetate, a common wastewater co-substrate, promoted cell growth and PET fragmentation. Of the multiple hydrolases encoded in the genome, intracellular proteomics detected only one, which was found in both acetate-only and PET-only conditions. Homology modeling of this hydrolase structure illustrated substrate binding analogous to reported PET hydrolases, despite dissimilar sequences. Mutants lacking this hydrolase gene were incapable of PET oligomer hydrolysis and had a 21% decrease in PET fragmentation; re-insertion of the gene restored both functions. Thus, we have identified constitutive production of a key PET-degrading hydrolase in wastewater Comamonas, which could be exploited for plastic bioconversion.
On the Rapid Calculation of Binding Affinities for Antigen and Antibody Design and Affinity Maturation Simulations
The accurate and efficient calculation of protein-protein binding affinities is an essential component in antibody and antigen design and optimization, and in computer modeling of antibody affinity maturation. Such calculations remain challenging despite advances in computer hardware and algorithms, primarily because proteins are flexible molecules, and thus, require explicit or implicit incorporation of multiple conformational states into the computational procedure. The astronomical size of the amino acid sequence space further compounds the challenge by requiring predictions to be computed within a short time so that many sequence variants can be tested. In this study, we compare three classes of methods for antibody/antigen (Ab/Ag) binding affinity calculations: (i) a method that relies on the physical separation of the Ab/Ag complex in equilibrium molecular dynamics (MD) simulations, (ii) a collection of 18 scoring functions that act on an ensemble of structures created using homology modeling software, and (iii) methods based on the molecular mechanics-generalized Born surface area (MM-GBSA) energy decomposition, in which the individual contributions of the energy terms are scaled to optimize agreement with the experiment. When applied to a set of 49 antibody mutations in two Ab/HIV gp120 complexes, all of the methods are found to have modest accuracy, with the highest Pearson correlations reaching about 0.6. In particular, the most computationally intensive method, i.e., MD simulation, did not outperform several scoring functions. The optimized energy decomposition methods provided marginally higher accuracy, but at the expense of requiring experimental data for parametrization. Within each method class, we examined the effect of the number of independent computational replicates, i.e., modeled structures or reinitialized MD simulations, on the prediction accuracy. We suggest using about ten modeled structures for scoring methods, and about five simulation replicates for MD simulations as a rule of thumb for obtaining reasonable convergence. We anticipate that our study will be a useful resource for practitioners working to incorporate binding affinity calculations within their protein design and optimization process.
Crystal structures of Salmonella enterica FraB deglycase reveal a conformational heterodimer with remarkable structural plasticity at the active site
Abstract Thefralocus ofSalmonella entericaencodes five genes for metabolism of fructose‐asparagine, an Amadori product formed by condensation of asparagine with glucose. In the last step of this pathway, the FraB deglycase cleaves 6‐phospho‐fructose‐aspartate into glucose‐6‐phosphate and aspartate. In homology models, FraB forms a homodimer with two equivalent active sites located at the dimer interface. E214 and H230, two invariant residues essential for catalysis, project into each active site cleft from opposing subunits of the dimer. Here, we have determined six crystal structures of FraB, three of a variant containing an N‐terminal His 6 tag and two mutations needed for crystallization (hereafter referred to as WT′), two with additional mutations to active site residues (E214A and P232A), and one of a variant with C‐terminal residues 313–325 deleted. Surprisingly, in the WT′ FraB structure, the two catalytic residues, E214 (general base) and H230 (general acid), are positioned ~22 Å apart. In the E214A and C‐terminus‐truncated FraB variants, however, a conformational change in the E214‐residing helix brings E214 and H230* to ~7 Å (* indicates residue from the second protomer that creates the inter‐subunit catalytic center). The loop bearing H230 also exhibits significant variation, ranging from being completely disordered to adopting open or closed states, with the nearby P232* residue being eithercisortrans. The C‐terminal residues 313–325 form a flexible “C‐tail” that can be fully disordered, bind in the active site to block access of substrate, or angle across the active site to wrap across the other subunit of the dimer and potentially close over substrate. Collectively, these structures reveal that FraB is a conformational heterodimer with two chemically identical subunits that are constrained to adopt different structures as they come together for catalysis. This plasticity likely involves correlated opening and closure of the two active sites for their respective binding and release of substrates and ligands.
Molecular dynamics analysis of the binding of human interleukin-6 with interleukin-6 α-receptor
Human interleukin-6 (hIL-6) is a multifunctional cytokine that regulates immune and inflammatory responses in addition to metabolic and regenerative processes and cancer. hIL-6 binding to the IL-6 receptor (IL-6Rα) induces homodimerization and recruitment of the glycoprotein (gp130) to form a hexameric signaling complex. Anti-IL-6 and IL-6R antibodies are clinically approved inhibitors of IL-6 signaling pathway for treating rheumatoid arthritis and Castleman's disease, respectively. There is a potential to develop novel small molecule IL-6 antagonists derived from understanding the structural basis for IL-6/IL-6Rα interactions. In this work, we combine homology modeling with extensive molecular dynamics (MD) simulations to examine the association of hIL-6 with IL-6Rα. A comparison with MD of apo hIL-6 reveals that the binding of hIL-6 to IL-6Rα induces structural and dynamic rearrangements in the AB loop region of hIL-6, disrupting intraprotein contacts and increasing the flexibility of residues 48 to 58 of the AB loop. In contrast, due to the involvement of residues 59 to 78 in forming contacts with the receptor, these residues of the AB loop are observed to rigidify in the presence of the receptor. The binary complex is primarily stabilized by two pairs of salt bridges, Arg181 (hIL-6)- Glu182 (IL-6Rα) and Arg184 (hIL-6)- Glu183 (IL-6Rα) as well as hydrophobic and aromatic stacking interactions mediated essentially by Phe residues in both proteins. An interplay of electrostatic, hydrophobic, hydrogen bonding, and aromatic stacking interactions facilitates the formation of the hIL-6/IL-6Rα complex.
Enhancement of anticancer drug sensitivity in multidrug resistance cells overexpressing ATP-binding cassette (ABC) transporter ABCC10 by CP55, a synthetic derivative of 5-cyano-6-phenylpyrimidin
Highlights: • Multidrug resistance (MDR) in cancer is one of the major obstacles in chemotherapy. • ABCC10 mediates resistance to wide spectrum of chemotherapeutic agents, but only limited number of studies are available. • Overexpression of ABCC10 has been linked to various types of cancers including colorectal, gastric and lung cancer. • CP55 is a potent ABCC10 inhibitor which enhances the anticancer drug sensitivity in ABCC10-overexpression MDR cells. ATP-binding cassette (ABC) transporter C10 (ABCC10), also named multidrug resistance protein 7 (MRP7), is a member of ABC transporter superfamily and has been revealed to transport a wide range of chemotherapeutic agents including taxanes, epothilone B, Vinca alkaloids, and anthracyclines. In our previous study, a 5-cyano-6-phenylpyrimidin derivative CP55 was synthesized and found significantly reversal effect of multidrug resistance (MDR) mediated by ABCB1. In this study, we found CP55 also efficiently reversed MDR mediated by ABCC10. Our in vitro study showed that co-treatment with CP55 significantly increased the efficacy of ABCC10-substrate anticancer drugs in MDR cells overexpressing ABCC10. Furthermore, we showed that treatment with CP55 increased the intracellular accumulation of [{sup 3}H]-labeled anticancer drugs and in-turn decreasing drug efflux by inhibiting the transport activity, without altering ABCC10 protein ex-pression level or cellular localization. Potential CP55-ABCC10 interactions were predicted via docking analysis using human ABCC10 homology model and obtained high docking score. Therefore, CP55 represents a promising therapeutic agent in the combinational treatment of chemo-resistant cancer related to ABCC10.