PyETV: a PyMOL evolutionary trace viewer to analyze functional site predictions in protein complexes

MS2/GmAMS1 encodes a bHLH transcription factor important for tapetum degeneration in soybean

KEY MESSAGE

GmAMS1 is the only functional AMS and works with GmTDF1-1 and GmMS3 to orchestrate the tapetum degeneration in soybean. Heterosis could significantly increase the production of major crops as well as soybean [Glycine max (L.) Merr.]. Stable male-sterile/female-fertile mutants including ms2 are useful resources to apply in soybean hybrid production. Here, we identified the detailed mutated sites of two classic mutants ms2 (Eldorado) and ms2 (Ames) in MS2/GmAMS1 via the high-throughput sequencing method. Subsequently, we verified that GmAMS1, a bHLH transcription factor, is the only functional AMS member in soybean through the complementary experiment in Arabidopsis; and elucidated the dysfunction of its homolog GmAMS2 is caused by the premature stop codon in the gene's coding sequence. Further qRT-PCR analysis and protein-protein interaction assays indicated GmAMS1 is required for expressing downstream members in the putative DYT1-TDF1-AMS-MYB80/MYB103/MS188-MS1 cascade module, and might regulate the upstream members in a feedback mechanism. GmAMS1 could interact with GmTDF1-1 and GmMS3 via different region, which contributes to dissect the mechanism in the tapetum degeneration process. Additionally, as a core member in the conserved cascade module controlling the tapetum development and degeneration, AMS is conservatively present in all land plant lineages, implying that AMS-mediated signaling pathway has been established before land plants diverged, which provides further insight into the tapetal evolution.

Conformational Dynamics, Energetics, and the Divergent Evolution of Allosteric Regulation: The Case of the Yeast MAPK Family

Allosteric mechanisms provide finely-tuned control over signalling proteins. Proteins of the same family may share high sequence identity and structural similarity but show distinct traits of allosteric control and evolutionary divergent regulation. Revealing the determinants of such properties may be important to understand the molecular bases of different regulatory pathways. Herein, we investigate whether and how evolutionarily-divergent traits of allosteric regulation in homologous proteins can be decoded in terms of internal dynamics and interaction networks that support functionally oriented conformations. In this framework, we start from the comparative analysis of the dynamics and energetics of the yeast MAP Kinases (MAPKs) Fus3 and Kss1 in their native basins. Importantly, distinctive dynamic and energetic stabilization features emerge, which can be related to the two proteins' differential ability to be phosphorylated and engage with the allosteric activator Ste5. We then expanded our study to other evolutionarily-related MAPKs. We show that the dynamical and energetical traits defining the distinct regulatory profiles of Fus3 and Kss1 can be traced along their evolutionary tree. Overall, our approach is able to reconnect (latent) allostery with the principal elements of protein structural stabilization and dynamics, showing how allosteric regulation was encrypted in MAPKs structure well before Ste5 appearance.

Variants in ZFX are associated with an X-linked neurodevelopmental disorder with recurrent facial gestalt

Pathogenic variants in multiple genes on the X chromosome have been implicated in syndromic and non-syndromic intellectual disability disorders. ZFX on Xp22.11 encodes a transcription factor that has been linked to diverse processes including oncogenesis and development, but germline variants have not been characterized in association with disease. Here, we present clinical and molecular characterization of 18 individuals with germline ZFX variants. Exome or genome sequencing revealed 11 variants in 18 subjects (14 males and 4 females) from 16 unrelated families. Four missense variants were identified in 11 subjects, with seven truncation variants in the remaining individuals. Clinical findings included developmental delay/intellectual disability, behavioral abnormalities, hypotonia, and congenital anomalies. Overlapping and recurrent facial features were identified in all subjects, including thickening and medial broadening of eyebrows, variations in the shape of the face, external eye abnormalities, smooth and/or long philtrum, and ear abnormalities. Hyperparathyroidism was found in four families with missense variants, and enrichment of different tumor types was observed. In molecular studies, DNA-binding domain variants elicited differential expression of a small set of target genes relative to wild-type ZFX in cultured cells, suggesting a gain or loss of transcriptional activity. Additionally, a zebrafish model of ZFX loss displayed an altered behavioral phenotype, providing additional evidence for the functional significance of ZFX. Our clinical and experimental data support that variants in ZFX are associated with an X-linked intellectual disability syndrome characterized by a recurrent facial gestalt, neurocognitive and behavioral abnormalities, and an increased risk for congenital anomalies and hyperparathyroidism.

Molecular design and synthesis of methoxy-substitued spiropyrans with photomodulated NIR-fluorescence

This study focuses on the molecular design and synthesis of salt spiropyrans with near-IR fluorescence. The structure of the obtained compounds was confirmed by NMR, IR and mass spectroscopy. In the course of studying the spectral and photoluminescent characteristics, it was possible to reveal the effect of some substituents in various positions on the properties of spiropyran dyes. Due to the structural similarity of one of the isomers to cyanine dyes, the obtained compounds are of interest as potential fluorescent probes for bioimagimg, in particular, for DNA studies. To reveal their ability of binding to DNA molecules molecular docking was carried out. Toxic effects of compounds demonstrating NIR fluorescence were studied on biofilms, as well as using bacterial lux-biosensors.

CovET: A covariation-evolutionary trace method that identifies protein structure-function modules

Measuring the relative effect that any two sequence positions have on each other may improve protein design or help better interpret coding variants. Current approaches use statistics and machine learning but rarely consider phylogenetic divergences which, as shown by Evolutionary Trace studies, provide insight into the functional impact of sequence perturbations. Here, we reframe covariation analyses in the Evolutionary Trace framework to measure the relative tolerance to perturbation of each residue pair during evolution. This approach (CovET) systematically accounts for phylogenetic divergences: at each divergence event, we penalize covariation patterns that belie evolutionary coupling. We find that while CovET approximates the performance of existing methods to predict individual structural contacts, it performs significantly better at finding structural clusters of coupled residues and ligand binding sites. For example, CovET found more functionally critical residues when we examined the RNA recognition motif and WW domains. It correlates better with large-scale epistasis screen data. In the dopamine D2 receptor, top CovET residue pairs recovered accurately the allosteric activation pathway characterized for Class A G protein-coupled receptors. These data suggest that CovET ranks highest the sequence position pairs that play critical functional roles through epistatic and allosteric interactions in evolutionarily relevant structure-function motifs. CovET complements current methods and may shed light on fundamental molecular mechanisms of protein structure and function.

The 5'UTR of HCoV-OC43 adopts a topologically constrained structure to intrinsically repress translation

The emergence of SARS-CoV-2, which is responsible for the COVID-19 pandemic, has highlighted the need for rapid characterization of viral mechanisms associated with cellular pathogenesis. Viral UTRs represent conserved genomic elements that contribute to such mechanisms. Structural details of most CoV UTRs are not available, however. Experimental approaches are needed to allow for the facile generation of high-quality viral RNA tertiary structural models, which can facilitate comparative mechanistic efforts. By integrating experimental and computational techniques, we herein report the efficient characterization of conserved RNA structures within the 5'UTR of the HCoV-OC43 genome, a lab-tractable model coronavirus. We provide evidence that the 5'UTR folds into a structure with well-defined stem-loops (SLs) as determined by chemical probing and direct detection of hydrogen bonds by NMR. We combine experimental base-pair restraints with global structural information from SAXS to generate a 3D model that reveals that SL1-4 adopts a topologically constrained structure wherein SLs 3 and 4 coaxially stack. Coaxial stacking is mediated by short linker nucleotides and allows SLs 1 to 2 to sample different cojoint orientations by pivoting about the SL3,4 helical axis. To evaluate the functional relevance of the SL3,4 coaxial helix, we engineered luciferase reporter constructs harboring the HCoV-OC43 5'UTR with mutations designed to abrogate coaxial stacking. Our results reveal that the SL3,4 helix intrinsically represses translation efficiency since the destabilizing mutations correlate with increased luciferase expression relative to wildtype without affecting reporter mRNA levels, thus highlighting how the 5'UTR structure contributes to the viral mechanism.

Investigating a Library of Flavonoids as Potential Inhibitors of a Cancer Therapeutic Target MEK2 Using in Silico Methods

The second leading cause of death in the world is cancer. Mitogen-activated protein kinase (MAPK) and extracellular signal-regulated protein kinase (ERK) 1 and 2 (MEK1/2) stand out among the different anticancer therapeutic targets. Many MEK1/2 inhibitors are approved and widely used as anticancer drugs. The class of natural compounds known as flavonoids is well-known for their therapeutic potential. In this study, we focus on discovering novel inhibitors of MEK2 from flavonoids using virtual screening, molecular docking analyses, pharmacokinetic prediction, and molecular dynamics (MD) simulations. A library of drug-like flavonoids containing 1289 chemical compounds prepared in-house was screened against the MEK2 allosteric site using molecular docking. The ten highest-scoring compounds based on docking binding affinity (highest score: -11.3 kcal/mol) were selected for further analysis. Lipinski's rule of five was used to test their drug-likeness, followed by ADMET predictions to study their pharmacokinetic properties. The stability of the best-docked flavonoid complex with MEK2 was examined for a 150 ns MD simulation. The proposed flavonoids are suggested as potential inhibitors of MEK2 and drug candidates for cancer therapy.

Potential Disruption of Systemic Hormone Transport by Tobacco Alkaloids Using Computational Approaches

Tobacco/nicotine is one of the most toxic and addictive substances and continues to pose a significant threat to global public health. The harmful effects of smoking/nicotine affect every system in the human body. Nicotine has been associated with effects on endocrine homeostasis in humans such as the imbalance of gonadal steroid hormones, adrenal corticosteroid hormones, and thyroid hormones. The present study was conducted to characterize the structural binding interactions of nicotine and its three important metabolites, cotinine, trans-3'-hydroxycotinine, and 5'-hydroxycotinine, against circulatory hormone carrier proteins, i.e., sex-hormone-binding globulin (SHBG), corticosteroid-binding globulin (CBG), and thyroxine-binding globulin (TBG). Nicotine and its metabolites formed nonbonded contacts and/or hydrogen bonds with amino acid residues of the carrier proteins. For SHBG, Phe-67 and Met-139 were the most important amino acid residues for nicotine ligand binding showing the maximum number of interactions and maximum loss in ASA. For CBG, Trp-371 and Asn-264 were the most important amino acid residues, and for TBG, Ser-23, Leu-269, Lys-270, Asn-273, and Arg-381 were the most important amino acid residues. Most of the amino acid residues of carrier proteins interacting with nicotine ligands showed a commonality with the interacting residues for the native ligands of the proteins. Taken together, the results suggested that nicotine and its three metabolites competed with native ligands for binding to their carrier proteins. Thus, nicotine and its three metabolites may potentially interfere with the binding of testosterone, estradiol, cortisol, progesterone, thyroxine, and triiodothyronine to their carrier proteins and result in the disbalance of their transport and homeostasis in the blood circulation.

Exploring the Natural Compounds in Flavonoids for Their Potential Inhibition of Cancer Therapeutic Target MEK1 Using Computational Methods

The Mitogen-Activated Protein Kinase (MAPK) signaling pathway plays an important role in cancer cell proliferation and survival. MAPKs' protein kinases MEK1/2 serve as important targets in drug designing against cancer. The natural compounds' flavonoids are known for their anticancer activity. This study aims to explore flavonoids for their inhibition ability, targeting MEK1 using virtual screening, molecular docking, ADMET prediction, and molecular dynamics (MD) simulations. Flavonoids (n = 1289) were virtually screened using molecular docking and have revealed possible inhibitors of MEK1. The top five scoring flavonoids based on binding affinity (highest score for MEK1 is -10.8 kcal/mol) have been selected for further protein-ligand interaction analysis. Lipinski's rule (drug-likeness) and absorption, distribution, metabolism, excretion, and toxicity predictions were followed to find a good balance of potency. The selected flavonoids of MEK1 have been refined with 30 (ns) molecular dynamics (MD) simulation. The five selected flavonoids are strongly suggested to be promising potent inhibitors for drug development as anticancer therapeutics of the therapeutic target MEK1.

Decoding Cancer Variants of Unknown Significance for Helicase-Nuclease-RPA Complexes Orchestrating DNA Repair During Transcription and Replication

All tumors have DNA mutations, and a predictive understanding of those mutations could inform clinical treatments. However, 40% of the mutations are variants of unknown significance (VUS), with the challenge being to objectively predict whether a VUS is pathogenic and supports the tumor or whether it is benign. To objectively decode VUS, we mapped cancer sequence data and evolutionary trace (ET) scores onto crystallography and cryo-electron microscopy structures with variant impacts quantitated by evolutionary action (EA) measures. As tumors depend on helicases and nucleases to deal with transcription/replication stress, we targeted helicase–nuclease–RPA complexes: (1) XPB-XPD (within TFIIH), XPF-ERCC1, XPG, and RPA for transcription and nucleotide excision repair pathways and (2) BLM, EXO5, and RPA plus DNA2 for stalled replication fork restart. As validation, EA scoring predicts severe effects for most disease mutations, but disease mutants with low ET scores not only are likely destabilizing but also disrupt sophisticated allosteric mechanisms. For sites of disease mutations and VUS predicted to be severe, we found strong co-localization to ordered regions. Rare discrepancies highlighted the different survival requirements between disease and tumor mutations, as well as the value of examining proteins within complexes. In a genome-wide analysis of 33 cancer types, we found correlation between the number of mutations in each tumor and which pathways or functional processes in which the mutations occur, revealing different mutagenic routes to tumorigenesis. We also found upregulation of ancient genes including BLM, which supports a non-random and concerted cancer process: reversion to a unicellular, proliferation-uncontrolled, status by breaking multicellular constraints on cell division. Together, these genes and global analyses challenge the binary “driver” and “passenger” mutation paradigm, support a gradient impact as revealed by EA scoring from moderate to severe at a single gene level, and indicate reduced regulation as well as activity. The objective quantitative assessment of VUS scoring and gene overexpression in the context of functional interactions and pathways provides insights for biology, oncology, and precision medicine.

Uncovering DNA-PKcs ancient phylogeny, unique sequence motifs and insights for human disease

DNA-dependent protein kinase catalytic subunit (DNA-PKcs) is a key member of the phosphatidylinositol-3 kinase-like (PIKK) family of protein kinases with critical roles in DNA-double strand break repair, transcription, metastasis, mitosis, RNA processing, and innate and adaptive immunity. The absence of DNA-PKcs from many model organisms has led to the assumption that DNA-PKcs is a vertebrate-specific PIKK. Here, we find that DNA-PKcs is widely distributed in invertebrates, fungi, plants, and protists, and that threonines 2609, 2638, and 2647 of the ABCDE cluster of phosphorylation sites are highly conserved amongst most Eukaryotes. Furthermore, we identify highly conserved amino acid sequence motifs and domains that are characteristic of DNA-PKcs relative to other PIKKs. These include residues in the Forehead domain and a novel motif we have termed YRPD, located in an α helix C-terminal to the ABCDE phosphorylation site loop. Combining sequence with biochemistry plus structural data on human DNA-PKcs unveils conserved sequence and conformational features with functional insights and implications. The defined generally progressive DNA-PKcs sequence diversification uncovers conserved functionality supported by Evolutionary Trace analysis, suggesting that for many organisms both functional sites and evolutionary pressures remain identical due to fundamental cell biology. The mining of cancer genomic data and germline mutations causing human inherited disease reveal that robust DNA-PKcs activity in tumors is detrimental to patient survival, whereas germline mutations compromising function are linked to severe immunodeficiency and neuronal degeneration. We anticipate that these collective results will enable ongoing DNA-PKcs functional analyses with biological and medical implications.

Targeting SARS-CoV-2 Nsp3 macrodomain structure with insights from human poly(ADP-ribose) glycohydrolase (PARG) structures with inhibitors

Arrival of the novel SARS-CoV-2 has launched a worldwide effort to identify both pre-approved and novel therapeutics targeting the viral proteome, highlighting the urgent need for efficient drug discovery strategies. Even with effective vaccines, infection is possible, and at-risk populations would benefit from effective drug compounds that reduce the lethality and lasting damage of COVID-19 infection. The CoV-2 MacroD-like macrodomain (Mac1) is implicated in viral pathogenicity by disrupting host innate immunity through its mono (ADP-ribosyl) hydrolase activity, making it a prime target for antiviral therapy. We therefore solved the structure of CoV-2 Mac1 from non-structural protein 3 (Nsp3) and applied structural and sequence-based genetic tracing, including newly determined A. pompejana MacroD2 and GDAP2 amino acid sequences, to compare and contrast CoV-2 Mac1 with the functionally related human DNA-damage signaling factor poly (ADP-ribose) glycohydrolase (PARG). Previously, identified targetable features of the PARG active site allowed us to develop a pharmacologically useful PARG inhibitor (PARGi). Here, we developed a focused chemical library and determined 6 novel PARGi X-ray crystal structures for comparative analysis. We applied this knowledge to discovery of CoV-2 Mac1 inhibitors by combining computation and structural analysis to identify PARGi fragments with potential to bind the distal-ribose and adenosyl pockets of the CoV-2 Mac1 active site. Scaffold development of these PARGi fragments has yielded two novel compounds, PARG-345 and PARG-329, that crystallize within the Mac1 active site, providing critical structure-activity data and a pathway for inhibitor optimization. The reported structural findings demonstrate ways to harness our PARGi synthesis and characterization pipeline to develop CoV-2 Mac1 inhibitors targeting the ADP-ribose active site. Together, these structural and computational analyses reveal a path for accelerating development of antiviral therapeutics from pre-existing drug optimization pipelines.

Structure and evolutionary trace-assisted screening of a residue swapping the substrate ambiguity and chiral specificity in an esterase

Our understanding of enzymes with high substrate ambiguity remains limited because their large active sites allow substrate docking freedom to an extent that seems incompatible with stereospecificity. One possibility is that some of these enzymes evolved a set of evolutionarily fitted sequence positions that stringently allow switching substrate ambiguity and chiral specificity. To explore this hypothesis, we targeted for mutation a serine ester hydrolase (EH3) that exhibits an impressive 71-substrate repertoire but is not stereospecific (e.e. 50%). We used structural actions and the computational evolutionary trace method to explore specificity-swapping sequence positions and hypothesized that position I244 was critical. Driven by evolutionary action analysis, this position was substituted to leucine, which together with isoleucine appears to be the amino acid most commonly present in the closest homologous sequences (max. identity, ca. 67.1%), and to phenylalanine, which appears in distant homologues. While the I244L mutation did not have any functional consequences, the I244F mutation allowed the esterase to maintain a remarkable 53-substrate range while gaining stereospecificity properties (e.e. 99.99%). These data support the possibility that some enzymes evolve sequence positions that control the substrate scope and stereospecificity. Such residues, which can be evolutionarily screened, may serve as starting points for further designing substrate-ambiguous, yet chiral-specific, enzymes that are greatly appreciated in biotechnology and synthetic chemistry.

Decoding Protein-protein Interactions: An Overview

Drug discovery has focused on the paradigm "one drug, one target" for a long time. However, small molecules can act at multiple macromolecular targets, which serves as the basis for drug repurposing. In an effort to expand the target space, and given advances in X-ray crystallography, protein-protein interactions have become an emerging focus area of drug discovery enterprises. Proteins interact with other biomolecules and it is this intricate network of interactions that determines the behavior of the system and its biological processes. In this review, we briefly discuss networks in disease, followed by computational methods for protein-protein complex prediction. Computational methodologies and techniques employed towards objectives such as protein-protein docking, protein-protein interactions, and interface predictions are described extensively. Docking aims at producing a complex between proteins, while interface predictions identify a subset of residues on one protein that could interact with a partner, and protein-protein interaction sites address whether two proteins interact. In addition, approaches to predict hot spots and binding sites are presented along with a representative example of our internal project on the chemokine CXC receptor 3 B-isoform and predictive modeling with IP10 and PF4.

An efficient chemical screening method for structure-based inhibitors to nucleic acid enzymes targeting the DNA repair-replication interface and SARS CoV-2

We present a Chemistry and Structure Screen Integrated Efficiently (CASSIE) approach (named for Greek prophet Cassandra) to design inhibitors for cancer biology and pathogenesis. CASSIE provides an effective path to target master keys to control the repair-replication interface for cancer cells and SARS CoV-2 pathogenesis as exemplified here by specific targeting of Poly(ADP-ribose) glycohydrolase (PARG) and ADP-ribose glycohydrolase ARH3 macrodomains plus SARS CoV-2 nonstructural protein 3 (Nsp3) Macrodomain 1 (Mac1) and Nsp15 nuclease. As opposed to the classical massive effort employing libraries with large numbers of compounds against single proteins, we make inhibitor design for multiple targets efficient. Our compact, chemically diverse, 5000 compound Goldilocks (GL) library has an intermediate number of compounds sized between fragments and drugs with predicted favorable ADME (absorption, distribution, metabolism, and excretion) and toxicological profiles. Amalgamating our core GL library with an approved drug (AD) library, we employ a combined GLAD library virtual screen, enabling an effective and efficient design cycle of ranked computer docking, top hit biophysical and cell validations, and defined bound structures using human proteins or their avatars. As new drug design is increasingly pathway directed as well as molecular and mechanism based, our CASSIE approach facilitates testing multiple related targets by efficiently turning a set of interacting drug discovery problems into a tractable medicinal chemistry engineering problem of optimizing affinity and ADME properties based upon early co-crystal structures. Optimization efforts are made efficient by a computationally-focused iterative chemistry and structure screen. Thus, we herein describe and apply CASSIE to define prototypic, specific inhibitors for PARG vs distinct inhibitors for the related macrodomains of ARH3 and SARS CoV-2 Nsp3 plus the SARS CoV-2 Nsp15 RNA nuclease.

PyVibMS: a PyMOL plugin for visualizing vibrations in molecules and solids

Visualizing vibrational motions calculated with different ab initio packages requires dedicated post-processing tools. Here, we present a PyMOL plugin called PyVibMS for visualizing the vibrational motions for both molecular and solid systems calculated by mainstream quantum chemical computer programs including Gaussian, Q-Chem, VASP, and CRYSTAL. Benefiting from the continuing development of the PyMOL platform, PyVibMS provides powerful functionalities and user-friendly interface. PyVibMS was written in Python and its open-source nature makes it flexible and sustainable. As an example, the motions of the Konkoli-Cremer local vibrational modes are shown in this work for the first time. PyVibMS is freely available at https://github.com/smutao/PyVibMS . Graphical abstract In this work, a PyMOL plugin named PyVibMS is developed to visualize molecular and lattice vibrations.

Maslinic acid differentially exploits the MAPK pathway in estrogen-positive and triple-negative breast cancer to induce mitochondrion-mediated, caspase-independent apoptosis

Breast cancer accounts for 1.4 million new cases every year. Triple-negative breast cancer (TNBC) is one the leading cause of mortality in developing countries and is associated with early age onset (under 40 years old). Chemotherapy has a poor success rate in patients with TNBC as compared to other types of breast cancers. It is due to the lack of expression of three validated molecular markers for breast cancer, the estrogen and progesterone receptors, and the amplification of HER-2/Neu. Therefore, a clear need exists for a greater understanding of TNBC at all levels and for the development of better therapies. We have studied the anti-tumor effects of a potential drug, maslinic acid, which can be extracted from olive oil industry waste. This natural product showed inhibitory effect at concentrations ranging from 30 to 50 µM within 24 h. It exhibited divergent effects in cell cycle progression for the MCF7 (estrogen positive) cell line when compared with TNBCs like MDA-MB-231 and MDA-MB-468. Also, maslinic acid treatment altered the mitochondrial membrane electrochemical potential and the reactive oxygen species (ROS) levels to cause a caspase-independent programmed cell death. In silico approaches and immunoblotting suggested the involvement of the MAPK pathway explaining the variability in cell cycle progression along with the apoptotic cell death caused by maslinic acid.

iPBAvizu: a PyMOL plugin for an efficient 3D protein structure superimposition approach

Background

Protein 3D structure is the support of its function. Comparison of 3D protein structures provides insight on their evolution and their functional specificities and can be done efficiently via protein structure superimposition analysis. Multiple approaches have been developed to perform such task and are often based on structural superimposition deduced from sequence alignment, which does not take into account structural features. Our methodology is based on the use of a Structural Alphabet (SA), i.e. a library of 3D local protein prototypes able to approximate protein backbone. The interest of a SA is to translate into 1D sequences into the 3D structures.

Results

We used Protein blocks (PB), a widely used SA consisting of 16 prototypes, each representing a conformation of the pentapeptide skeleton defined in terms of dihedral angles. Proteins are described using PB from which we have previously developed a sequence alignment procedure based on dynamic programming with a dedicated PB Substitution Matrix. We improved the procedure with a specific two-step search: (i) very similar regions are selected using very high weights and aligned, and (ii) the alignment is completed (if possible) with less stringent parameters. Our approach, iPBA, has shown to perform better than other available tools in benchmark tests. To facilitate the usage of iPBA, we designed and implemented iPBAvizu, a plugin for PyMOL that allows users to run iPBA in an easy way and analyse protein superimpositions.

Conclusions

iPBAvizu is an implementation of iPBA within the well-known and widely used PyMOL software. iPBAvizu enables to generate iPBA alignments, create and interactively explore structural superimposition, and assess the quality of the protein alignments.

Evolutionary Stability of Salt Bridges Hints Its Contribution to Stability of Proteins

The contribution of newly designed salt bridges to protein stabilization remains controversial even today. In order to solve this problem, we investigated salt bridges from two aspects: spatial distribution and evolutionary characteristics of salt bridges. Firstly, we analyzed spatial distribution of salt bridges in proteins, elucidating the basic requirements of forming salt bridges. Then, from an evolutionary point of view, the evolutionary characteristics of salt bridges as well as their neighboring residues were investigated in our study. The results demonstrate that charged residues appear more frequently than other neutral residues at certain positions of sequence even under evolutionary pressure, which are able to form electrostatic interactions that could increase the evolutionary stability of corresponding amino acid regions, enhancing their importance to stability of proteins. As a corollary, we conjectured that the newly designed salt bridges with more contribution to proteins, not only, are qualified spatial distribution of salt bridges, but also, are needed to further increase the evolutionary stability of corresponding amino acid regions. Based on analysis, the 8 mutations were accordingly constructed in the 1,4-α-glucan branching enzyme (EC 2.4.1.18, GBE) from Geobacillus thermoglucosidans STB02, of which 7 mutations improved thermostability of GBE. The enhanced thermostability of 7 mutations might be a result of additional salt bridges on residue positions that at least one of amino acids positions is conservative, improving their contribution of stabilization to proteins.

FARS2 deficiency; new cases, review of clinical, biochemical, and molecular spectra, and variants interpretation based on structural, functional, and evolutionary significance

An increasing number of mitochondrial diseases are found to be caused by pathogenic variants in nuclear encoded mitochondrial aminoacyl-tRNA synthetases. FARS2 encodes mitochondrial phenylalanyl-tRNA synthetase (mtPheRS) which transfers phenylalanine to its cognate tRNA in mitochondria. Since the first case was reported in 2012, a total of 21 subjects with FARS2 deficiency have been reported to date with a spectrum of disease severity that falls between two phenotypes; early onset epileptic encephalopathy and a less severe phenotype characterized by spastic paraplegia. In this report, we present an additional 15 individuals from 12 families who are mostly Arabs homozygous for the pathogenic variant Y144C, which is associated with the more severe early onset phenotype. The total number of unique pathogenic FARS2 variants known to date is 21 including three different partial gene deletions reported in four individuals. Except for the large deletions, all variants but two (one in-frame deletion of one amino acid and one splice-site variant) are missense. All large deletions and the single splice-site variant are in trans with a missense variant. This suggests that complete loss of function may be incompatible with life. In this report, we also review structural, functional, and evolutionary significance of select FARS2 pathogenic variants reported here.

Edema Induced by a Crotalus durissus terrificus Venom Serine Protease (Cdtsp 2) Involves the PAR Pathway and PKC and PLC Activation

Snake venom serine proteases (SVSPs) represent an essential group of enzymatic toxins involved in several pathophysiological effects on blood homeostasis. Some findings suggest the involvement of this class of enzymatic toxins in inflammation. In this paper, we purified and isolated a new gyroxin isoform from the Crotalus durissus terrificus (Cdt) venom, designated as Cdtsp 2, which showed significant proinflammatory effects in a murine model. In addition, we performed several studies to elucidate the main pathway underlying the edematogenic effect induced by Cdtsp 2. Enzymatic assays and structural analysis (primary structure analysis and three-dimensional modeling) were closely performed with pharmacological assays. The determination of edematogenic activity was performed using Cdtsp 2 isolated from snake venom, and was applied to mice treated with protein kinase C (PKC) inhibitor, phospholipase C (PLC) inhibitor, dexamethasone (Dexa), antagonists for protease-activated receptors (PARs), or saline (negative control). Additionally, we measured the levels of cyclooxygenase 2 (COX-2), malondialdehyde (MDA), and prostaglandin E2 (PGE2). Cdtsp 2 is characterized by an approximate molecular mass of 27 kDa, an isoelectric point (pI) of 4.5, and significant fibrinolytic activity, as well as the ability to hydrolyze Nα-benzoyl-l-arginine 4-nitroanilide (BAPNA). Its primary and three-dimensional structures revealed Cdtsp 2 as a typical snake venom serine protease that induces significant edema via the metabolism of arachidonic acid (AA), involving PARs, PKC, PLC, and COX-2 receptors, as well as inducing a significant increase in MDA levels. Our results showed that Cdtsp 2 is a serine protease with significant enzymatic activity, and it may be involved in the degradation of PAR1 and PAR2, which activate PLC and PKC to mobilize AA, while increasing oxidative stress. In this article, we provide a new perspective for the role of SVSPs beyond their effects on blood homeostasis.

Combinatorial inhibition of PTPN12-regulated receptors leads to a broadly effective therapeutic strategy in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is an aggressive subtype of breast cancer diagnosed in more than 200,000 women each year and is recalcitrant to targeted therapies. Although TNBCs harbor multiple hyperactive receptor tyrosine kinases (RTKs), RTK inhibitors have been largely ineffective in TNBC patients thus far. We developed a broadly effective therapeutic strategy for TNBC that is based on combined inhibition of receptors that share the negative regulator PTPN12. Previously, we and others identified the tyrosine phosphatase PTPN12 as a tumor suppressor that is frequently inactivated in TNBC. PTPN12 restrains several RTKs, suggesting that PTPN12 deficiency leads to aberrant activation of multiple RTKs and a co-dependency on these receptors. This in turn leads to the therapeutic hypothesis that PTPN12-deficient TNBCs may be responsive to combined RTK inhibition. However, the repertoire of RTKs that are restrained by PTPN12 in human cells has not been systematically explored. By methodically identifying the suite of RTK substrates (MET, PDGFRβ, EGFR, and others) inhibited by PTPN12, we rationalized a combination RTK-inhibitor therapy that induced potent tumor regression across heterogeneous models of TNBC. Orthogonal approaches revealed that PTPN12 was recruited to and inhibited these receptors after ligand stimulation, thereby serving as a feedback mechanism to limit receptor signaling. Cancer-associated mutation of PTPN12 or reduced PTPN12 protein levels diminished this feedback mechanism, leading to aberrant activity of these receptors. Restoring PTPN12 protein levels restrained signaling from RTKs, including PDGFRβ and MET, and impaired TNBC survival. In contrast with single agents, combined inhibitors targeting the PDGFRβ and MET receptors induced the apoptosis in TNBC cells in vitro and in vivo. This therapeutic strategy resulted in tumor regressions in chemo-refractory patient-derived TNBC models. Notably, response correlated with PTPN12 deficiency, suggesting that impaired receptor feedback may establish a combined addiction to these proto-oncogenic receptors. Taken together, our data provide a rationale for combining RTK inhibitors in TNBC and other malignancies that lack receptor-activating mutations.

S18 family of mitochondrial ribosomal proteins: evolutionary history and Gly132 polymorphism in colon carcinoma

S18 family of mitochondrial ribosomal proteins (MRPS18, S18) consists of three members, S18-1 to -3. Earlier, we found that overexpression of S18-2 protein resulted in immortalization and eventual transformation of primary rat fibroblasts. The S18-1 and -3 have not exhibited such abilities. To understand the differences in protein properties, the evolutionary history of S18 family was analyzed. The S18-3, followed by S18-1 and S18-2 emerged as a result of ancient gene duplication in the root of eukaryotic species tree, followed by two metazoan-specific gene duplications. However, the most conserved metazoan S18 homolog is the S18-1; it shares the most sequence similarity with S18 proteins of bacteria and of other eukaryotic clades. Evolutionarily conserved residues of S18 proteins were analyzed in various cancers. S18-2 is mutated at a higher rate, compared with S18-1 and -3 proteins. Moreover, the evolutionarily conserved residue, Gly132 of S18-2, shows genetic polymorphism in colon adenocarcinomas that was confirmed by direct DNA sequencing.Concluding, S18 family represents the yet unexplored important mitochondrial ribosomal proteins.

Peptide probes derived from pertuzumab by molecular dynamics modeling for HER2 positive tumor imaging

A high level of HER2 expression in breast cancer correlates with a higher tumor growth rate, high metastatic potential, and a poor long-term patient survival rate. Pertuzumab, a human monoclonal antibody, can reduce the effect of HER2 overexpression by preventing HER2 dimerization. In this study, a combination protocol of molecular dynamics modeling and MM/GBSA binding free energy calculations was applied to design peptides that interact with HER2 based on the HER2/pertuzumab crystal structure. Based on a β hairpin in pertuzumab from Glu46 to Lys65—which plays a key role in interacting with HER2—mutations were carried out in silico to improve the binding free energy of the hairpin that interacts with the Phe256-Lys314 of the HER2 protein. Combined the use of one-bead-one-compound library screening, among all the mutations, a peptide (58F63Y) with the lowest binding free energy was confirmed experimentally to have the highest affinity, and it may be used as a new probe in diagnosing and treating HER2-positive breast cancer.

Many therapeutic approaches, including the human monoclonal antibodies trastuzumab and pertuzumab, target the human epidermal growth factor receptor 2 (HER2) of any breast cancer that features HER2 overexpression. Compared to these antibodies, peptides have many advantages, including lower cost, easier synthesis, high affinity, and lower toxicity. Here, we first designed peptides that interact with HER2 protein based on the HER2/pertuzumab crystal structure (PDB entry: 1S78), using a combination protocol of molecular dynamics modeling, molecular mechanics/generalized Born solvent-accessible surface area (MM/GBSA) binding free energy calculations. Then, combined with the peptide library screening, six peptides were selected for further analysis and experimental validations. The results of ex vivo and in vivo experiments confirmed that one peptide (58F63Y) in particular has a strong affinity and high specificity to HER2-overexpressing tumors. This may due to more paired residues and lower binding free energy in peptide 58F63Y and HER2 complex based on free energy decomposition analysis and distances calculation. While both in silico and in vitro screenings point to the same high-affinity peptide, the findings suggest that in silico screening based on calculated binding free energies is rather reliable. Additionally, based on the calculation of binding free energies among mutants, we can reduce the library capacity of one-bead-one-compound screening. In summary, we present a rather simple and rapid means of deriving a peptide with a clear binding site to its target protein.

Structural Studies on the Extracellular Domain of Sensor Histidine Kinase YycG from Staphylococcus aureus and Its Functional Implications

Bacterial two-component signal transduction systems are used to adapt to fluctuations in the environment. YycG, a key two-component histidine kinase in Staphylococcus aureus, plays an essential role in cell viability and regulates cell wall metabolism, biofilm formation, virulence, and antibiotic resistance. For these reasons, YycG is considered a compelling target for the development of novel antibiotics. However, to date, the signaling mechanism of YycG and its stimulus are poorly understood mainly because of a lack of structural information on YycG. To address this deficiency, we determined the crystal structure of the extracellular domain of S. aureus YycG (YycGex) at 2.0-Å resolution. The crystal structure indicated two subunits with an extracellular Per-Arnt-Sim (PAS) topology packed into a dimer with interloop interactions. Disulfide scanning using cysteine-substituted mutants revealed that YycGex possessed dimeric interfaces not only in the loop but also in the helix α1. Cross-linking studies using intact YycG demonstrated that it was capable of forming high molecular weight oligomers on the cell membrane. Furthermore, we also observed that two auxiliary proteins of YycG, YycH and YycI, cooperatively interfered with the multimerization of YycG. From these results, we propose that signaling through YycG is regulated by multimerization and binding of YycH and YycI. These structural studies, combined with biochemical analyses, provide a better understanding of the signaling mechanism of YycG, which is necessary for developing novel antibacterial drugs targeting S. aureus.

The architecture of the Schizosaccharomyces pombe CCR4-NOT complex

CCR4-NOT is a large protein complex present both in cytoplasm and the nucleus of eukaryotic cells. Although it is involved in a variety of distinct processes related to expression of genetic information such as poly(A) tail shortening, transcription regulation, nuclear export and protein degradation, there is only fragmentary information available on some of its nine subunits. Here we show a comprehensive structural characterization of the native CCR4-NOT complex from Schizosaccharomyces pombe. Our cryo-EM 3D reconstruction of the complex, combined with techniques such as immunomicroscopy, RNA-nanogold labelling, docking of the available high-resolution structures and models of different subunits and domains, allow us to propose its full molecular architecture. We locate all functionally defined domains endowed with deadenylating and ubiquitinating activities, the nucleus-specific RNA-interacting subunit Mmi1, as well as surfaces responsible for protein–protein interactions. This information provides insight into cooperation of the different CCR4-NOT complex functions.

CCR4-NOT is a protein complex involved in a variety of important genetic processes. Here, the authors report the mid-resolution structure of this complex, and model the positions and contacts between the subunits, providing structural support for the previously reported functions of the complex.

UET: a database of evolutionarily-predicted functional determinants of protein sequences that cluster as functional sites in protein structures

The structure and function of proteins underlie most aspects of biology and their mutational perturbations often cause disease. To identify the molecular determinants of function as well as targets for drugs, it is central to characterize the important residues and how they cluster to form functional sites. The Evolutionary Trace (ET) achieves this by ranking the functional and structural importance of the protein sequence positions. ET uses evolutionary distances to estimate functional distances and correlates genotype variations with those in the fitness phenotype. Thus, ET ranks are worse for sequence positions that vary among evolutionarily closer homologs but better for positions that vary mostly among distant homologs. This approach identifies functional determinants, predicts function, guides the mutational redesign of functional and allosteric specificity, and interprets the action of coding sequence variations in proteins, people and populations. Now, the UET database offers pre-computed ET analyses for the protein structure databank, and on-the-fly analysis of any protein sequence. A web interface retrieves ET rankings of sequence positions and maps results to a structure to identify functionally important regions. This UET database integrates several ways of viewing the results on the protein sequence or structure and can be found at http://mammoth.bcm.tmc.edu/uet/.

Identifying a potential receptor for the antibacterial peptide of sponge Axinella donnani endosymbiont

Marine sponges and their associated bacteria are rich sources of novel secondary metabolites with therapeutic usefulness. In our earlier work, we have identified a novel antibacterial peptide from the marine sponge Axinella donnani endosymbiotic bacteria. In this work, we have carried out a comparative genomic analysis and identified a set of 60 proteins as probable receptor which is common in all the strains. The analysis on binding substrate showed that β barrel assembly machinery (BamA) of the outer membrane protein 85 (omp85) superfamily is a potential receptor protein for the antibacterial peptide. It plays a central role in OMP biogenesis, especially in cell viability. Further, the triplet and quartet motifs RGF and YGDG, respectively in L6 loop are conserved over all the strains and these conserved residues interact with antibacterial peptide to inhibit the BamA function, which is essential for OMP biogenesis.

Phylogenetic and 2D/3D Analysis of HCV 1a NS4A Gene/Protein in Pakistani Isolates

BACKGROUND

The nonstructural protein NS4A of hepatitis C virus is composed of 54 amino acids. This small size protein has vital role in many cellular functions. The most important reported function is being a cofactor of viral enzymes serine protease and helicase.

OBJECTIVES

The objective of this study was to analyze the phylogenetic variation, its impact in terms of translation and any functional change in protein structure at primary 2D/3D structure using computational tools from Pakistani patients isolates.

MATERIALS AND METHODS

Patient sera infected with Hepatitis C virus, genotype 1A, were obtained from Molecular Diagnostics lab, CEMB, University of the Punjab Lahore by using BD Vacutainer collection tubes (Becton Dickenson).

RESULTS

Phylogenetic analysis of the gene revealed that Pakistani 1a HCV strains are in the start of third cluster and there is a difference between inter Pakistani isolates at primary, secondary and tertiary levels.

CONCLUSIONS

Mutations were present in the central domain of NS4A (amino acids 21 - 34).

Structural basis for misfolding in myocilin-associated glaucoma

Olfactomedin (OLF) domain-containing proteins play roles in fundamental cellular processes and have been implicated in disorders ranging from glaucoma, cancers and inflammatory bowel disorder, to attention deficit disorder and childhood obesity. We solved crystal structures of the OLF domain of myocilin (myoc-OLF), the best studied such domain to date. Mutations in myoc-OLF are causative in the autosomal dominant inherited form of the prevalent ocular disorder glaucoma. The structures reveal a new addition to the small family of five-bladed β-propellers. Propellers are most well known for their ability to act as hubs for protein-protein interactions, a function that seems most likely for myoc-OLF, but they can also act as enzymes. A calcium ion, sodium ion and glycerol molecule were identified within a central hydrophilic cavity that is accessible via movements of surface loop residues. By mapping familial glaucoma-associated lesions onto the myoc-OLF structure, three regions sensitive to aggregation have been identified, with direct applicability to differentiating between neutral and disease-causing non-synonymous mutations documented in the human population worldwide. Evolutionary analysis mapped onto the myoc-OLF structure reveals conserved and divergent regions for possible overlapping and distinctive functional protein-protein or protein-ligand interactions across the broader OLF domain family. While deciphering the specific normal biological functions, ligands and binding partners for OLF domains will likely continue to be a challenging long-term experimental pursuit, atomic detail structural knowledge of myoc-OLF is a valuable guide for understanding the implications of glaucoma-associated mutations and will help focus future studies of this biomedically important domain family.

Differential effects of collagen prolyl 3-hydroxylation on skeletal tissues

Mutations in the genes encoding cartilage associated protein (CRTAP) and prolyl 3-hydroxylase 1 (P3H1 encoded by LEPRE1) were the first identified causes of recessive Osteogenesis Imperfecta (OI). These proteins, together with cyclophilin B (encoded by PPIB), form a complex that 3-hydroxylates a single proline residue on the α1(I) chain (Pro986) and has cis/trans isomerase (PPIase) activity essential for proper collagen folding. Recent data suggest that prolyl 3-hydroxylation of Pro986 is not required for the structural stability of collagen; however, the absence of this post-translational modification may disrupt protein-protein interactions integral for proper collagen folding and lead to collagen over-modification. P3H1 and CRTAP stabilize each other and absence of one results in degradation of the other. Hence, hypomorphic or loss of function mutations of either gene cause loss of the whole complex and its associated functions. The relative contribution of losing this complex's 3-hydroxylation versus PPIase and collagen chaperone activities to the phenotype of recessive OI is unknown. To distinguish between these functions, we generated knock-in mice carrying a single amino acid substitution in the catalytic site of P3h1 (Lepre1^H662A). This substitution abolished P3h1 activity but retained ability to form a complex with Crtap and thus the collagen chaperone function. Knock-in mice showed absence of prolyl 3-hydroxylation at Pro986 of the α1(I) and α1(II) collagen chains but no significant over-modification at other collagen residues. They were normal in appearance, had no growth defects and normal cartilage growth plate histology but showed decreased trabecular bone mass. This new mouse model recapitulates elements of the bone phenotype of OI but not the cartilage and growth phenotypes caused by loss of the prolyl 3-hydroxylation complex. Our observations suggest differential tissue consequences due to selective inactivation of P3H1 hydroxylase activity versus complete ablation of the prolyl 3-hydroxylation complex.

The prolyl 3-hydroxylase complex serves to hydroxylate a single residue in type I collagen and also serves as a collagen chaperone. The complex is comprised of prolyl 3-hydroxylase 1, cartilage associated protein, and cyclophilin B. Mutations have been identified in the genes encoding the complex members in patients with recessive Osteogenesis Imperfecta. Recent data suggest that prolyl 3-hydroxylation of collagen does not alter the stability of collagen but may rather mediate protein-protein interactions. Additionally, the collagen chaperoning function of the complex is an important rate limiting step in the modification of type I collagen. Irrespective of whether patients with mutations in the genes encoding the members of the prolyl 3-hydroxylase complex have hypomorphic or complete loss of function alleles, either circumstance leads to the loss of both functions of the prolyl 3-hydroxylase complex. Thus, it is unknown how collagen chaperoning and/or hydroxylation affect bone and cartilage homeostasis. In this study, we generated a mouse model lacking the prolyl 3-hydroxylation activity of the complex while maintaining the chaperoning function. We found that the hydroxylase mutant mice have normal cartilage and normal cortical bone but decreased trabecular bone, suggesting that there is a differential requirement for hydroxylation in different tissues.

xPyder: a PyMOL plugin to analyze coupled residues and their networks in protein structures

A versatile method to directly identify and analyze short- or long-range coupled or communicating residues in a protein conformational ensemble is of extreme relevance to achieve a complete understanding of protein dynamics and structural communication routes. Here, we present xPyder, an interface between one of the most employed molecular graphics systems, PyMOL, and the analysis of dynamical cross-correlation matrices (DCCM). The approach can also be extended, in principle, to matrices including other indexes of communication propensity or intensity between protein residues, as well as the persistence of intra- or intermolecular interactions, such as those underlying protein dynamics. The xPyder plugin for PyMOL 1.4 and 1.5 is offered as Open Source software via the GPL v2 license, and it can be found, along with the installation package, the user guide, and examples, at http://linux.btbs.unimib.it/xpyder/.

ETAscape: analyzing protein networks to predict enzymatic function and substrates in Cytoscape

Most proteins lack experimentally validated functions. To address this problem, we implemented the Evolutionary Trace Annotation (ETA) method in the Cytoscape network visualization environment. The result is the ETAscape plugin, which builds a structural genomics network based on local structural and evolutionary similarities among proteins and then globally diffuses known annotations across the resulting network. The plugin displays these novel functional annotations, their confidence, the molecular basis for individual matches and the set of matches that lead to a prediction.

AVAILABILITY

The ETA Network Plugin is available publicly for download at http://mammoth.bcm.tmc.edu/networks/.

IVSPlat 1.0: an integrated virtual screening platform with a molecular graphical interface

BACKGROUND

The virtual screening (VS) of lead compounds using molecular docking and pharmacophore detection is now an important tool in drug discovery. VS tasks typically require a combination of several software tools and a molecular graphics system. Thus, the integration of all the requisite tools in a single operating environment could reduce the complexity of running VS experiments. However, only a few freely available integrated software platforms have been developed.

RESULTS

A free open-source platform, IVSPlat 1.0, was developed in this study for the management and automation of VS tasks. We integrated several VS-related programs into a molecular graphics system to provide a comprehensive platform for the solution of VS tasks based on molecular docking, pharmacophore detection, and a combination of both methods. This tool can be used to visualize intermediate and final results of the VS execution, while also providing a clustering tool for the analysis of VS results. A case study was conducted to demonstrate the applicability of this platform.

CONCLUSIONS

IVSPlat 1.0 provides a plug-in-based solution for the management, automation, and visualization of VS tasks. IVSPlat 1.0 is an open framework that allows the integration of extra software to extend its functionality and modified versions can be freely distributed. The open source code and documentation are available at http://kyc.nenu.edu.cn/IVSPlat/.

Evolutionary trace for prediction and redesign of protein functional sites

The evolutionary trace (ET) is the single most validated approach to identify protein functional determinants and to target mutational analysis, protein engineering and drug design to the most relevant sites of a protein. It applies to the entire proteome; its predictions come with a reliability score; and its results typically reach significance in most protein families with 20 or more sequence homologs. In order to identify functional hot spots, ET scans a multiple sequence alignment for residue variations that correlate with major evolutionary divergences. In case studies this enables the selective separation, recoding, or mimicry of functional sites and, on a large scale, this enables specific function predictions based on motifs built from select ET-identified residues. ET is therefore an accurate, scalable and efficient method to identify the molecular determinants of protein function and to direct their rational perturbation for therapeutic purposes. Public ET servers are located at: http://mammoth.bcm.tmc.edu/.

Separation of recombination and SOS response in Escherichia coli RecA suggests LexA interaction sites

RecA plays a key role in homologous recombination, the induction of the DNA damage response through LexA cleavage and the activity of error-prone polymerase in Escherichia coli. RecA interacts with multiple partners to achieve this pleiotropic role, but the structural location and sequence determinants involved in these multiple interactions remain mostly unknown. Here, in a first application to prokaryotes, Evolutionary Trace (ET) analysis identifies clusters of evolutionarily important surface amino acids involved in RecA functions. Some of these clusters match the known ATP binding, DNA binding, and RecA-RecA homo-dimerization sites, but others are novel. Mutation analysis at these sites disrupted either recombination or LexA cleavage. This highlights distinct functional sites specific for recombination and DNA damage response induction. Finally, our analysis reveals a composite site for LexA binding and cleavage, which is formed only on the active RecA filament. These new sites can provide new drug targets to modulate one or more RecA functions, with the potential to address the problem of evolution of antibiotic resistance at its root.

In eubacteria, genome integrity is in large part orchestrated by RecA, which directly participates in recombination, induction of DNA damage response through LexA repressor cleavage and error-prone DNA synthesis. Yet, most of the interaction sites necessary for these vital processes are largely unknown. By comparing divergences among RecA sequences and computing putative functional regions, we discovered four functional sites of RecA. Targeted point-mutations were then tested for both recombination and DNA damage induction and reveal distinct RecA functions at each one of these sites. In particular, one new set of mutants is deficient in promoting LexA cleavage and yet maintains the ability to induce the DNA damage response. These results reveal specific amino acid determinants of the RecA–LexA interaction and suggest that LexA binds RecA_i and RecA_i+6 at a composite site on the RecA filament, which could explain the role of the active filament during LexA cleavage.