Modeling and Simulation of the NMDA Receptor at Coarse-Grained and Atomistic Levels.

N-Methyl-D-aspartate (NMDA) receptors are glutamate-gated excitatory channels that play essential roles in brain functions. While high-resolution structures were solved for an allosterically inhibited form of functional NMDA receptor, other key functional states (particularly the active open-channel state) have not yet been resolved at atomic resolutions. To decrypt the molecular mechanism of the NMDA receptor activation, structural modeling and simulation are instrumental in providing detailed information about the dynamics and energetics of the receptor in various functional states. In this chapter, we describe coarse-grained modeling of the NMDA receptor using an elastic network model and related modeling/analysis tools (e.g., normal mode analysis, flexibility and hotspot analysis, cryo-EM flexible fitting, and transition pathway modeling) based on available structures. Additionally, we show how to build an atomistic model of the active-state receptor with targeted molecular dynamics (MD) simulation and explore its energetics and dynamics with conventional MD simulation. Taken together, these modeling and simulation can offer rich structural and dynamic information which will guide experimental studies of the activation of this key receptor.

Nanomechanical collective vibration of SARS-CoV-2 spike proteins.

The development of effective therapeutics against COVID-19 requires a thorough understanding of the receptor recognition mechanism of the SARS-CoV-2 spike (S) protein. Here the multidomain collective dynamics on the trimer of the spike protein has been analyzed using normal mode analysis (NMA). A common nanomechanical profile was identified in the spike proteins of SARS-CoV-2 and its variants. The profile involves collective vibrations of the receptor-binding domain (RBD) and the N-terminal domain (NTD), which may mediate the physical interaction process. Quantitative analysis of the collective modes suggests a nanomechanical property involving large-scale conformational changes, which explains the difference in receptor binding affinity among different variants. These results support the use of intrinsic global dynamics as a valuable perspective for studying the allosteric and functional mechanisms of the S protein. This approach also provides a low-cost theoretical toolkit for screening potential pathogenic mutations and drug targets.

Design and computational evaluation of a novel multi-epitope hybrid vaccine against monkeypox virus: Potential targets and immunogenicity assessment for pandemic preparedness.

Monkeypox is a type of DNA-enveloped virus that belongs to the orthopoxvirus family, closely related to the smallpox virus. It can cause an infectious disease in humans known as monkeypox disease. Although there are multiple drugs and vaccines designed to combat orthopoxvirus infections, with a primary focus on smallpox, the recent spread of the monkeypox virus to over 50 countries have ignited a mounting global concern. This unchecked viral proliferation has raised apprehensions about the potential for a pandemic corresponding to the catastrophic impact of COVID-19. This investigation explored the structural proteins of monkeypox virus as potential candidates for designing a novel hybrid multi-epitope vaccine. The epitopes obtained from the selected proteins were screened to ensure their non-allergenicity, non-toxicity, and antigenicity to trigger T and B-cell responses. The interaction of the vaccine with toll-like receptor-3 (TLR-3) and major histocompatibility complexes (MHCs) was assessed using Cluspro 2.0. To establish the reliability of the docked complexes, a comprehensive evaluation was conducted using Immune and MD Simulations and Normal Mode Analysis. However, to validate the computational results of this study, additional in-vitro and in-vivo research is essential.

How alginate lyase produces quasi-monodisperse oligosaccharides: A normal-mode-based docking and molecular dynamics simulation study.

Oligosaccharide degradation products of alginate (AOS) hold significant potential in diverse fields, including pharmaceuticals, health foods, textiles, and agricultural production. Enzymatic alginate degradation is appealing due to its mild conditions, predictable activity, high yields, and controllability. However, the alginate degradation often results in a complex mixture of oligosaccharides, necessitating costly purification to isolate highly active oligosaccharides with a specific degree of polymerization (DP). Addressing this, our study centers on the alginate lyase AlyB from Vibrio Splendidus OU02, which uniquely breaks down alginate into mono-distributed trisaccharides. This enzyme features a polysaccharide lyase family 7 domain (PL-7) and a CBM32 carbohydrate-binding module connected by a helical structure. Through normal-mode-based docking and all-atom molecular simulations, we demonstrate that AlyB's substrate and product specificities are influenced by the spatial conformation of the catalytic pocket and the flexibility of its structure. The helically attached CBM is pivotal in releasing trisaccharides, which is crucial for avoiding further degradation. This study sheds light on AlyB's specificity and efficiency and contributes to the evolving field of enzyme design for producing targeted oligosaccharides, with significant implications for various bioindustries.

Elastic network model reveals distinct flexibilities of capping proteins bound to CARMIL and twinfilin-tail.

Capping protein (CP) binds to the barbed end of an actin-filament and inhibits its elongation. CARMIL binds CP and dissociates it from the barbed end of the actin-filament. The binding of CARMIL peptide alters the flexibility of CP, which is considered to facilitate the dissociation. Twinfilin also binds to CP through its C-terminal tail. The complex structures of the CP/twinfilin-tail (TW-tail) peptide indicate that the binding sites of CARMIL and TW-tail overlap. However, TW-tail binding does not facilitate the dissociation of CP from the barbed end. We extensively investigated the flexibilities of CP in the CP/TW-tail or CP/CARMIL complexes using an elastic network model and concluded that TW-tail binding does not alter the flexibility of CP. Our extensive analysis also highlighted that the strong contacts of peptides with the two domains of CP, that is, the CP-L and CP-S domains, are key to changing the flexibilities of CP. CARMIL peptides can interact strongly with both of the domains, while TW-tail peptides exclusively interact with the CP-S domain because the binding site of TW-tail on CP relatively shifts to the CP-S domain compared with that of CP/CARMIL. This result supports our hypothesis that the dissociation of CP from the barbed end is regulated by the flexibility of CP.

Screening of phytochemicals from Clerodendrum inerme (L.) Gaertn as potential anti-breast cancer compounds targeting EGFR: an in-silico approach.

Breast cancer (BC) is the most prevalent malignancy among women around the world. The epidermal growth factor receptor (EGFR) is a tyrosine kinase receptor (RTK) of the ErbB/HER family. It is essential for triggering the cellular signaling cascades that control cell growth and survival. However, perturbations in EGFR signaling lead to cancer development and progression. Hence, EGFR is regarded as a prominent therapeutic target for breast cancer. Therefore, in the current investigation, EGFR was targeted with phytochemicals from Clerodendrum inerme (L.) Gaertn (C. inerme). A total of 121 phytochemicals identified by gas chromatography-mass spectrometry (GC-MS) analysis were screened against EGFR through molecular docking, ADMET analyses (Absorption, Distribution, Metabolism, Excretion, and Toxicity), PASS predictions, and molecular dynamics simulation, which revealed three potential hit compounds with CIDs 10586 [i.e. alpha-bisabolol (-6.4 kcal/mol)], 550281 [i.e. 2,(4,4-Trimethyl-3-hydroxymethyl-5a-(3-methyl-but-2-enyl)-cyclohexene) (-6.5 kcal/mol)], and 161271 [i.e. salvigenin (-7.4 kcal/mol)]. The FDA-approved drug gefitinib was used to compare the inhibitory effects of the phytochemicals. The top selected compounds exhibited good ADMET properties and obeyed Lipinski's rule of five (ROF). The molecular docking analysis showed that salvigenin was the best among the three compounds and formed bonds with the key residue Met 793. Furthermore, the molecular mechanics generalized born surface area (MMGBSA) calculations, molecular dynamics simulation, and normal mode analysis validated the binding affinity of the compounds and also revealed the strong stability and compactness of phytochemicals at the docked site. Additionally, DFT and DOS analyses were done to study the reactivity of the compounds and to further validate the selected phytochemicals. These results suggest that the identified phytochemicals possess high inhibitory potential against the target EGFR and can treat breast cancer. However, further in vitro and in vivo investigations are warranted towards the development of these constituents into novel anti-cancer drugs.Communicated by Ramaswamy H. Sarma.

Structural and pKa Estimation of the Amphipathic HR1 in SARS-CoV-2: Insights from Constant pH MD, Linear vs. Nonlinear Normal Mode Analysis.

A comprehensive understanding of molecular interactions and functions is imperative for unraveling the intricacies of viral protein behavior and conformational dynamics during cellular entry. Focusing on the SARS-CoV-2 spike protein (SARS-CoV-2 sp), a Principal Component Analysis (PCA) on a subset comprising 131 A-chain structures in presence of various inhibitors was conducted. Our analyses unveiled a compelling correlation between PCA modes and Anisotropic Network Model (ANM) modes, underscoring the reliability and functional significance of low-frequency modes in adapting to diverse inhibitor binding scenarios. The role of HR1 in viral processing, both linear Normal Mode Analysis (NMA) and Nonlinear NMA were implemented. Linear NMA exhibited substantial inter-structure variability, as evident from a higher Root Mean Square Deviation (RMSD) range (7.30 Å), nonlinear NMA show stability throughout the simulations (RMSD 4.85 Å). Frequency analysis further emphasized that the energy requirements for conformational changes in nonlinear modes are notably lower compared to their linear counterparts. Using simulations of molecular dynamics at constant pH (cpH-MD), we successfully predicted the pKa order of the interconnected residues within the HR1 mutations at lower pH values, suggesting a transition to a post-fusion structure. The pKa determination study illustrates the profound effects of pH variations on protein structure. Key results include pKa values of 9.5179 for lys-921 in the D936H mutant, 9.50 for the D950N mutant, and a slightly higher value of 10.49 for the D936Y variant. To further understand the behavior and physicochemical characteristics of the protein in a biologically relevant setting, we also examine hydrophobic regions in the prefused states of the HR1 protein mutants D950N, D936Y, and D936H in our study. This analysis was conducted to ascertain the hydrophobic moment of the protein within a lipid environment, shedding light on its behavior and physicochemical properties in a biologically relevant context.

Comparative analysis of permanent and transient domain-domain interactions in multi-domain proteins.

Protein domains are structural, functional, and evolutionary units. These domains bring out the diversity of functionality by means of interactions with other co-existing domains and provide stability. Hence, it is important to study intra-protein inter-domain interactions from the perspective of types of interactions. Domains within a chain could interact over short timeframes or permanently, rather like protein-protein interactions (PPIs). However, no systematic study has been carried out between two classes, namely permanent and transient domain-domain interactions. In this work, we studied 263 two-domain proteins, belonging to either of these classes and their interfaces on the basis of several factors, such as interface area and details of interactions (number, strength, and types of interactions). We also characterized them based on residue conservation at the interface, correlation of residue motions across domains, its involvement in repeat formation, and their involvement in particular molecular processes. Finally, we could analyze the interactions arising from domains in two-domain monomeric proteins, and we observed significant differences between these two classes of domain interactions and a few similarities. This study will help to obtain a better understanding of structure-function and folding principles of multi-domain proteins.

Normal mode analysis and comparative study of intrinsic dynamics of alcohol oxidase enzymes from GMC protein family.

Glucose-Methanol-Choline (GMC) family enzymes are very important in catalyzing the oxidation of a wide range of structurally diverse substrates. Enzymes that constitute the GMC family, share a common tertiary fold but < 25% sequence identity. Cofactor FAD, FAD binding signature motif, and similar structural scaffold of the active site are common features of oxidoreductase enzymes of the GMC family. Protein functionality mainly depends on protein three-dimensional structures and dynamics. In this study, we used the normal mode analysis method to search the intrinsic dynamics of GMC family enzymes. We have explored the dynamical behavior of enzymes with unique substrate catabolism and active site characteristics from different classes of the GMC family. Analysis of individual enzymes and comparative ensemble analysis of enzymes from different classes has shown conserved dynamic motion at FAD binding sites. The present study revealed that GMC enzymes share a strong dynamic similarity (Bhattacharyya coefficient >90% and root mean squared inner product >52%) despite low sequence identity across the GMC family enzymes. The study predicts that local deformation energy between atoms of the enzyme may be responsible for the catalysis of different substrates. This study may help that intrinsic dynamics can be used to make meaningful classifications of proteins or enzymes from different organisms.Communicated by Ramaswamy H. Sarma.

Scipion-EM-ProDy: A Graphical Interface for the ProDy Python Package within the Scipion Workflow Engine Enabling Integration of Databases, Simulations and Cryo-Electron Microscopy Image Processing.

Macromolecular assemblies, such as protein complexes, undergo continuous structural dynamics, including global reconfigurations critical for their function. Two fast analytical methods are widely used to study these global dynamics, namely elastic network model normal mode analysis and principal component analysis of ensembles of structures. These approaches have found wide use in various computational studies, driving the development of complex pipelines in several software packages. One common theme has been conformational sampling through hybrid simulations incorporating all-atom molecular dynamics and global modes of motion. However, wide functionality is only available for experienced programmers with limited capabilities for other users. We have, therefore, integrated one popular and extensively developed software for such analyses, the ProDy Python application programming interface, into the Scipion workflow engine. This enables a wider range of users to access a complete range of macromolecular dynamics pipelines beyond the core functionalities available in its command-line applications and the normal mode wizard in VMD. The new protocols and pipelines can be further expanded and integrated into larger workflows, together with other software packages for cryo-electron microscopy image analysis and molecular simulations. We present the resulting plugin, Scipion-EM-ProDy, in detail, highlighting the rich functionality made available by its development.

Defensin-based therapeutic peptide design in attenuating V30M TTR-induced Familial Amyloid Polyneuropathy.

In the present study, we aimed to formulate an effective therapeutic candidate against V30M mutant transthyretin (TTR) protein to hinder its pathogenic misfolding. Nicotiana alata Defensin 1 (NaD1) Antimicrobial Peptide (AMP) was availed due to its tendency to aggregate, which may compete for aggregation-prone regions of pathogenic TTR protein. Based on NaD1's potential to bind to V30M TTR, we proposed NaD1-derived tetra peptides: CKTE and SKIL to be initial therapeutic candidates. Based on their association with mutant TTR protein, CKTE tetra peptide showed considerable interaction and curative potential as compared to SKIL tetra peptide. Further analyses from discrete molecular dynamics simulation corroborate CKTE tetra peptide's effectiveness as a 'beta-sheet breaker' against V30M TTR. Various post-simulation trajectory analyses suggested that CKTE tetra peptide alters the structural dynamics of pathogenic V30M TTR protein, thereby potentially attenuating its beta-sheets and impeding its aggregation. Normal mode analysis simulation corroborated that V30M TTR conformation is altered upon its interaction with CKTE peptide. Moreover, simulated thermal denaturation findings suggested that CKTE-V30M TTR complex is more susceptible to simulated denaturation, relative to pathogenic V30M TTR; further substantiating CKTE peptide's potential to alter V30M TTR's pathogenic conformation. Moreover, the residual frustration analysis augmented CKTE tetra peptide's proclivity in reorienting the conformation of V30M TTR. Therefore, we predicted that the tetra peptide, CKTE could be a promising therapeutic candidate in mitigating the amyloidogenic detrimental effects of V30M TTR-mediated familial amyloid polyneuropathy (FAP). Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-023-03646-4.

Network analysis of chromophore binding site in LOV domain.

Photoreceptor proteins are versatile toolbox for developing biosensors for optogenetic applications. These molecular tools get activated upon illumination of blue light, which in turn offers a non-invasive method for gaining high spatiotemporal resolution and precise control of cellular signal transduction. The Light-Oxygen-Voltage (LOV) domain family of proteins is a well-recognized system for constructing optogenetic devices. Translation of these proteins into efficient cellular sensors is possible by tuning their photochemistry lifetime. However, the bottleneck is the need for more understanding of the relationship between the protein environment and photocycle kinetics. Significantly, the effect of the local environment also modulates the electronic structure of chromophore, which perturbs the electrostatic and hydrophobic interaction within the binding site. This work highlights the critical factors hidden in the protein networks, linking with their experimental photocycle kinetics. It presents an opportunity to quantitatively examine the alternation in chromophore's equilibrium geometry and identify details which have substantial implications in designing synthetic LOV constructs with desirable photocycle efficiency.

Insight into the conserved structural dynamics of the C-terminus of mammal PrPC identifies structural core and possible structural role of pharmacological chaperones.

Misfolding of the prion protein is central to prion disease aetiology. Although understanding the dynamics of the native fold helps to decipher the conformational conversion mechanism, a complete depiction of distal but coupled prion protein sites common across species is lacking. To fill this gap, we used normal mode analysis and network analysis to examine a collection of prion protein structures deposited on the protein data bank. Our study identified a core of conserved residues that sustains the connectivity across the C-terminus of the prion protein. We propose how a well-characterized pharmacological chaperone may stabilize the fold. Also, we provide insight into the effect on the native fold of initial misfolding pathways identified by others using kinetics studies.

Computational studies on the design of NCI natural products as inhibitors to SARS-CoV-2 main protease.

The pandemic coronavirus disease (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in more than 5 million deaths globally. Currently there are no effective drugs available to treat COVID-19. The viral protease replication can be blocked by the inhibition of main protease that is encoded in polyprotein 1a and is therefore a potential protein target for drug discovery. We have carried out virtual screening of NCI natural compounds followed by molecular docking in order to identify hit molecules as probable SARS-CoV-2 main protease inhibitors. The molecular dynamics (MD) simulations of apo form in complex with N3, α-ketoamide and NCI natural products was used to validate the screened compounds. The MD simulations trajectories were analyzed using normal mode analysis and principal component analysis revealing dynamical nature of the protein. These findings aid in understanding the binding of natural products and molecular mechanisms of SARS-CoV-2 main protease inhibition.Communicated by Ramaswamy H. Sarma.

Scytonemin redox status in a filamentous cyanobacterium visualized by an excitation-laser-line-scanning spontaneous Raman scattering spectral microscope.

Some cyanobacteria produce a UVA-absorbing pigment, scytonemin, at extracellular sheaths. Although scytonemin-containing dark sheaths are recognizable through optical microscopes and its redox changes have been known for decades, there has been no report to obtain images of both oxidized and reduced scytonemins at subcellular resolution. Here, we show that a spontaneous Raman scattering spectral microscopy based on an excitation-laser-line-scanning method unveil 3D subcellular distributions of both the oxidized and reduced scytonemins in a filamentous cyanobacterium. The redox changes of scytonemin were supported by comparison in the Raman spectra between the cyanobacterial cells, solid-state scytonemin and quantum chemical normal mode analysis. Distributions of carotenoids, phycobilins, and the two redox forms of scytonemin were simultaneously visualized with an excitation wavelength at 1064 nm that is virtually free from the optical screening by the dark sheaths. The redox differentiation of scytonemin will advance our understanding of the redox homeostasis and secretion mechanisms of individual cyanobacteria as well as microscopic chemical environments in various microbial communities. The line-scanning Raman microscopy based on the 1064 nm excitation is thus a promising tool for exploring hitherto unreported Raman spectral features and distribution of nonfluorescent molecules embedded below nontransparent layers for visible light, while avoiding interference by autofluorescence.

Effect of phosphorylation on the structural dynamics, thermal stability of human dopamine transporter: A simulation study using normal modes, molecular dynamics and Markov State Model.

The Human Dopamine Transporter (hDAT) plays an essential role in modulating the Influx/Efflux of dopamine, and it is involved in the mechanism of certain neurodegenerative diseases such as Parkinson's disease. Several studies have reported important states for Dopamine transport: outward-facing open state (OFo), the outward-facing closed state (OFc), the holo-occluded state closed (holo), and the inward-facing open state (IFo). Furthermore, experimental assays have shown that different phosphorylation conditions in hDAT can affect the rate of dopamine absorption. We present a protocol using hybrid simulation methods to study the conformational dynamics and stability of states of hDAT under different phosphorylation sites. With this protocol, we explored the conformational space of hDAT, identified the states, and evaluated the free energy differences and the transition probabilities between them in each of the phosphorylation cases. We also presented the conformational changes and correlated them with those described in the literature. There is a thesis/hypothesis that the phosphorylation condition corresponding to NP-333 system (where all sites Ser/Thr from residue 2 to 62 and 254 to 613 are phosphorylated, except residue 333) would decrease the rate of dopamine transport from the extracellular medium to the intracellular medium by hDAT as previously described in the literature by Lin et al., 2003. Our results corroborated this thesis/hypothesis and the data reported. It is probably due to the affectation/changes/alteration of the conformational dynamics of this system that makes the intermediate states more likely and makes it difficult to initial states associated with the uptake of dopamine in the extracellular medium, corroborating the experimental results. Furthermore, our results showed that just single phosphorylation/dephosphorylation could alter intrinsic protein motions affecting the sampling of one or more states necessary for dopamine transport. In this sense, the modification of phosphorylation influences protein movements and conformational preferences, affecting the stability of states and the transition between them and, therefore, the transport.

Engineering Modified mRNA-Based Vaccine against Dengue Virus Using Computational and Reverse Vaccinology Approaches.

Dengue virus belonging to the family Flaviviridae and its four serotypes are responsible for dengue infections, which extend over 60 countries in tropical and subtropical areas of the world including Pakistan. During the ongoing dengue outbreak in Pakistan (2022), over 30,000 cases have been reported, and over 70 lives have been lost. The only commercialized vaccine against DENV, Dengvaxia, cannot be administered as a prophylactic measure to cure this infection due to various complications. Using machine learning and reverse vaccinology approaches, this study was designed to develop a tetravalent modified nucleotide mRNA vaccine using NS1, prM, and EIII sequences of dengue virus from Pakistani isolates. Based on high antigenicity, non-allergenicity, and toxicity profiling, B-cell epitope, cytotoxic T lymphocyte (CTL), and helper T lymphocyte (HTL) putative vaccine targets were predicted. Molecular docking confirmed favorable interactions between T-cell epitopes and their respective HLA alleles, while normal mode analysis validated high-affinity interactions of vaccine proteins with immune receptors. In silico immune simulations confirmed adequate immune responses to eliminate the antigen and generate memory. Codon optimization, physicochemical features, nucleotide modifications, and suitable vector availability further ensured better antigen expression and adaptive immune responses. We predict that this vaccine construct may prove to be a good vaccinal candidate against dengue virus in vitro as well.

Corrigendum: Deciphering the mechanism of inhibition of SERCA1a by sarcolipin using molecular simulations.

[This corrects the article DOI: 10.3389/fmolb.2020.606254.].

DeepHEMNMA: ResNet-based hybrid analysis of continuous conformational heterogeneity in cryo-EM single particle images.

Single-particle cryo-electron microscopy (cryo-EM) is a technique for biomolecular structure reconstruction from vitrified samples containing many copies of a biomolecular complex (known as single particles) at random unknown 3D orientations and positions. Cryo-EM allows reconstructing multiple conformations of the complexes from images of the same sample, which usually requires many rounds of 2D and 3D classifications to disentangle and interpret the combined conformational, orientational, and translational heterogeneity. The elucidation of different conformations is the key to understand molecular mechanisms behind the biological functions of the complexes and the key to novel drug discovery. Continuous conformational heterogeneity, due to gradual conformational transitions giving raise to many intermediate conformational states of the complexes, is both an obstacle for high-resolution 3D reconstruction of the conformational states and an opportunity to obtain information about multiple coexisting conformational states at once. HEMNMA method, specifically developed for analyzing continuous conformational heterogeneity in cryo-EM, determines the conformation, orientation, and position of the complex in each single particle image by image analysis using normal modes (the motion directions simulated for a given atomic structure or EM map), which in turn allows determining the full conformational space of the complex but at the price of high computational cost. In this article, we present a new method, referred to as DeepHEMNMA, which speeds up HEMNMA by combining it with a residual neural network (ResNet) based deep learning approach. The performance of DeepHEMNMA is shown using synthetic and experimental single particle images.

The Statistical Trends of Protein Evolution: A Lesson from AlphaFold Database.

The recent development of artificial intelligence provides us with new and powerful tools for studying the mysterious relationship between organism evolution and protein evolution. In this work, based on the AlphaFold Protein Structure Database (AlphaFold DB), we perform comparative analyses of the proteins of different organisms. The statistics of AlphaFold-predicted structures show that, for organisms with higher complexity, their constituent proteins will have larger radii of gyration, higher coil fractions, and slower vibrations, statistically. By conducting normal mode analysis and scaling analyses, we demonstrate that higher organismal complexity correlates with lower fractal dimensions in both the structure and dynamics of the constituent proteins, suggesting that higher functional specialization is associated with higher organismal complexity. We also uncover the topology and sequence bases of these correlations. As the organismal complexity increases, the residue contact networks of the constituent proteins will be more assortative, and these proteins will have a higher degree of hydrophilic-hydrophobic segregation in the sequences. Furthermore, by comparing the statistical structural proximity across the proteomes with the phylogenetic tree of homologous proteins, we show that, statistical structural proximity across the proteomes may indirectly reflect the phylogenetic proximity, indicating a statistical trend of protein evolution in parallel with organism evolution. This study provides new insights into how the diversity in the functionality of proteins increases and how the dimensionality of the manifold of protein dynamics reduces during evolution, contributing to the understanding of the origin and evolution of lives.