Self-association of an amphipathic helix peptide inhibitor of HIV-1 integrase assessed by electro spray ionization mass spectrometry in trifluoroethanol/water mixtures

Establishing the auto-associative properties of a molecule in solution can be important for determination of its structure and function. EAA26 (VESMNEELKKIIAQVRAQAEHLKTAY) has been designed to inhibit HIV-1 integrase via formation of a stable coiled-coil structure with a nearly homologous segment in the enzyme. The latter catalyzes the permanent incorporation of a DNA copy of the retrovirus genome into host cell DNA, and is thus essential to the life of the retrovirus. This makes integrase an obvious drug target in the therapy of AIDS. The present work has demonstrated, using electrospray ionization mass spectrometry (ESI-MS), that EAA26 is monomeric in pure water, and tetrameric and dimeric at respectively low and medium concentrations of 2,2,2-trifluoroethanol (TFE), and again monomeric at higher TFE concentrations. Thus, the apolar solvent TFE may contribute to either stabilization or disruption of the intermolecular hydrophobic contacts depending on its concentration in aqueous solution. Previous NMR and ultracentifugation results are thus confirmed, indicating the reliability of ESI-MS for defining the self-association state of biologically relevant peptides in both water and organic-water solutions.     

The role of plastic beta-hairpin and weak hydrophobic core in the stability and unfolding of a full sequence design protein

In this study, the thermal stability of a designed alpha/beta protein FSD (full sequence design) was studied by explicit solvent simulations at three moderate temperatures, 273 K, 300 K, and 330 K. The average properties of the ten trajectories at each temperature were analyzed. The thermal unfolding, as judged by backbone root-mean-square deviation and percentage of native contacts, was displayed with increased sampling outside of the native basin as the temperature was raised. The positional fluctuation of the hairpin residues was significantly higher than that of the helix residues at all three temperatures. The hairpin segment displayed certain plasticity even at 273 K. Apart from the terminal residues, the highest fluctuation was shown in the turn residues 7-9. Secondary structure analysis manifested the structural heterogeneity of the hairpin segment. It was also revealed by the simulation that the hydrophobic core was vulnerable to thermal denaturation. Consistent with the experiment, the I7Y mutation in the double mutant FSD-EY (FSD with mutations Q1E and I7Y) dramatically increased the protein stability in the simulation, suggesting that the plasticity of the hairpin can be partially compensated by a stronger hydrophobic core. As for the unfolding pathway, the breathing of the hydrophobic core and the separation of the two secondary structure elements (alpha helix and beta hairpin) was the initiation step of the unfolding. The loss of global contacts from the separation further destabilized the hairpin structure and also led to the unwinding of the helix.     

Oligomeric beta(beta)(alpha) miniprotein motifs: pivotal role of single hinge residue in determining the oligomeric state.

The role of a single glycine hinge residue in the structure of BBAT1, a beta(beta)(alpha) peptide that forms a discrete homotrimeric structure in solution, was evaluated with 11 new peptide sequences which differ only in the identity of the residue at the hinge position. The integrity of the structure and oligomeric state of the peptides was evaluated by using a combination of analytical ultracentrifugation and circular dichroism spectroscopy. Initially, it was discovered that the glycine hinge adopts backbone dihedral angles favored in D-amino acids and that incorporation of D-alanine at the hinge position stabilizes the trimer species. Subsequently, the effect of the side chains of different D-amino acids at the hinge position was evaluated. While incorporation of polar amino acids led to a destabilization of the oligomeric form of the peptide, only peptides including D-Ser or D-Asp at the hinge position were able to achieve a discrete trimer species. Incorporation of hydrophobic amino acids D-Leu and D-Phe led to oligomerization beyond a trimer to a tetrameric form. The dramatic differences among the thermodynamic stabilities and oligomeric states of these peptides illustrates the pivotal role of the hinge residue in the oligomerization of the beta(beta)(alpha) peptides.     

用分子对接方法研究HIV-1整合酶与病毒DNA的结合模式

用分子对接方法研究了HIV-1整合酶(Integrase, IN)二聚体与3’ 端加工(3’ Processing, 3’-P)前的8 bp及27 bp病毒DNA的相互作用, 并获得IN与27 bp病毒DNA的特异性结合模式. 模拟结果表明, IN有特异性DNA结合区和非特异性DNA结合区; IN二聚体B链的K14, R20, K156, K159, K160, K186, K188, R199和A链的K219, W243, K244, R262, R263是IN结合病毒DNA的关键残基; 并从结构上解释了能使IN发挥活性的病毒DNA的最小长度是15 bp. 通过分析结合能发现, IN与DNA稳定结合的主要因素是非极性相互作用, 而关键残基与病毒DNA相互识别主要依赖于极性相互作用. 模拟结果与实验数据较吻合.     

Stabilizing interactions between aromatic and basic side chains in alpha-helical peptides and proteins. Tyrosine effects on helix circular dichroism

Here we investigate the structures and energetics of interactions between aromatic (Phe or Tyr) and basic (Lys or Arg) amino acids in alpha-helices. Side chain interaction energies are measured using helical peptides, by quantifying their helicities with circular dichroism at 222 nm and interpreting the results with Lifson-Roig-based helix/coil theory. A difficulty in working with Tyr is that the aromatic ring perturbs the CD spectrum, giving an incorrect helicity. We calculated the effect of Tyr on the CD at 222 nm by deriving the intensities of the bands directly from the electronic and magnetic transition dipole moments through the rotational strengths corresponding to each excited state of the polypeptide. This gives an improved value of the helix preference of Tyr (from 0.48 to 0.35) and a correction to the helicity for the peptides containing Tyr. We find that Phe-Lys, Lys-Phe, Phe-Arg, Arg-Phe, and Tyr-Lys are all stabilizing by -0.10 to -0.18 kcal.mol-1 when placed i, i + 4 on the surface of a helix in aqueous solution, despite the great difference in polarity between these residues. Interactions between these side chains have previously been attributed to cation-pi bonds. A survey of protein structures shows that they are in fact predominantly hydrophobic interactions between the CH2 groups of Lys or Arg and the aromatic rings.     

Solvent dependence of PII conformation in model alanine peptides

Alanine residues in two model peptides, the pentapeptide AcGGAGGNH(2) and the 11mer AcO(2)A(7)O(2)NH(2), have been reported to have substantial PII conformation in water. The PII structure in both peptides is sensitive to solvent. In the presence of the organic solvent TFE, the conformation of the pentamer changes from PII to internally H-bonded gamma or beta turns, while the chain with seven alanines forms alpha helix. The PII structure in the 11mer is more stable than that in the shorter peptide as the TFE concentration increases. For the pentamer, a comparison of short-chain aliphatic alcohols to water shows that the PII content decreases in the order water > methanol > ethanol > 2-propanol, linearly according to empirical scales of solvent polarity. Thus, depending on the extent of local solvation as folding progresses, the peptide backbone as modeled by alanine oligomers shifts from PII to internally H-bonded (gamma or beta turn) conformations and to alpha helix in longer segments. On the other hand, the PII content of AcO(2)A(7)O(2)NH(2) increases significantly in the presence of guanidine, as does that of oligoproline peptides, while detergent sodium dodecyl sulfate (SDS) favors alpha helix in this peptide. The shorter peptide does not show a parallel increase in PII with guanidine.     

Examination of the folding of a short alanine-based helical peptide with salt bridges using molecular dynamics simulation

A molecular dynamics simulation of the folding of a short alanine-based helical peptide of 17 residues with three Glu...Lys (i, i + 4) salt bridge pairs, referred to as the AEK17 peptide, was carried out. The simulation gave an estimated simulation folding time of 2.5 ns, shorter than 12 ns for an alanine-based peptide of 16 residues with three Lys residues only, referred to as the AK16 peptide, simulated previously. After folded, the AEK17 peptide had a helical content of 77%, in excellent agreement with the experimentally determined value of 80%. An examination of the folding pathways of AEK17 indicated that the peptide proceeded via three-turn helix conformations more than the helix-turn-helix conformation in the folding pathways. An analysis of interactions indicated that the formation of hydrogen bonds between Lys residue side chains and backbone carbonyls is a major factor in the abundant conformation of the three-turn helix intermediate. The substitution of three Ala with Glu residues reduces the extent of hydrophobic interaction in alanine-based AK peptides with the result that the breaking of the interactions of Lys epsilon-NH3+(side chain)...C=O(backbone) is a major activation action for the AEK17 to achieve a complete fold, in contrast to the AK16 peptide, in which breaking non-native hydrophobic interaction is the rate-determining step.     

Interactions of peptide mimics of hyaluronic acid with the receptor for hyaluronan mediated motility (RHAMM)

Summary Using the hyaluronic acid (HA) binding region of the receptor for hyaluronan-mediated motility (RHAMM) as a model, a molecular perspective for peptide mimicry of the natural ligand was established by comparing the interaction sites of HA and unnatural peptide–ligands to RHAMM. This was accomplished by obtaining a series of octapeptide–ligands through screening experiments that bound to the HA binding domains of RHAMM (amino acids 517–576) and could be displaced by HA. These molecules were computationally docked onto a three-dimensional NMR based model of RHAMM. The NMR model showed that RHAMM(517–576) was a set of three helices, two of which contained the HA binding domains (HABDs) flanking a central groove. The structure was stabilized by hydrophobic interactions from four pairs of Val and Ile side chains extending into the groove. The presence of solvent exposed, positively charged side chains spaced 11Å apart matched the spacing of negative charges on HA. Docking experiments using flexible natural and artificial ligands demonstrated that HA and peptide–mimetics preferentially bound to the second helix that contains HABD-2. Three salt bridges between HA carboxylates and Lys548, Lys553 and Lys560 and two hydrophobic interactions involving Val538 and Val559 were predicted to stabilize the RHAMM-HA complex. The high affinity peptides and HA utilized the same charged residues, with additional contacts to other basic residues. However, hydrophobic contacts do not contribute to affinity for peptide ligand-RHAMM complexes. These results offer insight into how selectivity is achieved in the binding of HA to RHAMM, and how peptide competitors may compete for binding with HA on a single hyaladherin.     

某些环肽及其相应线型肽在不同溶剂中构象的CD谱研究

We used CD spectroscopy to study the conformations of three cyclic peptides (CP10E: cyclo[Glu(OBz1)-Pro-Gly-Glu(OBzl)-Gly]2, CP10K: cyclo[Lys(Z)-Pro-Gly-Lys(Z)-Gly]2, CP12K: cyclo[Phe-Lys(Z)-Pro-Gly-Lys(Z)-Gly]2 and their correspondent linear peptides (LP10E: Boc-[Glu(OBzl)-Pro-Gly-Glu(OBzl)-Gly]2-OPac, LP10K: Boc-[Lys(Z)-Pro-Lys(Z)-Pro]2-OMe, LP 12K: Bao- [-Lys(Z)-Pro-Gly-Lys(Z)-Gly]2- OMe) in three solvents of different polarity (chloroform, acetonitrile, 2,2,2-triliuroethanol), and it was found that all of linear and cyclicpeptides exists asγ-turn conformation in chloroform, however, in TFE＆ CH3CN solutions, the three linear peptides are inβ Ⅱ-turn conformations. CP10E isβI-turn conformation, CP10K ＆CP12K exists in more than one types of turn conformations. On the basis of our experiments, it was concluded: 1) In the presence of conformational constrained amino acids short linear peptides form obvious secondary structure; 2)The solvent polarity has influence on the peptide conformation and this influence on linear peptides is greater than that on cyclic peptides; 3)The backbone of cyclic peptide has constraint effect on its conformation and makes the secondary structure of cyclic peptide different from that of its relative linear peptide. This information might give some cules in the design of bioactive peptides with different receptor selectivity.     

Homology modeling of the receptor binding domain of human thrombopoietin

Platelet production in blood is regulated by a lineage specific humoral factor, thrombopoietin (TPO). The amino terminal domain of TPO (TPO-N) is responsible for the signal transduction mediated by the TPO receptor, c-mpl. From the predicted length of helices we found that TPO-N belongs to the long-chain subfamily of the four-helix bundle cytokine family. We built a three dimensional model of TPO-N by a comparative homology modeling procedure. The four helices of TPO-N with an up-up-down-down topology are stabilized by a tightly packed central hydrophobic core and the extended loop AB makes an additional hydrophobic core with helices B and D outside of the four helix bundle scaffold. An interpretation of the previous site directed mutageneses results in light of the model enabled us to identify two isolated receptor binding sites. The surface made of Lys 136, Lys 138 and Lys 140 in helix D, and Pro 42 and Glu 50 in loop AB forms the first receptor binding site, while the surface of Asp 8, Arg 10 and Lys14 in helix A represents the second binding site for the sequential receptor oligomerization.     

Critical residue that promotes protein dimerization: a story of partially exposed Phe25 in 14-3-3σ

Many proteins exist and function as oligomers. While hydrophobic interactions have been recognized as the major driving force for oligomerization, detailed molecular mechanisms for the assembly are unknown. Here, we used 14-3-3σ as a model protein and investigated the role of hydrophobic residues at the dimeric interface using MD simulations and coimmunoprecipitations. We found that a half-exposed and half-buried residue in the interface, Phe(25), plays a more important role in promoting homodimerization than the hydrophobic core residues by organizing both favorable hydrophobic and hydrophilic interactions. Phe(25) is critical in packing and stabilizing hydrophobic core residues. We conclude that the structural stability of hydrophobic cores is critical for a stable homodimer complex and this stable property can be bestowed by residues outside of hydrophobic core. The important organizing activity of Phe(25) for homodimerization of 14-3-3σ originates from its unique physical location, rigidity, size, and hydrophobicity. Thus, hydrophobic residues that are not deeply buried at the oligomeric interface may play important but different roles from the buried core residues and they may promote oligomerization by organizing co-operativity of core and other residues for favorable hydrophobic and electrostatic interactions.     

Denatured-state ensemble and the early-stage folding of the G29A mutant of the B-domain of protein A

The folding mechanism of the G29A mutant of the B-domain of protein A (BdpA) has been studied by all-atom molecular dynamics simulation using AMBER force field (ff03) and generalized Born continuum solvent model. Started from the extended chain conformation, a total of 16 simulations (400 ns each) at 300 K captured some early folding events of the G29A mutant of BdpA. In one of the 16 trajectories, the G29A mutant folded within 2.8 A (root mean square) of the wild-type NMR structure. We observed that the fast burial of hydrophobic residues was the driving force to bring the distant residues into close proximity. The initiation of the helix I and III occurred during the stage of hydrophobic collapse. The initiation and growth of the helix II was slow. Both the secondary structure formation and the development of the native tertiary contacts suggested a multistage folding process. Clustering analysis indicated that two helix species (helices I and III) could be intermediates. Further analysis revealed that the hydrophobic residues of partially folded helix II formed nativelike hydrophobic contacts with helices I and III that stabilized a nativelike state and delayed the completion of folding of the entire protein. The details of the early folding process were compared with other theoretical and experimental studies. It was found that a nativelike hydrophobic cluster was formed by residues including F(30), I(31), L(34), L(44), L(45), and A(48) that prevented further development of the native structures, and breaking the hydrophobic cluster like this one contributed to the rate-limiting step. This was in complete agreement with the recent kinetic measurements in which mutations of these residues to Gly and Ala substantially increased the folding rates by as much as 60 times. Apparently, destabilization of nonnative states dramatically enhanced the folding rates.     

Interplay of secondary structures and side-chain contacts in the denatured state of BBA1

The denatured state of a miniprotein BBA1 is studied under the native condition with the AMBER/Poisson-Boltzmann energy model and with the self-guided enhanced sampling technique. Forty independent trajectories are collected to sample the highly diversified denatured structures. Our simulation data show that the denatured BBA1 contains high percentage of native helix and native turn, but low percentage of native hairpin. Conditional population analysis indicates that the native helix formation and the native hairpin formation are not cooperative in the denatured state. Side-chain analysis shows that the native hydrophobic contacts are more preferred than the non-native hydrophobic contacts in the denatured BBA1. In contrast, the salt-bridge contacts are more or less nonspecific even if their populations are higher than those of hydrophobic contacts. Analysis of the trajectories shows that the native helix mostly initiates near the N terminus and propagates to the C terminus, and mostly forms from 3(10)-helix/turn to alpha helix. The same analysis shows that the native turn is important but not necessary in its formation in the denatured BBA1. In addition, the formations of the two strands in the native hairpin are rather asymmetric, demonstrating the likely influence of the protein environment. Energetic analysis shows that the native helix formation is largely driven by electrostatic interactions in denatured BBA1. Further, the native helix formation is associated with the breakup of non-native salt-bridge contacts and the accumulation of native salt-bridge contacts. However, the native hydrophobic contacts only show a small increase upon the native helix formation while the non-native hydrophobic contacts stay essentially the same, different from the evolution of hydrophobic contacts observed in an isolated helix folding.     

Structure and folding of glucagon-like peptide-1-(7–36)-amide in aqueous trifluoroethanol studied by NMR spectroscopy

The conformational changes of free, monomeric glucagon-like peptide-1-(7–36)-amide (GLP-1) in aqueous solution with increasing concentrations of 2,2,2-trifluoroethanol (TFE) were monitored by NMR spectroscopy. It was found that GLP-1 gradually assumes a stable, single-stranded helical structure in water solution when the TFE concentration is increased from 0 to 35% (v/v). No further structural changes were observed at higher TFE concentrations. The structure of GLP-1 in 35% TFE was determined from 292 distance restraints and 44 angle restraints by distance geometry, simulating annealing and restrained energy minimization. The helical structure extends from T7 to K28, with a less well-defined region around G16 and a disordered six-residue N-terminal domain. The folding process of GLP-1 from random coil (in water) to helix (in 35% TFE) is initiated by the formation of the C-terminal segment of the helix that is extended gradually towards the N-terminus of the peptide with increasing concentration of TFE. The exchange rates of the slow exchanging amide protons indicate that the C-terminal part of the helix is more stable than the N-terminal part. Copyright © 2001 John Wiley & Sons, Ltd.     

Effect of Phosphorylation on the Structural Behaviour of Peptides Derived from the Intrinsically Disordered C-Terminal Domain of Histone H1.0

To investigate the structural impact of phosphorylation on the human histone H1.0 C-terminal domain, we performed NMR structural studies of model peptides containing a single phosphorylation site: T¹¹⁸-H1.0 (T¹¹⁸PKK motif) and T¹⁴⁰-H1.0 (T¹⁴⁰PVK motif). Both model peptides are mainly disordered in aqueous solution in their non-phosphorylated and phosphorylated forms, but become structured in the presence of trifluoroethanol. The peptides T¹¹⁸-H1.0 and pT¹¹⁸-H1.0 contain two helical regions, a long amphipathic α helix spanning residues 104–115 and a short α/3₁₀ helix (residues 119–123), that are almost perpendicular in T¹¹⁸-H1.0 but have a poorly defined orientation in pT¹¹⁸-H1.0. Peptides T¹⁴⁰-H1.0 and pT¹⁴⁰-H1.0 form very similar α helices between residues 141–147. The TPKK and TPVK motifs show the same backbone conformation, but differ in their side-chain contacts; the Thr and pThr side chains interact with the i+2 Lys side chain in the TPKK motif, and with the i+3 Lys side chain in the TPVK motif. The pT phosphate group in pT¹¹⁸-H1.0 and pT¹⁴⁰-H1.0 has pK_a values below the intrinsic values, which can be explained by non-specific charge–charge interactions with nearby Lys. The non-polar Val in the TPVK motif accounts for the pT¹⁴⁰ pK_a being closer to the intrinsic pK_a value than the pT¹¹⁸ pK_a. Altogether, these results validate that minimalist strategies using model peptides can provide structural details difficult to obtain in short-lived intrinsically disordered proteins and domains.     

Macrocyclic beta-sheet peptides that mimic protein quaternary structure through intermolecular beta-sheet interactions

This paper reports the design, synthesis, and characterization of a family of cyclic peptides that mimic protein quaternary structure through beta-sheet interactions. These peptides are 54-membered-ring macrocycles comprising an extended heptapeptide beta-strand, two Hao beta-strand mimics [JACS 2000, 122, 7654] joined by one additional alpha-amino acid, and two delta-linked ornithine beta-turn mimics [JACS 2003, 125, 876]. Peptide 3a, as the representative of these cyclic peptides, contains a heptapeptide sequence (TSFTYTS) adapted from the dimerization interface of protein NuG2 [PDB ID: 1mio]. 1H NMR studies of aqueous solutions of peptide 3a show a partially folded monomer in slow exchange with a strongly folded oligomer. NOE studies clearly show that the peptide self-associates through edge-to-edge beta-sheet dimerization. Pulsed-field gradient (PFG) NMR diffusion coefficient measurements and analytical ultracentrifugation (AUC) studies establish that the oligomer is a tetramer. Collectively, these experiments suggest a model in which cyclic peptide 3a oligomerizes to form a dimer of beta-sheet dimers. In this tetrameric beta-sheet sandwich, the macrocyclic peptide 3a is folded to form a beta-sheet, the beta-sheet is dimerized through edge-to-edge interactions, and this dimer is further dimerized through hydrophobic face-to-face interactions involving the Phe and Tyr groups. Further studies of peptides 3b-3n, which are homologues of peptide 3a with 1-6 variations in the heptapeptide sequence, elucidate the importance of the heptapeptide sequence in the folding and oligomerization of this family of cyclic peptides. Studies of peptides 3b-3g show that aromatic residues across from Hao improve folding of the peptide, while studies of peptides 3h-3n indicate that hydrophobic residues at positions R3 and R5 of the heptapeptide sequence are important in oligomerization.     

Reversed-phase liquid chromatography as a tool in the determination of the hydrophilicity/hydrophobicity of amino acid side-chains at a ligand-receptor interface in the presence of different aqueous environments. II. Effect of varying peptide ligand hydrophobicity

The present study represents a continuation of our development of a chromatographic model for studying the hydrophobic interactions which characterize the way a ligand binds to its receptor. We have designed 18-residue amphipathic alpha-helical peptides (representing the hydrophobic binding domain of a ligand), where the non-polar face interacts with the non-polar face of a reversed-phase stationary phase (representing a receptor protein with a hydrophobic binding pocket). Two series of amphipathic alpha-helical peptides were subjected to reversed-phase liquid chromatography at pH 2.0, where the "native" Ala-face peptide contains seven Ala residues in its non-polar face and the "native" Leu-face series contains seven Leu residues in its non-polar face. Mutants of the two series were then prepared by replacing one residue in the centre of the non-polar face in both series of peptides, resulting in amino acid side-chains being exposed to a moderately non-polar environment (Ala series) or a very hydrophobic environment (Leu series) surrounding the substitution site. With this model, we have demonstrated that an increase in non-polarity of the ligand enhances hydrophilicity (decreases hydrophobicity) of all amino acids at the ligand-receptor interface, this effect being dependent on the intrinsic hydrophilicity/hydrophobicity of the side-chain. The addition of salt to the aqueous environment surrounding the binding site of the ligand and receptor was also shown to affect the hydrophilicity/hydrophobicity of amino acids in the binding interface. For the Ala-face mutants, the majority of the non-polar side-chains and the three positively charged residues (Arg, His, Lys) showed significant enhancement of hydrophobicity in the presence of salt; in contrast, in the much more hydrophobic environment of the Leu-face mutants, there was a trend of lesser hydrophobicity enhancement and/or significantly more hydrophilicity enhancement in the presence of salt. Our results should have major implications for the understanding of the hydrophilicity/ hydrophobicity of side-chains in varying hydrophobic and aqueous environments.     

大豆Em(LEA1)蛋白保守结构域的结构和聚集特性

采用圆二色谱(CD)和核磁共振波谱(NMR)方法研究了大豆Em(LEA1)蛋白保守基序Em-C和Em-2M多肽在不同环境中的结构及聚集行为.研究表明,在水和DMPG溶液中,两种多肽主要以无规结构形式存在.在50％ TFE溶液中,Em-C多肽折叠结构增加,含疏水残基的部分区域可能形成α-螺旋结构,且分子以二聚体形式存在;而Em-2M则以单体形式存在,且有序结构较少.以上结果表明,环境变化可能导致两种多肽的空间结构和聚集行为改变,这有助于理解Em蛋白在不同环境中的结构特点,及其重要区域在全长蛋白中所起的作用.     

A quantum mechanics/molecular mechanics study of the protein-ligand interaction for inhibitors of HIV-1 integrase

Human immunodeficiency virus type-1 integrase (HIV-1 IN) is an essential enzyme for effective viral replication. Diketo acids such as L-731,988 and S-1360 are potent and selective inhibitors of HIV-1 IN. In this study, we used molecular dynamics simulations, within the hybrid quantum mechanics/molecular mechanics (QM/MM) approach, to determine the protein-ligand interaction energy between HIV-1 IN and L-731,988 and 10 of its derivatives and analogues. This hybrid methodology has the advantage that it includes quantum effects such as ligand polarisation upon binding, which can be very important when highly polarisable groups are embedded in anisotropic environments, as for example in metal-containing active sites. Furthermore, an energy decomposition analysis was performed to determine the contributions of individual residues to the enzyme-inhibitor interactions on averaged structures obtained from rather extensive conformational sampling. Analysis of the results reveals first that there is a correlation between protein-ligand interaction energy and experimental strand transfer into human chromosomes and secondly that the Asn-155, Lys-156 and Lys-159 residues and the Mg(2+) ion are crucial to anti-HIV IN activity. These results may explain the available experimental data.     

Efficient Oligomerization of Aromatic Amino Acids Induced by Gaps in Four-Helix Bundles of DNA or RNA

The formation of peptides from amino acids is one of the processes associated with life. Because of the dominant role of translation in extant biology, peptide-forming processes that are RNA induced are of particular interest. We have previously reported the formation of phosphoramidate-linked peptido RNAs as the products of spontaneous condensation reactions between ribonucleotides and free amino acids in aqueous solution. We now asked whether four-helix bundle (4HB) DNA or RNA folding motifs with a single- or double-nucleotide gap next to a 5’-phosphate can act as reaction sites for phosphoramidate formation. For glycine, this was found to be the case, whereas phenylalanine and tryptophan showed accelerated formation of peptides without a covalent link to the nucleic acid. Free peptides with up to 11 tryptophan or phenylalanine residues were found in precipitates forming in the presence of gap-containing DNA or RNA 4HBs. Control experiments using motifs with just a nick or primer alone did not have the same effect. Because folded structures with a gap in a double helix are likely products of hybridization of strands formed in statistically controlled oligomerization reactions, our results are interesting in the context of prebiotic scenarios. Independent of a putative role in evolution, our findings suggest that for some aromatic amino acids an RNA-induced pathway for oligomerization exists that does not have a discernable link to translation.