Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

Structure and dynamics of the first archaeal parvulin reveal a new functionally important loop in parvulin-type prolyl isomerases

Parvulins are a group of peptidyl-prolyl isomerases (PPIases) responsible for important biological processes in all kingdoms of life. The PinA protein from the psychrophilic archaeon Cenarchaeum symbiosum is a parvulin-like PPIase. Due to its striking similarity to the human parvulins Pin1 and Par14, PinA constitutes an interesting subject for structural and functional studies. Here, we present the first high resolution NMR structure of an archaeal parvulin, PinA, based on 1798 conformational restraints. Structure calculation yields an ensemble of 20 convergent low energy structures with a backbone r.m.s.d. value of 0.6 Å within the secondary structure elements. The overall fold of PinA comprises the β-α₃-β-α-β₂ fold typical for all parvulin structures known so far, but with helix III being a short 3₁₀-helix. A detailed comparison of this high resolution structure of the first archaeal PinA protein with bacterial and eukaryotic parvulin PPIase structures reveals an atypically large catalytic binding site. This feature provides an explanation for cold-adapted protein function. Moreover, the residues in and around 3₁₀-helix III exhibit strong intramolecular dynamics on a microsecond to millisecond timescale and display structural heterogeneity within the NMR ensemble. A putative peptide ligand was found for PinA by phage display and was used for ¹H-¹⁵N-HSQC titrations. Again, the flexible region around 3₁₀-helix III as well as residues of the peptide binding pocket showed the strongest chemical shift perturbations upon peptide binding. The local flexibility of this region also was modulated by ligand binding. A glycine and two positively charged residues are conserved in most parvulin proteins in this flexible loop region, which may be of general functional importance for parvulin-type PPIases.     

A molecular dynamics study of solvent behavior around a protein

The solvent structure and behavior around a protein were examined by analyzing a trajectory of molecular dynamics simulation of thetrp-holorepressor in a periodic box of water. The calculated selfdiffusion coefficient indicated that the solvent within 10 Å of the protein had lower mobility. Examination of the solvent diffusion around different atoms of different kinds of residues showed no general tendency. Thisfact suggested that the solvent mobility is not influenced significantly bythe kind of the atom or residue they solvated. Distribution analysis aroundthe protein revealed two peaks of water oxygen: a sharp one at 2.8 Å around polar and charged atoms and a broad one at ～3.4 Å aroundapolar atoms. The former was stabilized by water–protein hydrogen bonds, and the latter was stabilized by water-lwater hydrogen bonds, suggesting the existence of a hydrophobic shell. An analysis of protein atom–water radial distribution functions confirmed these shell structures around polar or charged atoms and apolar ones. © 1993 Wiley-Liss, Inc.     

A new diagnostic tool for neurocysticercosis is a member of a cestode specific hydrophobic ligand binding protein family

A protein of unknown function has been identified as a key serological tool for diagnosis of human tapeworm neurocysticercosis, a major worldwide neurological disease. Our own sequence analysis predicts that this protein is a member of a newly identified cestode specific oligomeric hydrophobic ligand binding protein family. In this report, using a rat cestode model, we confirm that homologues of this protein can bind fatty acids and their derivatives, and thus suggest a biological function for this key diagnostic tool.     

A targeted molecular dynamics study of WPD loop movement in PTP1B

Targeted molecular dynamics was used to examine the mechanism of WPD loop closure in PTP1B, which is essential for the activity of the enzyme. Two important regions are identified: the R-loop (residues 113-118), which assists in substrate binding, and the S-loop (residues 198-209), which undergoes a conformational change that appears to be vital for the movement of the WPD loop. The S-loop is adjacent to the alpha3-helix, and its conformational change is coupled with a change of interactions between the alpha3- and alpha7-helices. This latter observation is of particular interest in connection with a novel class of allosteric inhibitors of PTP1B [Wiesmann et al., Nat. Struc. Mol. Biol. 11 (2004) 730-737]. These compounds prevent the closure of the WPD loop, forcing the enzyme to remain in a catalytically inactive conformation, by blocking the rearrangement of the alpha3-helix relative to the alpha7-helix.     

Folding profiles of antimicrobial scorpion venom-derived peptides on hydrophobic surfaces: a molecular dynamics study

Abstract

Most helical antimicrobial peptides (AMPs) are usually unfolded in aqueous solution; however they acquire their secondary structure in the presence of a hydrophobic environment such as lipid membranes. Being the biological membranes the main target of many AMPs it is necessary to understand their way of action. Pandinin 2 (Pin2) is an alpha-helical AMP isolated from the venom of the African scorpion Pandinus imperator which shows high antimicrobial activity against Gram-positive bacteria and it is less active against Gram-negative bacteria, nevertheless, it has strong hemolytic activity. Its chemically synthesized Pin2GVG analog has low hemolytic activity while keeping its antimicrobial activity. With the aim of exploring the partition and subsequent folding of these peptides, in this work we report the results of extensive molecular dynamics simulations of Pin2 and Pin2GVG peptides in the presence of 2 hydrophobic environments such as dodecyl-phosphocholine (DPC) micelle and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocoline (POPC) membrane. Our results indicate that Pin2 folds in DPC with a 79% of alpha-helical content, which is in agreement with the experimental results, while in POPC it has 62.5% of alpha-helical content. On the other hand, Pin2GVG presents a higher percentage of alpha-helical structure in POPC and a smaller content in DPC when compared with Pin2. These results can help to better choose the starting structures in future molecular dynamics simulations of AMPs, because these peptides can adopt slightly different conformations depending on the hydrophobic environment.

Communicated by Ramaswamy H. Sarma     

Crystal structures of the BtuF periplasmic-binding protein for vitamin B12 suggest a functionally important reduction in protein mobility upon ligand binding

BtuF is the periplasmic binding protein (PBP) for the vitamin B12 transporter BtuCD, a member of the ATP-binding cassette (ABC) transporter superfamily of transmembrane pumps. We have determined crystal structures of Escherichia coli BtuF in the apo state at 3.0 A resolution and with vitamin B12 bound at 2.0 A resolution. The structure of BtuF is similar to that of the FhuD and TroA PBPs and is composed of two alpha/beta domains linked by a rigid alpha-helix. B12 is bound in the "base-on" or vitamin conformation in a wide acidic cleft located between these domains. The C-terminal domain shares structural homology to a B12-binding domain found in a variety of enzymes. The same surface of this domain interacts with opposite surfaces of B12 when comparing ligand-bound structures of BtuF and the homologous enzymes, a change that is probably caused by the obstruction of the face that typically interacts with this domain by the base-on conformation of vitamin B12 bound to BtuF. There is no apparent pseudo-symmetry in the surface properties of the BtuF domains flanking its B12 binding site even though the presumed transport site in the previously reported crystal structure of BtuCD is located in an intersubunit interface with 2-fold symmetry. Unwinding of an alpha-helix in the C-terminal domain of BtuF appears to be part of conformational change involving a general increase in the mobility of this domain in the apo structure compared with the B12-bound structure. As this helix is located on the surface likely to interact with BtuC, unwinding of the helix upon binding to BtuC could play a role in triggering release of B12 into the transport cavity. Furthermore, the high mobility of this domain in free BtuF could provide an entropic driving force for the subsequent release of BtuF required to complete the transport cycle.     

Prevalent mutations of human prion protein: a molecular modeling and molecular dynamics study

Point mutations in the human prion protein gene, leading to amino acid substitutions in the human prion protein contribute to conversion of PrPC to PrPSc and amyloid formation, resulting in prion diseases such as familial Creutzfeldt-Jakob disease (CJD), Gerstmann-Straussler-Scheinker disease (GSS), and fatal familial insomnia. We have investigated impressions of prevalent mutations including Q217R, D202N, F198S, on the human prion protein and compared the mutant models with wild types. Structural analyses of models were performed with molecular modeling and molecular dynamics simulation methods. According to our results, frequently occurred mutations are observed in conserved and fully conserved sequences of human prion protein and the most fluctuation values occur in the Helix 1 around residues 144-152 and C-terminal end of the Helix 2. Our analysis of results obtained from MD simulation clearly shows that this long-range effect plays an important role in the conformational fluctuations in mutant structures of human prion protein. Results obtained from molecular modeling such as creation or elimination of some hydrogen bonds, increase or decrease of the accessible surface area and molecular surface, loss or accumulation of negative or positive charges on specific positions, and altering the polarity and pKa values, show that amino acid point mutations, though not urgently change the stability of PrP, might have some local impacts on the protein interactions which are required for oligomerization into fibrillar species.     

The molecular basis for ligand specificity in a mouse olfactory receptor: a network of functionally important residues

Sequence differences between members of the mouse olfac-tory receptor MOR42 subfamily (MOR42-3 and MOR42-1) are likely to be the basis for variation in ligand binding preference among these receptors. We investigated the specificity of MOR42-3 for a variety of dicarboxylic acids. We used site-directed mutagenesis, guided by homology modeling and ligand docking studies, to locate functionally important residues. Receptors were expressed in Xenopus oocytes and assayed using high throughput electrophysiology. The importance of the Val-113 residue, located deep within the receptor, was analyzed in the context of interhelical interactions. We also screened additional residues predicted to be involved in ligand binding site, based on comparison of ortholog/paralog pairs from the mouse and human olfactory receptor genomes (Man, O., Gilad, Y., and Lancet, D. (2004) Protein Sci. 13, 240-254). A network of 8 residues in transmembrane domains III, V, and VI was identified. These residues form part of the ligand binding pocket of MOR42-3. C12 dicarboxylic acid did not activate the receptor in our functional assay, yet our docking simulations predicted its binding site in MOR42-3. Binding without activation implied that C12 dicarboxylic acid might act as an antagonist. In our functional assay, C12 dicarboxylic acid did indeed act as an antagonist of MOR42-3, in agreement with molecular docking studies. Our results demonstrate a powerful approach based on the synergy between computational predictions and physiological assays.     

Quantum and molecular dynamics study for binding of macrocyclic inhibitors to human alpha-thrombin

Molecular dynamics simulations followed by quantum mechanical calculation and Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) analysis have been carried out to study binding of proline- and pyrazinone-based macrocyclic inhibitors (L86 and T76) to human alpha-thrombin. Detailed binding interaction energies between these inhibitors and individual protein fragments are calculated using DFT method based on a new quantum mechanical approach for computing protein-ligand interaction energy. The analysis of detailed interaction energies provides insight on the protein-ligand binding mechanism. Study shows that T76 and L86 bind to thrombin in a very similar "inhibition mode" except that T76 has relatively weaker binding interaction with Glu(217). The analysis from quantum calculation of binding interaction is consistent with the MM-PBSA calculation of binding free energy, and the calculated free energies for L86/T76-thrombin binding agree well with the experimental data.     

A network of ubiquitin ligases is important for the dynamics of misfolded protein aggregates in yeast

Quality control ubiquitin ligases promote degradation of misfolded proteins by the proteasome. If the capacity of the ubiquitin/proteasome system is exceeded, then misfolded proteins accumulate in aggregates that are cleared by the autophagic system. To identify components of the ubiquitin/proteasome system that protect against aggregation, we analyzed a GFP-tagged protein kinase, Ste11ΔN(K444R)-GFP, in yeast strains deleted for 14 different ubiquitin ligases. We show that deletion of almost all of these ligases affected the proteostatic balance in untreated cells such that Ste11ΔN(K444R)-GFP aggregation was changed significantly compared with the levels found in wild type cells. By contrast, aggregation was increased significantly in only six E3 deletion strains when Ste11ΔN(K444R)-GFP folding was impaired due to inhibition of the molecular chaperone Hsp90 with geldanamycin. The increase in aggregation of Ste11ΔN(K444R)-GFP due to deletion of UBR1 and UFD4 was partially suppressed by deletion of UBR2 due to up-regulation of Rpn4, which controls proteasome activity. Deletion of UBR1 in combination with LTN1, UFD4, or DOA10 led to a marked hypersensitivity to azetidine 2-carboxylic acid, suggesting some redundancy in the networks of quality control ubiquitin ligases. Finally, we show that Ubr1 promotes clearance of protein aggregates when the autophagic system is inactivated. These results provide insight into the mechanics by which ubiquitin ligases cooperate and provide feedback regulation in the clearance of misfolded proteins.     

Glycosylation is important for binding to human calcitonin receptors

Human calcitonin receptor (hCTR) subtypes contain three or four potential Asn-linked glycosylation sites in their extracellular amino termini. The role of glycosylation in hCTR function has not been identified, but it has been suggested that inhibition of glycosylation does not affect binding or signaling. To determine the role of glycosylation in hCTR biology, we studied the effects of inhibition of glycosylation and of substitution of Asn residues that are potential glycosylation sites. Native and mutated hCTRs were studied after transient expression in monkey kidney COS-1 cells. Tunicamycin, administered as part of a treatment protocol that inhibited glycosylation of all expressed receptors, decreased salmon calcitonin (sCT) binding affinities and signaling potencies at hCTRs with three or four potential glycosylation sites. In hCTR3, which contains three potential glycosylation sites at positions 26, 78, and 83, site-specific substitution of Asn-26 by Ala had no effect on sCT binding affinity or potency, whereas substitution of Asn-78 or Asn-83 lowered sCT affinity and potency. A mutant hCTR3 in which all three Asn residues were substituted with Ala exhibited no high-affinity sCT binding and potencies of several calcitonin analogues that were more than 100-fold lower than that of native hCTR3. Our data show that glycosylation is important for high-affinity binding and potency of calcitonin analogues at hCTRs.     

An alpha-subunit loop structure is required for GM2 activator protein binding by beta-hexosaminidase A

The alpha- and/or beta-subunits of human beta-hexosaminidase A (alphabeta) and B (betabeta) are approximately 60% identical. In vivo only beta-hexosaminidase A can utilize GM2 ganglioside as a substrate, but requires the GM2 activator protein to bind GM2 ganglioside and then interact with the enzyme, placing the terminal GalNAc residue in the active site of the alpha-subunit. A model for this interaction suggests that two loop structures, present only in the alpha-subunit, may be critical to this binding. Three amino acids in one of these loops are not encoded in the HEXB gene, while four from the other are removed posttranslationally from the pro-beta-subunit. Natural substrate assays with forms of hexosaminidase A containing mutant alpha-subunits demonstrate that only the site that is removed from the beta-subunit during its maturation is critical for the interaction. Our data suggest an unexpected biological role for such proteolytic processing events.     

A molecular dynamics study of Beta-Glucosidase B upon small substrate binding

Paenibacillus polymyxa β-glucosidase B (BglB), belongs to a GH family 1, is a monomeric enzyme that acts as an exo-β-glucosidase hydrolysing cellobiose and cellodextrins of higher degree of polymerization using retaining mechanism. A molecular dynamics (MD) simulation was performed at 300 K under periodic boundary condition for 5 ns using the complexes structure obtained from previous docking study, namely BglB-Beta-d-glucose and BglB-Cellobiose. From the root-mean-square deviation analysis, both enzyme complexes were reported to deviate from the initial structure in the early part of the simulation but it was stable afterwards. The root-mean-square fluctuation analysis revealed that the most flexible regions comprised of the residues from 26 to 29, 43 to 53, 272 to 276, 306 to 325 and 364 to 367. The radius of gyration analysis had shown the structure of BglB without substrate became more compact towards the end of the simulation compare to other two complexes. The residues His122 and Trp410 were observed to form stable hydrogen bond with occupancy higher than 10%. In conclusion, the behaviour of BglB enzyme towards the substrate binding was successfully explored via MD simulation approaches.     

Dynamic properties of individual water molecules in a hydrophobic pore lined with acyl chains: a molecular dynamics study

Recently, a certain class of synthetic molecules has been shown to form ion channels, the pore of which is lined with hydrophobic acyl chains [M. Sokabe, in: F. Oosawa, H. Hayashi, T. Yoshioka (Eds.), Transmembrane Signaling and Sensation, JSSP/VNU Science Press BV, Tokyo, 1984, p. 119; F. Hayashi, M. Sokabe, M. Takagi, K. Hayashi, U. Kishimoto, Biochim. Biophys. Acta, 510 (1978) 305; M.J. Pregel, L. Jullien, J. Canceill, L. Lacombe, J.M. Lehn, J. Chem. Soc. Perkin Trans., 2 (1995) 417; Y. Tanaka, Y. Kobuke, M. Sokabe, Angew. Chem. Int. Ed. Engl., 34 (1995) 693; M. Sokabe, Z. Qi, K. Donowaki, H. Ishida, K. Okubo, Biophys. J., 70 (1996) A201; H. Ishida, K. Donowaki, Y. Inoue, Z. Qi, M. Sokabe, Chem. Lett. (1997) p. 953]. As an initial step towards understanding the physical mechanisms of ion permeation across such a hydrophobic pore, systematic molecular dynamics simulations were performed to investigate dynamic and energetic properties of water molecules inside the pore using a dimer of alanine-N'-acylated cyclic peptide as a channel model. Dynamic energy profiles for water molecules indicated that the energy barrier at the middle region of the pore is approximately 2-3 kcal/mol higher than that in the cap water region which was defined as a vicinity region of the channel entrance. Energetics analyses demonstrated that the mutual interactions among intrapore water molecules are the major factor to give favorable interaction (negative energy contribution) for themselves. The pore, despite being lined with acyl chains, has a favorable van der Waals interaction with intrapore water molecules. These results may help to explain why water-filled channels can be formed by the hydrophobic helices in natural channels.     

Enantioselective substrate binding in a monooxygenase protein model by molecular dynamics and docking

The two-component flavoenzyme styrene monooxygenase (SMO) is an efficient alternative to several chemical epoxidation catalysts on a preparative scale. A first homology model of the catalytic domain (StyA) of SMO was constructed (Protein Data Bank ID 2HD8) based on the structure of para-hydroxybenzoate hydroxylase. The StyA protein structure was optimized by restrained molecular dynamics to reproduce specific pre-S binding orientations of styrene. Effects of all 10 point mutations examined were explained by the distance of the site to the styrene and FAD binding sites. Thirteen of 20 ligands could be accommodated in a catalytically active binding orientation, and predicted affinities correlated well with experimental turnover and inhibition. The binding cavity is almost completely hydrophobic except for a hydrogen-bonded network formed by three water molecules, the backbone of residues 300-302, and the flavin ribityl, similar to P293, and three crystal waters in para-hydroxybenzoate hydroxylase suggest that P302, T47, and the waters in StyA are a vital component of the catalytic mechanism. The current optimized and validated StyA model provides a good starting point for elucidation of the structural basis of StyA ligand binding and catalysis. Novel insights in the binding of ligands to SMO/StyA, provided by the current protein model, will aid the rational design of mutants with specific, altered enantioselective properties.     

The Erp protein is anchored at the surface by a carboxy-terminal hydrophobic domain and is important for cell-wall structure in Mycobacterium smegmatis

Erp (Exported Repetitive Protein), also known as P36, Pirg and Rv3810, is a member of a mycobacteria-specific family of extracellular proteins. In pathogenic species, the erp gene has been described as a virulence factor. The Erp proteins comprise three domains. The N- and C-terminal domains are similar in all mycobacterial species, while the central domain consists of a repeated module that differs considerably between species. Here we show that the Erp protein is loosely attached to the surface and that the carboxy-terminal domain, which displays hydrophobic features, anchors Erp at the surface of the bacillus. The hydrophobic region is not necessary for the complementation of the altered colony morphology of a Mycobacterium smegmatis erp- mutant but proved to be necessary to achieve resistance to detergent at wild-type levels.     

Conformational transitions in eosinophil cationic protein: a molecular dynamics study in aqueous environment

Extensive molecular dynamics simulations have been performed on eosinophil cationic protein (ECP). The two structures found in the crystallographic dimer (ECPA and ECPB) have been independently simulated. A small difference in the pattern of the sidechain hydrogen bonds in the starting structure has resulted in interesting differences in the conformations accessed during the simulations. In one simulation (ECPB), a stable equilibrium conformation was obtained and in the other (ECPA), conformational transitions at the level of sidechain interactions were observed. The conformational transitions exhibit the involvement of the solvent (water) molecules with a pore-like construct in the equilibrium conformation and an opening for a large number of water molecules during the transition phase. The details of these transitions are examined in terms of intra-protein hydrogen bonds, protein-water networks and the residence times of water molecules on the polar atoms of the protein. These properties show some significant differences in the region between the N-terminal helix and the loop before the C-terminal strand as a function of different conformations accessed during the simulations. However, the stable hydrogen bonds, the protein-water networks, and the hydration patterns in most part of the protein including the active site are very much similar in both the simulations, indicating the fact that these are intrinsic properties of proteins.     

Triadin binding to the C-terminal luminal loop of the ryanodine receptor is important for skeletal muscle excitation contraction coupling

Ca(2+) release from intracellular stores is controlled by complex interactions between multiple proteins. Triadin is a transmembrane glycoprotein of the junctional sarcoplasmic reticulum of striated muscle that interacts with both calsequestrin and the type 1 ryanodine receptor (RyR1) to communicate changes in luminal Ca(2+) to the release machinery. However, the potential impact of the triadin association with RyR1 in skeletal muscle excitation-contraction coupling remains elusive. Here we show that triadin binding to RyR1 is critically important for rapid Ca(2+) release during excitation-contraction coupling. To assess the functional impact of the triadin-RyR1 interaction, we expressed RyR1 mutants in which one or more of three negatively charged residues (D4878, D4907, and E4908) in the terminal RyR1 intraluminal loop were mutated to alanines in RyR1-null (dyspedic) myotubes. Coimmunoprecipitation revealed that triadin, but not junctin, binding to RyR1 was abolished in the triple (D4878A/D4907A/E4908A) mutant and one of the double (D4907A/E4908A) mutants, partially reduced in the D4878A/D4907A double mutant, but not affected by either individual (D4878A, D4907A, E4908A) mutations or the D4878A/E4908A double mutation. Functional studies revealed that the rate of voltage- and ligand-gated SR Ca(2+) release were reduced in proportion to the degree of interruption in triadin binding. Ryanodine binding, single channel recording, and calcium release experiments conducted on WT and triple mutant channels in the absence of triadin demonstrated that the luminal loop mutations do not directly alter RyR1 function. These findings demonstrate that junctin and triadin bind to different sites on RyR1 and that triadin plays an important role in ensuring rapid Ca(2+) release during excitation-contraction coupling in skeletal muscle.