Construction of a ferritin reactor: an efficient means for trapping various heavy metal ions in flowing seawater

An apparatus consisting of two pumps, a mixer, a ferritin reactor, and a spectrophotometer was constructed to study the ability to trap various heavy metal ions (M²⁺) and the dynamics of a reconstituted ferritin reactor in flowing seawater. Reconstituted pig spleen ferritin (PSF_r) is assembled from apo-protein shell to form a reconstituted iron core. The main components of the PSF_r are its core, which contains an Fe²⁺:P_i stoichiometry of 6.0±0.5, reconstituted from pig spleen apoferritin (apo PSF), Fe²⁺, inorganic phosphate (P_i), and O₂ (0.6 atm). The Fe³⁺—P_i clusters within the PSF_r core exhibit resistance to salt ranging from 1% to 6% NaCl. The ferritin reactor consists of PSF_r and an oscillating bag. Using the reactor, M²⁺ ions such as Cd²⁺, Zn²⁺, Co²⁺, and Mn²⁺ are directly trapped by the ferritin. We found a 1:2±0.2 stoichiometry of the trapped M²⁺ to the released iron as measured by chemical analysis or atomic absorption spectrometry; nontransient elements such as Na⁺, K⁺, Ca²⁺, etc., were scarcely trapped by the reactor. This study provides basic conditions for establishing a ferritin reactor and a convenient means for monitoring the pollution of heavy metal ions in seawater.     

Mechanistic insights into the dual inhibition strategy for checking Leishmaniasis

Leishmaniasis ^{1 1}These authors contributed equally.Communicated by Ramaswamy H. SarmaCommunicated by Ramaswamy H. Sarma is an endemic disease mainly caused by the protozoan Leishmania donovani (Ld). Polyamines have been identified as essential organic compounds for the growth and survival of Ld. These are synthesized in Ld by polyamine synthesis pathway comprising of many enzymes such as ornithine decarboxylase (ODC), spermidine synthase (SS), and S-adenosylmethionine decarboxylase. Inhibition of these enzymes in Ld offers a viable prospect to check its growth and development. In the present work, we used computational approaches to search natural inhibitors against ODC and SS enzymes. We predicted three-dimensional structures of ODC and SS using comparative modeling and molecular dynamics (MD) simulations. Thousands of natural compounds were virtually screened against target proteins using high throughput approach. MD simulations were then performed to examine molecular interactions between the screened compounds and functional residues of the active sites of the enzymes. Herein, we report two natural compounds of dual inhibitory nature active against the two crucial enzymes of polyamine pathway of Ld. These dual inhibitors have the potential to evolve as lead molecules in the development of antileishmanial drugs.     

TRPM7 provides an ion channel mechanism for cellular entry of trace metal ions

Trace metal ions such as Zn(2+), Fe(2+), Cu(2+), Mn(2+), and Co(2+) are required cofactors for many essential cellular enzymes, yet little is known about the mechanisms through which they enter into cells. We have shown previously that the widely expressed ion channel TRPM7 (LTRPC7, ChaK1, TRP-PLIK) functions as a Ca(2+)- and Mg(2+)-permeable cation channel, whose activity is regulated by intracellular Mg(2+) and Mg(2+).ATP and have designated native TRPM7-mediated currents as magnesium-nucleotide-regulated metal ion currents (MagNuM). Here we report that heterologously overexpressed TRPM7 in HEK-293 cells conducts a range of essential and toxic divalent metal ions with strong preference for Zn(2+) and Ni(2+), which both permeate TRPM7 up to four times better than Ca(2+). Similarly, native MagNuM currents are also able to support Zn(2+) entry. Furthermore, TRPM7 allows other essential metals such as Mn(2+) and Co(2+) to permeate, and permits significant entry of nonphysiologic or toxic metals such as Cd(2+), Ba(2+), and Sr(2+). Equimolar replacement studies substituting 10 mM Ca(2+) with the respective divalent ions reveal a unique permeation profile for TRPM7 with a permeability sequence of Zn(2+) approximately Ni(2+) > Ba(2+) > Co(2+) > Mg(2+) >/= Mn(2+) >/= Sr(2+) >/= Cd(2+) >/= Ca(2+), while trivalent ions such as La(3+) and Gd(3+) are not measurably permeable. With the exception of Mg(2+), which exerts strong negative feedback from the intracellular side of the pore, this sequence is faithfully maintained when isotonic solutions of these divalent cations are used. Fura-2 quenching experiments with Mn(2+), Co(2+), or Ni(2+) suggest that these can be transported by TRPM7 in the presence of physiological levels of Ca(2+) and Mg(2+), suggesting that TRPM7 represents a novel ion-channel mechanism for cellular metal ion entry into vertebrate cells.     

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.     

Mechanistic insights into transposon cleavage and integration by TnsB of ShCAST system

1. Download : Download high-res image (161KB)
2. Download : Download full-size image

     

Mechanistic insights into mode of action of potent natural antagonists of BACE-1 for checking Alzheimer’s plaque pathology

Alzheimer’s is a neurodegenerative disorder resulting in memory loss and decline in cognitive abilities. Accumulation of extracellular beta amyloidal plaques is one of the major pathology associated with this disease. β-Secretase or BACE-1 performs the initial and rate limiting step of amyloidic pathway in which 37–43 amino acid long peptides are generated which aggregate to form plaques. Inhibition of this enzyme offers a viable prospect to check the growth of these plaques. Numerous efforts have been made in recent years for the generation of BACE-1 inhibitors but many of them failed during the preclinical or clinical trials due to drug related or drug induced toxicity. In the present work, we have used computational methods to screen a large dataset of natural compounds to search for small molecules having BACE-1 inhibitory activity with low toxicity to normal cells. Molecular dynamics simulations were performed to analyze molecular interactions between the screened compounds and the active residues of the enzyme. Herein, we report two natural compounds of inhibitory nature active against β-secretase enzyme of amyloidic pathway and are potent lead molecules against Alzheimer’s disease.     

Mechanistic insights into RNase L through use of an MDMX-derived multi-functional protein domain

RNase L is part of the innate immune response to viral infection. It is activated by a small oligonucleotide (2–5A) whose synthesis is initiated as part of the interferon response. Binding of 2–5A to the N-terminal regulatory region, the ANK domain, of RNase L activates its ribonuclease activity and results in cleavage of RNA in the cell, which ultimately leads to apoptosis of the infected cell. The mechanism by which 2–5A activates the ribonuclease activity of RNase L is currently unclear but 2–5A has been shown to induce dimerization of RNase L. To investigate the importance of dimerization of RNase L, we developed a 15 kDa dimerization-inducing protein domain that was fused to the N-terminus of RNase L. From these studies we provide direct evidence that dimerization of RNase L occurs at physiologically relevant protein concentrations and correlates with activation of ribonuclease activity. We also show that the binding of 2–5A to RNase L promotes dimerization of the ANK domain and suggest how this could transmit a signal to the rest of the protein to activate ribonuclease activity. Finally, we show that the dimerization-inducing domain can be used as a general fusion partner to aid in protein expression and purification.     

Recent advances in ion channel research

The field of ion channels has entered into a rapid phase of development in the last few years, partly due to the breakthroughs in determination of the crystal structures of membrane proteins and advances in computer simulations of biomolecules. These advances have finally enabled the long-dreamed goal of relating function of a channel to its underlying molecular structure. Here we present simplified accounts of the competing permeation theories and then discuss their application to the potassium, gramicidin A and calcium channels.     

Ion binding properties and structure stability of the NaK channel

Ion distribution in the selectivity filter and ion-water and ion-protein interactions of NaK channel are systematically investigated by all-atom molecular dynamics simulations, with the tetramer channel protein being embedded in a solvated phospholipid bilayer. Analysis of the simulation results indicates that K⁺ ions prefer to bind within the sites formed by two adjacent planes of oxygen atoms from the selectivity filter, while Na⁺ ions are inclined to bind to a single plane of four oxygen atoms. At the same time, both K⁺ and Na⁺ ions can diffuse in the vestibule, accompanying with movements of the water molecules confined in a complex formed by the vestibule together with four small grottos connecting to it. As a result, K⁺ ions show a wide range of coordination numbers (6-8), while Na⁺ ions display a constant coordination number of ∼ 6 in the selectivity filter, which may result in the loss of selectivity of NaK. It is also found that a Ca²⁺ can bind at the extracellular site as reported in the crystal structure in a partially hydrated state, or at a higher site in a full hydration state. Furthermore, the carbonyl group of Asp66 can reorient to point towards the center pore when an ion exists in the vestibule, while that of Gly65 always aligns tangentially to the channel axis, as in the crystallographic structures.     

Constructing HIV-1 integrase tetramer and exploring influences of metal ions on forming integrase-DNA complex

HIV-1 integrase (IN) is essential for the replication of HIV-1 in human cells. At present, the complete structure of complex IN-DNA has not been resolved. In this paper, a HIV-1 IN tetramer model was built with homology modeling and molecular dynamics simulation approach, in which two Mg2+ ions were reasonably located in each catalytic core domain. Moreover, it was found that the AB and CD chains of HIV-1 IN tetramer were different in the structures and metal ions of HIV-1 IN tetramer might have great influences on DNA locating on IN. These findings may provide a more complete structural basis for guiding drug discovery and revealing integration mechanism.     

The passage of an ion through a synthetic ion channel

The simulated system consisted of a fatty acid bilayer membrane dividing two electrolyte layers each containing ions, and a channel composed of linked 15-crown-5 ether rings. The Na⁺ and F^?  ions in the aqueous electrolyte layers were too large to enter the channel, but the Li⁺ ions entered and were transported. Conditions that optimised the passive, electric-field-induced transport of Li⁺ ions through the channel were investigated. It was calculated and rationalised that the higher the numerical value of the electrostatic charge on the oxygen atoms of the crown ether rings, the more easily does the channel convey the Li⁺ ions.     

Ion-specific diffusion rates through transmembrane protein channels. A molecular dynamics study

The migration of different alkali metal cations through a transmembrane model channel is simulated by means of the molecular dynamics technique. The parameters of the model are chosen in close relation to the gramicidin A channel. Coulomb- and van der Waals-type potentials between the ions and flexible carbonyl groups of the pore-forming molecule are used to describe the ion channel interaction. The diffusion properties of the ions are obtained from three-dimensional trajectory calculations. The diffusion rates for the different ions Li+, Na+, K+ and Rb+ are affected not only by the mass of the particles but also very strongly by their size. The latter effect is more pronounced for rigid channels, i.e., for binding vibrational frequencies of the CO groups with v greater than 400 cm-1. In this range the selectivity sequence for the diffusion rates is the inverse of that expected from normal rate theory but agrees with that found in experiments for gramicidin A.     

Synthesis and structure-activity relationship of elastase inhibiting novel ethylated thiazole-triazole acetamide hybrids: Mechanistic insights through kinetics and computational contemplations

Keeping in mind the pharmacological importance of 2-aminothiazole and 1,2,4-triazole heterocyclic moieties, a series of novel ethylated bi-heterocyclic acetamide hybrids, 9a–p, was synthesized in a multi-step protocol. The structures of newly synthesized compounds were characterized by ¹H NMR, ¹³C NMR, IR and EI-MS spectral studies. The inhibitory effects of these bi-heterocyclic acetamides (9a–n) were evaluated against elastase and all these molecules were identified as potent inhibitors relative to the standard used. The Kinetics mechanism was analyzed by Lineweaver-Burk plots which revealed that, 9h, inhibited elastase competitively by forming an enzyme-inhibitor complex. The inhibition constants K_i calculated from Dixon plots for this compound was 0.9 µM. The computational study was articulate with the experimental results and these ligands unveiled good binding energy values (kcal/mol). So, these molecules can be considered as promising medicinal scaffolds for the treatment of skin melanoma, wrinkle formation, uneven pigmentation, and solar elastosis.     

Mixed modes in opening of KcsA potassium channel from a targeted molecular dynamics simulation

Potassium channels conduct K⁺ flow selectively across the membrane through a central pore. During a process called gating, the potassium channels undergo a conformational change that opens or closes the ion-conducting pore. The potassium channel KcsA has been structurally determined in its closed state. However, the dynamic mechanism of the gating transition of the KcsA channel is still being investigated. Here, a targeted molecular dynamics simulation up to 150 ns is performed to investigate the detailed opening process of the KcsA channel with an open Kv1.2 structure serving as the target. The channel arrived at a self-determined quasi-stable state within 60 ns. The rigid-body and hinge-bending modes are observed mixed together in the remaining 90 ns long quasi-stable state. The mixed-mode movement seems come from the competition between the helix rigidity and the biased-applied gating force.     

Proton channel models: Filling the gap between experimental data and the structural rationale

Voltage-gated proton channels are integral membrane proteins with the capacity to permeate elementary particles in a voltage and pH dependent manner. These proteins have been found in several species and are involved in various physiological processes. Although their primary topology is known, lack of details regarding their structures in the open conformation has limited analyses toward a deeper understanding of the molecular determinants of their function and regulation. Consequently, the function-structure relationships have been inferred based on homology models. In the present work, we review the existing proton channel models, their assumptions, predictions and the experimental facts that support them. Modeling proton channels is not a trivial task due to the lack of a close homolog template. Hence, there are important differences between published models. This work attempts to critically review existing proton channel models toward the aim of contributing to a better understanding of the structural features of these proteins.     

Molecular dynamics simulations of human and dog gastric lipases: insights into domain movements

Mammalian gastric lipases are stable and active under acidic conditions and also in the duodenal lumen. There has been considerable interest in acid stable lipases owing to their potential application in the treatment of pancreatic exocrine insufficiency. In order to gain insights into the domain movements of these enzymes, molecular dynamics simulations of human gastric lipase was performed at an acidic pH and under neutral conditions. For comparative studies, simulation of dog gastric lipase was also performed at an acidic pH. Analyses show, that in addition to the lid region, there is another region of high mobility in these lipases. The potential role of this novel region is discussed.     

Modelling the effect of a GHz electric field on the dynamics of K+ ions in KcsA potassium channel

The dynamics of potassium ions in a KcsA channel, located within a stochastically fluctuating medium, is modelled via the application of the molecular dynamics simulation method. We investigate the effect of presence and absence of an applied electric field, either constant or periodic, on the dynamics of the channel. It is found that the ions undergo a hopping motion when the channel is exposed to a constant electric field of strength 0.03 V/nm. Furthermore, an alternating electric field in the GHz range, normally present in the daily environment and encountered by most biological systems, is applied to the channel, showing that in this frequency range, the rigidity of the atomic bonds of the filter is increased, which in turn disturbs the ionic passage rate through the filter. Consequently, in this frequency range, the application of electric fields may affect the function of such channels.     

Mechanistic insights into the role of vitamin D and computational identification of potential lead compounds for Parkinson's disease

Parkinson's disease (PD) is the second most common neurodegenerative disorder that affects dopaminergic neurons in the midbrain. A recent study suggests that Orphan Nuclear Receptor 1 (NURR1) impairment may contribute to PD pathogenesis. Our study found three potent agonists for NURR1 protein based on structural and ligand-based screening methods. The pharmacophore is comprised of a hydrogen bond donor, a hydrophobic group, and two aromatic rings (DHRR). The Pharmacophore screening method screened 3142 compounds, of which 3 were screened using structure-based screening. An analysis of the molecules using Molecular Mechanics-Generalized Born Surface Area (binding free energy) revealed a range of −46.77 to −59.06 Kcal/mol. After that, chemical reactivity was investigated by density functional theory, and molecular dynamics simulation was performed (protein-ligand stability). Based on the computational studies, Lifechemical_16901310, Maybridge_2815310, and NPACT_392450 are promising agonists with respect to NURR1. To confirm the potency of the identified compounds, further validation and experiments must be conducted.     

Computer simulation of the KvAP voltage-gated potassium channel: steered molecular dynamics of the voltage sensor

The recent crystal structures of the voltage-gated potassium channel KvAP and its isolated voltage-sensing 'paddle' (composed of segments S1-S4) challenge existing models of voltage gating and raise a number of questions about the structure of the physiologically relevant state. We investigate a possible gating mechanism based on the crystal structures in a 10 ns steered molecular dynamics simulation of KvAP in a membrane-mimetic octane layer. The structure of the full KvAP protein has been modified by restraining the S2-S4 domain to the conformation of the isolated high-resolution paddle structure. After an initial relaxation, the paddle tips are pulled through the membrane from the intracellular to the extracellular side, corresponding to a putative change from closed to open. We describe the effect of this large-scale motion on the central pore domain, which remains largely unchanged, on the protein hydrogen-bonding network and on solvent. We analyze the motion of the S3b-S4 portion of the protein and propose a possible coupling mechanism between the paddle motion and the opening of the channel. Interactions between the arginine residues in S4, solvent and chloride ions are likely to play a role in the gating charge.