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.     

Biochemical and Structural Insights into RNA Binding by Ssh10b,a Member of the Highly Conserved Sac10b Protein Family in Archaea

Proteins of the Sac10b family are highly conserved in Archaea. Ssh10b, a member of the Sac10b family from the hyperthermophilic crenarchaeon Sulfolobus shibatae, binds to RNA in vivo. Here we show that binding by Ssh10b destabilizes RNA secondary structure. Structural analysis of Ssh10b in complex with a 25-bp RNA duplex containing local distortions reveals that Ssh10b binds the two RNA strands symmetrically as a tetramer with each dimer bound asymmetrically to a single RNA strand. Amino acid residues involved in double-stranded RNA binding are similar, but non-identical, to those in dsDNA binding. The dimer-dimer interaction mediated by the intermolecular β-sheet appears to facilitate the destabilization of base pairing in the secondary structure of RNA. Our results suggest that proteins of the Sac10b family may play important roles in RNA transactions requiring destabilization of RNA secondary structure in Sulfolobus.     

Thermal unfolding of the archaeal DNA and RNA binding protein Ssh10

The reversible thermal unfolding of the archaeal histone-like protein Ssh10b from the extremophile Sulfolobus shibatae was studied using differential scanning calorimetry and circular dichroism spectroscopy. Analytical ultracentrifugation and gel filtration showed that Ssh10b is a stable dimer in the pH range 2.5–7.0. Thermal denaturation data fit into a two-state unfolding model, suggesting that the Ssh10 dimer unfolds as a single cooperative unit with a maximal melting temperature of 99.9 °C and an enthalpy change of 134 kcal/mol at pH 7.0. The heat capacity change upon unfolding determined from linear fits of the temperature dependence of ΔH^cal is 2.55 kcal/(mol K). The low specific heat capacity change of 13 cal/(mol K residue) leads to a considerable flattening of the protein stability curve (ΔG (T)) and results in a maximal ΔG of only 9.5 kcal/mol at 320 K and a ΔG of only 6.0 kcal/mol at the optimal growth temperature of Sulfolobus.     

Ssh10b2与其同源蛋白Ssh10b在细胞含量及固定DNA超螺旋的能力上存在明显差异

极端嗜热古菌———芝田硫化叶菌(Sulfolobus shibatae)基因组含一对亲缘关系较远的同源基因,ssh10b和ssh10b2。这对同源基因编码的蛋白(Ssh10b和Ssh10b2)属于古菌Sac10b DNA结合蛋白家族。关于Ssh10b以及与其极为相似的硫矿硫化叶菌(S.solfataricus)Sso10b、嗜酸热硫化叶菌(S.acidocaldarius)Sac10b蛋白已有较多研究,推测这些蛋白可能在染色体组织和包装、DNA重组、基因表达调控等方面起作用。克隆并在大肠杆菌中表达了ssh10b2基因,纯化了重组Ssh10b2蛋白。免疫印迹定量分析表明,ssh10b2在芝田硫化叶菌中有表达,但其细胞含量仅相当于Ssh10b的约十分之一。重组Ssh10b2对双链DNA的亲和力低于Ssh10b。此外,Ssh10b2和Ssh10b在与双链DNA结合时表现出相似的凝胶阻滞模式。有意思的是,Ssh10b2固定DNA负超螺旋的能力明显低于Ssh10b。这些结果提示,Ssh10b和Ssh10b2可能具有不同的生理作用。     

Understanding the acylation mechanisms of active-site serine penicillin-recognizing proteins: a molecular dynamics simulation study

Beta-lactam antibiotics inhibit enzymes involved in the last step of peptidoglycan synthesis. These enzymes, also identified as penicillin-binding proteins (PBPs), form a long-lived acyl-enzyme complex with beta-lactams. Antibiotic resistance is mainly due to the production of beta-lactamases, which are enzymes that hydrolyze the antibiotics and so prevent them reaching and inactivating their targets, and to mutations of the PBPs that decrease their affinity for the antibiotics. In this study, we present a theoretical study of several penicillin-recognizing proteins complexed with various beta-lactam antibiotics. Hybrid quantum mechanical/molecular mechanical potentials in conjunction with molecular dynamics simulations have been performed to understand the role of several residues, and pK(a) calculations have also been done to determine their protonation state. We analyze the differences between the beta-lactamase TEM-1, the membrane-bound PBP2x of Streptococcus pneumoniae, and the soluble DD-transpeptidase of Streptomyces K15.     

Effect of molecular confinement on internal enzyme dynamics: frequency domain fluorometry and molecular dynamics simulation studies

The tryptophanyl emission decay of the mesophilic beta-galactosidase from Aspergillus oryzae free in buffer and entrapped in agarose gel is investigated as a function of temperature and compared to that of the hyperthermophilic enzyme from Sulfolobus solfataricus. Both enzymes are tetrameric proteins with a large number of tryptophanyl residues, so the fluorescence emission can provide information on the conformational dynamics of the overall protein structure rather than that of the local environment. The tryptophanyl emission decays are best fitted by bimodal Lorentzian distributions. The long-lived component is ascribed to close, deeply buried tryptophanyl residues with reduced mobility; the short-lived one arises from tryptophanyl residues located in more flexible external regions of each subunit, some of which are involved in forming the catalytic site. The center of both lifetime distribution components at each temperature increases when going from the free in solution mesophilic enzyme to the gel-entrapped and hyperthermophilic enzyme, thus indicating that confinement of the mesophilic enzyme in the agarose gel limits the freedom of the polypeptide chain. A more complex dependence is observed for the distribution widths. Computer modeling techniques are used to recognize that the catalytic sites are similar for the mesophilic and hyperthermophilic beta-galactosidases. The effect due to gel entrapment is considered in dynamic simulations by imposing harmonic restraints to solvent-exposed atoms of the protein with the exclusion of those around the active site. The temperature dependence of the tryptophanyl fluorescence emission decay and the dynamic simulation confirm that more rigid structures, as in the case of the immobilized and/or hyperthermophilic enzyme, require higher temperatures to achieve the requisite conformational dynamics for an effective catalytic action and strongly suggest a link between conformational rigidity and enhanced thermal stability.     

Proline in alpha-helix: stability and conformation studied by dynamics simulation

Free-energy simulations have been used to estimate the change in the conformational stability of short polyalanine alpha-helices when one of the alanines is replaced by a proline residue. For substituting proline in the middle of the helix the change in free energy of folding (delta delta G degrees) was calculated as 14 kJ/mol (3.4 kcal/mol), in excellent agreement with the one available experimental value. The helix containing proline was found to be strongly kinked; the free energy for reducing the angle of the kink from 40 degrees to 15 degrees was calculated, and found to be small. A tendency to alternate hydrogen bonding schemes was observed in the proline-containing helix. These observations for the oligopeptide agree well with the observation of a range of kink angles (18-35 degrees) and variety of hydrogen bonding schemes, in the rare instances where proline occurs in helices in globular proteins. For substituting proline at the N-terminus of the helix the change in free energy of folding (delta delta G degrees) was calculated as -4 kJ/mol in the first helical position (N1) and +6 kJ/mol in the second helical position (N2). The observed frequent occurrence of proline in position N1 in alpha-helices in proteins therefore has its origin in stability differences of secondary structure. The conclusion reached here that proline may be a better helix former in position N1 than (even) alanine, and thus be a helix initiator may be testable experimentally by measurements of fraction helical conformation of individual residues in oligopeptides of appropriate sequence. The relevance of these results in regards to the frequent occurrence of proline-containing helices in certain membrane proteins is discussed.     

Time-resolved fluorescence spectroscopy and molecular dynamics simulations point out the effects of pressure on the stability and dynamics of the porcine odorant-binding protein

The effects of hydrostatic pressure on the structure and stability of porcine odorant-binding protein (pOBP) in the presence and absence of the odorant molecule 2-isobutyl-3-methoxypyrazine (IBMP) were studied by steady-state and time-resolved fluorescence spectroscopy as well as by molecular dynamics simulation. The authors found that the application of moderate values of hydrostatic pressure to pOBP solutions perturbed the microenvironment of Trp(16) and disrupted its highly quenched complex with Met(39). In addition, compared with the protein in the absence of IBMP, the MD simulations experiments carried out at different pressures highlighted the role of this ligand in stabilizing the Trp(16)/Met(39) interaction even at 2000 bar. The obtained results will assist for the tailoring of this protein as specific sensing element in a new class of fluorescence-based biosensors for the detection of explosives.     

Exploring the effect of inhibitor AKB-9778 on VE-PTP by molecular docking and molecular dynamics simulation

Diabetic macular edema, also known as diabetic eye disease, is mainly caused by the overexpression of vascular endothelial protein tyrosine phosphatase (VE-PTP) at hypoxia/ischemic. AKB-9778 is a known VE-PTP inhibitor that can effectively interact with the active site of VE-PTP to inhibit the activity of VE-PTP. However, the binding pattern of VE-PTP with AKB-9778 and the dynamic implications of AKB-9778 on VE-PTP system at the molecular level are poorly understood. Through molecular docking, it was found that the AKB-9778 was docked well in the binding pocket of VE-PTP by the interactions of hydrogen bond and Van der Waals. Furthermore, after molecular dynamic simulations on VE-PTP system and VE-PTP ^AKB-9778 system, a series of postdynamic analyses found that the flexibility and conformation of the active site undergone an obvious transition after VE-PTP binding with AKB-9778. Moreover, by constructing the RIN, it was found that the different interactions in the active site were the detailed reasons for the conformational differences between these two systems. Thus, the finding here might provide a deeper understanding of AKB-9778 as VE-PTP Inhibitor.     

Exploring the mechanism how Marburg virus VP35 recognizes and binds dsRNA by molecular dynamics simulations and free energy calculations

Filoviruses often cause terrible infectious disease which has not been successfully dealt with pharmacologically. All filoviruses encode a unique protein termed VP35 which can mask doubled‐stranded RNA to deactivate interferon. The interface of VP35–dsRNA would be a feasible target for structure‐based antiviral agent design. To explore the essence of VP35–dsRNA interaction, molecular dynamics simulation combined with MM‐GBSA calculations were performed on Marburg virus VP35–dsRNA complex and several mutational complexes. The energetic analysis indicates that nonpolar interactions provide the main driving force for the binding process. Although the intermolecular electrostatic interactions play important roles in VP35–dsRNA interaction, the whole polar interactions are unfavorable for binding which result in a low binding affinity. Compared with wild type VP35, the studied mutants F228A, R271A, and K298A have obviously reduced binding free energies with dsRNA reflecting in the reduction of polar or nonpolar interactions. The results also indicate that the loss of binding affinity for one dsRNA strand would abolish the total binding affinity. Three important residues Arg271, Arg294, and Lys298 which makes the largest contribution for binding in VP35 lose their binding affinity significantly in mutants. The uncovering of VP35–dsRNA recognition mechanism will provide some insights for development of antiviral drug. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 849–860, 2014.     

Relationship between conformation and geometry as evidenced by molecular dynamics simulation of Cα,α-dialkylated glycines

The relationship between the local backbone conformation and bond angles at C^α of symmetrically substituted C^α,α-dialkylated glycines (C^α,α-dimethylglycine or α-aminoisobutyric acid, Aib; C^α,α-diethylglycine, Deg; C^α,α-di-n-propylglycine, Dpg) has been investigated by molecular dynamics (MD) simulation adopting flat bottom harmonic potentials, instead of the usual harmonic restraints, for the C^α bond angles. The MD simulations show that the C^α bond angles are related to the local backbone conformation, irrespectively of the side-chain length of Aib, Deg, and Dpg residues. Moreover, the N-C^α-C′ (τ) angle is the most sensitive conformational parameter and, in the folded form, is always larger and more flexible than in the extended one. © 1998 John Wiley & Sons, Inc. Biopoly 46: 239–244, 1998     

Insight into the intrinsic flexibility of the SL1 stem-loop from genomic RNA of HIV-1 as probed by molecular dynamics simulation

The SL1 stem-loop is the dimerization initiation site for linking the two copies of the RNA forming the HIV-1 genome. The 26 nucleotides stem contains a defect consisting on a highly conserved G-rich 1-3 asymmetrical internal loop, which is a major site for nucleocapsid protein binding. Several NMR attempts were undertaken to determine the internal loop structure in the SL1 monomer. However, the RNA constructs used in the different studies were largely mutated, in particular with replacement of the nine nucleotides apical loop by a tetraloop, and divergent results were obtained ranging from a rigid structure to a particularly large flexibility. To investigate the reasons for such discrepancies, we used molecular dynamics simulation of the SL1 monomer to probe the effect of mutations on the conformational stability of the internal loop and of the whole stem. It is found that in the wild-type sequence, the internal loop displays conformational variability originating mainly from the nine nucleotides apical loop flexibility that causes large conformational fluctuations (without changing the average structure) in the 7 bp duplex linking the apical and internal loops. The large amplitude atomic motions in the duplex are transmitted to the internal loop in which they induce conformational changes characterized by a labile hydrogen bond network such as G5 successively H-bonded to A29 and G30. On the contrary, with a four nucleotides apical loop, conformational fluctuations in the duplex are reduced by a factor of 2 and are not sufficiently energizing for promoting changes in the internal loop structure at the time scale of the simulations.     

Geometrical constraints limiting the poly(ADP‐ribose) conformation investigated by molecular dynamics simulation

Poly(ADP‐ribosylation) is a post‐transductional modification that regulates protein's function. Most of the proteins subjected to this control mechanism belong to machineries involved in DNA damage repair, or DNA interacting proteins. Poly(ADP‐ribose) polymers are long chains of even 100 monomer length that can be branched at several positions but, not withstanding its importance, nothing is known concerning its structure. To understand, which are the geometrical parameters that confer to the polymer the structural constraints that determine its interaction with the target proteins, we have performed molecular dynamics of three chains of different length, made by 5, 25, and 30 units, the last one being branched. Analysis of the simulations allowed us to identify the main intra‐ and inter‐monomer dihedral angles that govern the structure of the polymer that however, does not reach a unique definite conformation. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 78–86, 2014.     

Computational investigation of the HIV-1 Rev multimerization using molecular dynamics simulations and binding free energy calculations

The HIV Rev protein mediates the nuclear export of viral mRNA, and is thereby essential for the production of late viral proteins in the replication cycle. Rev forms a large organized multimeric protein-protein complex for proper functioning. Recently, the three-dimensional structures of a Rev dimer and tetramer have been resolved and provide the basis for a thorough structural analysis of the binding interaction. Here, molecular dynamics (MD) and binding free energy calculations were performed to elucidate the forces thriving dimerization and higher order multimerization of the Rev protein. It is found that despite the structural differences between each crystal structure, both display a similar behavior according to our calculations. Our analysis based on a molecular mechanics-generalized Born surface area (MM/GBSA) and a configurational entropy approach demonstrates that the higher order multimerization site is much weaker than the dimerization site. In addition, a quantitative hot spot analysis combined with a mutational analysis reveals the most contributing amino acid residues for protein interactions in agreement with experimental results. Additional residues were found in each interface, which are important for the protein interaction. The investigation of the thermodynamics of the Rev multimerization interactions performed here could be a further step in the development of novel antiretrovirals using structure based drug design. Moreover, the variability of the angle between each Rev monomer as measured during the MD simulations suggests a role of the Rev protein in allowing flexibility of the arginine rich domain (ARM) to accommodate RNA binding.     

Evolutionary conservation and DNA binding properties of the Ssh7 proteins fromSulfolobus shibatae

The thermoacidophilic archaeonSulfolobus shibatae synthesizes a large amount of the 7-ku DNA binding proteins known as Ssh7. Our hybridization experiments showed that two Ssh7-encoding genes existed in the genome of S.shibatae. These two genes, designatedssh7a andssh7b, have been cloned, sequenced and expressed inEscherichia coli. The two Ssh7 proteins differ only at three amino acid positions. In addition, thecis-regulatory sequences of thessh7a andssh7b genes are highly conserved. These results suggest the presence of a selective pressure to maintain not only the sequence but also the expression of the two genes. We have also found that there are two genes encoding the 7-ku protein inSulfolobus solfataricus. Based on this and other studies, we suggest that the gene encoding the 7-ku protein underwent duplication before the separation ofSulfolobus species. Binding of native Ssh7 and recombinant (r)Ssh7 to short duplex DNA fragments was analyzed by electrophoretic mobility shift assays. Both native and recombinant forms of the protein behaved in a similar fashion in the assays, suggesting that the interaction of Ssh7 with DNA is not affected either by specific lysine methylation found in the native Ssh7 proteins or by the difference between the two Ssh7 isomers in amino acid sequence. Our data show that Ssh7 binds duplex DNA fragments with a binding size of ∼ 6.6 base pairs and an apparent dissociation constant of (0.7–1.0) × 10^-7 mol/L under the assay conditions employed in the present study. In addition, Ssh7 binds more tightly to negatively supercoiled DNA than to linear or relaxed DNA.