Andrew M. Lynn

g_mmpbsa —A GROMACS Tool for High-Throughput MM-PBSA Calculations

Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), a method to estimate interaction free energies, has been increasingly used in the study of biomolecular interactions. Recently, this method has also been applied as a scoring function in computational drug design. Here a new tool g_mmpbsa, which implements the MM-PBSA approach using subroutines written in-house or sourced from the GROMACS and APBS packages is described. g_mmpbsa was developed as part of the Open Source Drug Discovery (OSDD) consortium. Its aim is to integrate high-throughput molecular dynamics (MD) simulations with binding energy calculations. The tool provides options to select alternative atomic radii and different nonpolar solvation models including models based on the solvent accessible surface area (SASA), solvent accessible volume (SAV), and a model which contains both repulsive (SASA-SAV) and attractive components (described using a Weeks-Chandler-Andersen like integral method). We showcase the effectiveness of the tool by comparing the calculated interaction energy of 37 structurally diverse HIV-1 protease inhibitor complexes with their experimental binding free energies. The effect of varying several combinations of input parameters such as atomic radii, dielectric constant, grid resolution, solute-solvent dielectric boundary definition, and nonpolar models was investigated. g_mmpbsa can also be used to estimate the energy contribution per residue to the binding energy. It has been used to identify those residues in HIV-1 protease that are most critical for binding a range of inhibitors.

Open source drug discovery– A new paradigm of collaborative research in tuberculosis drug development

It is being realized that the traditional closed-door and market driven approaches for drug discovery may not be the best suited model for the diseases of the developing world such as tuberculosis and malaria, because most patients suffering from these diseases have poor paying capacity. To ensure that new drugs are created for patients suffering from these diseases, it is necessary to formulate an alternate paradigm of drug discovery process. The current model constrained by limitations for collaboration and for sharing of resources with confidentiality hampers the opportunities for bringing expertise from diverse fields. These limitations hinder the possibilities of lowering the cost of drug discovery. The Open Source Drug Discovery project initiated by Council of Scientific and Industrial Research, India has adopted an open source model to power wide participation across geographical borders. Open Source Drug Discovery emphasizes integrative science through collaboration, open-sharing, taking up multi-faceted approaches and accruing benefits from advances on different fronts of new drug discovery. Because the open source model is based on community participation, it has the potential to self-sustain continuous development by generating a storehouse of alternatives towards continued pursuit for new drug discovery. Since the inventions are community generated, the new chemical entities developed by Open Source Drug Discovery will be taken up for clinical trial in a non-exclusive manner by participation of multiple companies with majority funding from Open Source Drug Discovery. This will ensure availability of drugs through a lower cost community driven drug discovery process for diseases afflicting people with poor paying capacity. Hopefully what LINUX the World Wide Web have done for the information technology, Open Source Drug Discovery will do for drug discovery.

Rationally Designed Transmembrane Peptide Mimics of the Multidrug Transporter Protein Cdr1 Act as Antagonists to Selectively Block Drug Efflux and Chemosensitize Azole-resistant Clinical Isolates of Candida albicans

Background: Transmembrane peptide mimics (TMPMs) may be useful to prevent helix-helix interaction of multidrug transporter proteins. Results: Rationally designed TMPMs chemosensitized azole-resistant clinical isolates of Candida by blocking drug efflux and in vivo improved the therapeutic efficacy of fluconazole. Conclusion: TMPMs as antagonists offer an alternative approach to chemosensitize azole-resistant Candida isolates. Significance: TMPM-based approach may be extended to other clinically relevant membrane proteins. Drug-resistant pathogenic fungi use several families of membrane-embedded transporters to efflux antifungal drugs from the cells. The efflux pump Cdr1 (Candida drug resistance 1) belongs to the ATP-binding cassette (ABC) superfamily of transporters. Cdr1 is one of the most predominant mechanisms of multidrug resistance in azole-resistant (AR) clinical isolates of Candida albicans. Blocking drug efflux represents an attractive approach to combat the multidrug resistance of this opportunistic human pathogen. In this study, we rationally designed and synthesized transmembrane peptide mimics (TMPMs) of Cdr1 protein (Cdr1p) that correspond to each of the 12 transmembrane helices (TMHs) of the two transmembrane domains of the protein to target the primary structure of the Cdr1p. Several FITC-tagged TMPMs specifically bound to Cdr1p and blocked the efflux of entrapped fluorescent dyes from the AR (Gu5) isolate. These TMPMs did not affect the efflux of entrapped fluorescent dye from cells expressing the Cdr1p homologue Cdr2p or from cells expressing a non-ABC transporter Mdr1p. Notably, the time correlation of single photon counting fluorescence measurements confirmed the specific interaction of FITC-tagged TMPMs with their respective TMH. By using mutant variants of Cdr1p, we show that these TMPM antagonists contain the structural information necessary to target their respective TMHs of Cdr1p and specific binding sites that mediate the interactions between the mimics and its respective helix. Additionally, TMPMs that were devoid of any demonstrable hemolytic, cytotoxic, and antifungal activities chemosensitize AR clinical isolates and demonstrate synergy with drugs that further improved the therapeutic potential of fluconazole in vivo.

Insight into Pleiotropic Drug Resistance ATP-binding Cassette Pump Drug Transport through Mutagenesis of Cdr1p Transmembrane Domains

Background: The Candida albicans efflux pump Cdr1p causes clinically significant antifungal resistance. Results: Biochemical mapping of Cdr1p transmembrane domain mutants reveals residues affecting drug transport. Conclusion: Functional characterization and homology modeling provide insight into the drug binding cavity and substrate promiscuity of Cdr1p. Significance: A platform is provided for systematic Cdr1p structure/function analysis and the rational design of transport modulators/inhibitors. The fungal ATP-binding cassette (ABC) transporter Cdr1 protein (Cdr1p), responsible for clinically significant drug resistance, is composed of two transmembrane domains (TMDs) and two nucleotide binding domains (NBDs). We have probed the nature of the drug binding pocket by performing systematic mutagenesis of the primary sequences of the 12 transmembrane segments (TMSs) found in the TMDs. All mutated proteins were expressed equally well and localized properly at the plasma membrane in the heterologous host Saccharomyces cerevisiae, but some variants differed significantly in efflux activity, substrate specificity, and coupled ATPase activity. Replacement of the majority of the amino acid residues with alanine or glycine yielded neutral mutations, but about 42% of the variants lost resistance to drug efflux substrates completely or selectively. A predicted three-dimensional homology model shows that all the TMSs, apart from TMS4 and TMS10, interact directly with the drug-binding cavity in both the open and closed Cdr1p conformations. However, TMS4 and TMS10 mutations can also induce total or selective drug susceptibility. Functional data and homology modeling assisted identification of critical amino acids within a drug-binding cavity that, upon mutation, abolished resistance to all drugs tested singly or in combinations. The open and closed Cdr1p models enabled the identification of amino acid residues that bordered a drug-binding cavity dominated by hydrophobic residues. The disposition of TMD residues with differential effects on drug binding and transport are consistent with a large polyspecific drug binding pocket in this yeast multidrug transporter.

HMM-ModE – Improved classification using profile hidden Markov models by optimising the discrimination threshold and modifying emission probabilities with negative training sequences

BackgroundProfile Hidden Markov Models (HMM) are statistical representations of protein families derived from patterns of sequence conservation in multiple alignments and have been used in identifying remote homologues with considerable success. These conservation patterns arise from fold specific signals, shared across multiple families, and function specific signals unique to the families. The availability of sequences pre-classified according to their function permits the use of negative training sequences to improve the specificity of the HMM, both by optimizing the threshold cutoff and by modifying emission probabilities to minimize the influence of fold-specific signals. A protocol to generate family specific HMMs is described that first constructs a profile HMM from an alignment of the familys sequences and then uses this model to identify sequences belonging to other classes that score above the default threshold (false positives). Ten-fold cross validation is used to optimise the discrimination threshold score for the model. The advent of fast multiple alignment methods enables the use of the profile alignments to align the true and false positive sequences, and the resulting alignments are used to modify the emission probabilities in the original model.ResultsThe protocol, called HMM-ModE, was validated on a set of sequences belonging to six sub-families of the AGC family of kinases. These sequences have an average sequence similarity of 63% among the group though each sub-group has a different substrate specificity. The optimisation of discrimination threshold, by using negative sequences scored against the model improves specificity in test cases from an average of 21% to 98%. Further discrimination by the HMM after modifying model probabilities using negative training sequences is provided in a few cases, the average specificity rising to 99%. Similar improvements were obtained with a sample of G-Protein coupled receptors sub-classified with respect to their substrate specificity, though the average sequence identity across the sub-families is just 20.6%. The protocol is applied in a high-throughput classification exercise on protein kinases.ConclusionThe protocol has the potential to maximise the contributions of discriminating residues to classify proteins based on their molecular function, using pre-classified positive and negative sequence training data. The high specificity of the method, and increasing availability of pre-classified sequence data holds the potential for its application in sequence annotation.

A key structural domain of the Candida albicans Mdr1 protein

A major multidrug transporter, MDR1 (multidrug resistance 1), a member of the MFS (major facilitator superfamily), invariably contributes to an increased efflux of commonly used azoles and thus corroborates their direct involvement in MDR in Candida albicans. The Mdr1 protein has two transmembrane domains, each comprising six transmembrane helices, interconnected with extracellular loops and ICLs (intracellular loops). The introduction of deletions and insertions through mutagenesis was used to address the role of the largest interdomain ICL3 of the MDR1 protein. Most of the progressive deletants, when overexpressed, eliminated the drug resistance. Notably, restoration of the length of the ICL3 by insertional mutagenesis did not restore the functionality of the protein. Interestingly, most of the insertion and deletion variants of ICL3 became amenable to trypsinization, yielding peptide fragments. The homology model of the Mdr1 protein showed that the molecular surface-charge distribution was perturbed in most of the ICL3 mutant variants. Taken together, these results provide the first evidence that the CCL (central cytoplasmic loop) of the fungal MFS transporter of the DHA1 (drug/proton antiporter) family is critical for the function of MDR. Unlike other homologous proteins, ICL3 has no apparent role in imparting substrate specificity or in the recruitment of the transporter protein.

Rational mutational analysis of a multidrug MFS transporter CaMdr1p of Candida albicans by employing a membrane environment based computational approach.

CaMdr1p is a multidrug MFS transporter of pathogenic Candida albicans. An over-expression of the gene encoding this protein is linked to clinically encountered azole resistance. In-depth knowledge of the structure and function of CaMdr1p is necessary for an effective design of modulators or inhibitors of this efflux transporter. Towards this goal, in this study, we have employed a membrane environment based computational approach to predict the functionally critical residues of CaMdr1p. For this, information theoretic scores which are variants of Relative Entropy (Modified Relative Entropy REM) were calculated from Multiple Sequence Alignment (MSA) by separately considering distinct physico-chemical properties of transmembrane (TM) and inter-TM regions. The residues of CaMdr1p with high REM which were predicted to be significantly important were subjected to site-directed mutational analysis. Interestingly, heterologous host Saccharomyces cerevisiae, over-expressing these mutant variants of CaMdr1p wherein these high REM residues were replaced by either alanine or leucine, demonstrated increased susceptibility to tested drugs. The hypersensitivity to drugs was supported by abrogated substrate efflux mediated by mutant variant proteins and was not attributed to their poor expression or surface localization. Additionally, by employing a distance plot from a 3D deduced model of CaMdr1p, we could also predict the role of these functionally critical residues in maintaining apparent inter-helical interactions to provide the desired fold for the proper functioning of CaMdr1p. Residues predicted to be critical for function across the family were also found to be vital from other previously published studies, implying its wider application to other membrane protein families.

Comparison of theoretical proteomes: Identification of COGs with conserved and variable pI within the multimodal pI distribution

BackgroundTheoretical proteome analysis, generated by plotting theoretical isoelectric points (pI) against molecular masses of all proteins encoded by the genome show a multimodal distribution for pI. This multimodal distribution is an effect of allowed combinations of the charged amino acids, and not due to evolutionary causes. The variation in this distribution can be correlated to the organisms ecological niche. Contributions to this variation maybe mapped to individual proteins by studying the variation in pI of orthologs across microorganism genomes.ResultsThe distribution of ortholog pI values showed trimodal distributions for all prokaryotic genomes analyzed, similar to whole proteome plots. Pairwise analysis of pI variation show that a few COGs are conserved within, but most vary between, the acidic and basic regions of the distribution, while molecular mass is more highly conserved. At the level of functional grouping of orthologs, five groups vary significantly from the population of orthologs, which is attributed to either conservation at the level of sequences or a bias for either positively or negatively charged residues contributing to the function. Individual COGs conserved in both the acidic and basic regions of the trimodal distribution are identified, and orthologs that best represent the variation in levels of the acidic and basic regions are listed.ConclusionThe analysis of pI distribution by using orthologs provides a basis for resolution of theoretical proteome comparison at the level of individual proteins. Orthologs identified that significantly vary between the major acidic and basic regions maybe used as representative of the variation of the entire proteome.

Lipid Induced Conformation of the Tachykinin Peptide Kassinin

Abstract Both the aqueous and lipid-induced structure of Kassinin, a dodecapeptide of amphibian origin, has been studied by two-dimensional proton nuclear magnetic resonance (2D 1H-NMR) spectroscopy and distance geometry calculations. Unambiguous NMR assignments of protons have been made with the aid of correlation spectroscopy (DQF-COSY and TOCSY) experiments and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The distance constraints obtained from the NMR data have been utilized in a distance geometry algorithm to generate a family of structures, which have been refined using restrained energy minimization and dynamics. These data show that, while in water Kassinin prefers to be in an extended chain conformation, in the presence of perdeuterated dodecylphosphocholine (DPC) micelles, a membrane model system, helical conformation is induced in the central core and C-terminal region (K4-M12) of the peptide. N-terminus though less defined also displays some degree of order and a possible turn structure. The conformation adopted by Kassinin in the presence of DPC micelles is consistent with the structural motif typical of neurokinin-1 selective agonists and with that reported for Eledoisin in hydrophobic environment.

What is the Functional Role of N-terminal Transmembrane Helices in the Metabolism Mediated by Liver Microsomal Cytochrome P450 and its Reductase?

We sought to clarify on the hitherto unresolved role of N-terminal transmembrane segments (TMS) of cytochrome P450 (CYP) and its’ reductase (CPR) in protein interaction/catalysis. TMS analyses show little evolutionary conservation in CYPs. The conserved CPR’s TMS poses limited scope for predictable/consistent hetero-recognition with the wide bevy of CYPs’ TMS, as evident from preliminary analyses and TMhit server predictions for inter-helical binding. Further, experimentations with four different CPR preparations (preps) and two liver microsomal CYPs (2C9 and 2E1) shows that the hydroxylated product formation rate is not quantitatively correlated to the extent of integrity of the CPR N-terms. Incorporation of cytochrome b5 in some reactions afforded similar rates while employing either fully intact or partially intact CPR. A survey of literature shows that liver microsomal CYPs function quite well even without the TMS or with significantly altered TMS. These observations negate the hypothesis that N-term TMS of CPR or CYP is obligatory for CYP–CPR interaction and catalysis. Also, in CYP2E1-mediated hydroxylation of para-nitrophenol, the extent of intactness or truncation did not significantly affect the CPR preps’ catalytic role at very low or high substrate concentrations. To interpret these results, we draw support from recently published research on reduced nicotinamide adenide dinucleotide phosphate oxidase (Takac et al., J Biol Chem, 286:13304–13313, 2011) and from our pertinent earlier works. We infer that CPR’ free TMS segment could alter the diffusible reactive oxygen species’ dynamics in the microenvironment, thereby altering the reaction outcome. Based on the evidence, we conclude that TMS merely facilitates “interaction/catalysis” by anchoring the CYP and CPR in the lipid interface.