FlexServ: an integrated tool for the analysis of protein flexibility.

SUMMARY: FlexServ is a web-based tool for the analysis of protein flexibility. The server incorporates powerful protocols for the coarse-grained determination of protein dynamics using different versions of Normal Mode Analysis (NMA), Brownian dynamics (BD) and Discrete Dynamics (DMD). It can also analyze user provided trajectories. The server allows a complete analysis of flexibility using a large variety of metrics, including basic geometrical analysis, B-factors, essential dynamics, stiffness analysis, collectivity measures, Lindemann's indexes, residue correlation, chain-correlations, dynamic domain determination, hinge point detections, etc. Data is presented through a web interface as plain text, 2D and 3D graphics. AVAILABILITY: http://mmb.pcb.ub.es/FlexServ; http://www.inab.org

Exploring the suitability of coarse-grained techniques for the representation of protein dynamics.

A systematic study of two coarse-grained techniques for the description of protein dynamics is presented. The two techniques exploit either Brownian or discrete molecular dynamics algorithms applied in the context of simple C(alpha)-C(alpha) potentials, like those used in coarse-grained normal mode analysis. Coarse-grained simulations of the flexibility of protein metafolds are compared to those computed with fully atomistic molecular dynamics simulations using state-of-the-art physical potentials and explicit solvent. Both coarse-grained models efficiently capture critical features of the protein dynamics.

GRID-MD-A tool for massive simulation of protein channels.

We present here a fast method for the exploration of channels in proteins based on molecular dynamics simulations of probe particles in a discrete grid space defined by an ensemble of protein conformations obtained either experimentally or by out-of-the-grid atomistic molecular dynamics simulations. The method is able to provide millisecond-long trajectories with a small computational effort, requires no human intervention in defining possible exit pathways and can detect both major and minor channels, giving a correct balance to the relative flux between them. The Grid-Molecular-Dynamics approach is then a suitable method for massive exploration of channels in proteins, even of those with unknown functional annotation.

Prediction and validation of protein intermediate states from structurally rich ensembles and coarse-grained simulations.

Protein conformational changes are at the heart of cell functions, from signalling to ion transport. However, the transient nature of the intermediates along transition pathways hampers their experimental detection, making the underlying mechanisms elusive. Here we retrieve dynamic information on the actual transition routes from principal component analysis (PCA) of structurally-rich ensembles and, in combination with coarse-grained simulations, explore the conformational landscapes of five well-studied proteins. Modelling them as elastic networks in a hybrid elastic-network Brownian dynamics simulation (eBDIMS), we generate trajectories connecting stable end-states that spontaneously sample the crystallographic motions, predicting the structures of known intermediates along the paths. We also show that the explored non-linear routes can delimit the lowest energy passages between end-states sampled by atomistic molecular dynamics. The integrative methodology presented here provides a powerful framework to extract and expand dynamic pathway information from the Protein Data Bank, as well as to validate sampling methods in general.

Fast Atomistic Molecular Dynamics Simulations from Essential Dynamics Samplings.

We present a new method for fast molecular dynamics simulations in cases where the new trajectories can be considered a perturbation or a combination of previously stored ones. The method is designed for the postgenomic scenario, where databases such as MoDEL will store curated equilibrium trajectories of all biomolecules (proteins, nucleic acids, etc.) of human interest. We demonstrate that the approach outlined here can, with accuracy and great computational efficiency, reproduce and extend original trajectories, describe dynamical effects due to perturbations (e.g., protein-ligand and protein-protein interactions and protein mutations) and predict the dynamics of large polymeric systems built up from previously studied fragments. The method can work simultaneously with low- and high-resolution pictures of the macromolecule, allowing the level of detail to be matched to that required for obtaining the information of biological interest.

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.

Coarse-grained representation of protein flexibility. Foundations, successes, and shortcomings.

Flexibility is the key magnitude to understand the variety of functions of proteins. Unfortunately, its experimental study is quite difficult, and in fact, most experimental procedures are designed to reduce flexibility and allow a better definition of the structure. Theoretical approaches have become then the alternative but face serious timescale problems, since many biologically relevant deformation movements happen in a timescale that is far beyond the possibility of current atomistic models. In this complex scenario, coarse-grained simulation methods have emerged as a powerful and inexpensive alternative. Along this chapter, we will review these coarse-grained methods, and explain their physical foundations and their range of applicability.

MoDEL (Molecular Dynamics Extended Library): a database of atomistic molecular dynamics trajectories.

More than 1700 trajectories of proteins representative of monomeric soluble structures in the protein data bank (PDB) have been obtained by means of state-of-the-art atomistic molecular dynamics simulations in near-physiological conditions. The trajectories and analyses are stored in a large data warehouse, which can be queried for dynamic information on proteins, including interactions. Here, we describe the project and the structure and contents of our database, and provide examples of how it can be used to describe the global flexibility properties of proteins. Basic analyses and trajectories stripped of solvent molecules at a reduced resolution level are available from our web server.

Structural plasticity and functional implications of internal cavities in distal mutants of type 1 non-symbiotic hemoglobin AHb1 from Arabidopsis thaliana.

The increasing number of nonsymbiotic plant hemoglobins discovered in genomic studies in the past decade raises intriguing questions about their physiological role. Among them, the nonsymbiotic hemoglobin AHb1 from Arabidopsis thaliana deserves particular attention, as it combines an extremely high oxygen affinity with an internal hexacoordination of the distal histidine HisE7 to the heme iron in the absence of exogenous ligands. In order to gain insight into the structure-function relationships of the protein, the ligand binding properties of mutants of two conserved residues of the distal cavity, HisE7 --> Leu and PheB10 --> Leu, were investigated by experimental and computational studies and compared to results determined for the wild type (wt) protein. The Fe(2+)-deoxy HisE7 --> Leu mutant exists, as expected, in the pentacoordinated form, while a mixture of penta- and hexacoordinated forms is found for the PheB10 --> Leu mutant, with an equilibrium shifted toward the pentacoordinated form with respect to the wt protein. Spectroscopic studies of the complexes of CO and CN(-) with AHb1 and its mutants show a subtle interplay of steric and electrostatic effects by distal residues on the ligand binding to the heme. Moreover, stopped-flow and flash photolysis experiments reveal substantial kinetic differences triggered by those mutations, which are particularly manifested in the enhanced geminate rebinding and bimolecular association rate. These findings are discussed in light of the drastic alterations found by molecular dynamics simulations in the nature and distribution of internal cavities in the protein matrix of the mutants, revealing an extremely large sensitivity of the protein structure to changes in distal HisE7 and PheB10 residues. Overall, data are consistent with the putative NO-dioxygenase activity attributed to AHb1.