Generating Triangulated Macromolecular Surfaces by Euclidean Distance Transform

Macromolecular surfaces are fundamental representations of their three-dimensional geometric shape. Accurate calculation of protein surfaces is of critical importance in the protein structural and functional studies including ligand-protein docking and virtual screening. In contrast to analytical or parametric representation of macromolecular surfaces, triangulated mesh surfaces have been proved to be easy to describe, visualize and manipulate by computer programs. Here, we develop a new algorithm of EDTSurf for generating three major macromolecular surfaces of van der Waals surface, solvent-accessible surface and molecular surface, using the technique of fast Euclidean Distance Transform (EDT). The triangulated surfaces are constructed directly from volumetric solids by a Vertex-Connected Marching Cube algorithm that forms triangles from grid points. Compared to the analytical result, the relative error of the surface calculations by EDTSurf is <2–4% depending on the grid resolution, which is 1.5–4 times lower than the methods in the literature; and yet, the algorithm is faster and costs less computer memory than the comparative methods. The improvements in both accuracy and speed of the macromolecular surface determination should make EDTSurf a useful tool for the detailed study of protein docking and structure predictions. Both source code and the executable program of EDTSurf are freely available at http://zhang.bioinformatics.ku.edu/EDTSurf

Construction, MD Simulation, and Hydrodynamic Validation of an All-Atom Model of a Monoclonal IgG Antibody

AbstractAt 150 kDa, antibodies of the IgG class are too large for their structure to be determined with current NMR methodologies. Because of hinge-region flexibility, it is difficult to obtain atomic-level structural information from the crystal, and questions regarding antibody structure and dynamics in solution remain unaddressed. Here we describe the construction of a model of a human IgG1 monoclonal antibody (trastuzumab) from the crystal structures of fragments. We use a combination of molecular-dynamics (MD) simulation, continuum hydrodynamics modeling, and experimental diffusion measurements to explore antibody behavior in aqueous solution. Hydrodynamic modeling provides a link between the atomic-level details of MD simulation and the size- and shape-dependent data provided by hydrodynamic measurements. Eight independent 40 ns MD trajectories were obtained with the AMBER program suite. The ensemble average of the computed transport properties over all of the MD trajectories agrees remarkably well with the value of the translational diffusion coefficient obtained with dynamic light scattering at 20°C and 27°C, and the intrinsic viscosity measured at 20°C. Therefore, our MD results likely represent a realistic sampling of the conformational space that an antibody explores in aqueous solution

How round is a protein? Exploring protein structures for globularity using conformal mapping.

We present a new algorithm that automatically computes a measure of the geometric difference between the surface of a protein and a round sphere. The algorithm takes as input two triangulated genus zero surfaces representing the protein and the round sphere, respectively, and constructs a discrete conformal map f between these surfaces. The conformal map is chosen to minimize a symmetric elastic energy E S (f) that measures the distance of f from an isometry. We illustrate our approach on a set of basic sample problems and then on a dataset of diverse protein structures. We show first that E S (f) is able to quantify the roundness of the Platonic solids and that for these surfaces it replicates well traditional measures of roundness such as the sphericity. We then demonstrate that the symmetric elastic energy E S (f) captures both global and local differences between two surfaces, showing that our method identifies the presence of protruding regions in protein structures and quantifies how these regions make the shape of a protein deviate from globularity. Based on these results, we show that E S (f) serves as a probe of the limits of the application of conformal mapping to parametrize protein shapes. We identify limitations of the method and discuss its extension to achieving automatic registration of protein structures based on their surface geometry

Improvements to the APBS biomolecular solvation software suite

The Adaptive Poisson-Boltzmann Solver (APBS) software was developed to solve the equations of continuum electrostatics for large biomolecular assemblages that has provided impact in the study of a broad range of chemical, biological, and biomedical applications. APBS addresses three key technology challenges for understanding solvation and electrostatics in biomedical applications: accurate and efficient models for biomolecular solvation and electrostatics, robust and scalable software for applying those theories to biomolecular systems, and mechanisms for sharing and analyzing biomolecular electrostatics data in the scientific community. To address new research applications and advancing computational capabilities, we have continually updated APBS and its suite of accompanying software since its release in 2001. In this manuscript, we discuss the models and capabilities that have recently been implemented within the APBS software package including: a Poisson-Boltzmann analytical and a semi-analytical solver, an optimized boundary element solver, a geometry-based geometric flow solvation model, a graph theory based algorithm for determining p$K_a$ values, and an improved web-based visualization tool for viewing electrostatics

Solvent accessible surface area approximations for rapid and accurate protein structure prediction

The burial of hydrophobic amino acids in the protein core is a driving force in protein folding. The extent to which an amino acid interacts with the solvent and the protein core is naturally proportional to the surface area exposed to these environments. However, an accurate calculation of the solvent-accessible surface area (SASA), a geometric measure of this exposure, is numerically demanding as it is not pair-wise decomposable. Furthermore, it depends on a full-atom representation of the molecule. This manuscript introduces a series of four SASA approximations of increasing computational complexity and accuracy as well as knowledge-based environment free energy potentials based on these SASA approximations. Their ability to distinguish correctly from incorrectly folded protein models is assessed to balance speed and accuracy for protein structure prediction. We find the newly developed “Neighbor Vector” algorithm provides the most optimal balance of accurate yet rapid exposure measures

Improving the resolution of interaction maps: A middleground between high-resolution complexes and genome-wide interactomes

Protein-protein interactions are ubiquitous in Biology and therefore central to understand living organisms. In recent years, large-scale studies have been undertaken to describe, at least partially, protein-protein interaction maps or interactomes for a number of relevant organisms including human. Although the analysis of interaction networks is proving useful, current interactomes provide a blurry and granular picture of the molecular machinery, i.e. unless the structure of the protein complex is known the molecular details of the interaction are missing and sometime is even not possible to know if the interaction between the proteins is direct, i.e. physical interaction or part of functional, not necessary, direct association. Unfortunately, the determination of the structure of protein complexes cannot keep pace with the discovery of new protein-protein interactions resulting in a large, and increasing, gap between the number of complexes that are thought to exist and the number for which 3D structures are available. The aim of the thesis was to tackle this problem by implementing computational approaches to derive structural models of protein complexes and thus reduce this existing gap. Over the course of the thesis, a novel modelling algorithm to predict the structure of protein complexes, V-D2OCK, was implemented. This new algorithm combines structure-based prediction of protein binding sites by means of a novel algorithm developed over the course of the thesis: VORFFIP and M-VORFFIP, data-driven docking and energy minimization. This algorithm was used to improve the coverage and structural content of the human interactome compiled from different sources of interactomic data to ensure the most comprehensive interactome. Finally, the human interactome and structural models were compiled in a database, V-D2OCK DB, that offers an easy and user-friendly access to the human interactome including a bespoken graphical molecular viewer to facilitate the analysis of the structural models of protein complexes. Furthermore, new organisms, in addition to human, were included providing a useful resource for the study of all known interactomes

Conformational Dynamics of Supramolecular Protein Assemblies in the EMDB

The Electron Microscopy Data Bank (EMDB) is a rapidly growing repository for the dissemination of structural data from single-particle reconstructions of supramolecular protein assemblies including motors, chaperones, cytoskeletal assemblies, and viral capsids. While the static structure of these assemblies provides essential insight into their biological function, their conformational dynamics and mechanics provide additional important information regarding the mechanism of their biological function. Here, we present an unsupervised computational framework to analyze and store for public access the conformational dynamics of supramolecular protein assemblies deposited in the EMDB. Conformational dynamics are analyzed using normal mode analysis in the finite element framework, which is used to compute equilibrium thermal fluctuations, cross-correlations in molecular motions, and strain energy distributions for 452 of the 681 entries stored in the EMDB at present. Results for the viral capsid of hepatitis B, ribosome-bound termination factor RF2, and GroEL are presented in detail and validated with all-atom based models. The conformational dynamics of protein assemblies in the EMDB may be useful in the interpretation of their biological function, as well as in the classification and refinement of EM-based structures.Comment: Associated online data bank available at: http://lcbb.mit.edu/~em-nmdb

VASCo: computation and visualization of annotated protein surface contacts

<p>Abstract</p> <p>Background</p> <p>Structural data from crystallographic analyses contain a vast amount of information on protein-protein contacts. Knowledge on protein-protein interactions is essential for understanding many processes in living cells. The methods to investigate these interactions range from genetics to biophysics, crystallography, bioinformatics and computer modeling. Also crystal contact information can be useful to understand biologically relevant protein oligomerisation as they rely in principle on the same physico-chemical interaction forces. Visualization of crystal and biological contact data including different surface properties can help to analyse protein-protein interactions.</p> <p>Results</p> <p>VASCo is a program package for the calculation of protein surface properties and the visualization of annotated surfaces. Special emphasis is laid on protein-protein interactions, which are calculated based on surface point distances. The same approach is used to compare surfaces of two aligned molecules. Molecular properties such as electrostatic potential or hydrophobicity are mapped onto these surface points. Molecular surfaces and the corresponding properties are calculated using well established programs integrated into the package, as well as using custom developed programs. The modular package can easily be extended to include new properties for annotation. The output of the program is most conveniently displayed in PyMOL using a custom-made plug-in.</p> <p>Conclusion</p> <p>VASCo supplements other available protein contact visualisation tools and provides additional information on biological interactions as well as on crystal contacts. The tool provides a unique feature to compare surfaces of two aligned molecules based on point distances and thereby facilitates the visualization and analysis of surface differences.</p

Structure and dynamics of biomacromolecules in solution: recent developments and future perspectives in SANS/SAXS and neutron spectroscopy

This French University Habilitation (“mémoire” to obtain the “Diplôme d’Habilitation à Diriger des Recherches”, DHDR) is divided into two parts: the first one deals with the results that I have obtained after my PhD thesis in 2003, the second one discusses open questions related to theseresults as well as mid- and long-term perspectives.Three different topics are presented in the result chapters: 1) a combination of small angle scattering (SAS) and nuclear magnetic resonance (NMR) for rigid-body modeling of biomacromolecular complexes, 2) the combined use of small angle X-ray (SAXS) and neutron (SANS) scattering for the study of unfolded proteins and 3) the study of biomacromolecular and solvent dynamics by neutron spectroscopy combining several instruments.Points 1) and 3) are discussed in great detail, both in the results and perspective sections, since they represent the most advanced projects of my research. Point 2) is dealt with more briefly in the results section, since few results are available so far. However, perspectives are discussed. Aspecial perspective chapter deals with applications on membrane proteins. Key publications for the different chapters are:Chapter 1: Gabel et al. (2006) A target function for quaternary structural refinement from small angle scattering and NMR orientational restraints. Eur. Biophys. J. 35(4), 313-327. Gabel et al. (2008) A structure refinement protocol combining NMR residual dipolar couplings and small angle scattering restraints. J. Biomol. NMR 41(4), 199-208.Chapter 2: Gabel et al. (2009) Quantitative Modelfree Analysis of Urea Binding to Unfolded Ubiquitin Using a Combination of Small Angle X-ray and Neutron Scattering. J. Am. Chem. Soc. 131(25), 8769-8771.Chapter 3: Gabel (2005) Protein dynamics in solution and powder measured by incoherent elastic neutron scattering: the influence of Q-range and energy resolution. Eur. Biophys. J. 31(1), 1-12. Gabel & Bellissent-Funel (2007) C-Phycocyanin Hydration Water Dynamics in the Presence of Trehalose: An Incoherent Elastic Neutron Scattering Study at Different Energy Resolutions. Biophys. J. 92(11), 4054-4063.The habilitation thesis is focused on methodological aspects and developments of small angle scattering and neutron spectroscopy. It is obvious that the approaches discussed here and their sophisticated levels of data analysis rely fundamentally on the quality of the sample, and inparticular on monodispersity (for SAS) and amount of material for spectroscopy. The paramount importance of good biochemistry and the use of complementary techniques for the characterization of samples can hardly be overestimated. They include, amongst others, gelfiltration, analytical ultracentrifugation, static and dynamic light scattering, NMR, etc. They are quite simply indispensable for doing good and accurate science with SAS, in particular in more complex systems (macromolecular complexes, membrane proteins …). If they are not presentedin more detail in this thesis, it is not out of ignorance of this fact but due to the lack of space