Stabilization of the integrase‐DNA complex by Mg2+ ions and prediction of key residues for binding HIV‐1 integrase inhibitors

The HIV‐1 integrase is an attractive target for the therapeutics development against AIDS, as no host homologue of this protein has been identified. The integrase strand transfer inhibitors (INSTIs), including raltegravir, specifically target the second catalytic step of the integration process by binding to the DDE motif of the catalytic site and coordinating Mg²⁺ ions. Recent X‐ray crystallographic structures of the integrase/DNA complex from prototype foamy virus allowed to investigate the role of the different partners (integrase, DNA, Mg²⁺ ions, raltegravir) in the complex stability using molecular dynamics (MD) simulations. The presence of Mg²⁺ ions is found to be essential for the stability, whereas the simultaneous presence of raltegravir and Mg²⁺ ions has a destabilizing influence. A homology model of HIV‐1 integrase was built on the basis of the X‐ray crystallographic information, and protein marker residues for the ligand binding were detected by clustering the docking poses of known HIV‐1 integrase inhibitors on the model. Interestingly, we had already identified some of these residues to be involved in HIV‐1 resistance mutations and in the stabilization of the catalytic site during the MD simulations. Classification of protein conformations along MD simulations, as well as of ligand docking poses, was performed by using an original learning method, based on self‐organizing maps. This allows us to perform a more in‐depth investigation of the free‐energy basins populated by the complex in MD simulations on the one hand, and a straightforward classification of ligands according to their binding residues on the other hand. Proteins 2014; 82:466–478. © 2013 Wiley Periodicals, Inc.     

Proline-induced hinges in transmembrane helices: possible roles in ion channel gating

A number of ion channels contain transmembrane (TM) alpha-helices that contain proline-induced molecular hinges. These TM helices include the channel-forming peptide alamethicin (Alm), the S6 helix from voltage-gated potassium (Kv) channels, and the D5 helix from voltage-gated chloride (CLC) channels. For both Alm and KvS6, experimental data implicate hinge-bending motions of the helix in an aspect of channel gating. We have compared the hinge-bending motions of these TM helices in bilayer-like environments by multi-nanosecond MD simulations in an attempt to describe motions of these helices that may underlie possible modes of channel gating. Alm is an alpha-helical channel-forming peptide, which contains a central kink associated with a Gly-x-x-Pro motif in its sequence. Simulations of Alm in a TM orientation for 10 ns in an octane slab indicate that the Gly-x-x-Pro motif acts as a molecular hinge. The S6 helix from Shaker Kv channels contains a Pro-Val-Pro motif. Modeling studies and recent experimental data suggest that the KvS6 helix may be kinked in the vicinity of this motif. Simulations (10 ns) of an isolated KvS6 helix in an octane slab and in a POPC bilayer reveal hinge-bending motions. A pattern-matching approach was used to search for possible hinge-bending motifs in the TM helices of other ion channel proteins. This uncovered a conserved Gly-x-Pro motif in TM helix D5 of CLC channels. MD simulations of a model of hCLC1-D5 spanning an octane slab suggest that this channel also contains a TM helix that undergoes hinge-bending motion. In conclusion, our simulations suggest a model in which hinge-bending motions of TM helices may play a functional role in the gating mechanisms of several different families of ion channels.     

New substituted aminopyrimidine derivatives as BACE1 inhibitors: in silico design,synthesis and biological assays

We report in this work new substituted aminopyrimidine derivatives acting as inhibitors of the catalytic site of BACE1. These compounds were obtained from a molecular modeling study. The theoretical and experimental study reported here was carried out in several steps: docking analysis, Molecular Dynamics (MD) simulations, Quantum Theory Atom in Molecules (QTAIM) calculations, synthesis and bioassays and has allowed us to propose some compounds of this series as new inhibitors of the catalytic site of BACE1. The QTAIM study has allowed us to obtain an excellent correlation between the electronic densities and the experimental data of IC₅₀. Also, using combined techniques (MD simulations and QTAIM calculations) enabled us to describe in detail the molecular interactions that stabilize the different L-R complexes. In addition, our results allowed us to determine what portion of these compounds should be changed in order to increase their affinity with the BACE1. Another interesting result is that a sort of synergism was observed when the effects of these new catalytic site inhibitors were combined with Ac-Tyr5-Pro6-Tyr7-Asp8-Ile9-Pro10-Leu11-NH₂, which we have recently reported as a modulator of BACE1 acting on its exosite.     

Assembling viral channel forming proteins: Vpu from HIV‐1

Different routes of assembly are probed for the transmembrane domain (TMD) of the bitopic membrane protein Vpu from HIV‐1. Vpu is responsible for the amplification of viral release from the host cell. The mode of action includes (i) heteroassembly with host factors and (ii) the formation of homo‐oligomers, which are able to conduct ions across the lipid membrane. Two different routes of assembling short sequences of the N terminus, including the TMD of Vpu, Vpu_1–32, and Vpu_8–26, are presented by using a combination of classical molecular dynamics (MD) simulations combined with a docking approach. The rim of alanines (Ala‐8, ‐11, ‐15, and ‐19) resembles an interlocking motif for the sequential assembly into a dimer and trimer. Simultaneous assembly results in oligomeric bundles (trimers to pentamers) with either tryptophans (Trp‐23) or purely hydrophobic residues facing the center. Bundles, with serines facing the pore (Ser‐24), are energetically not the lowest structures. For pentameric bundles with Ser‐24 facing the pore, no water column develops during a short 25 ns MD simulation. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 517–529, 2013.     

The homochiral metallosupramolecular column structure of rhenium (I) complex

The chiral construction of complex supramolecular systems is receiving much interest, especially with achiral building units. In this paper, we reported the helix supramolecular structure of an achiral rhenium complex crystal, [Re(CO)₃(C₁₆H₁₂N₂)(MeCN)] · ClO₄, which was formed by coordination of Re(CO)₅Cl with 2-methenepyridyl-1-aminonaphthalene and acetonitrile. X-ray diffraction analysis suggests that homohelical column metallosupramolecular architecture results from the array of the six-coordinated Re around and along a line: the helical axis in chiral P1211 (No. 4) space group. There is rare example of constructing helical architectures with both coordinated and hydrogen bonds from achiral units. The helical self-assembly process is induced by the intrastrand interaction and crystal growth. Here, the perchlorate acts as both the counter ion and the bridge linking the neighboring coordinated parts together by weak non-covalent interaction. Furthermore, the two-dimensional arrangement adopts double-edged axe-shaped motif, which is different from the common herringbone or brick-wall pattern in coordination polymers. Our investigation introduced the special arrangement model into complicated helical architecture, especially homochiral single-colony.     

Ab initio simulations of Cu binding sites on the N-terminal region of prion protein

The human prion protein binds Cu²⁺ ions in the octarepeat domain of the N-terminal tail up to full occupancy at pH 7.4. Recent experiments have shown that the HGGG octarepeat subdomain is responsible for holding the metal bound in a square-planar configuration. By using first principle ab initio molecular dynamics simulations of the Car–Parrinello type, the coordination of copper to the binding sites of the prion protein octarepeat region is investigated. Simulations are carried out for a number of structured binding sites. Results for the complexes Cu(HGGGW)(wat), Cu(HGGG), and [Cu(HGGG)]₂ are presented. While the presence of a Trp residue and a water molecule does not seem to affect the nature of the copper coordination, high stability of the bond between copper and the amide nitrogen of deprotonated Gly residues is confirmed in all cases. For the more interesting [Cu(HGGG)]₂ complex, a dynamically entangled arrangement of the two domains with exchange of amide nitrogen bonds between the two copper centers emerges, which is consistent with the short Cu–Cu distance observed in experiments at full copper occupancy.     

Human Rev1 relies on insert-2 to promote selective binding and accurate replication of stabilized G-quadruplex motifs

We previously reported that human Rev1 (hRev1) bound to a parallel-stranded G-quadruplex (G4) from the c-MYC promoter with high affinity. We have extended those results to include other G4 motifs, finding that hRev1 exhibited stronger affinity for parallel-stranded G4 than either anti-parallel or hybrid folds. Amino acids in the αE helix of insert-2 were identified as being important for G4 binding. Mutating E466 and Y470 to alanine selectively perturbed G4 binding affinity. The E466K mutant restored wild-type G4 binding properties. Using a forward mutagenesis assay, we discovered that loss of hRev1 increased G4 mutation frequency >200-fold compared to the control sequence. Base substitutions and deletions occurred around and within the G4 motif. Pyridostatin (PDS) exacerbated this effect, as the mutation frequency increased >700-fold over control and deletions upstream of the G4 site more than doubled. Mutagenic replication of G4 DNA (±PDS) was partially rescued by wild-type and E466K hRev1. The E466A or Y470A mutants failed to suppress the PDS-induced increase in G4 mutation frequency. These findings have implications for the role of insert-2, a motif conserved in vertebrates but not yeast or plants, in Rev1-mediated suppression of mutagenesis during G4 replication.     

A labile regulatory copper ion lies near the T1 copper site in the multicopper oxidase CueO

CueO, a multicopper oxidase, is part of the copper-regulatory cue operon in Escherichia coli, is expressed under conditions of copper stress and shows enhanced oxidase activity when additional copper is present. The 1.7-A resolution structure of a crystal soaked in CuCl2 reveals a Cu(II) ion bound to the protein 7.5 A from the T1 copper site in a region rich in methionine residues. The trigonal bipyramidal coordination sphere is unusual, containing two methionine sulfur atoms, two aspartate carboxylate oxygen atoms, and a water molecule. Asp-439 both ligates the labile copper and hydrogen-bonds to His-443, which ligates the T1 copper. This arrangement may mediate electron transfer from substrates to the T1 copper. Mutation of residues bound to the labile copper results in loss of oxidase activity and of copper tolerance, confirming a regulatory role for this site. The methionine-rich portion of the protein, which is similar to that of other proteins involved in copper homeostasis, does not display additional copper binding. The type 3 copper atoms of the trinuclear cluster in the structure are bridged by a chloride ion that completes a square planar coordination sphere for the T2 copper atom but does not affect oxidase activity.     

A Nanosecond Molecular Dynamics Study of Antiparallel d(G)7 Quadruplex Structures: Effect of the Coordinated Cations

Abstract

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G₇) quadruplex structures with different coordinated ions, namely Na⁺ and K⁺ ion, water and Na⁺ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of ～ 1.5 Å from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K⁺ ions are quite immobile. MD studies at 400K indicate that K⁺ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K⁺ ion coordinated structure is more favourable than that for the Na⁺ ion coordinated antiparallel quadruplex structure.     

Dynamics in an unusual acyl carrier protein from a ladderane lipid-synthesizing organism

Anaerobic ammonium-oxidizing (anammox) bacteria express a distinct acyl carrier protein implicated in the biosynthesis of the highly unusual “ladderane” lipids these organisms produce. This “anammox-specific” ACP, or amxACP, shows several unique features such as a conserved FF motif and an unusual sequence in the functionally important helix III. Investigation of the protein's structure and dynamics, both in the crystal by ensemble refinement and by MD simulations, reveals that helix III adopts a rare six-residue-long 3₁₀-helical conformation that confers a large degree of conformational and positional variability on this part of the protein. This way of introducing structural flexibility by using the inherent properties of 3₁₀-helices appears unique among ACPs. Moreover, the structure suggests a role for the FF motif in shielding the thioester linkage between the protein's prosthetic group and its acyl cargo from hydrolysis.     

Coordination of Ni2+ and Cu2+ to metal ion binding domains of E. coli SlyD protein

The C-terminal region of Escherichia coli SlyD is unstructured and extremely rich in potential metal-binding amino acids, especially in histidine residues. SlyD is able to bind two to seven nickel ions per molecule, in a variety of coordination geometries and coordination numbers. This protein contributes to the insertion of nickel into the hydrogenase precursor protein and it has a peptidyl-prolyl cis/trans-isomerase activity which can be regulated through nickel ions. This inspired us to undertake systematic studies on the coordination ability of two histidine-rich peptides from the C-terminus of the SlyD protein with nickel. Also, it is known that histidine-rich regions are part of a Cu^2 + binding domain involved in copper uptake under conditions of metal starvation in vivo in other bacteria. For this reason we decided to examine the complex formation of Ac-AHGHVHGAHDHHHD-NH₂ and Ac-GHGHDHGHEHG-NH₂ fragments with copper ions, which are also reference metal ions in this study. Experiments were performed in a DMSO/water 30:70 solvent. The Ac-AHGHVHGAHDHHHD-NH₂ and Ac-GHGHDHGHEHG-NH₂ fragments were synthesized and their interactions with Ni^2 + and Cu^2 + ions were studied by potentiometric, mass spectrometric, UV-vis, CD, EPR, and NMR spectroscopic techniques in solution. The results show that the Ac-GHGHDHGHEHG-NH₂ fragment forms equimolar complexes with both nickel and copper ions. At physiological pH, the metal ion is bound only through nitrogens from imidazole sidechain of histidine residues. On the contrary, Ac-AHGHVHGAHDHHHD-NH₂ binds 2 metal ions per molecule, at pH range 5 to 7, even if the 1:2 metal:peptide ratios were used. NMR studies indicate the involvement of all His residues in this pH-range in metal binding of the latter peptide. At higher pH, the stoichiometry changes to 1:1 and the His residues are displaced by amide nitrogens.     

Ion behavior in the selectivity filter of HCN1 channels

Hyperpolarization-activated cyclic-nucleotide gated channels (HCNs) are responsible for the generation of pacemaker currents (I_f or I_h) in cardiac and neuronal cells. Despite the overall structural similarity to voltage-gated potassium (Kv) channels, HCNs show much lower selectivity for K⁺ over Na⁺ ions. This increased permeability to Na⁺ is critical to their role in membrane depolarization. HCNs can also select between Na⁺ and Li⁺ ions. Here, we investigate the unique ion selectivity properties of HCNs using molecular-dynamics simulations. Our simulations suggest that the HCN1 pore is flexible and dilated compared with Kv channels with only one stable ion binding site within the selectivity filter. We also observe that ion coordination and hydration differ within the HCN1 selectivity filter compared with those in Kv and cyclic-nucleotide gated channels. Additionally, the C358T mutation further stabilizes the symmetry of the binding site and provides a more fit space for ion coordination, particularly for Li⁺.     

A comparative molecular dynamics study on BACE1 and BACE2 flap flexibility

Beta-amyloid precursor protein cleavage enzyme1 (BACE1) and beta-amyloid precursor protein cleavage enzyme2 (BACE2), members of aspartyl protease family, are close homologs and have high similarity in their protein crystal structures. However, their enzymatic properties are different, which leads to different clinical outcomes. In this study, we performed sequence analysis and all-atom molecular dynamic (MD) simulations for both enzymes in their ligand-free states in order to compare their dynamical flap behaviors. This is to enhance our understanding of the relationship between sequence, structure and the dynamics of this protein family. Sequence analysis shows that in BACE1 and BACE2, most of the ligand-binding sites are conserved, indicative of their enzymatic property as aspartyl protease members. The other conserved residues are more or less unsystematically localized throughout the structure. Herein, we proposed and applied different combined parameters to define the asymmetric flap motion; the distance, d₁, between the flap tip and the flexible region; the dihedral angle, φ, to account for the twisting motion and the TriCα angle, θ₂ and θ₁. All four combined parameters were found to appropriately define the observed “twisting” motion during the flaps different conformational states. Additional analysis of the parameters indicated that the flaps can exist in an ensemble of conformations, i.e. closed, semi-open and open conformations for both systems. However, the behavior of the flap tips during simulations is different between BACE1 and BACE2. The BACE1 active site cavity is more spacious as compared to that of BACE2. The analysis of 10S loop and 113S loop showed a similar trend to that of flaps, with the BACE1 loops being more flexible and less stable than those of BACE2. We believe that the results, methods and perspectives highlighted in this report would assist researchers in the discovery of BACE inhibitors as potential Alzheimer’s disease therapies.     

Contact pairs of RNA with magnesium ions-electrostatics beyond the Poisson-Boltzmann equation

The electrostatic interaction of RNA with its aqueous environment is most relevant for defining macromolecular structure and biological function. The attractive interaction of phosphate groups in the RNA backbone with ions in the water environment leads to the accumulation of positively charged ions in the first few hydration layers around RNA. Electrostatics of this ion atmosphere and the resulting ion concentration profiles have been described by solutions of the nonlinear Poisson-Boltzmann equation and atomistic molecular dynamics (MD) simulations. Much less is known on contact pairs of RNA phosphate groups with ions at the RNA surface, regarding their abundance, molecular geometry, and role in defining RNA structure. Here, we present a combined theoretical and experimental study of interactions of a short RNA duplex with magnesium (Mg²⁺) ions. MD simulations covering a microsecond time range give detailed hydration geometries as well as electrostatics and spatial arrangements of phosphate-Mg²⁺ pairs, including both pairs in direct contact and separated by a single water layer. The theoretical predictions are benchmarked by linear infrared absorption and nonlinear two-dimensional infrared spectra of the asymmetric phosphate stretch vibration which probes both local interaction geometries and electric fields. Contact pairs of phosphate groups and Mg²⁺ ions are identified via their impact on the vibrational frequency position and line shape. A quantitative analysis of infrared spectra for a range of Mg²⁺-excess concentrations and comparison with fluorescence titration measurements shows that on average 20–30% of the Mg²⁺ ions interacting with the RNA duplex form contact pairs. The experimental and MD results are in good agreement. In contrast, calculations based on the nonlinear Poisson-Boltzmann equation fail in describing the ion arrangement, molecular electrostatic potential, and local electric field strengths correctly. Our results underline the importance of local electric field mapping and molecular-level simulations to correctly account for the electrostatics at the RNA-water interface.     

Functional Analysis of the Na+,K+/H+ Antiporter PeNHX3 from the Tree Halophyte Populus euphratica in Yeast by Model-Guided Mutagenesis

Na⁺,K⁺/H⁺ antiporters are H⁺-coupled cotransporters that are crucial for cellular homeostasis. Populus euphratica, a well-known tree halophyte, contains six Na⁺/H⁺ antiporter genes (PeNHX1-6) that have been shown to function in salt tolerance. However, the catalytic mechanisms governing their ion transport remain largely unknown. Using the crystal structure of the Na⁺/H⁺ antiporter from the Escherichia coli (EcNhaA) as a template, we built the three-dimensional structure of PeNHX3 from P. euphratica. The PeNHX3 model displays the typical TM4-TM11 assembly that is critical for ion binding and translocation. The PeNHX3 structure follows the ‘positive-inside’ rule and exhibits a typical physicochemical property of the transporter proteins. Four conserved residues, including Tyr149, Asn187, Asp188, and Arg356, are indentified in the TM4-TM11 assembly region of PeNHX3. Mutagenesis analysis showed that these reserved residues were essential for the function of PeNHX3: Asn187 and Asp188 (forming a ND motif) controlled ion binding and translocation, and Tyr149 and Arg356 compensated helix dipoles in the TM4-TM11 assembly. PeNHX3 mediated Na⁺, K⁺ and Li⁺ transport in a yeast growth assay. Domain-switch analysis shows that TM11 is crucial to Li⁺ transport. The novel features of PeNHX3 in ion binding and translocation are discussed.     

Characterization of two tropomyosin isoforms from the fast skeletal muscle of bluefin tuna Thunnus thynnusorientalis

Fast skeletal muscle tropomyosin (TM) of tunas is composed of nearly equimolar amount of two isoforms designated α-TM and β-TM expediently based on their migration behavior in SDS-PAGE, whereas corresponding TMs from the other fish species are homogenous (α-type). The presence of β-TM is thus specific to tunas so far. The amino acid sequence of β-TM from bluefin tuna Thunnus thynnus orientalis, which has not been revealed to date unlike α-TM, was successfully obtained in this study by cDNA cloning. The coding region of β-TM cDNA comprised of an open reading frame of 855 bp encoding 284 amino acid residues, like most of the TMs. Unexpectedly, the sequence of β-TM showed high similarity to those of other vertebrate α-type TMs including tuna α-TM. Phylogenetic analysis also showed that β-TM has the closest relationship with α-TM of tuna. This fact was quite unlike the relation of mammalian α- and β-TMs. Based on the distribution of amino acid substitutions, it was suggested that tuna TM isoforms are the products of different genes. By thermodynamic analysis of native and reconstituted TMs, it was demonstrated that β-TM is less thermostable than α-TM. Proteolytic digestion also supported the lower stability of the former.     

Importance of polarization effect in the study of metalloproteins: Application of polarized protein specific charge scheme in predicting the reduction potential of azurin

Molecular dynamics (MD) simulation is commonly used in the study of protein dynamics, and in recent years, the extension of MD simulation to the study of metalloproteins is gaining much interest. Choice of force field is crucial in MD studies, and the inclusion of metal centers complicates the process of accurately describing the electrostatic environment that surrounds the redox centre. Herein, we would like to explore the importance of including electrostatic contribution from both protein and solvent in the study of metalloproteins. MD simulations with the implementation of thermodynamic integration will be conducted to model the reduction process of azurin from Pseudomonas aeruginosa. Three charge schemes will be used to derive the partial charges of azurin. These charge schemes differ in terms of the amount of immediate environment, respective to copper, considered during charge fitting, which ranges from the inclusion of copper and residues in the first coordination sphere during density functional theory charge fitting to the comprehensive inclusion of protein and solvent effect surrounding the metal centre using polarized protein‐specific charge scheme. From the simulations conducted, the relative reduction potential of the mutated azurins respective to that of wild‐type azurin (ΔE_cal) were calculated and compared with experimental values. The ΔE_cal approached experimental value with increasing consideration of environmental effect hence substantiating the importance of polarization effect in the study of metalloproteins. This study also attests the practicality of polarized protein‐specific charge as a computational tool capable of incorporating both protein environment and solvent effect into MD simulations. Proteins 2014; 82:2209–2219. © 2014 Wiley Periodicals, Inc.     

Molecular strategies to replace the structural metal site in the prokaryotic zinc finger domain

The specific arrangement of secondary elements in a local motif often totally relies on the formation of coordination bonds between metal ions and protein ligands. This is typified by the ~ 30 amino acid eukaryotic zinc finger motif in which a β-sheet and an α-helix are clustered around a zinc ion by various combinations of four ligands.     

Reversed micelles to mimic the active site of metalloenzymes

The structure and reactivity of cobalt(II), nickel(II), and copper(II) halides have been investigated in 0.20 M CTAX (X = Cl, Br) |CHCl₃ reversed micelles. The former two metal ions adopt a tetrahedral configuration at low water concentrations in the micelle. The tetrahedral complexes are converted to octahedral aqua complexes by increasing the water concentration (solvochromism) or by lowering the temperature (thermochromism). Upon reaction with imidazole, the tetrahedral cobalt and nickel halide complexes also undergo a structural transformation into an octahedral configuration with imidazole coordination. At low water concentrations, copper halides form a polynuclear complex bridged by halide ions and these halogen bridges are easily broken upon addition of water or imidazole. The copper complexes produced by reaction with imidazole were deduced to be CuIm₂X₂ and CuIm₄X₂ at intermediate and high ligand concentrations, respectively. It was also found that the cupric ion in reversed micelles is readily reduced to the cuprous ion with 2-mercaptoethanol, and the cuprous ion is oxidized to the cupric ion by reaction with hydrogen peroxide.