MicroRNAs (miRNAs) are endogenously produced ∼21-nt riboregulators that associate with Argonaute (Ago) proteins to direct mRNA cleavage or repress the translation of complementary RNAs. Capturing the molecular mechanisms of miRNA interacting with its target will not only reinforce the understanding of underlying RNA interference but also fuel the design of more effective small-interfering RNA strands. To address this, in the present work the RNA-bound (Ago-miRNA, Ago-miRNA-target) and RNA-free Ago forms were analyzed by performing both molecular dynamics simulations and thermodynamic analysis. Based on the principal component analysis results of the simulation trajectories as well as the correlation analysis in fluctuations of residues, we discover that: 1) three important (PAZ, Mid and PIWI) domains exist in Argonaute which define the global dynamics of the protein; 2) the interdomain correlated movements are so crucial for the interaction of Ago-RNAs that they not only facilitate the relaxation of the interactions between residues surrounding the RNA binding channel but also induce certain conformational changes; and 3) it is just these conformational changes that expand the cavity of the active site and open putative pathways for both the substrate uptake and product release. In addition, by thermodynamic analysis we also discover that for both the guide RNA 5′-end recognition and the facilitated site-specific cleavage of the target, the presence of two metal ions (of Mg²⁺) plays a predominant role, and this conclusion is consistent with the observed enzyme catalytic cleavage activity in the ternary complex (Ago-miRNA-mRNA). Our results find that it is the set of arginine amino acids concentrated in the nucleotide-binding channel in Ago, instead of the conventionally-deemed seed base-paring, that makes greater contributions in stabilizing the binding of the nucleic acids to Ago. One of the biggest surprises at the beginning of the ‘post-genome era’ was the discovery of numerous genes encoding microRNAs. The number of microRNA genes is estimated to be nearly 1% of that of protein-coding genes, which were found in genomes of such diverse organisms as Caenorhabditis elegans, Drosophila melanogaster, Arabidopsis thaliana, and Homo sapiens. Their products, tiny RNAs (miRNAs and siRNAs), are thought to bind to Argonaute (Ago) proteins and form effector complexes to direct mRNA cleavage or repress translation of complementary RNAs, during development, organogenesis, and very likely during many other processes. The cellular interactions between the miRNAs and their target RNAs associating with Ago are only beginning to be revealed, and details of this interaction mechanism at molecular level are still poorly understood. In this article we propose the possible mechanisms of miRNA-target interaction with special emphasis on their structural dynamic and thermodynamic aspects. The results of our model suggest the chemical and physical factors and effects that may be responsible for the miRNA-Ago assembly and miRNA-target recognition.

Dityrosine is the product of oxidation that has been linked to a number of serious pathological conditions. Evidence indicates that high amounts of dityrosine exist in oxidized milk powders and some milk related foodstuffs, further reducing the nutritional value of oxidized proteins. Therefore, we hypothesize that some receptors related to special diseases would be targets for dityrosine. However, the mechanisms of the interaction of dityrosine with probable targets are still unknown. In the present work, an inverse virtual screening approach was performed to screen possible novel targets for dityrosine. Molecular docking studies were performed on a panel of targets extracted from the potential drug target database (PDTD) to optimize and validate the screening results. Firstly, two different conformations cis- and trans- were found for dityrosine during minimization. Moreover, Tubulin (αT) (−11.0 kcal/mol) was identified as a target for cis-dityrosine (CDT), targets including αT (−11.2 kcal/mol) and thyroid hormone receptor beta-1 (−10.7 kcal/mol) presented high binding affinities for trans-dityrosine (TDT). Furthermore, in order to provide binding complexes with higher precision, the three docked systems were further refined by performing thermo dynamic simulations. A series of techniques for searching for the most stable binding pose and the calculation of binding free energy are elaborately provided in this work. The major interactions between these targets and dityrosine were hydrophobic, electrostatic and hydrogen bonding. The application of inverse virtual screening method may facilitate the prediction of unknown targets for known ligands, and direct future experimental assays.

Apical membrane antigen 1 (AMA1) interacts with RON2 to form a protein complex that plays a key role in the invasion of host cells by malaria parasites. Blocking this protein–protein interaction represents a potential route to controlling malaria and related parasitic diseases, but the polymorphic nature of AMA1 has proven to be a major challenge to vaccine-induced antibodies and peptide inhibitors exerting strain-transcending inhibitory effects. Here we present the X-ray crystal structure of AMA1 domains I and II from Plasmodium falciparum strain FVO. We compare our new structure to those of AMA1 from P. falciparum 3D7 and Plasmodium vivax. A combination of normalized B factor analysis and computational methods has been used to investigate the flexibility of the domain I loops and how this correlates with their roles in determining the strain specificity of human antibody responses and inhibitory peptides. We also investigated the domain II loop, a key region involved in inhibitor binding, by comparison of multiple AMA1 crystal structures. Collectively, these results provide valuable insights that should contribute to the design of strain-transcending agents targeting P. falciparum AMA1.

In order to explore the structure-activity correlation of benzothiadiazine series as inhibitors of genotype 1a HCV polymerase, a set of ligand- and receptor-based 3D-QSAR models were, for the first time, developed in the present work employing Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Index Analysis (CoMSIA) for 239 promising molecules. In addition, homology modeling, docking analysis, and molecular dynamics simulation (MD) were also applied to elucidate the probable binding modes of these inhibitors at the allosteric site of the enzyme. The statistical model validations assure the reliability of the obtained QSAR models. Changes in the binding affinity of the inhibitors attributing to modifications in the aromatic rings could be rationalized by the steric, electrostatic, hydrophobic, and hydrogen bond acceptor properties. (i) Hydrophobic substituents with similar size of benzo group like isosteres are preferential at positions 1 and 2 (ring B of benzothiadiazines). (ii) Substituents at position-3 containing a linear alkyl chain (four or five carbon atoms) or a branched alkyl chain (five-eight carbons) can increase the inhibitory activity by one to two orders of magnitude. (iii) A polar substituent like methanesulfonamide group at position-14 can enhance the activity of the drug by providing a hydrogen bonding interaction with the protein target. The results obtained from this work provide important guidelines in design of novel benzothiadiazine analogs as inhibitors of HCV genotype 1a NS5B.