Effects of intrinsically disordered regions in gp120 underlying HIV neutralization phenotypes.

HIV envelope protein gp120 is considered a primary molecular determinant of viral neutralization phenotype due to its critical role in viral entry and immune evasion. The intrinsically disordered regions (IDRs) in gp120 are responsible for their extensive sequence variations and significant structural rearrangements. Despite HIV neutralization phenotype and sequence/structural information of gp120 have been experimentally characterized, there remains a gap in our understanding of the correlation between the viral phenotype and IDRs in gp120. Here, we combined machine learning (ML) techniques and molecular dynamics (MD) simulations to gain data-driven and molecule-mechanism insights into relationships between viral sequence, structure, and phenotypes from the perspective of IDRs in gp120. ML models, trained only on the length and disorder score of IDRs, achieved equivalent performance to the best baseline model using amino acid sequences to discriminate HIV neutralization phenotype, indicating that the lengths or disorder of specific IDRs are strongly related to HIV neutralization phenotypes. Comparative MD analysis reveals that gp120 with extreme neutralization phenotypes in multiple conformational states, especially some IDRs, exhibit significantly distinct structural dynamics, conformational flexibility, and thermodynamic distributions. Taken together, our study provided insights into the role of IDRs in gp120 responding to HIV neutralization phenotypes, which will advance the understanding of molecular mechanisms underlying viral function associated with HIV neutralization phenotype and help develop antiviral vaccines or drugs.

Characterization of a GH10 extremely thermophilic xylanase from the metagenome of hot spring for prebiotic production.

A xylanase gene (named xyngmqa) was identified from the metagenomic data of the Gumingquan hot spring (92.5 °C, pH 9.2) in Tengchong City, Yunnan Province, southwest China. It showed the highest amino acid sequence identity (82.70%) to endo-1,4-beta-xylanase from Thermotoga caldifontis. A constitutive expression plasmid (denominated pSHY211) and double-layer plate (DLP) method were constructed for cloning, expression, and identification of the XynGMQA gene. The XynGMQA gene was synthesized and successfully expressed in Escherichia coli DH5α. XynGMQA exhibited optimal activity at 90 °C and pH 4.6, being thermostable by maintaining 100% of its activity after 2 h incubated at 80 °C. Interestingly, its enzyme activity was enhanced by high temperatures (70 and 80 °C) and low pH (3.0-6.0). About 150% enzyme activity was detected after incubation at 70 °C for 20 to 60 min or 80 °C for 10 to 40 min, and more than 140% enzyme activity after incubation at pH 3.0 to 6.0 for 12 h. Hydrolytic products of beechwood xylan with XynGMQA were xylooligosaccharides, including xylobiose (X2), xylotriose (X3), and xylotetraose (X4). These properties suggest that XynGMQA as an extremely thermophilic xylanase, may be exploited for biofuel and prebiotic production from lignocellulosic biomass.

Structural determinants underlying high-temperature adaptation of thermophilic xylanase from hot-spring microorganisms.

Thermophilic xylanases from hot-spring microorganisms play potential biological and industrial applications for renewable and sustainable social development. However, high-temperature adaptation mechanisms of these thermophilic xylanases remain elusive at the molecular and evolutionary levels. Here, two recently reported xylanases, named XynDRTY1 and XynM1, from hot springs were subjected to molecular dynamics (MD) simulations at a series of temperature gradients and comparatively analyzed in comparison with the evolutionary background of the xylanase family. Comparative analysis of MD trajectories revealed that the XynM1 exhibits smaller structural dynamics and greater thermal stability than the XynDRTY1, although both share a similar fold architecture with structural differences in the ßα_loops. Local regions whose conformational flexibility and regular secondary structure exhibited differences as temperature increases were closely related to the high-temperature adaptation of xylanase, implying that stabilization of these regions is a feasible strategy to improve the thermal stability of xylanases. Furthermore, coevolutionary information from the xylanase family further specified the structural basis of xylanases. Thanks to these results about the sequence, structure, and dynamics of thermophilic xylanases from hot springs, a series of high-temperature-related structural determinants were resolved to promote understanding of the molecular mechanism of xylanase high-temperature adaptation and to provide direct assistance in the improvement of xylanase thermal stability.

Study on molecular mechanisms of CD4 dependency and independency of HIV-1 gp120.

Different HIV-1 strains have different antibody neutralization phenotypes (or CD4-dependencies). However, the molecular mechanisms underlying these differences remain to be elucidated. In this study, we constructed gp120 structural models from the CD4-dependent, neutralization-resistant JR-FL strain and the CD4-independent, neutralization-sensitive R2 strain and carried out several conventional molecular dynamics (MD) simulations and free energy landscape (FEL) constructions. Comparative analyses of the MD simulations and FELs indicated that R2 gp120 had higher global structural flexibility and greater conformational diversity than JR-FL gp120. This provides the preconditions for R2 gp120 to adopt a more open conformation than JR-FL gp120. Essential dynamics (ED) analysis showed that the collective motions of R2 gp120 tend towards an open state while those of JR-FL gp120 tend to retain a closed state. Based on conformational selection theory, R2 gp120's more readily sampled open state makes it more sensitive to neutralizing antibodies (or more CD4-independent) than JR-FL gp120, which may explain why the HIV-1 R2 and JR-FL strains show CD4-independent and -dependent phenotypes, respectively. Our study provides thermodynamic and kinetic insights into the CD4-dependent and -independent molecular mechanisms of HIV-1 gp120 and helps shed light on HIV-1 immune evasion.

Electrostatic Interactions Are the Primary Determinant of the Binding Affinity of SARS-CoV-2 Spike RBD to ACE2: A Computational Case Study of Omicron Variants.

To explore the mechanistic origin that determines the binding affinity of SARS-CoV-2 spike receptor binding domain (RBD) to human angiotensin converting enzyme 2 (ACE2), we constructed the homology models of RBD-ACE2 complexes of four Omicron subvariants (BA.1, BA.2, BA.3 and BA.4/5), and compared them with wild type complex (RBDWT-ACE2) in terms of various structural dynamic properties by molecular dynamics (MD) simulations and binding free energy (BFE) calculations. The results of MD simulations suggest that the RBDs of all the Omicron subvariants (RBDOMIs) feature increased global structural fluctuations when compared with RBDWT. Detailed comparison of BFE components reveals that the enhanced electrostatic attractive interactions are the main determinant of the higher ACE2-binding affinity of RBDOMIs than RBDWT, while the weakened electrostatic attractive interactions determine RBD of BA.4/5 subvariant (RBDBA.4/5) lowest ACE2-binding affinity among all Omicron subvariants. The per-residue BFE decompositions and the hydrogen bond (HB) networks analyses indicate that the enhanced electrostatic attractive interactions are mainly through gain/loss of the positively/negatively charged residues, and the formation or destruction of the interfacial HBs and salt bridges can also largely affect the ACE2-binding affinity of RBD. It is worth pointing out that since Q493R plays the most important positive contribution in enhancing binding affinity, the absence of this mutation in RBDBA.4/5 results in a significantly weaker binding affinity to ACE2 than other Omicron subvariants. Our results provide insight into the role of electrostatic interactions in determining of the binding affinity of SARS-CoV-2 RBD to human ACE2.

Deciphering gp120 sequence variation and structural dynamics in HIV neutralization phenotype by molecular dynamics simulations and graph machine learning.

Human immunodeficiency virus (HIV) exploits the sequence variation and structural dynamics of the envelope glycoprotein gp120 to evade the immune attack of neutralization antibodies, contributing to various HIV neutralization phenotypes. Although the HIV neutralization phenotype has been experimentally characterized, the roles of rapid sequence variability and significant structural dynamics of gp120 are not well understood. Here, 45 prefusion gp120 from different HIV strains belong to three tiers of sensitive, moderate, and resistant neutralization phenotype are structurally modeled by homology modeling and then investigated by molecular dynamics (MD) simulations and graph machine learning (ML). Our results show that the structural deviations, population distribution, and conformational flexibility of gp120 are related to the HIV neutralization phenotype. Per-residue dynamics indicate the local regions especially in the second structural elements with high-flexibility, may be responsible for the HIV neutralization phenotype. Moreover, a graph ML model with the attention mechanism was trained to explore inherent representation related to the classification of the HIV neutralization phenotype, further distinguishing the strong related gp120 sequence variation together with structural dynamics in the HIV neutralization phenotype. Our study not only deciphers gp120 sequence variation and structural dynamics in the HIV neutralization phenotype but also explores complex relationships between the sequence, structure, and dynamics of protein by combining MD simulations and ML.

Characterization of a Cu2+, SDS, alcohol and glucose tolerant GH1 ß-glucosidase from Bacillus sp. CGMCC 1.16541.

A ß-glucosidase gene (bsbgl1a) from Bacillus sp. CGMCC 1.16541 was expressed in Escherichia coli BL21 and subsequently characterized. The amino acid sequence shared 83.64% identity with ß-glucosidase (WP_066390903.1) from Fictibacillus phosphorivorans. The recombinant ß-glucosidase (BsBgl1A) had a molecular weight of 52.2 kDa and could hydrolyze cellobiose, cellotriose, cellotetrose, p-nitrophenyl-ß-D-glucopyranoside (pNPG), and p-nitrophenyl-ß-D-xylopyranoside (pNPX). Optimal activity for BsBgl1A was recorded at 45 °C with a pH between 5.6 and 7.6, and 100% of its activity was maintained after a 24 h incubation between pH 4 and 9. Kinetic characterization revealed an enzymatic turnover (Kcat) of 616 ± 2 s-1 (with cellobiose) and 3.5 ± 0.1 s-1 (with p-nitrophenyl-ß-D-glucopyranoside). Interestingly, the recombinant enzyme showed cupric ion (Cu2+), sodium dodecyl sulfate (SDS) and alcohol tolerance at 10 mM for Cu2+ and 10% for both SDS and alcohol. Additionally, BsBgl1A had high tolerance for glucose (Ki = 2095 mM), which is an extremely desirable feature for industrial applications. Following the addition of BsBgl1A (0.05 mg/ml) to a commercial cellulase reaction system, glucose yields from sugarcane bagasse increased 100% after 1 day at 45 °C. This work identifies a Cu2+, SDS, alcohol, and glucose tolerant GH1 ß-glucosidase with potential applications in the hydrolysis of cellulose for the bioenergy industry.

New Insight into Mechanisms of Protein Adaptation to High Temperatures: A Comparative Molecular Dynamics Simulation Study of Thermophilic and Mesophilic Subtilisin-Like Serine Proteases.

In high-temperature environments, thermophilic proteins must possess enhanced thermal stability in order to maintain their normal biological functions. However, the physicochemical basis of the structural stability of thermophilic proteins at high temperatures remains elusive. In this study, we performed comparative molecular dynamics simulations on thermophilic serine protease (THM) and its homologous mesophilic counterpart (PRK). The comparative analyses of dynamic structural and geometrical properties suggested that THM adopted a more compact conformation and exhibited more intramolecular interactions and lower global flexibility than PRK, which could be in favor of its thermal stability in high-temperature environments. Comparison between protein solvent interactions and the hydrophobicity of these two forms of serine proteases showed that THM had more burial of nonpolar areas, and less protein solvent hydrogen bonds (HBs), indicating that solvent entropy maximization and mobility may play a significant role in THM's adaption to high temperature environments. The constructed funnel-like free energy landscape (FEL) revealed that, in comparison to PRK, THM had a relatively flat and narrow free energy surface, and a lower minimum free energy level, suggesting that the thermophilic form had lower conformational diversity and flexibility. Combining the FEL theory and our simulation results, we conclude that the solvent (entropy force) plays a significant role in protein adaption at high temperatures.

CD4-binding obstacles in conformational transitions and allosteric communications of HIV gp120.

As the only exposed viral protein at the membrane surface of HIV, envelope glycoprotein gp120 is responsible for recognizing host cells and mediating virus-cell membrane fusion. Available structures of gp120 indicate that it exhibits two distinct conformational states, called closed and open states. Although experimental data demonstrates that CD4 binding stabilizes open state of gp120, detailed structural dynamics and kinetics of gp120 during this process remain elusive. Here, two open-state gp120 simulation systems, one without any ligands (ligand-free) and the other complexed with CD4 (CD4-bound), were subjected to microsecond-scale molecular dynamics simulations following the conformational transitions and allosteric pathways of gp120 evaluated by using the Markov state model and a network-based method, respectively. Our results provide an atomic-resolution description of gp120 conformational transitions, suggesting that gp120 is intrinsically dynamic from the open state to closed state, whereas CD4 binding blocks these transitions. Consistent with experimental structures, five metastable conformations with different orientations of the V1/V2 region and V3 loop have been extracted. The binding of CD4 significantly enhances allosteric communications from the CD4-binding site to V3 loop and ß20-21 hairpin, resulting in high-affinity interactions with coreceptors and activation of the conformational transitions switcher, respectively. This study will facilitate the structural understanding of the CD4-binding effects on conformational transitions and allosteric pathways of gp120.

Anti-HIV drug repurposing against SARS-CoV-2.

A novel severe acute respiratory syndrome human coronavirus (SARS HCoV) was identified from respiratory illness patients (named SARS-CoV-2 by ICTV) in December 2019 and has recently emerged as a serious threat to world public health. However, no approved drugs have been found to effectively inhibit the virus. Since it has been reported that HIV protease inhibitors can be used as anti-SARS drugs by targeting SARS-CoV-1 3CLpro, we chose six approved anti-HIV drugs and investigated their binding interactions with 3CLpro to evaluate their potential to become clinical drugs for the new coronavirus pneumonia (COVID-19) caused by SARS-CoV-2 infection. The molecular docking results indicate that the 3CLpro of SARS-CoV-2 has a higher binding affinity for all the studied inhibitors than does SARS-CoV-1. Two docking complexes (indinavir and darunavir) with high docking scores were further subjected to MM-PBSA binding free energy calculations to detail the molecular interactions between these two protease inhibitors and SARS HCoV 3CLpro. Our results show that, among the inhibitors tested, darunavir has the highest binding affinity with SARS-CoV-2 and SARS-CoV-1 3CLpro, indicating that it may have the potential to be used as an anti-COVID-19 clinical drug. The mechanism behind the increased binding affinity of HIV protease inhibitors toward SARS-CoV-2 3CLpro (as compared to SARS-CoV-1) was investigated by MD simulations. Our study provides insight into the possible role of structural flexibility during interactions between SARS HCoV 3CLpro and inhibitors and sheds light on structure-based design of anti-COVID-19 drugs targeting SARS-CoV-2 3CLpro.

Probing intrinsic dynamics and conformational transition of HIV gp120 by molecular dynamics simulation.

The HIV envelope glycoprotein gp120 has evolved two distinct conformational states to balance viral infection and immune escape. One is a closed state resistant to most neutralization antibodies, and the other is an open state responsible for the binding of the receptor and coreceptors. Although the structures of gp120 in these two conformational states have been determined, a detailed molecular mechanism involving intrinsic dynamics and conformational transition is still elusive. In this study, µs-scale molecular dynamics simulation is performed to probe molecular dynamics and conformational transition away from the open state and approach the closed state. Our results reveal that open gp120 shows a larger structural deviation, higher conformational flexibility, and more conformational diversity than the form in the closed state, providing a structural explanation for receptor or coreceptor affinity at the open state and the neutralization resistance of closed conformation. Seven regions with greatly decreased coupled motions in the open states have been observed by dynamic cross-correlation analysis, indicating that conformational transition can be mainly attributed to the relaxation of intrinsic dynamics. Three conformations characterized by the structural orientations of the V1/V2 region and the V3 loop, suggesting gp120 is intrinsically dynamic from the open state to the closed state. Taken together, these findings shed light on the understanding of the conformational control mechanism of HIV.

Correction: Probing intrinsic dynamics and conformational transition of HIV gp120 by molecular dynamics simulation.

[This corrects the article DOI: 10.1039/D0RA06416E.].

Insight derived from molecular dynamics simulation into cold-adaptation mechanism of trypsins.

Communicated by Ramaswamy H. Sarma.

Effects of CD4 Binding on Conformational Dynamics, Molecular Motions, and Thermodynamics of HIV-1 gp120.

Human immunodeficiency virus type-1 (HIV-1) infection is triggered by its envelope (Env) glycoprotein gp120 binding to the host-cell receptor CD4. Although structures of Env/gp120 in the liganded state are known, detailed information about dynamics of the liganded gp120 has remained elusive. Two structural models, the CD4-free gp120 and the gp120-CD4 complex, were subjected to µs-scale multiple-replica molecular dynamics (MD) simulations to probe the effects of CD4 binding on the conformational dynamics, molecular motions, and thermodynamics of gp120. Comparative analyses of MD trajectories in terms of structural deviation and conformational flexibility reveal that CD4 binding effectively suppresses the overall conformational fluctuations of gp120. Despite the largest fluctuation amplitude of the V1/V2 region in both forms of gp120, the presence of CD4 prevents it from approaching the gp120 core. Comparison of the constructed free energy landscapes (FELs) shows that CD4 binding reduces the conformational entropy and conformational diversity while enhancing the stability of gp120. Further comparison of the representative structures extracted from free energy basins/minima of FELs reveals that CD4 binding weakens the reorientation ability of V1/V2 and hence hinders gp120 from transitioning out of the liganded state to the unliganded state. Therefore, locking gp120 conformation via restraining V1/V2 reorientation with small molecules seems to be a promising strategy to control HIV-1 infection. Our computer simulation results support the conformational selection mechanism for CD4 binding to gp120 and facilitate the understanding of HIV-1 immune evasion mechanisms.

Insight derived from molecular dynamics simulation into dynamics and molecular motions of cuticle-degrading serine protease Ver112.

Cuticle-degrading serine protease Ver112, which derived from a nematophagous fungus Lecanicillium psalliotae, has been exhibited to have high cuticle-degrading and nematicidal activities. We have performed molecular dynamics (MD) simulation based on the crystal structure of Ver112 to investigate its dynamic properties and large-scale concerted motions. The results indicate that the structural core of Ver112 shows a small fluctuation amplitude, whereas the substrate binding sites, and the regions close to and opposite the substrate binding sites experience significant conformational fluctuations. The large concerted motions obtained from essential dynamics (ED) analysis of MD trajectory can lead to open or close of the substrate binding sites, which are proposed to be linked to the functional properties of Ver112, such as substrate binding, orientation, catalytic, and release. The significant motion in the loop regions that is located opposite the binding sites are considered to play an important role in modulating the dynamics of the substrate binding sites. Furthermore, the bottom of free energy landscape (FEL) of Ver112 are rugged, which is mainly caused by the fluctuations of substrate binding regions and loops located opposite the binding site. In addition, the mechanism underlying the high flexibility and catalytic activity of Ver112 was also discussed. Our simulation study complements the biochemical and structural studies, and provides insight into the dynamics-function relationship of cuticle-degrading serine protease Ver112.

The effects of maternal cigarette smoking on pregnancy outcomes using assisted reproduction technologies: An updated meta-analysis.

There is strong evidence indicating that smoking has negative effects on female reproductive health. Studies to investigate the effects of female smoking on IVF outcomes have been conducted by several research groups, yet the results are controversial. To evaluate the impacts of female smoking on the outcomes of assisted reproduction, a meta-analysis was performed, which included studies published in English up to September 6, 2017 from the MEDLINE, EMBASE, and Cochrane library databases. Twenty-eight studies encompassing 5009 female smokers seeking assisted reproduction and 10,078 non-smokers were used in this meta-analysis. Significant negative outcomes were detected in the female smokers compared with non-smokers including decreases in live birth rate per cycle (OR=0.52, 95% CI 0.37-0.74), in clinical pregnancy rate per cycle (OR 0.59, 95% CI 0.51-0.68), in number of retrieved oocytes (MD=-0.87, 95% CI -1.39 to -0.25), and in average fertilization rate (MD=-4.80, 95% CI -8.49 to -2.02), as well as a significantly increased miscarriage rate per pregnancy (OR=2.48, 95% CI 1.79-3.43). In conclusion, the current meta-analysis provides compelling evidence that female smoking has a significantly negative impact on the outcomes of assisted reproductive technology (ART) and strongly recommends that female smokers will greatly benefit from a smoking cessation before employing ART to become pregnant.

Insights into the role of electrostatics in temperature adaptation: a comparative study of psychrophilic, mesophilic, and thermophilic subtilisin-like serine proteases.

To investigate the role of electrostatics in different temperature adaptations, we performed a comparative study on subtilisin-like serine proteases from psychrophilic Vibrio sp. PA-44 (VPR), mesophilic Engyodontium album (Tritirachium album) (PRK), and thermophilic Thermus aquaticus (AQN) using multiple-replica molecular dynamics (MD) simulations combined with continuum electrostatics calculations. The results reveal that although salt bridges are not a crucial factor in determining the overall thermostability of these three proteases, they on average provide the greatest, moderate, and least electrostatic stabilization to AQN, PRK, and VPR, respectively, at the respective organism growth temperatures. Most salt bridges in AQN are effectively stabilizing and thus contribute to maintaining the overall structural stability at 343 K, while nearly half of the salt bridges in VPR interconvert between being stabilizing and being destabilizing, likely aiding in enhancing the local conformational flexibility at 283 K. The individual salt bridges, salt-bridge networks, and calcium ions contribute differentially to local stability and flexibility of these three enzyme structures, depending on their spatial distributions and electrostatic strengths. The shared negatively charged surface potential at the active center of the three enzymes may provide the active-center flexibility necessary for nucleophilic attack and proton transfer. The differences in distributions of the electro-negative, electro-positive, and electro-neutral potentials, particularly over the back surfaces of the three proteases, may modulate/affect not only protein solubility and thermostability but also structural stability and flexibility/rigidity. These results demonstrate that electrostatics contributes to both heat and cold adaptation of subtilisin-like serine proteases through fine-tuning, either globally or locally, the structural stability and conformational flexibility/rigidity, thus providing a foundation for further engineering and mutagenesis studies.

Effect of the Solvent Temperatures on Dynamics of Serine Protease Proteinase K.

To obtain detailed information about the effect of the solvent temperatures on protein dynamics, multiple long molecular dynamics (MD) simulations of serine protease proteinase K with the solute and solvent coupled to different temperatures (either 300 or 180 K) have been performed. Comparative analyses demonstrate that the internal flexibility and mobility of proteinase K are strongly dependent on the solvent temperatures but weakly on the protein temperatures. The constructed free energy landscapes (FELs) at the high solvent temperatures exhibit a more rugged surface, broader spanning range, and higher minimum free energy level than do those at the low solvent temperatures. Comparison between the dynamic hydrogen bond (HB) numbers reveals that the high solvent temperatures intensify the competitive HB interactions between water molecules and protein surface atoms, and this in turn exacerbates the competitive HB interactions between protein internal atoms, thus enhancing the conformational flexibility and facilitating the collective motions of the protein. A refined FEL model was proposed to explain the role of the solvent mobility in facilitating the cascade amplification of microscopic motions of atoms and atomic groups into the global collective motions of the protein.

Insight derived from molecular dynamics simulations into molecular motions, thermodynamics and kinetics of HIV-1 gp120.

Although the crystal structures of the HIV-1 gp120 core bound and pre-bound by CD4 are known, the details of dynamics involved in conformational equilibrium and transition in relation to gp120 function have remained elusive. The homology models of gp120 comprising the N- and C-termini and loops V3 and V4 in the CD4-bound and CD4-unbound states were built and subjected to molecular dynamics (MD) simulations to investigate the differences in dynamic properties and molecular motions between them. The results indicate that the CD4-bound gp120 adopted a more compact and stable conformation than the unbound form during simulations. For both the unbound and bound gp120, the large concerted motions derived from essential dynamics (ED) analyses can influence the size/shape of the ligand-binding channel/cavity of gp120 and, therefore, were related to its functional properties. The differences in motion direction between certain structural components of these two forms of gp120 were related to the conformational interconversion between them. The free energy calculations based on the metadynamics simulations reveal a more rugged and complex free energy landscape (FEL) for the unbound than for the bound gp120, implying that gp120 has a richer conformational diversity in the unbound form. The estimated free energy difference of ∼-6.0 kJ/mol between the global minimum free energy states of the unbound and bound gp120 indicates that gp120 can transform spontaneously from the unbound to bound states, revealing that the bound state represents a high-probability "ground state" for gp120 and explaining why the unbound state resists crystallization. Our results provide insight into the dynamics-and-function relationship of gp120, and facilitate understandings of the thermodynamics, kinetics and conformational control mechanism of HIV-1 gp120.

Protein dynamics and motions in relation to their functions: several case studies and the underlying mechanisms.

Proteins are dynamic entities in cellular solution with functions governed essentially by their dynamic personalities. We review several dynamics studies on serine protease proteinase K and HIV-1 gp120 envelope glycoprotein to demonstrate the importance of investigating the dynamic behaviors and molecular motions for a complete understanding of their structure-function relationships. Using computer simulations and essential dynamic (ED) analysis approaches, the dynamics data obtained revealed that: (i) proteinase K has highly flexible substrate-binding site, thus supporting the induced-fit or conformational selection mechanism of substrate binding; (ii) Ca(2+) removal from proteinase K increases the global conformational flexibility, decreases the local flexibility of substrate-binding region, and does not influence the thermal motion of catalytic triad, thus explaining the experimentally determined decreased thermal stability, reduced substrate affinity, and almost unchanged catalytic activity upon Ca(2+) removal; (iii) substrate binding affects the large concerted motions of proteinase K, and the resulting dynamic pocket can be connected to substrate binding, orientation, and product release; (iv) amino acid mutations 375 S/W and 423 I/P of HIV-1 gp120 have distinct effects on molecular motions of gp120, facilitating 375 S/W mutant to assume the CD4-bound conformation, while 423 I/P mutant to prefer for CD4-unliganded state. The mechanisms underlying protein dynamics and protein-ligand binding, including the concept of the free energy landscape (FEL) of the protein-solvent system, how the ruggedness and variability of FEL determine protein's dynamics, and how the three ligand-binding models, the lock-and-key, induced-fit, and conformational selection are rationalized based on the FEL theory are discussed in depth.