Novel approach for structural identification of protein family: glyoxalase I

Glyoxalase is one of two enzymes of the glyoxalase detoxification system against methylglyoxal and other aldehydes, the metabolites derived from glycolysis. The glyoxalase system is found almost in all living organisms: bacteria, protozoa, plants, and animals, including humans, and is related to the class of ‘life essential proteins’. The enzyme belongs to the expanded Glyoxalase/Bleomycin resistance protein/Dioxygenase superfamily. At present the GenBank contains about 700 of amino acid sequences of this enzyme type, and the Protein Data Bank includes dozens of spatial structures. We have offered a novel approach for structural identification of glyoxalase I protein family, which is based on the selecting of basic representative proteins with known structures. On this basis, six new subfamilies of these enzymes have been derived. Most populated subfamilies A1 and A2 were based on representative human Homo sapiens and bacterial Escherichia coli enzymes. We have found that the principle feature, which defines the subfamilies’ structural differences, is conditioned by arrangement of N- and C-domains inside the protein monomer. Finely, we have deduced the structural classification for the glyoxalase I and assigned about 460 protein sequences distributed among six new subfamilies. Structural similarities and specific differences of all the subfamilies have been presented. This approach can be used for structural identification of thousands of the so-called hypothetical proteins with the known PDB structures allowing to identify many of already existing atomic coordinate entrees.     

The PHY Domain Dimer Interface of Bacteriophytochromes Mediates Cross-talk between Photosensory Modules and Output Domains

Protein dynamics play a major role for the catalytic function of enzymes, the interaction of protein complexes or signal integration in regulatory proteins. In the context of multi-domain proteins involved in light-regulation of enzymatic effectors, the central role of conformational dynamics is well established. Light activation of sensory modules is followed by long-range signal transduction to different effectors; rather than domino-style structural rearrangements, a complex interplay of functional elements is required to maintain functionality. One family of such sensor-effector systems are red-light-regulated phytochromes that control diguanylate cyclases involved in cyclic-dimeric-GMP formation. Based on structural and functional studies of one prototypic family member, the central role of the coiled-coil sensor-effector linker was established. Interestingly, subfamilies with different linker lengths feature strongly varying biochemical characteristics. The dynamic interplay of the domains involved, however, is presently not understood. Here we show that the PHY domain dimer interface plays an essential role in signal integration, and that a functional coupling with the coiled-coil linker element is crucial. Chimaeras of two biochemically different family members highlight the phytochrome-spanning helical spine as an essential structural element involved in light-dependent upregulation of enzymatic turnover. However, isolated structural elements can frequently not be assigned to individual characteristics, which further emphasises the importance of global conformational dynamics. Our results provide insights into the intricate processes at play during light signal integration and transduction in these photosensory systems and thus provide additional guidelines for a more directed design of novel sensor-effector combinations with potential applications as optogenetic tools.     

The integrase family of tyrosine recombinases: evolution of a conserved active site domain.

The integrases are a diverse family of tyrosine recombinases which rearrange DNA duplexes by means of conservative site-specific recombination reactions. Members of this family, of which the well-studied lambda Int protein is the prototype, were previously found to share four strongly conserved residues, including an active site tyrosine directly involved in transesterification. However, few additional sequence similarities were found in the original group of 27 proteins. We have now identified a total of 81 members of the integrase family deposited in the databases. Alignment and comparisons of these sequences combined with an evolutionary analysis aided in identifying broader sequence similarities and clarifying the possible functions of these conserved residues. This analysis showed that members of the family aggregate into subfamilies which are consistent with their biological roles; these subfamilies have significant levels of sequence similarity beyond the four residues previously identified. It was also possible to map the location of conserved residues onto the available crystal structures; most of the conserved residues cluster in the predicted active site cleft. In addition, these results offer clues into an apparent discrepancy between the mechanisms of different subfamilies of integrases.     

Functional diversity among spectrin isoforms

The purpose of this review on spectrin is to examine the functional properties of this ubiquitous family of membrane skeletal proteins. Major topics include spectrin-membrane linkages, spectrin-filament linkages, the subcellular localization of spectrins in various cell types and a discussion of major functional differences between erythroid and nonerythroid spectrins. This includes a summary of studies from our own laboratories on the functional and structural comparison of avian spectrin isoforms which are comprised of a common alpha subunit and a tissue-specific beta subunit. Consequently, the observed differences among these spectrins can be assigned to differences in the properties of the beta subunits.     

Cloning and sequencing of cDNA clones encoding chicken lamins A and B1 and comparison of the primary structures of vertebrate A- and B-type lamins

Nuclear lamins are intermediate-filament-type proteins forming a fibrillar meshwork underlying the inner nuclear membrane. The existence of multiple isoforms of lamin proteins in vertebrates is believed to reflect functional specializations during cell division and differentiation. Although biochemical criteria may be used to classify many lamin isoforms into A- and B-type subfamilies, the structural features distinguishing the members of these subfamilies remain to be characterized fully. Here, we report the complete primary structures of chicken lamins A and B1, as they are deduced from cloned cDNAs; in the accompanying paper we present the complete sequence of lamin B2, a second avian B-type lamin. Comparisons of the chicken lamin sequences with each other and with those of other lamins allow us to establish structural features that are common to members of both subfamilies. Conversely, multiple sequence alignments make it possible to identify a number of structural motifs that clearly differentiate B-type lamins from A-type lamins. With this information at hand, we attempt to correlate different biochemical properties of A- and B-type lamins with the presence or absence of specific sequence motifs.     

Enhanced functional and structural domain assignments using remote similarity detection procedures for proteins encoded in the genome of<Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis> H37Rv

The sequencing of theMycobacterium tuberculosis (MTB) H37Rv genome has facilitated deeper insights into the biology of MTB, yet the functions of many MTB proteins are unknown. We have used sensitive profile-based search procedures to assign functional and structural domains to infer functions of gene products encoded in MTB. These domain assignments have been made using a compendium of sequence and structural domain families. Functions are predicted for 78% of the encoded gene products. For 69% of these, functions can be inferred by domain assignments. The functions for the rest are deduced from their homology to proteins of known function. Superfamily relationships between families of unknown and known structures have increased structural information by ∼ 11%. Remote similarity detection methods have enabled domain assignments for 1325 ‘hypothetical proteins’. The most populated families in MTB are involved in lipid metabolism, entry and survival of the bacillus in host. Interestingly, for 353 proteins, which we refer to as MTB-specific, no homologues have been identified. Numerous, previously unannotated, hypothetical proteins have been assigned domains and some of these could perhaps be the possible chemotherapeutic targets. MTB-specific proteins might include factors responsible for virulence. Importantly, these assignments could be valuable for experimental endeavors. The detailed results are publicly available at http://hodgkin.mbu.iisc.ernet.in/∼dots. An erratum to this article is available at .     

Characterization of ligand-induced conformational states in the beta 2 adrenergic receptor

Drugs acting at G protein coupled receptors can be classified in biological assays as either agonists, partial agonists, neutral antagonists, or as inverse agonists. Very little is known about the actual molecular events and structural changes that occur in the receptor following ligand binding and during transmission of a signal across the membrane. Therefore, the structural basis for the biological classification of drug action remains unknown. To date, the conformational state of G protein coupled receptors has been inferred from the activity of the effector enzyme modulated by the G protein. We have used two different approaches to monitor conformational changes in beta 2 adrenergic receptor. Fluorescence spectroscopy can be used to directly monitor structural changes in purified beta 2 adrenergic receptor in real-time. The emission from many fluorescent molecules is strongly dependent on the polarity of the environment in which they are located. Thus, fluorescent probes covalently bound to proteins can be used as sensitive indicators of conformational changes and protein-protein interactions. In addition, we examined functional differences between agonists and partial agonists using fusion proteins between wild-type beta 2 receptor or a constitutively active beta 2 receptor mutant and Gs alpha. These receptor-G protein fusion proteins guarantee highly efficient coupling with a defined stoichiometry. The results of these experiments will be discussed in the context of current models of G protein coupled receptor activation.     

Analysis of the peroxiredoxin family: using active-site structure and sequence information for global classification and residue analysis

Peroxiredoxins (Prxs) are a widespread and highly expressed family of cysteine‐based peroxidases that react very rapidly with H₂O₂, organic peroxides, and peroxynitrite. Correct subfamily classification has been problematic because Prx subfamilies are frequently not correlated with phylogenetic distribution and diverge in their preferred reductant, oligomerization state, and tendency toward overoxidation. We have developed a method that uses the Deacon Active Site Profiler (DASP) tool to extract functional‐site profiles from structurally characterized proteins to computationally define subfamilies and to identify new Prx subfamily members from GenBank(nr). For the 58 literature‐defined Prx test proteins, 57 were correctly assigned, and none were assigned to the incorrect subfamily. The >3500 putative Prx sequences identified were then used to analyze residue conservation in the active site of each Prx subfamily. Our results indicate that the existence and location of the resolving cysteine vary in some subfamilies (e.g., Prx5) to a greater degree than previously appreciated and that interactions at the A interface (common to Prx5, Tpx, and higher order AhpC/Prx1 structures) are important for stabilization of the correct active‐site geometry. Interestingly, this method also allows us to further divide the AhpC/Prx1 into four groups that are correlated with functional characteristics. The DASP method provides more accurate subfamily classification than PSI‐BLAST for members of the Prx family and can now readily be applied to other large protein families. Proteins 2011. © 2010 Wiley‐Liss, Inc.     

Can functional regions of proteins be predicted from their coding sequences? The case study of G-protein coupled receptors

A filter based on a set of unsupervised neural networks trained with a winner-take-all strategy discloses signals along the coding sequences of G-protein coupled receptors. By comparing with the existing experimental data it appears that these signals correlate with putative functional domains of the proteins. After protein alignment within subfamilies, signals cluster in protein regions which, according to the presently available experimental results, are described as possible functional domains of the folded proteins. The mapping procedure reveals characteristic regions in the coding sequences common and/or characteristic of the receptor subtype. This is particularly noticeable for the third cytoplasmic loop, which is likely to be involved in the molecular coupling of all the subfamilies with G-proteins. The results indicate that our mapping can highlight intrinsic representative features of the coding sequences which, in the case of G-protein coupled receptors, are characteristic of protein functional regions and suggest a possible application of the filter for predicting functional determinants in proteins starting from the coding sequence.     

Similarity of different β-strands flanked in loops by glycines and prolines from distinct (α/β)8-barrel enzymes: Chance or a homology?

Many (alpha/beta)8-barrel enzymes contain their conserved sequence regions at or around the beta-strand segments that are often preceded and succeeded by glycines and prolines, respectively. alpha-Amylase is one of these enzymes. Its sequences exhibit a very low degree of similarity, but strong conservation is seen around its beta-strands. These conserved regions were used in the search for similarities with beta-strands of other (alpha/beta)8-barrel enzymes. The analysis revealed an interesting similarity between the segment around the beta 2-strand of alpha-amylase and the one around the beta 4-strand of glycolate oxidase that are flanked in loops by glycines and prolines. The similarity can be further extended on other members of the alpha-amylase and glycolate oxidase subfamilies, i.e., cyclodextrin glycosyltransferase and oligo-1,6-glucosidase, and flavocytochrome b2, respectively. Moreover, the alpha-subunit of tryptophan synthase, the (alpha/beta)8-barrel enzyme belonging to the other subfamily of (alpha/beta)8-barrels, has both investigated strands, beta 2 and beta 4, similar to beta 2 of alpha-amylase and beta 4 of glycolate oxidase. The possibilities of whether this similarity exists only by chance or is a consequence of some processes during the evolution of (alpha/beta)8-barrel proteins are briefly discussed.     

Solution structure of the RWD domain of the mouse GCN2 protein

GCN2 is the alpha-subunit of the only translation initiation factor (eIF2alpha) kinase that appears in all eukaryotes. Its function requires an interaction with GCN1 via the domain at its N-terminus, which is termed the RWD domain after three major RWD-containing proteins: RING finger-containing proteins, WD-repeat-containing proteins, and yeast DEAD (DEXD)-like helicases. In this study, we determined the solution structure of the mouse GCN2 RWD domain using NMR spectroscopy. The structure forms an alpha + beta sandwich fold consisting of two layers: a four-stranded antiparallel beta-sheet, and three side-by-side alpha-helices, with an alphabetabetabetabetaalphaalpha topology. A characteristic YPXXXP motif, which always occurs in RWD domains, forms a stable loop including three consecutive beta-turns that overlap with each other by two residues (triple beta-turn). As putative binding sites with GCN1, a structure-based alignment allowed the identification of several surface residues in alpha-helix 3 that are characteristic of the GCN2 RWD domains. Despite the apparent absence of sequence similarity, the RWD structure significantly resembles that of ubiquitin-conjugating enzymes (E2s), with most of the structural differences in the region connecting beta-strand 4 and alpha-helix 3. The structural architecture, including the triple beta-turn, is fundamentally common among various RWD domains and E2s, but most of the surface residues on the structure vary. Thus, it appears that the RWD domain is a novel structural domain for protein-binding that plays specific roles in individual RWD-containing proteins.     

Diversity, taxonomy and evolution of medium-chain dehydrogenase/reductase superfamily.

A comprehensive, structural and functional, in silico analysis of the medium-chain dehydrogenase/reductase (MDR) superfamily, including 583 proteins, was carried out by use of extensive database mining and the blastp program in an iterative manner to identify all known members of the superfamily. Based on phylogenetic, sequence, and functional similarities, the protein members of the MDR superfamily were classified into three different taxonomic categories: (a) subfamilies, consisting of a closed group containing a set of ideally orthologous proteins that perform the same function; (b) families, each comprising a cluster of monophyletic subfamilies that possess significant sequence identity among them and might share or not common substrates or mechanisms of reaction; and (c) macrofamilies, each comprising a cluster of monophyletic protein families with protein members from the three domains of life, which includes at least one subfamily member that displays activity related to a very ancient metabolic pathway. In this context, a superfamily is a group of homologous protein families (and/or macrofamilies) with monophyletic origin that shares at least a barely detectable sequence similarity, but showing the same 3D fold. The MDR superfamily encloses three macrofamilies, with eight families and 49 subfamilies. These subfamilies exhibit great functional diversity including noncatalytic members with different subcellular, phylogenetic, and species distributions. This results from constant enzymogenesis and proteinogenesis within each kingdom, and highlights the huge plasticity that MDR superfamily members possess. Thus, through evolution a great number of taxa-specific new functions were acquired by MDRs. The generation of new functions fulfilled by proteins, can be considered as the essence of protein evolution. The mechanisms of protein evolution inside MDR are not constrained to conserve substrate specificity and/or chemistry of catalysis. In consequence, MDR functional diversity is more complex than sequence diversity. MDR is a very ancient protein superfamily that existed in the last universal common ancestor. It had at least two (and probably three) different ancestral activities related to formaldehyde metabolism and alcoholic fermentation. Eukaryotic members of this superfamily are more related to bacterial than to archaeal members; horizontal gene transfer among the domains of life appears to be a rare event in modern organisms.     

DOM-fold: a structure with crossing loops found in DmpA, ornithine acetyltransferase, and molybdenum cofactor-binding domain

Understanding relationships between sequence, structure, and evolution is important for functional characterization of proteins. Here, we define a novel DOM-fold as a consensus structure of the domains in DmpA (L-aminopeptidase D-Ala-esterase/amidase), OAT (ornithine acetyltransferase), and MocoBD (molybdenum cofactor-binding domain), and discuss possible evolutionary scenarios of its origin. As shown by a comprehensive structure similarity search, DOM-fold distinguished by a two-layered beta/alpha architecture of a particular topology with unusual crossing loops is unique to those three protein families. DmpA and OAT are evolutionarily related as indicated by their sequence, structural, and functional similarities. Structural similarity between the DmpA/OAT superfamily and the MocoBD domains has not been reported before. Contrary to previous reports, we conclude that functional similarities between DmpA/OAT proteins and N-terminal nucleophile (Ntn) hydrolases are convergent and are unlikely to be inherited from a common ancestor.     

Sequence and hydropathy profile analysis of two classes of secondary transporters

A structural class in the MemGen classification of membrane proteins is a set of evolutionary related proteins sharing a similar global fold. A structural class contains both closely related pairs of proteins for which homology is clear from sequence comparison and very distantly related pairs, for which it is not possible to establish homology based on sequence similarity alone. In the latter case the evolutionary link is based on hydropathy profile analysis. Here, we use these evolutionary related sets of proteins to analyze the relationship between E-values in BLAST searches, sequence similarities in multiple sequence alignments and structural similarities in hydropathy profile analyses. Two structural classes of secondary transporters termed ST[3], which includes the Ion Transporter (IT) superfamily and ST[4], which includes the DAACS family (TC# 2.A.23) were extracted from the NCBI protein database. ST[3] contains 2051 unique sequences distributed over 32 families and 59 subfamilies. ST[4] is a smaller class containing 399 unique sequences distributed over 2 families and 7 subfamilies. One subfamily in ST[4] contains a new class of binding protein dependent secondary transporters. Comparison of the averaged hydropathy profiles of the subfamilies in ST[3] and ST[4] revealed that the two classes represent different folds. Divergence of the sequences in ST[4] is much smaller than observed in ST[3], suggesting different constraints on the proteins during evolution. Analysis of the correlation between the evolutionary relationship of pairs of proteins in a class and the BLAST E-value revealed that: (i) the BLAST algorithm is unable to pick up the majority of the links between proteins in structural class ST[3], (ii) ‘low complexity filtering’ and ‘composition based statistics’ improve the specificity, but strongly reduce the sensitivity of BLAST searches for distantly related proteins, indicating that these filters are too stringent for the proteins analyzed, and (iii) the E-value cut-off, which may be used to evaluate evolutionary significance of a hit in a BLAST search is very different for the two structural classes of membrane proteins.     

Dynamics alignment: comparison of protein dynamics in the SCOP database

A novel methodology for comparison of protein dynamics is presented. Protein dynamics is calculated using the Gaussian network model and the modes of motion are globally aligned using the dynamic programming algorithm of Needleman and Wunsch, commonly used for sequence alignment. The alignment is fast and can be used to analyze large sets of proteins. The methodology is applied to the four major classes of the SCOP database: "all alpha proteins," "all beta proteins," "alpha and beta proteins," and "alpha/beta proteins". We show that different domains may have similar global dynamics. In addition, we report that the dynamics of "all alpha proteins" domains are less specific to structural variations within a given fold or superfamily compared with the other classes. We report that domain pairs with the most similar and the least similar global dynamics tend to be of similar length. The significance of the methodology is that it suggests a new and efficient way of mapping between the global structural features of protein families/subfamilies and their encoded dynamics.     

Blue copper proteins: a comparative analysis of their molecular interaction properties

Blue copper proteins are type-I copper-containing redox proteins whose role is to shuttle electrons from an electron donor to an electron acceptor in bacteria and plants. A large amount of experimental data is available on blue copper proteins; however, their functional characterization is hindered by the complexity of redox processes in biological systems. We describe here the application of a semiquantitative method based on a comparative analysis of molecular interaction fields to gain insights into the recognition properties of blue copper proteins. Molecular electrostatic and hydrophobic potentials were computed and compared for a set of 33 experimentally-determined structures of proteins from seven blue copper subfamilies, and the results were quantified by means of similarity indices. The analysis provides a classification of the blue copper proteins and shows that (I) comparison of the molecular electrostatic potentials provides useful information complementary to that highlighted by sequence analysis; (2) similarities in recognition properties can be detected for proteins belonging to different subfamilies, such as amicyanins and pseudoazurins, that may be isofunctional proteins; (3) dissimilarities in interaction properties, consistent with experimentally different binding specificities, may be observed between proteins belonging to the same subfamily, such as cyanobacterial and eukaryotic plastocyanins; (4) proteins with low sequence identity, such as azurins and pseudoazurins, can have sufficient similarity to bind to similar electron donors and acceptors while having different binding specificity profiles.     

Elucidation of subfamily segregation and intramolecular coevolution of the olfactomedin-like proteins by comprehensive phylogenetic analysis and gene expression pattern assessment

The categorization of genes by structural distinctions relevant to biological characteristics is very important for understanding of gene functions and predicting functional implications of uncharacterized genes. It was absolutely necessary to deploy an effective and efficient strategy to deal with the complexity of the large olfactomedin-like (OLF) gene family sharing sequence similarity but playing diversified roles in many important biological processes, as the simple highest-hit homology analysis gave incomprehensive results and led to inappropriate annotation for some uncharacterized OLF members. In light of evolutionary information that may facilitate the classification of the OLF family and proper association of novel OLF genes with characterized homologs, we performed phylogenetic analysis on all 116 OLF proteins currently available, including two novel members cloned by our group. The OLF family segregated into seven subfamilies and members with similar domain compositions or functional properties all fell into relevant subfamilies. Furthermore, our Northern blot analysis and previous studies revealed that the typical human OLF members in each subfamily exhibited tissue-specific expression patterns, which in turn supported the segregation of the OLF subfamilies with functional divergence. Interestingly, the phylogenetic tree topology for the OLF domains alone was almost identical with that of the full-length tree representing the unique phylogenetic feature of full-length OLF proteins and their particular domain compositions. Moreover, each of the major functional domains of OLF proteins kept the same phylogenetic feature in defining similar topology of the tree. It indicates that the OLF domain and the various domains in flanking non-OLF regions have coevolved and are likely to be functionally interdependent. Expanded by a plausible gene duplication and domain couplings scenario, the OLF family comprises seven evolutionarily and functionally distinct subfamilies, in which each member shares similar structural and functional characteristics including the composition of coevolved and interdependent domains. The phylogenetically classified and preliminarily assessed subfamily framework may greatly facilitate the studying on the OLF proteins. Furthermore, it also demonstrated a feasible and reliable strategy to categorize novel genes and predict the functional implications of uncharacterized proteins based on the comprehensive phylogenetic classification of the subfamilies and their relevance to preliminary functional characteristics.     

Sequence and hydropathy profile analysis of two classes of secondary transporters

A structural class in the MemGen classification of membrane proteins is a set of evolutionary related proteins sharing a similar global fold. A structural class contains both closely related pairs of proteins for which homology is clear from sequence comparison and very distantly related pairs, for which it is not possible to establish homology based on sequence similarity alone. In the latter case the evolutionary link is based on hydropathy profile analysis. Here, we use these evolutionary related sets of proteins to analyze the relationship between E-values in BLAST searches, sequence similarities in multiple sequence alignments and structural similarities in hydropathy profile analyses. Two structural classes of secondary transporters termed ST[3], which includes the Ion Transporter (IT) superfamily and ST[4], which includes the DAACS family (TC# 2.A.23) were extracted from the NCBI protein database. ST[3] contains 2051 unique sequences distributed over 32 families and 59 subfamilies. ST[4] is a smaller class containing 399 unique sequences distributed over 2 families and 7 subfamilies. One subfamily in ST[4] contains a new class of binding protein dependent secondary transporters. Comparison of the averaged hydropathy profiles of the subfamilies in ST[3] and ST[4] revealed that the two classes represent different folds. Divergence of the sequences in ST[4] is much smaller than observed in ST[3], suggesting different constraints on the proteins during evolution. Analysis of the correlation between the evolutionary relationship of pairs of proteins in a class and the BLAST E-value revealed that: (i) the BLAST algorithm is unable to pick up the majority of the links between proteins in structural class ST[3], (ii) "low complexity filtering" and "composition based statistics" improve the specificity, but strongly reduce the sensitivity of BLAST searches for distantly related proteins, indicating that these filters are too stringent for the proteins analyzed, and (iii) the E-value cut-off, which may be used to evaluate evolutionary significance of a hit in a BLAST search is very different for the two structural classes of membrane proteins.     

The crystal structure of murine coronin-1: a regulator of actin cytoskeletal dynamics in lymphocytes

Mammalian coronin-1 is preferentially expressed in hematopoietic cells and plays a poorly understood role in the dynamic reorganization of the actin cytoskeleton. Sequence analysis of coronin-1 revealed five WD40 repeats that were predicted to form a beta propeller. They are followed by a 130 residue extension and a 30 residue leucine zipper domain that is responsible for multimerization of the protein. Here, we present the crystal structure of murine coronin-1 without the leucine zipper at 1.75 A resolution. Coronin-1 forms a seven-bladed beta propeller composed of the five predicted WD40 repeats and two additional blades that lack any homology to the canonical WD40 motif. The C-terminal extension adopts an extended conformation, packs tightly against the bottom surface of the propeller, and is likely to be required for the structural stability of the propeller. Analysis of charged and conserved surface residues delineate possible binding sites for F-actin on the beta propeller.     

Sequence and structural analysis of COVID-19 E and M proteins with MERS virus E and M proteins—A comparative study

The outbreak of SARS in 2003, MERS in 2012, and now COVID-19 in 2019 has demonstrated that Coronaviruses are capable of causing primary lethal infections in humans, and the pandemic is now a global concern. The COVID-19 belongs to the beta coronavirus family encoding 29 proteins, of which four are structural, the Spike, Membrane, Envelope, and Nucleocapsid proteins. Here we have analyzed and compared the Membrane (M) and Envelope (E) proteins of COVID-19 and MERS with SARS and Bat viruses. The sequence analysis of conserved regions of both E and M proteins revealed that many regions of COVID-19 are similar to Bat and SARS viruses while the MERS virus showed variations. The essential binding motifs found in SARS appeared in COVID-19. Besides, the M protein of COVID-19 showed a distinct serine phosphorylation site in the C-terminal domain, which looked like a catalytic triad seen in serine proteases. A Dileucine motif occurred many times in the sequence of the M protein of all the four viruses compared. Concerning the structural part, the COVID-19 E protein showed more similarity to Bat while MERS shared similarity with the SARS virus. The M protein of both COVID-19 and MERS displayed variations in the structure. The interaction between M and E proteins was also studied to know the additional binding regions. Our study highlights the critical motifs and structural regions to be considered for further research to design better inhibitors for the infection caused by these viruses.