Multidrug-resistant (MDR) gram-negative bacteria have increased the prevalence of fatal sepsis in modern times. Colistin is a cationic antimicrobial peptide (CAMP) antibiotic that permeabilizes the bacterial outer membrane (OM) and has been used to treat these infections. The OM outer leaflet is comprised of endotoxin containing lipid A, which can be modified to increase resistance to CAMPs and prevent clearance by the innate immune response. One type of lipid A modification involves the addition of phosphoethanolamine to the 1 and 4′ headgroup positions by phosphoethanolamine transferases. Previous structural work on a truncated form of this enzyme suggested that the full-length protein was required for correct lipid substrate binding and catalysis. We now report the crystal structure of a full-length lipid A phosphoethanolamine transferase from Neisseria meningitidis, determined to 2.75-Å resolution. The structure reveals a previously uncharacterized helical membrane domain and a periplasmic facing soluble domain. The domains are linked by a helix that runs along the membrane surface interacting with the phospholipid head groups. Two helices located in a periplasmic loop between two transmembrane helices contain conserved charged residues and are implicated in substrate binding. Intrinsic fluorescence, limited proteolysis, and molecular dynamics studies suggest the protein may sample different conformational states to enable the binding of two very different-sized lipid substrates. These results provide insights into the mechanism of endotoxin modification and will aid a structure-guided rational drug design approach to treating multidrug-resistant bacterial infections.lipid modification | multidrug resistance | molecular dynamics | Neisseria | membrane protein structure

The multidrug transporter P-glycoprotein (P-gp) is central to the development of multidrug resistance in cancer. While residues essential for transport and binding have been identified, the location, composition, and specificity of potential drug binding sites are uncertain. Here molecular dynamics simulations are used to calculate the free energy profile for the binding of morphine and nicardipine to P-gp. We show that morphine and nicardipine primarily interact with key residues implicated in binding and transport from mutational studies, binding at different but overlapping sites within the transmembrane pore. Their permeation pathways were distinct but involved overlapping sets of residues. The results indicate that the binding location and permeation pathways of morphine and nicardipine are not well separated and cannot be considered as unique. This has important implications for our understanding of substrate uptake and transport by P-gp. Our results are independent of the choice of starting structure and consistent with a range of experimental studies.

Despite decades of research, the mechanism of action of the ABC multidrug transporter P-glycoprotein (P-gp) remains elusive. Due to experimental limitations, many researchers have turned to molecular dynamics simulation studies in order to investigate different aspects of P-gp function. However, such studies are challenging and caution is required when interpreting the results. P-gp is highly flexible and the time scale on which it can be simulated is limited. There is also uncertainty regarding the accuracy of the various crystal structures available, let alone the structure of the protein in a physiologically relevant environment. In this study, three alternative structural models of mouse P-gp (3G5U, 4KSB, 4M1M), all resolved to 3.8 Å, were used to initiate sets of simulations of P-gp in a membrane environment in order to determine: a) the sensitivity of the results to differences in the starting configuration; and b) the extent to which converged results could be expected on the times scales commonly simulated for this system. The simulations suggest that the arrangement of the nucleotide binding domains (NBDs) observed in the crystal structures is not stable in a membrane environment. In all simulations, the NBDs rapidly associated (within 10 ns) and changes within the transmembrane helices were observed. The secondary structure within the transmembrane domain was best preserved in the 4M1M model under the simulation conditions used. However, the extent to which replicate simulations diverged on a 100 to 200 ns timescale meant that it was not possible to draw definitive conclusions as to which structure overall was most stable, or to obtain converged and reliable results for any of the properties examined. The work brings into question the reliability of conclusions made in regard to the nature of specific interactions inferred from previous simulation studies on this system involving similar sampling times. It also highlights the need to demonstrate the statistical significance of any results obtained in simulations of large flexible proteins, especially where the initial structure is uncertain.

A family of flavin/deazaflavin-dependent oxidoreductases (FDORs) from mycobacteria has been recently characterized and found to play a variety of catalytic roles, including the activation of prodrugs such as the candidate anti-tuberculosis drug pretomanid (PA-824). However, our understanding of the catalytic mechanism used by these enzymes is relatively limited. To address this, we have used a combination of quantum mechanics and molecular dynamics calculations to study the catalytic mechanism of the activation of pretomanid by the deazaflavin-dependent nitroreductase (Ddn) from Mycobacterium tuberculosis. The preferred pathway involves an initial hydride transfer step from the deprotonated cofactor (i.e., FH), with subsequent protonation, before a series of spontaneous intramolecular reactions to form the final reactive nitrogen species. The most likely proton source is a hydroxonium ion within the solvent accessible active site. Intriguingly, catalysis of the rate-determining hydride transfer step is aided by three tyrosine residues that form a hydrophobic barrier around the active site that, upon reaction, is then disrupted to allow increased water accessibility to facilitate the subsequent proton transfer step. The catalytic mechanism we propose is consistent with previous experimental observations of the Ddn enzyme and will inform the design of improved prodrugs in the future.

Bilirubin is a potent antioxidant that is produced from the reduction of the heme degradation product biliverdin. In mammalian cells and Cyanobacteria, NADH/NADPH-dependent biliverdin reductases (BVRs) of the Rossmann-fold have been shown to catalyze this reaction. Here, we describe the characterization of Rv2074 from Mycobacterium tuberculosis, which belongs to a structurally and mechanistically distinct family of F 420 H 2 -dependent BVRs (F-BVRs) that are exclusively found in Actinobacteria. We have solved the crystal structure of Rv2074 bound to its cofactor, F 420 , and used this alongside molecular dynamics simulations, site-directed mutagenesis and NMR spectroscopy to elucidate its catalytic mechanism. The production of bilirubin by Rv2074 could exploit the anti-oxidative properties of bilirubin and contribute to the range of immunoevasive mechanisms that have evolved in M. tuberculosis to allow persistent infection.

The utility of enzymes as industrial catalysts is reduced by their low stability and limited operating range. Connal and co-workers have developed a simple, onestep procedure to afford an enzyme-inspired catalyst that contains the common functional motif present at the active site of many hydrolases-the catalytic triadincorporated into a single unit. This enzyme-inspired catalyst is then easily immobilized onto a resin support and hydrophobically tuned to perform enzymelike hydrolytic catalysis. Computational modeling provides clues to a concerted, two-step mechanism analogous to native enzymes.

The active sites of the (6-4) photolyases contain two conserved histidine residues, which, in the Drosophila melanogaster enzyme, correspond to His365 and His369. While there are nine combinations in which the three possible protonation states of the two histidines (with protons on Nδ (HID), Nε (HIE), or both Nδ and Nε (HIP)) can be paired, there is presently no consensus as to which of these states is present, let alone mechanistically relevant. EPR hyperfine couplings for selected protons of the FADH(•) radical have previously been used to address this issue. Our QM/MM calculations show, however, that the experimental couplings are equally well reproduced by each of the nine combinations. Since the EPR results seemingly cannot be used to unequivocally assign the protonation states, the pKa values of the two histidines were calculated using the popular PROPKA, H++, and APBS approaches, in various environments and for several lesions. These techniques consistently indicate that, at pH = 7, both His365 and His369 should be neutral, although His369 is found to be more prone to becoming protonated. In a comparative approach, a series of molecular dynamics simulations was performed with all nine combinations, employing various reference crystal structures and different oxidation states of the FAD cofactor. The overall result of this approach is in agreement with our pKa results. Consequently, although the introduction of the reduced cofactor results in an increased stability for selected protonated states, particularly the His365═HID and His369═HIP combination, the neutral combination His365═HID and His365═HIE stands out as the most relevant state for the activity of the enzyme.

We have investigated the conformational phase spaces of both Met-enkephalin and Ada-enkephalin in 2,2,2-trifluoroethanol in order to connect them to their respective CD spectra. To this end, we have characterized the conformational preferences of the zwitterionic and neutral forms of Met-enkephalin and of both the R- and the S-epimers of Ada-enkephalin, as obtained by classical molecular dynamics. The CD spectrum for each peptide was subsequently obtained with a procedure of successive averaging, which accounts for the behavior of the solvent, the side chains, and the backbone variations of the peptides. To make an appropriate comparison with experimental results, we have produced composite spectra that account for the appropriate contributions of the zwitterionic and neutral forms of the peptides as well as the expected epimeric ratio. Such a procedure results in theoretically obtained CD spectra that show significant promise in terms of reproducing their experimentally measured counterparts.