Theoretical studies on protein-nucleic acid interactions. I. Interaction of positively charged amino acids with nucleic acid fragments

Coulombic interactions between the side chains of charged amino acids (Arg⁺, Lys⁺, and His⁺) and negatively charged phosphate groups of nucleic acid fragments have been studied theoretically. Diribose monophosphate and dideoxyribose monophosphate are chosen as model systems for single-stranded RNA and DNA, respectively. The interaction energies have been calculated by second-order perturbation theory using simplified formulas for individual terms. The interaction energy in this formalism is a sum of electrostatic, polarization, dispersion, and repulsive energies. Our results show that about 90% of the total interaction energy is contributed by the electrostatic term alone. Contribution from the repulsive term exceeds that from the dispersion term. Calculated interaction energies suggest that Lys⁺ and His⁺ form more stable complexes with RNA than with single-stranded DNA. On the other hand, Arg⁺ has a higher affinity for DNA than for RNA. The affinity of nucleic acids for the three amino acids is in the order Lys⁺ > His⁺ > Arg⁺. Further, the basic amino acid residues form more stable complexes with A-DNA than with B-DNA. The role of the Coulombic interactions in the specific recognition of nucleic acids by proteins is discussed.     

Theoretical studies on protein-nucleic acid interactions. II. Hydrogen bonding of amino acid side chains with bases and base pairs of nucleic acids

With a view to understanding the role of hydrogen bonds in the recognition of nucleic acids by proteins, hydrogen bonding between the bases and base pairs of nucleic acids and the amino acids (Asn, Gln, Asp and Glu, and charged residues Arg⁺, Glu^?, and Asp^?) has been studied by a second-order perturbation theory. Binding energies have been calculated for all possible configurations involving a pair of hydrogen bonds between the base (or base pair) and the amino acid residue. Our results show that the hydrogen bonding in these cases has a large contribution from electrostatic interaction. In general, the charged amino acids, compared to the uncharged ones, form more stable complexes with bases or base pairs. The hydrogen-bond energies are an order of magnitude smaller than the Coulombic interaction energies between basic amino acids (Lys⁺, Arg⁺, and His⁺) and the phosphate groups of nucleic acids. The stabilities of the complexes of amino acids Asn, Gln, Asp, and Glu with bases are in the order: G–X > C–X > A–X U–X or T–X, and G · C–X > A · T(U)–X, where X is one of these amino acid residues. It has been shown that Glu^? and Asp^? can recognize guanine in single-stranded nucleic acids; Arg⁺ can recognize G · C base pairs from A · T base pairs in double-stranded structures.     

Conformational Activation of Antithrombin by Heparin Involves an Altered Exosite Interaction with Protease

Heparin allosterically activates antithrombin as an inhibitor of factors Xa and IXa by enhancing the initial Michaelis complex interaction of inhibitor with protease through exosites. Here, we investigate the mechanism of this enhancement by analyzing the effects of alanine mutations of six putative antithrombin exosite residues and three complementary protease exosite residues on antithrombin reactivity with these proteases in unactivated and heparin-activated states. Mutations of antithrombin Tyr²⁵³ and His³¹⁹ exosite residues produced massive 10–200-fold losses in reactivity with factors Xa and IXa in both unactivated and heparin-activated states, indicating that these residues made critical attractive interactions with protease independent of heparin activation. By contrast, mutations of Asn²³³, Arg²³⁵, Glu²³⁷, and Glu²⁵⁵ exosite residues showed that these residues made both repulsive and attractive interactions with protease that depended on the activation state and whether the critical Tyr²⁵³/His³¹⁹ residues were mutated. Mutation of factor Xa Arg¹⁴³, Lys¹⁴⁸, and Arg¹⁵⁰ residues that interact with the exosite in the x-ray structure of the Michaelis complex confirmed the importance of all residues for heparin-activated antithrombin reactivity and Arg¹⁵⁰ for native serpin reactivity. These results demonstrate that the exosite is a key determinant of antithrombin reactivity with factors Xa and IXa in the native as well as the heparin-activated state and support a new model of allosteric activation we recently proposed in which a balance between attractive and repulsive exosite interactions in the native state is shifted to favor the attractive interactions in the activated state through core conformational changes induced by heparin binding.     

Structure-activity relationship of gelatinase biosynthesis-activating pheromone of Enterococcus faecalis

The expression of pathogenicity-related extracellular proteases, namely, gelatinase and serine protease, in Enterococcus faecalis is positively regulated by a quorum-sensing system mediated by an autoinducing peptide called gelatinase biosynthesis-activating pheromone (GBAP). GBAP is an 11-amino-acid-residue cyclic peptide containing a lactone linkage. To study the structure-activity relationship of GBAP, we synthesized a series of GBAP analogues and evaluated their activities by a gelatinase-inducing assay and newly developed receptor-binding assays in which fluorescence-labeled peptides bound onto the FsrC-overexpressing Lactococcus lactis cell surface were observed by fluorescent microscopy and quantified by using a fluorophotometer. Alanine-scanning analysis of GBAP showed that the entire ring region was involved in the GBAP agonist activity, while side chains of the tail region were not strictly recognized. The alanine substitution of Phe⁷ or Trp¹⁰ almost abolished their receptor-binding abilities and GBAP agonist activities, suggesting that these two aromatic side chains are strongly involved in receptor interaction and activation. Furthermore, the Trp¹⁰ substitution with natural and unnatural aromatic amino acids, except pentafluorophenylalanine, caused no loss of agonist activity. This suggested the importance of a negative electrostatic potential created by an π-electron cloud on the aromatic ring surface. Structural analysis of GBAP with nuclear magnetic resonance spectroscopy revealed that the ring region adopted a hairpin-like fold and was tightly packed into a compact form. The side chain of Trp¹⁰ was partially buried in the core structure, contributing to the stabilization of the compact form, while that of Phe⁷ was extended from the core structure into the solvent and was probably directly involved in receptor binding.     

Side-chain dynamics are critical for water permeation through aquaporin-1

Molecular dynamics simulations of aquaporin-1 embedded in a solvated lipid bilayer were carried out to investigate the mechanism of water permeation. The 2.2 Å resolution crystal structure of the bovine protein was used for five independent trajectories. During the equilibration and preparatory steps in which the protein was held fixed, water molecules inside the water channel adopted the same positions as observed in the crystal structure but they did not pass through the channel, suggesting that the dynamic motion of the protein is critical for water permeation. When the protein atoms were allowed to move, the side chains of the two asparagines in the two conserved Asn-Pro-Ala motifs near the center of the channel formed hydrogen bonds with water and helped water molecules move through the channel by actively aligning them for transport. The main-chain oxygen atoms, which were exposed to the pore surface in the crystal structure, also contributed to water transfer. Besides the constriction region observed in the crystal structure (Arg¹⁹⁷, Phe⁵⁸, His¹⁸², and Cys¹⁹¹), we found that His⁷⁶ and Val¹⁵⁵ act as a valve by dynamically blocking water permeation and helping control flow.     

Exploring Strong Interactions in Proteins with Quantum Chemistry and Examples of Their Applications in Drug Design

Objectives

Three strong interactions between amino acid side chains (salt bridge, cation-π, and amide bridge) are studied that are stronger than (or comparable to) the common hydrogen bond interactions, and play important roles in protein-protein interactions.

Methods

Quantum chemical methods MP2 and CCSD(T) are used in calculations of interaction energies and structural optimizations.

Results

The energies of three types of amino acid side chain interactions in gaseous phase and in aqueous solutions are calculated using high level quantum chemical methods and basis sets. Typical examples of amino acid salt bridge, cation-π, and amide bridge interactions are analyzed, including the inhibitor design targeting neuraminidase (NA) enzyme of influenza A virus, and the ligand binding interactions in the HCV p7 ion channel. The inhibition mechanism of the M2 proton channel in the influenza A virus is analyzed based on strong amino acid interactions.

Conclusion

(1) The salt bridge interactions between acidic amino acids (Glu^- and Asp^-) and alkaline amino acids (Arg⁺, Lys⁺ and His⁺) are the strongest residue-residue interactions. However, this type of interaction may be weakened by solvation effects and broken by lower pH conditions. (2) The cation- interactions between protonated amino acids (Arg⁺, Lys⁺ and His⁺) and aromatic amino acids (Phe, Tyr, Trp and His) are 2.5 to 5-fold stronger than common hydrogen bond interactions and are less affected by the solvation environment. (3) The amide bridge interactions between the two amide-containing amino acids (Asn and Gln) are three times stronger than hydrogen bond interactions, which are less influenced by the pH of the solution. (4) Ten of the twenty natural amino acids are involved in salt bridge, or cation-, or amide bridge interactions that often play important roles in protein-protein, protein-peptide, protein-ligand, and protein-DNA interactions.     

Syntheses and biological evaluation of analogs of luteinizing hormone-releasing hormone (LH-RH) modified in position 2, 3, 4 or 5

Syntheses by the conventional methods as well as the chemical, physical and biological properties are described of the following analogs of the LH-releasing hormone (LH-RH): [Leu³]-LH-RH, [Phe³]-LH-RH, [Trp²] [His³]-LH-RH, Des-Trp³-LH-RH, Des-His²-[Phe⁵]-LH-RH, [Ala⁴]-LH-RH, [Phe⁵]-LH-RH and [Ala⁴] [Phe⁵]-LH-RH. $\text{In vivo}$ assays showed that [Leu³]-LH-RH did not release LH in doses as high as 5 – 25 μg, having less than 0.0008% of LH-RH activity, while [Phe³]-LH-RH had 0.43% of the LH-RH activity of natural LH-RH. The LH-RH activities of [Trp²] [His³]-LH-RH, Des-Trp³-LH-RH and Des-His²-[Phe⁵]-LH-RH were extremely low. On the other hand, [Ala⁴]-LH-RH, [Phe⁵]-LH-RH and [Ala⁴] [Phe⁵]-LH-RH had significant LH-RH activity. The structure-activity relationship of LH-RH is discussed on the basis of these findings.     

Conformational analysis of cyclic angiotensin II analogues

Summary Conformationally restricted cyclic analogues of angiotensin II (ANG II), Asp¹-Arg²-Val³-Tyr⁴-Val⁵-His⁶-Pro⁷-Phe⁸, with a link between positions 3 and 5 have considerable biological activity. It is proposed that the spatial arrangement of the pharmacophore groups of Tyr⁴, His⁶ and Phe⁸ side chains and the C-terminal carboxyl group in ANG II and active analogues is similar. Conformational analysis of ANG II and two cyclic analogues c[Sar¹, Lys³,Glu⁵]ANG II and c[Sar¹,Hcy³,Mpt⁵]ANG II was performed, and a geometrical comparison of the low-energy conformations of these compounds allowed one to propose a model of receptor-bound conformation in terms of the spatial arrangement of the pharmacophore groups. This model is characterised by the close spatial location of the His⁶-Phe⁸ side chains and the Tyr⁴ C-terminal carboxyl group and is stabilised by the electrostatic interaction of Arg² and the C-terminal carboxyl group.Abbreviations ANG II angiotensin II - Hcy homocysteine - Mpt trans-4-mercaptoproline     

Homology modeling of 3D structure of human chitinase-like protein CHI3L2

The human genome encodes six proteins of family 18 glycosyl hydrolases, two active chitinases and four chitinase-like lectins (chi-lectins) lacking catalytic activity. The present article is dedicated to homology modeling of 3D structure of human chitinase 3-like 2 protein (CHI3L2), which is overexpressed in glial brain tumors, and its structural comparison with homologous chi-lectin CHI3L1. Two crystal structures of CHI3L1 in free state (Protein Data Bank codes 1HJX and 1NWR) were used as structural templates for the homology modeling by Modeller 9.7 program, and the best quality model structure was selected from the obtained model ensemble. Analysis of potential oligosaccharide-binding groove structures of CHI3L1 and CHI3L2 revealed significant differences between these two homologous proteins. 8 of 19 amino acid residues important for ligand binding are substituted in CHI3L2: Tyr³⁴/Asp³⁹, Trp⁶⁹/Lys⁷⁴, Trp⁷¹/Lys⁷⁶, Trp⁹⁹/Tyr¹⁰⁴, Asn¹⁰⁰/Leu¹⁰⁵, Met²⁰⁴/Leu²¹⁰, Tyr²⁰⁶/Phe²¹² and Arg²⁶³/His²⁷¹. The differences between these residues could influence the structure of the ligand-binding groove and substantially change the ability of CHI3L2 to bind oligosaccharide ligands.     

Free‐energy profiles for ions in the influenza M2‐TMD channel

M₂ transmembrane domain channel (M₂‐TMD) permeation properties are studied using molecular dynamics simulations of M₂‐TMD (1NYJ) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl or KCl saline solution. This study allows examination of spontaneous cation and anion entry into the selectivity filter. Three titration states of the M₂‐TMD tetramer are modeled for which the four His³⁷ residues, forming the selectivity filter, are net uncharged, +2 charged, or +3 charged. M₂‐TMD structural properties from our simulations are compared with the properties of other models extracted from NMR and X‐ray studies. During 10 ns simulations, chloride ions occasionally occupy the positively‐charged selectivity filter region, and from umbrella sampling simulations, Cl^? has a lower free‐energy barrier in the selectivity‐filter region than either Na⁺ or NH, and NH has a lower free‐energy barrier than Na⁺. For Na⁺ and Cl^?, the free‐energy barriers are less than 5 kcal/mol, suggesting that the 1NYJ conformation would probably not be exquisitely proton selective. We also point out a rotameric configuration of Trp⁴¹ that could fully occlude the channel. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

The amino acid sequence of human APOA-I, an apolipoprotein isolated from high density lipoproteins.

The complete amino acid sequence of human A-I has been determined by manual and automated Edman degradation of intact and peptide fragments of A-I. A-I is a single chain protein of 243 residues with the following amino acid composition: Asp₁₆, Asn₅, Thr₁₀, Ser₁₅, Glu₂₇, Gln₁₉, Pro₁₀, Gly₁₀, Ala₁₉, Val₁₃, Met₃, Leu₃₇, Tyr₇, Phe₆, Trp₄, Lys₂₁, His₅, and Arg₁₆. The amino acid sequence contains no linear segments of hydrophobic or hydrophilic residues. A detailed correlation of the amino acid sequence, conformation, and self association of A-I will add further insight into the molecular mechanisms involved in protein-protein and protein-lipid interactions.     

Regulatory role of the extreme C‐terminal end of the S6 inner helix in C‐terminal‐truncated Kv1.2 channel activation

The transmembrane part of the S6 inner helix of the Kv1.2 potassium channel is a pivotal part in sustaining channel activity. However, the role of its extreme C‐terminal end, which is located on the cytoplasmic side of the channel, is largely unknown. Here, we investigated the role of the extreme C‐terminal end of the S6 inner helix (the HRET region) by constructing a series of C‐terminal‐truncated mutations related to this region in the C‐terminal‐truncated Kv1.2 channel. Mutations on Thr⁴²¹ or Glu⁴²⁰ significantly altered the activation of the truncated channel. Mutations on Arg⁴¹⁹ demonstrated that neutral or basic, but not acidic amino acid, is essential at the position for the truncated channel activation, and no functional channel was observed when the channel was truncated from His⁴¹⁸. Hence, our results indicate that the extreme C‐terminal end of the S6 inner helix plays an important regulatory role in the activation of the C‐terminal‐truncated Kv1.2 channel.     

The lipid-associated 3D structure of SPA, a broad-spectrum neuropeptide antagonist with anticancer properties

[D-Arg¹, D-Trp^5,7,9, Leu¹¹] substance P (SPA) belongs to a family of peptides including antagonist G and SpD that act as broad-spectrum neuropeptide antagonists at several peripheral receptors. The lipid-induced structure of these peptides may be important for the receptor interactions of these analogs. Thus we describe the tertiary structure of SPA in the presence of sodium dodecylsulfate micelles at pH 5.0, and 25°C as determined from two-dimensional ¹H-NMR data recorded at 500 MHz. The resulting three-dimensional structure can be generally described as two type IV nonstandard turns around Arg¹*, Pro², Lys³, and Pro⁴ and Gln⁶, Trp⁷*, Phe⁸, and Trp⁹* residues, respectively, inserted into the interfacial region of the micelles (the asterisks denote D-form amino acid). These turns juxtapose the N- and C-termini of SPA and may form the basis of this peptide's unique ability to inhibit peptide receptor interactions at multiple receptor types.     

Determination of β-glucosidase aggregating factor (BGAF) binding and polymerization regions on the maize β-glucosidase isozyme Glu1

β-Glucosidases (Glu1 and Glu2) in maize specifically interact with a lectin called β-glucosidase aggregating factor (BGAF). We have shown that the N-terminal (Glu⁵⁰–Val¹⁴⁵) and the C-terminal (Phe⁴⁶⁶–Ala⁵¹²) regions of maize Glu1 are involved in binding to BGAF. Sequence comparison between sorghum β-glucosidases (dhurrinases, which do not bind to BGAF) and maize β-glucosidases, and the 3D-structure of Glu1 suggested that the BGAF-binding site on Glu1 is much smaller than predicted previously. To define more precisely the BGAF-binding site, we constructed additional chimeric β-glucosidases. The results showed that a region spanning 11 amino acids (Ile⁷²–Thr⁸²) on Glu1 is essential and sufficient for BGAF binding, whereas the extreme N-terminal region Ser¹–Thr²⁹, together with C-terminal region Phe⁴⁶⁶–Ala⁵¹², affects the size of Glu1–BGAF complexes. The dissociation constants (K_d) of chimeric β-glucosidase–BGAF interactions also demonstrated that the extreme N-terminal and C-terminal regions are important but not essential for binding. To confirm the importance of Ile⁷²–Thr⁸² on Glu1 for BGAF binding, we constructed a chimeric sorghum β-glucosidase, Dhr2 (C-11, Dhr2 whose Val⁷²–Glu⁸² region was replaced with the Ile⁷²–Thr⁸² region of Glu1). C-11 binds to BGAF, indicating that the Ile⁷²–Thr⁸² region is indeed a major interaction site on Glu1 involved in BGAF binding.     

Subunit Symmetry at the Extracellular Domain-Transmembrane Domain Interface in Acetylcholine Receptor Channel Gating

Transmitter molecules bind to synaptic acetylcholine receptor channels (AChRs) to promote a global channel-opening conformational change. Although the detailed mechanism that links ligand binding and channel gating is uncertain, the energy changes caused by mutations appear to be more symmetrical between subunits in the transmembrane domain compared with the extracellular domain. The only covalent connection between these domains is the pre-M1 linker, a stretch of five amino acids that joins strand β10 with the M1 helix. In each subunit, this linker has a central Arg (Arg^3′), which only in the non-α-subunits is flanked by positively charged residues. Previous studies showed that mutations of Arg^3′ in the α-subunit alter the gating equilibrium constant and reduce channel expression. We recorded single-channel currents and estimated the gating rate and equilibrium constants of adult mouse AChRs with mutations at the pre-M1 linker and the nearby residue Glu⁴⁵ in non-α-subunits. In all subunits, mutations of Arg^3′ had similar effects as in the α-subunit. In the ϵ-subunit, mutations of the flanking residues and Glu⁴⁵ had only small effects, and there was no energy coupling between ϵGlu⁴⁵ and ϵArg^3′. The non-α-subunit Arg^3′ residues had Φ-values that were similar to those for the α-subunit. The results suggest that there is a general symmetry between the AChR subunits during gating isomerization in this linker and that the central Arg is involved in expression more so than gating. The energy transfer through the AChR during gating appears to mainly involve Glu⁴⁵, but only in the α-subunits.     

Purification of Alkaline Proteinase from Aspergillus Candidus

An alkaline proteinase of Aspergillus Candidus was purified from wheat bran solid culture by batchwise treatment with Amberlite IRC–50 and sequential chromatography on DEAE-cellulose, hydroxylapatite and Sephadex G–100 gel. This purification results in a 18-fold increase of proteolytic activity and the enzyme preparation was homogeneous in sedimentation analysis of the ultracentrifuge and polyacrylamide gel disc electrophoresis. The molecular weight was estimated to be about 23,000 by gel glltration and 22,000 by calculation from the amino acid composition. The enzyme consisted of Lys₁₄, His₄, Arg₃, Asp₂₅, Thr₁₅, Ser₂₃, Glu₁₅, Pro₇, Gly₂₂, Ala₂₄, Met₂, Val₁₆, Ile₁₁, Leu₁₀, Tyr₆, Phe₇, Trp₂ and amide ammonia₁₄ and did not contain cysteine or cystine.     

α‐MSH tripeptide analogs activate the melanocortin 1 receptor and reduce UV‐induced DNA damage in human melanocytes

One skin cancer prevention strategy that we are developing is based on synthesizing and testing melanocortin analogs that reduce and repair DNA damage resulting from exposure to solar ultraviolet (UV) radiation, in addition to stimulating pigmentation. Previously, we reported the effects of tetrapeptide analogs of α‐melanocortin (α‐MSH) that were more potent and stable than the physiological α‐MSH, and mimicked its photoprotective effects against UV‐induced DNA damage in human melanocytes. Here, we report on a panel of tripeptide analogs consisting of a modified α‐MSH core His⁶‐d ‐Phe⁷‐Arg⁸, which contained different N‐capping groups, C‐terminal modifications, or arginine mimics. The most potent tripeptides in activating cAMP formation and tyrosinase of human melanocytes were three analogs with C‐terminal modifications. The most effective C‐terminal tripeptide mimicked α‐MSH in reducing hydrogen peroxide generation and enhancing nucleotide excision repair following UV irradiation. The effects of these three analogs required functional MC1R, as they were absent in human melanocytes that expressed non‐functional receptor. These results demonstrate activation of the MC1R by tripeptide melanocortin analogs. Designing small analogs for topical delivery should prove practical and efficacious for skin cancer prevention.     

Exploring a potential palonosetron allosteric binding site in the 5-HT3 receptor

Palonosetron (Aloxi) is a potent second generation 5-HT₃ receptor antagonist whose mechanism of action is not yet fully understood. Palonosetron acts at the 5-HT₃ receptor binding site but recent computational studies indicated other possible sites of action in the extracellular domain. To test this hypothesis we mutated a series of residues in the 5-HT3A receptor subunit (Tyr⁷³, Phe¹³⁰, Ser¹⁶³, and Asp¹⁶⁵) and in the 5-HT3B receptor subunit (His⁷³, Phe¹³⁰, Glu¹⁷⁰, and Tyr¹⁴³) that were previously predicted by in silico docking studies to interact with palonosetron. Homomeric (5-HT₃A) and heteromeric (5-HT₃AB) receptors were then expressed in HEK293 cells to determine the potency of palonosetron using both fluorimetric and radioligand methods to test function and ligand binding, respectively. The data show that the substitutions have little or no effect on palonosetron inhibition of 5-HT-evoked responses or binding. In contrast, substitutions in the orthosteric binding site abolish palonosetron binding. Overall, the data support a binding site for palonosetron at the classic orthosteric binding pocket between two 5-HT₃A receptor subunits but not at allosteric sites previously identified by in silico modelling and docking.     

Activation and Proton Transport Mechanism in Influenza A M2 Channel

Molecular dynamics trajectories 2 μs in length have been generated for the pH-activated, tetrameric M2 proton channel of the influenza A virus in all protonation states of the pH sensor located at the His³⁷ tetrad. All simulated structures are in very good agreement with high-resolution structures. Changes in the channel caused by progressive protonation of His³⁷ provide insight into the mechanism of proton transport. The channel is closed at both His³⁷ and Trp⁴¹ sites in the singly and doubly protonated states, but it opens at Trp⁴¹ upon further protonation. Anions access the charged His³⁷ and by doing so stabilize the protonated states of the channel. The narrow opening at the His³⁷ site, further blocked by anions, is inconsistent with the water-wire mechanism of proton transport. Instead, conformational interconversions of His³⁷ correlated with hydrogen bonding to water molecules indicate that these residues shuttle protons in high-protonation states. Hydrogen bonds between charged and uncharged histidines are rare. The valve at Val²⁷ remains on average quite narrow in all protonation states but fluctuates sufficiently to support water and proton transport. A proton transport mechanism in which the channel, depending on pH, opens at either the histidine or valine gate is only partially supported by the simulations.