Using mechanism similarity to understand enzyme evolution

Enzyme reactions take place in the active site through a series of catalytic steps, which are collectively termed the enzyme mechanism. The catalytic step is thereby the individual unit to consider for the purposes of building new enzyme mechanisms — i.e. through the mix and match of individual catalytic steps, new enzyme mechanisms and reactions can be conceived. In the case of natural evolution, it has been shown that new enzyme functions have emerged through the tweaking of existing mechanisms by the addition, removal, or modification of some catalytic steps, while maintaining other steps of the mechanism intact. Recently, we have extracted and codified the information on the catalytic steps of hundreds of enzymes in a machine-readable way, with the aim of automating this kind of evolutionary analysis. In this paper, we illustrate how these data, which we called the “rules of enzyme catalysis”, can be used to identify similar catalytic steps across enzymes that differ in their overall function and/or structural folds. A discussion on a set of three enzymes that share part of their mechanism is used as an exemplar to illustrate how this approach can reveal divergent and convergent evolution of enzymes at the mechanistic level.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12551-022-01022-9.     

Using microarrays and functional genomics to understand and target cancer

Conclusion In our work microarrays are of importance for screening the whole genome to sort out important genes and also for constructing pathways systems for the understanding of the links between genes. Our studies demonstrate the potential clinical usefulness of microarrays for the classification of patients and for understanding cancer biology and identifying new drug targets.     

Novel aspects of calmodulin target recognition and activation.

Several crystal and NMR structures of calmodulin (CaM) in complex with fragments derived from CaM-regulated proteins have been reported recently and reveal novel ways for CaM to interact with its targets. This review will discuss and compare features of the interaction between CaM and its target domains derived from the plasma membrane Ca2+-pump, the Ca2+-activated K+-channel, the Ca2+/CaM-dependent kinase kinase and the anthrax exotoxin. Unexpected aspects of CaM/target interaction observed in these complexes include: (a) binding of the Ca2+-pump domain to only the C-terminal part of CaM (b) dimer formation with fragments of the K+-channel (c) insertion of CaM between two domains of the anthrax exotoxin (d) binding of Ca2+ ions to only one EF-hand pair and (e) binding of CaM in an extended conformation to some of its targets. The mode of interaction between CaM and these targets differs from binding conformations previously observed between CaM and peptides derived from myosin light chain kinase (MLCK) and CaM-dependent kinase IIalpha (CaMKIIalpha). In the latter complexes, CaM engulfs the CaM-binding domain peptide with its two Ca2+-binding lobes and forms a compact, ellipsoid-like complex. In the early 1990s, a model for the activation of CaM-regulated proteins was developed based on this observation and postulated activation through the displacement of an autoinhibitory or regulatory domain from the target protein upon binding of CaM. The novel structures of CaM-target complexes discussed here demonstrate that this mechanism of activation may be less general than previously believed and seems to be not valid for the anthrax exotoxin, the CaM-regulated K+-channel and possibly also not for the Ca2+-pump.     

A metadynamic approach to understand the recognition mechanism of the histone H3 tail with the ATRXADD domain

The binding affinity between the histone 3 (H3) tail and the ADD domain of ATRX (ATRX_ADD) increases with the subsequent addition of methyl groups on lysine 9 on H3. To improve our understanding of how the difference in methylation state affects binding between H3 and the ATRX_ADD, we adopted a metadynamic approach to explore the recognition mechanism between the two proteins and identify the key intermolecular interactions that mediate this protein-peptide interaction (PPI). The non-methylated H3 peptide is recognized only by the PHD finger of ATRX_ADD while mono-, di-, and trimethylated H3 is recognized by both the PHD and GATA-like zinc finger of the domain. Furthermore, water molecules play an important role in orienting the lysine 9 anchor towards the GATA-like zinc finger, which results in stabilizing the lysine 9 binding pocket on ATRX_ADD. We compared our computational results against experimentally determined NMR and X-ray structures by demonstrating the RMSD, order parameter S² and hydration number of the complex. The metadynamics data provide new insight into roles of water-bridges and the mechanisms through which K9 hydration stabilizes the H3K9me3:ATRX_ADD PPI, providing context for the high affinity demonstrated between this protein and peptide.     

Modelling of lysozyme binding to a cation exchange surface at atomic detail: the role of flexibility

Different approaches were made to predict the adsorbed orientation based on rigid, flexible, or a mixture of both models. To determine the role of flexibility during adsorption, the orientation of lysozyme adsorbed to a negatively charged ligand surface was predicted by a rigid and a flexible model based on two differing protein structures at atomic resolution. For the rigid model, the protein structures were placed at different distances from the ligand surface and the electrostatic interaction energy was calculated for all possible orientations. The results were compared to a flexible model where the binding to the ligand surface was modeled by multiple molecular dynamics simulations starting with 14 initial orientations. Different aspects of the adsorption process were not covered by the rigid model and only detectable by the flexible model. Whereas the results of the rigid model depended sensitively on the protein-surface distance and the protein structure, the preferred orientation obtained by the flexible model was closer to a previous experimental determined orientation, robust toward the initial orientation and independent of the initial protein structure. Additionally, it was possible to obtain insights into the preferred binding process of lysozyme on a negatively charged surface by the flexible model.     

Using the recognition code to swap homeodomain target specificity in cell culture

The homeodomain (HD) is a 60 amino acid-long DNA-binding domain. A large fraction of HDs binds with high affinity sequences containing the 5′-TAAT-3′ core motif. However, NK-2 class HDs recognizes sequences containing the 5′-CAAG-3′ core motif. By using a cell transfection approach, here we show that modification of residues located in the N-terminal arm (at positions 6, 7 and 8) and in the recognition helix (at position 54) is enough to swap the “in vivo” binding specificity of TTF-1 HD (which is a member of the NK-2 class HD) from 5′-CAAG-3′ to 5′-TAAT-3′-containing targets. The role of residue at position 54 is also supported by data obtained with the HD of the Drosophila engrailed protein. These data support the notion that DNA-binding specificity “in vivo” is dictated by few critical residues.     

Using mathematical models to understand the AIDS epidemic

The most urgent public-health problem today is to devise effective strategies to minimize the destruction caused by the AIDS epidemic. This complex problem will involve medical advances and new public-health and education initiatives. Mathematical models based on the underlying transmission mechanisms of the AIDS virus can help the medical/scientific community understand and anticipate its spread in different populations and evaluate the potential effectiveness of different approaches for bringing the epidemic under control. Before we can use models to predict the future, we must carefully test them against the past spread of the infection and for sensitivity to parameter changes. The long and extremely variable incubation period and the low probability of transmitting the AIDS virus in a single contact imply that population structure and variations in infectivity both play an important role in its spread. The population structure is caused by differences between people in numbers of sexual partners and the use of intravenous drugs and because of the way in which people mix among age, ethnic, and social groups. We use a simplified approach to investigate the effects of variation in incubation periods and infectivity specific to the AIDS virus, and we compare a model of random partner choices with a model in which partners both come from similar behavior groups.     

Using evolutionary information for the query and target improves fold recognition

In this study, we show that it is possible to increase the performance over PSI-BLAST by using evolutionary information for both query and target sequences. This information can be used in three different ways: by sequence linking, profile-profile alignments, and by combining sequence-profile and profile-sequence searches. If only PSI-BLAST is used, 16% of superfamily-related protein domains can be detected at 90% specificity, but if a sequence-profile and a profile-sequence search are combined, this is increased to 20%, profile-profile searches detects 19%, whereas a linking procedure identifies 22% of these proteins. All three methods show equal performance, but the best combination of speed and accuracy seems to be obtained by the combined searches, because this method shows a good performance even at high specificity and the lowest computational cost. In addition, we show that the E-values reported by all these methods, including PSI-BLAST, underestimate the true rate of false positives. This behavior is seen even if a very strict E-value cutoff and a limited number of iterations are used. However, the difference is more pronounced with a looser E-value cutoff and more iterations.     

Interaction of calmodulin with the phosphofructokinase target sequence

Ca4.calmodulin (Ca4.CaM) inhibits the glycolytic enzyme phosphofructokinase, by preventing formation of its active tetramer. Fluorescence titrations show that the affinity of complex formation of Ca4.CaM with the key 21-residue target peptide increases 1000-fold from pH 9.0 to 4.8, suggesting the involvement of histidine and carboxylic acid residues. 1H NMR pH titration indicates a marked increase in pKa of the peptide histidine on complex formation and HSQC spectra show related pH-dependent changes in the conformation of the complex. This unusually strong sensitivity of a CaM-target complex to pH suggests a potential functional role for Ca4.CaM in regulation of the glycolytic pathway.     

Using genetics to understand the dynamics of wild primate populations

While much can be learned about primates by means of observation, the slow life history of many primates means that even decades of dedicated effort cannot illuminate long-term evolutionary processes. For example, while the size of a contemporary population can be estimated from field censuses, it is often desirable to know whether a population has been constant or changing in size over a time frame of hundreds or thousands of years. Even the nature of “a population” is open to question, and the extent to which individuals successfully disperse among defined populations is also difficult to estimate by using observational methods alone. Researchers have thus turned to genetic methods to examine the size, structure, and evolutionary histories of primate populations. Many results have been gained by study of sequence variation of maternally inherited mitochondrial DNA, but in recent years researchers have been increasingly focusing on analysis of short, highly variable microsatellite segments in the autosomal genome for a high-resolution view of evolutionary processes involving both sexes. In this review we describe some of the insights thus gained, and discuss the likely impact on this field of new technologies such as high-throughput DNA sequencing.     

Apocalmodulin binds to the myosin light chain kinase calmodulin target site.

The interaction of a 20-residue-long peptide derived from the calmodulin-binding domain of the smooth muscle myosin light chain kinase with calcium-free calmodulin (apocalmodulin) was studied using a combination of isothermal titration calorimetry and differential scanning calorimetry. We showed that: (i) a significant binding between apocalmodulin and the target peptide (RS20) exists in the absence of salt (Ka = 10(6) M-1), (ii) the peptide interacts with the C-terminal lobe of calmodulin and adopts a partly helical conformation, and (iii) the presence of salt weakens the affinity of the peptide for apocalmodulin, emphasizing the importance of electrostatic interactions in the complex. Based on these results and taking into account the work of Bayley et al. (Bayley, P. M., Findlay, W.A., and Martin, S. R. (1996) Protein Sci. 5, 1215-1228), we suggest a physiological role for apocalmodulin.     

Genetics and bone. Using the mouse to understand man

The rationale for use of inbred strains of mice in bone research is well recognized and includes: a) practical factors (economics of scale, rapid development of adult status, pre-existing knowledge, down-sized technologies) and b) proven methodologies for genetic studies (polygenic trait analyses, mapping tools, genomic sequencing, methods for gene manipulation). Initial investigations of inbred strains of mice showed that femoral and lumbar vertebral volumetric bone mineral density (BMD, mg/mm(3)) by pQCT varied in excess of 50% for femurs and 9% in vertebral BMD. Two strains - low BMD C57BL/6J (B6) mice and high BMD C3H/HeJ (C3H) - were investigated for insights to their BMD diversity. B6C3F2 females derived from intercrossing B6C3F1s were raised to adult skeletal status at 4 months, then necropsied for phenotyping of bone and genotyping of genomic DNA. 1000 F2 females were genotyped for PCR product polymorphisms on all 19 autosomes at approximately 15 cM. Genome wide analyses for genotype-phenotype correlations showed 10 chromosomes (Chrs) carried genes for femoral and 7 Chrs for vertebral BMD. LOD scores ranged from 2.90 to 24.4, and percent of F2 variance accounted for ranged from 1 to 10%. Analyses of main effects revealed both dominant-recessive and additive inheritance patterns. Both progenitor strains carried alleles with positive and negative effects on BMD of each bone sites. A remarkable array of additonal skeletal phenotypes (femur and vertebral geometry, strength measures, serum markers) also proved polygenic in nature, with complex segregation patterns. Verification of BMD quantitative trait loci (QTLs) was undertaken by creating congenic B6 strains carrying individual QTL regions from C3H. Following 6 cycles of backcrossing a QTL-containing region from C3H to the B6 strain, N6F2 congenic strain mice were aged to 4 months, then genotyped for the QTL region and phenotyped for skeletal traits. Comparison of mice homozygous for C3H alleles versus homozygous for B6 alleles in the QTL regions showed that femoral BMD increased or decreased significantly in congenic strains, as was predicted from F2 data. Gender differences specific to BMD QTLs have been revealed, as have more than 30 additional phenotypes associated with cortical and trabecular structural parameters and biomechanical properties.     

Dissection of the pathway of molecular recognition by calmodulin

Amide hydrogen exchange has been used to examine the structural dynamics and energetics of the interaction of a peptide corresponding to the calmodulin-binding domain of smooth muscle myosin light chain kinase (smMLCKp) with calcium-saturated calmodulin. Heteronuclear NMR (15)N-(1)H correlation spectroscopy was used to quantify amide proton exchange rates of the uniformly (15)N-labeled domain bound to calmodulin. A key feature of a proposed model for molecular recognition by calmodulin [Ehrhardt et al. (1995) Biochemistry 34, 2731-2738] is tested by examination of the dependence of amide hydrogen exchange on applied hydrostatic pressure. Hydrogen exchange rates and corresponding protection factors (1/K(op)) for individual amide protons of the bound smMLCKp domain span 5 orders of magnitude at ambient pressure. Individual protection factors decrease significantly in a linear fashion with increasing hydrostatic pressure. A common pressure dependence is revealed by a constant large negative volume change across the residues comprising the core of the bound helical domain. The pattern of protection factors and their response to hydrostatic pressure is consistent with a structural reorganization that results in the concerted disruption of ion pairs between calmodulin and the bound domain. These observations reinforce a model for the molecular recognition pathway where formation of the initial encounter complex is followed by helix-coil transitions in the bound state and subsequent concerted formation of the extensive ion pair network defining the intermolecular contact surface between CaM and the target domain in the final, compact complex structure.     

Using zebrafish as the model organism to understand organ regeneration

The limited regenerative capacity of several organs, such as central nervous system(CNS), heart and limb in mammals makes related major diseases quite difficult to recover. Therefore, dissection of the cellular and molecular mechanisms underlying organ regeneration is of great scientific and clinical interests. Tremendous progression has already been made after extensive investigations using several model organisms for decades. Unfortunately, distance to the final achievement of the goal still remains. Recently, zebrafish became a popular model organism for the deep understanding of regeneration based on its powerful regenerative capacity, in particular the organs that are limitedly regenerated in mammals. Additionally, zebrafish are endowed with other advantages good for the study of organ regeneration. This review summarizes the recent progress in the study of zebrafish organ regeneration, in particular regeneration of fin, heart, CNS, and liver as the representatives. We also discuss reasons of the reduced regenerative capacity in higher vertebrate, the roles of inflammation during regeneration, and the difference between organogenesis and regeneration.     

Engineering a peptide aptamer to target calmodulin for the inhibition of Magnaporthe oryzae

To develop an antimicrobial agent for preventing the devasting damage caused by rice blast, a novel peptide aptamer was identified to interact with calmodulin (CaM) for the inhibition of the spore development in the pathogen Magnaporthe oryzae. A peptide aptamer designated as SNP-D4, consisted of the scaffold protein Staphylococcus aureus nuclease (SN) and an exposed surface loop of 16 random amino acids, was screened from the constructed peptide aptamer libraries by bacterial two-hybrid system using CaM of M. oryzae as the bait. The preliminary inhibition in the sporulation development was observed after treating with the crude extracts expressing SNP-D4. The inhibition efficacies of the purified SNP-D4 were quantified at the stages of conidial germination, germ tube elongation, and appressorium formation in M. oryzae. The binding affinity analysis revealed that SNP-D4 interacted with CaM at a dissociation constant (K_d) of about 20 μM. Moreover, the N-terminus of CaM was identified as the key binding region.     

Studies of the reduction and protonation behavior of tetraheme cytochromes using atomic detail

A comparative study of tetraheme cytochrome c3 molecules from several species was carried out using recently developed theoretical methods based on continuum electrostatics. The binding joint equilibrium of electrons and protons was simulated, revealing the complete thermodynamic aspects of electron-proton coupling in these molecules. The method yields excellent accuracy in terms of midpoint potentials, giving the correct reduction orders in all molecules examined, except for one heme site. The coupling between electrons and protons is shown to be present and significant at physiological pH in all cases. This phenomenon, known as the redox-Bohr effect, though of thermodynamic nature, is shown to have an intrinsic "dynamic" character at the molecular level (in the sense of the empty/occupied fluctuations at the microscopic level), with the binding states of redox and protonatable sites displaying both correlated averages and correlated fluctuations. The protonatable sites more directly involved in the redox-Bohr effect are identified using, among other properties, the statistical correlation between pairs of sites, which automatically reflects indirect effects mediated by other sites. Several sites are identified in this analysis. Propionate D of heme I seems to be the most interesting, generally showing a high correlation not only with its own heme, but also with heme II, corresponding to an indirect stabilization of the reduced forms of both hemes. Other interesting sites are the free histidines of two of the cytochromes and propionate D of heme IV, the latter being potentially associated with redox-induced structural changes. Among the set of cytochromes c3 analyzed in this study, significant differences are observed for several properties of the acidic cytochrome included in the set, from Desulfovibrio africanus, supporting the hypothesis of a different functional role.