Binding Sites of theE. ColiDNA Recombinase Protein to the ssDNA: A Computational Study

Modeling the Homologous Recombination Process: Methods, Successes and Challenges

Homologous recombination (HR) is a fundamental process common to all species. HR aims to faithfully repair DNA double strand breaks. HR involves the formation of nucleoprotein filaments on DNA single strands (ssDNA) resected from the break. The nucleoprotein filaments search for homologous regions in the genome and promote strand exchange with the ssDNA homologous region in an unbroken copy of the genome. HR has been the object of intensive studies for decades. Because multi-scale dynamics is a fundamental aspect of this process, studying HR is highly challenging, both experimentally and using computational approaches. Nevertheless, knowledge has built up over the years and has recently progressed at an accelerated pace, borne by increasingly focused investigations using new techniques such as single molecule approaches. Linking this knowledge to the atomic structure of the nucleoprotein filament systems and the succession of unstable, transient intermediate steps that takes place during the HR process remains a challenge; modeling retains a very strong role in bridging the gap between structures that are stable enough to be observed and in exploring transition paths between these structures. However, working on ever-changing long filament systems submitted to kinetic processes is full of pitfalls. This review presents the modeling tools that are used in such studies, their possibilities and limitations, and reviews the advances in the knowledge of the HR process that have been obtained through modeling. Notably, we will emphasize how cooperative behavior in the HR nucleoprotein filament enables modeling to produce reliable information.

Binding mechanism of CDK5 with roscovitine derivatives based on molecular dynamics simulations and MM/PBSA methods

Roscovitine derivatives are potent inhibitors of cyclin-dependent kinase 5 (CDK5), but they exhibit different activities, which has not been understood clearly up to now. On the other hand, the task of drug design is difficult because of the fuzzy binding mechanism. In this context, the methods of molecular docking, molecular dynamics (MD) simulation, and binding free energy analysis are applied to investigate and reveal the detailed binding mechanism of four roscovitine derivatives with CDK5. The electrostatic and van der Waals interactions of the four inhibitors with CDK5 are analyzed and discussed. The calculated binding free energies in terms of MM-PBSA method are consistent with experimental ranking of inhibitor effectiveness for the four inhibitors. The hydrogen bonds of the inhibitors with Cys83 and Lys33 can stabilize the inhibitors in binding sites. The van der Waals interactions, especially the pivotal contacts with Ile10 and Leu133 have larger contributions to the binding free energy and play critical roles in distinguishing the variant bioactivity of four inhibitors. In terms of binding mechanism of the four inhibitors with CDK5 and energy contribution of fragments of each inhibitor, two new CDK5 inhibitors are designed and have stronger inhibitory potency.

Molecular dynamics of zinc-finger ubiquitin binding domains: a comparative study of histone deacetylase 6 and ubiquitin-specific protease 5

HDAC6 is a unique cytoplasmic histone deacetylase characterized by two deacetylase domains, and by a zinc-finger ubiquitin binding domain (ZnF-UBP) able to recognize ubiquitin (Ub). The latter has recently been demonstrated to be involved in the progression of neurodegenerative diseases and in mediating infection by the influenza A virus. Nowadays, understanding the dynamic and energetic features of HDAC6 ZnF-UBP-Ub recognition is considered as a crucial step for the conception of HDAC6 potential modulators. In this study, the atomic, solvent-related, and thermodynamic features behind HDAC6 ZnF-UBP-Ub recognition have been analyzed through molecular dynamics simulations. The behavior was then compared to the prototypical ZnF-UBP from ubiquitin-specific protease 5 (USP5) in order to spot relevant differences useful for selective drug design. Principal component analysis highlighted flapping motions of the L2A loop which were lowered down upon Ub binding in both systems. While polar and nonpolar interactions involving Ub G75 and G76 residues were also common features stabilizing both complexes, salt bridges showed a different pattern, more significant in HDAC6 ZnF-UBP-Ub, whose energetic contribution in USP5 ZnF-UBP-Ub was compensated by the presence of a more stable bridging water molecule. Whereas molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) free energies of binding were comparable for both systems, in agreement with experiments, computational alanine scanning and free energy decomposition data revealed that HDAC6 E1141 and D1178 are potential hotspots for the design of selective HDAC6 modulators.

Modelling and enhanced molecular dynamics to steer structure-based drug discovery

The ever-increasing gap between the availabilities of the genome sequences and the crystal structures of proteins remains one of the significant challenges to the modern drug discovery efforts. The knowledge of structure-dynamics-functionalities of proteins is important in order to understand several key aspects of structure-based drug discovery, such as drug-protein interactions, drug binding and unbinding mechanisms and protein-protein interactions. This review presents a brief overview on the different state of the art computational approaches that are applied for protein structure modelling and molecular dynamics simulations of biological systems. We give an essence of how different enhanced sampling molecular dynamics approaches, together with regular molecular dynamics methods, assist in steering the structure based drug discovery processes.

Potential human cholesterol esterase inhibitor design: benefits from the molecular dynamics simulations and pharmacophore modeling studies

Human pancreatic cholesterol esterase (hCEase) is one of the lipases found to involve in the digestion of large and broad spectrum of substrates including triglycerides, phospholipids, cholesteryl esters, etc. The presence of bile salts is found to be very important for the activation of hCEase. Molecular dynamic simulations were performed for the apoform and bile salt complexed form of hCEase using the co-ordinates of two bile salts from bovine CEase. The stability of the systems throughout the simulation time was checked and two representative structures from the highly populated regions were selected using cluster analysis. These two representative structures were used in pharmacophore model generation. The generated pharmacophore models were validated and used in database screening. The screened hits were refined for their drug-like properties based on Lipinski's rule of five and ADMET properties. The drug-like compounds were further refined by molecular docking simulation using GOLD program based on the GOLD fitness score, mode of binding, and molecular interactions with the active site amino acids. Finally, three hits of novel scaffolds were selected as potential leads to be used in novel and potent hCEase inhibitor design. The stability of binding modes and molecular interactions of these final hits were re-assured by molecular dynamics simulations.

Mechanism of antisense oligonucleotide interaction with natural RNAs

Oligonucleotides find several numbers of applications: as diagnostic probes, RT and PCR primers and antisense agents due to their ability of forming specific interactions with complementary nucleotide sequences within nucleic acids. These interactions are strongly affected by accessibility of the target sequence in the RNA structure. In the present work the mechanism of invasion of RNA structure by oligonucleotide was investigated using a model system: yeast tRNA(Phe) and oligonucleotides complementary to the 3'-part of this molecule. Kinetics of interaction of oligonucleotides with in vitro transcript of yeast tRNAPhe was studied using stopped-flow technique with fluorescence quenching detection, 5'-DABCYL labeled oligonucleotide was hybridized with 3'-fluorescein labeled tRNA(Phe). The results evidence for a four-step invasion process of the oligonucleotide-RNA complex formation. The process is initiated by formation of transition complexes with nucleotides in the T-loop and ACCA sequence. This complex formation is followed by RNA unfolding and formation of an extended heteroduplex with the oligonucleotide via strand displacement process. Computer modeling of oligonucleotide-tRNA(Phe) interaction revealed potential factors that could favor transition complexes formation and confirmed the proposed mechanism, showing the oligonucleotide to be a molecular "wedge". Our data evidence that oligonucleotide invasion into structured RNA is initiated by loop-single strand interactions, similar to the initial step of the antisense RNA-RNA interactions. The obtained results can be used for choosing efficient oligonucleotide probes.

Using computer-aided drug design and medicinal chemistry strategies in the fight against diabetes

The aim of this work is to present a simple, practical and efficient protocol for drug design, in particular Diabetes, which includes selection of the illness, good choice of a target as well as a bioactive ligand and then usage of various computer aided drug design and medicinal chemistry tools to design novel potential drug candidates in different diseases. We have selected the validated target dipeptidyl peptidase IV (DPP-IV), whose inhibition contributes to reduce glucose levels in type 2 diabetes patients. The most active inhibitor with complex X-ray structure reported was initially extracted from the BindingDB database. By using molecular modification strategies widely used in medicinal chemistry, besides current state-of-the-art tools in drug design (including flexible docking, virtual screening, molecular interaction fields, molecular dynamics, ADME and toxicity predictions), we have proposed 4 novel potential DPP-IV inhibitors with drug properties for Diabetes control, which have been supported and validated by all the computational tools used herewith.

Screening from the world's largest TCM database against H1N1 virus

The swine influenza virus (H1N1) 2009 pandemic highlights the importance of having effective anti-viral strategies. Recently, oseltamivir (Tamiflu) resistant influenza viruses are identified; which further emphasizes the urgency in developing new antiviral agents. In influenza virus replication cycle, viral surface glycoprotein, hemagglutinin, is responsible for viral entry into host cells. Hence, a potentially effective antiviral strategy is to inhibit viral entry mechanism. To develop novel antiviral agent that inhibits viral entry, we analyzed 20,000 traditional Chinese medicine (TCM) ingredients in hemagglutinin subtype H1 sialic acid binding site found on H1N1 virus. We then performed molecular dynamics simulations to investigate receptor-ligand interaction of the candidates obtained from docking. Here, we report three TCM derivatives that have high binding affinities to H1 sialic acid binding site residues based on structure-based calculations. The top three derivatives, xylopine_2, rosmaricine_14 and rosmaricine_15, all have an amine group that interact with Glu83 and a pyridinium group that interact with Asp103. Molecular dynamics simulations show that these derivatives form strong hydrogen bonding with Glu83 but interact transiently with Asp103. We therefore suggest that an enhanced hemagglutinin inhibitor, based on our scaffold, should be designed to bind both Glu83 and Asp103 with high affinity.

Thermal stability and unfolding pathways of Sso7d and its mutant F31A: insight from molecular dynamics simulation

The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the protein's stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.