Discontinuous instabilities in disordered solids.

Under a sufficiently large load, a solid material will flow via rearrangements, where particles change neighbors. Such plasticity is most easily described in the athermal, quasistatic limit of zero temperature and infinitesimal loading rate, where rearrangements occur only when the system becomes mechanically unstable. For disordered solids, the instabilities marking the onset of rearrangements have long been believed to be fold instabilities, in which an energy barrier disappears and the frequency of a normal mode of vibration vanishes continuously. Here, we report that there exists another, anomalous, type of instability caused by the breaking of a "stabilizing bond," whose removal creates an unstable vibrational mode. For commonly studied systems, such as those with harmonic finite-range interparticle interactions, such "discontinuous instabilities" are not only inevitable, they often dominate the modes of failure. Stabilizing bonds are a subset of all the bonds in the system and are prevalent in disordered solids generally. Although they do not trigger discontinuous instabilities in systems with vanishing stiffness at the interaction cutoff, they are, even in those cases, local indicators of incipient mechanical failure. They therefore provide an accurate structural predictor of instabilities not only of the discontinuous type but of the fold type as well.

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.

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.

Tunnel-like region observed as a potential allosteric site in Staphylococcus aureus Glyceraldehyde-3-phosphate dehydrogenase.

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) catalyzing the sixth step of glycolysis has been investigated for allosteric features that might be used as potential target for specific inhibition of Staphylococcus aureus (S.aureus). X-ray structure of bacterial enzyme for which a tunnel-like opening passing through the center previously proposed as an allosteric site has been subjected to six independent 500 ns long Molecular Dynamics simulations. Harmonic bond restraints were employed at key residues to underline the allosteric feature of this region. A noticeable reduction was observed in the mobility of NAD+ binding domains when restrictions were applied. Also, a substantial decrease in cross-correlations between distant Cα fluctuations was detected throughout the structure. Mutual information (MI) analysis revealed a similar decrease in the degree of correspondence in positional fluctuations in all directions everywhere in the receptor. MI between backbone and side-chain torsional variations changed its distribution profile and decreased considerably around the catalytic sites when restraints were employed. Principal component analysis clearly showed that the restrained state sampled a narrower range of conformations than apo state, especially in the first principal mode due to restriction in the conformational flexibility of NAD+ binding domain. Clustering the trajectory based on catalytic site residues displayed a smaller repertoire of conformations for restrained state compared to apo. Representative snapshots subjected to k-shortest pathway analysis revealed the impact of bond restraints on the allosteric communication which displayed distinct optimal and suboptimal pathways for two states, where observed frequencies of critical residues Gln51 and Val283 at the proposed site changed considerably.

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.

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.

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.

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.

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.

Low-frequency collective motion of DNA-binding domain defines p53 function.

Most mutations in the DNA-binding domain (DBD) of p53 inactivate or rescue the protein function interacting with the minor groove of DNA. However, how the conformation changes propagating from the mutation sites result in distinct molecular recognition is still not well understood. As the protein mobility is an intrinsic property encrypted in its primary structure, we examined if different structures of wild-type and mutant p53 core domains display any unique patterns of intrinsic mobility. Normal mode calculation was employed to characterize the collective dynamics of DBD in p53 monomer and tetramer as well as their mutants. Intriguingly, the low-frequency collective motions of DBD show similar patterns between the wild-type protein and the rescued mutants. The analysis on atomic backbone fluctuations and low-frequency vibration mode statistics does further support the correlation between the intrinsic collective motion of DBD and the p53 protein function. The mutations in the DBD influence the low-frequency vibration of the p53 tetramer via the change of the collective motions among its four monomers. These findings thus provide new insights for understanding the physical mechanism of p53 protein structure-function relationship and help find the small molecule drug to modulate protein dynamic for disease therapy.

Pseudokinases repurpose flexibility signatures associated with the protein kinase fold for noncatalytic roles.

The bilobal protein kinase-like fold in pseudokinases lack one or more catalytic residues, conserved in canonical protein kinases, and are considered enzymatically deficient. Tertiary structures of pseudokinases reveal that their loops topologically equivalent to activation segments of kinases adopt contracted configurations, which is typically extended in active conformation of kinases. Herein, anisotropic network model based normal mode analysis (NMA) was conducted on 51 active conformation structures of protein kinases and 26 crystal structures of pseudokinases. Our observations indicate that although backbone fluctuation profiles are similar for individual kinase-pseudokinase families, low intensity mean square fluctuations in pseudo-activation segment and other sub-structures impart rigidity to pseudokinases. Analyses of collective motions from functional modes reveal that pseudokinases, compared to active kinases, undergo distinct conformational transitions using the same structural fold. All-atom NMA of protein kinase-pseudokinase pairs from each family, sharing high amino acid sequence identities, yielded distinct community clusters, partitioned by residues exhibiting highly correlated fluctuations. It appears that atomic fluctuations from equivalent activation segments guide community membership and network topologies for respective kinase and pseudokinase. Our findings indicate that such adaptations in backbone and side-chain fluctuations render pseudokinases competent for catalysis-independent roles.

The effect of dimerization and ligand binding on the dynamics of Kaposi's sarcoma-associated herpesvirus protease.

The Kaposi's sarcoma-associated herpesvirus protease is essential for virus maturation. This protease functions under allosteric regulation that establishes its enzymatic activity upon dimerization. It exists in equilibrium between an inactive monomeric state and an active, weakly associating, dimeric state that is stabilized upon ligand binding. The dynamics of the protease dimer and its monomer were studied using the Gaussian network model and the anisotropic network model , and its role in mediating the allosteric regulation is demonstrated. We show that the dimer is composed of five dynamical domains. The central domain is formed upon dimerization and composed of helix five of each monomer, in addition to proximal and distal domains of each monomer. Dimerization reduces the mobility of the central domains and increases the mobility of the distal domains, in particular the binding site within them. The three slowest ANM modes of the dimer assist the protease in ligand binding, motion of the conserved Arg142 and Arg143 toward the oxyanion, and reducing the activation barrier for the tetrahedral transition state by stretching the bond that is cleaved by the protease. In addition, we show that ligand binding reduces the motion of helices α1 and α5 at the interface and explain how ligand binding can stabilize the dimer.

Analysis of changes of cavity volumes in predefined directions of protein motions and cavity flexibility.

Dynamics of protein cavities associated with protein fluctuations and conformational plasticity is essential for their biological function. NMR ensembles, molecular dynamics (MD) simulations, and normal mode analysis (NMA) provide appropriate frameworks to explore functionally relevant protein dynamics and cavity changes relationships. Within this context, we have recently developed analysis of null areas (ANA), an efficient method to calculate cavity volumes. ANA is based on a combination of algorithms that guarantees its robustness against numerical differentiations. This is a unique feature with respect to other methods. Herein, we present an updated and improved version that expands it use to quantify changes in cavity features, like volume and flexibility, due to protein structural distortions performed on predefined biologically relevant directions, for example, directions of largest contribution to protein fluctuations (principal component analysis [PCA modes]) obtained by MD simulations or ensembles of NMR structures, collective NMA modes or any other direction of motion associated with specific conformational changes. A web page has been developed where its facilities are explained in detail. First, we show that ANA can be useful to explore gradual changes of cavity volume and flexibility associated with protein ligand binding. Secondly, we perform a comparison study of the extent of variability between protein backbone structural distortions, and changes in cavity volumes and flexibilities evaluated for an ensemble of NMR active and inactive conformers of the epidermal growth factor receptor structures. Finally, we compare changes in size and flexibility between sets of NMR structures for different homologous chains of dynein.

Integration of network models and evolutionary analysis into high-throughput modeling of protein dynamics and allosteric regulation: theory, tools and applications.

Proteins are dynamical entities that undergo a plethora of conformational changes, accomplishing their biological functions. Molecular dynamics simulation and normal mode analysis methods have become the gold standard for studying protein dynamics, analyzing molecular mechanism and allosteric regulation of biological systems. The enormous amount of the ensemble-based experimental and computational data on protein structure and dynamics has presented a major challenge for the high-throughput modeling of protein regulation and molecular mechanisms. In parallel, bioinformatics and systems biology approaches including genomic analysis, coevolution and network-based modeling have provided an array of powerful tools that complemented and enriched biophysical insights by enabling high-throughput analysis of biological data and dissection of global molecular signatures underlying mechanisms of protein function and interactions in the cellular environment. These developments have provided a powerful interdisciplinary framework for quantifying the relationships between protein dynamics and allosteric regulation, allowing for high-throughput modeling and engineering of molecular mechanisms. Here, we review fundamental advances in protein dynamics, network theory and coevolutionary analysis that have provided foundation for rapidly growing computational tools for modeling of allosteric regulation. We discuss recent developments in these interdisciplinary areas bridging computational biophysics and network biology, focusing on promising applications in allosteric regulations, including the investigation of allosteric communication pathways, protein-DNA/RNA interactions and disease mutations in genomic medicine. We conclude by formulating and discussing future directions and potential challenges facing quantitative computational investigations of allosteric regulatory mechanisms in protein systems.

LARMD: integration of bioinformatic resources to profile ligand-driven protein dynamics with a case on the activation of estrogen receptor.

Protein dynamics is central to all biological processes, including signal transduction, cellular regulation and biological catalysis. Among them, in-depth exploration of ligand-driven protein dynamics contributes to an optimal understanding of protein function, which is particularly relevant to drug discovery. Hence, a wide range of computational tools have been designed to investigate the important dynamic information in proteins. However, performing and analyzing protein dynamics is still challenging due to the complicated operation steps, giving rise to great difficulty, especially for nonexperts. Moreover, there is a lack of web protocol to provide online facility to investigate and visualize ligand-driven protein dynamics. To this end, in this study, we integrated several bioinformatic tools to develop a protocol, named Ligand and Receptor Molecular Dynamics (LARMD, http://chemyang.ccnu.edu.cn/ccb/server/LARMD/ and http://agroda.gzu.edu.cn:9999/ccb/server/LARMD/), for profiling ligand-driven protein dynamics. To be specific, estrogen receptor (ER) was used as a case to reveal ERß-selective mechanism, which plays a vital role in the treatment of inflammatory diseases and many types of cancers in clinical practice. Two different residues (Ile373/Met421 and Met336/Leu384) in the pocket of ERß/ERα were the significant determinants for selectivity, especially Met336 of ERß. The helix H8, helix H11 and H7-H8 loop influenced the migration of selective agonist (WAY-244). These computational results were consistent with the experimental results. Therefore, LARMD provides a user-friendly online protocol to study the dynamic property of protein and to design new ligand or site-directed mutagenesis.

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.

Variation of Structural and Dynamical Flexibility of Myelin Basic Protein in Response to Guanidinium Chloride.

Myelin basic protein (MBP) is intrinsically disordered in solution and is considered as a conformationally flexible biomacromolecule. Here, we present a study on perturbation of MBP structure and dynamics by the denaturant guanidinium chloride (GndCl) using small-angle scattering and neutron spin-echo spectroscopy (NSE). A concentration of 0.2 M GndCl causes charge screening in MBP resulting in a compact, but still disordered protein conformation, while GndCl concentrations above 1 M lead to structural expansion and swelling of MBP. NSE data of MBP were analyzed using the Zimm model with internal friction (ZIF) and normal mode (NM) analysis. A significant contribution of internal friction was found in compact states of MBP that approaches a non-vanishing internal friction relaxation time of approximately 40 ns at high GndCl concentrations. NM analysis demonstrates that the relaxation rates of internal modes of MBP remain unaffected by GndCl, while structural expansion due to GndCl results in increased amplitudes of internal motions. Within the model of the Brownian oscillator our observations can be rationalized by a loss of friction within the protein due to structural expansion. Our study highlights the intimate coupling of structural and dynamical plasticity of MBP, and its fundamental difference to the behavior of ideal polymers in solution.

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex.

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.

Using normal mode analysis on protein structural models. How far can we go on our predictions?

Normal mode analysis (NMA) is a fast and inexpensive approach that is largely used to gain insight into functional protein motions, and more recently to create conformations for further computational studies. However, when the protein structure is unknown, the use of computational models is necessary. Here, we analyze the capacity of NMA in internal coordinate space to predict protein motion, its intrinsic flexibility, and atomic displacements, using protein models instead of native structures, and the possibility to use it for model refinement. Our results show that NMA is quite insensitive to modeling errors, but that calculations are strictly reliable only for very accurate models. Our study also suggests that internal NMA is a more suitable tool for the improvement of structural models, and for integrating them with experimental data or in other computational techniques, such as protein docking or more refined molecular dynamics simulations.

Predicting cryptic ligand binding sites based on normal modes guided conformational sampling.

To greatly expand the druggable genome, fast and accurate predictions of cryptic sites for small molecules binding in target proteins are in high demand. In this study, we have developed a fast and simple conformational sampling scheme guided by normal modes solved from the coarse-grained elastic models followed by atomistic backbone refinement and side-chain repacking. Despite the observations of complex and diverse conformational changes associated with ligand binding, we found that simply sampling along each of the lowest 30 modes is near optimal for adequately restructuring cryptic sites so they can be detected by existing pocket finding programs like fpocket and concavity. We further trained machine-learning protocols to optimize the combination of the sampling-enhanced pocket scores with other dynamic and conservation scores, which only slightly improved the performance. As assessed based on a training set of 84 known cryptic sites and a test set of 14 proteins, our method achieved high accuracy of prediction (with area under the receiver operating characteristic curve >0.8) comparable to the CryptoSite server. Compared with CryptoSite and other methods based on extensive molecular dynamics simulation, our method is much faster (1-2 hours for an average-size protein) and simpler (using only pocket scores), so it is suitable for high-throughput processing of large datasets of protein structures at the genome scale.