Search PubMed⌕ Search

Biomedical subjects

Enrico O Purisima

Publications and source records attributed to Enrico O Purisima.

16 recordsLinked to original sources
PubMed2006

MicroRNAs preferentially target the genes with high transcriptional regulation complexity.

Over the past few years, microRNAs (miRNAs) have emerged as a new prominent class of gene regulatory factors that negatively regulate expression of approximately one-third of the genes in animal genomes at post-transcriptional level. However, it is still unclear why some genes are regulated by miRNAs but others are not, i.e. what principles govern miRNA regulation in animal genomes. In this study, we systematically analyzed the relationship between transcription factors (TFs) and miRNAs in gene regulation. We found that the genes with more TF-binding sites have a higher probability of being targeted by miRNAs and have more miRNA-binding sites on average. This observation reveals that the genes with higher cis-regulation complexity are more coordinately regulated by TFs at the transcriptional level and by miRNAs at the post-transcriptional level. This is a potentially novel discovery of mechanism for coordinated regulation of gene expression. Gene ontology analysis further demonstrated that such coordinated regulation is more popular in the developmental genes.

Chromosome Mapping↗
PubMed2006

Principles of microRNA regulation of a human cellular signaling network.

MicroRNAs (miRNAs) are endogenous approximately 22-nucleotide RNAs, which suppress gene expression by selectively binding to the 3'-noncoding region of specific messenger RNAs through base-pairing. Given the diversity and abundance of miRNA targets, miRNAs appear to functionally interact with various components of many cellular networks. By analyzing the interactions between miRNAs and a human cellular signaling network, we found that miRNAs predominantly target positive regulatory motifs, highly connected scaffolds and most downstream network components such as signaling transcription factors, but less frequently target negative regulatory motifs, common components of basic cellular machines and most upstream network components such as ligands. In addition, when an adaptor has potential to recruit more downstream components, these components are more frequently targeted by miRNAs. This work uncovers the principles of miRNA regulation of signal transduction networks and implies a potential function of miRNAs for facilitating robust transitions of cellular response to extracellular signals and maintaining cellular homeostasis.

Computational Biology↗
PubMed2006

Binding site-based classification of coronaviral papain-like proteases.

The coronavirus replicase gene encodes one or two papain-like proteases (termed PL1pro and PL2pro) implicated in the N-terminal processing of the replicase polyprotein and thus contributing to the formation of the viral replicase complex that mediates genome replication. Using consensus fold recognition with the 3D-JURY meta-predictor followed by model building and refinement, we developed a structural model for the single PLpro present in the severe acute respiratory syndrome coronavirus (SCoV) genome, based on significant structural relationships to the catalytic core domain of HAUSP, a ubiquitin-specific protease (USP). By combining the SCoV PLpro model with comparative sequence analyses we show that all currently known coronaviral PLpros can be classified into two groups according to their binding site architectures. One group includes all PL2pros and some of the PL1pros, which are characterized by a restricted USP-like binding site. This group is designated the R-group. The remaining PL1pros from some of the coronaviruses form the other group, featuring a more open papain-like binding site, and is referred to as the O-group. This two-group, binding site-based classification is consistent with experimental data accumulated to date for the specificity of PLpro-mediated polyprotein processing and PLpro inhibition. It also provides an independent evaluation of the similarity-based annotation of PLpro-mediated cleavage sites, as well as a basis for comparison with previous groupings based on phylogenetic analyses.

Amino Acid Sequence↗
PubMed2006

Impaired helix 12 dynamics due to proline 892 substitutions in the androgen receptor are associated with complete androgen insensitivity.

Structural studies of the ligand-binding domain (LBD) of several steroid receptors have revealed that the dynamic properties of the C-terminal helix 12 (H12) are the major determinant of the activation mode of these receptors. H12 exhibits high mobility and different conformations in the absence of ligand. Upon ligand binding, H12 is stabilized in a precise position to seal the ligand-binding pocket and finalize the assembly of the activation function (AF-2) domain. In this study, we investigated the role of the conserved proline 892 of the androgen receptor (AR) in directing the dynamic location and orientation of the AR-H12. We used a combined approach including kinetic and biochemical assays with molecular dynamic simulations to analyze two substitutions (P892A and P892L) identified in individuals with complete androgen insensitivity syndrome. Our analyses revealed distinct mechanisms by which these substitutions impair H12 function resulting in severely defective receptors. The AR-P892A receptor exhibited reduced ligand binding and transactivational potential because of an increased flexibility in H12. The AR-P892L substitution renders the receptor inactive due to a distorted, unstructured and misplaced H12. To confirm the mutants' inability to stabilize H12 in an active position, we have developed a novel in vivo assay to evaluate the accessibility of the H12-docking site on the AR-LBD surface. An extrinsic AR-H12 peptide was able to interact with wild-type and mutant LBDs in the absence of ligand. Ligand-induced proper positioning of the intrinsic H12 of wild-type AR prevented these interactions, whereas the misplacement of the mutants' H12 did not. Proline at this position may be critical for H12 dynamics not only in the AR, but also in other nuclear receptors where this proline is conserved.

Amino Acid Sequence↗
PubMed2006

Molecular surface generation using a variable-radius solvent probe.

Protein-ligand binding occurs through interactions at the molecular surface. Hence, a proper description of this surface is essential to our understanding of the process of molecular recognition. Recent studies have noted the inadequacy of using a fixed 1.4 A solvent probe radius to generate the molecular surface. This assumes that water molecules approach all surface atoms at an equal distance irrespective of polarity, which is not the case. To adequately model the protein-water boundary requires that the solvent probe radius change according to the polarity of its contacting atoms, smaller near polar atoms and larger near apolar atoms. To our knowledge, no method currently exists to generate the molecular surface of a protein in this manner. Using a modification of the marching tetrahedra algorithm, we present a method to generate molecular surfaces using a variable radius solvent probe. The resulting surface lacks many of the unrealistic small crevices in nonpolar regions that are found when utilizing an invariant 1.4 A solvent probe, while maintaining the fine detail of the surface at polar regions. On application of the method on a test set of 20 protein structures taken from the Protein Data Bank (PDB), we also find far fewer empty unsolvated cavities that are present when using only a 1.4 A solvent probe, while the majority of solvated and polar cavities is retained. This suggests that the majority of empty cavities previously observed in protein structures might simply be artifacts of the surfacing method. We also find that the variable probe surface can have significant effects on electrostatic calculations by generating a better tuned description of the protein-water boundary. We also examined the binding interfaces of a diverse set of 55 protein-protein complexes. We find that using a variable probe results in an increase in perceived shape complementarity at these sites compared to using a 1.4 A solvent probe. The molecular volume and surface area are geometric values that determine various important properties for macromolecules, and the altered description afforded by a variable solvent probe molecular surface can have significant implications in protein recognition, energetics, folding, and stability calculations.

Entropy↗
PubMed2006

Coupled atomic charge selectivity for optimal ligand-charge distributions at protein binding sites.

Charge optimization as a tool for both analyzing and enhancing binding electrostatics has become an attractive approach over the past few years. An interesting feature of this method for molecular design is that it provides not only the optimal charge magnitudes, but also the selectivity of a particular atomic center for its optimal charge. The current approach to compute the charge selectivity at a given atomic center of a ligand in a particular binding process is based on the binding-energy cost incurred upon the perturbation of the optimal charge distribution by a unit charge at the given atomic center, while keeping the other atomic partial charges at their optimal values. A limitation of this method is that it does not take into account the possible concerted changes in the other atomic charges that may incur a lower energetic cost than perturbing a single charge. Here, we describe a novel approach for characterizing charge selectivity in a concerted manner, taking into account the coupling between the ligand charge centers in the binding process. We apply this novel charge selectivity measure to the celecoxib molecule, a nonsteroidal anti-inflammatory agent binding to cyclooxygenase 2 (COX2), which has been recently shown to also exhibit cross-reactivity toward carbonic anhydrase II (CAII), to which it binds with nanomolar affinity. The uncoupled and coupled charge selectivity profiles over the atomic centers of the celecoxib ligand, binding independently to COX2 and CAII, are analyzed comparatively and rationalized with respect to available experimental data. Very different charge selectivity profiles are obtained for the uncoupled versus coupled selectivity calculations.

Binding Sites↗
PubMed2005

Inhibition of a cathepsin L-like cysteine protease by a chimeric propeptide-derived inhibitor.

Like other papain-related cathepsins, congopain from Trypanosoma congolense is synthesized as a zymogen. We have previously identified a proregion-derived peptide (Pcp27), acting as a weak and reversible inhibitor of congopain. Pcp27 contains a 5-mer YHNGA motif, which is essential for selectivity in the inhibition of its mature form [Lalmanach, G., Lecaille, F., Chagas, J. R., Authié, E., Scharfstein, J., Juliano, M. A., and Gauthier, F. (1998) J. Biol. Chem. 273, 25112-25116]. In the work presented here, a homology model of procongopain was generated and subsequently used to model a chimeric 50-mer peptide (called H3-Pcp27) corresponding to the covalent linkage of an unrelated peptide (H3 helix from Antennapedia) to Pcp27. Molecular simulations suggested that H3-Pcp27 (pI = 9.99) maintains an N-terminal helical conformation, and establishes more complementary electrostatic interactions (E(coul) = -25.77 kcal/mol) than 16N-Pcp27, the 34-mer Pcp27 sequence plus the 16 native residues upstream from the proregion (E(coul) = 0.20 kcal/mol), with the acid catalytic domain (pI = 5.2) of the mature enzyme. In silico results correlated with the significant improvement of congopain inhibition by H3-Pcp27 (K(i) = 24 nM), compared to 16N-Pcp27 (K(i) = 1 microM). In addition, virtual alanine scanning of H3 and 16N identified the residues contributing most to binding affinity. Both peptides did not inhibit human cathepsins B and L. In conclusion, these data support the notion that the positively charged H3 helix favors binding, without modifying the selectivity of Pcp27 for congopain.

Amino Acid Sequence↗
PubMed2005

Carboxy-monopeptidase substrate specificity of human cathepsin X.

Cathepsin X is a papain-like cysteine protease with restricted positional specificity, acting primarily as a carboxy-monopeptidase. We mapped the specificities at the S2, S1, and S1' subsites of human cathepsin X by systematically and independently substituting the P2, P1, and P1' positions of the carboxy-monopeptidase substrate Abz-FRF(4NO(2)) with natural amino acids. Human cathepsin X has broad S2, S1, and S1' specificities within two orders of magnitude in k(cat)/K(M), excluding proline that is not tolerated at these subsites. Glycine is not favored in S2, but is among the preferred residues in S1 and S1', which highlights S2 as the affinity-determinant subsite. The presence of peculiar residues at several binding site positions (Asp76, His234, Asn75, and Glu72) does not translate into a markedly different sequence specificity profile relative to other human cathepsins. These findings suggest that a specific function of human cathepsin X is unlikely to result from sequence specificity, but rather from a combination of its unique positional specificity and the co-localization of enzyme and substrate in a specific cellular environment.

Amino Acid Sequence↗
PubMed2003

Bridging structural biology and genetics by computational methods: an investigation into how the R774C mutation in the AR gene can result in complete androgen insensitivity syndrome.

Recent structural studies of the ligand-binding domain (LBD) of the androgen receptor (AR) have raised more questions than answers, as most of the known pathogenic mutations of the AR gene causing androgen insensitivity syndrome (AIS) are not in the ligand-binding pocket. In this study, we have investigated one such pathogenic mutation, by examining details of its altered atomic structure using a computational technique of molecular dynamics (MD) simulations extended over 4 ns, effectively creating a 4D structural model. The mutation R774C, which is in the LBD of the AR gene, causes complete AIS (CAIS), producing ARs that have a unique thermolabile profile, being thermostable at 22 degrees C but thermolabile at 37 degrees C. We have therefore investigated this mutation by MD simulations at 293 K (20 degrees C), 300 K (27 degrees C), and 310 K (37 degrees C). The MD simulations indicate that: 1) the mutation causes local structural distortions, which result in changes in the shape of the ligand-binding pocket; 2) the mutation alters the dynamic nature of the protein and results in a more diverse conformational distribution of the ligand-binding pocket; and 3) the effect of the mutation on AR structure could be largely reversed by lowering the temperature at which the MD simulations were conducted. These results therefore strongly support the biochemical data, e.g., the mutants' inability to form AR-ligand complexes at 37 degrees C and its characteristic reversible thermolability, clearly indicating the value of such computational methods.

Androgen-Insensitivity Syndrome↗
PubMed2003

On the transferability of hydration-parametrized continuum electrostatics models to solvated binding calculations.

Using molecular mechanics force field partial atomic charges, we show the nonuniqueness of the parametrization of continuum electrostatics models with respect to solute atomic radii and interior dielectric constant based on hydration (vacuum-to-water transfer) free energy data available for small molecules. Moreover, parameter sets that are optimal and equivalent for hydration free energy calculations lead to large variations of calculated absolute and relative electrostatic binding free energies. Hence, parametrization of solvation effects based on hydration data, although a necessary condition, is not sufficient to guarantee its transferability to the calculation of binding free energies in solution.

Journal Article↗
PubMed2003

Profiling charge complementarity and selectivity for binding at the protein surface.

A novel analysis and representation of the protein surface in terms of electrostatic binding complementarity and selectivity is presented. The charge optimization methodology is applied in a probe-based approach that simulates the binding process to the target protein. The molecular surface is color coded according to calculated optimal charge or according to charge selectivity, i.e., the binding cost of deviating from the optimal charge. The optimal charge profile depends on both the protein shape and charge distribution whereas the charge selectivity profile depends only on protein shape. High selectivity is concentrated in well-shaped concave pockets, whereas solvent-exposed convex regions are not charge selective. This suggests the synergy of charge and shape selectivity hot spots toward molecular selection and recognition, as well as the asymmetry of charge selectivity at the binding interface of biomolecular systems. The charge complementarity and selectivity profiles map relevant electrostatic properties in a readily interpretable way and encode information that is quite different from that visualized in the standard electrostatic potential map of unbound proteins.

Binding Sites↗
PubMed2002

Design of noncovalent inhibitors of human cathepsin L. From the 96-residue proregion to optimized tripeptides.

A novel series of noncovalent inhibitors of cathepsin L have been designed to mimic the mode of autoinhibition of procathepsin L. Just like the propeptide, these peptide-based inhibitors have a reverse-binding mode relative to a substrate and span both the S' and S subsites of the enzyme active site. In contrast to previous studies in which even moderate truncation of the full-length propeptide led to rapid reduction in potency, these blocked tripeptide-sized inhibitors maintain nanomolar potency. Moreover, these short peptides show higher selectivity (up to 310-fold) for inhibiting cathepsin L over K versus only 2-fold selectivity of the 96-residue propeptide of cathepsin L. A 1.9 A X-ray crystallographic structure of the complex of cathepsin L with one of the inhibitors confirms the designed reverse-binding mode of the inhibitor as well as its noncovalent nature. Enzymatic analysis also shows the inhibitors to be resistant to hydrolysis at elevated concentrations of the enzyme. The mode of inhibition of these molecules provides a general strategy for inhibiting other cathepsins as well as other proteases.

Binding Sites↗
PubMed2002

Implications of the reactive thiol and the proximal non-proline cis-peptide bond in the Structure and function of Vibrio harveyi luciferase.

The role of a highly reactive cysteine residue, Cys106, in Vibrio harveyi luciferase in modulating the substrate-enzyme interactions and in turn affecting the enzyme activity has been extensively investigated over the past three decades. Replacing Cys106 with valine dramatically hinders the ability of luciferase to stabilize the C4a-hydroperoxyflavin intermediate [Abu-Soud, H. M., Clark, A. C., Francisco, W. A., Baldwin, T. O., and Raushel, F. M. (1993) J. Biol. Chem. 268, 7699-7706] and consume aldehyde substrate [Xi, L., Cho, K.-W., Herndon, M. E., and Tu, S.-C. (1990) J. Biol. Chem. 265, 4200-4203], therefore markedly decreasing enzyme activity. On the basis of the structure-activity relationship of flavin analogues and the location of the phosphate binding site of flavin mononucleotide (FMN) coupled with molecular modeling, the functional part of the isoalloxazine ring of FMN, the thiol side chain of Cys106, the methyl group of Ala75, and the unique non-prolyl cis-peptide bond between Ala74 and Ala75 were found to be closely packed [Lin, L. Y., Sulea, T., Szittner, R., Vassilyev, V., Purisima, E. O., and Meighen, E. A. (2001) Protein Sci. 10, 1563-1571]. Here, by mutating Ala75 to Gly, we restored key wild-type properties to the C106V mutant, in particular, high enzyme activity and a stable C4a-hydroperoxyflavin intermediate, demonstrating that the primary reason for the dark phenotype of the C106V mutant was the unfavorable steric interaction between Val106 and Ala75 side chains, which could in turn disturb the cis-oriented amide linkage of Ala74 and Ala75. Moreover, significant red shifts in light emission of 3-10 nm were measured for luciferases carrying Val106 with the spectrum of the double mutant C106V/A75G now red shifted to that of Photobacterium phosphoreum luciferase, which also has Val and Gly at positions 106 and 75, respectively. These results strengthen the validity of the binding geometry of the modeled flavin with the re-face of the pyrimidine end of the isoalloxazine ring next to Cys106 and implicate the Ala74-Ala75 cis-peptide as a key component in the bioluminescence reaction.

Enzyme Stability↗
PubMed2002

A fast pairwise evaluation of molecular surface area.

A fast and general analytical approach was developed for the calculation of the approximate van der Waals and solvent-accessible surface areas. The method is based on three basic ideas: the use of the Lorentz transformation formula, a rigid-geometry approximation, and a single fitting parameter that can be refitted on the fly during a simulation. The Lorentz transformation equation is used for the summation of the areas of an atom buried by its neighboring contacting atoms, and implies that a sum of the buried pairwise areas cannot be larger than the surface area of the isolated spherical atom itself. In a rigid-geometry approximation we numerically calculate and keep constant the surface of each atom buried by the atoms involved in 1-2 and 1-3 interactions. Only the contributions from the nonbonded atoms (1-4 and higher interactions) are considered in terms of the pairwise approximation. The accuracy and speed of the method is competitive with other pairwise algorithms. A major strength of the method is the ease of parametrization.

Journal Article↗
PubMedDate not supplied

Effect of limited solid-state glycation on the conformation of lysozyme by ESI-MSMS peptide mapping and molecular modeling.

Although protein glycation has been implicated in the alteration of protein functionality, both in vivo (in biological systems) and in vitro (in food systems), the effect of the protein-bound glycan moiety on the structure/conformation of proteins that result in the modification of functionality is not clear. In this article, we report a study of the conformational changes of glycated lysozyme using LC-ESI-MSMS peptide mapping, and molecular modeling. A comparison of the RP-HPLC of the tryptic digests of unglycated and glycated lysozyme showed markedly different chromatographic profiles. Analysis of the peptide composition of the chromatographic fractions of the tryptic digests revealed that glycation of lysozyme resulted in the modification of its conformation. Glycation-induced changes in the conformation of lysozyme resulted in the exposure of its active site region to increased proteolytic activity of trypsin. Molecular simulation of triglycated lysozyme also showed that limited glycation of lysozyme caused reorientation of the active site residues (Arg 45, Arg 68, Asn 44, and Trp 62) and increased solvent accessibility into the active site region of the protein. The results of the modeling experiment corroborated the results of the RP-HPLC and ESI-MSMS peptide mapping.

Animals↗