rDock: A Fast,Versatile and Open Source Program for Docking Ligands to Proteins and Nucleic Acids

Identification of chemical compounds with specific biological activities is an important step in both chemical biology and drug discovery. When the structure of the intended target is available, one approach is to use molecular docking programs to assess the chemical complementarity of small molecules with the target; such calculations provide a qualitative measure of affinity that can be used in virtual screening (VS) to rank order a list of compounds according to their potential to be active. rDock is a molecular docking program developed at Vernalis for high-throughput VS (HTVS) applications. Evolved from RiboDock, the program can be used against proteins and nucleic acids, is designed to be computationally very efficient and allows the user to incorporate additional constraints and information as a bias to guide docking. This article provides an overview of the program structure and features and compares rDock to two reference programs, AutoDock Vina (open source) and Schrödinger''s Glide (commercial). In terms of computational speed for VS, rDock is faster than Vina and comparable to Glide. For binding mode prediction, rDock and Vina are superior to Glide. The VS performance of rDock is significantly better than Vina, but inferior to Glide for most systems unless pharmacophore constraints are used; in that case rDock and Glide are of equal performance. The program is released under the Lesser General Public License and is freely available for download, together with the manuals, example files and the complete test sets, at http://rdock.sourceforge.net/

This is a PLOS Computational Biology Software Article.

     

La desoxycytidine aminohydrolase des feuilles de Zea mays

Deoxycytidine aminohydrolase was extracted from Zea Mays leaves and purified 22-fold. The extract deaminates deoxycytidine, cytidine, 5-methyldeoxycytidine, 5-bromodeoxycytidine and arabinosylcytosine. Its molecular weight was estimated to be 77 000-78 000 based on the elution pattern from Sephadex G100 and G200. Kinetic studies indicated that the hydrolysis of the substrates proceeds at rates which can formally be described by Michaelis-Menten kinetics. p-Chloromercuribenzoate inhibits deoxycytidine aminohydrolase activity, as do mercaptoethanol and dithiothreitol. The activity is also inhibited by deoxyuridine (mixed type), uridine (non competitive), 5-bromodeoxyuridine (competitive).     

TAPASS: Tool for annotation of protein amyloidogenicity in the context of other structural states

Numerous studies have demonstrated that the propensity of a protein to form amyloids or amorphous aggregates is encoded by its amino acid sequence. This led to the emergence of several computational programs to predict amyloidogenicity from amino acid sequences. However, a growing number of studies indicate that an accurate prediction of the protein aggregation can only be achieved when also accounting for the overall structural context of the protein, and the likelihood of transition between the initial state and the aggregate. Here, we describe a computational pipeline called TAPASS, which was designed to do just that. The pipeline assigns each residue of a protein as belonging to a structured region or an intrinsically disordered region (IDR). For this purpose, TAPASS uses either several state-of-the-art programs for prediction of IDRs, of transmembrane regions and of structured domains or the artificial intelligence program AlphaFold. In the next step, this assignment is crossed with amyloidogenicity prediction. As a result, TAPASS allows the detection of Exposed Amyloidogenic Regions (EARs) located within intrinsically disordered regions (IDRs) and carrying high amyloidogenic potential. TAPASS can substantially improve the prediction of amyloids and be used in proteome-wide analysis to discover new amyloid-forming proteins. Its results, combined with clinical data, can create individual risk profiles for different amyloidoses, opening up new opportunities for personalised medicine. The architecture of the pipeline is designed so that it makes it easy to add new individual predictors as they become available. TAPASS can be used through the web interface (https://bioinfo.crbm.cnrs.fr/index.php?route=tools&tool=32).     

MacSECPROT: Two programs in basic to evaluate protein secondary structure

Two interactive programs in BASIC are described, which provide useful tools to evaluate protein secondary structure. Output is given in two formats: (1) graphics are displayed on screen, which can be printed immediately, and (2) textfiles are saved to disk as permanent records and can be printed with a word-processing program. The programs are fast and easy to use and could be a valuable teaching aid in biochemical and molecular biology courses. Program lists are written in Microsoft® BASIC for the Apple® Macintosh™, but can be adapted to other machines accepting graphic commands.     

ClustScan: an integrated program package for the semi-automatic annotation of modular biosynthetic gene clusters and in silico prediction of novel chemical structures

The program package ‘ClustScan’ (Cluster Scanner) is designed for rapid, semi-automatic, annotation of DNA sequences encoding modular biosynthetic enzymes including polyketide synthases (PKS), non-ribosomal peptide synthetases (NRPS) and hybrid (PKS/NRPS) enzymes. The program displays the predicted chemical structures of products as well as allowing export of the structures in a standard format for analyses with other programs. Recent advances in understanding of enzyme function are incorporated to make knowledge-based predictions about the stereochemistry of products. The program structure allows easy incorporation of additional knowledge about domain specificities and function. The results of analyses are presented to the user in a graphical interface, which also allows easy editing of the predictions to incorporate user experience. The versatility of this program package has been demonstrated by annotating biochemical pathways in microbial, invertebrate animal and metagenomic datasets. The speed and convenience of the package allows the annotation of all PKS and NRPS clusters in a complete Actinobacteria genome in 2–3 man hours. The open architecture of ClustScan allows easy integration with other programs, facilitating further analyses of results, which is useful for a broad range of researchers in the chemical and biological sciences.     

The DNA inspector II a program for analyzing and manipulating DNA sequence on the Apple Macintosh

This paper describes a new set of linked programs for analyzing and manipulating DNA sequences using the Apple Macintosh computer. The Macintosh is ideally suited for this purpose because of its extensive use of icons, pull-down menus, and dialog boxes, which form a very intuitive user interface, and its high-resolution graphics, which produce brilliant displays. Because of the extreme ease of use, the learning time for using the program is literally on the order of 5–10 minutes, thus making available computer analyses for many molecular geneticists who previously were not willing to invest the larger amount of time needed to learn how to run the analysis routines in other programs. The DNA Inspector II carries out all of the day-to-day routines needed in a molecular genetics lab. Results of the analyses can be saved to disk, printed, or transferred directly to other applications (e.g., word processing) through the “clipboard, ” an area of memory shared by all applications. The availability of Apple's relatively low-cost LaserWriter printer allows publication-quality figures to be generated very easily from within the program.     

Resveratrol: distribución,propiedades y perspectivas

Resveratrol is a natural polyphenol which can be found in many plants and fruits, such as peanuts, mulberries, blueberries and, above all, in grapes and red wine. Its synthesis is regulated by the presence of stressful factors, such as fungal contamination and ultra-violet radiation. In plants, it plays a role as a phytoalexin, showing a capacity to inhibit the development of certain infections. Plant extracts which contain resveratrol have been employed by traditional medicine for more than 2000 years. Resveratrol was first isolated, and its properties were initially studied with scientific methods, thirty years ago. Its in vitro properties have been extensively studied and demonstrated. It is worth highlighting its activity as an anti-cancer agent, platelet anti-aggregation agent, anti-inflammatory, antiallergenic, etc. The activity of its in vivo properties are not so clear. There are many studies that report benefits on the cardiovascular system, illnesses such as diabetes, and in longevity. However, other authors did not find any agreement between in vitro and in vivo studies. This discrepancy is due to the bioavailability of resveratrol. After an oral dose, it has been demonstrated that the absorption is very high, but the metabolic pathways leave just a little free resveratrol in blood, therefore the bioavailability in the target tissues is very low and the concentrations used in in vitro studies are not found in these tissues. Thus, resveratrol is a very active molecule for maintaining health, but due to the low bioavailability not all the in vitro effects can be translated to in vivo. This opens a new potential approach, seeking derivatives of resveratrol that can be measured in the desired tissues     

fastruct: model‐based clustering made faster

Bayesian model‐based clustering programs have gained increased popularity in studies of population structure since the publication of the software structure . These programs are generally acknowledged as performing well, but their running‐time may be prohibitive. fastruct is a non‐Bayesian implementation of the classical model with no‐admixture uncorrelated allele frequencies. This new program relies on the expectation–maximization principle, and produces assignment rivalling other model‐based clustering programs. In addition, it can be manyfold faster than Bayesian implementations. The software consists of a command‐line engine, which is suitable for batch analysis of data, and a graphical interface, which is convenient for exploring data.     

GONOME: measuring correlations between GO terms and genomic positions

Background

The design of oligonucleotides and PCR primers for studying large genomes is complicated by the redundancy of sequences. The eukaryotic genomes are particularly difficult to study due to abundant repeats. The speed of most existing primer evaluation programs is not sufficient for large-scale experiments.

Results

In order to improve the efficiency and success rate of automatic primer/oligo design, we created a novel method which allows rapid masking of repeats in large sequence files, for example in eukaryotic genomes. It also allows the detection of all alternative binding sites of PCR primers and the prediction of PCR products. The new method was implemented in a collection of efficient programs, the GENOMEMASKER package. The performance of the programs was compared to other similar programs. We also modified the PRIMER3 program, to be able to design primers from lowercase-masked sequences.

Conclusion

The GENOMEMASKER package is able to mask the entire human genome for non-unique primers within 6 hours and find locations of all binding sites for 10 000 designed primer pairs within 10 minutes. Additionally, it predicts all alternative PCR products from large genomes for given primer pairs.     

Development of GMDR-GPU for Gene-Gene Interaction Analysis and Its Application to WTCCC GWAS Data for Type 2 Diabetes

Although genome-wide association studies (GWAS) have identified a significant number of single-nucleotide polymorphisms (SNPs) associated with many complex human traits, the susceptibility loci identified so far can explain only a small fraction of the genetic risk. Among other possible explanations, the lack of a comprehensive examination of gene–gene interaction (G×G) is often considered a source of the missing heritability. Previously, we reported a model-free Generalized Multifactor Dimensionality Reduction (GMDR) approach for detecting G×G in both dichotomous and quantitative phenotypes. However, the computational burden and less efficient implementation of the original programs make them impossible to use for GWAS. In this study, we developed a graphics processing unit (GPU)-based GMDR program (named GWAS-GPU), which is able not only to analyze GWAS data but also to run much faster than the earlier version of the GMDR program. As a demonstration of the program, we used the GMDR-GPU software to analyze a publicly available GWAS dataset on type 2 diabetes (T2D) from the Wellcome Trust Case Control Consortium. Through an exhaustive search of pair-wise interactions and a selected search of three- to five-way interactions conditioned on significant pair-wise results, we identified 24 core SNPs in six genes (FTO: rs9939973, rs9940128, rs9922047, rs1121980, rs9939609, rs9930506; TSPAN8: rs1495377; TCF7L2: rs4074720, rs7901695, rs4506565, rs4132670, rs10787472, rs11196205, rs10885409, rs11196208; L3MBTL3: rs10485400, rs4897366; CELF4: rs2852373, rs608489; RUNX1: rs445984, rs1040328, rs990074, rs2223046, rs2834970) that appear to be important for T2D. Of these core SNPs, 11 in FTO, TSPAN8, and TCF7L2 have been reported to be associated with T2D, obesity, or both, providing an independent replication of previously reported SNPs. Importantly, we identified three new susceptibility genes; i.e., L3MBTL3, CELF4, and RUNX1, for T2D, a finding that warrants further investigation with independent samples.     

Improved analysis of bacterial CGH data beyond the log-ratio paradigm

Background

While the theory of enzyme kinetics is fundamental to analyzing and simulating biochemical systems, the derivation of rate equations for complex mechanisms for enzyme-catalyzed reactions is cumbersome and error prone. Therefore, a number of algorithms and related computer programs have been developed to assist in such derivations. Yet although a number of algorithms, programs, and software packages are reported in the literature, one or more significant limitation is associated with each of these tools. Furthermore, none is freely available for download and use by the community.

Results

We have implemented an algorithm based on the schematic method of King and Altman (KA) that employs the topological theory of linear graphs for systematic generation of valid reaction patterns in a GUI-based stand-alone computer program called KAPattern. The underlying algorithm allows for the assumption steady-state, rapid equilibrium-binding, and/or irreversibility for individual steps in catalytic mechanisms. The program can automatically generate MathML and MATLAB output files that users can easily incorporate into simulation programs.

Conclusion

A computer program, called KAPattern, for generating rate equations for complex enzyme system is a freely available and can be accessed at http://www.biocoda.org.     

A greedier approach for finding tag SNPs

MOTIVATION: Recent studies have shown that a small subset of Single Nucleotide Polymorphisms (SNPs) (called tag SNPs) is sufficient to capture the haplotype patterns in a high linkage disequilibrium region. To find the minimum set of tag SNPs, exact algorithms for finding the optimal solution could take exponential time. On the other hand, approximation algorithms are more efficient but may fail to find the optimal solution. RESULTS: We propose a hybrid method that combines the ideas of the branch-and-bound method and the greedy algorithm. This method explores larger solution space to obtain a better solution than a traditional greedy algorithm. It also allows the user to adjust the efficiency of the program and quality of solutions. This algorithm has been implemented and tested on a variety of simulated and biological data. The experimental results indicate that our program can find better solutions than previous methods. This approach is quite general since it can be used to adapt other greedy algorithms to solve their corresponding problems. AVAILABILITY: The program is available upon request.     

Purification and properties of three intracellular proteinases from Candida albicans

Three intracellular proteinases termed A, B and C were purified to homogeneity from the unicellular form of the yeast Candida albicans. Enzyme A is an aspartic proteinase that acts on a variety of proteins. Its optimal pH is around 5 and it is displaced to 6.5 by KSCN. It is not significantly inhibited by PMSF, TLCK (Tos-Lys-CHCl₂) or soybean trypsin inhibitor but it is inhibited by pepstatin. Its molecular weight is 60 000. Enzyme B is a dipeptidase that acts on esters or on dipeptides without blocks in either the carboxyl or amino ends. Its pH optimum is around 7.5 and the molecular weight is 57 000. It is inhibited by PMSF, TLCK and DANME (N₂Ac-Nle-OMe). Proteinase C is an aminopeptidase with an optimum pH around 8. Its molecular weight was 67 000 when determined by SDS gel electrophoresis and 243 000 when determined by gel filtration. It is active towards dipeptides in which at least one amino acid is apolar and is not active when the N-terminal amino acid is blocked. It is inhibited by EDTA or o-phenanthroline and activated by several divalent cations.     

A Silurian short-great-appendage arthropod

A new arthropod, Enalikter aphson gen. et sp. nov., is described from the Silurian (Wenlock Series) Herefordshire Lagerstätte of the UK. It belongs to the Megacheira (=short-great-appendage group), which is recognized here, for the first time, in strata younger than mid-Cambrian age. Discovery of this new Silurian taxon allows us to identify a Devonian megacheiran representative, Bundenbachiellus giganteus from the Hunsrück Slate of Germany. The phylogenetic position of megacheirans is controversial: they have been interpreted as stem chelicerates, or stem euarthropods, but when Enalikter and Bundenbachiellus are added to the most comprehensive morphological database available, a stem euarthropod position is supported. Enalikter represents the only fully three-dimensionally preserved stem-group euarthropod, it falls in the sister clade to the crown-group euarthropods, and it provides new insights surrounding the origin and early evolution of the euarthropods. Recognition of Enalikter and Bundenbachiellus as megacheirans indicates that this major arthropod group survived for nearly 100 Myr beyond the mid-Cambrian.     

Action de la chlorhexidine sur l'expression de la virulence de Candida albicans

Candida albicans is an opportunistic yeast. Its pathogenicity is linked to the susceptibility of the host surface as well as to particular factors of the strain: adhesion, filamentous growth and secretion of proteolytic enzymes. Chlorhexidine digluconate is an antiseptic with fungicidal properties. The action of the antiseptic on the growth of the yeast shows a minimal inhibitory concentration (MIC) at 50 μg·mL⁻¹ and a minimal fungicidal concentration (MFC) at 100 μg·mL⁻¹. The consequences of antiseptic treatment are studied using two indicators of pathogenicity: filamentation and the secretion of acid proteinase. Concerning the morphological indicator, a complete inhibition of filamentation for antiseptic concentration far below the MIC is observed; concerning the biochemical indicator, the results show a maintenance of proteolytic activity for concentration equal to the MIC and a disappearance of this activity between the MIC and the MFC. Concerning the proteolytic activity, Skim Milk (Difco) was substituted for bovine serum albumin with success. This mixture of casein allows to use a cheaper and more easily available substrate.     

Bugs in computational chemistry software and their consequences: the importance of the source code

The increase in computer power and the development of new mathematical concepts implemented in software have allowed computational chemistry to emerge as a new research field. Although programs were freely distributed during the "golden age" of this discipline, today they are usually copyrighted and have become easier and easier to use through sophisticated graphical interfaces. This "democratization" is a vector of success for this discipline. Nowadays, non-theoreticians can use such programs more easily and solve chemistry-related problems with the computer. The number of program offerings has rapidly grown and private companies specialized in molecular modeling have appeared and compete to sell their products. Thus, numerous software packages, often presenting similar capabilities, are now available on the market. Within this context, the availability of the program source code remains, in our opinion, an important criterion for program selection.     

Controlled in meso phase crystallization--a method for the structural investigation of membrane proteins

We investigated in meso crystallization of membrane proteins to develop a fast screening technology which combines features of the well established classical vapor diffusion experiment with the batch meso phase crystallization, but without premixing of protein and monoolein. It inherits the advantages of both methods, namely (i) the stabilization of membrane proteins in the meso phase, (ii) the control of hydration level and additive concentration by vapor diffusion. The new technology (iii) significantly simplifies in meso crystallization experiments and allows the use of standard liquid handling robots suitable for 96 well formats. CIMP crystallization furthermore allows (iv) direct monitoring of phase transformation and crystallization events. Bacteriorhodopsin (BR) crystals of high quality and diffraction up to 1.3 Å resolution have been obtained in this approach. CIMP and the developed consumables and protocols have been successfully applied to obtain crystals of sensory rhodopsin II (SRII) from Halobacterium salinarum for the first time.     

Increased genetic gains in sheep,beef and dairy breeding programs from using female reproductive technologies combined with optimal contribution selection and genomic breeding values

Background

Female reproductive technologies such as multiple ovulation and embryo transfer (MOET) and juvenile in vitro embryo production and embryo transfer (JIVET) can boost rates of genetic gain but they can also increase rates of inbreeding. Inbreeding can be managed using the principles of optimal contribution selection (OCS), which maximizes genetic gain while placing a penalty on the rate of inbreeding. We evaluated the potential benefits and synergies that exist between genomic selection (GS) and reproductive technologies under OCS for sheep and cattle breeding programs.

Methods

Various breeding program scenarios were simulated stochastically including: (1) a sheep breeding program for the selection of a single trait that could be measured either early or late in life; (2) a beef breeding program with an early or late trait; and (3) a dairy breeding program with a sex limited trait. OCS was applied using a range of penalties (severe to no penalty) on co-ancestry of selection candidates, with the possibility of using multiple ovulation and embryo transfer (MOET) and/or juvenile in vitro embryo production and embryo transfer (JIVET) for females. Each breeding program was simulated with and without genomic selection.

Results

All breeding programs could be penalized to result in an inbreeding rate of 1 % increase per generation. The addition of MOET to artificial insemination or natural breeding (AI/N), without the use of GS yielded an extra 25 to 60 % genetic gain. The further addition of JIVET did not yield an extra genetic gain. When GS was used, MOET and MOET + JIVET programs increased rates of genetic gain by 38 to 76 % and 51 to 81 % compared to AI/N, respectively.

Conclusions

Large increases in genetic gain were found across species when female reproductive technologies combined with genomic selection were applied and inbreeding was managed, especially for breeding programs that focus on the selection of traits measured late in life or that are sex-limited. Optimal contribution selection was an effective tool to optimally allocate different combinations of reproductive technologies. Applying a range of penalties to co-ancestry of selection candidates allows a comprehensive exploration of the inbreeding vs. genetic gain space.