Biotite: new tools for a versatile Python bioinformatics library

Sequence modeling and design from molecular to genome scale with Evo

The genome is a sequence that encodes the DNA, RNA, and proteins that orchestrate an organism's function. We present Evo, a long-context genomic foundation model with a frontier architecture trained on millions of prokaryotic and phage genomes, and report scaling laws on DNA to complement observations in language and vision. Evo generalizes across DNA, RNA, and proteins, enabling zero-shot function prediction competitive with domain-specific language models and the generation of functional CRISPR-Cas and transposon systems, representing the first examples of protein-RNA and protein-DNA codesign with a language model. Evo also learns how small mutations affect whole-organism fitness and generates megabase-scale sequences with plausible genomic architecture. These prediction and generation capabilities span molecular to genomic scales of complexity, advancing our understanding and control of biology.

Drude2019IDPC polarizable force field reveals structure-function relationship of insulin

Intrinsically disordered proteins (IDPs) lack stable tertiary structures under physiological conditions, yet play key roles in biological processes and associated with human complex diseases. Their conformational characteristics and high content of charged residues make the use of polarizable force fields an advantageous for simulating IDPs. The Drude2019IDP polarizable force field, previously introduced, has demonstrated comprehensive enhancements and improvements in dipeptides, short peptides, and IDPs, achieving a balanced sampling between IDPs and structured proteins. However, the performance in simulating 5 dipeptides was found to be underestimate. Therefore, we individually performed reweighting and grid-based energy correction map (CMAP) optimization for these 5 dipeptides, resulting in the enhanced Drude2019IDPC force field. The performance of Drude2019IDPC was evaluated with 5 dipeptides, 5 disordered short peptides, and a representative IDP. The results demonstrated a marked improvement comparing with original Drude2019IDP. To further substantiate the capabilities of Drude2019IDPC, MD simulation and Markov state model (MSM) were applied to wild type and mutant for insulin, to elucidate the difference of conformational characteristics and transition path. The findings reveal that mutation can maintain the monomorphic characteristics, providing insights for engineered insulin development. These results indicate that Drude2019IDPC could be used to reveal the structure-function relationship for other proteins.

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

Anisotropic Network Analysis of Open/Closed HCN4 Channel Advocates Asymmetric Subunit Cooperativity in cAMP Modulation of Gating

Hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are opened in an allosteric manner by membrane hyperpolarization and cyclic nucleotides such as cAMP. Because of conflicting reports from experimental studies on whether cAMP binding to the four available binding sites in the channel tetramer operates cooperatively in gating, we employ here a computational approach as a promising route to examine ligand-induced conformational changes after binding to individual sites. By combining an elastic network model (ENM) with linear response theory (LRT) for modeling the apo-holo transition of the cyclic nucleotide-binding domain (CNBD) in HCN channels, we observe a distinct pattern of cooperativity matching the "positive-negative-positive" cooperativity reported from functional studies. This cooperativity pattern is highly conserved among HCN subtypes (HCN4, HCN1), but only to a lesser extent visible in structurally related channels, which are only gated by voltage (KAT1) or cyclic nucleotides (TAX4). This suggests an inherent cooperativity between subunits in HCN channels as part of a ligand-triggered gating mechanism in these channels.

Linking Gene Fusions to Bone Marrow Failure and Malignant Transformation in Dyskeratosis Congenita

Dyskeratosis Congenita (DC) is a multisystem disorder intrinsically associated with telomere dysfunction, leading to bone marrow failure (BMF). Although the pathology of DC is largely driven by mutations in telomere-associated genes, the implications of gene fusions, which emerge due to telomere-induced genomic instability, remain unexplored. We meticulously analyzed gene fusions in RNA-Seq data from DC patients to provide deeper insights into DC's progression. The most significant DC-specific gene fusions were subsequently put through in silico assessments to ascertain biophysical and structural attributes, including charge patterning, inherent disorder, and propensity for self-association. Selected candidates were then analyzed using deep learning-powered structural predictions and molecular dynamics simulations to gauge their potential for forming higher-order oligomers. Our exploration revealed that genes participating in fusion events play crucial roles in upholding genomic stability, facilitating hematopoiesis, and suppressing tumors. Notably, our analysis spotlighted a particularly disordered polyampholyte fusion protein that exhibits robust higher-order oligomerization dynamics. To conclude, this research underscores the potential significance of several high-confidence gene fusions in the progression of BMF in DC, particularly through the dysregulation of genomic stability, hematopoiesis, and tumor suppression. Additionally, we propose that these fusion proteins might hold a detrimental role, specifically in inducing proteotoxicity-driven hematopoietic disruptions.

SBILib: a handle for protein modeling and engineering

SUMMARY

The SBILib Python library provides an integrated platform for the analysis of macromolecular structures and interactions. It combines simple 3D file parsing and workup methods with more advanced analytical tools. SBILib includes modules for macromolecular interactions, loops, super-secondary structures, and biological sequences, as well as wrappers for external tools with which to integrate their results and facilitate the comparative analysis of protein structures and their complexes. The library can handle macromolecular complexes formed by proteins and/or nucleic acid molecules (i.e. DNA and RNA). It is uniquely capable of parsing and calculating protein super-secondary structure and loop geometry. We have compiled a list of example scenarios which SBILib may be applied to and provided access to these within the library.

AVAILABILITY AND IMPLEMENTATION

SBILib is made available on Github at https://github.com/structuralbioinformatics/SBILib.