Hydrophobic Gating and 1/f Noise of the Anthrax Toxin Channel

Synthesis, computational docking and molecular dynamics studies of a new class of spiroquinoxalinopyrrolidine embedded chromanone hybrids as potent anti-cholinesterase agents

Novel structurally intriguing heterocycles embedded with spiropyrrolidine, quinoxaline and chromanone units were synthesized in good yields using a [Bmim]Br accelerated multicomponent reaction strategy. The key step of the reaction is 1,3-dipolar cycloaddition involving highly functionalized dipolarophile, viz. 3-benzylidenechroman-4-one, to afford spiroquinoxalinopyrrolidine embedded chromanone hybrid heterocycles. The formation of spiro products occurs via two C-C, two N-C and one C-N bonds possessing four adjoining stereogenic centers, two of which are spiro carbons. The newly synthesized spiro compounds showed potent acetylcholinesterase and butyrylcholinesterase inhibitory activities. Moreover, compounds with fluorine displayed the highest AChE (3.20 ± 0.16 μM) and BChE (18.14 ± 0.06 μM) inhibitory activities. Further, docking studies, followed by all-atom molecular dynamics, showed results that are consistent with in vitro experimental findings. Although docking scores for the synthesized derivatives were higher than those of the standard drug, MD MMPBSA results showed better binding of synthesized derivatives (-93.5 ± 11.9 kcal mol-1) compared to the standard drug galantamine (-66.2 ± 12.3 kcal mol-1) for AChE but exhibited similar values (-98.1 ± 11.2 and -97.9 ± 11.5 kcal mol-1) for BChE. These differences observed in drug binding with AChE/BChE are consistent with RMSD, RMSF, LIG plots, and FEL structural analysis. Taken together, these derivatives could be potential candidates as inhibitors of AChE and BChE.

In silico energetic and molecular dynamic simulations studies demonstrate potential effect of the point mutations with implications for protein engineering in BDNF

Protein engineering by directed evolution is time-consuming. Hence, in silico techniques like FoldX-Yasara for ∆∆G calculation, and SNPeffect for predicting propensity for aggregation, amyloid formation, and chaperone binding are employed to design proteins. Here, we used in silico techniques to engineer BDNF-NTF3 interaction and validated it using mutations with known functional implications for NGF dimer. The structures of three mutants representing a positive, negative, or neutral ∆∆G involving two interface residues in BDNF and two mutations representing a neutral and positive ∆∆G in NGF, which is aligned with BDNF, were selected for molecular dynamics (MD) simulation. Our MD results conclude that the secondary structure of individual protomers of the positive and negative mutants displayed a similar or different conformation from the NTF3 monomer, respectively. The positive mutants showed fewer hydrophobic interactions and higher hydrogen bonds compared to the wild-type, negative, and neutral mutants with similar SASA, suggesting solvent-mediated disruption of hydrogen-bonded interactions. Similar results were obtained for mutations with known functional implications for NGF and BDNF. The results suggest that mutations with known effects in homologous proteins could help in validation, and in silico directed evolution experiments could be a viable alternative to the experimental technique used for protein engineering.

Computational analysis of sodium-dependent phosphate transporter SLC20A1/PiT1 gene identifies missense variations C573F, and T58A as high-risk deleterious SNPs

SLC20A1/PiT1 is a sodium-dependent inorganic phosphate transporter, initially recognized as the retroviral receptor for Gibbon Ape Leukemia Virus in humans. SNPs in SLC20A1 is associated with Combined Pituitary Hormone Deficiency and Sodium Lithium Counter transport. Using in silico techniques, we have screened the nsSNPs for their deleterious effect on the structure and function of SLC20A1. Screening with sequence and structure-based tools on 430 nsSNPs, filtered 17 nsSNPs which are deleterious. To evaluate the role of these SNPs, protein modeling and MD simulations were performed. A comparative analysis of model generated with SWISS-MODEL and AlphaFold shows that many residues are in the disallowed region of Ramachandran plot. Since SWISS-MODEL structure has a 25-residue deletion, the AlphaFold structure was used to perform MD simulation for equilibration and structure refinement. Further, to understand perturbation of energetics, we performed in silico mutagenesis and ΔΔG calculation using FoldX on MD refined structures, which yielded SNPs that are neutral (3), destabilizing (12) and stabilizing (2) on protein structure. Furthermore, to elucidate the impact of SNPs on structure, we performed MD simulations to discern the changes in RMSD, Rg, RMSF and LigPlot of interacting residues. RMSF profiles of representative SNPs revealed that A114V (neutral) and T58A (positive) were more flexible & C573F (negative) was more rigid compared to wild type, which is also reflected in the changes in number of local interacting residues in LigPlot and ΔΔG. Taken together, our results show that SNPs can lead to structural perturbations and impact the function of SLC20A1 with potential implications for disease.Communicated by Ramaswamy H. Sarma.

Anthrax Toxin: Model System for Studying Protein Translocation

Dedicated translocase channels are nanomachines that often, but not always, unfold and translocate proteins through narrow pores across the membrane. Generally, these molecular machines utilize external sources of free energy to drive these reactions, since folded proteins are thermodynamically stable, and once unfolded they contain immense diffusive configurational entropy. To catalyze unfolding and translocate the unfolded state at appreciable timescales, translocase channels often utilize analogous peptide-clamp active sites. Here we describe how anthrax toxin has been used as a biophysical model system to study protein translocation. The tripartite bacterial toxin is composed of an oligomeric translocase channel, protective antigen (PA), and two enzymes, edema factor (EF) and lethal factor (LF), which are translocated by PA into mammalian host cells. Unfolding and translocation are powered by the endosomal proton gradient and are catalyzed by three peptide-clamp sites in the PA channel: the α clamp, the ϕ clamp, and the charge clamp. These clamp sites interact nonspecifically with the chemically complex translocating chain, serve to minimize unfolded state configurational entropy, and work cooperatively to promote translocation. Two models of proton gradient driven translocation have been proposed: (i) an extended-chain Brownian ratchet mechanism and (ii) a proton-driven helix-compression mechanism. These models are not mutually exclusive; instead the extended-chain Brownian ratchet likely operates on β-sheet sequences and the helix-compression mechanism likely operates on α-helical sequences. Finally, we compare and contrast anthrax toxin with other related and unrelated translocase channels.

Chloroquine-analogues block anthrax protective antigen channels in steady-state and kinetic studies

The tripartite anthrax toxin from Bacillus anthracis represents the prototype of A-B type of toxins, where the effector A (an enzymatic subunit) is transported with the help of a binding component B into a target cell. Anthrax toxin consists of three different molecules, two effectors, lethal factor (LF) and edema factor (EF) and the binding component also known as protective antigen (PA). PA forms heptamers or octamers following binding to host cell's receptors and mediates the translocation of the effectors into the cytosol via the endosomal pathway. The cation-selective PA₆₃-channel is able to reconstitute in lipid membranes and can be blocked by chloroquine and other heterocyclic compounds. This suggests that the PA₆₃-channel contains a binding site for quinolines. In this study, we investigated the structure-function relationship of different quinolines for the block of the PA₆₃-channel. The affinity of the different chloroquine analogues to the PA₆₃-channel as provided by the equilibrium dissociation constant was measured using titrations. Some quinolines had a much higher affinity to the PA₆₃-channel than chloroquine itself. We also performed ligand-induced current noise measurements using fast Fourier transformation to get insight in the kinetics of the binding of some quinolines to the PA₆₃-channel. The on-rate constants of ligand binding were around 10⁸M^-1·s^-1 at 150mM KCl and were only little dependent on the individual quinoline. The off-rates varied between 4s^-1 and 160s^-1 and depended much more on the structure of the molecules than the on-rate constants. The possible use of the 4-aminoquinolines as a therapy is discussed.

Anthrax toxin channel: What we know based on over 30 years of research

Protective antigen channel is the central component of the deadly anthrax exotoxin responsible for binding and delivery of the toxin's enzymatic lethal and edema factor components into the cytosol. The channel, which is more than three times longer than the lipid bilayer membrane thickness and has a 6-Å limiting diameter, is believed to provide a sophisticated unfoldase and translocase machinery for the foreign protein transport into the host cell cytosol. The tripartite toxin can be reengineered, one component at a time or collectively, to adapt it for the targeted cancer therapeutic treatments. In this review, we focus on the biophysical studies of the protective antigen channel-forming activity, small ion transport properties, enzymatic factor translocation, and blockage comparing it with the related clostridial binary toxin channels. We address issues linked to the anthrax toxin channel structural dynamics and lipid dependence, which are yet to become generally recognized as parts of the toxin translocation machinery.