Search NASA⌕ Search

SEARCH · Search NASA

Results for “protein flexibility”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A bacterial surface layer protein exploits multistep crystallization for rapid self-assembly

Surface layers (S-layers) are crystalline protein coats surrounding microbial cells. S-layer proteins (SLPs) regulate their extracellular self-assembly by crystallizing when exposed to an environmental trigger. However, molecular mechanisms governing rapid protein crystallization in vivo or in vitro are largely unknown. Here, we demonstrate that the Caulobacter crescentus SLP readily crystallizes into sheets in vitro via a calcium-triggered multistep assembly pathway. This pathway involves 2 domains serving distinct functions in assembly. The C-terminal crystallization domain forms the physiological 2-dimensional (2D) crystal lattice, but full-length protein crystallizes multiple orders of magnitude faster due to the N-terminal nucleation domain. Observing crystallization using a time course of electron cryo-microscopy (Cryo-EM) imaging reveals a crystalline intermediate wherein N-terminal nucleation domains exhibit motional dynamics with respect to rigid lattice-forming crystallization domains. Dynamic flexibility between the 2 domains rationalizes efficient S-layer crystal nucleation on the curved cellular surface. Rate enhancement of protein crystallization by a discrete nucleation domain may enable engineering of kinetically controllable self-assembling 2D macromolecular nanomaterials.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The Epidermal Microbiome Within an Aggregation of Leopard Sharks ( Triakis semifasciata ) Has Taxonomic Flexibility with Gene Functional Stability Across Three Time-points

The epidermis of Chondrichthyan fishes consists of dermal denticles with production of minimal but protein-rich mucus that collectively, influence the attachment and biofilm development of microbes, facilitating a unique epidermal microbiome. Here, we use metagenomics to provide the taxonomic and functional characterization of the epidermal microbiome of the Triakis semifasciata (leopard shark) at three time-points collected across 4 years to identify links between microbial groups and host metabolism. Our aims include (1) describing the variation of microbiome taxa over time and identifying recurrent microbiome members (present across all time-points); (2) investigating the relationship between the recurrent and flexible taxa (those which are not found consistently across time-points); (3) describing the functional compositions of the microbiome which may suggest links with the host metabolism; and (4) identifying whether metabolic processes are shared across microbial genera or are unique to specific taxa. Microbial members of the microbiome showed high similarity between all individuals (Bray–Curtis similarity index = 82.7, where 0 = no overlap, 100 = total overlap) with the relative abundance of those members varying across sampling time-points, suggesting flexibility of taxa in the microbiome. One hundred and eighty-eight genera were identified as recurrent, including Pseudomonas, Erythrobacter, Alcanivorax, Marinobacter, and Sphingopxis being consistently abundant across time-points, while Limnobacter and Xyella exhibited switching patterns with high relative abundance in 2013, Sphingobium and Sphingomona in 2015, and Altermonas, Leeuwenhoekiella, Gramella, and Maribacter in 2017. Of the 188 genera identified as recurrent, the top 19 relatively abundant genera formed three recurrent groups. The microbiome also displayed high functional similarity between individuals (Bray–Curtis similarity index = 97.6) with gene function composition remaining consistent across all time-points. These results show that while the presence of microbial genera exhibits consistency across time-points, their abundances do fluctuate. Microbial functions however remain stable across time-points; thus, we suggest the leopard shark microbiomes exhibit functional redundancy. We show coexistence of microbes hosted in elasmobranch microbiomes that encode genes involved in utilizing nitrogen, but not fixing nitrogen, degrading urea, and resistant to heavy metal.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Domain architecture and catalysis of the Staphylococcus aureus fatty acid kinase

Fatty acid kinase (Fak) is a two-component enzyme that generates acyl-phosphate for phospholipid synthesis. Fak consists of a kinase domain protein (FakA) that phosphorylates a fatty acid enveloped by a fatty acid binding protein (FakB). The structural basis for FakB function has been established, but little is known about FakA. Here, we used limited proteolysis to define three separate FakA domains: the amino terminal FakA_N, the central FakA_L, and the carboxy terminal FakA_C. The isolated domains lack kinase activity, but activity is restored when FakA_N and FakA_L are present individually or connected as FakA_NL. The X-ray structure of the monomeric FakA_N captures the product complex with ADP and two Mg 2+ ions bound at the nucleotide site. The FakA_L domain encodes the dimerization interface along with conserved catalytic residues Cys240, His282, and His284. AlphaFold analysis of FakA_L predicts the catalytic residues are spatially clustered and pointing away from the dimerization surface. Furthermore, the X-ray structure of FakA_C shows that it consists of two subdomains that are structurally related to FakB. Analytical ultracentrifugation demonstrates that FakA_C binds FakB, and site-directed mutagenesis confirms that a positively charged wedge on FakB meshes with a negatively charged groove on FakA_C. Finally, small angle X-ray scattering analysis is consistent with freely rotating FakA_N and FakA_C domains tethered by flexible linkers to FakA_L. These data reveal specific roles for the three independently folded FakA protein domains in substrate binding and catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

In Situ Proteolysis Condition-Induced Crystallization of the XcpVWX Complex in Different Lattices

Although prevalent in the determination of protein structures; crystallography always has the bottleneck of obtaining high-quality protein crystals for characterizing a wide range of proteins; especially large protein complexes. Stable fragments or domains of proteins are more readily to crystallize; which prompts the use of in situ proteolysis to remove flexible or unstable structures for improving crystallization and crystal quality. In this work; we investigated the effects of in situ proteolysis by chymotrypsin on the crystallization of the XcpVWX complex from the Type II secretion system of Pseudomonas aeruginosa. Different proteolysis conditions were found to result in two distinct lattices in the same crystallization solution. With a shorter chymotrypsin digestion at a lower concentration; the crystals exhibited a P3 hexagonal lattice that accommodates three complex molecules in one asymmetric unit. By contrast; a longer digestion with chymotrypsin of a 10-fold higher concentration facilitated the formation of a compact P212121 orthorhombic lattice with only one complex molecule in each asymmetric unit. The molecules in the hexagonal lattice have shown high atomic displacement parameter values compared with the ones in the orthorhombic lattice. Taken together; our results clearly demonstrate that different proteolysis conditions can result in the generation of distinct lattices in the same crystallization solution; which can be exploited in order to obtain different crystal forms of a better quality

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

In situ ligand restraints from quantum-mechanical methods

In macromolecular crystallographic structure refinement, ligands present challenges for the generation of geometric restraints due to their large chemical variability, their possible novel nature and their specific interaction with the binding pocket of the protein. Quantum-mechanical approaches are useful for providing accurate ligand geometries, but can be plagued by the number of minima in flexible molecules. In an effort to avoid these issues, the Quantum Mechanical Restraints (QMR) procedure optimizes the ligand geometry in situ , thus accounting for the influence of the macromolecule on the local energy minima of the ligand. The optimized ligand geometry is used to generate target values for geometric restraints during the crystallographic refinement. As demonstrated using a sample of >2330 ligand instances in >1700 protein–ligand models, QMR restraints generally result in lower deviations from the target stereochemistry compared with conventionally generated restraints. In particular, the QMR approach provides accurate torsion restraints for ligands and other entities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

X‐ray scattering reveals disordered linkers and dynamic interfaces in complexes and mechanisms for DNA double‐strand break repair impacting cell and cancer biology

Abstract Evolutionary selection ensures specificity and efficiency in dynamic metastable macromolecular machines that repair DNA damage without releasing toxic and mutagenic intermediates. Here we examine non‐homologous end joining (NHEJ) as the primary conserved DNA double‐strand break (DSB) repair process in human cells. NHEJ has exemplary key roles in networks determining the development, outcome of cancer treatments by DSB‐inducing agents, generation of antibody and T‐cell receptor diversity, and innate immune response for RNA viruses. We determine mechanistic insights into NHEJ structural biochemistry focusing upon advanced small angle X‐ray scattering (SAXS) results combined with X‐ray crystallography (MX) and cryo‐electron microscopy (cryo‐EM). SAXS coupled to atomic structures enables integrated structural biology for objective quantitative assessment of conformational ensembles and assemblies in solution, intra‐molecular distances, structural similarity, functional disorder, conformational switching, and flexibility. Importantly, NHEJ complexes in solution undergo larger allosteric transitions than seen in their cryo‐EM or MX structures. In the long‐range synaptic complex, X‐ray repair cross‐complementing 4 (XRCC4) plus XRCC4‐like‐factor (XLF) form a flexible bridge and linchpin for DNA ends bound to KU heterodimer (Ku70/80) and DNA‐PKcs (DNA‐dependent protein kinase catalytic subunit). Upon binding two DNA ends, auto‐phosphorylation opens DNA‐PKcs dimer licensing NHEJ via concerted conformational transformations of XLF‐XRCC4, XLF–Ku80, and LigIV BRCT –Ku70 interfaces. Integrated structures reveal multifunctional roles for disordered linkers and modular dynamic interfaces promoting DSB end processing and alignment into the short‐range complex for ligation by LigIV. Integrated findings define dynamic assemblies fundamental to designing separation‐of‐function mutants and allosteric inhibitors targeting conformational transitions in multifunctional complexes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The yeast ALA synthase C‐terminus positively controls enzyme structure and function

Abstract 5‐Aminolevulinic acid synthase (ALAS) is a pyridoxal 5′‐phosphate (PLP)‐dependent enzyme that catalyzes the first and rate‐limiting step of heme biosynthesis in α‐proteobacteria and several non‐plant eukaryotes. All ALAS homologs contain a highly conserved catalytic core, but eukaryotes also have a unique C‐terminal extension that plays a role in enzyme regulation. Several mutations in this region are implicated in multiple blood disorders in humans. In Saccharomyces cerevisiae ALAS (Hem1), the C‐terminal extension wraps around the homodimer core to contact conserved ALAS motifs proximal to the opposite active site. To determine the importance of these Hem1 C‐terminal interactions, we determined the crystal structure of S. cerevisiae Hem1 lacking the terminal 14 amino acids (Hem1 ΔCT). With truncation of the C‐terminal extension, we show structurally and biochemically that multiple catalytic motifs become flexible, including an antiparallel β‐sheet important to Fold‐Type I PLP‐dependent enzymes. The changes in protein conformation result in an altered cofactor microenvironment, decreased enzyme activity and catalytic efficiency, and ablation of subunit cooperativity. These findings suggest that the eukaryotic ALAS C‐terminus has a homolog‐specific role in mediating heme biosynthesis, indicating a mechanism for autoregulation that can be exploited to allosterically modulate heme biosynthesis in different organisms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Conformational motions and ligand-binding underlying gating and regulation in IP 3 R channel

Inositol-1,4,5-trisphosphate receptors (IP3Rs) are activated by IP 3 and Ca 2+ and their gating is regulated by various intracellular messengers that finely tune the channel activity. Here, using single particle cryo-EM analysis we determined 3D structures of the nanodisc-reconstituted IP 3 R1 channel in two ligand-bound states. These structures provide unprecedented details governing binding of IP 3 , Ca 2+ and ATP, revealing conformational changes that couple ligand-binding to channel opening. Using a deep-learning approach and 3D variability analysis we extracted molecular motions of the key protein domains from cryo-EM density data. We find that IP 3 binding relies upon intrinsic flexibility of the ARM2 domain in the tetrameric channel. Our results highlight a key role of dynamic side chains in regulating gating behavior of IP 3 R channels. This work represents a stepping-stone to developing mechanistic understanding of conformational pathways underlying ligand-binding, activation and regulation of the channel.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

The crystal structure of methanogen McrD , a methyl‐coenzyme M reductase‐associated protein

Methyl‐coenzyme M reductase (MCR) is a multi‐subunit (α 2 β 2 γ 2 ) enzyme responsible for methane formation via its unique F 430 cofactor. The genes responsible for producing MCR ( mcrA , mcrB and mcrG ) are typically colocated with two other highly conserved genes mcrC and mcrD . We present here the high‐resolution crystal structure for McrD from a human gut methanogen Methanomassiliicoccus luminyensis strain B10. The structure reveals that McrD comprises a ferredoxin‐like domain assembled into an α + β barrel‐like dimer with conformational flexibility exhibited by a functional loop. The description of the M. luminyensis McrD crystal structure contributes to our understanding of this key conserved methanogen protein typically responsible for promoting MCR activity and the production of methane, a greenhouse gas.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Rational design of photosynthetic reaction center protein maquettes

New technologies for efficient solar-to-fuel energy conversion will help facilitate a global shift from dependence on fossil fuels to renewable energy. Nature uses photosynthetic reaction centers to convert photon energy into a cascade of electron-transfer reactions that eventually produce chemical fuel. The design of new reaction centers de novo deepens our understanding of photosynthetic charge separation and may one day allow production of biofuels with higher thermodynamic efficiency than natural photosystems. Recently, we described the multi-step electron-transfer activity of a designed reaction center maquette protein (the RC maquette), which can assemble metal ions, tyrosine, a Zn tetrapyrrole, and heme into an electron-transport chain. Here, we detail our modular strategy for rational protein design and show that the intended RC maquette design agrees with crystal structures in various states of assembly. A flexible, dynamic apo-state collapses by design into a more ordered holo-state upon cofactor binding. Crystal structures illustrate the structural transitions upon binding of different cofactors. Spectroscopic assays demonstrate that the RC maquette binds various electron donors, pigments, and electron acceptors with high affinity. We close with a critique of the present RC maquette design and use electron-tunneling theory to envision a path toward a designed RC with a substantially higher thermodynamic efficiency than natural photosystems.

09 BIOMASS FUELS↗
DOE OSTI2024

ZMPY3D: accelerating protein structure volume analysis through vectorized 3D Zernike moments and Python-based GPU integration

Abstract Motivation Volumetric 3D object analyses are being applied in research fields such as structural bioinformatics, biophysics, and structural biology, with potential integration of artificial intelligence/machine learning (AI/ML) techniques. One such method, 3D Zernike moments, has proven valuable in analyzing protein structures (e.g., protein fold classification, protein–protein interaction analysis, and molecular dynamics simulations). Their compactness and efficiency make them amenable to large-scale analyses. Established methods for deriving 3D Zernike moments, however, can be inefficient, particularly when higher order terms are required, hindering broader applications. As the volume of experimental and computationally-predicted protein structure information continues to increase, structural biology has become a “big data” science requiring more efficient analysis tools. Results This application note presents a Python-based software package, ZMPY3D, to accelerate computation of 3D Zernike moments by vectorizing the mathematical formulae and using graphical processing units (GPUs). The package offers popular GPU-supported libraries such as CuPy and TensorFlow together with NumPy implementations, aiming to improve computational efficiency, adaptability, and flexibility in future algorithm development. The ZMPY3D package can be installed via PyPI, and the source code is available from GitHub. Volumetric-based protein 3D structural similarity scores and transform matrix of superposition functionalities have both been implemented, creating a powerful computational tool that will allow the research community to amalgamate 3D Zernike moments with existing AI/ML tools, to advance research and education in protein structure bioinformatics. Availability and implementation ZMPY3D, implemented in Python, is available on GitHub (https://github.com/tawssie/ZMPY3D) and PyPI, released under the GPL License.

Lai, Jhih-Siang (ORCID:0000000156775890)↗
DOE OSTI2024

Adjusting the Energy Profile for CH–O Interactions Leads to Improved Stability of RNA Stem-Loop Structures in MD Simulations

The role of ribonucleic acid (RNA) in biology continues to grow, but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Further many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, the role of weak CH···O interactions is explored, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, it is shown that CH···O interactions are overly repulsive in the widely used Amber RNA force fields. A simple, targeted adjustment of CH···O repulsion that leaves the remainder of the force field unchanged was developed. Then, the standard and modified force fields were tested using molecular dynamics (MD) simulations with explicit water and salt, amassing over 300 μs of data for multiple RNA systems containing important features such as the presence of loops, base stacking interactions as well as canonical and noncanonical base pairing. In this work and others, standard force fields lead to reproducible unfolding of the NMR-based structures. Including a targeted CH···O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Structural basis of the amidase ClbL central to the biosynthesis of the genotoxin colibactin

Colibactin is a genotoxic natural product produced by select commensal bacteria in the human gut microbiota. The compound is a bis-electrophile that is predicted to form interstrand DNA cross-links in target cells, leading to double-strand DNA breaks. The biosynthesis of colibactin is carried out by a mixed NRPS–PKS assembly line with several noncanonical features. An amidase, ClbL, plays a key role in the pathway, catalyzing the final step in the formation of the pseudodimeric scaffold. ClbL couples α-aminoketone and β-ketothioester intermediates attached to separate carrier domains on the NRPS–PKS assembly. Here, the 1.9 Å resolution structure of ClbL is reported, providing a structural basis for this key step in the colibactin biosynthetic pathway. The structure reveals an open hydrophobic active site surrounded by flexible loops, and comparison with homologous amidases supports its unusual function and predicts macromolecular interactions with pathway carrier-protein substrates. Modeling protein–protein interactions supports a predicted molecular basis for enzyme–carrier domain interactions. Overall, the work provides structural insight into this unique enzyme that is central to the biosynthesis of colibactin.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

RCSB Protein Data Bank 1D tools and services

Abstract Motivation Interoperability between polymer sequences and structural data is essential for providing a complete picture of protein and gene features and helping to understand biomolecular function. Results Herein, we present two resources designed to improve interoperability between the RCSB Protein Data Bank, the NCBI and the UniProtKB data resources and visualize integrated data therefrom. The underlying tools provide a flexible means of mapping between the different coordinate spaces and an interactive tool allows convenient visualization of the 1-dimensional data over the web. Availabilityand implementation https://1d-coordinates.rcsb.org and https://rcsb.github.io/rcsb-saguaro. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Crystal structure of Methanococcus jannaschii dihydroorotase with substrate bound

Here, we report the X-ray structural analysis of dihydroorotase from Methanococcus jannaschii co-crystallized with dihydroorotate at pH 6.5. The crystals are isomorphous with the crystals of the apoenzyme, with space group P3 2 21 and unit-cell dimensions a = b = 111.4, c = 101.2 Å. The structure was refined to R = 0.169 and R free = 0.186 at a resolution of 1.87 Å at room temperature. During crystallization the degradative reaction took place, and the electron density in the active site corresponds to the substrate carbamoyl aspartate. Carbamoyl aspartate interacts with the protein in the active site in a manner similar to that observed for Escherichia coli and human dihydroorotases. However, the flexible loop (residues 140–151) adopts both conformations in the crystal, loop-out and loop-in, with the loop-out conformation having higher occupancy. This contrasts with our expectations for the flexible loop to be exclusively in the loop-in conformation, which is the conformation that it adopts to stabilize the binding of the substrate in the known systems. Additional studies with substrate analogs that resemble carbamoyl aspartate and different crystallization conditions would provide further insight into the conformation of the flexible loop in this system.

Methanococcus jannaschii↗
DOE OSTI2021

OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations

We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

tinyIFD: A High-Throughput Binding Pose Refinement Workflow Through Induced-Fit Ligand Docking

A critical step in structure-based drug discovery is predicting whether and how a candidate molecule binds to a model of a therapeutic target. However, substantial protein side chain movements prevent current screening methods, such as docking, from accurately predicting the ligand conformations and require expensive refinements to produce viable candidates. Here, we present the development of a high-throughput and flexible ligand pose refinement workflow, called “tinyIFD”. The main features of the workflow include the use of specialized high-throughput, small-system MD simulation code mdgx.cuda and an actively learning model zoo approach. We show the application of this workflow on a large test set of diverse protein targets, achieving 66% and 76% success rates for finding a crystal-like pose within the top-2 and top-5 poses, respectively. We also applied this workflow to the SARS-CoV-2 main protease (M pro ) inhibitors, where we demonstrate the benefit of the active learning aspect in this workflow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Heterogeneous flexibility can contribute to chromatin segregation in the cell nucleus

The highly and slightly condensed forms of chromatin, heterochromatin and euchromatin, respectively, segregate in the cell nucleus. Heterochromatin is more abundant in the nucleus periphery. Here, in this work, we study the mechanism of heterochromatin segregation by modeling interphase chromosomes as diblock ring copolymers confined in a rigid spherical shell using molecular dynamics simulations. In our model, heterochromatin and euchromatin are distinguished by their bending stiffnesses only, while an interaction potential between the spherical shell and chromatin is used to model lamin-associated proteins. Our simulations indicate that in the absence of attractive interactions between the nuclear shell and the chromatin, most heterochromatin segregates towards the nuclear interior due to the depletion of less flexible heterochromatin segments from the nuclear periphery. This inverted chromatin distribution, which is opposite to the conventional case with heterochromatin dominating at the periphery, is in accord with experimental observations in rod cells. This “inversion” is also found to be independent of the heterochromatin concentration and chromosome number. The chromatin distribution at the periphery found in vivo can be recovered by further increasing the bending stiffness of heterochromatin segments or by turning on attractive interactions between the nuclear shell and heterochromatin. Our results indicate that the bending stiffness of chromatin could be a contributor to chromosome organization along with differential effects of HP ⁢1⁢𝛼-driven phase segregation and of loop extruders and interactions with the nuclear envelope and topological constraints.

biomolecular & subcellular processes↗