Protein promiscuity: drug resistance and native functions--HIV-1 case

OBI: A computational tool for the analysis and systematization of the positive selection in proteins

There are multiple tools for positive selection analysis, including vaccine design and detection of variants of circulating drug-resistant pathogens in population selection. However, applying these tools to analyze a large number of protein families or as part of a comprehensive phylogenomics pipeline could be challenging. Since many standard bioinformatics tools are only available as executables, integrating them into complex Bioinformatics pipelines may not be possible.

We have developed OBI, an open-source tool aimed to facilitate positive selection analysis on a large scale. It can be used as a stand-alone command-line app that can be easily installed and used as a Conda package.

Some advantages of using OBI are:•

It speeds up the analysis by automating the entire process

•

It allows multiple starting points and customization for the analysis

•

It allows the retrieval and linkage of structural and evolutive data for a protein through

We hope to provide with OBI a solution for reliably speeding up large-scale protein evolutionary and structural analysis.

Substrate-Assisted Catalysis in the Reaction Catalyzed by Salicylic Acid Binding Protein 2 (SABP2), a Potential Mechanism of Substrate Discrimination for Some Promiscuous Enzymes

Although one of an enzyme's hallmarks is the high specificity for their natural substrates, substrate promiscuity has been reported more frequently. It is known that promiscuous enzymes generally show different catalytic efficiencies to different substrates, but our understanding of the origin of such differences is still lacking. Here we report the results of quantum mechanical/molecular mechanical simulations and an experimental study of salicylic acid binding protein 2 (SABP2). SABP2 has promiscuous esterase activity toward a series of substrates but shows a high activity toward its natural substrate, methyl salicylate (MeSA). Our results demonstrate that this enzyme may use substrate-assisted catalysis involving the hydroxyl group from MeSA to enhance the activity and achieve substrate discrimination.

Molecular dynamics studies on both bound and unbound renin protease

The aspartic protease renin (REN) catalyses the rate-limiting step in the Renin-Angiotensin-Aldosterone System (RAAS), which regulates cardiovascular and renal homoeostasis in living organisms. Renin blockage is therefore an attractive therapeutic strategy for the treatment of hypertension. Herein, computational approaches were used to provide a structural characterization of the binding site, flap opening and dynamic rearrangements of REN in the key conserved residues and water molecules, with the binding of a dodecapeptide substrate or different inhibitors. All these structural insights during catalysis may assist future studies in developing novel strategies for REN inactivation. Our molecular dynamics simulations of several unbound-REN and bound-REN systems indicate similar flexible-segments plasticity with larger fluctuations in those belonging to the C-domain (exposed to the solvent). These segments are thought to assist the flap opening and closure to allow the binding of the substrate and catalytic water molecules. The unbound-REN simulation suggests that the flap can acquire three different conformations: closed, semi-open and open. Our results indicate that the semi-open conformation is already sufficient and appropriate for the binding of the angiotensinogen (Ang) tail, thus contributing to the high specificity of REN, and that both semi-open and open flap conformations are present in free and complexed enzymes. We additionally observed that the Tyr75-Trp39 H-bond has an important role in assisting flap movement, and we highlight several conserved water molecules and amino acids that are essential for the proper catalytic activity of REN.

A measure of the promiscuity of proteins and characteristics of residues in the vicinity of the catalytic site that regulate promiscuity

Promiscuity, the basis for the evolution of new functions through 'tinkering' of residues in the vicinity of the catalytic site, is yet to be quantitatively defined. We present a computational method Promiscuity Indices Estimator (PROMISE)--based on signatures derived from the spatial and electrostatic properties of the catalytic residues, to estimate the promiscuity (PromIndex) of proteins with known active site residues and 3D structure. PromIndex reflects the number of different active site signatures that have congruent matches in close proximity of its native catalytic site, the quality of the matches and difference in the enzymatic activity. Promiscuity in proteins is observed to follow a lognormal distribution (μ = 0.28, σ = 1.1 reduced chi-square = 3.0E-5). The PROMISE predicted promiscuous functions in any protein can serve as the starting point for directed evolution experiments. PROMISE ranks carboxypeptidase A and ribonuclease A amongst the more promiscuous proteins. We have also investigated the properties of the residues in the vicinity of the catalytic site that regulates its promiscuity. Linear regression establishes a weak correlation (R(2)∼0.1) between certain properties of the residues (charge, polar, etc) in the neighborhood of the catalytic residues and PromIndex. A stronger relationship states that most proteins with high promiscuity have high percentages of charged and polar residues within a radius of 3 Å of the catalytic site, which is validated using one-tailed hypothesis tests (P-values∼0.05). Since it is known that these characteristics are key factors in catalysis, their relationship with the promiscuity index cross validates the methodology of PROMISE.

Probing the sequence space available for HIV-1 evolution

We designed a novel experimental approach to probe the sequence space available for HIV-1 evolution. Selective pressure was put on conserved HIV-1 genomic sequences by means of RNA interference (RNAi). Virus escape was monitored in many parallel cultures, and we scored the mutations selected in the RNAi target sequences. The experimentally induced sequence variation closely resembles the sequence variation of natural HIV-1 strains. This indicates that we actually mapped a restricted area of sequence space compatible with virus replication.

Predicting the complex structure and functional motions of the outer membrane transporter and signal transducer FecA

Escherichia coli requires an efficient transport and signaling system to successfully sequester iron from its environment. FecA, a TonB-dependent protein, serves a critical role in this process: first, it binds and transports iron in the form of ferric citrate, and second, it initiates a signaling cascade that results in the transcription of several iron transporter genes in interaction with inner membrane proteins. The structure of the plug and barrel domains and the periplasmic N-terminal domain (NTD) are separately available. However, the linker connecting the plug and barrel and the NTD domains is highly mobile, which may prevent the determination of the FecA structure as a whole assembly. Here, we reduce the conformation space of this linker into most probable structural models using the modeling tool CABS, then apply normal-mode analysis to investigate the motions of the whole structure of FecA by using elastic network models. We relate the FecA domain motions to the outer-inner membrane communication, which initiates transcription. We observe that the global motions of FecA assign flexibility to the TonB box and the NTD, and control the exposure of the TonB box for binding to the TonB inner membrane protein, suggesting how these motions relate to FecA function. Our simulations suggest the presence of a communication between the loops on both ends of the protein, a signaling mechanism by which a signal could be transmitted by conformational transitions in response to the binding of ferric citrate.

Latent evolutionary potentials under the neutral mutational drift of an enzyme

Biological systems exhibit mutational robustness, or neutrality, whereby the impact of mutations is minimized. Does neutrality hamper their ability to adapt in the face of changing environments? We monitored changes in genotype and phenotype that occur within a neutral mutational network of an enzyme, experimentally and computationally (see accompanying article). Using the enzyme PON1 as a model, we performed random mutagenesis and purifying selection to purge deleterious mutations. We characterized approximately 300 variants that are apparently neutral, or close to neutral, with respect to PON1's levels of expression and native lactonase activity. Their activities with promiscuous substrates and ligands indicated significant changes in adaptive potentials. Almost half of the variants exhibited changes in promiscuous activities, specificities, or inhibition, and several of these were found to be one or two mutations, closer to potentially new phenotypes: aryl esterase, thiolactonase, phosphotriesterase, or drug resistance. This empirical measure of phenotypic changes under neutrality supports the notion that sequence changes that are neutral, i.e., non-adaptive, in a current context can facilitate adaptation under changing circumstances, by both expanding the activity range of existing enzymes and thus providing an immediate advantage, and by reducing the number of mutations required for divergence of new functions.

Functional clustering of yeast proteins from the protein-protein interaction network

BACKGROUND

The abundant data available for protein interaction networks have not yet been fully understood. New types of analyses are needed to reveal organizational principles of these networks to investigate the details of functional and regulatory clusters of proteins.

RESULTS

In the present work, individual clusters identified by an eigenmode analysis of the connectivity matrix of the protein-protein interaction network in yeast are investigated for possible functional relationships among the members of the cluster. With our functional clustering we have successfully predicted several new protein-protein interactions that indeed have been reported recently.

CONCLUSION

Eigenmode analysis of the entire connectivity matrix yields both a global and a detailed view of the network. We have shown that the eigenmode clustering not only is guided by the number of proteins with which each protein interacts, but also leads to functional clustering that can be applied to predict new protein interactions.

Comparative pharmacogenomics of antiretroviral and cytotoxic treatments

Few genetic markers are used routinely to predict clinical effectiveness and toxic effects despite the fact that physicians and their patients are consistently confronted with this balance. Because one of the goals of pharmacogenomics is to identify individuals and target populations that might have adverse outcomes, pharmaceutical companies have been reluctant to use a strategy that might identify patients who are not eligible for a particular treatment. This view is changing because drug-discovery programmes and treatments that target specific pathways, are showing improvements in surrogate and survival endpoints. HIV and cancer are now regarded as chronic diseases, which commonly need life-long systemic treatment from the time of diagnosis. HIV and cancer medicine have used pharmacogenomics to some extent in clinical care. Common and classic features of pharmacogenomics that are related to both antiretroviral treatment and to cytotoxic treatment are discussed in this review, providing a framework for individual treatment of these diseases.

Protein promiscuity: drug resistance and native functions--HIV-1 case

The association of a drug with its target protein has the effect of blocking the protein activity and is termed a promiscuous function to distinguish from the protein's native function (Tawfik and associates, Nat. Genet. 37, 73-6, 2005). Obviously, a protein has not evolved naturally for drug association or drug resistance. Promiscuous protein functions exhibit unique traits of evolutionary adaptability, or evolvability, which is dependent on the induction of novel phenotypic traits by a small number of mutations. These mutations might have small effects on native functions, but large effects on promiscuous function; for example, an evolving protein could become increasingly drug resistant while maintaining its original function. Ariel Fernandez, in his opinion piece, notes that drug-binding "promiscuity" can hardly be dissociated from native functions; a dominant approach to drug discovery is the protein-native-substrate transition-state mimetic strategy. Thus, man-made ligands (e.g. drugs) have been successfully crafted to restrain enzymatic activity by focusing on the very same structural features that determine the native function. Using the successful inhibition of HIV-1 protease as an example, Fernandez illustrates how drug designers have employed naturally evolved features of the protein to suppress its activity. Based on these arguments, he dismisses the notion that drug binding is quintessentially promiscuous, even though in principle, proteins did not evolve to associate with man made ligands. In short, Fernandez argues that there may not be separate protein domains that one could term promiscuous domains. While acknowledging that drugs may bind promiscuously or in a native-like manner a la Fernandez, Tawfik maintains the role of evolutionary adaptation, even when a drug binds native-like. In the case of HIV-1 protease, drugs bind natively, and the initial onset of mutations results in drug resistance in addition to a dramatic decline in enzymatic activity and fitness of the virus. A chain of compensatory mutations follows this, and then the virus becomes fully fit and drug resistant. Ben Berkhout and Rogier Sanders subscribe to the evolution of new protein functions through gene duplication. With two identical protein domains, one domain can be released from a constraint imposed by the original function and it is thus free to move in sequence space toward a new function without loss of the original function. They emphasize that the forced evolution of drug-resistance differs significantly from the spontaneous evolution of an additional protein function. For instance, the latter process could proceed gradually on an evolutionary time scale, whereas the acquisition of drug-resistance is an all or nothing process for a virus, leading to the failure or success of therapy. They find no evidence to the thesis that resistance-mutations appear more rapidly in promiscuous domains than native domains. Berkhout and Sanders illustrate the genetic plasticity of HIV-1 by citing examples in which well-conserved amino acid residues of catalytic domains are forced to mutate under drug-pressure. HIV drug resistance biology is very complex. Instead of a viral protein, a drug can be targeted at a cellular protein. For example, Berkhout and Sanders claim, a drug targeted at the cellular protein CCR5 inhibits the binding of the viral envelope glycoprotein (Env) to CCR5. However, Env mutates so that it binds to the CCR5-drug complex and develops drug resistance. Interestingly, CCR5 has not evolved to bind to Env, but to a series of chemokines. Andrzej Kloczkowski, Taner Sen, and Bob Jernigan point out the importance of protein motions for binding. They believe it is likely that different ligands can bind to the diverse protein conformations sampled in the course of normal protein conformational fluctuations. They have been applying simple elastic network models to extract the motions as normal modes, which yield relatively small numbers of conformations that are useful for developing protein mechanisms; while these are typically small motions, for some proteins they can be quite large in scale. One of the major advantages of the approach is that only relatively small numbers of modes are important contributors to the overall motion -- so the approach provides a way to systematically map out a protein's motions. These models successfully represent the conformational fluctuations manifested in the crystallographic B-factors, and often suggest motions related to protein functional behaviors, such as those observed for reverse transcriptase, where two dominant hinges clearly relate to the processing steps -- one showing anti-correlation between the polymerase and ribonuclease H sites related to the translation and positioning of the nucleic acid chain, and another for opening and closing the polymerase site. Disordered proteins represent a more extreme case where the set of accessible conformations is much larger; thus they could offer up a broader range of possible binding forms. Whether evolution controls the functional motions for proteins remains little studied. Intriguingly, buried in the existing databases of protein-protein interactions may be information that can shed light on the extent of promiscuous binding among proteins themselves. Within these data there are cases where large numbers of diverse proteins have been shown to interact with a single protein; some of these could represent promiscuous protein-protein binding. Uncovering these promiscuous behaviors could be important for comprehending the details of how proteins can bind promiscuously to one another, and can exhibit even greater promiscuity in their binding to small molecules. The evolutionary routes, the dynamics of the target protein, and the many other aspects that need to be addressed while designing a drug that may dodge drug resistance, indicate the complexity and multi-disciplinary nature of the issue of drug resistance.