Julio A. Kovacs

ADP_EM: fast exhaustive multi-resolution docking for high-throughput coverage

MOTIVATION Efficient fitting tools are needed to take advantage of a fast growth of atomic models of protein domains from crystallography or comparative modeling, and low-resolution density maps of larger molecular assemblies. Here, we report a novel fitting algorithm for the exhaustive and fast overlay of partial high-resolution models into a low-resolution density map. The method incorporates a fast rotational search based on spherical harmonics (SH) combined with a simple translational scanning. RESULTS This novel combination makes it possible to accurately dock atomic structures into low-resolution electron-density maps in times ranging from seconds to a few minutes. The high-efficiency achieved with simulated and experimental test cases preserves the exhaustiveness needed in these heterogeneous-resolution merging tools. The results demonstrate its efficiency, robustness and high-throughput coverage. AVAILABILITY http://sbg.cib.csic.es/Software/ADP_EM. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Predictions of Protein Flexibility: First-Order Measures

The normal modes of a molecule are utilized, in conjunction with classical conformal vector field theory, to define a function that measures the capability of the molecule to deform at each of its residues. An efficient algorithm is presented to calculate the local chain deformability from the set of normal modes of vibration. This is done by considering each mode as an off‐grid sample of a deformation vector field. Predictions of deformability are compared with experimental data in the form of dihedral angle differences between two conformations of ten kinases by using a modified correlation function. Deformability calculations correlate well with experimental results and validate the applicability of this method to protein flexibility predictions. Proteins 2004.

Topology representing neural networks reconcile biomolecular shape, structure, and dynamics

Topology-representing networks (TRNs) generate reduced models of biomolecules and thereby facilitate the 5tting of molecular fragments into large macromolecular complexes. The components of such complexes undergo a wide range of motions, and shapes observed at low resolution often deviate from the known atomic structures. What is required for the modeling of such motions is a combination of global shape constraints based on the low-resolution data with a local modeling of atomic interactions. We present a novel Motion Capture Network that freezes inessential degrees of freedom to maintain the stereochemistryof an atomic model. TRN-based deformable models retain much of the mechanical properties of biological macromolecules. The elastic models yield a decomposition of the predicted motion into vibrational normal modes and are amenable to interactive manipulation with haptic rendering software. c

DFprot: A webtool for predicting local chain deformability

UNLABELLED 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.

Damped-Dynamics Flexible Fitting

In fitting atomic structures into EM maps, it often happens that the map corresponds to a different conformation of the structure. We have developed a new methodology to handle these situations that preserves the covalent geometry of the structure and allows the modeling of large deformations. The first goal is achieved by working in generalized coordinates (positional and internal coordinates), and the second by avoiding harmonic potentials. Instead, we use dampers (shock absorbers) between every pair of atoms, combined with a force field that attracts the atomic structure toward incompletely occupied regions of the EM map. The trajectory obtained by integrating the resulting equations of motion converges to a conformation that, in our validation cases, was very close to the target atomic structure. Compared to current methods, our approach is more efficient and robust against wrong solutions and to overfitting, and does not require user intervention or subjective decisions. Applications to the computation of transition pathways between known conformers, homology and loop modeling, as well as protein docking, are also discussed.

Comparing an Atomic Model or Structure to a Corresponding Cryo-electron Microscopy Image at the Central Axis of a Helix

Three-dimensional density maps of biological specimens from cryo-electron microscopy (cryo-EM) can be interpreted in the form of atomic models that are modeled into the density, or they can be compared to known atomic structures. When the central axis of a helix is detectable in a cryo-EM density map, it is possible to quantify the agreement between this central axis and a central axis calculated from the atomic model or structure. We propose a novel arc-length association method to compare the two axes reliably. This method was applied to 79 helices in simulated density maps and six case studies using cryo-EM maps at 6.4-7.7 Å resolution. The arc-length association method is then compared to three existing measures that evaluate the separation of two helical axes: a two-way distance between point sets, the length difference between two axes, and the individual amino acid detection accuracy. The results show that our proposed method sensitively distinguishes lateral and longitudinal discrepancies between the two axes, which makes the method particularly suitable for the systematic investigation of cryo-EM map-model pairs.

Tracing Actin Filament Bundles in Three-Dimensional Electron Tomography Density Maps of Hair Cell Stereocilia

Cryo-electron tomography (cryo-ET) is a powerful method of visualizing the three-dimensional organization of supramolecular complexes, such as the cytoskeleton, in their native cell and tissue contexts. Due to its minimal electron dose and reconstruction artifacts arising from the missing wedge during data collection, cryo-ET typically results in noisy density maps that display anisotropic XY versus Z resolution. Molecular crowding further exacerbates the challenge of automatically detecting supramolecular complexes, such as the actin bundle in hair cell stereocilia. Stereocilia are pivotal to the mechanoelectrical transduction process in inner ear sensory epithelial hair cells. Given the complexity and dense arrangement of actin bundles, traditional approaches to filament detection and tracing have failed in these cases. In this study, we introduce BundleTrac, an effective method to trace hundreds of filaments in a bundle. A comparison between BundleTrac and manually tracing the actin filaments in a stereocilium showed that BundleTrac accurately built 326 of 330 filaments (98.8%), with an overall cross-distance of 1.3 voxels for the 330 filaments. BundleTrac is an effective semi-automatic modeling approach in which a seed point is provided for each filament and the rest of the filament is computationally identified. We also demonstrate the potential of a denoising method that uses a polynomial regression to address the resolution and high-noise anisotropic environment of the density map.