A Metalgel with Liquid Metal Continuum Immobilized in Polymer Network.

Gels are formed by fluids that expand throughout the whole volume of 3D polymer networks. To unlock unprecedented properties, exploring new fluids immobilized in polymer networks is crucial. Here, a new liquid metal-polymer gel material termed "metalgel" is introduced via fluid replacement strategy, featuring 92.40% vol liquid metal fluid as a continuum immobilized by interconnected nanoscale polymer network. The unique structure endows metalgel with high electrical conductivity (up to 3.18 × 106 S·m‒1), tissue-like softness (Young's modulus as low as 70 kPa), and low gas permeability (4.50 × 10‒22 m2·s‒1·Pa‒1). Besides, metalgel demonstrates electrical stability under extreme deformations, such as being run over by a 4.5-metric-tonne truck, and maintains its integrity in various environments for up to 180 days. The immobilization of high-volume-fraction liquid metal fluid is realized by electrostatic interactions is further revealed. Potential applications for metalgel are diverse and include soft electromagnetic shielding, hermetic sealing, and stimulating/sensing electrodes in implantable bioelectronics, underscoring its broad applicability.

Carbon chain elongation characterizations of electrode-biofilm microbes in electro-fermentation.

The higher efficiency of electro-fermentation in synthesizing medium-chain fatty acids (MCFAs) compared to traditional fermentation has been acknowledged. However, the functional mechanisms of electrode-biofilm enhancing MCFAs synthesis remain research gaps. To address this, this study proposed a continuous flow electrode-biofilm reactor for chain elongation (CE). After 225 days of operation, stable electrode-biofilms formed and notably improved caproate yield by more than 38 %. The electrode-biofilm was enriched with more CE microorganisms and electroactive bacteria compared to the suspended sludge microorganisms, including Caproicibacterium, Oscillibacter and Pseudoramibacter. Besides, the upregulated CE pathways were evaluated by metagenomic analysis, and the results indicated that the pathways such as acetyl-CoA and malonyl-[acp] formation, reverse beta-oxidation, and fatty acid biosynthesis pathway were all markedly enhanced in cathodic biofilm, more than anodic biofilm and suspended microorganisms. Moreover, microbial community regulated processes like bacterial chemotaxis, flagellar assembly and quorum sensing, crucial for electrode-biofilm formation. Electron transfer, energy metabolism, and microbial interactions were found to be prominently upregulated in the cathodic biofilm, surpassing levels observed in anodic biofilm and suspended sludge microorganisms, which further enhanced CE efficiency. In addition, the statistical analyses further highlighted key microbial functions and interactions within the cathodic biofilm. Oscillospiraceae_bacterium was identified to be the most active microbe, alongside pivotal roles played by Caproiciproducens_sp._NJN-50, Clostridiales_bacterium, Prevotella_sp. and Pseudoclavibacter_caeni. Eventually, the proposed microbial collaboration mechanisms of cathodic biofilm were ascertained. Overall, this study uncovered the biological effects of the electrode-biofilm on MCFAs electrosynthesis, thereby advancing biochemicals production and filling the knowledge gaps in CE electroactive biofilm reactors.

Fusion Dynamics and Size-Dependence of Droplet Microstructure in ssDNA-Mediated Protein Phase Separation.

Biomolecular condensation involving proteins and nucleic acids has been recognized to play crucial roles in genome organization and transcriptional regulation. However, the biophysical mechanisms underlying the droplet fusion dynamics and microstructure evolution during the early stage of liquid-liquid phase separation (LLPS) remain elusive. In this work, we study the phase separation of linker histone H1, which is among the most abundant chromatin proteins, in the presence of single-stranded DNA (ssDNA) capable of forming a G-quadruplex by using molecular simulations and experimental characterization. We found that droplet fusion is a rather stochastic and kinetically controlled process. Productive fusion events are triggered by the formation of ssDNA-mediated electrostatic bridges within the droplet contacting zone. The droplet microstructure is size-dependent and evolves driven by maximizing the number of electrostatic contacts. We also showed that the folding of ssDNA to the G-quadruplex promotes LLPS by increasing the multivalency and strength of protein-DNA interactions. These findings provide deep mechanistic insights into the growth dynamics of biomolecular droplets and highlight the key role of kinetic control during the early stage of ssDNA-protein condensation.

Predicting protein conformational motions using energetic frustration analysis and AlphaFold2.

Proteins perform their biological functions through motion. Although high throughput prediction of the three-dimensional static structures of proteins has proved feasible using deep-learning-based methods, predicting the conformational motions remains a challenge. Purely data-driven machine learning methods encounter difficulty for addressing such motions because available laboratory data on conformational motions are still limited. In this work, we develop a method for generating protein allosteric motions by integrating physical energy landscape information into deep-learning-based methods. We show that local energetic frustration, which represents a quantification of the local features of the energy landscape governing protein allosteric dynamics, can be utilized to empower AlphaFold2 (AF2) to predict protein conformational motions. Starting from ground state static structures, this integrative method generates alternative structures as well as pathways of protein conformational motions, using a progressive enhancement of the energetic frustration features in the input multiple sequence alignment sequences. For a model protein adenylate kinase, we show that the generated conformational motions are consistent with available experimental and molecular dynamics simulation data. Applying the method to another two proteins KaiB and ribose-binding protein, which involve large-amplitude conformational changes, can also successfully generate the alternative conformations. We also show how to extract overall features of the AF2 energy landscape topography, which has been considered by many to be black box. Incorporating physical knowledge into deep-learning-based structure prediction algorithms provides a useful strategy to address the challenges of dynamic structure prediction of allosteric proteins.

Feedback of chain elongation microorganisms on iron-based conductive materials: Enhanced microbial functions and biotoxicity adaptation mechanisms.

Despite the increased research efforts aimed at understanding iron-based conductive materials (CMs) for facilitating chain elongation (CE) to produce medium chain fatty acids (MCFAs), the impact of these materials on microbial community functions and the adaptation mechanisms to their biotoxicity remain unclear. This study found that the supply of zero-valent iron (ZVI) and magnetite enhanced the MCFAs carbon-flow distribution by 26 % and 52 %, respectively. Metagenomic analysis revealed the upregulation of fatty acid metabolism, pyruvate metabolism and ABC transporters with ZVI and magnetite. The predominant functional microorganisms were Massilibacterium and Tidjanibacter with ZVI, and were Petrimonas and Candidatus_Microthrix with magnetite. Furthermore, it was demonstrated that CE microorganisms respond and adapt to the biotoxicity of iron-based CMs by adjusting Two-component system and Quorum sensing for the first time. In summary, this study provided a new deep-insight on the feedback mechanisms of CE microorganisms on iron-based CMs.

Upgrading soybean dreg to caproate via intermediate of lactate and mediator of biochar.

To address the environmental hazards posed by high-yield soybean dreg (SD), a high-value strategy is firstly proposed by synthesizing caproate through chain elongation (CE). Optimized conditions for lactate-rich broth as intermediate, utilizing 50 % inoculum ratio, 40 g/L substrate concentration, and pH 5, resulting in 2.05 g/L caproate from direct fermentation. Leveraging lactate-rich broth supplemented with ethanol, caproate was optimized to 2.76 g/L under a refined electron donor to acceptor of 2:1. Furthermore, incorporating 20 g/L biochar elevated caproate production to 3.05 g/L and significantly shortened the lag phase. Mechanistic insights revealed that biochar's surface-existed quinone and hydroquinone groups exhibit potent redox characteristics, thereby facilitating electron transfer. Moreover, biochar up-regulated the abundance of key genes involved in CE process (especially fatty acids biosynthesis pathway), also enriching Lysinibacillus and Pseudomonas as an unrecognized cooperation to CE. This study paves a way for sustainable development of SD by upgrading to caproate.

Deciphering the Multi-state Conformational Equilibrium of HDM2 in the Regulation of p53 Binding: Perspectives from Molecular Dynamics Simulation and NMR Analysis.

HDM2 negatively regulates the activity of the tumor suppressor p53. Previous NMR studies have shown that apo-HDM2 interconverts between an "open" state in which the N-terminal "lid" is disordered and a "closed" state in which the lid covers the p53-binding site in the core region. Molecular dynamics (MD) simulation studies have been performed to elucidate the conformational dynamics of HDM2, but the direct relevance of the experimental and computational analyses is unclear. In addition, how the phosphorylation of S17 in the lid contributes to the inhibition of p53 binding remains controversial. Here, we used both NMR and MD simulations to investigate the conformational dynamics of apo-HDM2. The NMR analysis revealed that apo-HDM2 exists in a fast-exchanging equilibrium within two closed states, closed 1 and closed 2, in addition to a previously demonstrated slow-exchanging "open-closed" equilibrium. MD simulations visualized two characteristic closed states, where the spatial orientation of the key residues corresponds well to the chemical shift changes of the NMR spectra. Furthermore, the phosphorylation of S17 induced an equilibrium shift toward closed 1, thereby suppressing the binding of p53 to HDM2. This study reveals a multi-state equilibrium of apo-HDM2 and provides new insights into the regulation mechanism of HDM2-p53 interactions.

The Statistical Trends of Protein Evolution: A Lesson from AlphaFold Database.

The recent development of artificial intelligence provides us with new and powerful tools for studying the mysterious relationship between organism evolution and protein evolution. In this work, based on the AlphaFold Protein Structure Database (AlphaFold DB), we perform comparative analyses of the proteins of different organisms. The statistics of AlphaFold-predicted structures show that, for organisms with higher complexity, their constituent proteins will have larger radii of gyration, higher coil fractions, and slower vibrations, statistically. By conducting normal mode analysis and scaling analyses, we demonstrate that higher organismal complexity correlates with lower fractal dimensions in both the structure and dynamics of the constituent proteins, suggesting that higher functional specialization is associated with higher organismal complexity. We also uncover the topology and sequence bases of these correlations. As the organismal complexity increases, the residue contact networks of the constituent proteins will be more assortative, and these proteins will have a higher degree of hydrophilic-hydrophobic segregation in the sequences. Furthermore, by comparing the statistical structural proximity across the proteomes with the phylogenetic tree of homologous proteins, we show that, statistical structural proximity across the proteomes may indirectly reflect the phylogenetic proximity, indicating a statistical trend of protein evolution in parallel with organism evolution. This study provides new insights into how the diversity in the functionality of proteins increases and how the dimensionality of the manifold of protein dynamics reduces during evolution, contributing to the understanding of the origin and evolution of lives.

Simultaneous medium chain fatty acids production and process carbon emissions reduction in a continuous-flow reactor: Re-understanding of carbon flow distribution.

Due to its wide application and high value, the production of medium chain fatty acids (MCFAs) from waste biomass has become one of the worldwide research hotspots. Increasing the carbon element participation from short-chain fatty acids to the form of MCFAs is also conductive to reduce the release of biogas from biological treatment process, because carbon is in the form of MCFAs instead of biogas which directly contribute to process carbon emissions reduction. However, many barriers limiting MCFAs production and application remain to be resolved. Aiming continuous MCFAs production from lactate-rich waste biomass, this study optimized the operation conditions and clarified the main limiting factors and possible mechanisms. The maximum caproic acid concentration of 2.757 g/L were obtained at the Upflow Velocity (ULV) of 1.15 m/h and pH 4.9-5.1. Caproiciproducens, Pseudoramibacter, norank_f_Eubacteriaceae, and Oscillibacter were identified to be the dominant microbial genus responsible for MCFAs production from lactate. The reduction of carbon emissions calculation was also studied in the present processes.

Effect of substrate structure on medium chain fatty acids production and reactor microbiome.

The medium chain fatty acids (MCFAs) produced from organic wastes can replace part fossil-fuel-based products to promote the sustainable development of economy and environment. However, the selection and collocation of feedstocks for MCFAs production are lack of reference basis. This study thereby aimed to investigate how the commonly used electron donor (ED) and substrate configuration affect MCFAs synthesis and then obtain the optimal substrate composition. It was found that the optimized ratios for ethanol/acetate, lactate/acetate, and ethanol/lactate/acetate were 3/1, 2/1, and 2/1/1, respectively, and the optimal substrate concentration was 400 mM C. Combining ethanol and lactate as co-EDs effectively concentrated substrate-carbon-flow (increased by 20-28% than sole ED) on MCFAs synthesis by promoting the elongation of butyrate and reutilization of by-products. As a result, the higher MCFAs yield of 646.22 mg COD/g COD and selectivity of 67.72% were obtained from co-EDs than those from sole ED. Moreover, the key functional bacteria enriched under different ED were also discrepant, which were Clostridium sensu stricto for ethanol, Corynebacterium for lactate, and Veillonella and Oscillibacter for ethanol-lactate, respectively. This study provided a basic but significant reference for the scale-up MCFAs production.

Unraveling the Coupling between Conformational Changes and Ligand Binding in Ribose Binding Protein Using Multiscale Molecular Dynamics and Free-Energy Calculations.

Conformational changes of proteins upon ligand binding are usually explained in terms of several mechanisms including the induced fit, conformational selection, or their mixtures. Due to the slow time scales, conventional molecular dynamics (cMD) simulations based on the atomistic models cannot easily simulate the open-to-closed conformational transition in proteins. In our previous study, we have developed an enhanced sampling scheme (generalized replica exchange with solute tempering selected surface charged residues: gREST_SSCR) for multidomain proteins and applied it to ligand-mediated conformational changes in the G134R mutant of ribose-binding protein (RBPG134R) in solution. The free-energy landscape (FEL) of RBPG134R in the presence of a ribose at the binding site included the open and closed states and two intermediates, open-like and closed-like forms. Only the open and open-like forms existed in the FEL without a ribose. In the current study, the coupling between the conformational changes and ligand binding is further investigated using coarse-grained MD, multiple atomistic cMD, and free-energy calculations. The ribose is easily dissociated from the binding site of wild-type RBP and RBPG134R in the cMD simulations starting from the open and open-like forms. In contrast, it is stable at the binding site in the simulations from the closed and closed-like forms. The free-energy calculations provide the binding affinities of different structures, supporting the results of cMD simulations. Importantly, cMD simulations from the closed-like structures reveal transitions toward the closed one in the presence of a bound ribose. On the basis of the computational results, we propose a molecular mechanism in which conformational selection and induced fit happen in the first and second halves of the open-to-closed transition in RBP, respectively.

The role of gelsolin domain 3 in familial amyloidosis (Finnish type).

In the disease familial amyloidosis, Finnish type (FAF), also known as AGel amyloidosis (AGel), the mechanism by which point mutations in the calcium-regulated actin-severing protein gelsolin lead to furin cleavage is not understood in the intact protein. Here, we provide a structural and biochemical characterization of the FAF variants. X-ray crystallography structures of the FAF mutant gelsolins demonstrate that the mutations do not significantly disrupt the calcium-free conformations of gelsolin. Small-angle X-ray-scattering (SAXS) studies indicate that the FAF calcium-binding site mutants are slower to activate, whereas G167R is as efficient as the wild type. Actin-regulating studies of the gelsolins at the furin cleavage pH (6.5) show that the mutant gelsolins are functional, suggesting that they also adopt relatively normal active conformations. Deletion of gelsolin domains leads to sensitization to furin cleavage, and nanobody-binding protects against furin cleavage. These data indicate instability in the second domain of gelsolin (G2), since loss or gain of G2-stabilizing interactions impacts the efficiency of cleavage by furin. To demonstrate this principle, we engineered non-FAF mutations in G3 that disrupt the G2-G3 interface in the calcium-activated structure. These mutants led to increased furin cleavage. We carried out molecular dynamics (MD) simulations on the FAF and non-FAF mutant G2-G3 fragments of gelsolin. All mutants showed an increase in the distance between the center of masses of the 2 domains (G2 and G3). Since G3 covers the furin cleavage site on G2 in calcium-activated gelsolin, this suggests that destabilization of this interface is a critical step in cleavage.

Exploration of the folding dynamics of human telomeric G-quadruplex with a hybrid atomistic structure-based model.

Structure-based models or Go-like models, which are built from one or multiple particular experimental structures, have been successfully applied to the folding of proteins and RNAs. Recently, a variant termed the hybrid atomistic model advances the description of backbone and side chain interactions of traditional structure-based models, by borrowing the description of local interactions from classical force fields. In this study, we assessed the validity of this model in the folding problem of human telomeric DNA G-quadruplex, where local dihedral terms play important roles. A two-state model was developed and a set of molecular dynamics simulations was conducted to study the folding dynamics of sequence Htel24, which was experimentally validated to adopt two different (3 + 1) hybrid G-quadruplex topologies in K+ solution. Consistent with the experimental observations, the hybrid-1 conformation was found to be more stable and the hybrid-2 conformation was kinetically more favored. The simulations revealed that the hybrid-2 conformation folded in a higher cooperative manner, which may be the reason why it was kinetically more accessible. Moreover, by building a Markov state model, a two-quartet G-quadruplex state and a misfolded state were identified as competing states to complicate the folding process of Htel24. Besides, the simulations also showed that the transition between hybrid-1 and hybrid-2 conformations may proceed an ensemble of hairpin structures. The hybrid atomistic structure-based model reproduced the kinetic partitioning folding dynamics of Htel24 between two different folds, and thus can be used to study the complex folding processes of other G-quadruplex structures.

Metal cofactor modulated folding and target recognition of HIV-1 NCp7.

The HIV-1 nucleocapsid 7 (NCp7) plays crucial roles in multiple stages of HIV-1 life cycle, and its biological functions rely on the binding of zinc ions. Understanding the molecular mechanism of how the zinc ions modulate the conformational dynamics and functions of the NCp7 is essential for the drug development and HIV-1 treatment. In this work, using a structure-based coarse-grained model, we studied the effects of zinc cofactors on the folding and target RNA(SL3) recognition of the NCp7 by molecular dynamics simulations. After reproducing some key properties of the zinc binding and folding of the NCp7 observed in previous experiments, our simulations revealed several interesting features in the metal ion modulated folding and target recognition. Firstly, we showed that the zinc binding makes the folding transition states of the two zinc fingers less structured, which is in line with the Hammond effect observed typically in mutation, temperature or denaturant induced perturbations to protein structure and stability. Secondly, We showed that there exists mutual interplay between the zinc ion binding and NCp7-target recognition. Binding of zinc ions enhances the affinity between the NCp7 and the target RNA, whereas the formation of the NCp7-RNA complex reshapes the intrinsic energy landscape of the NCp7 and increases the stability and zinc affinity of the two zinc fingers. Thirdly, by characterizing the effects of salt concentrations on the target RNA recognition, we showed that the NCp7 achieves optimal balance between the affinity and binding kinetics near the physiologically relevant salt concentrations. In addition, the effects of zinc binding on the inter-domain conformational flexibility and folding cooperativity of the NCp7 were also discussed.

Consequences of Energetic Frustration on the Ligand-Coupled Folding/Dimerization Dynamics of Allosteric Protein S100A12.

Allosteric proteins are featured by energetic degeneracy of two (or more) functionally relevant conformations, therefore their energy landscapes are often locally frustrated. How such frustration affects the protein folding/binding dynamics is not well understood. Here, by using molecular simulations we study the consequences of local frustration in the dimerization dynamics of allosteric proteins based on a homodimer protein S100A12. Despite of the structural symmetry of the two EF-hand motifs in the three-dimensional structures, the S100A12 homodimer shows allosteric behaviors and local frustration only in half of its structural elements, i.e., the C-terminal EF-hand. We showed that such spatially asymmetric location of frustration leads to asymmetric dimerization pathways, in which the dimerization is dominantly initiated by the interchain binding of the minimally frustrated N-terminal EF-hands, achieving optimal balance between the requirements of rapid conformational switching and interchain assembling to the energy landscapes. We also showed that the local frustration, as represented by the double-basin topography of the energy landscape, gives rise to multiple cross-linked dimerization pathways, in which the dimerization is coupled with the allosteric motions of the C-terminal EF-hands. Binding of metal ions tends to reshape the energy landscape and modulate the dimerization pathways. In addition, by employing the frustratometer method, we showed that the highly frustrated residue-pairs in the C-terminal EF-hand are partially unfolded during the conformational transitions of the native homodimer, leading to lowing of free energy barrier. Our results revealed tight interplay between the local frustration of the energy landscape and the dimerization dynamics for allosteric proteins.

A novel missense mutation in the ALPL gene causes dysfunction of the protein.

Hypophosphatasia (HP) is a rare genetic disease caused by mutation in the alkaline phosphatase, liver/bone/kidney (ALPL) gene with highly variable clinical manifestations. Efforts have been made to collect cases with novel mutations and to examine how a missense mutation affects ALPL protein function, which remains difficult to predict. The present study investigated the underlying mechanism of ALPL dysfunction in a patient diagnosed with HP. Bidirectional sequencing of the ALPL gene was conducted in a 5­year­old Chinese girl preliminary diagnosed with childhood HP. Sorting Intolerant from Tolerant (SIFT) and Polymorphism Phenotyping v2 (PolyPhen­2) tools were used to forecast the impact of the mutation on protein function. Site­directed mutagenesis was performed and transfected into cells to verify the role of the specific mutation. Furthermore, the mechanism of the impact was investigated via all­atom molecular dynamics (MD) simulation. The patient demonstrated a compound heterozygote with two missense mutations in the ALPL gene, p.Trp29Arg and p.Ile395Val. Trp29 and Ile395 were determined to be 'tolerable' by SIFT, whereas they were 'possibly damaging' by PolyPhen­2 in terms of conservation. Additionally, HEK293 cells were transfected with plasmids expressing wild type and/or mutated ALPL. Only 4.1% of ALP activity remained when Trp29 was substituted by Arg, whereas 19.1, 33.7, 50.1 and 7.6% ALP activity remained in cells expressing p.Ile395Val, wild type+p.Trp29Arg, wild type+p.Ile395Val and p.Trp29Arg+p.Ile395Val substitutions, respectively. All­atom MD simulation demonstrated that the N­terminal helix of mutated ALPL, where Trp29 is located, separated from the main body of the protein after 30 nsec, and moved freely. These results demonstrated that p.Trp29Arg, as a novel missense mutation in the ALPL gene, reduced the enzymatic activity of ALPL. This effect may be associated with an uncontrolled N­terminal helix. These results provide novel information about the genetic basis of HP, and may facilitate the development of future therapies.

Piecing Together the Allosteric Patterns of Chaperonin GroEL.

Despite considerable efforts, elucidating the allostery of large macromolecular assemblies at a molecular level in solution remains technically challenging due to its structural complexity. Here we have employed an approach combining amide backbone hydrogen/deuterium exchange coupled with mass spectrometry, fluorescence spectroscopy, and molecular simulations to characterize allosteric patterns of chaperonin GroEL, an ∼800 kDa tetradecamer from E. coli. Using available crystal structures of GroEL, we quantitatively map out GroEL allosteric changes in solution by resolving exchange behaviors of 133 overlapping proteolytic peptides with more than 95% sequence coverage. This comprehensive analysis gives a refined resolution down to five residues to pilot the GroEL allosteric determinants, of which the localized dynamics is monitored by tryptophan-mutated GroEL. Furthermore, the GroEL conformational transition is evaluated by molecular dynamics simulations with an atomic-interaction-based coarse-grained model. Collectively, we provide a practical methodology to analyze GroEL allostery in solution.