In Vitro Analysis of PDZ-dependent CFTR Macromolecular Signaling Complexes

Cystic fibrosis transmembrane conductance regulator (CFTR), a chloride channel located primarily at the apical membranes of epithelial cells, plays a crucial role in transepithelial fluid homeostasis^1-3. CFTR has been implicated in two major diseases: cystic fibrosis (CF)⁴ and secretory diarrhea⁵. In CF, the synthesis or functional activity of the CFTR Cl- channel is reduced. This disorder affects approximately 1 in 2,500 Caucasians in the United States⁶. Excessive CFTR activity has also been implicated in cases of toxin-induced secretory diarrhea (e.g., by cholera toxin and heat stable E. coli enterotoxin) that stimulates cAMP or cGMP production in the gut⁷.Accumulating evidence suggest the existence of physical and functional interactions between CFTR and a growing number of other proteins, including transporters, ion channels, receptors, kinases, phosphatases, signaling molecules, and cytoskeletal elements, and these interactions between CFTR and its binding proteins have been shown to be critically involved in regulating CFTR-mediated transepithelial ion transport in vitro and also in vivo^8-19. In this protocol, we focus only on the methods that aid in the study of the interactions between CFTR carboxyl terminal tail, which possesses a protein-binding motif [referred to as PSD95/Dlg1/ZO-1 (PDZ) motif], and a group of scaffold proteins, which contain a specific binding module referred to as PDZ domains. So far, several different PDZ scaffold proteins have been reported to bind to the carboxyl terminal tail of CFTR with various affinities, such as NHERF1, NHERF2, PDZK1, PDZK2, CAL (CFTR-associated ligand), Shank2, and GRASP^20-27. The PDZ motif within CFTR that is recognized by PDZ scaffold proteins is the last four amino acids at the C terminus (i.e., 1477-DTRL-1480 in human CFTR)²⁰. Interestingly, CFTR can bind more than one PDZ domain of both NHERFs and PDZK1, albeit with varying affinities²². This multivalency with respect to CFTR binding has been shown to be of functional significance, suggesting that PDZ scaffold proteins may facilitate formation of CFTR macromolecular signaling complexes for specific/selective and efficient signaling in cells^16-18.Multiple biochemical assays have been developed to study CFTR-involving protein interactions, such as co-immunoprecipitation, pull-down assay, pair-wise binding assay, colorimetric pair-wise binding assay, and macromolecular complex assembly assay^16-19,28,29. Here we focus on the detailed procedures of assembling a PDZ motif-dependent CFTR-containing macromolecular complex in vitro, which is used extensively by our laboratory to study protein-protein or domain-domain interactions involving CFTR^16-19,28,29.     

In Vitro HIV-1 Evolution in Response to Triple Reverse Transcriptase Inhibitors & In Silico Phenotypic Analysis

Background

Effectiveness of ART regimens strongly depends upon complex interactions between the selective pressure of drugs and the evolution of mutations that allow or restrict drug resistance.

Methods

Four clinical isolates from NRTI-exposed, NNRTI-naive subjects were passaged in increasing concentrations of NVP in combination with 1 µM 3 TC and 2 µM ADV to assess selective pressures of multi-drug treatment. A novel parameter inference procedure, based on a stochastic viral growth model, was used to estimate phenotypic resistance and fitness from in vitro combination passage experiments.

Results

Newly developed mathematical methods estimated key phenotypic parameters of mutations arising through selective pressure exerted by 3 TC and NVP. Concentrations of 1 µM 3 TC maintained the M184V mutation, which was associated with intrinsic fitness deficits. Increasing NVP concentrations selected major NNRTI resistance mutations. The evolutionary pathway of NVP resistance was highly dependent on the viral genetic background, epistasis as well as stochasticity. Parameter estimation indicated that the previously unrecognized mutation L228Q was associated with NVP resistance in some isolates.

Conclusion

Serial passage of viruses in the presence of multiple drugs may resemble the selection of mutations observed among treated individuals and populations in vivo and indicate evolutionary preferences and restrictions. Phenotypic resistance estimated here “in silico” from in vitro passage experiments agreed well with previous knowledge, suggesting that the unique combination of “wet-” and “dry-lab” experimentation may improve our understanding of HIV-1 resistance evolution in the future.     

In Silico Analysis Reveals Sequential Interactions and Protein Conformational Changes during the Binding of Chemokine CXCL-8 to Its Receptor CXCR1

Chemokine CXCL-8 plays a central role in human immune response by binding to and activate its cognate receptor CXCR1, a member of the G-protein coupled receptor (GPCR) family. The full-length structure of CXCR1 is modeled by combining the structures of previous NMR experiments with those from homology modeling. Molecular docking is performed to search favorable binding sites of monomeric and dimeric CXCL-8 with CXCR1 and a mutated form of it. The receptor-ligand complex is embedded into a lipid bilayer and used in multi ns molecular dynamics (MD) simulations. A multi-steps binding mode is proposed: (i) the N-loop of CXCL-8 initially binds to the N-terminal domain of receptor CXCR1 driven predominantly by electrostatic interactions; (ii) hydrophobic interactions allow the N-terminal Glu-Leu-Arg (ELR) motif of CXCL-8 to move closer to the extracellular loops of CXCR1; (iii) electrostatic interactions finally dominate the interaction between the N-terminal ELR motif of CXCL-8 and the EC-loops of CXCR1. Mutation of CXCR1 abrogates this mode of binding. The detailed binding process may help to facilitate the discovery of agonists and antagonists for rational drug design.     

In Silico Prediction and In Vivo Validation of Daphnia pulex Micrornas

Daphnia pulex, the crustacean with the first sequenced genome, is an important organism that has been widely used in ecological and toxicological research. MicroRNAs (miRNAs) are 21–25 nucleotide small non-coding RNAs that are involved in a myriad of physiological processes. In this research, we predicted 75 D. pulex miRNAs by sequence homology and secondary structure identification from the full genome sequence. Fourteen predicted miRNAs were selected for quantitative real time polymerase chain reaction (RT-PCR) validation. Out of these, eight (mir-8, mir-9, mir-12, mir-92, mir-100, mir-133, mir-153 and mir-283) were successfully amplified and validated. Next, expression levels were quantified at three different life stages (days 4, 8 and 12 of age) using U6 spliceosomal RNA as a reference gene. The expression of mir-8, mir-9, mir-12, mir-92 and mir-100 significantly differed across time suggesting these microRNAs might play a critical role during D. pulex development. This is the first study to identify and validate miRNAs in D. pulex, which is an important first step in further studies that evaluate their roles in development and response to environmental and ecological stimuli.     

In Silico Analysis of Putative Paralytic Shellfish Poisoning Toxins Export Proteins in Cyanobacteria

Paralytic shellfish poisoning toxins (PSTs) are a family of more than 30 natural alkaloids synthesized by dinoflagellates and cyanobacteria whose toxicity in animals is mediated by voltage-gated Na⁺ channel blocking. The export of PST analogues may be through SxtF and SxtM, two putative MATE (multidrug and toxic compound extrusion) family transporters encoded in PSTs biosynthetic gene cluster (sxt). sxtM is present in every sxt cluster analyzed; however, sxtF is only present in the Cylindrospermopsis-Raphidiopsis clade. These transporters are energetically coupled with an electrochemical gradient of proton (H⁺) or sodium (Na⁺) ions across membranes. Because the functional role of PSTs remains unknown and methods for genetic manipulation in PST-producing organisms have not yet been developed, protein structure analyses will allow us to understand their function. By analyzing the sxt cluster of eight PST-producing cyanobacteria, we found no correlation between the presence of sxtF or sxtM and a specific PSTs profile. Phylogenetic analyses of SxtF/M showed a high conservation of SxtF in the Cylindrospermopsis-Raphidiopsis clade, suggesting conserved substrate affinity. Two domains involved in Na⁺ and drug recognition from NorM proteins (MATE family) of Vibrio parahaemolyticus and V. cholerae are present in SxtF/M. The Na⁺ recognition domain was conserved in both SxtF/M, indicating that Na⁺ can maintain the role as a cation anti-transporter. Consensus motifs for toxin binding differed between SxtF and SxtM implying differential substrate binding. Through protein modeling and docking analysis, we found that there is no marked affinity between the recognition domain and a specific PST analogue. This agrees with our previous results of PST export in R. brookii D9, where we observed that the response to Na⁺ incubation was similar to different analogues. These results reassert the hypothesis regarding the involvement of Na⁺ in toxin export, as well as the motifs L³⁹⁸XGLQD⁴⁰³ (SxtM) and L³⁹⁰VGLRD³⁹⁵ (SxtF) in toxin recognition.     

In Silico Analysis of Functional Single Nucleotide Polymorphisms in the Human TRIM22 Gene

Tripartite motif protein 22 (TRIM22) is an evolutionarily ancient protein that plays an integral role in the host innate immune response to viruses. The antiviral TRIM22 protein has been shown to inhibit the replication of a number of viruses, including HIV-1, hepatitis B, and influenza A. TRIM22 expression has also been associated with multiple sclerosis, cancer, and autoimmune disease. In this study, multiple in silico computational methods were used to identify non-synonymous or amino acid-changing SNPs (nsSNP) that are deleterious to TRIM22 structure and/or function. A sequence homology-based approach was adopted for screening nsSNPs in TRIM22, including six different in silico prediction algorithms and evolutionary conservation data from the ConSurf web server. In total, 14 high-risk nsSNPs were identified in TRIM22, most of which are located in a protein interaction module called the B30.2 domain. Additionally, 9 of the top high-risk nsSNPs altered the putative structure of TRIM22''s B30.2 domain, particularly in the surface-exposed v2 and v3 regions. These same regions are critical for retroviral restriction by the closely-related TRIM5α protein. A number of putative structural and functional residues, including several sites that undergo post-translational modification, were also identified in TRIM22. This study is the first extensive in silico analysis of the highly polymorphic TRIM22 gene and will be a valuable resource for future targeted mechanistic and population-based studies.     

In Silico Classification of Proteins from Acidic and Neutral Cytoplasms

Protein acidostability is a common problem in biopharmaceutical and other industries. However, it remains a great challenge to engineer proteins for enhanced acidostability because our knowledge of protein acidostabilization is still very limited. In this paper, we present a comparative study of proteins from bacteria with acidic (AP) and neutral cytoplasms (NP) using an integrated statistical and machine learning approach. We construct a set of 393 non-redundant AP-NP ortholog pairs and calculate a total of 889 sequence based features for these proteins. The pairwise alignments of these ortholog pairs are used to build a residue substitution propensity matrix between APs and NPs. We use Gini importance provided by the Random Forest algorithm to rank the relative importance of these features. A scoring function using the 10 most significant features is developed and optimized using a hill climbing algorithm. The accuracy of the score function is 86.01% in predicting AP-NP ortholog pairs and is 76.65% in predicting non-ortholog AP-NP pairs, suggesting that there are significant differences between APs and NPs which can be used to predict relative acidostability of proteins. The overall trends uncovered in the study can be used as general guidelines for designing acidostable proteins. To best of our knowledge, this work represents the first systematic comparative study of the acidostable proteins and their non-acidostable orthologs.     

Effect of Starvation on the Endocytic Pathway in Dictyostelium Cells

Dictyostelium discoideum amoebae have been used extensively to study the structure and dynamics of the endocytic pathway. Here, we show that while the general structure of the endocytic pathway is maintained in starved cells, its dynamics rapidly slow down. In addition, analysis of apm3 and lvsB mutants reveals that the functional organization of the endocytic pathway is profoundly modified upon starvation. Indeed, in these mutant cells, some of the defects observed in rich medium persist in starved cells, notably an abnormally slow transfer of endocytosed material between endocytic compartments. Other parameters, such as endocytosis of the fluid phase or the rate of fusion of postlysosomes to the cell surface, vary dramatically upon starvation. Studying the endocytic pathway in starved cells can provide a different perspective, allowing the primary (invariant) defects resulting from specific mutations to be distinguished from their secondary (conditional) consequences.Dictyostelium discoideum is a widely used model organism for studying the organization and function of the endocytic pathway. In Dictyostelium, the organization of the endocytic pathway is similar to that in higher eukaryotes. The pathway in Dictyostelium can be divided into four steps (see Fig. S1 in the supplemental material): uptake at the plasma membrane of particles and medium, transfer through early acidic endocytic compartments (lysosomes), passage into less acidic postlysosomes (PLs), and finally, exocytosis of undigested materials (17, 20). Thus, Dictyostelium recapitulates many of the functions of the endocytic pathway in mammalian cells, including some features observed in most cell types (lysosome biogenesis) and some observed only in specialized cells (phagocytosis, macropinocytosis, and lysosome secretion).Dictyostelium amoebae live in the soil, where they feed by ingesting and digesting other microorganisms. In addition, axenic laboratory strains can macropinocytose medium to ensure their growth. Accordingly, both in natural situations and in laboratory settings, the endocytic pathway plays a key role in the acquisition of nutrients by Dictyostelium cells. In agreement with this notion, several observations suggest that the physiology of the endocytic pathway is sensitive to nutrient availability. In particular, starvation induces secretion of lysosomal enzymes by an unknown mechanism (11). The morphology of the endocytic pathway is also sensitive to nutritional cues, as shown for example by the observation that formation of multilamellar endosomes is enhanced in cells fed with bacteria (18).Here, we analyzed the effect of starvation on the organization as well as the dynamics of the endocytic pathway. We found that, while the overall organization was not extensively modified in starved cells, the dynamics of endocytic compartments were altered. Moreover, analysis of two specific knockout mutants, the apm3 (6) and lvsB (8) strains, revealed that their phenotype was profoundly altered upon starvation, providing further insight about the role of Apm3 and LvsB in the endocytic pathway.     

Differential Behavior of Missense Mutations in the Intersubunit Contact Domain of the Human Pyruvate Kinase M2 Isozyme

In this study, we attempted to understand the mechanism of regulation of the activity and allosteric behavior of the pyruvate kinase M₂ enzyme and two of its missense mutations, H391Y and K422R, found in cells from Bloom syndrome patients, prone to develop cancer. Results show that despite the presence of mutations in the intersubunit contact domain, the K422R and H391Y mutant proteins maintained their homotetrameric structure, similar to the wild-type protein, but showed a loss of activity of 75 and 20%, respectively. Interestingly, H391Y showed a 6-fold increase in affinity for its substrate phosphoenolpyruvate and behaved like a non-allosteric protein with compromised cooperative binding. However, the affinity for phosphoenolpyruvate was lost significantly in K422R. Unlike K422R, H391Y showed enhanced thermal stability, stability over a range of pH values, a lesser effect of the allosteric inhibitor Phe, and resistance toward structural alteration upon binding of the activator (fructose 1,6-bisphosphate) and inhibitor (Phe). Both mutants showed a slight shift in the pH optimum from 7.4 to 7.0. Although this study signifies the importance of conserved amino acid residues in long-range communications between the subunits of multimeric proteins, the altered behavior of mutants is suggestive of their probable role in tumor-promoting growth and metabolism in Bloom syndrome patients with defective pyruvate kinase M₂.Pyruvate kinase (PK³; EC 2.7.1.40), a pacemaker of the glycolytic pathway, catalyzes irreversibly the transphosphorylation from P-enolpyruvate to ADP, generating pyruvate and ATP (1, 2). There are four different isozymes (L, R, M₁, and M₂) in mammalian tissues, which differ in their regulatory properties. These isozymes are allosteric in nature with the exception of the M₁ form, present in skeletal muscle and brain (3–7). PKM₂ is a ubiquitous prototype enzyme present in all tissues during the embryonic stage and is gradually replaced by other isozymic forms in specific tissues during development. The M₂, L, and R isozymes show homotropic cooperative activation with P-enolpyruvate and heterotropic cooperative activation with Fru-1,6-P₂ (8–10). The M₁ isozyme is regulated by neither P-enolpyruvate nor Fru-1,6-P₂ because of its intrinsic active conformation in the R-state (5, 6). Under unfavorable conditions such as hypoxia and lack of glucose supply, the anaerobic tissues and tumor cells rely heavily on PKM₂ for ATP production (7). Therefore, stringent control of PK activity is of great importance not only for cell metabolism but also for tumorigenic proliferation.The M₁ and M₂ isozymes are produced from a single gene locus by mutually exclusive alternative splicing (11–14). In the human M₁ and M₂ isozymes, the exon that is exchanged because of alternative splicing encodes 56 amino acids, in which a total of 22 amino acids differ within a length of 45 residues. The residues located in this region form the major intersubunit contact domain (8). The distinguishable kinetic properties of the M₁ and M₂ isozymes are attributed to these amino acid substitutions. It has been shown by x-ray crystallographic analyses and computer modeling that the corresponding regions of their polypeptides participate directly in the intersubunit contact, which is responsible for the intersubunit communication required for allosteric cooperativity (8, 15).PK has been largely conserved throughout evolution. The enzyme is usually a homotetramer composed of four identical subunits, and each subunit consists of four domains: the A-, B-, and C-domains and the N-terminal domain. The structure of human PKM₂ was recently determined in complex with inhibitors (16). In mammalian cells, PK activity is regulated by two different mechanisms: one at the level of expression and the other through allosteric regulation. The catalytic site usually composes a small part of the enzyme, but allosteric control is transmitted over a long range, thus increasing the number of possible residues involved in regulation. The allosteric transition in PK involves mutual rotations of the A- and C-domains within each subunit and the subunit within the tetramer (14). The residues at the subunit interfaces have the critical function of relaying the allosteric signal from and to the catalytic and regulatory sites. This region also transmits the allosteric signal between P-enolpyruvate- and Fru-1,6-P₂-binding sites. Despite the availability of structural details of several PK isozymes, it is difficult to identify the structural elements that play an important role in PK regulation and propagation of the allosteric signals. Although the role of some of the PK residues (positions 340, 389, 398, 401, 402, 408, 423, and 427) has been studied in allosteric regulation (10, 17–19) by in vitro site-directed mutagenesis, the absence of these mutations in any naturally occurring condition presents limitations in attributing a biological role to the introduced changes.The natural mutations H391Y and K422R (reported previously as K421R) were reported by us for the first time in the PKM₂ gene in a Bloom syndrome cell line and in the lymphocytes of an Indian Bloom syndrome patient, respectively (20). The two missense mutations, located in the region of the intersubunit contact domain (Fig. 1, A and B), presented with the biochemical phenotype of down-regulated enzyme activity to different extents (20) and were expected to influence the allosteric nature of the enzyme. The regulatory behavior of allosteric PK has been described by a two-state model that proposes an active (R) and an inactive (T) form of the macromolecule with differential affinity for ligands (15). Upon binding of the substrate or its analogs, the enzyme undergoes a transition from a low activity/low affinity conformation (T state) to a high activity/high affinity conformation (R state). The binding of phenylalanine produces a global structural change and exhibits reduced affinity for substrate P-enolpyruvate in the T state (21–23). Previous studies have demonstrated that each individual domain acts as a rigid body and that, upon transition from the T to the R state, the domain of the functional tetramer modifies its relative orientation by 29°. These movements bring conformational change to the active site, which, upon transition to the T state, undergoes a distortion of the P-enolpyruvate-binding site (24).Open in a separate windowFIGURE 1.A, ribbon diagram of the overall structure of PK showing the positions of the two mutations, H391Y and K422R, along with the active site and Fru-1,6-P₂-binding site. B, intersubunit contact domain of PK. The major amino acid residues and side chains at the tetramer interface region are shown.Because the mutations observed by us previously (20) are located at highly conserved positions not only in different isozymic forms but also across the species (supplemental Fig. S1) and are observed in the genetic background of a syndrome prone to cancer in early age, a study related to the structure-function correlations of these mutations is likely to provide insight into their possible biological importance, especially in the context of recent research highlighting the importance of PKM₂ in tumor promotion and growth. In this study, we investigated the role of the two natural missense mutations, after site-directed mutagenesis in the PKM₂ gene, in the regulation of allosteric properties as well as their effects on the secondary and tertiary structures in comparison with wild-type PKM₂ (PK-WT). An attempt has also been made to understand the effects of these mutations at the interface of the subunits on the signal transmission pathway within the protein.     

In Silico Determination of the Effect of Multi-Target Drugs on Calcium Dynamics Signaling Network Underlying Sea Urchin Spermatozoa Motility

The motility of spermatozoa of both Lytechinus pictus and Strongylocentrotus purpuratus sea urchin species is modulated by the egg-derived decapeptide speract via an oscillatory [Ca²⁺]-dependent signaling pathway. Comprehension of this pathway is hence directly related to the understanding of regulated sperm swimming. Niflumic acid (NFA), a nonsteroidal anti-inflammatory drug alters several ion channels. Though unspecific, NFA profoundly affects how sea urchin sperm respond to speract, increasing the [Ca²⁺]_i oscillation period, amplitude, peak and average level values of the responses in immobilized and swimming cells. A previous logical network model we developed for the [Ca²⁺] dynamics of speract signaling cascade in sea urchin sperm allows integrated dissection of individual and multiple actions of NFA. Among the channels affected by NFA are: hyperpolarization-activated and cyclic nucleotide gated Na⁺ channels (HCN), [Ca²⁺]-dependent Cl⁻ channels (CaCC) and [Ca²⁺]-dependent K⁺ channels (CaKC), all present in the sea urchin genome. Here, using our model we investigated the effect of blocking in silico HCN and CaCC channels suggested by experiments. Regarding CaKC channels, arguments can be provided for either their blockage or activation by NFA. Our study yielded two scenarios compliant with experimental observations: i) under CaKC inhibition, this [Ca²⁺]-dependent K⁺ channel should be different from the Slo1 channel and ii) under activation of the CaKC channel, another [Ca²⁺] channel not considered previously in the network is required, such as the pH-dependent CatSper channel. Additionally, our findings predict cause-effect relations resulting from a selective inhibition of those channels. Knowledge of these relations may be of consequence for a variety of electrophysiological studies and have an impact on drug related investigations. Our study contributes to a better grasp of the network dynamics and suggests further experimental work.     

In Silico Target-Specific siRNA Design Based on Domain Transfer in Heterogeneous Data

RNA interference via exogenous small interference RNAs (siRNA) is a powerful tool in gene function study and disease treatment. Designing efficient and specific siRNA on target gene remains the key issue in RNAi. Although various in silico models have been proposed for rational siRNA design, most of them focus on the efficiencies of selected siRNAs, while limited effort has been made to improve their specificities targeted on specific mRNAs, which is related to reducing off-target effects (OTEs) in RNAi. In our study, we propose for the first time that the enhancement of target specificity of siRNA design can be achieved computationally by domain transfer in heterogeneous data sources from different siRNA targets. A transfer learning based method i.e., heterogeneous regression (HEGS) is presented for target-specific siRNA efficacy modeling and feature selection. Based on the model, (1) the target regression model can be built by extracting information from related data in other targets/experiments, thus increasing the target specificity in siRNA design with the help of information from siRNAs binding to other homologous genes, and (2) the potential features correlated to the current siRNA design can be identified even when there is lack of experimental validated siRNA affinity data on this target. In summary, our findings present useful instructions for a better target-specific siRNA design, with potential applications in genome-wide high-throughput screening of effective siRNA, and will provide further insights on the mechanism of RNAi.     

In Silico Screening of the Key Cellular Remodeling Targets in Chronic Atrial Fibrillation

Chronic atrial fibrillation (AF) is a complex disease with underlying changes in electrophysiology, calcium signaling and the structure of atrial myocytes. How these individual remodeling targets and their emergent interactions contribute to cell physiology in chronic AF is not well understood. To approach this problem, we performed in silico experiments in a computational model of the human atrial myocyte. The remodeled function of cellular components was based on a broad literature review of in vitro findings in chronic AF, and these were integrated into the model to define a cohort of virtual cells. Simulation results indicate that while the altered function of calcium and potassium ion channels alone causes a pronounced decrease in action potential duration, remodeling of intracellular calcium handling also has a substantial impact on the chronic AF phenotype. We additionally found that the reduction in amplitude of the calcium transient in chronic AF as compared to normal sinus rhythm is primarily due to the remodeling of calcium channel function, calcium handling and cellular geometry. Finally, we found that decreased electrical resistance of the membrane together with remodeled calcium handling synergistically decreased cellular excitability and the subsequent inducibility of repolarization abnormalities in the human atrial myocyte in chronic AF. We conclude that the presented results highlight the complexity of both intrinsic cellular interactions and emergent properties of human atrial myocytes in chronic AF. Therefore, reversing remodeling for a single remodeled component does little to restore the normal sinus rhythm phenotype. These findings may have important implications for developing novel therapeutic approaches for chronic AF.     

In Silico Adoption of an Orphan Nuclear Receptor NR4A1

A 4.1μs molecular dynamics simulation of the NR4A1 (hNur77) apo-protein has been undertaken and a previously undetected druggable pocket has become apparent that is located remotely from the ‘traditional’ nuclear receptor ligand-binding site. A NR4A1/bis-indole ligand complex at this novel site has been found to be stable over 1 μs of simulation and to result in an interesting conformational transmission to a remote loop that has the capacity to communicate with a NBRE within a RXR-α/NR4A1 heterodimer. Several features of the simulations undertaken indicate how NR4A1 can be affected by alternate-site modulators.     

In Silico Prediction of Mutant HIV-1 Proteases Cleaving a Target Sequence

HIV-1 protease represents an appealing system for directed enzyme re-design, since it has various different endogenous targets, a relatively simple structure and it is well studied. Recently Chaudhury and Gray (Structure (2009) 17: 1636–1648) published a computational algorithm to discern the specificity determining residues of HIV-1 protease. In this paper we present two computational tools aimed at re-designing HIV-1 protease, derived from the algorithm of Chaudhuri and Gray. First, we present an energy-only based methodology to discriminate cleavable and non cleavable peptides for HIV-1 proteases, both wild type and mutant. Secondly, we show an algorithm we developed to predict mutant HIV-1 proteases capable of cleaving a new target substrate peptide, different from the natural targets of HIV-1 protease. The obtained in silico mutant enzymes were analyzed in terms of cleavability and specificity towards the target peptide using the energy-only methodology. We found two mutant proteases as best candidates for specificity and cleavability towards the target sequence.