|KOSMOS: a universal morph server for nucleic acids, proteins and their complexes.|
|Jump to Full Text|
|PMID: 22669912 Owner: NLM Status: MEDLINE|
|KOSMOS is the first online morph server to be able to address the structural dynamics of DNA/RNA, proteins and even their complexes, such as ribosomes. The key functions of KOSMOS are the harmonic and anharmonic analyses of macromolecules. In the harmonic analysis, normal mode analysis (NMA) based on an elastic network model (ENM) is performed, yielding vibrational modes and B-factor calculations, which provide insight into the potential biological functions of macromolecules based on their structural features. Anharmonic analysis involving elastic network interpolation (ENI) is used to generate plausible transition pathways between two given conformations by optimizing a topology-oriented cost function that guarantees a smooth transition without steric clashes. The quality of the computed pathways is evaluated based on their various facets, including topology, energy cost and compatibility with the NMA results. There are also two unique features of KOSMOS that distinguish it from other morph servers: (i) the versatility in the coarse-graining methods and (ii) the various connection rules in the ENM. The models enable us to analyze macromolecular dynamics with the maximum degrees of freedom by combining a variety of ENMs from full-atom to coarse-grained, backbone and hybrid models with one connection rule, such as distance-cutoff, number-cutoff or chemical-cutoff. KOSMOS is available at http://bioengineering.skku.ac.kr/kosmos.|
|Sangjae Seo; Moon Ki Kim|
|Type: Journal Article; Research Support, Non-U.S. Gov't Date: 2012-06-04|
|Title: Nucleic acids research Volume: 40 ISSN: 1362-4962 ISO Abbreviation: Nucleic Acids Res. Publication Date: 2012 Jul|
|Created Date: 2012-06-28 Completed Date: 2012-12-12 Revised Date: 2013-04-01|
Medline Journal Info:
|Nlm Unique ID: 0411011 Medline TA: Nucleic Acids Res Country: England|
|Languages: eng Pagination: W531-6 Citation Subset: IM|
|Department of Nano Science and Technology and School of Mechanical Engineering, Sungkyunkwan University, 300, Cheoncheon-dong, Jangan-gu, Suwon 440-746, Korea.|
|APA/MLA Format Download EndNote Download BibTex|
Lactoferrin / chemistry
Multiprotein Complexes / chemistry
Nucleic Acid Conformation
Proteins / chemistry*
RNA / chemistry*
Ribosomes / chemistry
|0/Multiprotein Complexes; 0/Proteins; 63231-63-0/RNA; 9007-49-2/DNA; EC 3.4.21.-/Lactoferrin|
Journal ID (nlm-ta): Nucleic Acids Res
Journal ID (iso-abbrev): Nucleic Acids Res
Journal ID (publisher-id): nar
Journal ID (hwp): nar
Publisher: Oxford University Press
© The Author(s) 2012. Published by Oxford University Press.
Received Day: 26 Month: 2 Year: 2012
Revision Received Day: 4 Month: 5 Year: 2012
Accepted Day: 11 Month: 5 Year: 2012
collection publication date: Month: 7 Year: 2012
Print publication date: Month: 7 Year: 2012
Electronic publication date: Day: 4 Month: 6 Year: 2012
pmc-release publication date: Day: 4 Month: 6 Year: 2012
Volume: 40 Issue: Web Server issue
First Page: W531 Last Page: W536
PubMed Id: 22669912
Publisher Id: gks525
|KOSMOS: a universal morph server for nucleic acids, proteins and their complexes|
|Moon Ki Kim12*|
1Department of Nano Science and Technology and 2School of Mechanical Engineering, Sungkyunkwan University, 300, Cheoncheon-dong, Jangan-gu, Suwon 440-746, Korea
|Correspondence: *To whom correspondence should be addressed. Tel: +82 31 299 4840; Fax: +82 31 290 5889; Email: firstname.lastname@example.org
A great deal of effort has been devoted to studying macromolecular dynamics to elucidate the relationship between structure and function. An earlier attempt with all-atom normal mode analysis (NMA) using empirical potential showed that only a few of the lowest normal modes of a protein are needed to adequately capture the protein’s collective motions, which are believed to be strongly correlated with that protein’s biological functions (1,2). Although this all-atom-based approach is a rigorous method for analyzing molecular vibrations, its computational cost is high for large proteins. Alternatively, the coarse-grained elastic network model (ENM), which employs a single parameter force field, was introduced as a well-balanced approach for understanding the intrinsic vibration features of proteins, with consideration given to the modeling accuracy and computation time (3–7). In particular, this simplified NMA has successfully revealed that large-amplitude domain motions play critical roles in ligand binding, signal transduction and allosteric responses in proteins (8–11).
There are already several web services that provide an easy-to-use interface for the calculation of the normal modes of proteins. ANM (12), oGNM (13) and WEBnm@ (14) similarly adopt point-mass-based models, such as a Cα coarse-grained model and an all-atom model, whereas NOMAD-Ref (15) and ElNémo (16) utilize the building block approach to reduce the computational burden by grouping several consecutive amino acids along the backbone of proteins (17). Obviously, these coarse-graining methods are computationally efficient and easy to implement, but such an automatic approach without any consideration of the intrinsic rigidity (and/or flexibility) of proteins is not suitable for capturing large-scale domain motions effectively (18,19). To overcome this drawback, rigid-cluster-based ENMs have been proposed. Domain-NMA assigns a much higher spring constant to the intra-domain contacts (18). Cluster-NMA (19) and a hybrid elastic network model (HENM) (20) fully utilize mathematical details on the translation and rotation motions of rigid domains.
In addition to the coarse-graining levels, the network connectivity is also one of the most important key factors in determining the quality of ENM-based analysis. Currently, the distance-cutoff method is widely used for the ENM because of its simplicity. However, a small cutoff value reduced the stability of the ENM, whereas a large cutoff value makes the system too stiff so that we cannot observe functional hinge motions, especially from closed to open conformations, in some cases. Alternatively, the chemical-cutoff method has been proposed (21). In this method, different spring constants are assigned to each pair of representative atoms, depending on the types of chemical interactions involved.
KOSMOS (Korea Skku Morph Server) was recently launched, emphasizing its versatility regarding target systems, coarse-graining levels and connection rules. First, KOSMOS accepts any macromolecular structure, including DNA/RNA, proteins and even their complexes. Second, it provides a variety of choices for ENMs, ranging from all-atom to a hybrid model. Information of the rigidity of the macromolecules can be incorporated through the HENM technique. Third, a connection rule can be selected for the ENM, from simple distance-cutoff rule to the more sophisticated chemical-cutoff rule. In addition, KOSMOS can also generate transition pathways between two known conformations using elastic network interpolation (ENI) (22).
Figure 1 demonstrates the entire procedure of KOSMOS. The server is implemented using HTML, Java Script and PHP for the input and output web environments. KOSMOS accepts only standard PDB files. Users can either upload their own file in PDB format or type a PDB ID. When a PDB ID is entered, to prevent duplication, KOSMOS automatically searches the database and then lists all previous results. There are two different query options in KOSMOS, basic and advanced. A basic query requires only a PDB ID as an input. The calculations will be performed automatically using the default parameters and options. An advanced query consists of several consecutive steps to determine more detailed analysis conditions. The user is asked to choose the coarse-graining level, the connection rule, the cutoff value, etc. When a query is completed, a new window is automatically opened, in which the user can check the status of the computation and access to the results page directly.
KOSMOS runs in the MATLAB environment. When the inquiry is accepted, MATLAB starts to run the preprogrammed codes and passes the results to a MySQL database. The user can not only sort the database by individual parameters, but also search the entire database based on various categories provided. The result page provides a 3D visualization tool using a JMOL plug-in so that one can interactively animate the NMA or ENI results. Analysis conditions and downloadable results files are also available on this page.
As a default, KOSMOS performs the macromolecular dynamics analysis based on coarse-grained ENMs. Only the atomic coordinates of Cαs and phosphates are sampled from an input file to represent amino acids and nucleotides, respectively. For ENI computations, KOSMOS automatically aligns the two given PDB files by only taking common residues available in both structures. All of the representative atoms are then networked with virtual linear springs. Depending on the user’s choice of coarse-graining levels and connection rules, the connectivity and spring constants should be determined separately. For example, if the basic query is selected with the distance-cutoff method as a default, a uniform spring constant is assigned to each pair of representative atoms within the cutoff distance of 12 or 24 Å for proteins and DNA/RNA, respectively. These default values are carefully determined by considering the system stability of the ENMs (21). In contrast to the distance-cutoff method, the chemical-cutoff method assigns different spring constants that depend on the type of chemical interaction (e.g. covalent bonds, disulfide bonds, hydrogen bonds, salt bridges and Van der Waals interactions). Once the construction of the ENM is complete, KOSMOS develops a simple Hookean potential energy function:. Here, we assume that δi is small fluctuation of xi(0). Because the standard NMA procedure has been widely reported on in the literature, we will not describe it in detail here. A full description can be obtained from Ref. (22).
ENI is a topology-oriented optimization method to generate transition pathways between two known conformations of the same macromolecule (23,24). ENI starts from standardized coordinates by adjusting the center of mass and principal axis between two conformations because these conformations may have different initial positions and orientations. One can iteratively generate intermediate conformations by solving the following potential-like cost function, in which the initial conformation is systematically deformed toward the target conformation to reduce the root mean square deviation (RMSD) between these conformations:
Unlike the point mass-based models, hybrid models are necessary to decompose the structure into rigid clusters and point masses. One can use various web servers, such as FlexWeb, KINAR-Web, RigidFinder and PiSQRD, to determine rigid clusters in proteins (25–28). Once the rigid cluster information is loaded into KOSMOS, the position vector of residue i within rigid cluster b at time t is defined as:is a translation vector of rigid cluster b and is a rotation matrix formed by:
Lactoferrin is a transferrin protein that is responsible for iron ion transport. Ligand binding induces a conformational change in lactoferrin that involves a rigid-body rotation of the N2 domain relative to the N1 domain (29). To observe such a ‘tweezers-like’ motion, we performed the ENI from the holo form (PDB ID: 1LFG) to the Apo form (PDB ID: 1LFH). We evaluated intermediate conformations with cost and RMSD values to verify the plausibility of the computed pathway, which primarily consists of vertical bending and tweezers-like motion (Figure 2A). Both the distance-cutoff NMA and the chemical-cutoff NMA of the closed forms of lactoferrin were tested. As shown in Figure 2B, the conventional distance-cutoff method causes over-constraints between two close lobes, whereas the chemical-cutoff method properly rules out such inter-domain interactions. As a result, Figure 2C demonstrates that the chemical-cutoff NMA captures both the bending motion and the tweezers-like motion at Mode 1 and Mode 4, respectively. The distance-cutoff NMA shows only a bending motion in the low-frequency modes. This result shows the superior performance of the chemical-cutoff method over the simpler distance-cutoff method in terms of the reality of the ENM, especially for the closed forms of macromolecules.
The ribosome is responsible for the protein synthesis process, called translation, which decodes the genetic information of mRNA to synthesize proteins. The 70S ribosome consists of a small 30S subunit containing the 16S rRNA and 21 proteins and a large 50S subunit composed of the 23S rRNA, the 5S rRNA and 34 proteins. Frank and Agrawal (30) have experimentally observed that elongation factor G induces a ratchet-like rotation in the 30S subunit relative to the 50S subunit during ribosomal translocation. There have been several computational efforts to determine the large conformational transition of the ribosome in atomic detail. Specifically, Wang and Tama have individually reproduced the ratchet-like motion by NMA based on a coarse-grained representation of the ribosome (31,32). The 70S ribosome structure was built by fitting the 30S and 50S subunits (PDB ID: 1GIX and 1GIY) and was tested on various NMA web servers. Most servers failed to analyze this protein–RNA complex. Although WEBnm@ and ElNémo returned NMA results, WEBnm@ missed all of the RNA information during computation and ElNémo showed a limitation in capturing the collective motion evidenced by small collectivity values <0.02. Figure 3A illustrates the HENM of the ribosome constructed by a mixture of five rigid clusters and coarse-grained point masses representing tRNA, rRNA and proteins near the flexible interface between the two subunits. As shown in Figure 3B, the predicted motions on KOSMOS demonstrated good agreement with the current hypothesis (30). In addition, unlike the previous NMA results, KOSMOS returned the main ratchet-like motion as the lowest mode of the ribosome, which was found at the third lowest mode elsewhere (31,32). This result implies that the proposed HENM can capture the intrinsic dynamics of the ribosome more precisely than other modeling tools.
A fully automated ENM-based morph server, called KOSMOS, was recently launched to provide both harmonic and anharmonic analyses represented by NMA and ENI, respectively. Although several webs services already exist that deal with either NMA or transition pathway generation in proteins, KOSMOS is the first automated web server handling analyses for both DNA/RNA and proteins, with an extensive choice of various coarse-graining types and connection rules. In particular, the HENM is a computationally efficient method for studying supramolecular dynamics. The chemical-cutoff method provides more realistic NMA results for macromolecules in closed forms. Regarding visualization, KOSMOS adopts a user-oriented JMol plug-in, which gives us another degree of freedom in 3D interactive visualization. Consequently, the versatility and reliability of KOSMOS will have a significant impact on the area of computational structural biology and bioinformatics by making open access to analysis results, which is optional in KOSMOS, and source codes available to the public.
National Research Foundation of Korea funded by the Ministry of Education, Science and Technology [R33-10079, 2011-0014584]; Faculty Research Funds from Sungkyunkwan University. Funding for open access charge: National Research Foundation of Korea funded by the Ministry of Education, Science and Technology [R33-10079].
Conflict of interest statement. None declared.
|1.||Go N,Noguti T,Nishikawa T. Dynamics of a small globular protein in terms of low-frequency vibrational modesProc. Natl Acad. Sci. USAYear: 198380369637006574507|
|2.||Brooks B,Karplus M. Normal modes for specific motions of macromolecules: application to the hinge-bending mode of lysozymeProc. Natl Acad. Sci. USAYear: 198582499549993860838|
|3.||Tirion MM. Large amplitude elastic motions in proteins from a single-parameter, atomic analysisPhys. Rev. Lett.Year: 1996771905190810063201|
|4.||Bahar I,Atilgan AR,Erman B. Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potentialFold. Des.Year: 199721731819218955|
|5.||Atilgan AR,Durell SR,Jernigan RL,Demirel MC,Keskin O,Bahar I. Anisotropy of fluctuation dynamics of proteins with an elastic network modelBiophys. J.Year: 20018050551511159421|
|6.||Bahar I,Rader AJ. Coarse-grained normal mode analysis in structural biologyCurr. Opin. Struct. Biol.Year: 20051558659216143512|
|7.||Ma J. Usefulness and limitations of normal mode analysis in modeling dynamics of biomolecular complexesStructureYear: 20051337338015766538|
|8.||Thomas A,Field MJ,Perahia D. Analysis of the low-frequency normal modes of the R state of aspartate transcarbamylase and a comparison with the T state modesJ. Mol. Biol.Year: 19962614905068780788|
|9.||Ma J,Karplus M. The allosteric mechanism of the chaperonin GroEL: a dynamic analysisProc. Natl Acad. Sci. USAYear: 199895850285079671707|
|10.||Li G,Cui Q. A coarse-grained normal mode approach for macromolecules: an efficient implementation and application to Ca(2+)-ATPaseBiophys. J.Year: 2002832457247412414680|
|11.||Xu C,Tobi D,Bahar I. Allosteric changes in protein structure computed by a simple mechanical model: hemoglobin T<–>R2 transitionJ. Mol. Biol.Year: 200333315316814516750|
|12.||Eyal E,Yang LW,Bahar I. Anisotropic network model: systematic evaluation and a new web interfaceBioinformaticsYear: 2006222619262716928735|
|13.||Yang LW,Rader AJ,Liu X,Jursa CJ,Chen SC,Karimi HA,Bahar I. oGNM: online computation of structural dynamics using the Gaussian Network ModelNucleic Acids Res.Year: 200634W24W3116845002|
|14.||Hollup SM,Salensminde G,Reuter N. WEBnm@: a web application for normal mode analyses of proteinsBMC Bioinform.Year: 2005652|
|15.||Lindahl E,Azuara C,Koehl P,Delarue M. NOMAD-Ref: visualization, deformation and refinement of macromolecular structures based on all-atom normal mode analysisNucleic Acids Res.Year: 200634W52W5616845062|
|16.||Suhre K,Sanejouand YH. ElNemo: a normal mode web server for protein movement analysis and the generation of templates for molecular replacementNucleic Acids Res.Year: 200432W610W61415215461|
|17.||Tama F,Gadea FX,Marques O,Sanejouand YH. Building-block approach for determining low-frequency normal modes of macromoleculesProteinsYear: 2000411710944387|
|18.||Song G,Jernigan RL. An enhanced elastic network model to represent the motions of domain-swapped proteinsProteinsYear: 20066319720916447281|
|19.||Schuyler AD,Chirikjian GS. Normal mode analysis of proteins: a comparison of rigid cluster modes with C(alpha) coarse grainingJ. Mol. Graph. Model.Year: 20042218319314629977|
|20.||Kim MK,Jernigan RL,Chirikjian GS. Rigid-cluster models of conformational transitions in macromolecular machines and assembliesBiophys. J.Year: 200589435515833998|
|21.||Jeong JI,Jang Y,Kim MK. A connection rule for alpha-carbon coarse-grained elastic network models using chemical bond informationJ. Mol. Graph. Model.Year: 20062429630616289973|
|22.||Kim MK,Chirikjian GS,Jernigan RL. Elastic models of conformational transitions in macromoleculesJ. Mol. Graph. Model.Year: 20022115116012398345|
|23.||Kim MK,Jernigan RL,Chirikjian GS. Efficient generation of feasible pathways for protein conformational transitionsBiophys. J.Year: 2002831620163012202386|
|24.||Kim MK,Jang Y,Jeong JI. Using harmonic analysis and optimization to study macromolecular dynamicsInt. J. Control Autom. Syst.Year: 20064382393|
|25.||Jacobs DJ,Rader AJ,Kuhn LA,Thorpe MF. Protein flexibility predictions using graph theoryProteinsYear: 20014415016511391777|
|26.||Aleksiev T,Potestio R,Pontiggia F,Cozzini S,Micheletti C. PiSQRD: a web server for decomposing proteins into quasi-rigid dynamical domainsBioinformaticsYear: 2009252743274419696046|
|27.||Abyzov A,Bjornson R,Felipe M,Gerstein M. RigidFinder: a fast and sensitive method to detect rigid blocks in large macromolecular complexesProteinsYear: 20107830932419705487|
|28.||Fox N,Jagodzinski F,Li Y,Streinu I. KINARI-Web: a server for protein rigidity analysisNucleic Acids Res.Year: 201139W177W18321693559|
|29.||Gerstein M,Anderson BF,Norris GE,Baker EN,Lesk AM,Chothia C. Domain closure in lactoferrin. Two hinges produce a see-saw motion between alternative close-packed interfacesJ. Mol. Biol.Year: 19932343573728230220|
|30.||Frank J,Agrawal RK. A ratchet-like inter-subunit reorganization of the ribosome during translocationNatureYear: 200040631832210917535|
|31.||Tama F,Valle M,Frank J,Brooks CL. Dynamic reorganization of the functionally active ribosome explored by normal mode analysis and cryo-electron microscopyProc. Natl Acad. Sci. USAYear: 20031009319932312878726|
|32.||Wang Y,Rader AJ,Bahar I,Jernigan RL. Global ribosome motions revealed with elastic network modelJ. Struct. Biol.Year: 200414730231415450299|
Previous Document: Quantitative proteomics profiling of the poly(ADP-ribose)-related response to genotoxic stress.
Next Document: CMWeb: an interactive on-line tool for analysing residue-residue contacts and contact prediction met...