Plasmodium Topoisomerase VIB and Spo11 Constitute Functional Type IIB Topoisomerase in Malaria Parasite: Its Possible Role in Mitochondrial DNA Segregation.

The human malaria parasite undergoes a noncanonical cell division, namely, endoreduplication, where several rounds of nuclear, mitochondrial, and apicoplast replication occur without cytoplasmic division. Despite its importance in Plasmodium biology, the topoisomerases essential for decatenation of replicated chromosome during endoreduplication remain elusive. We hypothesize that the topoisomerase VI complex, containing Plasmodium falciparum topiosomerase VIB (PfTopoVIB) and catalytic P. falciparum Spo11 (PfSpo11), might be involved in the segregation of the Plasmodium mitochondrial genome. Here, we demonstrate that the putative PfSpo11 is the functional ortholog of yeast Spo11 that can complement the sporulation defects of the yeast Δspo11 strain, and the catalytic mutant Pfspo11Y65F cannot complement such defects. PfTopoVIB and PfSpo11 display a distinct expression pattern compared to the other type II topoisomerases of Plasmodium and are induced specifically at the late schizont stage of the parasite, when the mitochondrial genome segregation occurs. Furthermore, PfTopoVIB and PfSpo11 are physically associated with each other at the late schizont stage, and both subunits are localized in the mitochondria. Using PfTopoVIB- and PfSpo11-specific antibodies, we immunoprecipitated the chromatin of tightly synchronous early, mid-, and late schizont stage-specific parasites and found that both the subunits are associated with the mitochondrial genome during the late schizont stage of the parasite. Furthermore, PfTopoVIB inhibitor radicicol and atovaquone show synergistic interaction. Accordingly, atovaquone-mediated disruption of mitochondrial membrane potential reduces the import and recruitment of both subunits of PfTopoVI to mitochondrial DNA (mtDNA) in a dose-dependent manner. The structural differences between PfTopoVIB and human TopoVIB-like protein could be exploited for development of a novel antimalarial agent. IMPORTANCE This study demonstrates a likely role of topoisomerase VI in the mitochondrial genome segregation of Plasmodium falciparum during endoreduplication. We show that PfTopoVIB and PfSpo11 remain associated and form the functional holoenzyme within the parasite. The spatiotemporal expression of both subunits of PfTopoVI correlates well with their recruitment to the mitochondrial DNA at the late schizont stage of the parasite. Additionally, the synergistic interaction between PfTopoVI inhibitor and the disruptor of mitochondrial membrane potential, atovaquone, supports that topoisomerase VI is the mitochondrial topoisomerase of the malaria parasite. We propose that topoisomerase VI may act as a novel target against malaria.

De novo design of anti-tuberculosis agents using a structure-based deep learning method.

Mycobacterium tuberculosis (Mtb) is a pathogen of major concern due to its ability to withstand both first- and second-line antibiotics, leading to drug resistance. Thus, there is a critical need for identification of novel anti-tuberculosis agents targeting Mtb-specific proteins. The ceaseless search for novel antimicrobial agents to combat drug-resistant bacteria can be accelerated by the development of advanced deep learning methods, to explore both existing and uncharted regions of the chemical space. The adaptation of deep learning methods to under-explored pathogens such as Mtb is a challenging aspect, as most of the existing methods rely on the availability of sufficient target-specific ligand data to design novel small molecules with optimized bioactivity. In this work, we report the design of novel anti-tuberculosis agents targeting the Mtb chorismate mutase protein using a structure-based drug design algorithm. The structure-based deep learning method relies on the knowledge of the target protein's binding site structure alone for conditional generation of novel small molecules. The method eliminates the need for curation of a high-quality target-specific small molecule dataset, which remains a challenge even for many druggable targets, including Mtb chorismate mutase. Novel molecules are proposed, that show high complementarity to the target binding site. The graph attention model could identify the probable key binding site residues, which influenced the conditional molecule generator to design new molecules with pharmacophoric features similar to the known inhibitors.

Ion Selectivity and Permeation Mechanism in a Cyclodextrin-Based Channel.

Synthetic ion channels are a promising technology in the medical and materials sciences because of their ability to conduct ions. Channels based on cyclodextrin, a cyclic oligomer of glucose, are of particular interest because of their nontoxicity and biocompatibility. Using molecular dynamics-based free energy calculations, this study identifies cyclodextrin channel types that are best suited to serve as synthetic ion channels. Free energy profiles show that the connectivity in the channel determines whether the channel is cation-selective or anion-selective. Furthermore, the energy barrier for ion transport is governed by the number of glucose molecules making up the cyclodextrin units of the channel. A detailed mechanism is proposed for ion transport through these channels. Findings from this study will aid in designing cyclodextrin-based channels that could be either cation-selective or anion-selective, by modifying the linkages of the channel or the number of glucose molecules in the cyclodextrin rings.

Elucidation of DNA Repair Function of PfBlm and Potentiation of Artemisinin Action by a Small-Molecule Inhibitor of RecQ Helicase.

Artemisinin (ART)-based combination therapies are recommended as first- and second-line treatments for Plasmodium falciparum malaria. Here, we investigated the impact of the RecQ inhibitor ML216 on the repair of ART-mediated damage in the genome of P. falciparumPfBLM and PfWRN were identified as members of the RecQ helicase family in P. falciparum However, the role of these RecQ helicases in DNA double-strand break (DSB) repair in this parasite has not been explored. Here, we provide several lines of evidence to establish the involvement of PfBlm in DSB repair in P. falciparum First, we demonstrate that PfBlm interacts with two well-characterized DSB repair proteins of this parasite, namely, PfRad51 and PfalMre11. Second, we found that PfBLM expression was upregulated in response to DNA-damaging agents. Third, through yeast complementation studies, we demonstrated that PfBLM could complement the DNA damage sensitivity of a Δsgs1 mutant of Saccharomyces cerevisiae, in contrast to the helicase-dead mutant PfblmK83R Finally, we observe that the overexpression of PfBLM induces resistance to DNA-damaging agents and offers a survival advantage to the parasites. Most importantly, we found that the RecQ inhibitor ML216 inhibits the repair of DSBs and thereby renders parasites more sensitive to ART. Such synergism between ART and ML216 actions was observed for both drug-sensitive and multidrug-resistant strains of P. falciparum Taken together, these findings establish the implications of PfBlm in the Plasmodium DSB repair pathway and provide insights into the antiparasitic activity of the ART-ML216 combination.IMPORTANCE Malaria continues to be a serious threat to humankind not only because of the morbidity and mortality associated with the disease but also due to the huge economic burden that it imparts. Resistance to all available drugs and the unavailability of an effective vaccine cry for an urgent discovery of newer drug targets. Here, we uncovered a role of the PfBlm helicase in Plasmodium DNA double-strand break repair and established that the parasitic DNA repair mechanism can be targeted to curb malaria. The small-molecule inhibitor of PfBlm tested in this study acts synergistically with two first-line malaria drugs, artemisinin (ART) and chloroquine, in both drug-sensitive and multidrug-resistant strains of P. falciparum, thus qualifying this chemical as a potential partner in ART-based combination therapy. Additionally, the identification of this new specific inhibitor of the Plasmodium homologous recombination (HR) mechanism will now allow us to investigate the role of HR in Plasmodium biology.

Human aminolevulinate synthase structure reveals a eukaryotic-specific autoinhibitory loop regulating substrate binding and product release.

5'-aminolevulinate synthase (ALAS) catalyzes the first step in heme biosynthesis, generating 5'-aminolevulinate from glycine and succinyl-CoA. Inherited frameshift indel mutations of human erythroid-specific isozyme ALAS2, within a C-terminal (Ct) extension of its catalytic core that is only present in higher eukaryotes, lead to gain-of-function X-linked protoporphyria (XLP). Here, we report the human ALAS2 crystal structure, revealing that its Ct-extension folds onto the catalytic core, sits atop the active site, and precludes binding of substrate succinyl-CoA. The Ct-extension is therefore an autoinhibitory element that must re-orient during catalysis, as supported by molecular dynamics simulations. Our data explain how Ct deletions in XLP alleviate autoinhibition and increase enzyme activity. Crystallography-based fragment screening reveals a binding hotspot around the Ct-extension, where fragments interfere with the Ct conformational dynamics and inhibit ALAS2 activity. These fragments represent a starting point to develop ALAS2 inhibitors as substrate reduction therapy for porphyria disorders that accumulate toxic heme intermediates.

Selectivity and transport in aquaporins from molecular simulation studies.

The transport of water through aquaporins is a dynamic process that involves rapid movement of a chain of water molecules through the pore of the aquaporin. Structures of aquaporins solved using X-ray crystallography have provided some insights into how water is transported through these channels, and how certain structural features of the pore might help exclude other solutes from passing through the pore. However, such techniques provide only a static picture, and a dynamic picture of the transport and selectivity mechanism at work in aquaporins is possible with molecular dynamics (MD) simulations. In MD simulations, the forces between the different atoms in a system are computed, and the atoms are then allowed to move under the influence of these forces. This allows the sampling of different conformations of the molecule being studied, including conformations that are crucial in driving biological phenomena like water transport. Simulation studies have provided insights into a number of aspects of aquaporins, including the role of the asparagine-proline-alanine (NPA) motif and the aromatic/arginine (ar/R) constriction, water transport mechanism, mechanisms defining the selectivity of the channel, interaction with lipids, response to external electric field, and binding of putative drug molecules. This chapter provides a brief review of the current status of computational modeling of aquaporins using MD simulations. Initially, a brief account of force fields and MD simulations is presented followed by an account of how MD simulations have contributed to further our understanding of different aspects of aquaporins.

TPP riboswitch aptamer: Role of Mg2+ ions, ligand unbinding, and allostery.

Riboswitches are non-coding RNAs that regulate gene expression in response to the binding of metabolites. Their abundance in bacteria makes them ideal drug targets. The prokaryotic thiamine pyrophosphate (TPP) riboswitch regulates gene expression in a wide range of bacteria by undergoing conformational changes in response to the binding of TPP. Although an experimental structure for the aptamer domain of the riboswitch is now available, details of the conformational changes that occur during the binding of the ligand, and the factors that govern these conformational changes, are still not clear. This study employs microsecond-scale molecular dynamics simulations to provide insights into the functioning of the riboswitch aptamer in atomistic detail. A mechanism for the transmission of conformational changes from the ligand-binding site to the P1 switch helix is proposed. Mg2+ ions in the binding site play a critical role in anchoring the ligand to the riboswitch. Finally, modeling the egress of TPP from the binding site reveals a two-step mechanism for TPP unbinding. Findings from this study can motivate the design of future studies aimed at modulating the activity of this drug target.

Cholic Acid-Derived Amphiphile which Combats Gram-Positive Bacteria-Mediated Infections via Disintegration of Lipid Clusters.

Inappropriate and uncontrolled use of antibiotics results in the emergence of antibiotic resistance, thereby threatening the present clinical regimens to treat infectious diseases. Therefore, new antimicrobial agents that can prevent bacteria from developing drug resistance are urgently needed. Selective disruption of bacterial membranes is the most effective strategy for combating microbial infections as accumulation of genetic mutations will not allow for the emergence of drug resistance against these antimicrobials. In this work, we tested cholic acid (CA) derived amphiphiles tethered with different alkyl chains for their ability to combat Gram-positive bacterial infections. In-depth biophysical and biomolecular simulation studies suggested that the amphiphile with a hexyl chain (6) executes more effective interactions with Gram-positive bacterial membranes as compared to other hydrophobic counterparts. Amphiphile 6 is effective against multidrug resistant Gram-positive bacterial strains as well and does not allow the adherence of S. aureus on amphiphile 6 coated catheters implanted in mice. Further, treatment of wound infections with amphiphile 6 clears the bacterial infections. Therefore, the current study presents strategic guidelines in design and development of CA-derived membrane-targeting antimicrobials for Gram-positive bacterial infections.

Microsecond simulation of human aquaporin 2 reveals structural determinants of water permeability and selectivity.

Human aquaporin 2 (AQP2) from the family of aquaporins assumes great physiological importance, owing to its association with nephrogenic diabetes insipidus (NDI). The present study provides detailed insights into the transport properties of AQP2 with the use of microsecond-scale molecular dynamics simulations, and explains how these channels conduct water molecules while at the same time excluding other molecules. Water transport is seen to be diffusion-limited, with a barrier of only 1.6kcalmol-1, and the channel is more water-permeable than other known aquaporins. A constriction site with a pore-facing phenylalanine and arginine is proposed to serve as a selectivity filter as well as a gate modulating the conductance state of the channel. Water molecules form a continuous single-file in the pore lumen, and the orientation of water molecules in this chain is governed by water-protein interactions. A mutant is designed that exhibits different orientation of water molecules, leading to altered permeability. The study complements experimental studies by revealing details of the transport mechanism, energetics, and kinetics. Furthermore, insights obtained into the regulation of permeability in the channel offer the promise of devising new strategies for altering the permeability of the channel under diseased conditions.

Urea Mimics Nucleobases by Preserving the Helical Integrity of B-DNA Duplexes via Hydrogen Bonding and Stacking Interactions.

Urea lesions are formed in DNA because of free radical damage of the thymine base, and their occurrence in DNA blocks DNA polymerases, which has deleterious consequences. Recently, it has been shown that urea is capable of forming hydrogen bonding and stacking interactions with nucleobases, which are responsible for the unfolding of RNA in aqueous urea. Base pairing and stacking are inherent properties of nucleobases; because urea is able to form both, this study attempts to investigate if urea can mimic nucleobases in the context of nucleic acid structures by examining the effect of introducing urea lesions complementary to the four different nucleobases on the overall helical integrity of B-DNA duplexes and their thermodynamic stabilities using molecular dynamics (MD) simulations. The MD simulations resulted in stable duplexes without significant changes in the global B-DNA conformation. The urea lesions occupy intrahelical positions by forming hydrogen bonds with nitrogenous nucleobases, in agreement with experimental results. Furthermore, these urea lesions form hydrogen bonding and stacking interactions with other nucleobases of the same and partner strands, analogous to nucleobases in typical B-DNA duplexes. Direct hydrogen bond interactions are observed for the urea-purine pairs within DNA duplexes, whereas two different modes of pairing, namely, direct hydrogen bonds and water-mediated hydrogen bonds, are observed for the urea-pyrimidine pairs. The latter explains the complexities involved in interpreting the experimental nuclear magnetic resonance data reported previously. Binding free energy calculations were further performed to confirm the thermodynamic stability of the urea-incorporated DNA duplexes with respect to pure duplexes. This study suggests that urea mimics nucleobases by pairing opposite all four nucleobases and maintains the overall structure of the B-DNA duplexes.

Urea-Aromatic Stacking and Concerted Urea Transport: Conserved Mechanisms in Urea Transporters Revealed by Molecular Dynamics.

Urea transporters are membrane proteins that selectively allow urea molecules to pass through. It is not clear how these transporters allow rapid conduction of urea, a polar molecule, in spite of the presence of a hydrophobic constriction lined by aromatic rings. The current study elucidates the mechanism that is responsible for this rapid conduction by performing free energy calculations on the transporter dvUT with a cumulative sampling time of about 1.3 µs. A parallel arrangement of aromatic rings in the pore enables stacking of urea with these rings, which, in turn, lowers the energy barrier for urea transport. Such interaction of the rings with urea is proposed to be a conserved mechanism across all urea-conducting proteins. The free energy landscape for the permeation of multiple urea molecules reveals an interplay between interurea interaction and the solvation state of the urea molecules. This is for the first time that multiple molecule permeation through any small molecule transporter has been modeled.

Cooperation of Hydrophobic Gating, Knock-on Effect, and Ion Binding Determines Ion Selectivity in the p7 Channel.

Ion channels selectively allow certain ions to pass through at much higher rates than others, and thereby modulate ionic concentrations across cell membranes. The current molecular dynamics study elucidates the intricate mechanisms that render ion selectivity to the viral channel p7 by employing free energy calculations. Free energy barriers of 5.4 and 19.4 kcal mol(-1) for K(+) and Ca(2+), respectively, explain the selectivity of the channel reported in experiments. Initially, the permeating ions encounter a hydrophobic barrier followed by stabilization in an ion-binding site. Electrostatic repulsion between the permeating ions propels one of the ions out of the binding site to complete the process of permeation. K(+) and Ca(2+) are seen to exhibit different modes of binding toward a ring of asparagine residues, which serves as the binding site. The findings illustrate how the overall selectivity of a channel can be achieved by a combination of subtle differences.

Ion hydration dynamics in conjunction with a hydrophobic gating mechanism regulates ion permeation in p7 viroporin from hepatitis C virus.

The selectivity of the p7 channel from hepatitis C virus (HCV) toward K(+) over Ca(2+) has made the channel an intriguing system for investigating ion permeation. The present study employs umbrella sampling free energy calculations to investigate the atomistic details of cation conduction through the channel. The free energy profiles suggest that the energy barrier for Ca(2+) conduction is higher than that for K(+) conduction by about 4.5 kcal/mol, thus explaining the selectivity exhibited by the channel toward K(+). A hydrophobic stretch in the channel is proposed to be the primary factor that discriminates K(+) from Ca(2+), and the ion solvation dynamics in this stretch reveals interesting insights into the atomistic mechanisms involved. Two-dimensional free energy landscapes for the ion permeation reveal differences in the lateral motions of K(+) and Ca(2+) with respect to the pore axis, and provide additional details of ion-protein interactions that govern selectivity.

Prediction of the structures of helical membrane proteins based on a minimum unfavorable contacts approach.

An understanding of structure-function relationships of membrane proteins continues to be a challenging problem, owing to the difficulty in obtaining their structures experimentally. This study suggests a method for modeling membrane protein structures that can be used to generate a reliable initial conformation prior to the use of other approaches for sampling conformations. It involves optimizing the orientation of hydrophilic residues so as to minimize unfavorable contacts with the hydrophobic tails of the lipid bilayer. Starting with the optimized initial conformation for three different proteins modeled based on this method, two independent approaches have been used for sampling the conformational space of the proteins. Both approaches are able to predict structures reasonably close to experimental structures, indicating that the initial structure enables the sampling of conformations that are close to the native structure. Possible improvements in the method for making it broadly applicable to helical membrane proteins are discussed.

Modeling the structure of SARS 3a transmembrane protein using a minimum unfavorable contact approach.

3a is an accessory protein from SARS coronavirus that is known to play a significant role in the proliferation of the virus by forming tetrameric ion channels. Although the monomeric units are known to consist of three transmembrane (TM) domains, there are no solved structures available for the complete monomer. The present study proposes a structural model for the transmembrane region of the monomer by employing our previously tested approach, which predicts potential orientations of TM α-helices by minimizing the unfavorable contact surfaces between the different TM domains. The best model structure comprising all three α-helices has been subjected to MD simulations to examine its quality. The TM bundle was found to form a compact and stable structure with significant intermolecular interactions. The structural features of the proposed model of 3a account for observations from previous experimental investigations on the activity of the protein. Further analysis indicates that residues from the TM2 and TM3 domains are likely to line the pore of the ion channel, which is in good agreement with a recent experimental study. In the absence of an experimental structure for the protein, the proposed structure can serve as a useful model for inferring structure-function relationships about the protein. Graphical AbstractThe structure of the membrane protein 3a from SARS coronavirus is modeled using an approach that minimizes unfavorable contacts between transmembrane domains. A structure for a complete monomeric form of the protein thereby proposed is able to account for the behavior of the protein reported in previous experimental studies. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s12039-015-0982-z) contains supplementary material, which is available to authorized users.

Atomistic detailed mechanism and weak cation-conducting activity of HIV-1 Vpu revealed by free energy calculations.

The viral protein U (Vpu) encoded by HIV-1 has been shown to assist in the detachment of virion particles from infected cells. Vpu forms cation-specific ion channels in host cells, and has been proposed as a potential drug target. An understanding of the mechanism of ion transport through Vpu is desirable, but remains limited because of the unavailability of an experimental structure of the channel. Using a structure of the pentameric form of Vpu--modeled and validated based on available experimental data--umbrella sampling molecular dynamics simulations (cumulative simulation time of more than 0.4 µs) were employed to elucidate the energetics and the molecular mechanism of ion transport in Vpu. Free energy profiles corresponding to the permeation of Na+ and K+ were found to be similar to each other indicating lack of ion selection, consistent with previous experimental studies. The Ser23 residue is shown to enhance ion transport via two mechanisms: creating a weak binding site, and increasing the effective hydrophilic length of the channel, both of which have previously been hypothesized in experiments. A two-dimensional free energy landscape has been computed to model multiple ion permeation, based on which a mechanism for ion conduction is proposed. It is shown that only one ion can pass through the channel at a time. This, along with a stretch of hydrophobic residues in the transmembrane domain of Vpu, explains the slow kinetics of ion conduction. The results are consistent with previous conductance studies that showed Vpu to be a weakly conducting ion channel.

Molecular dynamics simulations reveal the HIV-1 Vpu transmembrane protein to form stable pentamers.

The human immunodeficiency virus type I (HIV-1) Vpu protein is 81 residues long and has two cytoplasmic and one transmembrane (TM) helical domains. The TM domain oligomerizes to form a monovalent cation selective ion channel and facilitates viral release from host cells. Exactly how many TM domains oligomerize to form the pore is still not understood, with experimental studies indicating the existence of a variety of oligomerization states. In this study, molecular dynamics (MD) simulations were performed to investigate the propensity of the Vpu TM domain to exist in tetrameric, pentameric, and hexameric forms. Starting with an idealized α-helical representation of the TM domain, a thorough search for the possible orientations of the monomer units within each oligomeric form was carried out using replica-exchange MD simulations in an implicit membrane environment. Extensive simulations in a fully hydrated lipid bilayer environment on representative structures obtained from the above approach showed the pentamer to be the most stable oligomeric state, with interhelical van der Waals interactions being critical for stability of the pentamer. Atomic details of the factors responsible for stable pentamer structures are presented. The structural features of the pentamer models are consistent with existing experimental information on the ion channel activity, existence of a kink around the Ile17, and the location of tetherin binding residues. Ser23 is proposed to play an important role in ion channel activity of Vpu and possibly in virus propagation.