Shu Qun Liu

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) Ca2+ 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 Ca2+ 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.

Characterizing structural features of cuticle-degrading proteases from fungi by molecular modeling

BackgroundSerine proteases secreted by nematode and insect pathogenic fungi are bio-control agents which have commercial potential for developing into effective bio-pesticides. A thorough understanding of the structural and functional features of these proteases would significantly assist with targeting the design of efficient bio-control agents.ResultsStructural models of serine proteases PR1 from entomophagous fungus, Ver112 and VCP1 from nematophagous fungi, have been modeled using the homology modeling technique based on the crystal coordinate of the proteinase K. In combination with multiple sequence alignment, these models suggest one similar calcium-binding site and two common disulfide bridges in the three cuticle-degrading enzymes. In addition, the predicted models of the three cuticle-degrading enzymes present an essentially identical backbone topology and similar geometric properties with the exception of a limited number of sites exhibiting relatively large local conformational differences only in some surface loops and the N-, C termini. However, they differ from each other in the electrostatic surface potential, in hydrophobicity and size of the S4 substrate-binding pocket, and in the number and distribution of hydrogen bonds and salt bridges within regions that are part of or in close proximity to the S2-loop.ConclusionThese differences likely lead to variations in substrate specificity and catalytic efficiency among the three enzymes. Amino acid polymorphisms in cuticle-degrading enzymes were discussed with respect to functional effects and host preference. It is hoped that these structural models would provide a further basis for exploitation of these serine proteases from pathogenic fungi as effective bio-control agents.

Molecular motions of human HIV-1 gp120 envelope glycoproteins

The HIV-1 gp120 exterior envelope glycoprotein undergoes a series of conformational rearrangements while sequentially interacting with the receptor CD4 and the coreceptor CCR5 or CXCR4 on the surface of host cells to initiate virus entry. Both the crystal structures of the HIV-1 gp120 core bound by CD4 and antigen 17b, and the SIV gp120 core pre-bound by CD4 are known. Despite the wealth of knowledge on these static snapshots of molecular conformations, the details of molecular motions crucial to intervention remain elusive. We presented a comprehensive comparative analysis of dynamic behavior of gp120 in its CD4-complexed, CD4-free and CD4-unliganded states based on the homology models with modeled V3 and V4 loops. CONCOORD computer simulation was utilized to generate ensembles of feasible protein structures, which were subsequently analyzed by essential dynamics technique to identify preferred concerted motions. The revealed collective fluctuations are dominated by complex motional modes such as rotation/twisting, flexing/closing, and shortness/elongation between or within the inner, outer, and bridging-sheet domains. An attempt has been made to relate these modes to receptor/coreceptor association and neutralization avoidance. Covariance web analysis revealed four subdomains that undergo concerted motion in gp120. The structural components in gp120 that move in concert with CD4 were also identified, which may be the suitable target for inhibitor design to interrupt CD4-gp120 interaction. The differences in B-factors between the three gp120 states revealed certain structural regions that could be related either to CD4 association or to subsequent dissociation of gp120 from gp41. These dynamics data provide new insights into the structure-function relationship of gp120 and may aid in structure-based anti-HIV vaccine design.

Insights derived from molecular dynamics simulation into the molecular motions of serine protease proteinase K

Serine protease proteinase K, a member of the subtilisin family of enzymes, is of significant industrial, agricultural and biotechnological importance. Despite the wealth of structural information about proteinase K provided by static X-ray structures, a full understanding of the enzymatic mechanism requires further insight into the dynamic properties of this enzyme. Molecular dynamics simulations and essential dynamics (ED) analysis were performed to investigate the molecular motions in proteinase K. The results indicate that the internal core of proteinase K is relatively rigid, whereas the surface-exposed loops, most notably the substrate-binding regions, exhibit considerable conformational fluctuations. Further ED analysis reveals that the large concerted motions in the substrate-binding regions cause opening/closing of the substrate-binding pockets, thus supporting the proposed induced-fit mechanism of substrate binding. The distinct electrostatic/hydrogen-bonding interactions between Asp39 and His69 and between His69 and Ser224 within the catalytic triad lead to different thermal motions and orientations of these three catalytic residues, which can be related to their different functional roles in the catalytic process. Statistical analyses of the geometrical/functional properties as well as evolutionary conservation of the glycines in proteinase K-like proteins reveal that glycines may play an important role in determining the folding architecture and structural flexibility of this class of enzymes. Our simulation study complements the biochemical and structural studies and provides new insights into the dynamic structural basis of the functional properties of this class of enzymes.

Bacteria can mobilize nematode-trapping fungi to kill nematodes

In their natural habitat, bacteria are consumed by bacterivorous nematodes; however, they are not simply passive preys. Here we report a defensive mechanism used by certain bacteria to mobilize nematode-trapping fungi to kill nematodes. These bacteria release urea, which triggers a lifestyle switch in the fungus Arthrobotrys oligospora from saprophytic to nematode–predatory form; this predacious form is characterized by formation of specialized cellular structures or ‘traps’. The bacteria significantly promote the elimination of nematodes by A. oligospora. Disruption of genes involved in urea transport and metabolism in A. oligospora abolishes the urea-induced trap formation. Furthermore, the urea metabolite ammonia functions as a signal molecule in the fungus to initiate the lifestyle switch to form trap structures. Our findings highlight the importance of multiple predator–prey interactions in prey defense mechanisms.

The effect of calciums on molecular motions of proteinase K

The native serine protease proteinase K binds two calcium cations. It has been reported that Ca2+ removal decreased the enzyme’s thermal stability and to some extent the substrate affinity, but has discrepant effects on catalytic activity of the enzyme. Molecular dynamics simulations were performed on the Ca2+-bound and Ca2+-free proteases to investigate the mechanism by which the calciums affect the structural stability, molecular motions, and catalytic activity of proteinase K. Very similar structural properties were observed between these two forms of proteinase K during simulations; and several long-lived hydrogen bonds and salt bridges common to both forms of proteinase K were found to be crucial in maintaining the local conformations around these two Ca2+ sites. Although Ca2+ removal enhanced the overall flexibility of proteinase K, the flexibility in a limited number of segments surrounding the substrate-binding pockets decreased. The largest differences in the equilibrium structures of the two simulations indicate that, upon the removal of Ca2+, the large concerted motion originating from the Ca1 site can transmit to the substrate-binding regions but not to the catalytic triad residues. In conjunction with the large overlap of the essential subspaces between the two simulations, these results not only provide insight into the dynamics of the underlying molecular mechanism responsible for the unchanged enzymatic activity as well as the decreased thermal stability and substrate affinity of proteinase K upon Ca2+ removal, but also complement the experimentally determined structural and biochemical data.

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.

Pattern of Mutation Rates in the Germline of Drosophila melanogaster Males from a Large-Scale Mutation Screening Experiment

The sperm or eggs of sexual organisms go through a series of cell divisions from the fertilized egg; mutations can occur at each division. Mutations in the lineage of cells leading to the sperm or eggs are of particular importance because many such mutations may be shared by somatic tissues and also may be inherited, thus having a lasting consequence. For decades, little has been known about the pattern of the mutation rates along the germline development. Recently it was shown from a small portion of data that resulted from a large-scale mutation screening experiment that the rates of recessive lethal or nearly lethal mutations differ dramatically during the germline development of Drosophila melanogaster males. In this paper the full data set from the experiment and its analysis are reported by taking advantage of a recent methodologic advance. By analyzing the mutation patterns with different levels of recessive lethality, earlier published conclusions based on partial data are found to remain valid. Furthermore, it is found that for most nearly lethal mutations, the mutation rate at the first cell division is even greater than previous thought compared with those at other divisions. There is also some evidence that the mutation rate at the second division decreases rapidly but is still appreciably greater than those for the rest of the cleavage stage. The mutation rate at spermatogenesis is greater than late cleavage and stem-cell stages, but there is no evidence that rates are different among the five cell divisions of the spermatogenesis. We also found that a modestly biased sampling, leading to slightly more primordial germ cells after the eighth division than those reported in the literature, provides the best fit to the data. These findings provide conceptual and numerical basis for exploring the consequences of differential mutation rates during individual development.

Comparative thermal unfolding study of psychrophilic and mesophilic subtilisin-like serine proteases by molecular dynamics simulations

Molecular dynamics (MD) simulations of a subtilisin-like serine protease VPR from the psychrophilic marine bacterium Vibrio sp. PA-44 and its mesophilic homologue, proteinase K (PRK), have been performed for 20 ns at four different temperatures (300, 373, 473, and 573 K). The comparative analyses of MD trajectories reveal that at almost all temperatures, VPR exhibits greater structural fluctuations/deviations, more unstable regular secondary structural elements, and higher global flexibility than PRK. Although these two proteases follow similar unfolding pathways at high temperatures, VPR initiates unfolding at a lower temperature and unfolds faster at the same high temperatures than PRK. These observations collectively indicate that VPR is less stable and more heat-labile than PRK. Analyses of the structural/geometrical properties reveal that, when compared to PRK, VPR has larger radius of gyration (Rg), less intramolecular contacts and hydrogen bonds (HBs), more protein-solvent HBs, and smaller burial of nonpolar area and larger exposure of polar area. These suggest that the increased flexibility of VPR would be most likely caused by its reduced intramolecular interactions and more favourable protein-solvent interactions arising from the larger exposure of the polar area, whereas the enhanced stability of PRK could be ascribed to its increased intramolecular interactions arising from the better optimized hydrophobicity. The factors responsible for the significant differences in local flexibility between these two proteases were also analyzed and ascertained. This study provides insights into molecular basis of thermostability of homologous serine proteases adapted to different temperatures.

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.