Characterization of alternate encounter assemblies of SARS-CoV-2 main protease.

The assembly of two monomeric constructs spanning segments 1-199 (MPro1-199) and 10-306 (MPro10-306) of SARS-CoV-2 main protease (MPro) was examined to assess the existence of a transient heterodimer intermediate in the N-terminal autoprocessing pathway of MPro model precursor. Together, they form a heterodimer population accompanied by a 13-fold increase in catalytic activity. Addition of inhibitor GC373 to the proteins increases the activity further by ∼7-fold with a 1:1 complex and higher order assemblies approaching 1:2 and 2:2 molecules of MPro1-199 and MPro10-306 detectable by analytical ultracentrifugation and native mass estimation by light scattering. Assemblies larger than a heterodimer (1:1) are discussed in terms of alternate pathways of domain III association, either through switching the location of helix 201 to 214 onto a second helical domain of MPro10-306 and vice versa or direct interdomain III contacts like that of the native dimer, based on known structures and AlphaFold 3 prediction, respectively. At a constant concentration of MPro1-199 with molar excess of GC373, the rate of substrate hydrolysis displays first order dependency on the MPro10-306 concentration and vice versa. An equimolar composition of the two proteins with excess GC373 exhibits half-maximal activity at ∼6 µM MPro1-199. Catalytic activity arises primarily from MPro1-199 and is dependent on the interface interactions involving the N-finger residues 1 to 9 of MPro1-199 and E290 of MPro10-306. Importantly, our results confirm that a single N-finger region with its associated intersubunit contacts is sufficient to form a heterodimeric MPro intermediate with enhanced catalytic activity.

Comparison of structural networks across homologous proteins.

Protein sequence determines its structure and function. The indirect relationship between protein function and structure lies deep-rooted in the structural topology that has evolved into performing optimal function. The evolution of structure and its interconnectivity has been conventionally studied by comparing the root means square deviation between protein structures at the backbone level. Two factors that are necessary for the quantitative comparison of non-covalent interactions are (a) explicit inclusion of the coordinates of side-chain atoms and (b) consideration of multiple structures from the conformational landscape to account for structural variability. We have recently addressed these fundamental issues by investigating the alteration of inter-residue interactions across an ensemble of protein structure networks through a graph spectral approach. In this study, we have developed a rigorous method to compare the structure networks of homologous proteins, with a wide range of sequence identity percentages. A range of dissimilarity measures that show the extent of change in the network across homologous structures are generated, which also includes the comparison of the protein structure variability. We discuss in detail, scenarios where the variation of structure is not accompanied by loss or gain of the overall network and its vice versa. The sequence-based phylogeny among the homologs is also compared with the lineage obtained from information from such a robust structure comparison. In summary, we can obtain a quantitative comparison score for the structure networks of homologous proteins, which also enables us to study the evolution of protein function based on the variation of their topologies.

Molecular investigation of the tandem Tudor domain and plant homeodomain histone binding domains of the epigenetic regulator UHRF2.

Ubiquitin-like containing PHD and ring finger (UHRF)1 and UHRF2 are multidomain epigenetic proteins that play a critical role in bridging crosstalk between histone modifications and DNA methylation. Both proteins contain two histone reader domains, called tandem Tudor domain (TTD) and plant homeodomain (PHD), which read the modification status on histone H3 to regulate DNA methylation and gene expression. To shed light on the mechanism of histone binding by UHRF2, we have undergone a detailed molecular investigation with the TTD, PHD and TTD-PHD domains and compared the binding activity to its UHRF1 counterpart. We found that unlike UHRF1 where the PHD is the primary binding contributor, the TTD of UHRF2 has modestly higher affinity toward the H3 tail, while the PHD has a weaker binding interaction. We also demonstrated that like UHRF1, the aromatic amino acids within the TTD are important for binding to H3K9me3 and a conserved aspartic acid within the PHD forms an ionic interaction with R2 of H3. However, while the aromatic amino acids in the TTD of UHRF1 contribute to selectivity, the analogous residues in UHRF2 contribute to both selectivity and affinity. We also discovered that the PHD of UHRF2 contains a distinct asparagine in the H3R2 binding pocket that lowers the binding affinity of the PHD by reducing a potential electrostatic interaction with the H3 tail. Furthermore, we demonstrate the PHD and TTD of UHRF2 cooperate to interact with the H3 tail and that dual domain engagement with the H3 tail relies on specific amino acids. Lastly, our data indicate that the unique stretch region in the TTD of UHRF2 can decrease the melting temperature of the TTD-PHD and represents a disordered region. Thus, these subtle but important mechanistic differences are potential avenues for selectively targeting the histone binding interactions of UHRF1 and UHRF2 with small molecules.

A study of a hierarchical structure of proteins and ligand binding sites of receptors using the triangular spatial relationship-based structure comparison method and development of a size-filtering feature designed for comparing different sizes of protein structures.

The presence of receptors and the specific binding of the ligands determine nearly all cellular responses. Binding of a ligand to its receptor causes conformational changes of the receptor that triggers the subsequent signaling cascade. Therefore, systematically studying structures of receptors will provide insight into their functions. We have developed the triangular spatial relationship (TSR)-based method where all possible triangles are constructed with Cα atoms of a protein as vertices. Every triangle is represented by an integer denoted as a "key" computed through the TSR algorithm. A structure is thereby represented by a vector of integers. In this study, we have first defined substructures using different types of keys. Second, using different types of keys represents a new way to interpret structure hierarchical relations and differences between structures and sequences. Third, we demonstrate the effects of sequence similarity as well as sample size on the structure-based classifications. Fourth, we show identification of structure motifs, and the motifs containing multiple triangles connected by either an edge or a vertex are mapped to the ligand binding sites of the receptors. The structure motifs are valuable resources for the researchers in the field of signal transduction. Next, we propose amino-acid scoring matrices that capture "evolutionary closeness" information based on BLOSUM62 matrix, and present the development of a new visualization method where keys are organized according to evolutionary closeness and shown in a 2D image. This new visualization opens a window for developing tools with the aim of identification of specific and common substructures by scanning pixels and neighboring pixels. Finally, we report a new algorithm called as size filtering that is designed to improve structure comparison of large proteins with small proteins. Collectively, we provide an in-depth interpretation of structure relations through the detailed analyses of different types of keys and their associated key occurrence frequencies, geometries, and labels. In summary, we consider this study as a new computational platform where keys are served as a bridge to connect sequence and structure as well as structure and function for a deep understanding of sequence, structure, and function relationships of the protein family.

Analyses of displacements resulting from a point mutation in proteins.

The effects of a single residue substitution on the protein backbone are frequently quite small and there are many other potential sources of structural variation for protein. We present here a methodology considering different sources of distortions in order to isolate the very effect of the mutation. To validate our methodology, we consider a well-studied family with many single mutants: the human lysozyme. Most of the perturbations are expected to be at the very localisation of the mutation, but in many cases the effects are propagated at long range. We show that the distances between the mutated residue and the 5% most disturbed residues exponentially decreases. One third of the affected residues are in direct contact with the mutated position; the remaining two thirds are potential allosteric effects. We confirm the reliability of the residues identified as significantly perturbed by comparing our results to experimental studies. We confirm with the present method all the previously identified perturbations. This study shows that mutations have long-range impact on protein backbone that can be detected, although the displacement of the affected atoms is small.

Molecular characterization and expression of TLR7 and TLR8 in barbel chub (Squaliobarbus curriculus): Responses to stimulation of grass carp reovirus and lipopolysaccharide.

The barbel chub (Squaliobarbus curriculus) is a kind of small size commercial fish species that is widely spread in Asia and has shown significant resistance to disease. In this study, the full-length cDNA sequences of Toll-like receptor (TLR) 7 and 8 from S. curriculus, designated as ScTLR7 and ScTLR8, were cloned, and their differences in the structure and the responses to the grass carp (GCRV) infection and lipopolysaccharide stimulation were investigated. The full-length 3715 base pair (bp) cDNA of ScTLR7 contained a complete open reading frame of 3162 bp and encoded a putative polypeptide of 1053 amino acid residues. The full-length 4624 base pair (bp) cDNA of ScTLR8 contained a complete open reading frame of 3072 bp and encoded a putative polypeptide of 1023 amino acid residues. ScTLR7 and ScTLR8 consisted of N-terminal signal peptide, leucine-rich repeats (LRRs), and Toll/IL-1 Receptors domain. LRR motifs of ScTLR7 and ScTLR8 bend into horseshoe-like solenoid structure, while the number of LRRs between the two genes is different. Phylogenetic analysis showed that both the ScTLR7 and ScTLR8 were closely clustered with Ctenopharyngodon idellus and Megalobrama amblycephala. Quantitative real-time polymerase chain reaction analysis showed that the expression levels of ScTLR7 in S. curriculus were most abundant in the brain followed by the spleen and heart, and the lowest in the intestine. The highest expression level of ScTLR8 was observed in spleen and the lowest in the liver. After LPS stimulation, the relative expression levels of both ScTLR7 and ScTLR8 exhibited an overall trend of up-regulation. The expression levels of type I-IFN showed an overall trend of down-regulation at time points of 12, 24, 72 and 168 h compared to that of 6 h after LPS stimulation. Compared to 6 h post GCRV infection, the transcription level of ScTLR7 was up-regulated from 12 to 168 h, and transcription levels of ScTLR8, MyD88, and type I-IFN were firstly up-regulated from 12 to 72 h, and then down-regulated from 72 to 168 h. Correlation analysis showed that expression level of ScTLR7 in the spleen was significantly positively correlated with that of MyD88 (Pearson correlation coefficient: 0.909, P: 0.033), and a significantly positive correlation was also observed between expression levels of MyD88 and type I IFN (Pearson correlation coefficient: 0.962, P: 0.009), after GCRV stimulation. These results indicate that ScTLR7 and ScTLR8 may play important roles in the responses to the grass carp (GCRV) infection and lipopolysaccharide stimulation and trigger different downstream immune signal pathways.

LCS-TA to identify similar fragments in RNA 3D structures.

BACKGROUND: In modern structural bioinformatics, comparison of molecular structures aimed to identify and assess similarities and differences between them is one of the most commonly performed procedures. It gives the basis for evaluation of in silico predicted models. It constitutes the preliminary step in searching for structural motifs. In particular, it supports tracing the molecular evolution. Faced with an ever-increasing amount of available structural data, researchers need a range of methods enabling comparative analysis of the structures from either global or local perspective. RESULTS: Herein, we present a new, superposition-independent method which processes pairs of RNA 3D structures to identify their local similarities. The similarity is considered in the context of structure bending and bonds' rotation which are described by torsion angles. In the analyzed RNA structures, the method finds the longest continuous segments that show similar torsion within a user-defined threshold. The length of the segment is provided as local similarity measure. The method has been implemented as LCS-TA algorithm (Longest Continuous Segments in Torsion Angle space) and is incorporated into our MCQ4Structures application, freely available for download from http://www.cs.put.poznan.pl/tzok/mcq/ . CONCLUSIONS: The presented approach ties torsion-angle-based method of structure analysis with the idea of local similarity identification by handling continuous 3D structure segments. The first method, implemented in MCQ4Structures, has been successfully utilized in RNA-Puzzles initiative. The second one, originally applied in Euclidean space, is a component of LGA (Local-Global Alignment) algorithm commonly used in assessing protein models submitted to CASP. This unique combination of concepts implemented in LCS-TA provides a new perspective on structure quality assessment in local and quantitative aspect. A series of computational experiments show the first results of applying our method to comparison of RNA 3D models. LCS-TA can be used for identifying strengths and weaknesses in the prediction of RNA tertiary structures.

RNA-TVcurve: a Web server for RNA secondary structure comparison based on a multi-scale similarity of its triple vector curve representation.

BACKGROUND: RNAs have been found to carry diverse functionalities in nature. Inferring the similarity between two given RNAs is a fundamental step to understand and interpret their functional relationship. The majority of functional RNAs show conserved secondary structures, rather than sequence conservation. Those algorithms relying on sequence-based features usually have limitations in their prediction performance. Hence, integrating RNA structure features is very critical for RNA analysis. Existing algorithms mainly fall into two categories: alignment-based and alignment-free. The alignment-free algorithms of RNA comparison usually have lower time complexity than alignment-based algorithms. RESULTS: An alignment-free RNA comparison algorithm was proposed, in which novel numerical representations RNA-TVcurve (triple vector curve representation) of RNA sequence and corresponding secondary structure features are provided. Then a multi-scale similarity score of two given RNAs was designed based on wavelet decomposition of their numerical representation. In support of RNA mutation and phylogenetic analysis, a web server (RNA-TVcurve) was designed based on this alignment-free RNA comparison algorithm. It provides three functional modules: 1) visualization of numerical representation of RNA secondary structure; 2) detection of single-point mutation based on secondary structure; and 3) comparison of pairwise and multiple RNA secondary structures. The inputs of the web server require RNA primary sequences, while corresponding secondary structures are optional. For the primary sequences alone, the web server can compute the secondary structures using free energy minimization algorithm in terms of RNAfold tool from Vienna RNA package. CONCLUSION: RNA-TVcurve is the first integrated web server, based on an alignment-free method, to deliver a suite of RNA analysis functions, including visualization, mutation analysis and multiple RNAs structure comparison. The comparison results with two popular RNA comparison tools, RNApdist and RNAdistance, showcased that RNA-TVcurve can efficiently capture subtle relationships among RNAs for mutation detection and non-coding RNA classification. All the relevant results were shown in an intuitive graphical manner, and can be freely downloaded from this server. RNA-TVcurve, along with test examples and detailed documents, are available at: http://ml.jlu.edu.cn/tvcurve/ .

Validation of protein structure models using network similarity score.

Accurate structural validation of proteins is of extreme importance in studies like protein structure prediction, analysis of molecular dynamic simulation trajectories and finding subtle changes in very similar structures. The benchmarks for today's structure validation are scoring methods like global distance test-total structure (GDT-TS), TM-score and root mean square deviations (RMSD). However, there is a lack of methods that look at both the protein backbone and side-chain structures at the global connectivity level and provide information about the differences in connectivity. To address this gap, a graph spectral based method (NSS-network similarity score) which has been recently developed to rigorously compare networks in diverse fields, is adopted to compare protein structures both at the backbone and at the side-chain noncovalent connectivity levels. In this study, we validate the performance of NSS by investigating protein structures from X-ray structures, modeling (including CASP models), and molecular dynamics simulations. Further, we systematically identify the local and the global regions of the structures contributing to the difference in NSS, through the components of the score, a feature unique to this spectral based scoring scheme. It is demonstrated that the method can quantify subtle differences in connectivity compared to a reference protein structure and can form a robust basis for protein structure comparison. Additionally, we have also introduced a network-based method to analyze fluctuations in side chain interactions (edge-weights) in an ensemble of structures, which can be an useful tool for the analysis of MD trajectories.

A comparative analysis of the foamy and ortho virus capsid structures reveals an ancient domain duplication.

BACKGROUND: The Spumaretrovirinae (foamy viruses) and the Orthoretrovirinae (e.g. HIV) share many similarities both in genome structure and the sequences of the core viral encoded proteins, such as the aspartyl protease and reverse transcriptase. Similarity in the gag region of the genome is less obvious at the sequence level but has been illuminated by the recent solution of the foamy virus capsid (CA) structure. This revealed a clear structural similarity to the orthoretrovirus capsids but with marked differences that left uncertainty in the relationship between the two domains that comprise the structure. METHODS: We have applied protein structure comparison methods in order to try and resolve this ambiguous relationship. These included both the DALI method and the SAP method, with rigorous statistical tests applied to the results of both methods. For this, we employed collections of artificial fold 'decoys' (generated from the pair of native structures being compared) to provide a customised background distribution for each comparison, thus allowing significance levels to be estimated. RESULTS: We have shown that the relationship of the two domains conforms to a simple linear correspondence rather than a domain transposition. These similarities suggest that the origin of both viral capsids was a common ancestor with a double domain structure. In addition, we show that there is also a significant structural similarity between the amino and carboxy domains in both the foamy and ortho viruses. CONCLUSIONS: These results indicate that, as well as the duplication of the double domain capsid, there may have been an even more ancient gene-duplication that preceded the double domain structure. In addition, our structure comparison methodology demonstrates a general approach to problems where the components have a high intrinsic level of similarity.

Development of a motif-based topology-independent structure comparison method to identify evolutionarily related folds.

Structure conservation, functional similarities, and homologous relationships that exist across diverse protein topologies suggest that some regions of the protein fold universe are continuous. However, the current structure classification systems are based on hierarchical organizations, which cannot accommodate structural relationships that span fold definitions. Here, we describe a novel, super-secondary-structure motif-based, topology-independent structure comparison method (SmotifCOMP) that is able to quantitatively identify structural relationships between disparate topologies. The basis of SmotifCOMP is a systematically defined super-secondary-structure motif library whose representative geometries are shown to be saturated in the Protein Data Bank and exhibit a unique distribution within the known folds. SmotifCOMP offers a robust and quantitative technique to compare domains that adopt different topologies since the method does not rely on a global superposition. SmotifCOMP is used to perform an exhaustive comparison of the known folds and the identified relationships are used to produce a nonhierarchical representation of the fold space that reflects the notion of a continuous and connected fold universe. The current work offers insight into previously hypothesized evolutionary relationships between disparate folds and provides a resource for exploring novel ones. Proteins 2016; 84:1859-1874. © 2016 Wiley Periodicals, Inc.

Evaluation of free modeling targets in CASP11 and ROLL.

We present an assessment of 'template-free modeling' (FM) in CASP11and ROLL. Community-wide server performance suggested the use of automated scores similar to previous CASPs would provide a good system of evaluating performance, even in the absence of comprehensive manual assessment. The CASP11 FM category included several outstanding examples, including successful prediction by the Baker group of a 256-residue target (T0806-D1) that lacked sequence similarity to any existing template. The top server model prediction by Zhang's Quark, which was apparently selected and refined by several manual groups, encompassed the entire fold of target T0837-D1. Methods from the same two groups tended to dominate overall CASP11 FM and ROLL rankings. Comparison of top FM predictions with those from the previous CASP experiment revealed progress in the category, particularly reflected in high prediction accuracy for larger protein domains. FM prediction models for two cases were sufficient to provide functional insights that were otherwise not obtainable by traditional sequence analysis methods. Importantly, CASP11 abstracts revealed that alignment-based contact prediction methods brought about much of the CASP11 progress, producing both of the functionally relevant models as well as several of the other outstanding structure predictions. These methodological advances enabled de novo modeling of much larger domain structures than was previously possible and allowed prediction of functional sites. Proteins 2016; 84(Suppl 1):51-66. © 2015 Wiley Periodicals, Inc.

Determining and visualizing flexibility in protein structures.

How to compare the structures of an ensemble of protein conformations is a fundamental problem in structural biology. As has been previously observed, the widely used RMSD measure due to Kabsch, in which a rigid-body superposition minimizing the least-squares positional deviations is performed, has its drawbacks when comparing and visualizing a set of flexible protein structures. Here, we develop a method, fleximatch, of protein structure comparison that takes flexibility into account. Based on a distance matrix measure of flexibility, a weighted superposition of distance matrices rather than of atomic coordinates is performed. Subsequently, this allows a consistent determination of (a) a superposition of structures for visualization, (b) a partitioning of the protein structure into rigid molecular components (core atoms), and (c) an atomic mobility measure. The method is suitable for highlighting both particularly flexible and rigid parts of a protein from structures derived from NMR, X-ray diffraction or molecular simulation.

Bh3 induced conformational changes in Bcl-Xl revealed by crystal structure and comparative analysis.

Apoptosis or programmed cell death is a regulatory process in cells in response to stimuli perturbing physiological conditions. The Bcl-2 family of proteins plays an important role in regulating homeostasis during apoptosis. In the process, the molecular interactions among the three members of this family, the pro-apoptotic, anti-apoptotic and BH3-only proteins at the mitochondrial outer membrane define the fate of a cell. Here, we report the crystal structures of the human anti-apoptotic protein Bcl-XL in complex with BH3-only BID(BH3) and BIM(BH3) peptides determined at 2.0 Å and 1.5 Å resolution, respectively. The BH3 peptides bind to the canonical hydrophobic pocket in Bcl-XL and adopt an alpha helical conformation in the bound form. Despite a similar structural fold, a comparison with other BH3 complexes revealed structural differences due to their sequence variations. In the Bcl-XL-BID(BH3) complex we observed a large pocket, in comparison with other BH3 complexes, lined by residues from helices α1, α2, α3, and α5 located adjacent to the canonical hydrophobic pocket. These results suggest that there are differences in the mode of interactions by the BH3 peptides that may translate into functional differences in apoptotic regulation.

Structure alignment of membrane proteins: Accuracy of available tools and a consensus strategy.

Protein structure alignment methods are used for the detection of evolutionary and functionally related positions in proteins. A wide array of different methods are available, but the choice of the best method is often not apparent to the user. Several studies have assessed the alignment accuracy and consistency of structure alignment methods, but none of these explicitly considered membrane proteins, which are important targets for drug development and have distinct structural features. Here, we compared 13 widely used pairwise structural alignment methods on a test set of homologous membrane protein structures (called HOMEP3). Each pair of structures was aligned and the corresponding sequence alignment was used to construct homology models. The model accuracy compared to the known structures was assessed using scoring functions not incorporated in the tested structural alignment methods. The analysis shows that fragment-based approaches such as FR-TM-align are the most useful for aligning structures of membrane proteins. Moreover, fragment-based approaches are more suitable for comparison of protein structures that have undergone large conformational changes. Nevertheless, no method was clearly superior to all other methods. Additionally, all methods lack a measure to rate the reliability of a position within a structure alignment. To solve both of these problems, we propose a consensus-type approach, combining alignments from four different methods, namely FR-TM-align, DaliLite, MATT, and FATCAT. Agreement between the methods is used to assign confidence values to each position of the alignment. Overall, we conclude that there remains scope for the improvement of structural alignment methods for membrane proteins.

Molecular mechanisms of protein kinase regulation by calcium/calmodulin.

Many human protein kinases are regulated by the calcium-sensor protein calmodulin, which binds to a short flexible segment C-terminal to the enzyme's catalytic kinase domain. Our understanding of the molecular mechanism of kinase activity regulation by calcium/calmodulin has been advanced by the structures of two protein kinases-calmodulin kinase II and death-associated protein kinase 1-bound to calcium/calmodulin. Comparison of these two structures reveals a surprising level of diversity in the overall kinase-calcium/calmodulin arrangement and functional readout of activity, as well as complementary mechanisms of kinase regulation such as phosphorylation.

Function-based assessment of structural similarity measurements using metal co-factor orientation.

Structure comparison is widely used to quantify protein relationships. Although there are several approaches to calculate structural similarity, specifying significance thresholds for similarity metrics is difficult due to the inherent likeness of common secondary structure elements. In this study, metal co-factor location is used to assess the biological relevance of structural alignments. The distance between the centroids of bound co-factors adds a chemical and function-relevant constraint to the structural superimposition of two proteins. This additional dimension can be used to define cut-off values for discriminating valid and spurious alignments in large alignment sets. The hypothesis underlying our approach is that metal coordination sites constrain structural evolution, thus revealing functional relationships between distantly related proteins. A comparison of three related nitrogenases shows the sequence and fold constraints imposed on the protein structures up to 18 Å away from the centers of their bound metal clusters.

A common structural scaffold in CTD phosphatases that supports distinct catalytic mechanisms.

The phosphorylation and dephosphorylation of the carboxyl-terminal domain (CTD) of the largest RNA polymerase II (RNAPII) subunit is a critical regulatory checkpoint for transcription and mRNA processing. This CTD is unique to eukaryotic organisms and it contains multiple tandem-repeats with the consensus sequence Tyr(1) -Ser(2) -Pro(3) -Thr(4) -Ser(5) -Pro(6) -Ser(7) . Traditionally, CTD phosphatases that use metal-ion-independent (cysteine-based) and metal-ion-assisted (aspartate-based) catalytic mechanisms have been considered to belong to two independent groups. However, using structural comparisons we have identified a common structural scaffold in these two groups of CTD phosphatases. This common scaffold accommodates different catalytic processes with the same substrate specificity, in this case phospho-serine/threonine residues flanked by prolines. Furthermore, this scaffold provides a structural connection between two groups of protein tyrosine phosphatases (PTPs): Cys-based (classes I, II, and III) and Asp-based (class IV) PTPs. Redundancy in catalytic mechanisms is not infrequent and may arise in specific biological settings. To better understand the activity of the CTD phosphatases, we combined our structural analyses with data on CTD phosphatase expression in different human and mouse tissues. The results suggest that aspartate- and cysteine-based CTD-dephosphorylation acts in concert during cellular stress, when high levels of reactive oxygen species can inhibit the nucleophilic function of the catalytic cysteine, as occurs in mental and neurodegenerative disorders like schizophrenia, Alzheimer's and Parkinson's diseases. Moreover, these findings have significant implications for the study of the RNAPII-CTD dephosphorylation in eukaryotes.

Rapid determination of RMSDs corresponding to macromolecular rigid body motions.

Finding the root mean sum of squared deviations (RMSDs) between two coordinate vectors that correspond to the rigid body motion of a macromolecule is an important problem in structural bioinformatics, computational chemistry, and molecular modeling. Standard algorithms compute the RMSD with time proportional to the number of atoms in the molecule. Here, we present RigidRMSD, a new algorithm that determines a set of RMSDs corresponding to a set of rigid body motions of a macromolecule in constant time with respect to the number of atoms in the molecule. Our algorithm is particularly useful for rigid body modeling applications, such as rigid body docking, and also for high-throughput analysis of rigid body modeling and simulation results. We also introduce a constant-time rotation RMSD as a similarity measure for rigid molecules. A C++ implementation of our algorithm is available at http://nano-d.inrialpes.fr/software/RigidRMSD.

A design of a nature-like fishway to solve the fractured river connectivity caused by small hydropower based on hydrodynamics and fish behaviors.

Achieving a harmonious alignment between the biological characteristics of fish and hydrodynamics patterns is crucial for ensuring the efficacy of fish passage facilities. In this study, based on the hydrodynamic characteristics of the river and the biological characteristics of fish, we evaluated the internal flow field in the nature-like fishway of Congen II hydropower station located along the Chabao river and explored methods to improve the operation efficiency. Based on comprehensive considerations of the flow field, turbulent kinetic energy, and the migration pathways of fish, it is found that the implementation of a continuous oblique bottom slope represents a more cost-effective and operationally convenient solution. The influence of different permutation of bulkheads in the nature-like fishway on operational efficiency was further examined. Our investigation revealed that the nature-like fishway with the continuous slope of 2% and the arrangement of three bulkheads in each row (model 3) exhibited a relatively simple velocity distribution and linear flow line, which poses challenges for fish in locating resting areas. In addition, the distribution of low turbulence kinetic energy area in the mainstream made it less favorable for fish to transition from the mainstream to the rest area within the fishway. The nature-like fishway with the continuous slope of 2% and the arrangement of two or three bulkheads in staggered rows (model 4) demonstrated better performance. Several potential fish migration routes for both model 3 and model 4 were proposed based on the numerical simulation results. In model 3, fish exhibited a continuous sprint through the concentrated high-speed area, which was less favorable for fish to rest and forage. In contrast, model 4 exhibited a diversified flow velocity distribution, enabling fish to make timely changes in their direction during migration. This feather proved to be advantageous in enhancing fish migration within the passage. The design of nature-like fishway in this study provides an important reference and technical support for the construction and optimization of the nature-like fishway for low dams, and is of great significance for restoring river connectivity destroyed by small hydropower construction and improving fish migration.