Source author record

Rafael Najmanovich

Rafael Najmanovich appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Biomolecules Genomics Populations and Evolution

Catalog footprint

What is connected

4works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

The DynaSig-ML Python package: automated learning of biomolecular dynamics-function relationships

Summary: The DynaSig-ML (Dynamical Signatures - Machine Learning) Python package allows the efficient, user-friendly exploration of 3D dynamics-function relationships in biomolecules, using datasets of experimental measures from large numbers of sequence variants. The DynaSig-ML package is built around the Elastic Network Contact Model (ENCoM), the first and only sequence-sensitive coarse-grained NMA model, which is used to generate the input Dynamical Signatures. Starting from in silico mutated structures, the whole pipeline can be run with just a few lines of Python and modest computational resources. The compute-intensive steps can also easily be parallelized in the case of either large biomolecules or vast amounts of sequence variants. As an example application, we use the DynaSig-ML package to predict the evolutionary fitness of the bacterial enzyme VIM-2 lactamase from deep mutational scan data. Availability and implementation: DynaSig-ML is open source software available at https: //github.com/gregorpatof/dynasigml package Contact: rafael.najmanovich@umontreal.ca

preprint2020arXiv

The NRGTEN Python package: an extensible toolkit for coarse-grained normal mode analysis of proteins, nucleic acids, small molecules and their complexes

Summary: Coarse-grained normal mode analysis (NMA) is a fast computational technique to study the dynamics of biomolecules. Here we present the Najmanovich Research Group Toolkit for Elastic Networks (NRGTEN). NRGTEN is a Python toolkit that implements four different NMA models in addition to popular and novel metrics to benchmark and measure properties from these models. Furthermore, the toolkit is available as a public Python package and is easily extensible for the development or implementation of additional NMA models. The inclusion of the ENCoM model (Elastic Network Contact Model) developed in our group within NRGTEN is noteworthy, owing to its account for the specific chemical nature of atomic interactions. This makes possible some unique predictions of the effect of mutations, such as on stability (via changes in vibrational entropy differences), on the transition probability between different conformational states or on the flexibility profile of the whole macromolecule/complex (to study allostery and signalling). In addition, all NMA models can be used to generate conformational ensembles from a starting structure to aid in protein-protein, protein-ligand or other docking studies among applications. NRGTEN is freely available via a public Python package which can be easily installed on any modern machine and includes a detailed user guide hosted online. Availability and implementation: https://github.com/gregorpatof/nrgten_package/ Contact: rafael.najmanovich@umontreal.ca

preprint2013arXiv

Protein flexibility upon ligand binding: Docking predictions and statistical analysis

Side chain flexibility is an important factor in ligand binding. In order to determine the extent to which side chain flexibility is involved in ligand binding, a knowledge-based approach was taken. A database composed of examples of protein structures in the presence or absence of a given ligand is used to analyze which side chains undergo side chain conformational changes. Such an analysis has determined that up to 40% of binding site do not present side chain conformational changes. A total of three residues undergoing side chain conformational changes encompass approximately 85% of the binding sites studied. When analyzing the propensities of different amino acids to undergo side chain conformational changes we find that there are considerable differences between different amino acids. A support vector machine learning approach was used to create a classifier system utilizing information about the solvent accessible area as well as flexibility scale value of each specific side chain to be predicted together with its neighboring side chains. An accuracy level of 70% is reached using this approach. The fact that a small number of residues undergo side chain conformational changes in the majority of binding sites makes it feasible to introduce side chain flexibility in docking simulations. An algorithm has been developed for introducing side chain flexibility utilizing a hybrid genetic-algorithm/exhaustive-search procedure and a surface complementarity based scoring function. This approach is implemented in the software tool FlexAID. FlexAID utilizes a rotamer library to create alternative conformations for a list of residues that are exhaustively searched during the docking simulation. The performance of FlexAID in rigid local as well as global simulations falls in the 70-80% range for both local and global simulations.

preprint2012arXiv

Deep-sequencing of the Peach Latent Mosaic Viroid Reveals New Aspects of Population Heterogeneity

Viroids are small circular single-stranded infectious RNAs that are characterized by a relatively high mutation level. Knowledge of their sequence heterogeneity remains largely elusive, and, as yet, no strategy attempting to address this question from a population dynamics point of view is in place. In order to address these important questions, a GF305 indicator peach tree was infected with a single variant of the Avsunviroidae family member Peach latent mosaic viroid (PLMVd). Six months post-inoculation, full-length circular conformers of PLMVd were isolated, deep-sequenced and the resulting sequences analyzed using an original bioinformatics scheme specifically designed and developed in order to evaluate the richness of a given the sequence's population. Two distinct libraries were analyzed, and yielded 1125 and 1061 different PLMVd variants respectively, making this study the most productive to date (by more than an order of magnitude) in terms of the reporting of novel viroid sequences. Sequence variants exhibiting up to ~20% of mutations relative to the inoculated viroid were retrieved, clearly illustrating the high divergence dynamic inside a unique population. Using a novel hierarchical clustering algorithm, the different variants obtained were grouped into either 7 or 8 clusters depending on the library being analyzed. Most of the sequences contained, on average, between 4.6 and 6.3 mutations relative to the variant used initially to inoculate the plant. Interestingly, it was possible to reconstitute the sequence evolution between these clusters. On top of providing a reliable pipeline for the treatment of viroid deep-sequencing, this study sheds new light on the importance of the sequence variation that may take place in a viroid population and which may result in the formation of a quasi-species.