The energetics and evolution of oxidoreductases in deep time.

The core metabolic reactions of life drive electrons through a class of redox protein enzymes, the oxidoreductases. The energetics of electron flow is determined by the redox potentials of organic and inorganic cofactors as tuned by the protein environment. Understanding how protein structure affects oxidation-reduction energetics is crucial for studying metabolism, creating bioelectronic systems, and tracing the history of biological energy utilization on Earth. We constructed ProtReDox (https://protein-redox-potential.web.app), a manually curated database of experimentally determined redox potentials. With over 500 measurements, we can begin to identify how proteins modulate oxidation-reduction energetics across the tree of life. By mapping redox potentials onto networks of oxidoreductase fold evolution, we can infer the evolution of electron transfer energetics over deep time. ProtReDox is designed to include user-contributed submissions with the intention of making it a valuable resource for researchers in this field.

Small protein folds at the root of an ancient metabolic network.

Life on Earth is driven by electron transfer reactions catalyzed by a suite of enzymes that comprise the superfamily of oxidoreductases (Enzyme Classification EC1). Most modern oxidoreductases are complex in their structure and chemistry and must have evolved from a small set of ancient folds. Ancient oxidoreductases from the Archean Eon between ca. 3.5 and 2.5 billion years ago have been long extinct, making it challenging to retrace evolution by sequence-based phylogeny or ancestral sequence reconstruction. However, three-dimensional topologies of proteins change more slowly than sequences. Using comparative structure and sequence profile-profile alignments, we quantify the similarity between proximal cofactor-binding folds and show that they are derived from a common ancestor. We discovered that two recurring folds were central to the origin of metabolism: ferredoxin and Rossmann-like folds. In turn, these two folds likely shared a common ancestor that, through duplication, recruitment, and diversification, evolved to facilitate electron transfer and catalysis at a very early stage in the origin of metabolism.

Anoxic photogeochemical oxidation of manganese carbonate yields manganese oxide.

The oxidation states of manganese minerals in the geological record have been interpreted as proxies for the evolution of molecular oxygen in the Archean eon. Here we report that an Archean manganese mineral, rhodochrosite (MnCO3), can be photochemically oxidized by light under anoxic, abiotic conditions. Rhodochrosite has a calculated bandgap of about 5.4 eV, corresponding to light energy centering around 230 nm. Light at that wavelength would have been present on Earth's surface in the Archean, prior to the formation of stratospheric ozone. We show experimentally that the photooxidation of rhodochrosite in suspension with light centered at 230 nm produced H2 gas and manganite (γ-MnOOH) with an apparent quantum yield of 1.37 × 10-3 moles hydrogen per moles incident photons. Our results suggest that manganese oxides could have formed abiotically on the surface in shallow waters and on continents during the Archean eon in the absence of molecular oxygen.

Biophysical analysis of the structural evolution of substrate specificity in RuBisCO.

Ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO) is the most abundant enzyme on Earth. However, its catalytic rate per molecule of protein is extremely slow and the binding of the primary substrate, CO2, is competitively displaced by O2. Hence, carbon fixation by RuBisCO is highly inefficient; indeed, in higher C3 plants, about 30% of the time the enzyme mistakes CO2 for O2 Using genomic and structural analysis, we identify regions around the catalytic site that play key roles in discriminating between CO2 and O2 Our analysis identified positively charged cavities directly around the active site, which are expanded as the enzyme evolved with higher substrate specificity. The residues that extend these cavities have recently been under selective pressure, indicating that larger charged pockets are a feature of modern RuBisCOs, enabling greater specificity for CO2 This paper identifies a key structural feature that enabled the enzyme to evolve improved CO2 sequestration in an oxygen-rich atmosphere and may guide the engineering of more efficient RuBisCOs.

Bioactive compounds of corn silk and their role in management of glycaemic response.

Management of glycaemic response is perhaps the most critical part of antidiabetic therapy. Hypoglycaemia is an avoidable complication caused by conventional drugs used in the treatment of diabetes. It triggers commonly during the intensification of anti-hyperglycemic therapy used to render glycemic control in diabetic patients. The commercial oral hypoglycaemic drugs, insulin, herbal medicines and plant extracts are therefore used as a part of the treatment of diabetes. The demand for treating diabetes, through herbal and plant resources is due to their lesser adverse reactions and better phytochemical benefits. Corn silk has been shown to have anti-allergic, anti-inflammatory, and anti-hypertensive effects when extracted in various solvents. Corn silk has medicinal characteristics and has long been used as a traditional medicine in many nations, although the mechanism of action is unknown. The hypoglycaemic effects of corn silk are investigated in this review. The phytochemical components present in corn silk-like flavonoids, phenolics, terpenoids, tannins, sterols, and alkaloids are phytochemical components that have hypoglycemic activity and a mechanism for lowering blood glucose levels. There is a lack of a homogenized database on the hypoglycemic properties of corn silk thus the present review attempts to critically analyse it and provide specific recommendations of its doses.

Evaluating Mineral Lattices as Evolutionary Proxies for Metalloprotein Evolution.

Protein coordinated iron-sulfur clusters drive electron flow within metabolic pathways for organisms throughout the tree of life. It is not known how iron-sulfur clusters were first incorporated into proteins. Structural analogies to iron-sulfide minerals present on early Earth, suggest a connection in the evolution of both proteins and minerals. The availability of large protein and mineral crystallographic structure data sets, provides an opportunity to explore co-evolution of proteins and minerals on a large-scale using informatics approaches. However, quantitative comparisons are confounded by the infinite, repeating nature of the mineral lattice, in contrast to metal clusters in proteins, which are finite in size. We address this problem using the Niggli reduction to transform a mineral lattice to a finite, unique structure that when translated reproduces the crystal lattice. Protein and reduced mineral structures were represented as quotient graphs with the edges and nodes corresponding to bonds and atoms, respectively. We developed a graph theory-based method to calculate the maximum common connected edge subgraph (MCCES) between mineral and protein quotient graphs. MCCES can accommodate differences in structural volumes and easily allows additional chemical criteria to be considered when calculating similarity. To account for graph size differences, we use the Tversky similarity index. Using consistent criteria, we found little similarity between putative ancient iron-sulfur protein clusters and iron-sulfur mineral lattices, suggesting these metal sites are not as evolutionarily connected as once thought. We discuss possible evolutionary implications of these findings in addition to suggesting an alternative proxy, mineral surfaces, for better understanding the coevolution of the geosphere and biosphere.

De novo design of symmetric ferredoxins that shuttle electrons in vivo.

A symmetric origin for bacterial ferredoxins was first proposed over 50 y ago, yet, to date, no functional symmetric molecule has been constructed. It is hypothesized that extant proteins have drifted from their symmetric roots via gene duplication followed by mutations. Phylogenetic analyses of extant ferredoxins support the independent evolution of N- and C-terminal sequences, thereby allowing consensus-based design of symmetric 4Fe-4S molecules. All designs bind two [4Fe-4S] clusters and exhibit strongly reducing midpoint potentials ranging from -405 to -515 mV. One of these constructs efficiently shuttles electrons through a designed metabolic pathway in Escherichia coli These finding establish that ferredoxins consisting of a symmetric core can be used as a platform to design novel electron transfer carriers for in vivo applications. Outer-shell asymmetry increases sequence space without compromising electron transfer functionality.

Biallelic hypomorphic variants in ALDH1A2 cause a novel lethal human multiple congenital anomaly syndrome encompassing diaphragmatic, pulmonary, and cardiovascular defects.

This study shows a causal association between ALDH1A2 variants and a novel, severe multiple congenital anomaly syndrome in humans that is neonatally lethal due to associated pulmonary hypoplasia and respiratory failure. In two families, exome sequencing identified compound heterozygous missense variants in ALDH1A2. ALDH1A2 is involved in the conversion of retinol (vitamin A) into retinoic acid (RA), which is an essential regulator of diaphragm and cardiovascular formation during embryogenesis. Reduced RA causes cardiovascular, diaphragmatic, and associated pulmonary defects in several animal models, matching the phenotype observed in our patients. In silico protein modeling showed probable impairment of ALDH1A2 for three of the four substitutions. In vitro studies show a reduction of RA. Few pathogenic variants in genes encoding components of the retinoic signaling pathway have been described to date, likely due to embryonic lethality. Thus, this study contributes significantly to knowledge of the role of this pathway in human diaphragm and cardiovascular development and disease. Some clinical features in our patients are also observed in Fryns syndrome (MIM# 229850), syndromic microphthalmia 9 (MIM# 601186), and DiGeorge syndrome (MIM# 188400). Patients with similar clinical features who are genetically undiagnosed should be tested for recessive ALDH1A2-deficient malformation syndrome.

Design Strategies to Tune the Structural and Mechanical Properties of Synthetic Collagen Hydrogels.

As an important component of biomaterials, collagen provides three-dimensional scaffolds and biological cues for cell adhesion and proliferation in tissue engineering. Recombinant collagen-like proteins, which were initially discovered in Streptococcus pyogenes and produced in heterologous hosts, have been chemically and genetically engineered for biomaterial applications. However, existing collagen-like proteins do not form gels, limiting their utility as biomaterials. Here, we present a series of rationally designed collagen-like proteins composed of a trimerization domain, triple-helical domains with various lengths, and a pair of heterotrimeric coiled-coil sequences attached to the N- and C-termini as adhesive ends. These designed proteins fold into triple helices and form self-supporting gels. As the triple-helical domains are lengthened, the gels become less stiff, pore sizes increase, and structural anisotropy decreases. Moreover, cell-culture assay confirms that the designed proteins are noncytotoxic. This study provides a design strategy for collagen-based biomaterials. The sequence variations reveal a relationship between the protein primary structure and material properties, where variations in the cross-linking density and association energies define the gelation of the protein network.

Lysosomal protein thermal stability does not correlate with cellular half-life: global observations and a case study of tripeptidyl-peptidase 1.

Late-infantile neuronal ceroid lipofuscinosis (LINCL) is a neurodegenerative lysosomal storage disorder caused by mutations in the gene encoding the protease tripeptidyl-peptidase 1 (TPP1). Progression of LINCL can be slowed or halted by enzyme replacement therapy, where recombinant human TPP1 is administered to patients. In this study, we utilized protein engineering techniques to increase the stability of recombinant TPP1 with the rationale that this may lengthen its lysosomal half-life, potentially increasing the potency of the therapeutic protein. Utilizing multiple structure-based methods that have been shown to increase the stability of other proteins, we have generated and evaluated over 70 TPP1 variants. The most effective mutation, R465G, increased the melting temperature of TPP1 from 55.6°C to 64.4°C and increased its enzymatic half-life at 60°C from 5.4 min to 21.9 min. However, the intracellular half-life of R465G and all other variants tested in cultured LINCL patient-derived lymphoblasts was similar to that of WT TPP1. These results provide structure/function insights into TPP1 and indicate that improving in vitro thermal stability alone is insufficient to generate TPP1 variants with improved physiological stability. This conclusion is supported by a proteome-wide analysis that indicates that lysosomal proteins have higher melting temperatures but also higher turnover rates than proteins of other organelles. These results have implications for similar efforts where protein engineering approaches, which are frequently evaluated in vitro, may be considered for improving the physiological properties of proteins, particularly those that function in the lysosomal environment.

Marginal protein stability drives subcellular proteome isoelectric point.

There exists a positive correlation between the pH of subcellular compartments and the median isoelectric point (pI) for the associated proteomes. Proteins in the human lysosome-a highly acidic compartment in the cell-have a median pI of ∼6.5, whereas proteins in the more basic mitochondria have a median pI of ∼8.0. Proposed mechanisms reflect potential adaptations to pH. For example, enzyme active site general acid/base residue pKs are likely evolved to match environmental pH. However, such effects would be limited to a few residues on specific proteins, and might not affect the proteome at large. A protein model that considers residue burial upon folding recapitulates the correlation between proteome pI and environmental pH. This correlation can be fully described by a neutral evolution process; no functional selection is included in the model. Proteins in acidic environments incur a lower energetic penalty for burying acidic residues than basic residues, resulting in a net accumulation of acidic residues in the protein core. The inverse is true under alkaline conditions. The pI distributions of subcellular proteomes are likely not a direct result of functional adaptations to pH, but a molecular spandrel stemming from marginal stability.

How electrostatic networks modulate specificity and stability of collagen.

One-quarter of the 28 types of natural collagen exist as heterotrimers. The oligomerization state of collagen affects the structure and mechanics of the extracellular matrix, providing essential cues to modulate biological and pathological processes. A lack of high-resolution structural information limits our mechanistic understanding of collagen heterospecific self-assembly. Here, the 1.77-Å resolution structure of a synthetic heterotrimer demonstrates the balance of intermolecular electrostatics and hydrogen bonding that affects collagen stability and heterospecificity of assembly. Atomistic simulations and mutagenesis based on the solved structure are used to explore the contributions of specific interactions to energetics. A predictive model of collagen stability and specificity is developed for engineering novel collagen structures.

Modular origins of biological electron transfer chains.

Oxidoreductases catalyze electron transfer reactions that ultimately provide the energy for life. A limited set of ancestral protein-metal modules are presumably the building blocks that evolved into this diverse protein family. However, the identity of these modules and their path to modern oxidoreductases is unknown. Using a comparative structural analysis approach, we identify a set of fundamental electron transfer modules that have evolved to form the extant oxidoreductases. Using transition metal-containing cofactors as fiducial markers, it is possible to cluster cofactor microenvironments into as few as four major modules: bacterial ferredoxin, cytochrome c, symerythrin, and plastocyanin-type folds. From structural alignments, it is challenging to ascertain whether modules evolved from a single common ancestor (homology) or arose by independent convergence on a limited set of structural forms (analogy). Additional insight into common origins is contained in the spatial adjacency network (SPAN), which is based on proximity of modules in oxidoreductases containing multiple cofactor electron transfer chains. Electron transfer chains within complex modern oxidoreductases likely evolved through repeated duplication and diversification of ancient modular units that arose in the Archean eon.

Comparative dynamics of tropomyosin in vertebrates and invertebrates.

Tropomyosin (Tpm) is an extended α-helical coiled-coil homodimer that regulates actinomyosin interactions in muscle. Molecular simulations of four Tpms, two from the vertebrate class Mammalia (rat and pig), and two from the invertebrate class Malacostraca (shrimp and lobster), showed that despite extensive sequence and structural homology across metazoans, dynamic behavior-particularly long-range structural fluctuations-were clearly distinct. Vertebrate Tpms were more flexible and sampled complex, multi-state conformational landscapes. Invertebrate Tpms were more rigid, sampling a highly constrained harmonic landscape. Filtering of trajectories by principle component analysis into essential subspaces showed significant overlap within but not between phyla. In vertebrate Tpms, hinge-regions decoupled long-range interhelical motions and suggested distinct domains. In contrast, crustacean Tpms did not exhibit long-range dynamic correlations-behaving more like a single rigid rod on the nanosecond time scale. These observations suggest there may be divergent mechanisms for Tpm binding to actin filaments, where conformational flexibility in mammalian Tpm allows a preorganized shape complementary to the filament surface, and where rigidity in the crustacean Tpm requires concerted bending and binding.

Zeta Potential Prediction from Protein Structure in General Aqueous Electrolyte Solutions.

The ζ, or electrokinetic, potential is the effective charge energy of a molecule in a solution, defining its electrostatic interactions in the solution. A computational protocol for computing ζ potential from the high-resolution structures of proteins (ZPRED) is described. This model considers both protein and solution components and incorporates a number of electrokinetic models that account for many of the complexities of protein electrophoresis. Experimental observations of electrophoretic mobilities using a benchtop light scattering instrument match computed mobilities for different proteins over a wide range of aqueous solution conditions. ZPRED is a tool for optimizing protein sequence and solution conditions (pH, ionic composition and strength, temperature) to disperse molecules by charge repulsion, preventing aggregation. This is an important factor in enhancing the stability of engineered biologics or industrial protein catalysts.

Design of a Fe4 S4 cluster into the core of a de novo four-helix bundle.

We explore the capacity of the de novo protein, S824, to incorporate a multinuclear iron-sulfur cluster within the core of a single-chain four-helix bundle. This topology has a high intrinsic designability because sequences are constrained largely by the pattern of hydrophobic and hydrophilic amino acids, thereby allowing for the extensive substitution of individual side chains. Libraries of novel proteins based on these constraints have surprising functional potential and have been shown to complement the deletion of essential genes in E. coli. Our structure-based design of four first-shell cysteine ligands, one per helix, in S824 resulted in successful incorporation of a cubane Fe4 S4 cluster into the protein core. A number of challenges were encountered during the design and characterization process, including nonspecific metal-induced aggregation and the presence of competing metal-cluster stoichiometries. The introduction of buried iron-sulfur clusters into the helical bundle is an initial step toward converting libraries of designed structures into functional de novo proteins with catalytic or electron-transfer functionalities.

Sugar profile and rheological behaviour of four different Indian honey varieties.

The aim of the study was to quantify sugar profile and rheological behaviour of four Indian honey varieties (Cotton, Coriander, Dalbergia, and Murraya). The effect of temperature (5, 10, 20 and 30 °C) on rheological behaviour of these honey varieties was also studied. Fourteen sugars (three monosaccharides, six disaccharides, four trisaccharides, and one oligosaccharide) were quantified. The concentration of glucose and fructose varied from 33.40-34.06% to 36.86-41.15%, respectively. Result indicated that monosaccharides were the dominant sugars among all the honey samples. Low amounts of turanose, trehalose, melibiose, and raffinose were present in the range of 0.02-0.03%, 0.11-0.26%, 0.09-0.18% and 0.06-0.12%, respectively in the analyzed honey varieties. The rheological behaviour of all four honey varieties was analysed at different temperatures i.e. 5, 10, 20 and 30 °C followed an Arrhenius model. In analysed honey varieties storage modulus (G') was less than loss modulus (G″) which confirmed the Newtonian behaviour of all honey samples.

Catalytic linkage between caspase activity and proteostasis in Archaea.

The model haloarchaeon, Haloferax volcanii possess an extremely high, and highly specific, basal caspase activity in exponentially growing cells that closely resembles caspase-4. This activity is specifically inhibited by the pan-caspase inhibitor, z-VAD-FMK, and has no cross-reactivity with other known protease families. Although it is one of the dominant cellular proteolytic activities in exponentially growing H. volcanii cells, the interactive cellular roles remain unknown and the protein(s) responsible for this activity remain elusive. Here, biochemical purification and in situ trapping with caspase targeted covalent inhibitors combined with genome-enabled proteomics, structural analysis, targeted gene knockouts and treatment with canavanine demonstrated a catalytic linkage between caspase activity and thermosomes, proteasomes and cdc48b, a cell division protein and proteasomal degradation facilitating ATPase, as part of an 'interactase' of stress-related protein complexes with an established link to the unfolded protein response (UPR). Our findings provide novel cellular and biochemical context for the observed caspase activity in Archaea and add new insight to understanding the role of this activity, implicating their possible role in the establishment of protein stress and ER associated degradation pathways in Eukarya.

Optimization of process parameters on hydroxymethylfurfural content, diastase and invertase activity of coriander honey.

Floral authenticity of coriander (Coriandrum sativum) honey samples was confirmed by melissopalynology. Effect of temperature, time and pH on quality parameters i.e. hydroxymethylfurfural (HMF) content, diastase and invertase activity of coriander honey was analysed using response surface methodology. Central composite rotatable design was adopted for optimization of process variables. An increased in HMF content was observed with increase in temperature and pH whereas diastase activity decreased with increase in temperature and with a pH value other than the optimum value of 4.6-5.6. Invertase activity was maximum at 4.8 pH. Interaction effect of temperature and pH was significant for HMF whereas interaction effect of temperature and time was significant for HMF, diastase and invertase activity. Optimization of variables was done by the mathematical method, and optimized values of HMF content, diastase, and invertase activity were obtained as 7.78 (mg/kg), 17.95 DN and 13.96 IN, respectively at 47.5 °C (temperature), 4.7 (pH) and 9 min (time).

Hydrodynamic radius coincides with the slip plane position in the electrokinetic behavior of lysozyme.

The zeta potential (ζ) is the effective charge energy of a solvated protein, describing the magnitude of electrostatic interactions in solution. It is commonly used in the assessment of adsorption processes and dispersion stability. Predicting ζ from molecular structure would be useful to the structure-based molecular design of drugs, proteins, and other molecules that hold charge-dependent function while remaining suspended in solution. One challenge in predicting ζ is identifying the location of the slip plane (XSP ), a distance from the protein surface where ζ is theoretically defined. This study tests the hypothesis that the XSP can be estimated by the Stokes-Einstein hydrodynamic radius (Rh ), using globular hen egg white lysozyme as a model system. Although the XSP and Rh differ in their theoretical definitions, with the XSP being the position of the ζ during electrokinetic phenomena (e.g., electrophoresis) and the Rh being a radius pertaining to the edge of solvation during diffusion, they both represent the point where water and ions no longer adhere to a molecule. This work identifies the limited range of ionic strengths in which the XSP can be determined using diffusivity measurements and the Stokes-Einstein equation. In addition, a computational protocol is developed for determining the ζ from a protein crystal structure. At low ionic strengths, a hyperdiffusivity regime exists, requiring direct measurement of electrophoretic mobility to determine ζ. This work, therefore, supports a basic tenant of EDL theory that the electric double layer during diffusion and electrophoresis are equivalent in the Stokes-Einstein regime.