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.

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.

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.

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.

pH and the Breast Cancer Recurrent Mutation D538G Affect the Process of Activation of Estrogen Receptor α.

Estrogen receptor α (ERα) is a regulatory protein that can access a set of distinct structural configurations. ERα undergoes extensive remodeling as it interacts with different agonists and antagonists, as well as transcription activation and repression factors. Moreover, breast cancer tumors resistant to hormone therapy have been associated with the imbalance between the active and inactive ERα states. Cancer-activating mutations in ERα play a crucial role in this imbalance and can promote the progression of cancer. However, the rate of this progression can also be increased by dysregulated pH in the tumor microenvironment. Many molecular aspects of the process of activation of ERα that can be affected by these pH changes and mutations are still unclear. Thus, we applied computational and experimental techniques to explore the activation process dynamics of ER for environments with different pHs and in the presence of one of the most recurrent cancer-activating mutations, D538G. Our results indicated that the effect of the pH increase associated with the D538G mutation promoted a robust stabilization of the active state of ER. We were also able to determine the main protein regions that have the most potential to influence the activation process under different pH conditions, which may provide targets of future therapeutics for the treatment of hormone-resistant breast cancer tumors. Finally, the approach used here can be applied for proteins associated with the proliferation of other cancer types, which can also have their function affected by small pH changes.

In silico identification of noncompetitive inhibitors targeting an uncharacterized allosteric site of falcipain-2.

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme.

Water Bridges Play a Key Role in Affinity and Selectivity for Malarial Protease Falcipain-2.

Falcipain-2 (FP-2) is hemoglobinase considered an attractive drug target of Plasmodium falciparum. Recently, it has been shown that peptidomimetic nitriles containing a 3-pyridyl (3Pyr) moiety at P2 display high affinity and selectivity for FP-2 with respect to human cysteine cathepsins (hCats), outperforming other P2-Pyr isomers and analogs. Further characterization demonstrated that certain P3 variants of these compounds possess micromolar inhibition of parasite growth in vitro and no cytotoxicity against human cell lines. However, the structural determinants underlying the selectivity of the 3Pyr-containing nitriles for FP-2 remain unknown. In this work, we conduct a thorough computational study combining MD simulations and free energy calculations to decipher the bases of the selectivity of the aforementioned nitriles. Our results reveal that water bridges involving the nitrogen and one carboxyl oxygen of I85 and D234 of FP-2, respectively, and the nitrogen of the neutral 3Pyr moiety, which are either less prevalent or nonexistent in the other complexes, explain the experimental activity profiles. The presence of crystallographic waters close to the bridging water positions in the studied proteases strongly supports the occurrence of such interactions. Overall, our findings suggest that selective FP-2 inhibitors can be designed by promoting water bridge formation at the bottom of the S2 subsite and/or by introducing complementary groups that displace the bridging water.

Rational Design of Chymotrypsin Inhibitor 2 by Optimizing Non-Native Interactions.

Rational design of proteins via mutagenesis is crucial for several biotechnological applications. A significant challenge of the computational strategies used to predict optimized mutations is to understand the influence of each amino acid during the folding process. In the present work, chymotrypsin inhibitor 2 (CI2) and several of its designed mutants have been simulated using a non-native hydrophobic and electrostatic potential as a structure-based Cα model. Through these simulations, we could identify the most critical folding stage to accelerate CI2 and also the charged residues responsible for providing its thermostability. The replacement of ionizable residues for hydrophobic ones tended to promote the formation of the CI2 secondary structure in the early transition state, which speeds up folding. However, this same replacement destabilized the native structure, and there was a decrease in the protein thermostability. Such a simple method proved to be capable of providing valuable information about thermodynamics and kinetics of CI2 and its mutations, thus being a fast alternative to the study of rational protein design.

Diffusion of tRNA inside the ribosome is position-dependent.

In recent years, there has been a growing interest to quantify the energy landscape that governs ribosome dynamics. However, in order to quantitatively integrate theoretical predictions and experimental measurements, it is essential that one has a detailed understanding of the associated diffusive properties. Here, all-atom explicit-solvent simulations (50 µs of aggregate sampling) predict that the diffusion coefficient of a tRNA molecule will depend on its position within the ribosome. Specifically, during aa-tRNA accommodation (i.e., the process by which tRNA enters the ribosome), the apparent diffusion coefficient decreases by approximately an order of magnitude. By comparing these to values obtained with an energetically "smooth" model, we show that the observed nonuniform behavior likely arises from electrostatic and solvation interactions between the tRNA and ribosome. These calculations also reveal the hierarchical character of ribosomal energetics, where steric interactions induce a large-scale free-energy barrier, and short-scale roughness determines the rate of diffusive movement across the landscape.

TKSA-MC: A web server for rational mutation through the optimization of protein charge interactions.

The TKSAMC is a web server which calculates protein charge-charge interactions via the Tanford-Kirkwood Surface Accessibility model with the Monte Carlo method for sampling different protein protonation states. The optimization of charge-charge interactions via directed mutations has successfully enhanced the thermal stability of different proteins and could be a key to protein engineering improvement. The server presents the electrostatic free energy contribution of each polar-charged residue to the protein native state stability. Specific residues are suggested to be mutated for improving thermal stability. The choice of a residue is based on its fraction of side chain exposed to solvent and its positive free energy contribution, which tends to destabilize the protein native state. Any residue energy contribution can be shown as a function of pH condition. The web server is freely available at UNESP (São Paulo State University - DF/IBILCE): http://tksamc.df.ibilce.unesp.br and also on GitHub https://github.com/contessoto/tksamc.

Topography of funneled landscapes determines the thermodynamics and kinetics of protein folding.

The energy landscape approach has played a fundamental role in advancing our understanding of protein folding. Here, we quantify protein folding energy landscapes by exploring the underlying density of states. We identify three quantities essential for characterizing landscape topography: the stabilizing energy gap between the native and nonnative ensembles δE, the energetic roughness ΔE, and the scale of landscape measured by the entropy S. We show that the dimensionless ratio between the gap, roughness, and entropy of the system Λ=δE/(ΔE√(2S)) accurately predicts the thermodynamics, as well as the kinetics of folding. Large Λ implies that the energy gap (or landscape slope towards the native state) is dominant, leading to more funneled landscapes. We investigate the role of topological and energetic roughness for proteins of different sizes and for proteins of the same size, but with different structural topologies. The landscape topography ratio Λ is shown to be monotonically correlated with the thermodynamic stability against trapping, as characterized by the ratio of folding temperature versus trapping temperature. Furthermore, Λ also monotonically correlates with the folding kinetic rates. These results provide the quantitative bridge between the landscape topography and experimental folding measurements.

Capturing transition paths and transition states for conformational rearrangements in the ribosome.

To reveal the molecular determinants of biological function, one seeks to characterize the interactions that are formed in conformational and chemical transition states. In other words, what interactions govern the molecule's energy landscape? To accomplish this, it is necessary to determine which degrees of freedom can unambiguously identify each transition state. Here, we perform simulations of large-scale aminoacyl-transfer RNA (aa-tRNA) rearrangements during accommodation on the ribosome and project the dynamics along experimentally accessible atomic distances. From this analysis, we obtain evidence for which coordinates capture the correct number of barrier-crossing events and accurately indicate when the aa-tRNA is on a transition path. Although a commonly used coordinate in single-molecule experiments performs poorly, this study implicates alternative coordinates along which rearrangements are accurately described as diffusive movements across a one-dimensional free-energy profile. From this, we provide the theoretical foundation required for single-molecule techniques to uncover the energy landscape governing aa-tRNA selection by the ribosome.

Assembling a xylanase-lichenase chimera through all-atom molecular dynamics simulations.

Multifunctional enzyme engineering can improve enzyme cocktails for emerging biofuel technology. Molecular dynamics through structure-based models (SB) is an effective tool for assessing the tridimensional arrangement of chimeric enzymes as well as for inferring the functional practicability before experimental validation. This study describes the computational design of a bifunctional xylanase-lichenase chimera (XylLich) using the xynA and bglS genes from Bacillus subtilis. In silico analysis of the average solvent accessible surface area (SAS) and the root mean square fluctuation (RMSF) predicted a fully functional chimera, with minor fluctuations and variations along the polypeptide chains. Afterwards, the chimeric enzyme was built by fusing the xynA and bglS genes. XylLich was evaluated through small-angle X-ray scattering (SAXS) experiments, resulting in scattering curves with a very accurate fit to the theoretical protein model. The chimera preserved the biochemical characteristics of the parental enzymes, with the exception of a slight variation in the temperature of operation and the catalytic efficiency (kcat/Km). The absence of substantial shifts in the catalytic mode of operation was also verified. Furthermore, the production of chimeric enzymes could be more profitable than producing a single enzyme separately, based on comparing the recombinant protein production yield and the hydrolytic activity achieved for XylLich with that of the parental enzymes.

Inferring a weighted elastic network from partial unfolding with coarse-grained simulations.

A number of studies have demonstrated that simple elastic network models can reproduce experimental B-factors, providing insights into the structure-function properties of proteins. Here, we report a study on how to improve an elastic network model and explore its performance by predicting the experimental B-factors. Elastic network models are built on the experimental Cα coordinates, and they only take the pairs of Cα atoms within a given cutoff distance rc into account. These models describe the interactions by elastic springs with the same force constant. We have developed a method based on numerical simulations with a simple coarse-grained force field, to attribute weights to these spring constants. This method considers the time that two Cα atoms remain connected in the network during partial unfolding, establishing a means of measuring the strength of each link. We examined two different coarse-grained force fields and explored the computation of these weights by unfolding the native structures.

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.

Analyzing the effect of homogeneous frustration in protein folding.

The energy landscape theory has been an invaluable theoretical framework in the understanding of biological processes such as protein folding, oligomerization, and functional transitions. According to the theory, the energy landscape of protein folding is funneled toward the native state, a conformational state that is consistent with the principle of minimal frustration. It has been accepted that real proteins are selected through natural evolution, satisfying the minimum frustration criterion. However, there is evidence that a low degree of frustration accelerates folding. We examined the interplay between topological and energetic protein frustration. We employed a Cα structure-based model for simulations with a controlled nonspecific energetic frustration added to the potential energy function. Thermodynamics and kinetics of a group of 19 proteins are completely characterized as a function of increasing level of energetic frustration. We observed two well-separated groups of proteins: one group where a little frustration enhances folding rates to an optimal value and another where any energetic frustration slows down folding. Protein energetic frustration regimes and their mechanisms are explained by the role of non-native contact interactions in different folding scenarios. These findings strongly correlate with the protein free-energy folding barrier and the absolute contact order parameters. These computational results are corroborated by principal component analysis and partial least square techniques. One simple theoretical model is proposed as a useful tool for experimentalists to predict the limits of improvements in real proteins.

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.

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.

Coarse-Grained Simulations of Protein Folding: Bridging Theory and Experiments.

Computational coarse-grained models play a fundamental role as a research tool in protein folding, and they are important in bridging theory and experiments. Folding mechanisms are generally discussed using the energy landscape framework, which is well mapped within a class of simplified structure-based models. In this chapter, simplified computer models are discussed with special focus on structure-based ones.

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.