Yifei Kong

How to describe protein motion without amino acid sequence and atomic coordinates

This paper reports a computational method, the quantized elastic deformational model, that can reliably describe the conformational flexibility of a protein in the absence of the amino acid sequence and atomic coordinates. The essence of this method lies in the fact that, in modeling the functionally important conformational changes such as domain movements, it is possible to abandon the traditional concepts of protein structure (bonds, angles, dihedrals, etc.) and treat the protein as an elastic object. The shape and mass distribution of the object are described by the electron density maps, at various resolutions, from methods such as x-ray diffraction or cryo-electron microscopy. The amplitudes and directionality of the elastic deformational modes of a protein, whose patterns match the biologically relevant conformational changes, can then be derived solely based on the electron density map. The method yields an accurate description of protein dynamics over a wide range of resolutions even as low as 15–20 Å at which there is nearly no visually distinguishable internal structures. Therefore, this method dramatically enhances the capability of studying protein motions in structural biology. It is also expected to have ample applications in related fields such as bioinformatics, structural genomics, and proteomics, in which ones ability to extract functional information from the not-so-well-defined structural models is vitally important.

Signaling pathways of PDZ2 domain: A molecular dynamics interaction correlation analysis

PDZ domains are found in many signaling proteins. One of their functions is to provide scaffolds for forming membrane‐associated protein complexes by binding to the carboxyl termini of their partners. PDZ domains are thought also to play a signal transduction role by propagating the information that binding has occurred to remote sites. In this study, a molecular dynamics (MD) simulation‐based approach, referred to as an interaction correlation analysis, is applied to the PDZ2 domain to identify the possible signal transduction pathways. A residue correlation matrix is constructed from the interaction energy correlations between all residue pairs obtained from the MD simulations. Two continuous interaction pathways, starting at the ligand binding pocket, are identified by a hierarchical clustering analysis of the residue correlation matrix. One pathway is mainly localized at the N‐terminal side of helix α1 and the adjacent C‐terminus of loop β1‐β2. The other pathway is perpendicular to the central β‐sheet and extends toward the side of PDZ2 domain opposite to the ligand binding pocket. The results complement previous studies based on multiple sequence analysis, NMR, and MD simulations. Importantly, they reveal the energetic origin of the long‐range coupling. The PDZ2 results, as well as the earlier rhodopsin analysis, show that the interaction correlation analysis is a robust approach for determining pathways of intramolecular signal transduction. Proteins 2009.

Conformational pathways in the gating of Escherichia coli mechanosensitive channel

The pathway of the gating conformational transition of Escherichia coli mechanosensitive channel was simulated, using the recently modeled open and closed structures, by targeted molecular dynamics method. The transition can be roughly viewed as a four-stage process. The initial motion under a lower tension load is predominantly elastic deformation. The opening of the inner hydrophobic pore on a higher tension load takes place after the major expansion of the outer channel dimension. The hypothetical N-terminal S1 helical bundle has been confirmed to form the hydrophobic gate, together with the M1 helices. The sequential breaking of the tandem hydrophobic constrictions on the M1 and S1 helices makes the two parts of the gate strictly coupled, acting as a single gate. The simulation also revealed that there is no significant energetic coupling between the inner S1 bundle and the outer M2 transmembrane helices. The molten-globular-like structural features of the S1 bundle in its intermediate open states may account for the observed multiple subconductance states. Moreover, the intermediate open states of mechanosensitive channels are not symmetric, i.e., the opening does not follow iris-like motion, which sharply contrasts to the potassium channel KcsA.

A Structural-informatics Approach for Mining β-Sheets: Locating Sheets in Intermediate-resolution Density Maps

Here, we report a new computational method, called sheetminer, for mining beta-sheets in the density maps at intermediate resolutions of 6 to 10A. The method employs a multi-step ad hoc morphological analysis of density maps to identify the unique characteristics of beta-sheets. It was tested on density maps from 12 protein crystal structures that were artificially blurred to intermediate resolutions. There are a total of 35 independent beta-sheets with a wide distribution of morphology. The method successfully located 34 of them and missed only one. The method was also applied to an experimental 9A electron cryomicroscopic structure and an 8A X-ray density map. In both cases, the sheet-searching results were found to agree very well with known high-resolution crystal structures. Collectively, these results demonstrate clearly the robustness of sheetminer in locating the regions belonging to beta-sheets in the intermediate-resolution density maps. Furthermore, sheetminer is completely complementary to all other existing computational methods, including helixhunter and threading algorithms. Their combined usage has the potential to significantly enhance the computational modeling capacity for a much more complete interpretation of structural data at intermediate resolutions, from which extraction of functional information would be more effective. This is particularly important in the field of structural genomics, in which the fast screening approach may not always yield crystals that diffract to atomic resolution. An exciting future application of sheetminer is as a valuable tool for revealing the structures of amyloid fibrils that are rich in beta-motifs.

Domain movements in human fatty acid synthase by quantized elastic deformational model

This paper reports the results of applying a computational method called the quantized elastic deformational model, to the determination of conformational flexibility of the supermolecular complex of human fatty acid synthase. The essence of this method is the ability to model large-scale conformational changes such as domain movements by treating the protein as an elastic object without the knowledge of protein primary sequence and atomic coordinates. The calculation was based on the electron density maps of the synthase at 19 Å. The results suggest that the synthase is a very flexible molecule. Two types of flexible hinges in the structure were identified. One is an intersubunit hinge formed by the intersubunit connection and the other is an intrasubunit hinge located between domains I and II. Despite the fact that the dimeric synthase has a chemically symmetric structure, large domain movements around the hinge region occur in various directions and allow the molecule to adopt a wide range of conformations. These domain movements are likely to be important in facilitating and regulating the entire palmitate synthesis by coordinating the communication between components of the molecule, for instance, adjusting the distance between various active sites inside the catalytic reaction center. Finally, the ability to describe protein motions of a supermolecular complex, without the information of protein sequence and atomic coordinates, is a major advance in computational modeling of protein dynamics. The method provides an unprecedented ability to model protein motions at such a low resolution of structure.

Intrinsic flexibility and gating mechanism of the potassium channel KcsA

The gating mechanism of the potassium channel KcsA was studied by normal mode analysis. The results provided an atomic description of the locations of the pivot points and the motional features of key structural elements in the gating process. Two pivot points were found in the motions of the inner TM2 helical bundle that directly modulate the size of the central channel pore. One point is an intrasubunit hinge point that sharply divides the structural flexibility between the more rigid selectivity filter and the more mobile peripheral transmembrane helices. Such a division is vital for KcsA because it permits the large-scale motions of transmembrane helices required for the gating and, in the meantime, maintains the rigidity of the filter region essential for the selectivity. The other pivot point is an intersubunit one at which all four TM2 helices are bundled together. During the gating process, each TM2 helix undergoes a lever-like swinging motion pivoting on the intrasubunit hinge, and the entire TM2 bundle undergoes a concerted rotational motion around the central channel axis constrained around the intersubunit bundle point. This series of motions leads to a dramatic enlargement of the intracellular gate without loosening up the structural integrity.

Substructure synthesis method for simulating large molecular complexes

This paper reports a computational method for describing the conformational flexibility of very large biomolecular complexes using a reduced number of degrees of freedom. It is called the substructure synthesis method, and the basic concept is to treat the motions of a given structure as a collection of those of an assemblage of substructures. The choice of substructures is arbitrary and sometimes quite natural, such as domains, subunits, or even large segments of biomolecular complexes. To start, a group of low-frequency substructure modes is determined, for instance by normal mode analysis, to represent the motions of the substructure. Next, a desired number of substructures are joined together by a set of constraints to enforce geometric compatibility at the interface of adjacent substructures, and the modes for the assembled structure can then be synthesized from the substructure modes by applying the Rayleigh–Ritz principle. Such a procedure is computationally much more desirable than solving the full eigenvalue problem for the whole assembled structure. Furthermore, to show the applicability to biomolecular complexes, the method is used to study F-actin, a large filamentous molecular complex involved in many cellular functions. The results demonstrate that the method is capable of studying the motions of very large molecular complexes that are otherwise completely beyond the reach of any conventional methods.

Conformational Flexibility of Pyruvate Dehydrogenase Complexes: A Computational Analysis by Quantized Elastic Deformational Model

Pyruvate dehydrogenase complex (PDC) is one of the largest multienzyme complexes known and consists of a dodecahedral E2 core to which other components are attached. We report the results of applying a new computational method, quantized elastic deformational model, to simulating the conformational fluctuations of the truncated E2 core, using low-resolution electron cryomicroscopy density maps. The motional features are well reproduced; especially, the symmetric breathing mode revealed in simulation is nearly identical with what was observed experimentally. Structural details of the motions of the trimeric building blocks, which are critical to facilitating the global expansion and contraction of the complex, were revealed. Using the low-resolution maps from electron cryomicroscopy reconstructions, the simulations showed a picture of the motional mechanism of the PDC core, which is an example without precedent of thermally activated global dynamics. Moreover, the current results support an earlier suggestion that, at low resolution and without the use of amino acid sequence and atomic coordinates, it is possible for computer simulations to provide an accurate description of protein dynamics.

Simulation of F-Actin Filaments of Several Microns

Here we report the results of applying substructure synthesis method to the simulation of F-actin filaments of several microns in length. The elastic deformational modes of long F-actin filaments were generated from the vibrational modes of the 13-subunit repeat of F-actin using a hierarchical synthesis scheme. The computationally synthesized deformational modes, in the very low-frequency regime, are in good agreement with theoretical solutions for long homogeneous elastic rods, which confirmed the usefulness of substructure synthesis method. Other low-frequency modes carry rich local deformational features that are unique to F-actins. All these modes thus provide a theoretical basis set for a description of spontaneously occurring thermal deformations, such as undulations, of the filaments. The results demonstrate that substructure synthesis method, as a method for computational modal analysis, is capable of scaling up the microscopic dynamic information, obtained from atomistic simulations, to a wide range of macroscopic length scale. Moreover, the combination of substructure synthesis method and hierarchical synthesis scheme provides an effective way in dealing with complex systems of periodic repeats that are abundant in cells.