Comparative analysis of HIV-1 Tat variants

HIV-1 Tat protein is a crucial element for viral replication; therefore, its inhibition might be exploited against the AIDS infection. To gain insights on the natural variability of this protein, we present a comparative investigation on the relationship between the primary sequences and the experimentally available three-dimensional structures from the HIV-1 Tat variants Z2, BRU, and MAL. Our computational tools include sequence conservation algorithms, structural analysis, electrostatic modeling, and molecular dynamics (MD) simulations. We find that two regions located between residues 10-18 and 41-52 display the highest primary sequence conservation, while the most conserved region among the available structures corresponds approximately to the segment between positions approximately 44 and 50. Furthermore, in spite of their large structural divergence, Tat variants share a common mode for long-range intramolecular interactions. Finally, the flexibility of the Z2, BRU, and MAL variants, as emerging from multinanosecond MD simulations, is rather similar. Based on this work, we conclude that the turnlike region between amino acids 44 and 50 is structurally most conserved, emerging as an important motif for pharmaceutical targeting aimed toward inhibiting Tat action.     

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.     

Peptide models of electrostatic interactions in proteins: NMR studies on two beta-turn tetrapeptides containing Asp-His and Asp-Lys salt bridges

Two model peptides Boc-Asp-Pro-Aib-X-NHMe [X = His (1) and X = Lys (2)] were synthesized to simulate intramolecular electrostatic interactions between ionizable side chains. Conformational analysis by 270-MHz 1H NMR in (CD3)2SO reveals that the backbone secondary structures of these two peptides are stabilized by two strong intramolecular hydrogen bonds, involving the consecutive carboxy-terminal NH groups. 1H NMR chemical shifts were measured in 1, 2, and a protected derivative, Boc-Asp(OBzl)-Pro-Aib-His-NHMe (3). These shifts were also measured for the model compounds Ac-Lys-NHMe, Boc-Asp-NHMe, and Boc-His-NHMe in their different states of ionization. An analysis of the chemical shifts of the ionization-sensitive reporter resonances suggests the formation of a strong intramolecular salt bridge in the lysyl peptide 2 and a bridge of moderate strength in the histidyl peptide 1. A comparison of the temperature dependence of chemical shifts in peptides 1-3 suggests that intramolecular salt bridge formation results in diminished backbone flexibility. The results establish that proximity effects confer far greater stability to intramolecular ion pair interactions vis-a-vis their intermolecular counterparts. The salt bridge interaction in peptide 1 displays a remarkable sensitivity to the dielectric constant of the solvent medium. The results suggest that these peptides are good simulators of the role of salt bridges in the structural dynamics of proteins.     

Protein molecular dynamics with electrostatic force entirely determined by a single Poisson-Boltzmann calculation

The electrostatic force including the intramolecular Coulombic interactions and the electrostatic contribution of solvation effect were entirely calculated by using the finite difference Poisson-Boltzmann method (FDPB), which was incorporated into the GROMOS96 force field to complete a new finite difference stochastic dynamics procedure (FDSD). Simulations were performed on an insulin dimer. Different relative dielectric constants were successively assigned to the protein interior; a value of 17 was selected as optimal for our system. The simulation data were analyzed and compared with those obtained from 500-ps molecular dynamics (MD) simulation with explicit water and a 500-ps conventional stochastic dynamics (SD) simulation without the mean solvent force. The results indicate that the FDSD method with GROMOS96 force field is suitable to study the dynamics and structure of proteins in solution if used with the optimal protein dielectric constant.     

Multisite phosphorylation disrupts arginine-glutamate salt bridge networks required for binding of cytoplasmic linker-associated protein 2 (CLASP2) to end-binding protein 1 (EB1)

A group of diverse proteins reversibly binds to growing microtubule plus ends through interactions with end-binding proteins (EBs). These +TIPs control microtubule dynamics and microtubule interactions with other intracellular structures. Here, we use cytoplasmic linker-associated protein 2 (CLASP2) binding to EB1 to determine how multisite phosphorylation regulates interactions with EB1. The central, intrinsically disordered region of vertebrate CLASP proteins contains two SXIP EB1 binding motifs that are required for EB1-mediated plus-end-tracking in vitro. In cells, both EB1 binding motifs can be functional, but most of the binding free energy results from nearby electrostatic interactions. By employing molecular dynamics simulations of the EB1 interaction with a minimal CLASP2 plus-end-tracking module, we find that conserved arginine residues in CLASP2 form extensive hydrogen-bond networks with glutamate residues predominantly in the unstructured, acidic C-terminal tail of EB1. Multisite phosphorylation of glycogen synthase kinase 3 (GSK3) sites near the EB1 binding motifs disrupts this electrostatic "molecular Velcro." Molecular dynamics simulations and (31)P NMR spectroscopy indicate that phosphorylated serines participate in intramolecular interactions with and sequester arginine residues required for EB1 binding. Multisite phosphorylation of these GSK3 motifs requires priming phosphorylation by interphase or mitotic cyclin-dependent kinases (CDKs), and we find that CDK- and GSK3-dependent phosphorylation completely disrupts CLASP2 microtubule plus-end-tracking in mitosis.     

Charge grouping approaches to calculation of electrostatic forces in molecular dynamics of macromolecules.

Calculation of long-range electrostatic interactions is the most time-consuming step in theoretical simulation of the structure and dynamics of macromolecules. In practice very short cutoff distances are used, which may distort the behavior of the model system. We describe two accurate approaches to calculation of electrostatic forces based on hierarchical grouping of charges into cubes. The first is similar to the O(NlogN) algorithm developed by Barnes, J. and Hut, P., Nature (London) 324, 446-449 (1986), for simulation of a gravitational motion of N bodies. The second approach we formulate for a system with periodic boundary conditions in the nearest image approximation. The calculation of electrostatic interactions and a charge grouping procedure are faster than O(N2). The average inaccuracy in the force introduced by the grouping does not exceed 1%. We describe a small modification of the same approach which makes it suitable for long strongly charged polymers as well. This accurate approach to calculation of electrostatic interactions is illustrated with an example of the dynamics of ions near DNA. Quick equilibration of the ionic distribution is observed during molecular dynamics simulation if electrostatic forces are properly calculated, while the behavior and distribution of ions are less realistic when the conventional cutoff distances are used.     

Electronic structure and conformational properties of <Emphasis Type="Italic">1H</Emphasis>-indole-3-acetic acid

The conformational space of 1H-Indole-3-Acetic Acid (IAA) was scanned using molecular dynamics at semiempirical level, and complemented with functional density calculations at B3LYP/6-31G** level, 14 conformers of lowest energy were obtained. Electronic distributions were analyzed at a higher calculation level, thus improving the basis set (B3LYP/6-311++G**). A topological study based on Bader’s theory (AIM: atoms in molecules) and natural bond orbital (NBO) framework performed with the aim to analyze the stability and reactivity of the conformers allowed the understanding of electronic aspects relevant in the study of the antioxidant properties of IAA. Intramolecular hydrogen bonds were found and were characterized as blue-shifting hydrogen bonding interactions. Furthermore, molecular electrostatic potential maps (MEPs) were obtained and analyzed in the light of AIM and NBO results, thus showing subtle but essential features related not only to reactivity but also with intramolecular weak interactions, charge delocalization and structure stabilization.     

Probing Tat peptide-TAR RNA interactions by psoralen photo-cross-linking

Replication of human immunodeficiency virus type 1 (HIV-1) requires specific interactions of Tat protein with the trans-activation responsive region (TAR) RNA, a 59-base stem-loop structure located at the 5'-end of all HIV mRNAs. We have used a site-specific cross-linking method based on psoralen photochemistry to determine the effect of core residues from the Tat sequence on the protein orientation in the Tat-TAR complex and on the specificity of Tat-TAR binding. We synthesized two Tat fragments, Tat(42-72) and Tat(37-72), and incorporated a psoralen-modified amino acid at position 41 during solid-phase assembly of the peptides. We used these psoralen-Tat conjugates to form specific complexes with TAR RNA. Upon near-ultraviolet irradiation (360 nm), psoralen-Asp41-Tat(37-72) cross-linked to a single site in the TAR RNA sequence. The RNA-protein complex was purified and the cross-link site on TAR RNA was determined by primer extension analysis, which revealed that Asp41 of Tat is close to U42 of the lower stem region of TAR RNA. Specificity of the RNA-peptide cross-linking reactions was determined by competition experiments. Our results show that the addition of only four residues (Cys37-Thr40) from the Tat core region significantly enhanced the specificity of the Tat peptide-TAR interactions without altering the site or chemical nature of the cross-link. These studies provide new insights into RNA-protein recognition that could be useful in designing peptidomimetics for RNA targeting. Such psoralen-peptide conjugates provide a new class of probes for sequence-specific protein-nucleic acid interactions and could be used to selectively control gene expression or to induce site-directed mutations.     

Theoretical studies of the electrostatic interactions in aspartic proteinases, intramolecular interactions in pepsin and penicillopepsin

A semi-empirical approach has been used to estimate the intramolecular electrostatic interactions in pepsin and penicillopepsin. The pH-dependence of the free energy electrostatic term was calculated, and the pH-dependence of the domain interactions has been estimated. As it was shown, the contribution of electrostatic interactions is rather small for the stabilization of the native structure. At the same time the electrostatic repulsion between domains increases with the increase of pH. The later can be the cause of the alkaline denaturation of pepsin and domain mobility.     

Structural study of an RNA aptamer for a Tat protein complexed with ligands.

An RNA aptamer for an HIV Tat protein has been isolated by the in vitro SELEX method. The RNA aptamer binds to the Tat protein 50-100 times more strongly than native TAR RNA does. Here, we have investigated the structure of the RNA aptamer complexed with ligands, partial peptide fragments of the Tat protein or argininamide, by multidimensional 1H/13C/15N NMR. It is strongly suggested that two U:A:U base triples are formed in the RNA aptamer upon binding of ligands. Specific hydrogen bonds between arginine side chains of ligands and guanine bases located adjacent to the base triples are identified. On the basis of many intramolecular and intermolecular NOEs, a structural model of the complex has been constructed.