Molecular dynamics simulation studies of induced fit and conformational capture in U1A–RNA binding: Do molecular substates code for specificity?

Molecular dynamics (MD) simulations on stem loop 2 of U1 small nuclear RNA and a construct of the U1A protein were carried out to obtain predictions of the structures for the unbound forms in solution and to elucidate dynamical aspects of induced fit upon binding. A crystal structure of the complex between the U1A protein and stem loop 2 RNA and an NMR structure for the uncomplexed form of the U1A protein are available from Oubridge et al. (Nature, 1994, Vol. 372, pp. 432-438) and Avis et al. (Journal of Molecular Biology, 1996, Vol. 257, pp. 398-411), respectively. As a consequence, U1A-RNA binding is a particularly attractive case for investigations of induced fit in protein-nucleic acid complexation. When combined with the available structural data, the results from simulations indicate that structural adaptation of U1A protein and RNA define distinct mechanisms for induced fit. For the protein, the calculations indicate that induced fit upon binding involves a non-native thermodynamic substate in which the structure is preorganized for binding. In contrast, induced fit of the RNA involves a distortion of the native structure in solution to an unstable form. However, the RNA solution structures predicted from simulation show evidence that structures in which groups of bases are favorably oriented for binding the U1A protein are thermally accessible. These results, which quantify with computational modeling recent proposals on induced fit and conformational capture by Leuillot and Varani (Biochemistry, 2001, Vol. 40, pp. 7947-7956) and by Williamson (Nature Structural Biology, 2000, Vol. 7, pp. 834-837) suggest an important role for intrinsic molecular architecture and substates other than the native form in the specificity of protein-RNA interactions.     

Molecular dynamics simulation studies for DNA sequence recognition by reactive metabolites of anticancer compounds

The discovery of novel anticancer molecules 5F‐203 (NSC703786) and 5‐aminoflavone (5‐AMF, NSC686288) has addressed the issues of toxicity and reduced efficacy by targeting over expressed Cytochrome P450 1A1 (CYP1A1) in cancer cells. CYP1A1 metabolizes these compounds into their reactive metabolites, which are proven to mediate their anticancer effect through DNA adduct formation. However, the drug metabolite–DNA binding has not been explored so far. Hence, understanding the binding characteristics and molecular recognition for drug metabolites with DNA is of practical and fundamental interest. The present study is aimed to model binding preference shown by reactive metabolites of 5F‐203 and 5‐AMF with DNA in forming DNA adducts. To perform this, three different DNA crystal structures covering sequence diversity were selected, and 12 DNA‐reactive metabolite complexes were generated. Molecular dynamics simulations for all complexes were performed using AMBER 11 software after development of protocol for DNA‐reactive metabolite system. Furthermore, the MM‐PBSA/GBSA energy calculation, per‐nucleotide energy decomposition, and Molecular Electrostatic Surface Potential analysis were performed. The results obtained from present study clearly indicate that minor groove in DNA is preferable for binding of reactive metabolites of anticancer compounds. The binding preferences shown by reactive metabolites were also governed by specific nucleotide sequence and distribution of electrostatic charges in major and minor groove of DNA structure. Overall, our study provides useful insights into the initial step of mechanism of reactive metabolite binding to the DNA and the guidelines for designing of sequence specific DNA interacting anticancer agents. Copyright © 2014 John Wiley & Sons, Ltd.     

Molecular dynamics simulation of interleukin-2 and its complex and determination of the binding free energy

Interleukin-2 (IL-2) protein belongs to the signal modulator cytokine's family and therefore it is prevalent for immunological responses. It has been identified as a centrally important potential drug target for the inhibition of protein-protein interactions; so as to suppress the immunological responses associated with autoimmune, inflammatory and immunological diseases, and cancer. In the present work, we have performed two independent 100?ns of molecular dynamics (MD) simulations on the apo IL-2 protein and its ligand-bound complex (with a potent inhibitor FRG), to study the effect of inhibitor binding on the dynamics and stability of the protein. The calculation of binding free energy via post-processing end state method of Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA) and Molecular Mechanics Generalised Born Surface Area (MM-GBSA) has inferred a good correlation in accordance with the already reported experimental data, demonstrating that the free energy of binding calculated by the two methods has no significant difference. The investigation of individual components of free energy revealed that the association of IL-2 protein with FRG ligand is primarily driven by the van der Waals energy contribution that represents the non-polar/hydrophobic energy contribution as dominant in this case of ligand binding.     

Molecular dynamics,thermodynamic, and mutational binding studies for tumor-specific LyP-1 in complex with p32

Recent studies in tumor homing peptides have shown the specificity of LyP-1 (CGNKRTRGC) to tumor lymphatics. In this present work, we evaluated the possible interactions between cyclic LyP-1 and its receptor, p32, with molecular dynamics and docking studies in order to lead the design of novel LyP-1 derivatives, which could bind to p32 more effectively and perform enhanced antitumor effect. The total binding enthalpy energies have been obtained by MM-PBSA thermodynamic computations and the favorability of p32.LyP-1 complex in water has been shown by explicit water MD computations. The last 30 ns of molecular dynamics trajectory have shown the strong interaction of LyP-1 with the inner surface chains of p32, especially with chains B and C. ALA-SCAN mutagenesis studies have indicated the considerable influence of Asn3, Lys4, Arg5, and Arg7 amino acid residues on the specific binding of LyP-1. Within the knowledge of the critical role of p32 receptor in cancer cell metabolism, this study can lead to further developments in anticancer therapy by targeting p32 with LyP-1 derivatives as active targeting moiety. This data can also be applied for the development of new drug delivery systems in which LyP-1 can be used for its targeting and anticancer properties.     

Molecular dynamics simulation of the induced‐fit binding process of DNA aptamer and L‐argininamide

Aptamers are rare functional nucleic acids with binding affinity to and specificity for target ligands. Recent experiments have lead to the proposal of an induced‐fit binding mechanism for L ‐argininamide (Arm) and its binding aptamer. However, at the molecular level, this mechanism between the aptamer and its coupled ligand is still poorly understood. The present study used explicit solvent molecular dynamics (MD) simulations to examine the critical bases involved in aptamer‐Arm binding and the induced‐fit binding process at atomic resolution. The simulation results revealed that the Watson‐Crick pair (G10‐C16), C9, A12, and C17 bases play important roles in aptamer‐Arm binding, and that binding of Arm results in an aptamer conformation optimized through a general induced‐fit process. In an aqueous solution, the mechanism has the following characteristic stages: (a) adsorption stage, the Arm anchors to the binding site of aptamer with strong electrostatic interaction; (b) binding stage, the Arm fits into the binding site of aptamer by hydrogen‐bond formation; and (c) complex stabilization stage, the hydrogen bonding and electrostatic interactions cooperatively stabilize the complex structure. This study provides dynamics information on the aptamer‐ligand induced‐fit binding mechanism. The critical bases in aptamer‐ligand binding may provide a guideline in aptamer design for molecular recognition engineering.     

Targeted molecular dynamics simulation studies of binding and conformational changes in E. coli MurD

Enzymes involved in the biosynthesis of bacterial peptidoglycan, an essential cell wall polymer unique to prokaryotic cells, represent a highly interesting target for antibacterial drug design. Structural studies of E. coli MurD, a three-domain ATP hydrolysis driven muramyl ligase revealed two inactive open conformations of the enzyme with a distinct C-terminal domain position. It was hypothesized that the rigid body rotation of this domain brings the enzyme to its closed active conformation, a structure, which was also determined experimentally. Targeted molecular dynamics 1 ns-length simulations were performed in order to examine the substrate binding process and gain insight into structural changes in the enzyme that occur during the conformational transitions into the active conformation. The key interactions essential for the conformational transitions and substrate binding were identified. The results of such studies provide an important step toward more powerful exploitation of experimental protein structures in structure-based inhibitor design.     

Molecular dynamics simulation of intrinsically disordered proteins

Intrinsically disordered proteins are biomolecules that do not have a definite 3D structure; therefore, their dynamical simulation cannot start from a known list of atomistic positions, such as a Protein Data Bank file. We describe a method to start a computer simulation of these proteins. The first step of the procedure is the creation of a multi-rod configuration of the molecule, derived from its primary sequence. This structure is dynamically evolved in vacuo until its gyration radius reaches the experimental average value; at this point solvent molecules, in explicit or implicit implementation, are added to the protein and a regular molecular dynamics simulation follows. We have applied this procedure to the simulation of tau, one of the largest totally disordered proteins.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

Molecular dynamics simulation of truncated bovine adrenodoxin

The 128 amino acid long soluble protein adrenodoxin (Adx) is a typical member of the ferredoxin protein family that are electron carrier proteins with an iron-sulfur cofactor. Adx carries electrons from adrenodoxin reductase (AdR) to cytochrome P450s. Its binding modes to these proteins were previously characterized by site-directed mutagenesis, by X-ray crystallography for the complex Adx:AdR, and by NMR. However, no clear evidence has been provided for the driving force that promotes Adx detachment from AdR upon reduction. Here, we characterized the conformational dynamics of unbound Adx in the oxidized and reduced forms using 2-20 ns long molecular dynamics simulations. The most noticeable difference between both forms is the enhanced flexibility of the loop (47-51) surrounding the iron-sulfur cluster in the reduced form. Together with several structural displacements at the binding interface, this increased flexibility may be the key factor promoting unbinding of reduced Adx from AdR. This points to an intrinsic property of reduced Adx that drives dissociation.     

Trax is a component of the Translin-containing RNA binding complex

Translin is a nucleic acid binding protein that has been implicated in regulating the targeting and translation of dendritic RNA. In previous studies, we found that Translin and its partner protein, Trax, are components of a gel-shift complex that is highly enriched in brain extracts. In those studies, we employed a DNA oligonucleotide, GS1, as a probe to label the complex. Translin has also been identified as a component of a gel-shift complex detected using an RNA oligonucleotide probe, derived from the 3' UTR of protamine-2 mRNA. Although we had assumed that these probes labeled the same complex, recent studies indicate that association of Trax with Translin suppresses its RNA binding activity. As these findings challenge this assumption and suggest that the native RNA binding complex does not contain Trax, we have re-examined this issue. We have found that the gel-shift complexes labeled with either GS1 or protamine-2 probes are "supershifted" by addition of Trax antibodies, indicating that both are heteromeric Translin/Trax complexes. In addition, cross-competition studies provide additional evidence that these probes label the same complex. Furthermore, analysis of recombinant Translin/Trax complexes generated by co-transfection of Trax with Translin in hEK293T demonstrates that they are labeled with either probe. Although recombinant Translin forms a homomeric nucleic acid binding complex in vitro, our findings indicate that both Trax and Translin are components of the native gel-shift complex labeled with either GS1 or protamine-2 probes.     

Molecular simulation of N-acetylneuraminic acid analogs and molecular dynamics studies of cholera toxin-Neu5Gc complex

Cholera toxin (CT) is an AB₅ protein complex secreted by the pathogen Vibrio cholera, which is responsible for cholera infection. N-acetylneuraminic acid (NeuNAc) is a derivative of neuraminic acid with nine-carbon backbone. NeuNAc is distributed on the cell surface mainly located in the terminal components of glycoconjugates, and also plays an important role in cell–cell interaction. In our current study, molecular docking and molecular dynamic (MD) simulations were implemented to identify the potent NeuNAc analogs with high-inhibitory activity against CT protein. Thirty-four NeuNAc analogs, modified in different positions C-1/C-2/C-4/C-5/C-7/C-8/C-9, were modeled and docked against the active site of CT protein. Among the 34 NeuNAc analogs, the analog Neu5Gc shows the least extra precision glide score of ?9.52 and glide energy of ?44.71?kcal/mol. NeuNAc analogs block the CT active site residues HIS:13, ASN:90, LYS:91, GLN:56, GLN:61, and TRP:88 through intermolecular hydrogen bonding. The MD simulation for CT-Neu5Gc docking complex was performed using Desmond. MD simulation of CT-Neu5Gc complex reveals the stable nature of docking interaction.     

Molecular dynamics simulations of apo and holo forms of fatty acid binding protein 5 and cellular retinoic acid binding protein II reveal highly mobile protein,retinoic acid ligand,and water molecules

Structural and dynamic properties from a series of 300 ns molecular dynamics, MD, simulations of two intracellular lipid binding proteins, iLBPs, (Fatty Acid Binding Protein 5, FABP5, and Cellular Retinoic Acid Binding Protein II, CRABP-II) in both the apo form and when bound with retinoic acid reveal a high degree of protein and ligand flexibility. The ratio of FABP5 to CRABP-II in a cell may determine whether it undergoes natural apoptosis or unrestricted cell growth in the presence of retinoic acid. As a result, FABP5 is a promising target for cancer therapy. The MD simulations presented here reveal distinct differences in the two proteins and provide insight into the binding mechanism. CRABP-II is a much larger, more flexible protein that closes upon ligand binding, where FABP5 transitions to an open state in the holo form. The traditional understanding obtained from crystal structures of the gap between two β-sheets of the β-barrel common to iLBPs and the α-helix cap that forms the portal to the binding pocket is insufficient for describing protein conformation (open vs. closed) or ligand entry and exit. When the high degree of mobility between multiple conformations of both the ligand and protein are examined via MD simulation, a new mode of ligand motion that improves understanding of binding dynamics is revealed.     

Exploring the selectivity of a ligand complex with CDK2/CDK1: a molecular dynamics simulation approach

Cyclin‐dependent kinases (CDKs) are core components of the cell cycle machinery that govern the transition between phases during cell cycle progression. Abnormalities in CDKs activity and regulation are common features of cancer, making CDK family members attractive targets for the development of anticancer drugs. Their inhibitors have entered in clinical trials to treat cancer. Very recently, Heathcote et al. (J. Med. Chem. 2010, 53:8508–8522) have found a ligand BS194 that has a high affinity with CDK2 (IC₅₀ = 3 nm ) but shows low affinity with CDK1 (IC₅₀ = 30 nm ). To understand the selectivity, we used homology modeling, molecular docking, molecular dynamics, and free‐energy calculation to analyze the interactions. A rational three‐dimensional model of the CDK1/BS194 complex is built. We found that Leu83 is a key residue that recognizes BS194 more effectively with CDK2 with good binding free energies rather than CDK1. Energetic analysis reveals that van der Waals interaction and non‐polar contributions to solvent are favorable in the formation of complexes and amine group of the ligand, which plays a crucial role for binding selectivity between CDK2 and CDK1. Copyright © 2012 John Wiley & Sons, Ltd.     

Targeted molecular dynamics simulation studies of calcium binding and conformational change in the C-terminal half of gelsolin

Gelsolin consists of six related domains (G1-G6) and the C-terminal half (G4-G6) acts as a calcium sensor during the activation of the whole molecule, a process that involves large domain movements. In this study, we used targeted molecular dynamics simulations to elucidate the conformational transitions of G4-G6 at an atomic level. Domains G4 and G6 are initially ruptured, followed by a rotation of G6 by approximately 90 degrees , which is the dominant conformational change. During this period, local conformational changes occur at the G4 and G5 calcium-binding sites, facilitating large changes in interdomain distances. Alterations in the binding affinities of the calcium ions in these three domains appear to be related to local conformational changes at their binding sites. Analysis of the relative stabilities of the G4-G6-bound calcium ions suggests that they bind first to G6, then to G4, and finally to G5.     

Structural studies of SNARE complex and its interaction with complexin by molecular dynamics simulation

Since neurotransmitter releasing into the synaptic space delivers electrical signals from presynaptic neural cell to the postsynaptic cell, neurotransmitter secretion must be much orchestrated. Crowded intracellular vesicles involving neurotransmitters present a question of the how secretory vesicles fuse onto the plasma membrane in a fast synchronized fashion. Complexin is one of the most experimentally studied proteins that regulate assembly of fusogenic four‐helix SNARE complex to synchronized neurotransmitter secretion. We used MD simulation to investigate the interaction of complexin with the neural SNARE complex in detail. Our results show that the SNARE complex interacts with the complexin central helix by forming salt bridges and hydrogen bonds. Complexin also can interact with the Q‐SNARE complex instead of synaptobrevin to decrease the Q‐SNARE flexibility. The complexin alpha‐accessory helix and the C‐terminal region of synaptobrevin can interact with the same region of syntaxin. Although the alpha‐accessory helix aids the tight binding of the central helix to the SNARE complex, its proximity with synaptobrevin causes the destabilization of syntaxin and Sn1 helices. This study suggests that the alpha‐accessory helix of complexin can be an inhibiting factor for membrane fusion by competing with synaptobrevin for binding to the Q‐SNARE complex. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 560–570, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

RBPMS不同剪接体的真核表达和亚细胞定位

目的：构建具有多种剪接形式的RNA结合蛋白（RBPMS）基因的真核表达载体,并在真核细胞中表达,确定不同形式的RBPMS在细胞中的定位。方法：采用PCR技术从人卵巢cDNA文库中扩增RBPMS基因的几种完整编码序列（命名为RBPl～RBP4）,克隆到带绿色荧光蛋白标签的pEGFP-C1表达载体上,转染人胚肾细胞293T,Western印迹鉴定RBPMS的表达,并利用激光共聚焦显微镜观察RBPMS不同剪接体在细胞中的定位。结果：限制性内切酶分析和DNA序列测定表明构建的重组表达载体正确,Western印迹实验证明RBP1～RBP4表达成功。通过激光共聚焦显微镜观察,RBP1／4围绕胞核在核膜的周围呈聚集状分布;RBP2则在细胞质和细胞核中均有分布,但会出现斑点状聚集;RBP3呈半月状紧密分布在细胞核周围;RBPMS中的RNA识别基序缺失后,这种现象消失,与空载体对照类似,在细胞核和细胞质中均有分布。结论：构建并表达了RBPMS基因的真核表达载体,RBPMS不同剪接体及RNA识别基序缺失后具有不同的亚细胞分布模式,提示具有不同的功能。