WEBnm@: a web application for normal mode analyses of proteins

Normal mode analysis based on an elastic network model for biomolecules in the Protein Data Bank, which uses dihedral angles as independent variables

We have developed a computer program, named PDBETA, that performs normal mode analysis (NMA) based on an elastic network model that uses dihedral angles as independent variables. Taking advantage of the relatively small number of degrees of freedom required to describe a molecular structure in dihedral angle space and a simple potential-energy function independent of atom types, we aimed to develop a program applicable to a full-atom system of any molecule in the Protein Data Bank (PDB). The algorithm for NMA used in PDBETA is the same as the computer program FEDER/2, developed previously. Therefore, the main challenge in developing PDBETA was to find a method that can automatically convert PDB data into molecular structure information in dihedral angle space. Here, we illustrate the performance of PDBETA with a protein-DNA complex, a protein-tRNA complex, and some non-protein small molecules, and show that the atomic fluctuations calculated by PDBETA reproduce the temperature factor data of these molecules in the PDB. A comparison was also made with elastic-network-model based NMA in a Cartesian-coordinate system.

Localized lipid packing of transmembrane domains impedes integrin clustering

Integrin clustering plays a pivotal role in a host of cell functions. Hetero-dimeric integrin adhesion receptors regulate cell migration, survival, and differentiation by communicating signals bidirectionally across the plasma membrane. Thus far, crystallographic structures of integrin components are solved only separately, and for some integrin types. Also, the sequence of interactions that leads to signal transduction remains ambiguous. Particularly, it remains controversial whether the homo-dimerization of integrin transmembrane domains occurs following the integrin activation (i.e. when integrin ectodomain is stretched out) or if it regulates integrin clustering. This study employs molecular dynamics modeling approaches to address these questions in molecular details and sheds light on the crucial effect of the plasma membrane. Conducting a normal mode analysis of the intact αllbβ3 integrin, it is demonstrated that the ectodomain and transmembrane-cytoplasmic domains are connected via a membrane-proximal hinge region, thus merely transmembrane-cytoplasmic domains are modeled. By measuring the free energy change and force required to form integrin homo-oligomers, this study suggests that the β-subunit homo-oligomerization potentially regulates integrin clustering, as opposed to α-subunit, which appears to be a poor regulator for the clustering process. If α-subunits are to regulate the clustering they should overcome a high-energy barrier formed by a stable lipid pack around them. Finally, an outside-in activation-clustering scenario is speculated, explaining how further loading the already-active integrin affects its homo-oligomerization so that focal adhesions grow in size.

Focal adhesions are complex, dynamic structures of multiple proteins that act as the cell's mechanical anchorage to its surrounding. Integrins are proteins linking the cell inner and outer environments, which act as a bridge that crosses the cell membrane. Integrins respond to mechanical loads exerted to them by changing their conformations. Several diseases, such as atherosclerosis and different types of cancer, are caused by altered function of integrins. Essential to the formation of focal adhesions is the process of integrin clustering. Bidirectional integrin signaling involves conformational changes in this protein, clustering, and finally the assembly of a large intracellular adhesion complex. Integrin clustering is defined as the interaction of integrins to form lateral assemblies that eventually lead to focal adhesion formation. The effect of the plasma membrane on formation of integrin clusters has been largely neglected in current literature; subsequently some apparently contradictory data has been reported by a number of researchers in the field. Using a molecular dynamics modeling approach, a computational method that simulates systems in a full-atomic scale, we probe the role of the plasma membrane in integrin clustering and hypothesize a clustering scenario that explains the relationship between integrin activation and focal adhesion growth.

Molecular modelling and simulations in cancer research

The complexity of cancer and the vast amount of experimental data available have made computer-aided approaches necessary. Biomolecular modelling techniques are becoming increasingly easier to use, whereas hardware and software are becoming better and cheaper. Cross-talk between theoretical and experimental scientists dealing with cancer-research from a molecular approach, however, is still uncommon. This is in contrast to other fields, such as amyloid-related diseases, where molecular modelling studies are widely acknowledged. The aim of this review paper is therefore to expose some of the more common approaches in molecular modelling to cancer scientists in simple terms, illustrating success stories while also revealing the limitations of computational studies at the molecular level.

Structure-based molecular modeling approaches to GPCR oligomerization

Classical structure-based drug design techniques using G-protein-coupled receptors (GPCRs) as targets focus nearly exclusively on binding at the orthosteric site of a single receptor. Dimerization and oligomerization of GPCRs, proposed almost 30 years ago, have, however, crucial relevance for drug design. Targeting these complexes selectively or designing small molecules that affect receptor-receptor interactions might provide new opportunities for novel drug discovery. In order to study the mechanisms and dynamics that rule GPCRs oligomerization, it is essential to understand the dynamic process of receptor-receptor association and to identify regions that are suitable for selective drug binding, which may be determined with experimental methods such as Förster resonance energy transfer (FRET) or Bioluminescence resonance energy transfer (BRET) and computational sequence- and structure-based approaches. The aim of this chapter is to provide a comprehensive description of the structure-based molecular modeling methods for studying GPCR dimerization, that is, protein-protein docking, molecular dynamics, normal mode analysis, and electrostatics studies.

Computer-Aided Protein Directed Evolution: a Review of Web Servers, Databases and other Computational Tools for Protein Engineering

The combination of computational and directed evolution methods has proven a winning strategy for protein engineering. We refer to this approach as computer-aided protein directed evolution (CAPDE) and the review summarizes the recent developments in this rapidly growing field. We will restrict ourselves to overview the availability, usability and limitations of web servers, databases and other computational tools proposed in the last five years. The goal of this review is to provide concise information about currently available computational resources to assist the design of directed evolution based protein engineering experiment.

Application of Drug-Perturbed Essential Dynamics/Molecular Dynamics (ED/MD) to Virtual Screening and Rational Drug Design

We present here the first application of a new algorithm, essential dynamics/molecular dynamics (ED/MD), to the field of small molecule docking. The method uses a previously existing molecular dynamics (MD) ensemble of a protein or protein-drug complex to generate, with a very small computational cost, perturbed ensembles which represent ligand-induced binding site flexibility in a more accurate way than the original trajectory. The use of these perturbed ensembles in a standard docking program leads to superior performance than the same docking procedure using the crystal structure or ensembles obtained from conventional MD simulations as templates. The simplicity and accuracy of the method opens up the possibility of introducing protein flexibility in high-throughput docking experiments.

KOSMOS: a universal morph server for nucleic acids, proteins and their complexes

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.

Wordom: a user-friendly program for the analysis of molecular structures, trajectories, and free energy surfaces

Wordom is a versatile, user-friendly, and efficient program for manipulation and analysis of molecular structures and dynamics. The following new analysis modules have been added since the publication of the original Wordom paper in 2007: assignment of secondary structure, calculation of solvent accessible surfaces, elastic network model, motion cross correlations, protein structure network, shortest intra-molecular and inter-molecular communication paths, kinetic grouping analysis, and calculation of mincut-based free energy profiles. In addition, an interface with the Python scripting language has been built and the overall performance and user accessibility enhanced. The source code of Wordom (in the C programming language) as well as documentation for usage and further development are available as an open source package under the GNU General Purpose License from http://wordom.sf.net.

Computational methods for studying serpin conformational change and structural plasticity

Currently, over a hundred high-resolution structures of serpins are available, exhibiting a wide range of conformations. However, our understanding of serpin dynamics and conformational change is still limited, mainly due to challenges of monitoring structural changes and characterizing transient conformations using experimental methods. Insight can be provided, however, by employing theoretical and computational approaches. In this chapter, we present an overview of such methods, focusing on molecular dynamics and simulation. As serpin conformational dynamics span a wide range of timescales, we discuss the relative merits of each method and suggest which method is suited to specific conformational phenomena.

Conformational dynamics data bank: a database for conformational dynamics of proteins and supramolecular protein assemblies

The conformational dynamics data bank (CDDB, http://www.cdyn.org) is a database that aims to provide comprehensive results on the conformational dynamics of high molecular weight proteins and protein assemblies. Analysis is performed using a recently introduced coarse-grained computational approach that is applied to the majority of structures present in the electron microscopy data bank (EMDB). Results include equilibrium thermal fluctuations and elastic strain energy distributions that identify rigid versus flexible protein domains generally, as well as those associated with specific functional transitions, and correlations in molecular motions that identify molecular regions that are highly coupled dynamically, with implications for allosteric mechanisms. A practical web-based search interface enables users to easily collect conformational dynamics data in various formats. The data bank is maintained and updated automatically to include conformational dynamics results for new structural entries as they become available in the EMDB. The CDDB complements static structural information to facilitate the investigation and interpretation of the biological function of proteins and protein assemblies essential to cell function.

A molecular dynamics investigation of vinculin activation

Vinculin activation plays a critical role in focal adhesion initiation and formation. In its native state, vinculin is in an autoinhibitory conformation in which domain 1 prevents interaction of the vinculin tail domain with actin by steric hindrance. Once activated, vinculin is able to interact with both actin and talin. Several hypotheses have been put forth addressing the mechanisms of vinculin activation. One set of studies suggests that vinculin interaction with talin is sufficient to cause activation, whereas another set of studies suggests that a simultaneous interaction with several binding partners is necessary to achieve vinculin activation. Using molecular-dynamics (MD) simulations, we investigate the mechanisms of vinculin activation and suggest both a trajectory of conformational changes leading to vinculin activation, and key structural features that are likely involved in stabilizing the autoinhibited conformation. Assuming that the simultaneous interaction of vinculin with both actin and talin causes a stretching force on vinculin, and that vinculin activation results from a removal of steric hindrance blocking the actin-binding sites, we simulate with MD the stretching and activation of vinculin. The MD simulations are further confirmed by normal-mode analysis and simulation after residue modification. Taken together, the results of these simulations suggest that bending of the vinculin-binding-site region in vinculin away from the vinculin tail is the likely trajectory of vinculin activation.

Coarse-grained models to study dynamics of nanoscale biomolecules and their applications to the ribosome

Biopolymers are of dynamic nature and undergo functional motions spanning a large spectrum of timescales. To study the internal dynamics of nano-sized molecular complexes that exceed hundred thousands of atoms with atomic detail is computationally inefficient. Therefore, to achieve both the spatial and temporal scales of biological interest coarse-grained models of macromolecules are often used. By uniting groups of atoms into single interacting centers one decreases the resolution of the system and gets rid of the irrelevant degrees of freedom. This simplification, even though it requires parameterization, makes the studies of biomolecular dynamics computationally tractable and allows us to reach beyond the microsecond time frame. Here, I review the coarse-grained models of macromolecules composed of proteins and nucleic acids. I give examples of one-bead models that were developed to investigate the internal dynamics and focus on their applications to the ribosome--the nanoscale protein synthesis machine.

Integrated prediction of one-dimensional structural features and their relationships with conformational flexibility in helical membrane proteins

Background

Many structural properties such as solvent accessibility, dihedral angles and helix-helix contacts can be assigned to each residue in a membrane protein. Independent studies exist on the analysis and sequence-based prediction of some of these so-called one-dimensional features. However, there is little explanation of why certain residues are predicted in a wrong structural class or with large errors in the absolute values of these features. On the other hand, membrane proteins undergo conformational changes to allow transport as well as ligand binding. These conformational changes often occur via residues that are inherently flexible and hence, predicting fluctuations in residue positions is of great significance.

Results

We performed a statistical analysis of common patterns among selected one-dimensional equilibrium structural features (ESFs) and developed a method for simultaneously predicting all of these features using an integrated system. Our results show that the prediction performance can be improved if multiple structural features are trained in an integrated model, compared to the current practice of developing individual models. In particular, the performance of the solvent accessibility and bend-angle prediction improved in this way. The well-performing bend-angle prediction can be used to predict helical positions with severe kinks at a modest success rate. Further, we showed that single-chain conformational dynamics, measured by B-factors derived from normal mode analysis, could be predicted from observed and predicted ESFs with good accuracy. A web server was developed (http://tardis.nibio.go.jp/netasa/htmone/) for predicting the one-dimensional ESFs from sequence information and analyzing the differences between the predicted and observed values of the ESFs.

Conclusions

The prediction performance of the integrated model is significantly better than that of the models performing the task separately for each feature for the solvent accessibility and bend-angle predictions. The predictability of the features also plays a role in determining flexible positions. Although the dynamics studied here concerns local atomic fluctuations, a similar analysis in terms of global structural features will be helpful in predicting large-scale conformational changes, for which work is in progress.

Multiple interests in structural models of DARC transmembrane protein

Duffy Antigen Receptor for Chemokines (DARC) is an unusual transmembrane chemokine receptor which (i) binds the two main chemokine families and (ii) does not transduct any signal as it lacks the DRY consensus sequence. It is considered as silent chemokine receptor, a tank useful for chemiotactism. DARC had been particularly studied as a major actor of malaria infection by Plasmodium vivax. It is also implicated in multiple chemokine inflammation, inflammatory diseases, in cancer and might play a role in HIV infection and AIDS. In this review, we focus on the interest to build structural model of DARC to understand more precisely its abilities to bind its physiological ligand CXCL8 and its malaria ligand. We also present innovative development on VHHs able to bind DARC protein. We underline difficulties and limitations of such bioinformatics approaches and highlight the crucial importance of biological data to conduct these kinds of researches.

In silico studies on DARC

The Duffy Antigen/Receptor for Chemokine (DARC) is a seven segment transmembrane protein. It was firstly discovered as a blood group antigen and was the first specific gene locus assigned to a specific autosome in man. It became more famous as an erythrocyte receptor for malaria parasites (Plasmodium vivax and Plasmodium knowlesi), and finally for chemokines. DARC is an unorthodox chemokine receptor as (i) it binds chemokines of both CC and CXC classes and (ii) it lacks the Asp-Arg-Tyr consensus motif in its second cytoplasmic loop hence cannot couple to G proteins and activate their signaling pathways. DARC had also been associated to cancer progression, numerous inflammatory diseases, and possibly to AIDS. In this review, we will summarize important biological data on DARC. Then we shall focus on recent development of the elaboration and analyzes of structural models of DARC. We underline the difficulty to propose pertinent structural models of transmembrane protein using comparative modeling process, and other dedicated approaches as the Protein Blocks. The chosen structural models encompass most of the biochemical data known to date. Finally, we present recent development of protein-protein docking between DARC structural models and CXCL-8 structures. We propose a hierarchical search based on separated rigid and flexible docking.

Molecular mechanics of the alpha-actinin rod domain: bending, torsional, and extensional behavior

α-Actinin is an actin crosslinking molecule that can serve as a scaffold and maintain dynamic actin filament networks. As a crosslinker in the stressed cytoskeleton, α-actinin can retain conformation, function, and strength. α-Actinin has an actin binding domain and a calmodulin homology domain separated by a long rod domain. Using molecular dynamics and normal mode analysis, we suggest that the α-actinin rod domain has flexible terminal regions which can twist and extend under mechanical stress, yet has a highly rigid interior region stabilized by aromatic packing within each spectrin repeat, by electrostatic interactions between the spectrin repeats, and by strong salt bridges between its two anti-parallel monomers. By exploring the natural vibrations of the α-actinin rod domain and by conducting bending molecular dynamics simulations we also predict that bending of the rod domain is possible with minimal force. We introduce computational methods for analyzing the torsional strain of molecules using rotating constraints. Molecular dynamics extension of the α-actinin rod is also performed, demonstrating transduction of the unfolding forces across salt bridges to the associated monomer of the α-actinin rod domain.

The cell interacts with its environment in both biochemical and mechanical ways. In this study we explore one of the ways in which the cell interacts mechanically with its environment. α-Actinin is a cytoskeletal crosslinker: it functions to scaffold the cytoskeletal actin filaments that provide mechanical reinforcement to the cell. In its functional environment α-actinin is exposed to a multitude of mechanical stresses as it attaches itself to a dynamic network of actin filaments. The actin filaments extend, rotate, and bend the α-actinin crosslinkers. In this study we employ molecular dynamics techniques to understand the structural characteristics of α-actinin that underlie its ability to provide a scaffold in such a stressed environment. We analyzed the natural frequencies of α-actinin and simulated force-induced bending, extension, and twisting. Our results suggest that α-actinin has structural flexibility facilitating crosslinking in a dynamic environment and also structural rigidity stabilizing the linkage in the stressed environment. We have discovered novel natural bending movements of the rod domain that enhance its function as a crosslinker. We have also demonstrated the specific structural characteristics of α-actinin that give it the previously suggested property of having partial flexibility. Our results enhance the understanding of structural mechanics in the cytoskeletal molecules.

HingeProt: automated prediction of hinges in protein structures

Proteins are highly flexible molecules. Prediction of molecular flexibility aids in the comprehension and prediction of protein function and in providing details of functional mechanisms. The ability to predict the locations, directions, and extent of molecular movements can assist in fitting atomic resolution structures to low-resolution EM density maps and in predicting the complex structures of interacting molecules (docking). There are several types of molecular movements. In this work, we focus on the prediction of hinge movements. Given a single protein structure, the method automatically divides it into the rigid parts and the hinge regions connecting them. The method employs the Elastic Network Model, which is very efficient and was validated against a large data set of proteins. The output can be used in applications such as flexible protein-protein and protein-ligand docking, flexible docking of protein structures into cryo-EM maps, and refinement of low-resolution EM structures. The web server of HingeProt provides convenient visualization of the results and is available with two mirror sites at http://www.prc.boun.edu.tr/appserv/prc/HingeProt3 and http://bioinfo3d.cs.tau.ac.il/HingeProt/.

Dealing with structural variability in molecular replacement and crystallographic refinement through normal-mode analysis

Normal-mode analysis (NMA) can be used to generate multiple structural variants of a given template model, thereby increasing the chance of finding the molecular-replacement solution. Here, it is shown that it is also possible to directly refine the amplitudes of the normal modes against experimental data (X-ray or cryo-EM), generalizing rigid-body refinement methods by adding just a few additional degrees of freedom that sample collective and large-amplitude movements. It is also argued that the situation where several (conformations of) models are present simultaneously in the crystal can be studied with adjustable occupancies using techniques derived from statistical thermodynamics and already used in molecular modelling.

Two-level approach to efficient visualization of protein dynamics

Proteins are highly flexible and large amplitude deformations of their structure, also called slow dynamics, are often decisive to their function. We present a two-level rendering approach that enables visualization of slow dynamics of large protein assemblies. Our approach is aligned with a hierarchical model of large scale molecules. Instead of constantly updating positions of large amounts of atoms, we update the position and rotation of residues, i.e., higher level building blocks of a protein. Residues are represented by one vertex only indicating its position and additional information defining the rotation. The atoms in the residues are generated on-the-fly on the GPU, exploiting the new graphics hardware geometry shader capabilities. Moreover, we represent the atoms by billboards instead of tessellated spheres. Our representation is then significantly faster and pixel precise. We demonstrate the usefulness of our new approach in the context of our collaborative bioinformatics project.

Thorough validation of protein normal mode analysis: a comparative study with essential dynamics

The deformation patterns of a large set of representative proteins determined by essential dynamics extracted from atomistic simulations and coarse-grained normal mode analysis are compared. Our analysis shows that the deformational space obtained with both approaches is quite similar when taking into account a representative number of modes. The results provide not only a comprehensive validation of the use of a low-frequency modal spectrum to describe protein flexibility, but also a complete picture of normal mode limitations.

Determinants for DNA target structure selectivity of the human LINE-1 retrotransposon endonuclease

The human LINE-1 endonuclease (L1-EN) is the targeting endonuclease encoded by the human LINE-1 (L1) retrotransposon. L1-EN guides the genomic integration of new L1 and Alu elements that presently account for ∼28% of the human genome. L1-EN bears considerable technological interest, because its target selectivity may ultimately be engineered to allow the site-specific integration of DNA into defined genomic locations. Based on the crystal structure, we generated L1-EN mutants to analyze and manipulate DNA target site recognition. Crystal structures and their dynamic and functional analysis show entire loop grafts to be feasible, resulting in altered specificity, while individual point mutations do not change the nicking pattern of L1-EN. Structural parameters of the DNA target seem more important for recognition than the nucleotide sequence, and nicking profiles on DNA oligonucleotides in vitro are less well defined than the respective integration site consensus in vivo. This suggests that additional factors other than the DNA nicking specificity of L1-EN contribute to the targeted integration of non-LTR retrotransposons.

TMM@: a web application for the analysis of transmembrane helix mobility

Background

To understand the mechanism by which a protein transmits a signal through the cell membrane, an understanding of the flexibility of its transmembrane (TM) region is essential. Normal Mode Analysis (NMA) has become the method of choice to investigate the slowest motions in macromolecular systems. It has been widely used to study transmembrane channels and pumps. It relies on the hypothesis that the vibrational normal modes having the lowest frequencies (also named soft modes) describe the largest movements in a protein and are the ones that are functionally relevant. In particular NMA can be used to study dynamics of TM regions, but no tool making this approach available for non-experts, has been available so far.

Results

We developed the web-application TMM@ (TransMembrane α-helical Mobility analyzer). It uses NMA to characterize the propensity of transmembrane α-helices to be displaced. Starting from a structure file at the PDB format, the server computes the normal modes of the protein and identifies which helices in the bundle are the most mobile. Each analysis is performed independently from the others and results can be visualized using only a web browser. No additional plug-in or software is required. For users who would like to further analyze the output data with their favourite software, raw results can also be downloaded.

Conclusion

We built a novel and unique tool, TMM@, to study the mobility of transmembrane α-helices. The tool can be applied to for example membrane transporters and provides biologists studying transmembrane proteins with an approach to investigate which α-helices are likely to undergo the largest displacements, and hence which helices are most likely to be involved in the transportation of molecules in and out of the cell.

Crystal structure and allosteric regulation of the cytoplasmic Escherichia coli L-asparaginase I

AnsA is the cytoplasmic asparaginase from Escherichia coli involved in intracellular asparagine utilization. Analytical ultracentifugation and X-ray crystallography reveal that AnsA forms a tetrameric structure as a dimer of two intimate dimers. Kinetic analysis of the enzyme reveals that AnsA is positively cooperative, displaying a sigmoidal substrate dependence curve with an [S](0.5) of 1 mM L-asparagine and a Hill coefficient (n(H)) of 2.6. Binding of L-asparagine to an allosteric site was observed in the crystal structure concomitant with a reorganization of the quarternary structure, relative to the apo enzyme. The carboxyl group of the bound asparagine makes salt bridges and hydrogen bonds to Arg240, while the N(delta2) nitrogen interacts with Thr162. Mutation of Arg240 to Ala increases the [S](0.5) value to 5.9 mM, presumably by reducing the affinity of the site for L-asparagine, although the enzyme retains cooperativity. Mutation of Thr162 to Ala results in an active enzyme with no cooperativity. Transmission of the signal from the allosteric site to the active site appears to involve subtle interactions at the dimer-dimer interface and relocation of Gln118 into the vicinity of the active site to position the probable catalytic water molecule. These data define the structural basis for the cooperative regulation of the intracellular asparaginase that is required for proper functioning within the cell.

Structure of phenylalanine hydroxylase from Colwellia psychrerythraea 34H, a monomeric cold active enzyme with local flexibility around the active site and high overall stability

The characteristic of cold-adapted enzymes, high catalytic efficiency at low temperatures, is often associated with low thermostability and high flexibility. In this context, we analyzed the catalytic properties and solved the crystal structure of phenylalanine hydroxylase from the psychrophilic bacterium Colwellia psychrerythraea 34H (CpPAH). CpPAH displays highest activity with tetrahydrobiopterin (BH(4)) as cofactor and at 25 degrees C (15 degrees C above the optimal growth temperature). Although the enzyme is monomeric with a single L-Phe-binding site, the substrate binds cooperatively. In comparison with PAH from mesophilic bacteria and mammalian organisms, CpPAH shows elevated [S(0.5)](L-Phe) (= 1.1 +/- 0.1 mm) and K(m)(BH(4))(= 0.3 +/- 0.1 mm), as well as high catalytic efficiency at 10 degrees C. However, the half-inactivation and denaturation temperature is only slightly lowered (T(m) approximately 52 degrees C; where T(m) is half-denaturation temperature), in contrast to other cold-adapted enzymes. The crystal structure shows regions of local flexibility close to the highly solvent accessible binding sites for BH(4) (Gly(87)/Phe(88)/Gly(89)) and l-Phe (Tyr(114)-Pro(118)). Normal mode and COREX analysis also detect these and other areas with high flexibility. Greater mobility around the active site and disrupted hydrogen bonding abilities for the cofactor appear to represent cold-adaptive properties that do not markedly affect the thermostability of CpPAH.

DFprot: a webtool for predicting local chain deformability

DFprot is a web-based server for predicting main-chain deformability from a single protein conformation. The server automatically performs a normal-mode analysis (NMA) of the uploaded structure and calculates its capability to deform at each of its residues. Non-specialists can easily and rapidly obtain a quantitative first approximation of the flexibility of their structures with a simple and efficient interface.

AVAILABILITY

http://sbg.cib.csic.es/Software/DFprot.

Impact of EGFR point mutations on the sensitivity to gefitinib: insights from comparative structural analyses and molecular dynamics simulations

Emergence of resistant mutations in drug targets represents a serious problem in the targeted chemotherapy. One challenging issue is to understand the atomic-detailed effect of the mutation on the target. Another intriguing issue is how to predict specific mutations that would show up in the clinical setting, leading to drug resistance. By computational approaches, we have investigated structural, dynamics and energetic effects of a series of EGFR mutations identified from the lung cancer patients. We demonstrated mutation L858R caused gefitinib move closer to the hinge region, whereas T790M caused the ligand escape from the binding pocket. In particular, the T790M decreased the size of the hydrophobic slot formed by L718 and G796. This suggests that, to be effective against the T790M mutant, the inhibitors should avoid interactions with the hydrophobic slot. Mutations T790M, L858R, and their combinations are found to cause different conformational redistribution and to perturb the electrostatic potential at the ATP-binding pocket. Normal mode analysis revealed the mutations resulted in changes in the correlated movements in the protein. In an attempt to develop a computational descriptor for predicting the functional effect of EGFR mutations, we have developed a Plarm algorithm, and the Plarm score was found to be an excellent predictor of the functional impact of six clinical relevant mutations in EGFR tyrosine kinase domains, including T790M, L858R, G719C, L861Q, T790M + L858R double mutant, and delL747-P753insS. The Plarm algorithm could be readily extended to investigate other drug targets.

Conformational Flexibility, Hydration and State Parameter Fluctuations of Fibroblast Growth Factor-10: Effects of Ligand Binding

Differential effects of ligand binding on local and global fibroblast growth factor-10 (FGF-10) flexibility and stability have been investigated utilizing a variety of experimental and computational techniques. Normal mode analysis was used to predict the low frequency motions and regional flexibility of FGF-10. Similarly, regional variations in local folding/unfolding equilibria were characterized with the COREX/BEST algorithm. Experimental adiabatic and isothermal compressibilities of FGF-10 alone and in the presence of polyanions are compared. Furthermore, the effect of polyanions on the coefficient of thermal expansion is compared. Measurements of density, heat capacity, compressibility, and expansibility were combined to calculate experimentally determined volume and enthalpy fluctuations. Global effects of polyanions on FGF-10 flexibility, thermodynamic fluctuations, and hydration vary depending on the size and charge density of the polyanion. Local effects of polyanions were investigated utilizing time-resolved fluorescence spectroscopy and red edge excitation spectroscopy (REES). Increased rigidity of the protein matrix or an increased solvent response surrounding the Trp residues is observed in the presence of polyanions. Similarly, time-resolved spectroscopy reveals increased ground state heterogeneity and increased dipole relaxation on the time scale of fluorescence for FGF-10 in the presence of polyanions. These polyanions increase heterogeneity, global flexibility, and fluctuations while increasing the melting temperature (Tm) of FGF-10.

NOMAD-Ref: visualization, deformation and refinement of macromolecular structures based on all-atom normal mode analysis

Normal mode analysis (NMA) is an efficient way to study collective motions in biomolecules that bypasses the computational costs and many limitations associated with full dynamics simulations. The NOMAD-Ref web server presented here provides tools for online calculation of the normal modes of large molecules (up to 100 000 atoms) maintaining a full all-atom representation of their structures, as well as access to a number of programs that utilize these collective motions for deformation and refinement of biomolecular structures. Applications include the generation of sets of decoys with correct stereochemistry but arbitrary large amplitude movements, the quantification of the overlap between alternative conformations of a molecule, refinement of structures against experimental data, such as X-ray diffraction structure factors or Cryo-EM maps and optimization of docked complexes by modeling receptor/ligand flexibility through normal mode motions. The server can be accessed at the URL .

oGNM: online computation of structural dynamics using the Gaussian Network Model

An assessment of the equilibrium dynamics of biomolecular systems, and in particular their most cooperative fluctuations accessible under native state conditions, is a first step towards understanding molecular mechanisms relevant to biological function. We present a web-based system, oGNM that enables users to calculate online the shape and dispersion of normal modes of motion for proteins, oligonucleotides and their complexes, or associated biological units, using the Gaussian Network Model (GNM). Computations with the new engine are 5–6 orders of magnitude faster than those using conventional normal mode analyses. Two cases studies illustrate the utility of oGNM. The first shows that the thermal fluctuations predicted for 1250 non-homologous proteins correlate well with X-ray crystallographic data over a broad range [7.3–15 Å] of inter-residue interaction cutoff distances and the correlations improve with increasing observation temperatures. The second study, focused on 64 oligonucleotides and oligonucleotide–protein complexes, shows that good agreement with experiments is achieved by representing each nucleotide by three GNM nodes (as opposed to one-node-per-residue in proteins) along with uniform interaction ranges for all components of the complexes. These results open the way to a rapid assessment of the dynamics of DNA/RNA-containing complexes. The server can be accessed at .

UMMS: constrained harmonic and anharmonic analyses of macromolecules based on elastic network models

UMass Morph Server (UMMS) has been developed for the broad impact on the study of molecular dynamics (MD). The elastic network model (ENM) of a given macromolecule has been proven as a useful tool for analyzing thermal behaviors locally and predicting folding pathways globally. UMMS utilizes coarse-grained ENMs at various levels. These simplifications remarkably save computation time compared with all-atom MD simulations so that one can bring down massive computational problems from a supercomputer to a PC. To improve computational efficiency and physical reality of ENMs, the symmetry-constrained, rigid-cluster, hybrid and chemical-bond ENMs have been developed and implemented at UMMS. One can request both harmonic normal mode analysis of a single macromolecule and anharmonic pathway generation between two conformations of a same molecule using elastic network interpolation at http://biomechanics.ecs.umass.edu/umms.html.

Coarse-grained normal mode analysis in structural biology

The realization that experimentally observed functional motions of proteins can be predicted by coarse-grained normal mode analysis has renewed interest in applications to structural biology. Notable applications include the prediction of biologically relevant motions of proteins and supramolecular structures driven by their structure-encoded collective dynamics; the refinement of low-resolution structures, including those determined by cryo-electron microscopy; and the identification of conserved dynamic patterns and mechanically key regions within protein families. Additionally, hybrid methods that couple atomic simulations with deformations derived from coarse-grained normal mode analysis are able to sample collective motions beyond the range of conventional molecular dynamics simulations. Such applications have provided great insight into the underlying principles linking protein structures to their dynamics and their dynamics to their functions.

Characterization of protein matrix motions in the Rb. sphaeroides photosynthetic reaction center

We use Normal Mode Analysis to investigate motions in the photosynthetic reaction center (RC) protein. We identify the regions involved in concerted fluctuations of the protein matrix and analyze the normalized amplitudes and the directionality of the first few dominant modes. We also seek to quantify the coupling of normal modes to long-range electron transfer (ET). We find that a quasi-continuous spectrum of protein motions rather than one individual mode contributes to light-driven electron transfer. This is consistent with existing theoretical models (e.g. the spin-boson/dispersed polaron model) for the coupling of the protein and solvent "bath" to charge separation events. [Figure: see text].