bims-lances Biomed News
on Landscapes from Cryo-EM and Simulations
Issue of 2024–03–24
seven papers selected by
James M. Krieger, National Centre for Biotechnology



  1. Biochemistry. 2024 Mar 19.
      Ras-related nuclear protein (Ran) is a member of the Ras superfamily of small guanosine triphosphatases (GTPases) and a regulator of multiple cellular processes. In healthy cells, the GTP-bound form of Ran is concentrated at chromatin, creating a Ran•GTP gradient that provides the driving force for nucleocytoplasmic transport, mitotic spindle assembly, and nuclear envelope formation. The Ran•GTP gradient is maintained by the regulator of chromatin condensation 1 (RCC1), a guanine nucleotide exchange factor that accelerates GDP/GTP exchange in Ran. RCC1 interacts with nucleosomes, which are the fundamental repeating units of eukaryotic chromatin. Here, we present a cryo-EM analysis of a trimeric complex composed of the nucleosome core particle (NCP), RCC1, and Ran. While the contacts between RCC1 and Ran in the complex are preserved compared with a previously determined structure of RCC1-Ran, our study reveals that RCC1 and Ran interact dynamically with the NCP and undergo rocking motions on the nucleosome surface. Furthermore, the switch 1 region of Ran, which plays an important role in mediating conformational changes associated with the substitution of GDP and GTP nucleotides in Ras family members, appears to undergo disorder-order transitions and forms transient contacts with the C-terminal helix of histone H2B. Nucleotide exchange assays performed in the presence and absence of NCPs are not consistent with an active role for nucleosomes in nucleotide exchange, at least in vitro. Instead, the nucleosome stabilizes RCC1 and serves as a hub that concentrates RCC1 and Ran to promote efficient Ran•GDP to Ran•GTP conversion.
    DOI:  https://doi.org/10.1021/acs.biochem.3c00724
  2. bioRxiv. 2024 Mar 06. pii: 2024.03.02.583144. [Epub ahead of print]
      α-Klotho (KLA) is a type-1 membranous protein that can associate with fibroblast growth factor receptor (FGFR) to form co-receptor for FGF23. The ectodomain of unassociated KLA is shed as soluble KLA (sKLA) to exert FGFR/FGF23-independent pleiotropic functions. The previously determined X-ray crystal structure of the extracellular region of sKLA in complex with FGF23 and FGFR1c suggests that sKLA functions solely as an on-demand coreceptor for FGF23. To understand the FGFR/FGF23-independent pleiotropic functions of sKLA, we investigated biophysical properties and structure of apo-sKLA. Mass photometry revealed that sKLA can form a stable structure with FGFR and/or FGF23 as well as sKLA dimer in solution. Single particle cryogenic electron microscopy (cryo-EM) supported the dimeric structure of sKLA. Cryo-EM further revealed a 3.3Å resolution structure of apo-sKLA that overlays well with its counterpart in the ternary complex with several distinct features. Compared to the ternary complex, the KL2 domain of apo-sKLA is more flexible. 3D variability analysis revealed that apo-sKLA adopts conformations with different KL1-KL2 interdomain bending and rotational angles. The potential multiple forms and shapes of sKLA support its role as FGFR-independent hormone with pleiotropic functions. A comprehensive understanding of the sKLA conformational landscape will provide the foundation for developing klotho-related therapies for diseases.
    DOI:  https://doi.org/10.1101/2024.03.02.583144
  3. bioRxiv. 2024 Feb 20. pii: 2024.02.15.580591. [Epub ahead of print]
      The groundbreaking achievements of AlphaFold2 (AF2) approaches in protein structure modeling marked a transformative era in structural biology. Despite the success of AF2 tools in predicting single protein structures, these methods showed intrinsic limitations in predicting multiple functional conformations of allosteric proteins and fold-switching systems. The recent NMR-based structural determination of the unbound ABL kinase in the active state and two inactive low-populated functional conformations that are unique for ABL kinase presents an ideal challenge for AF2 approaches. In the current study we employ several implementations of AF2 methods to predict protein conformational ensembles and allosteric states of the ABL kinase including (a) multiple sequence alignments (MSA) subsampling approach; (b) SPEACH_AF approach in which alanine scanning is performed on generated MSAs; and (c) introduced in this study randomized full sequence mutational scanning for manipulation of sequence variations combined with the MSA subsampling. We show that the proposed AF2 adaptation combined with local frustration mapping of conformational states enable accurate prediction of the ABL active and intermediate structures and conformational ensembles, also offering a robust approach for interpretable characterization of the AF2 predictions and limitations in detecting hidden allosteric states. We found that the large high frustration residue clusters are uniquely characteristic of the low-populated, fully inactive ABL form and can define energetically frustrated cracking sites of conformational transitions, presenting difficult targets for AF2 methods. This study uncovered previously unappreciated, fundamental connections between distinct patterns of local frustration in functional kinase states and AF2 successes/limitations in detecting low-populated frustrated conformations, providing a better understanding of benefits and limitations of current AF2-based adaptations in modeling of conformational ensembles.
    DOI:  https://doi.org/10.1101/2024.02.15.580591
  4. Methods Mol Biol. 2024 ;2754 77-90
      The electron microscopy metainference integrative structural biology method enables the combination of cryo-electron microscopy electron density maps with molecular modeling techniques such as molecular dynamics to unveil the atomistic biomolecular structural ensemble and the error in the map data in an efficient manner. Here we illustrate the electron microscopy metainference protocol and analysis used to elucidate the atomistic structural ensemble of the microtubule-associated protein tau bound to microtubules by using state-of-the-art molecular mechanic force field and the electron density map.
    Keywords:  Alzheimer’s disease; Bayesian inference for cryo-EM data modeling; Integrative structural biology; Microtubules; Protein functional dynamics; Structural ensemble determination; Tau
    DOI:  https://doi.org/10.1007/978-1-0716-3629-9_4
  5. Biophys Physicobiol. 2023 ;20(2): e200022
      Protein functions associated with biological activity are precisely regulated by both tertiary structure and dynamic behavior. Thus, elucidating the high-resolution structures and quantitative information on in-solution dynamics is essential for understanding the molecular mechanisms. The main experimental approaches for determining tertiary structures include nuclear magnetic resonance (NMR), X-ray crystallography, and cryogenic electron microscopy (cryo-EM). Among these procedures, recent remarkable advances in the hardware and analytical techniques of cryo-EM have increasingly determined novel atomic structures of macromolecules, especially those with large molecular weights and complex assemblies. In addition to these experimental approaches, deep learning techniques, such as AlphaFold 2, accurately predict structures from amino acid sequences, accelerating structural biology research. Meanwhile, the quantitative analyses of the protein dynamics are conducted using experimental approaches, such as NMR and hydrogen-deuterium mass spectrometry, and computational approaches, such as molecular dynamics (MD) simulations. Although these procedures can quantitatively explore dynamic behavior at high resolution, the fundamental difficulties, such as signal crowding and high computational cost, greatly hinder their application to large and complex biological macromolecules. In recent years, machine learning techniques, especially deep learning techniques, have been actively applied to structural data to identify features that are difficult for humans to recognize from big data. Here, we review our approach to accurately estimate dynamic properties associated with local fluctuations from three-dimensional cryo-EM density data using a deep learning technique combined with MD simulations.
    Keywords:  3D-CNN; allostery; big data; machine learning; single particle analysis
    DOI:  https://doi.org/10.2142/biophysico.bppb-v20.0022
  6. J Chem Inf Model. 2024 Mar 20.
      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.
    DOI:  https://doi.org/10.1021/acs.jcim.4c00034
  7. Biophys J. 2024 Mar 21. pii: S0006-3495(24)00204-2. [Epub ahead of print]
      Comparative methods in molecular evolution and structural biology rely heavily upon the site-wise analysis of DNA sequence and protein structure, both static forms of information. However, it is widely accepted that protein function results from nanoscale non-random machine-like motions induced by evolutionarily conserved molecular interactions. Comparisons of molecular dynamics (MD) simulations conducted between homologous sites representative of different functional or mutational states can potentially identify local effects on binding interaction and protein evolution. Additionally, comparisons of different (i.e. non-homologous) sites within MD simulations could be employed to identify functional shifts in local time-coordinated dynamics indicative of logic-gating within proteins. However, comparative MD analysis is challenged by the large fraction of protein motion caused by random thermal noise in the surrounding solvent. Therefore, properly de-noised MD comparisons could reveal functional sites involving these machine-like dynamics with good accuracy. Here, we introduce ATOMDANCE, a user-interfaced suite of comparative machine learning based de-noising tools designed for identifying functional sites and the patterns of coordinated motion they can create within MD simulations. ATOMDANCE-maxDemon4.0 employs Gaussian kernel functions to compute site-wise maximum mean discrepancy (MMD) between learned features of motion, thereby assessing de-noised differences in the non-random motions between functional or evolutionary states (e.g. ligand bound vs. unbound, wild-type vs. mutant). ATOMDANCE-maxDemon4.0 also employs MMD to analyze potential random amino-acid replacements allowing for a site-wise test of neutral vs. non-neutral evolution on the divergence of dynamic function in protein homologs. Lastly, ATOMDANCE-Choreograph2.0 employs mixed-model ANOVA and graph network to detect regions where time synchronized shifts in dynamics occur. Here, we demonstrate ATOMDANCE's utility for identifying key sites involved in dynamic responses during functional binding interactions involving DNA, small molecule drugs, and virus-host recognition, as well as understanding shifts in global and local site coordination occurring during allosteric activation of a pathogenic protease.
    Keywords:  allostery; comparative method; coordinated dynamics; machine learning; molecular dynamics simulation; molecular evolution
    DOI:  https://doi.org/10.1016/j.bpj.2024.03.024