Roland H. Stote

Structural basis for hijacking of cellular LxxLL motifs by papillomavirus E6 oncoproteins.

Targeting HPV Papillomaviruses infect mammalian epithelial cells and induce cancers, including cervical cancer in humans. Vaccines against human papillomavirus (HPV) can prevent, but not cure, infection. A key viral oncoprotein, E6, acts by binding and inactivating many host proteins. Zanier et al. (p. 694) determined high-resolution crystal structures of bovine papillomavirus bound to a peptide from the focal adhesion protein, paxillin, and of HPV bound to a peptide from the ubiquitin ligase E6AP. The structures show that the peptide binds in a pocket formed by two zinc domains and a linker helix, which represents a promising target for therapeutics. Crystal structures show how a key oncoprotein in human papillomavirus binds host proteins. E6 viral oncoproteins are key players in epithelial tumors induced by papillomaviruses in vertebrates, including cervical cancer in humans. E6 proteins target many host proteins by specifically interacting with acidic LxxLL motifs. We solved the crystal structures of bovine (BPV1) and human (HPV16) papillomavirus E6 proteins bound to LxxLL peptides from the focal adhesion protein paxillin and the ubiquitin ligase E6AP, respectively. In both E6 proteins, two zinc domains and a linker helix form a basic-hydrophobic pocket, which captures helical LxxLL motifs in a way compatible with other interaction modes. Mutational inactivation of the LxxLL binding pocket disrupts the oncogenic activities of both E6 proteins. This work reveals the structural basis of both the multifunctionality and the oncogenicity of E6 proteins.

Protein-protein recognition and interaction hot spots in an antigen-antibody complex: free energy decomposition identifies "efficient amino acids".

The molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) method was applied to the study of the protein–protein complex between a camelid single chain variable domain (cAb‐Lys3) and hen egg white lysozyme (HEL), and between cAb‐Lys3 and turkey egg white lysozyme (TEL). The electrostatic energy was estimated by solving the linear Poisson–Boltzmann equation. A free energy decomposition scheme was developed to determine binding energy hot spots of each complex. The calculations identified amino acids of the antibody that make important contributions to the interaction with lysozyme. They further showed the influence of small structural variations on the energetics of binding and they showed that the antibody amino acids that make up the hot spots are organized in such a way as to mimic the lysozyme substrate. Through further analysis of the results, we define the concept of “efficient amino acids,” which can provide an assessment of the binding potential of a particular hot spot interaction. This information, in turn, can be useful in the rational design of small molecules that mimic the antibody. The implications of using free energy decomposition to identify regions of a protein–protein complex that could be targeted by small molecules inhibitors are discussed. Proteins 2007.

Critical Role of Desolvation in the Binding of 20-Hydroxyecdysone to the Ecdysone Receptor

The insect steroid hormone 20-hydroxyecdysone (20E) binds to its cognate nuclear receptor composed of the ecdysone receptor (EcR) and Ultraspiracle (USP) and triggers the main developmental transitions, in particular molting and metamorphosis. We present the crystal structure of the ligand-binding domains of EcR/USP in complex with 20E at 2.4Å resolution and compare it with published structures of EcR/USP bound to ponasterone A (ponA). ponA is essentially identical to 20E but lacks the 25-OH group of 20E. The structure of 20E-bound EcR indicates that an additional hydrogen bond is formed compared with the ponA-bound receptor, yet, paradoxically, ponA has a significantly higher affinity for EcR than 20E. Theoretical studies based on docking and free energy methods lead to a rationale for understanding the difference in binding affinities between 20E and ponA. Results of the calculations indicate that the favorable contribution from the extra H-bond made by 25-OH of 20E is counterbalanced by its larger desolvation cost compared with that of ponA. The contribution of 25-OH to the binding affinity is further compared with those of 20- and 22-OH groups. Ligands that lack the 20- or 22-OH group are indeed known to bind less favorably to EcR than 20E, an effect opposite to that observed for ponA. The results indicate that their respective contributions to receptor-ligand complex stability reside mostly in their different contributions to solvation/desolvation. Together, the data demonstrate the critical role of ligand desolvation in determining binding affinity, with general implications for the binding of hormones to their cognate nuclear receptors.

Calculation of the aqueous solvation energy and entropy, as well as free energy, of simple polar solutes

With the advent of more powerful computers, the question of calculating thermodynamic quantities, such as the energy and the entropy, in solute-solvent systems is revisited. The calculation of these thermodynamic quantitites was limited in the past by their slow convergence relative to the free energy. Using molecular dynamics simulations, the energy, entropy, and free energy of solvation of NMA and CH(3)NH(2), as well as their relative values, have been determined. Three different methods (the thermodynamic perturbation method, the thermodynamic integration method, and a finite-difference method) are compared. The finite difference method gives the best results and accurate values for the entropy and energy were obtained using a reasonable amount to computer time. The results suggest that a meaningful thermodynamic description of biomolecular processes can be realized with present methods and the available computer time.

Structural insights into the molecular mechanism of vitamin d receptor activation by lithocholic Acid involving a new mode of ligand recognition.

The vitamin D receptor (VDR), an endocrine nuclear receptor for 1α,25-dihydroxyvitamin D3, acts also as a bile acid sensor by binding lithocholic acid (LCA). The crystal structure of the zebrafish VDR ligand binding domain in complex with LCA and the SRC-2 coactivator peptide reveals the binding of two LCA molecules by VDR. One LCA binds to the canonical ligand-binding pocket, and the second one, which is not fully buried, is anchored to a site located on the VDR surface. Despite the low affinity of the alternative site, the binding of the second molecule promotes stabilization of the active receptor conformation. Biological activity assays, structural analysis, and molecular dynamics simulations indicate that the recognition of two ligand molecules is crucial for VDR agonism by LCA. The unique binding mode of LCA provides clues for the development of new chemical compounds that target alternative binding sites for therapeutic applications.

Histone H3 Tails Containing Dimethylated Lysine and Adjacent Phosphorylated Serine Modifications Adopt a Specific Conformation during Mitosis and Meiosis

ABSTRACT Condensation of chromatin, mediated in part by posttranslational modifications of histones, is essential for cell division during mitosis. Histone H3 tails are dimethylated on lysine (Kme2) and become phosphorylated on serine (Sp) residues during mitosis. We have explored the possibility that these double modifications are involved in the establishment of H3 tail conformations during the cell cycle. Here we describe a specific chromatin conformation occurring at Kme2 and adjacently phosphorylated S of H3 tails upon formation of a hydrogen bond. This conformation appears exclusively between early prophase and early anaphase of the mitosis, when chromatin condensation is highest. Moreover, we observed that the conformed H3Kme2Sp tail is present at the diplotene and metaphase stages in spermatocytes and oocytes. Our data together with results obtained by cryoelectron microscopy suggest that the conformation of Kme2Sp-modified H3 tails changes during mitosis and meiosis. This is supported by biostructural modeling of a modified histone H3 tail bound by an antibody, indicating that Kme2Sp-modified H3 tails can adopt at least two different conformations. Thus, the H3K9me2S10p and the H3K27me2S28p sites are involved in the acquisition of specific chromatin conformations during chromatin condensation for cell division.

Molecular Basis for Bcl-2 Homology 3 Domain Recognition in the Bcl-2 Protein Family IDENTIFICATION OF CONSERVED HOT SPOT INTERACTIONS

The proteins of the Bcl-2 family are important regulators of apoptosis, or programmed cell death. These proteins regulate this fundamental biological process via the formation of heterodimers involving both pro- and anti-apoptotic family members. Disruption of the balance between anti- and pro-apoptotic Bcl-2 proteins is the cause of numerous pathologies. Bcl-xl, an anti-apoptotic protein of this family, is known to form heterodimers with multiple pro-apoptotic proteins, such as Bad, Bim, Bak, and Bid. To elucidate the molecular basis of this recognition process, we used molecular dynamics simulations coupled with the Molecular Mechanics/Poisson-Boltzmann Surface Area approach to identify the amino acids that make significant energetic contributions to the binding free energy of four complexes formed between Bcl-xl and pro-apoptotic Bcl-2 homology 3 peptides. A fifth protein-peptide complex composed of another anti-apoptotic protein, Bcl-w, in complex with the peptide from Bim was also studied. The results identified amino acids of both the anti-apoptotic proteins as well as the Bcl-2 homology 3 (BH3) domains of the pro-apoptotic proteins that make strong, recurrent interactions in the protein complexes. The calculations show that the two anti-apoptotic proteins, Bcl-xl and Bcl-w, share a similar recognition mechanism. Our results provide insight into the molecular basis for the promiscuous nature of this molecular recognition process by members of the Bcl-2 protein family. These amino acids could be targeted in the design of new mimetics that serve as scaffolds for new antitumoral molecules.

Far infrared spectra of solid state aliphatic amino acids in different protonation states

Far infrared spectra of zwitterionic, cationic, and anionic forms of aliphatic amino acids in solid state have been studied experimentally. Measurements were done on glycine, L-alanine, L-valine, L-leucine, and L-isoleucine powder samples and film samples obtained from dried solutions prepared at pH ranging from 1 to 13. Solid state density functional theory calculations were also performed, and detailed potential energy distributions were obtained from normal mode results. A good correspondence between experimental and simulated spectra was achieved and this allowed us to propose an almost complete band assignment for the far infrared spectra of zwitterionic forms. In the 700-50 cm(-1) range, three regions were identified, each corresponding to a characteristic set of normal modes. A first region between 700 and 450 cm(-1) mainly contained the carboxylate bending, rocking, and wagging modes as well as the ammonium torsional mode. The 450-250 cm(-1) region was representative of backbone and sidechain skeletal bending modes. At last, the low wavenumber zone, below 250 cm(-1), was characteristic of carboxylate and skeletal torsional modes and of lattice modes. Assignments are also proposed for glycine cationic and anionic forms, but could not be obtained for all aliphatic amino acids due to the lack of structural data. This work is intended to provide fundamental information for the understanding of peptides vibrational properties.

Phosphorylation of the retinoic acid receptor alpha induces a mechanical allosteric regulation and changes in internal dynamics.

Nuclear receptor proteins constitute a superfamily of proteins that function as ligand dependent transcription factors. They are implicated in the transcriptional cascades underlying many physiological phenomena, such as embryogenesis, cell growth and differentiation, and apoptosis, making them one of the major signal transduction paradigms in metazoans. Regulation of these receptors occurs through the binding of hormones, and in the case of the retinoic acid receptor (RAR), through the binding of retinoic acid (RA). In addition to this canonical scenario of RAR activity, recent discoveries have shown that RAR regulation also occurs as a result of phosphorylation. In fact, RA induces non-genomic effects, such as the activation of kinase signaling pathways, resulting in the phosphorylation of several targets including RARs themselves. In the case of RARα, phosphorylation of Ser369 located in loop L9–10 of the ligand-binding domain leads to an increase in the affinity for the protein cyclin H, which is part of the Cdk-activating kinase complex of the general transcription factor TFIIH. The cyclin H binding site in RARα is situated more than 40 Å from the phosphorylated serine. Using molecular dynamics simulations of the unphosphorylated and phosphorylated forms of the receptor RARα, we analyzed the structural implications of receptor phosphorylation, which led to the identification of a structural mechanism for the allosteric coupling between the two remote sites of interest. The results show that phosphorylation leads to a reorganization of a local salt bridge network, which induces changes in helix extension and orientation that affects the cyclin H binding site. This results in changes in conformation and flexibility of the latter. The high conservation of the residues implicated in this signal transduction suggests a mechanism that could be applied to other nuclear receptor proteins.

Force field parameters for the simulation of modified histone tails

We describe the development of force field parameters for methylated lysines and arginines, and acetylated lysine for the CHARMM all‐atom force field. We also describe a CHARMM united‐atom force field for modified sidechains suitable for use with fragment‐based docking methods. The development of these parameters is based on results of ab initio quantum mechanics calculations of model compounds with subsequent refinement and validation by molecular mechanics and molecular dynamics simulations. The united‐atom parameters are tested by fragment docking to target proteins using the MCSS procedure. The all‐atom force field is validated by molecular dynamics simulations of multiple experimental structures. In both sets of calculations, the computational predictions using the force field were compared to the corresponding experimental structures. We show that the parameters yield an accurate reproduction of experimental structures. Together with the existing CHARMM force field, these parameters will enable the general modeling of post‐translational modifications of histone tails.