Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein folding”

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 541 records · Page 30

Ceruloplasmin copper induces oxidant damage by a redox process utilizing cell-derived superoxide as reductant

Oxidative damage by transition metals bound to proteins may be an important pathogenic mechanism. Ceruloplasmin (Cp) is a Cu-containing plasma protein thought to be involved in oxidative modification of lipoproteins. We have previously shown that Cp increased cell-mediated low-density lipoprotein (LDL) oxidation by a process requiring cell-derived superoxide, but the underlying chemical mechanism(s) is (are) unknown. We now show that superoxide reduction of Cp Cu is a critical reaction in cellular LDL oxidation. By bathocuproine disulfonate (BCS) binding and by superoxide utilization, we showed that exogenous superoxide reduces a single Cp Cu atom, the same Cu required for LDL oxidation. The Cu atom remained bound to Cp during the redox cycle. Three avenues of evidence showed that vascular cells reduce Cp Cu by a superoxide-dependent process. The 2-fold higher rate of Cp Cu reduction by smooth muscle cells (SMC) compared to endothelial cells (EC) was consistent with their relative rates of superoxide release. Furthermore, Cp Cu reduction by cells was blocked by Cu,Zn superoxide dismutase (SOD1). Finally, the level of superoxide produced by EC and SMC was sufficient to cause the amount of Cu reduction observed. An important role of Cp Cu reduction in LDL oxidation was suggested by results showing that SOD1 inhibited Cp Cu reduction and LDL oxidation by SMC with equal potency, while tumor necrosis factor-alpha stimulated both processes. In summary, these results show that superoxide is a critical cellular reductant of divalent transition metals involved in oxidation, and that protein-bound Cu is a substrate for this reaction. The role of these mechanisms in oxidative processes in vivo has yet to be defined.

NASA Discipline Regulatory Physiology↗

Nitric oxide delivery and heme-assisted S-nitrosation by the bedbug nitrophorin

Here, nitrophorins are heme proteins used by blood feeding insects to deliver nitric oxide (NO) to a victim, leading to vasodilation and antiplatelet activity. Cimex lectularius (bedbug) nitrophorin (cNP) accomplishes this with a cysteine ligated ferric (Fe(III)) heme. In the acidic environment of the insect's salivary glands, NO binds tightly to cNP. During a blood meal, cNP-NO is delivered to the feeding site where dilution and increased pH lead to NO release. In a previous study, cNP was shown to not only bind heme, but to also nitrosate the proximal cysteine, leading to Cys-NO (SNO) formation. SNO formation requires oxidation of the proximal cysteine, which was proposed to be metal-assisted through accompanying reduction of ferric heme and formation of Fe(II)-NO. Here, we report the 1.6 Å crystal structure of cNP first chemically reduced and then exposed to NO, and show that Fe(II)-NO is formed but SNO is not, supporting a metal-assisted SNO formation mechanism. Crystallographic and spectroscopic studies of mutated cNP show that steric crowding of the proximal site inhibits SNO formation while a sterically relaxed proximal site enhances SNO formation, providing insight into specificity for this poorly understood modification. Experiments examining the pH dependence for NO implicate direct protonation of the proximal cysteine as the underlying mechanism. At lower pH, thiol heme ligation predominates, leading to a smaller trans effect and 60-fold enhanced NO affinity (K d = 70 nM). Unexpectedly, we find that thiol formation interferes with SNO formation, suggesting cNP-SNO is unlikely to form in the insect salivary glands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Light‐Driven Metabolic Pathways in Non‐Photosynthetic Biohybrid Bacteria

Biomanufacturing via microorganisms relies on carbon substrates for molecular feedstocks and a source of energy to carry out enzymatic reactions. This creates metabolic bottlenecks and lowers the efficiency for substrate conversion. Nanoparticle biohybridization with proteins and whole cell surfaces can bypass the need for redox cofactor regeneration for improved secondary metabolite production in a non-specific manner. Here we propose using nanobiohybrid organisms (Nanorgs), intracellular protein-nanoparticle hybrids formed through the spontaneous coupling of core-shell quantum dots (QDs) with histidine-tagged enzymes in non-photosynthetic bacteria, for light-mediated control of bacterial metabolism. Furthermore, this proved to eliminate metabolic constrictions and replace glucose with light as the source of energy in Escherichia coli, with an increase in growth by 1.7-fold in 75 % reduced nutrient media. Metabolomic tracking through carbon isotope labeling confirmed flux shunting through targeted pathways, with accumulation of metabolites downstream of respective targets. Finally, application of Nanorgs with the Ehrlich pathway improved isobutanol titers/yield by 3.9-fold in 75 % less sugar from E. coli strains with no genetic alterations. These results demonstrate the promise of Nanorgs for metabolic engineering and low-cost biomanufacturing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mining metatranscriptomes reveals a vast world of viroid-like circular RNAs

Viroids and viroid-like covalently closed circular (ccc) RNAs are minimal replicators that typically encode no proteins and hijack cellular enzymes for replication. The extent and diversity of viroid-like agents are poorly understood. We developed a computational pipeline to identify viroid-like cccRNAs and applied it to 5,131 metatranscriptomes and 1,344 plant transcriptomes. The search yielded 11,378 viroid-like cccRNAs spanning 4,409 species-level clusters, a 5-fold increase compared to the previously identified viroid-like elements. Within this diverse collection, we discovered numerous putative viroids, satellite RNAs, retrozymes, and ribozy-like viruses. Diverse ribozyme combinations and unusual ribozymes within the cccRNAs were identified. Self-cleaving ribozymes were identified in ambiviruses, some mito-like viruses and capsid-encoding satellite virus-like cccRNAs. The broad presence of viroid-like cccRNAs in diverse transcriptomes and ecosystems implies that their host range is far broader than currently known, and matches to CRISPR spacers suggest that some cccRNAs replicate in prokaryotes.

59 BASIC BIOLOGICAL SCIENCES↗

A Chimeric LBT-GFP Biosensor Exhibits Antithetical Fluorescence Responses to Ca 2+ and Dy 3+ Binding

Rare earth elements (REEs) are critical components in emerging technologies, but their mining and refining processes are often laborious, costly, and environmentally damaging. Developing green and efficient separation methods for REEs is crucial. Biomolecular approaches using lanthanide-binding proteins and peptides show promise for selective REE extraction and separation. In this study, we present the design and characterization of a genetically encoded fluorescence indicator (GEFI) construct that combines a superfolder green fluorescent protein (sfGFP) with a dual lanthanide-binding tag (2×dLBT). The 2×dLBT insert induces conformational changes in sfGFP upon lanthanide binding, modulating the fluorescence intensity. The sfGFP-2×dLBT biosensor exhibited distinct fluorescence responses to different lanthanide ions, with the highest dynamic range observed for heavy REEs like dysprosium (Dy 3+ ). Interestingly, the sensor displayed an antithetical response, where low concentrations of lanthanides initially quenched the fluorescence, but higher concentrations led to a significant fluorescence increase (1.5-fold). The Ca 2+ ion on the other hand showed only a dose-dependent quenching of the fluorescence response. Based on these observations, the biphasic response of the biosensor to lanthanides was eliminated by pretreating the sensor with calcium, which further expanded the dynamic range up to 3-fold for Dy 3+ . The lanthanide-selective and concentration-dependent fluorescence changes of the sfGFP-2×dLBT biosensor demonstrate its potential as a platform for developing specific sensors for various REEs. These sensors could enable rapid and cost-effective determination of REE composition in complex mixtures, facilitating the separation and recovery of critical REEs from electronic waste and other REE-containing sources.

59 BASIC BIOLOGICAL SCIENCES↗

Dueling Backbones: Comparing Peptoid and Peptide Analogues of a Mussel Adhesive Protein

Ensembles of amino acid side chains often dominate the interfacial interactions of intrinsically disordered proteins; however, backbone contributions are far from negligible. Using a combination of nanoscale force measurements and molecular dynamics simulations, we demonstrated with analogous mussel-mimetic adhesive peptides and peptoids both 34 residues long that highly divergent adhesive/cohesive outcomes can be achieved on mica surfaces by altering backbone chemistry only. The Phe, Tyr, and Dopa containing peptoid variants used in this study deposited as dehydrated and incompressible films that facilitated analysis of peptoid side chain contributions to adhesion and cohesion. For example, whereas Phe and Dopa peptoids exhibited similar cohesion, Dopa peptoids were ~3 times more adhesive than Phe peptoids on mica. Compared with the peptides, Phe peptoid achieved only ~20% of Phe containing peptide adhesion, but the Dopa peptoids were >2-fold more adhesive than the Dopa peptides. Cation-π interactions accounted for some but not all of the cohesive interactions. Finally, our results were corroborated by molecular dynamics simulations and highlight the importance of backbone chemistry and the potential of peptoids or peptoid/peptide hybrids as wet adhesives and primers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Small molecule biosensors

Biosensors for small molecules can be used in applications that range from metabolic engineering to orthogonal control of transcription. Biosensors are produced based on a ligand-binding domain (LBD) using a method that, in principle, can be applied for any target molecule. The LBD is fused to either a fluorescent protein or a transcriptional activator and is destabilized by mutation such that the fusion accumulates only in cells containing the target ligand. The power of this method is illustrated by developing biosensors for digoxin and progesterone. Addition of ligand to cells expressing a biosensor activates transcription in yeast, mammalian cells and plants, with a dynamic range of up to about 100-fold or more. The biosensors are used to improve the biotransformation of pregnenolone to progesterone in yeast and to regulate CRISPR activity in mammalian cells. This work provides a general methodology to develop biosensors for a broad range of molecules.

Mandell, Daniel J.↗

Computational Prediction of Coiled–Coil Protein Gelation Dynamics and Structure

Protein hydrogels represent an important and growing biomaterial for a multitude of applications, including diagnostics and drug delivery. We have previously explored the ability to engineer the thermoresponsive supramolecular assembly of coiled–coil proteins into hydrogels with varying gelation properties, where we have defined important parameters in the coiled–coil hydrogel design. Using Rosetta energy scores and Poisson–Boltzmann electrostatic energies, we iterate a computational design strategy to predict the gelation of coiled–coil proteins while simultaneously exploring five new coiled–coil protein hydrogel sequences. Provided this library, we explore the impact of in silico energies on structure and gelation kinetics, where we also reveal a range of blue autofluorescence that enables hydrogel disassembly and recovery. As a result of this library, we identify the new coiled–coil hydrogel sequence, Q5, capable of gelation within 24 h at 4 °C, a more than 2-fold increase over that of our previous iteration Q2. The fast gelation time of Q5 enables the assessment of structural transition in real time using small-angle X-ray scattering (SAXS) that is correlated to coarse-grained and atomistic molecular dynamics simulations revealing the supramolecular assembling behavior of coiled–coils toward nanofiber assembly and gelation. This work represents the first system of hydrogels with predictable self-assembly, autofluorescent capability, and a molecular model of coiled–coil fiber formation.

36 MATERIALS SCIENCE↗

Hydrodynamic Effects on Biofilm Development and Recombinant Protein Expression

Hydrodynamics play an important role in the rate of cell attachment and nutrient and oxygen transfer, which can affect biofilm development and the level of recombinant protein production. In the present study, the effects of different flow conditions on the development of Escherichia coli biofilms and the expression of a model recombinant protein (enhanced green fluorescent protein, eGFP) were examined. Planktonic and biofilm cells were grown at two different flow rates in a recirculating flow cell system for 7 days: 255 and 128 L h−1 (corresponding to a Reynolds number of 4600 and 2300, respectively). The fluorometric analysis showed that the specific eGFP production was higher in biofilms than in planktonic cells under both hydrodynamic conditions (3-fold higher for 255 L h−1 and 2-fold higher for 128 L h−1). In the biofilm cells, the percentage of eGFP-expressing cells was on average 52% higher at a flow rate of 255 L h−1. Furthermore, a higher plasmid copy number (PCN) was obtained for the highest flow rate for both planktonic (244 PCN/cell versus 118 PCN/cell) and biofilm cells (43 PCN/cell versus 29 PCN/cell). The results suggested that higher flow velocities promoted eGFP expression in E. coli biofilms.

Soares, Alexandra↗

Differential expression of members of the annexin multigene family in Arabidopsis

Although in most plant species no more than two annexin genes have been reported to date, seven annexin homologs have been identified in Arabidopsis, Annexin Arabidopsis 1-7 (AnnAt1--AnnAt7). This establishes that annexins can be a diverse, multigene protein family in a single plant species. Here we compare and analyze these seven annexin gene sequences and present the in situ RNA localization patterns of two of these genes, AnnAt1 and AnnAt2, during different stages of Arabidopsis development. Sequence analysis of AnnAt1--AnnAt7 reveals that they contain the characteristic four structural repeats including the more highly conserved 17-amino acid endonexin fold region found in vertebrate annexins. Alignment comparisons show that there are differences within the repeat regions that may have functional importance. To assess the relative level of expression in various tissues, reverse transcription-PCR was carried out using gene-specific primers for each of the Arabidopsis annexin genes. In addition, northern blot analysis using gene-specific probes indicates differences in AnnAt1 and AnnAt2 expression levels in different tissues. AnnAt1 is expressed in all tissues examined and is most abundant in stems, whereas AnnAt2 is expressed mainly in root tissue and to a lesser extent in stems and flowers. In situ RNA localization demonstrates that these two annexin genes display developmentally regulated tissue-specific and cell-specific expression patterns. These patterns are both distinct and overlapping. The developmental expression patterns for both annexins provide further support for the hypothesis that annexins are involved in the Golgi-mediated secretion of polysaccharides.

Non-NASA Center↗

Prostanoid-induced expression of matrix metalloproteinase-1 messenger ribonucleic acid in rat osteosarcoma cells

Individual prostanoids have distinct potencies in activating intracellular signaling pathways and regulating gene expression in osteoblastic cells. The E-series prostaglandins (PGs) are known to stimulate matrix metalloproteinase-1 (MMP-1) synthesis and secretion in certain rodent and human osteoblastic cells, yet the intracellular events involved remain unclear. To further characterize this response and its signal transduction pathway(s), we examined prostanoid-induced expression of the MMP-1 gene in the rat osteoblastic osteosarcoma cell line UMR 106-01. Northern blot analysis demonstrated that prostaglandin E2 (PGE2) and PGE1 were very potent stimulators (40-fold) of MMP-1 transcript abundance, PGF2 alpha and prostacyclin were weak stimulators (4-fold), and thromboxane-B2 had no effect. The marked increase in MMP-1 transcript abundance after PGE2 treatment was first detected at 2 h, became maximal at 4 h, and persisted beyond 24 h. This response was dose dependent and elicited maximal and half-maximal effects with concentrations of 10(-6) and 0.6 x 10(-7) M, respectively. Cycloheximide, a protein synthesis inhibitor, completely blocked this effect of PGE2, suggesting that the expression of other genes is required. Nuclear run-on experiments demonstrated that PGE2 rapidly activates MMP-1 gene transcription, with a maximal increase at 2-4 h. The second messenger analog, 8-bromo-cAMP, mimicked the effects of PGE2 by stimulating a dose-dependent increase in MMP-1 messenger RNA (mRNA) levels, with a maximal effect quantitatively similar to that observed with PGE2. Thus, in UMR 106-01 cells, different prostanoids have distinct potencies in stimulating MMP-1 mRNA abundance. Our data suggest that PGE2 stimulation of MMP-1 synthesis is due to activation of MMP-1 gene transcription and a subsequent marked increase in MMP-1 mRNA abundance. This effect is dependent on de novo protein synthesis and is mimicked by protein kinase-A activation.

NASA Discipline Musculoskeletal↗

Structural, biophysical, and biochemical insights into C–S bond cleavage by dimethylsulfone monooxygenase

Sulfur is an essential element for life. Bacteria can obtain sulfur from inorganic sulfate; but in the sulfur starvation–induced response,Pseudomonadsemploy two-component flavin-dependent monooxygenases (TC-FMOs) from themsuandsfnoperons to assimilate sulfur from environmental compounds including alkanesulfonates and dialkylsulfones. Here, we report binding studies of oxidized FMN to enzymes involved within theP. fluorescensenzymatic pathway responsible for converting dimethylsulfone (DMSO 2 ) to sulfite. In this catabolic pathway, SfnG serves as the initial TC-FMO for sulfur assimilation, which is investigated in detail by solving the 2.6-Å resolution crystal structure of unliganded SfnG and the 1.75-Å resolution crystal structure of the SfnG ternary complex containing FMN and DMSO 2 . We find that SfnG adopts a (β/α) 8 barrel fold with a distinct quaternary configuration from other tetrameric class C TC-FMOs. To probe the unexpected tetramer arrangement, structural heterogeneity is assessed by chromatography and light scattering to confirm ligand binding correlates with a tetramer. Binding of FMN and DMSO 2 accompanies ordering of the active site, with DMSO 2 bound on thesi-face of the flavin. A previously unobserved protein backbone conformation is found within the oxygen-binding site on there-face of the flavin. Functional assays and the positioning of ligands with respect to the oxygen-binding site are consistent with use of an N5-(hydro)peroxyflavin pathway. Biochemical endpoint assays and docking studies reveal SfnG breaks the C–S bond of a range of dialkylsulfones.

Science & Technology - Other Topics↗

Rare ribosomal RNA sequences from archaea stabilize the bacterial ribosome

The ribosome serves as the universally conserved translator of the genetic code into proteins and supports life across diverse temperatures ranging from below freezing to above 120°C. Ribosomes are capable of functioning across this wide range of temperatures even though the catalytic site for peptide bond formation, the peptidyl transferase center, is nearly universally conserved. Here we find that Thermoproteota, a phylum of thermophilic Archaea, substitute cytidine for uridine at large subunit rRNA positions 2554 and 2555 (Escherichia coli numbering) in the A loop, immediately adjacent to the binding site for the 3'-end of A-site tRNA. We show by cryo-EM that E. coli ribosomes with uridine to cytidine mutations at these positions retain the proper fold and post-transcriptional modification of the A loop. Additionally, these mutations do not affect cellular growth, protect the large ribosomal subunit from thermal denaturation, and increase the mutational robustness of nucleotides in the peptidyl transferase center. This work identifies sequence variation across archaeal ribosomes in the peptidyl transferase center that likely confers stabilization of the ribosome at high temperatures and develops a stable mutant bacterial ribosome that can act as a scaffold for future ribosome engineering efforts.

59 BASIC BIOLOGICAL SCIENCES↗

MARTINI-Compatible Coarse-Grained Model for the Mesoscale Simulation of Peptoids

Peptoids (poly-N-substituted glycines) are a class of synthetic polymers that are regioisomers of peptides (poly-C-substituted glycines), in which the point of side-chain connectivity is shifted from the backbone C to the N atom. Peptoids have found diverse applications as peptidomimetic drugs, protein mimetic polymers, surfactants, and catalysts. Computational modeling is valuable in the understanding and design of peptoid-based nanomaterials. In this work, we report the bottom-up parameterization of coarse-grained peptoid force fields based on the MARTINI peptide force field against all-atom peptoid simulation data. Our parameterization pipeline iteratively refits coarse-grained bonded interactions using iterative Boltzmann inversion and nonbonded interactions by matching the potential of mean force for chain extension. We assure good sampling of the amide bond cis/trans isomerizations in the all-atom simulation data using parallel bias metadynamics. We develop coarse-grained models for two representative peptoids—polysarcosine (poly(N-methyl glycine)) and poly(N-((4-bromophenyl)ethyl)glycine)—and show their structural and thermodynamic properties to be in excellent accord with all-atom calculations but up to 25-fold more efficient and compatible with MARTINI force fields. Here, this work establishes a new rigorously parameterized coarse-grained peptoid force field for the understanding and design of peptoid nanomaterials at length and time scales inaccessible to all-atom calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The structural basis for 2′−5′/3′−5′-cGAMP synthesis by cGAS

Abstract cGAS activates innate immune responses against cytosolic double-stranded DNA. Here, by determining crystal structures of cGAS at various reaction stages, we report a unifying catalytic mechanism. apo-cGAS assumes an array of inactive conformations and binds NTPs nonproductively. Dimerization-coupled double-stranded DNA-binding then affixes the active site into a rigid lock for productive metal•substrate binding. A web-like network of protein•NTP, intra-NTP, and inter-NTP interactions ensures the stepwise synthesis of 2′−5′/3′−5′-linked cGAMP while discriminating against noncognate NTPs and off-pathway intermediates. One divalent metal is sufficient for productive substrate binding, and capturing the second divalent metal is tightly coupled to nucleotide and linkage specificities, a process which manganese is preferred over magnesium by 100-fold. Additionally, we elucidate how mouse cGAS achieves more stringent NTP and linkage specificities than human cGAS. Together, our results reveal that an adaptable, yet precise lock-and-key-like mechanism underpins cGAS catalysis.

59 BASIC BIOLOGICAL SCIENCES↗

Rational Design of Lanmodulin Variants for Size-Based Selectivity of Individual Rare Earth Elements

Rare earth elements (REEs) are essential to modern technologies, yet their high physical and chemical similarity makes separation of individual REEs difficult and environmentally taxing. Metalloproteins offer a promising alternative for selective REE binding, as they tend to have high metal ion affinity and specificity. Lanmodulin (LanM), in particular, has arisen as a potential candidate for REE separation as it exhibits picomolar affinity for elements in the REE family. Prior work has shown that the single point mutation D9N can shift LanM’s preference away from lanthanides toward actinides, motivating efforts to tune selectivity of LanM through targeted mutagenesis. Here, we tested the hypothesis that introducing selective aspartic acid to glutamic acid substitutions in the metal coordinating EF hands of LanM would impose steric constraints that would drive LanM affinity away from larger ions, such as La3+, to smaller ions, such as Y3+. To test this hypothesis, a combination of computational and experimental approaches were employed to evaluate the signal mutations LanM D5E and LanM D3E and the double mutants LanM D1ED5E and LanM D3ED9E. Surprisingly, increasing the number of mutations within the metal center did not enhance affinity for smaller REEs, or decrease affinity for larger ions. Only the single point mutation LanM D5E weakened La3+ binding by one order of magnitude relative to LanM wild type (WT), and pairing it with a second mutation to produce LanM D1ED5E drove La3+ affinity to be stronger than that seen for LanM WT. The D3E mutation alone prevented proper expression and folding, but paring it with D9E to produce LanM D3ED9E rescued expression and yielded La3+ affinities comparable to LanM WT. All variants that expressed (LanM D5E, LanM D1ED5E, LanM D3ED9E) displayed Y3+ affinities comparable to LanM WT. Overall, these results highlight the tunability of LanM’s metal-binding environment but also expose current limitations in predicting structural responses to point mutations within a protein sequence. This work establishes a foundation that can be used for refining computational and experimental strategies to engineer metalloproteins with tailored REE selectivity.

Close, Emily [Pacific Northwest National Laborator↗

Effect of Increased Cyclic AMP Concentration on Muscle Protein Synthesis and Beta-Adrenergic Receptor Expression in Chicken Skeletal Muscle Cells in Culture

Analogies of epinephrine are known to cause hypertrophy of skeletal muscle when fed to animals. These compounds presumably exert their physiological action through interaction with the P-adrenergic receptor. Since the intracellular signal generated by the Beta-adrenergic receptor is cyclic AMP (cAMP), experiments were initiated in cell culture to determine if artificial elevation of cAMP by treatment with forskolin would alter muscle protein metabolism and P-adrenergic receptor expression. Chicken skeletal muscle cells after 7 days in culture were treated with 0.2-30 micrometers forskolin for a total of three days. At the end of the treatment period, both the concentration of cAMP and the quantity of myosin heavy chain (MHC) were measured. Concentration of cAMP in forskolin-treated cells increased up to 10-fold in a dose dependent manner. In contrast, the quantity of MHC was increased approximately 50% above control cells at 0.2 micrometers forskolin, but exhibited a gradual decline at higher levels of forskolin so that the quantity of MHC in cells treated with 30 micrometers forskolin was not significantly different from controls. Curiously, the intracellular concentration of cAMP which elicited the maximum increase in the quantity of MHC was only 40% higher than cAMP concentration in control cells.

Young, R. B.↗

MULTICOM2 open-source protein structure prediction system powered by deep learning and distance prediction

Protein structure prediction is an important problem in bioinformatics and has been studied for decades. However, there are still few open-source comprehensive protein structure prediction packages publicly available in the field. In this paper, we present our latest open-source protein tertiary structure prediction system—MULTICOM2, an integration of template-based modeling (TBM) and template-free modeling (FM) methods. The template-based modeling uses sequence alignment tools with deep multiple sequence alignments to search for structural templates, which are much faster and more accurate than MULTICOM1. The template-free (ab initio or de novo) modeling uses the inter-residue distances predicted by DeepDist to reconstruct tertiary structure models without using any known structure as template. In the blind CASP14 experiment, the average TM-score of the models predicted by our server predictor based on the MULTICOM2 system is 0.720 for 58 TBM (regular) domains and 0.514 for 38 FM and FM/TBM (hard) domains, indicating that MULTICOM2 is capable of predicting good tertiary structures across the board. It can predict the correct fold for 76 CASP14 domains (95% regular domains and 55% hard domains) if only one prediction is made for a domain. The success rate is increased to 3% for both regular and hard domains if five predictions are made per domain. Moreover, the prediction accuracy of the pure template-free structure modeling method on both TBM and FM targets is very close to the combination of template-based and template-free modeling methods. This demonstrates that the distance-based template-free modeling method powered by deep learning can largely replace the traditional template-based modeling method even on TBM targets that TBM methods used to dominate and therefore provides a uniform structure modeling approach to any protein. Finally, on the 38 CASP14 FM and FM/TBM hard domains, MULTICOM2 server predictors (MULTICOM-HYBRID, MULTICOM-DEEP, MULTICOM-DIST) were ranked among the top 20 automated server predictors in the CASP14 experiment. After combining multiple predictors from the same research group as one entry, MULTICOM-HYBRID was ranked no. 5. The source code of MULTICOM2 is freely available at https://github.com/multicom-toolbox/multicom/tree/multicom_v2.0 .

97 MATHEMATICS AND COMPUTING↗