Distinguishing Biomolecular Pathways and Metastable States

Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample

Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.

Understanding the Energy Landscape of Intrinsically Disordered Protein Ensembles

A substantial portion of various organisms' proteomes comprises intrinsically disordered proteins (IDPs) that lack a defined three-dimensional structure. These IDPs exhibit a diverse array of conformations, displaying remarkable spatiotemporal heterogeneity and exceptional conformational flexibility. Characterizing the structure or structural ensemble of IDPs presents significant conceptual and methodological challenges owing to the absence of a well-defined native structure. While databases such as the Protein Ensemble Database (PED) provide IDP ensembles obtained through a combination of experimental data and molecular modeling, the absence of reaction coordinates poses challenges in comprehensively understanding pertinent aspects of the system. In this study, we leverage the energy landscape visualization method (JCTC, 6482, 2019) to scrutinize four IDP ensembles sourced from PED. ELViM, a methodology that circumvents the need for a priori reaction coordinates, aids in analyzing the ensembles. The specific IDP ensembles investigated are as follows: two fragments of nucleoporin (NUL: 884-993 and NUS: 1313-1390), yeast sic 1 N-terminal (1-90), and the N-terminal SH3 domain of Drk (1-59). Utilizing ELViM enables the comprehensive validation of ensembles, facilitating the detection of potential inconsistencies in the sampling process. Additionally, it allows for identifying and characterizing the most prevalent conformations within an ensemble. Moreover, ELViM facilitates the comparative analysis of ensembles obtained under diverse conditions, thereby providing a powerful tool for investigating the functional mechanisms of IDPs.

Helical Content Correlations and Hydration Structures of the Folding Ensemble of the B Domain of Protein A

The B domain of protein A (BdpA), a small three-helix bundle, folds on a time scale of a few microseconds with heterogeneous native and unfolded states. It is widely used as a model for understanding protein folding mechanisms. In this work, we use structure-based models (SBMs) and atomistic simulations to comprehensively investigate how BdpA folding is associated with the formation of its secondary structure. The energy landscape visualization method (ELViM) was used to characterize the pathways that connect the folded and unfolded states of BdpA as well as the sets of structures displaying specific ellipticity patterns. We show that the native state conformational diversity is due mainly to the conformational variability of helix I. Helices I, II, and III occur in a weakly correlated manner, with Spearman's rank correlation coefficients of 0.1539 (I and II), 0.1259 (I and III), and 0.2561 (II and III). These results, therefore, suggest the highest cooperativity between helices II and III. Our results allow the clustering of partially folded structures of folding of the B domain of protein A on the basis of its secondary structure, paving the way to an understanding of environmental factors in the relative stability of the basins of the folding ensemble, which are illustrated by the structural dependency of the protein hydration structures, as computed with minimum-distance distribution functions.

ELViM: Exploring Biomolecular Energy Landscapes through Multidimensional Visualization

Molecular dynamics (MD) simulations provide a powerful means of exploring the dynamic behavior of biomolecular systems at the atomic level. However, analyzing the vast data sets generated by MD simulations poses significant challenges. This article discusses the energy landscape visualization method (ELViM), a multidimensional reduction technique inspired by the energy landscape theory. ELViM transcends one-dimensional representations, offering a comprehensive analysis of the effective conformational phase space without the need for predefined reaction coordinates. We apply the ELViM to study the folding landscape of the antimicrobial peptide Polybia-MP1, showcasing its versatility in capturing complex biomolecular dynamics. Using dissimilarity matrices and a force-scheme approach, the ELViM provides intuitive visualizations, revealing structural correlations and local conformational signatures. The method is demonstrated to be adaptable, robust, and applicable to various biomolecular systems.

Exploring the folding landscape of leptin: Insights into threading pathways

The discovery of new protein topologies with entanglements and loop-crossings have shown the impact of local amino acid arrangement and global three-dimensional structures. This phenomenon plays a crucial role in understanding how protein structure relates to folding and function, affecting the global stability, and biological activity. Protein entanglements encompassing knots and non-trivial topologies add complexity to their folding free energy landscapes. However, the initial native contacts driving the threading event for entangled proteins remains elusive. The Pierced Lasso Topology (PLT) represents an entangled topology where a covalent linker creates a loop in which the polypeptide backbone is threaded through. Compared to true knotted topologies, PLTs are simpler topologies where the covalent-loop persists in all conformations. In this work, the PLT protein leptin, is used to visualize and differentiate the preference for slipknotting over plugging transition pathways along the folding route. We utilize the Energy Landscape Visualization Method (ELViM), a multidimensional projection technique, to visualize and distinguish early threaded conformations that cannot be observed in an in vitro experiment. Critical contacts for the leptin threading mechanisms were identified where the competing pathways are determined by the formation of a hairpin loop in the unfolded basin. Thus, prohibiting the dominant slipknotting pathway. Furthermore, ELViM offers insights into distinct folding pathways associated with slipknotting and plugging providing a novel tool for de novo design and in vitro experiments with residue specific information of threading events in silico.

Probing Mastoparan-like Antimicrobial Peptides Interaction with Model Membrane Through Energy Landscape Analysis

Antimicrobial Peptides (AMPs) have emerged as promising alternatives to conventional antibiotics due to their capacity to disrupt the lipid packing of bacterial cell membranes. This mechanism of action may prevent the development of resistance by bacteria. Understanding their role in lipid packing disruption and their structural properties upon interaction with bacterial membranes is highly desirable. In this study, we employed Molecular Dynamics simulations and the Energy Landscape Visualization Method (ELViM) to characterize and compare the conformational ensembles of mastoparan-like Polybia-MP1 and its analogous H-MP1, in which histidines replace lysine residues. Two situations were analyzed: (i) the peptides in their free state in an aqueous solution containing water and ions and (ii) the peptides spontaneously adsorbing onto an anionic lipid bilayer, used as a bacteria membrane mimetic. ELViM was used to project a single effective conformational phase space for both peptides, providing a comparative analysis. This projection enabled us to map the conformational ensembles of each peptide in an aqueous solution and assess the structural effects of substituting lysines with histidines in H-MP1. Furthermore, a single conformational phase space analysis was employed to describe structural changes during the adsorption process using the same framework. We show that ELViM provides a comprehensive analysis, able to identify discrepancies in the conformational ensembles of these peptides that may affect their affinity to the membrane and adsorption kinetics.

Unveiling Metastable Ensembles of GRB2 and the Relevance of Interdomain Communication during Folding

The folding process of multidomain proteins is a highly intricate phenomenon involving the assembly of distinct domains into a functional three-dimensional structure. During this process, each domain may fold independently while interacting with others. The folding of multidomain proteins can be influenced by various factors, including their composition, the structure of each domain, or the presence of disordered regions, as well as the surrounding environment. Misfolding of multidomain proteins can lead to the formation of nonfunctional structures associated with a range of diseases, including cancers or neurodegenerative disorders. Understanding this process is an important step for many biophysical analyses such as stability, interaction, malfunctioning, and rational drug design. One such multidomain protein is growth factor receptor-bound protein 2 (GRB2), an adaptor protein that is essential in regulating cell survival. GRB2 consists of one central Src homology 2 (SH2) domain flanked by two Src homology 3 (SH3) domains. The SH2 domain interacts with phosphotyrosine regions in other proteins, while the SH3 domains recognize proline-rich regions on protein partners during cell signaling. Here, we combined computational and experimental techniques to investigate the folding process of GRB2. Through computational simulations, we sampled the conformational space and mapped the mechanisms involved by the free energy profiles, which may indicate possible intermediate states. From the molecular dynamics trajectories, we used the energy landscape visualization method (ELViM), which allowed us to visualize a three-dimensional (3D) representation of the overall energy surface. We identified two possible parallel folding routes that cannot be seen in a one-dimensional analysis, with one occurring more frequently during folding. Supporting these results, we used differential scanning calorimetry (DSC) and fluorescence spectroscopy techniques to confirm these intermediate states in vitro. Finally, we analyzed the deletion of domains to compare our model outputs to previously published results, supporting the presence of interdomain modulation. Overall, our study highlights the significance of interdomain communication within the GRB2 protein and its impact on the formation, stability, and structural plasticity of the protein, which are crucial for its interaction with other proteins in key signaling pathways.

Characterizing the Folding Transition-State Ensembles in the Energy Landscape of an RNA Tetraloop

Molecular dynamics (MD) simulations have become increasingly powerful and can now describe the folding/unfolding of small biomolecules in atomic detail. However, a major challenge in MD simulations is to represent the complex energy landscape of biomolecules using a small number of reaction coordinates. In this study, we investigate the folding pathways of an RNA tetraloop, gcGCAAgc, using five classical MD simulations with a combined simulation time of approximately 120 μs. Our approach involves analyzing the tetraloop dynamics, including the folding transition state ensembles, using the energy landscape visualization method (ELViM). The ELViM is an approach that uses internal distances to compare any two conformations, allowing for a detailed description of the folding process without requiring root mean square alignment of structures. This method has previously been applied to describe the energy landscape of disordered β-amyloid peptides and other proteins. The ELViM results in a non-linear projection of the multidimensional space, providing a comprehensive representation of the tetraloop's energy landscape. Our results reveal four distinct transition-state regions and establish the paths that lead to the folded tetraloop structure. This detailed analysis of the tetraloop's folding process has important implications for understanding RNA folding, and the ELViM approach can be used to study other biomolecules.

Probing the Energy Landscape of Spectrin R15 and R16 and the Effects of Non-native Interactions

Understanding the details of a protein folding mechanism can be a challenging and complex task. One system with an interesting folding behavior is the α-spectrin domain, where the R15 folds three-orders of magnitude faster than its homologues R16 and R17, despite having similar structures. The molecular origins that explain these folding rate differences remain unclear, but our previous work revealed that a combined effect produced by non-native interactions could be a reasonable cause for these differences. In this study, we explore further the folding process by identifying the molecular paths, metastable states, and the collective motions that lead these unfolded proteins to their native state conformation. Our results uncovered the differences between the folding pathways for the wild-type R15 and R16 and an R16 mutant. The metastable ensembles that speed down the folding were identified using an energy landscape visualization method (ELViM). These ensembles correspond to similar experimentally reported configurations. Our observations indicate that the non-native interactions are also associated with secondary structure misdocking. This computational methodology can be used as a fast, straightforward protocol for shedding light on systems with unclear folding or conformational traps.

Exploring the Folding Mechanism of Dimeric Superoxide Dismutase

The Cu/Zn Human Superoxide Dismutase (SOD1) is a dimeric metalloenzyme whose genetic mutations are directly related to amyotrophic lateral sclerosis (ALS), so understanding its folding mechanism is of fundamental importance. Currently, the SOD1 dimer formation is studied via molecular dynamics simulations using a simplified structure-based model and an all-atom model. Results from the simplified model reveal a mechanism dependent on distances between monomers, which are limited by constraints to mimic concentration dependence. The stability of intermediates (during the int state) is significantly affected by this distance, as well as by the presence of two folded monomers prior to dimer formation. The kinetics of interface formation are also highly dependent on the separation distance. The folding temperature of the dimer is about 4.2% higher than that of the monomer, a value not too different from experimental data. All-atom simulations on the apo dimer give binding free energy between monomers similar to experimental values. An intermediate state is evident for the apo form at a separation distance between monomers slightly larger than the native distance which has little formed interface between monomers. We have shown that this intermediate is stabilized by non-native intra- and intercontacts.

Multi-state Energy Landscape for Photoreaction of Stilbene and Dimethyl-stilbene

We have recently developed the reaction space projector (ReSPer) method, which constructs a reduced-dimensionality reaction space uniquely determined from reference reaction paths for a polyatomic molecular system and projects classical trajectories into the same reaction space. In this paper, we extend ReSPer to the analysis of photoreaction dynamics and relaxation processes of stilbene and present the concept of a "multi-state energy landscape," incorporating the ground- and excited-state reaction subspaces. The multi-state energy landscape successfully explains the previously established photoreaction processes of cis-stilbene, such as the cis-trans photoisomerization and photocyclization. In addition, we discuss the difference in the excited-state reaction dynamics between stilbene and 1,1'-dimethyl stilbene based on a common reaction subspace determined from the framework part of reference structures with different number of atoms. This approach allows us to target any molecule with a common framework, greatly expanding the applicability of the ReSPer analysis. The multi-state energy landscape provides fruitful insight into photochemical reactions, exploring the excited- and ground-state potential energy surfaces, as well as comprehensive reaction processes with nonradiative transitions between adiabatic states, within the stage of a reduced-dimensionality reaction space.

Reaction Space Projector (ReSPer) for Visualizing Dynamic Reaction Routes Based on Reduced-Dimension Space

To analyze chemical reaction dynamics based on a reaction path network, we have developed the "Reaction Space Projector" (ReSPer) method with the aid of the dimensionality reduction method. This program has two functions: the construction of a reduced-dimensionality reaction space from a molecular structure dataset, and the projection of dynamic trajectories into the low-dimensional reaction space. In this paper, we apply ReSPer to isomerization and bifurcation reactions of the Au₅ cluster and succeed in analyzing dynamic reaction routes involved in multiple elementary reaction processes, constructing complicated networks (called "closed islands") of nuclear permutation-inversion (NPI) isomerization reactions, and elucidating dynamic behaviors in bifurcation reactions with reference to bundles of trajectories. Interestingly, in the second application, we find a correspondence between the contribution ratios in the ability to visualize and the symmetry of the morphology of closed islands. In addition, the third application suggests the existence of boundaries that determine the selectivity in bifurcation reactions, which was discussed in the phase space. The ReSPer program is a versatile and robust tool to clarify dynamic reaction mechanisms based on the reduced-dimensionality reaction space without prior knowledge of target reactions.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

Resolving the fine structure in the energy landscapes of repeat proteins

Ankyrin (ANK) repeat proteins are coded by tandem occurrences of patterns with around 33 amino acids. They often mediate protein–protein interactions in a diversity of biological systems. These proteins have an elongated non-globular shape and often display complex folding mechanisms. This work investigates the energy landscape of representative proteins of this class made up of 3, 4 and 6 ANK repeats using the energy-landscape visualisation method (ELViM). By combining biased and unbiased coarse-grained molecular dynamics AWSEM simulations that sample conformations along the folding trajectories with the ELViM structure-based phase space, one finds a three-dimensional representation of the globally funnelled energy surface. In this representation, it is possible to delineate distinct folding pathways. We show that ELViMs can project, in a natural way, the intricacies of the highly dimensional energy landscapes encoded by the highly symmetric ankyrin repeat proteins into useful low-dimensional representations. These projections can discriminate between multiplicities of specific parallel folding mechanisms that otherwise can be hidden in oversimplified depictions.

Examining the Ensembles of Amyloid-β Monomer Variants and Their Propensities to Form Fibers Using an Energy Landscape Visualization Method

The amyloid-β (Aβ) monomer, an intrinsically disordered peptide, is produced by the cleavage of the amyloid precursor protein, leading to Aβ-40 and Aβ-42 as major products. These two isoforms generate pathological aggregates, whose accumulation correlates with Alzheimer's disease (AD). Experiments have shown that even though the natural abundance of Aβ-42 is smaller than that for Aβ-40, the Aβ-42 is more aggregation-prone compared to Aβ-40. Moreover, several single-point mutations are associated with early onset forms of AD. This work analyzes coarse-grained associative-memory, water-mediated, structure and energy model (AWSEM) simulations of normal Aβ-40 and Aβ-42 monomers, along with six single-point mutations associated with early onset disease. We analyzed the simulations using the energy landscape visualization method (ELViM), a reaction-coordinate-free approach suited to explore the frustrated energy landscapes of intrinsically disordered proteins. ELViM is shown to distinguish the monomer ensembles of variants that rapidly form fibers from those that do not form fibers as readily. It also delineates the amino acid contacts characterizing each ensemble. The results shed light on the potential of ELViM to probe intrinsically disordered proteins.

Exploring Energy Landscapes of Intrinsically Disordered Proteins: Insights into Functional Mechanisms

Intrinsically disordered proteins (IDPs) lack a rigid three-dimensional structure and populate a polymorphic ensemble of conformations. Because of the lack of a reference conformation, their energy landscape representation in terms of reaction coordinates presents a daunting challenge. Here, our newly developed energy landscape visualization method (ELViM), a reaction coordinate-free approach, shows its prime application to explore frustrated energy landscapes of an intrinsically disordered protein, prostate-associated gene 4 (PAGE4). PAGE4 is a transcriptional coactivator that potentiates the oncogene c-Jun. Two kinases, namely, HIPK1 and CLK2, phosphorylate PAGE4, generating variants phosphorylated at different serine/threonine residues (HIPK1-PAGE4 and CLK2-PAGE4, respectively) with opposing functions. While HIPK1-PAGE4 predominantly phosphorylates Thr51 and potentiates c-Jun, CLK2-PAGE4 hyperphosphorylates PAGE4 and attenuates transactivation. To understand the underlying mechanisms of conformational diversity among different phosphoforms, we have analyzed their atomistic trajectories simulated using AWSEM forcefield, and the energy landscapes were elucidated using ELViM. This method allows us to identify and compare the population distributions of different conformational ensembles of PAGE4 phosphoforms using the same effective phase space. The results reveal a predominant conformational ensemble with an extended C-terminal segment of WT PAGE4, which exposes a functional residue Thr51, implying its potential of undertaking a fly-casting mechanism while binding to its cognate partner. In contrast, for HIPK1-PAGE4, a compact conformational ensemble enhances its population sequestering phosphorylated-Thr51. This clearly explains the experimentally observed weaker affinity of HIPK1-PAGE4 for c-Jun. ELViM appears as a powerful tool, especially to analyze the highly frustrated energy landscape representation of IDPs where appropriate reaction coordinates are hard to apprehend.

pH and Charged Mutations Modulate Cold Shock Protein Folding and Stability: A Constant pH Monte Carlo Study

The folding and stability of proteins is a fundamental problem in several research fields. In the present paper, we have used different computational approaches to study the effects caused by changes in pH and for charged mutations in cold shock proteins from Bacillus subtilis (Bs-CspB). First, we have investigated the contribution of each ionizable residue for these proteins to their thermal stability using the TKSA-MC, a Web server for rational mutation via optimizing the protein charge interactions. Based on these results, we have proposed a new mutation in an already optimized Bs-CspB variant. We have evaluated the effects of this new mutation in the folding energy landscape using structure-based models in Monte Carlo simulation at constant pH, SBM-CpHMC. Our results using this approach have indicated that the charge rearrangements already in the unfolded state are critical to the thermal stability of Bs-CspB. Furthermore, the conjunction of these simplified methods was able not only to predict stabilizing mutations in different pHs but also to provide essential information about their effects in each stage of protein folding.