Unveiling the N-Terminal Homodimerization of BCL11B by Hybrid Solvent Replica-Exchange Simulations

UMSBP2 is chromatin remodeler that functions in regulation of gene expression and suppression of antigenic variation in trypanosomes

Universal Minicircle Sequence binding proteins (UMSBPs) are CCHC-type zinc-finger proteins that bind the single-stranded G-rich UMS sequence, conserved at the replication origins of minicircles in the kinetoplast DNA, the mitochondrial genome of kinetoplastids. Trypanosoma brucei UMSBP2 has been recently shown to colocalize with telomeres and to play an essential role in chromosome end protection. Here we report that TbUMSBP2 decondenses in vitro DNA molecules, which were condensed by core histones H2B, H4 or linker histone H1. DNA decondensation is mediated via protein-protein interactions between TbUMSBP2 and these histones, independently of its previously described DNA binding activity. Silencing of the TbUMSBP2 gene resulted in a significant decrease in the disassembly of nucleosomes in T. brucei chromatin, a phenotype that could be reverted, by supplementing the knockdown cells with TbUMSBP2. Transcriptome analysis revealed that silencing of TbUMSBP2 affects the expression of multiple genes in T. brucei, with a most significant effect on the upregulation of the subtelomeric variant surface glycoproteins (VSG) genes, which mediate the antigenic variation in African trypanosomes. These observations suggest that UMSBP2 is a chromatin remodeling protein that functions in the regulation of gene expression and plays a role in the control of antigenic variation in T. brucei.

Fundamental Redesign of the TIGER2hs Kernel to Address Severe Parameter Sensitivity

Replica exchange molecular dynamics simulations are one of the most popular approaches to enhance conformational sampling of molecular systems. Applications range from protein folding to protein-protein or other host-guest interactions, as well as binding free energy calculations. While these methods are computationally expensive, highly accurate results can be obtained. We recently developed TIGER2hs, an improved version of the temperature intervals with global exchange of replicas (TIGER2) algorithm. This method combines the replica-based enhanced sampling in an explicit solvent with a hybrid solvent energy evaluation. During the exchange attempts, bulk water is replaced by an implicit solvent model, allowing sampling with significantly less replicas than parallel tempering (REMD). This enables accurate enhanced sampling calculations with only a fraction of computational resources compared to REMD. Our latest results highlight several issues with sampling imbalance and parameter sensitivity within the original TIGER2 exchange algorithms that affect the overall state populations. A high sensitivity on replica number and maximum temperature is eliminated by changing to a pairwise exchange kernel (PE) without additional sorting. Simulations are controlled by adjusting the average temperature change per exchange ⟨ΔT/χ⟩ to below 30 K to mimic a controlled temperature mixing of replicas similar to REMD. Thus, this parameter provides an applicable property for selecting combinations of replica number and maximum temperature to adjust simulations for best accuracy, with flexible resource investment. This increases the robustness of the method and ensures results in excellent agreement with REMD, as demonstrated for three different peptides.