Mixed, nonclassical behavior in a classic allosteric protein.

Allostery is a major driver of biological processes requiring coordination. Thus, it is one of the most fundamental and remarkable phenomena in nature, and there is motivation to understand and manipulate it to a multitude of ends. Today, it is often described in terms of two phenomenological models proposed more than a half-century ago involving only T(tense) or R(relaxed) conformations. Here, methyl-based NMR provides extensive detail on a dynamic T to R switch in the classical dimeric allosteric protein, yeast chorismate mutase (CM), that occurs in the absence of substrate, but only with the activator bound. Switching of individual subunits is uncoupled based on direct observation of mixed TR states in the dimer. This unique finding excludes both classic models and solves the paradox of a coexisting hyperbolic binding curve and highly skewed substrate-free T-R equilibrium. Surprisingly, structures of the activator-bound and effector-free forms of CM appear the same by NMR, providing another example of the need to account for dynamic ensembles. The apo enzyme, which has a sigmoidal activity profile, is shown to switch, not to R, but to a related high-energy state. Thus, the conformational repertoire of CM does not just change as a matter of degree depending on the allosteric input, be it effector and/or substrate. Rather, the allosteric model appears to completely change in different contexts, which is only consistent with modern ensemble-based frameworks.

BigBind: Learning from Nonstructural Data for Structure-Based Virtual Screening.

Deep learning methods that predict protein-ligand binding have recently been used for structure-based virtual screening. Many such models have been trained using protein-ligand complexes with known crystal structures and activities from the PDBBind data set. However, because PDBbind only includes 20K complexes, models typically fail to generalize to new targets, and model performance is on par with models trained with only ligand information. Conversely, the ChEMBL database contains a wealth of chemical activity information but includes no information about binding poses. We introduce BigBind, a data set that maps ChEMBL activity data to proteins from the CrossDocked data set. BigBind comprises 583 K ligand activities and includes 3D structures of the protein binding pockets. Additionally, we augmented the data by adding an equal number of putative inactives for each target. Using this data, we developed Banana (basic neural network for binding affinity), a neural network-based model to classify active from inactive compounds, defined by a 10 µM cutoff. Our model achieved an AUC of 0.72 on BigBind's test set, while a ligand-only model achieved an AUC of 0.59. Furthermore, Banana achieved competitive performance on the LIT-PCBA benchmark (median EF1% 1.81) while running 16,000 times faster than molecular docking with Gnina. We suggest that Banana, as well as other models trained on this data set, will significantly improve the outcomes of prospective virtual screening tasks.

STOPLIGHT: A Hit Scoring Calculator.

We introduce STOPLIGHT, a web portal to assist medicinal chemists in prioritizing hits from screening campaigns and the selection of compounds for optimization. STOPLIGHT incorporates services to assess 6 physiochemical and structural properties, 6 assay liabilities, and 11 pharmacokinetic properties, for any small molecule represented by its SMILES string. We briefly describe each service and illustrate the utility of this portal with a case study. The STOPLIGHT portal provides a user-friendly tool to guide hit selection in early drug discovery campaigns, whereby compounds with unfavorable properties can be quickly recognized and eliminated.

Accurate ligand-protein docking in CASP15 using the ClusPro LigTBM server.

In the ligand prediction category of CASP15, the challenge was to predict the positions and conformations of small molecules binding to proteins that were provided as amino acid sequences or as models generated by the AlphaFold2 program. For most targets, we used our template-based ligand docking program ClusPro ligTBM, also implemented as a public server available at https://ligtbm.cluspro.org/. Since many targets had multiple chains and a number of ligands, several templates, and some manual interventions were required. In a few cases, no templates were found, and we had to use direct docking using the Glide program. Nevertheless, ligTBM was shown to be a very useful tool, and by any ranking criteria, our group was ranked among the top five best-performing teams. In fact, all the best groups used template-based docking methods. Thus, it appears that the AlphaFold2-generated models, despite the high accuracy of the predicted backbone, have local differences from the x-ray structure that make the use of direct docking methods more challenging. The results of CASP15 confirm that this limitation can be frequently overcome by homology-based docking.

Tyrosinase-Catalyzed Peptide Macrocyclization for mRNA Display.

mRNA display of macrocyclic peptides has proven itself to be a powerful technique to discover high-affinity ligands for a protein target. However, only a limited number of cyclization chemistries are known to be compatible with mRNA display. Tyrosinase is a copper-dependent oxidase that oxidizes tyrosine phenol to an electrophilic o-quinone, which is readily attacked by cysteine thiol. Here we show that peptides containing tyrosine and cysteine are rapidly cyclized upon tyrosinase treatment. Characterization of the cyclization reveals it to be widely applicable to multiple macrocycle sizes and scaffolds. We combine tyrosinase-mediated cyclization with mRNA display to discover new macrocyclic ligands targeting melanoma-associated antigen A4 (MAGE-A4). These macrocycles potently inhibit the MAGE-A4 binding axis with nanomolar IC50 values. Importantly, macrocyclic ligands show clear advantage over noncyclized analogues with ∼40-fold or greater decrease in IC50 values.

Fluorescent-Tagged Antiscalants-The New Materials for Scale Inhibition Mechanism Studies, Antiscalant Traceability and Antiscaling Efficacy Optimization during CaCO3 and CaSO4·2H2O Scale Formation.

Equipment scaling leads to reduced production efficiency in a wide range of industrial applications worldwide. Various antiscaling agents are currently commonly used to mitigate this problem. However, irrespective of their long and successful application in water treatment technologies, little is known about the mechanisms of scale inhibition, particularly the localization of scale inhibitors on scale deposits. The lack of such knowledge is a limiting factor in the development of applications for antiscalants. Meanwhile, fluorescent fragments integrated into scale inhibitor molecules have provided a successful solution to the problem. The focus of this study is, therefore, on the synthesis and investigation of a novel fluorescent antiscalant: (2-(6-morpholino-1,3-dioxo-1H-benzo[de]isoquinolin-2(3H)yl)ethylazanediyl)bis(methylenephosphonic acid) (ADMP-F) which is an analog of the commercial antiscalant: aminotris(methylenephosphonic acid) (ATMP). ADMP-F has been found to effectively control the precipitation of CaCO3 and CaSO4 in solution and is a promising tracer for organophosphonate scale inhibitors. ADMP-F was compared with two other fluorescent antiscalants-polyacrylate (PAA-F1) and bisphosphonate (HEDP-F)-and was found to be highly effective: PAA-F1 > ADMP-F >> HEDP-F (CaCO3) and PAA-F1 > ADMP-F > HEDP-F (CaSO4·2H2O). The visualization of the antiscalants on the deposits provides unique information on their location and reveals differences in the "antiscalant-deposit" interactions for scale inhibitors of different natures. For these reasons, a number of important refinements to the mechanisms of scale inhibition are proposed.

Mechanisms of Phototoxic Effects of Cationic Porphyrins on Human Cells In Vitro.

The toxic effects of four cationic porphyrins on various human cells were studied in vitro. It was found that, under dark conditions, porphyrins are almost nontoxic, while, under the action of light, the toxic effect was observed starting from nanomolar concentrations. At a concentration of 100 nM, porphyrins caused inhibition of metabolism in the MTT test in normal and cancer cells. Furthermore, low concentrations of porphyrins inhibited colony formation. The toxic effect was nonlinear; with increasing concentrations of various porphyrins, up to about 1 µM, the effect reached a plateau. In addition to the MTT test, this was repeated in experiments examining cell permeability to trypan blue, as well as survival after 24 h. The first visible manifestation of the toxic action of porphyrins is blebbing and swelling of cells. Against the background of this process, permeability to porphyrins and trypan blue appears. Subsequently, most cells (even mitotic cells) freeze in this swollen state for a long time (24 and even 48 h), remaining attached. Cellular morphology is mostly preserved. Thus, it is clear that the cells undergo mainly necrotic death. The hypothesis proposed is that the concentration dependence of membrane damage indicates a limited number of porphyrin targets on the membrane. These targets may be any ion channels, which should be considered in photodynamic therapy.

Mapping co-regulatory interactions among ligand-binding sites in ryanodine receptor 1.

Ryanodine receptor 1 (RyR1) is an intracellular calcium ion (Ca2+ ) release channel required for skeletal muscle contraction. Although cryo-electron microscopy identified binding sites of three coactivators Ca2+ , ATP, and caffeine (CFF), the mechanism of co-regulation and synergy of these activators is unknown. Here, we report allosteric connections among the three ligand-binding sites and pore region in (i) Ca2+ bound-closed, (ii) ATP/CFF bound-closed, (iii) Ca2+ /ATP/CFF bound-closed, and (iv) Ca2+ /ATP/CFF bound-open RyR1 states. We identified two dominant networks of interactions that mediate communication between the Ca2+ -binding site and pore region in Ca2+ bound-closed state, which partially overlapped with the pore communications in ATP/CFF bound-closed RyR1 state. In Ca2+ /ATP/CFF bound-closed and -open RyR1 states, co-regulatory interactions were analogous to communications in the Ca2+ bound-closed and ATP/CFF bound-closed states. Both ATP- and CFF-binding sites mediate communication between the Ca2+ -binding site and the pore region in Ca2+ /ATP/CFF bound-open RyR1 structure. We conclude that Ca2+ , ATP, and CFF propagate their effects to the pore region through a network of overlapping interactions that mediate allosteric control and molecular synergy in channel regulation.

Harmonically mode-locked laser pulse accumulation in a self-resonating optical cavity.

Optical enhancement cavities enabling laser pulses to be coherently stacked in free space are used in several applications to enhance accessible optical power. In this study, we develop an optical cavity that accumulates harmonically mode-locked laser pulses with a self-resonating mechanism for X-ray sources based on laser-Compton scattering. In particular, a Fabry-Perot cavity composed of 99% reflectance mirrors maintained the optical resonance in a feedback-free fashion for more than half an hour and automatically resumed the accumulation even if the laser oscillation was suspended. In contrast to conventional optical enhancement cavity systems with a dedicated feedback controller, this characteristic is highly beneficial in practical applications, such as for laser-Compton scattering X-ray sources. Lastly, upscaling and adoption of the proposed system might improve the operability and equipment use of laser Compton-scattering X-ray sources.

Integrated approach to elucidate metal-implant related adverse outcome pathways.

Exogenous metal particles and ions from implant devices are known to cause severe toxic events with symptoms ranging from adverse local tissue reactions to systemic toxicities, potentially leading to the development of cancers, heart conditions, and neurological disorders. Toxicity mechanisms, also known as Adverse Outcome Pathways (AOPs), that explain these metal-induced toxicities are severely understudied. Therefore, we deployed in silico structure- and knowledge-based approaches to identify proteome-level perturbations caused by metals and pathways that link these events to human diseases. We captured 177 structure-based, 347 knowledge-based, and 402 imputed metal-gene/protein relationships for chromium, cobalt, molybdenum, nickel, and titanium. We prioritized 72 proteins hypothesized to directly contact implant surfaces and contribute to adverse outcomes. Results of this exploratory analysis were formalized as structured AOPs. We considered three case studies reflecting the following possible situations: (i) the metal-protein-disease relationship was previously known; (ii) the metal-protein, protein-disease, and metal-disease relationships were individually known but were not linked (as a unified AOP); and (iii) one of three relationships was unknown and was imputed by our methods. These situations were illustrated by case studies on nickel-induced allergy/hypersensitivity, cobalt-induced heart failure, and titanium-induced periprosthetic osteolysis, respectively. All workflows, data, and results are freely available in https://github.com/DnlRKorn/Knowledge_Based_Metallomics/. An interactive view of select data is available at the ROBOKOP Neo4j Browser at http://robokopkg.renci.org/browser/.

Structure of prion ß-oligomers as determined by short-distance crosslinking constraint-guided discrete molecular dynamics simulations.

The conversion of the native monomeric cellular prion protein (PrPC ) into an aggregated pathological ß-oligomeric form (PrPß ) and an infectious form (PrPSc ) is the central element in the development of prion diseases. The structure of the aggregates and the molecular mechanisms of the conformational changes involved in the conversion are still unknown. We applied mass spectrometry combined with chemical crosslinking, hydrogen/deuterium exchange, limited proteolysis, and surface modification for the differential characterization of the native and the urea+acid-converted prion ß-oligomer structures to obtain insights into the mechanisms of conversion and aggregation. For the determination of the structure of the monomer and the dimer unit of the ß-oligomer, we applied a recently-developed approach for de novo protein structure determination which is based on the incorporation of zero-length and short-distance crosslinking data as intra- and inter-protein constraints in discrete molecular dynamics simulations (CL-DMD). Based on all of the structural-proteomics experimental data and the computationally predicted structures of the monomer units, we propose the potential mode of assembly of the ß-oligomer. The proposed ß-oligomer assembly provides a clue on the ß-sheet nucleation site, and how template-based conversion of the native prion molecule occurs, growth of the prion aggregates, and maturation into fibrils may occur.

Conformational ensemble of native α-synuclein in solution as determined by short-distance crosslinking constraint-guided discrete molecular dynamics simulations.

Combining structural proteomics experimental data with computational methods is a powerful tool for protein structure prediction. Here, we apply a recently-developed approach for de novo protein structure determination based on the incorporation of short-distance crosslinking data as constraints in discrete molecular dynamics simulations (CL-DMD) for the determination of conformational ensemble of the intrinsically disordered protein α-synuclein in the solution. The predicted structures were in agreement with hydrogen-deuterium exchange, circular dichroism, surface modification, and long-distance crosslinking data. We found that α-synuclein is present in solution as an ensemble of rather compact globular conformations with distinct topology and inter-residue contacts, which is well-represented by movements of the large loops and formation of few transient secondary structure elements. Non-amyloid component and C-terminal regions were consistently found to contain ß-structure elements and hairpins.

Inter-Active Site Communication Mediated by the Dimer Interface ß-Sheet in the Half-the-Sites Enzyme, Thymidylate Synthase.

Thymidylate synthase (TS) is a dimeric enzyme conserved in all life forms that exhibits the allosteric feature of half-the-sites activity. Neither the reason for nor the mechanism of this phenomenon is understood. We used a combined nuclear magnetic resonance (NMR) and molecular dynamics approach to study a stable intermediate preceding hydride transfer, which is the rate-limiting and half-the-sites step. In NMR titrations with ligands leading to this intermediate, we measured chemical shifts of the apoenzyme (lig0), the saturated holoenzyme (lig2), and the typically elusive singly bound (lig1) states. Approximately 40 amides showed quartet patterns providing direct NMR evidence of coupling between the active site and probes >30 Å away in the distal subunit. Quartet peak patterns have symmetrical character, indicating reciprocity in communicating the first and second binding events to the distal protomer. Quartets include key catalytic residues and map to the dimer interface ß-sheet, which also represents the shortest path between the two active sites. Simulations corroborate the coupling observed in solution in that there is excellent overlap between quartet residues and main-chain atoms having intersubunit cross-correlated motions. Simulations identify five hot spot residues, three of which lie at the kink in the unique ß-bulge abutting the active sites on either end of the sheet. Interstrand cross-correlated motions become more organized and pronounced as the enzyme progresses from lig0 to lig1 and ultimately lig2. Coupling in the apparently symmetrical complex has implications for half-the-sites reactivity and potentially resolves the paradox of inequivalent TS active sites despite the vast majority of X-ray structures appearing to be symmetrical.

Mapping allosteric linkage to channel gating by extracellular domains in the human epithelial sodium channel.

The epithelial sodium channel (ENaC) mediates sodium absorption in lung, kidney, and colon epithelia. Channels in the ENaC/degenerin family possess an extracellular region that senses physicochemical changes in the extracellular milieu and allosterically regulates the channel opening. Proteolytic cleavage activates the ENaC opening, by the removal of specific segments in the finger domains of the α- and γ ENaC-subunits. Cleavage causes perturbations in the extracellular region that propagate to the channel gate. However, it is not known how the channel structure mediates the propagation of activation signals through the extracellular sensing domains. Here, to identify the structure-function determinants that mediate allosteric ENaC activation, we performed MD simulations, thiol modification of residues substituted by cysteine, and voltage-clamp electrophysiology recordings. Our simulations of an ENaC heterotetramer, α1ßα2γ, in the proteolytically cleaved and uncleaved states revealed structural pathways in the α-subunit that are responsible for ENaC proteolytic activation. To validate these findings, we performed site-directed mutagenesis to introduce cysteine substitutions in the extracellular domains of the α-, ß-, and γ ENaC-subunits. Insertion of a cysteine at the α-subunit Glu557 site, predicted to stabilize a closed state of ENaC, inhibited ENaC basal activity and retarded the kinetics of proteolytic activation by 2-fold. Our results suggest that the lower palm domain of αENaC is essential for ENaC activation. In conclusion, our integrated computational and experimental approach suggests key structure-function determinants for ENaC proteolytic activation and points toward a mechanistic model for the allosteric communication in the extracellular domains of the ENaC/degenerin family channels.

High-speed atomic force microscopy reveals structural dynamics of α-synuclein monomers and dimers.

α-Synuclein (α-syn) is the major component of the intraneuronal inclusions called Lewy bodies, which are the pathological hallmark of Parkinson's disease. α-Syn is capable of self-assembly into many different species, such as soluble oligomers and fibrils. Even though attempts to resolve the structures of the protein have been made, detailed understanding about the structures and their relationship with the different aggregation steps is lacking, which is of interest to provide insights into the pathogenic mechanism of Parkinson's disease. Here we report the structural flexibility of α-syn monomers and dimers in an aqueous solution environment as probed by single-molecule time-lapse high-speed AFM. In addition, we present the molecular basis for the structural transitions using discrete molecular dynamics (DMD) simulations. α-Syn monomers assume a globular conformation, which is capable of forming tail-like protrusions over dozens of seconds. Importantly, a globular monomer can adopt fully extended conformations. Dimers, on the other hand, are less dynamic and show a dumbbell conformation that experiences morphological changes over time. DMD simulations revealed that the α-syn monomer consists of several tightly packed small helices. The tail-like protrusions are also helical with a small ß-sheet, acting as a "hinge". Monomers within dimers have a large interfacial interaction area and are stabilized by interactions in the non-amyloid central (NAC) regions. Furthermore, the dimer NAC-region of each α-syn monomer forms a ß-rich segment. Moreover, NAC-regions are located in the hydrophobic core of the dimer.

Migratory birds along the Mediterranean - Black Sea Flyway as carriers of zoonotic pathogens.

At the crossroad between Europe, Asia, and Africa, Bulgaria is part of the Mediterranean - Black Sea Flyway (MBSF) used by millions of migratory birds. In this study, bird species migrating through Bulgaria were investigated as carriers of zoonotic pathogens. In total, 706 birds belonging to 46 species were checked for the presence of various bacterial pathogens (Campylobacter, Yersinia, Salmonella, Listeria, Escherichia coli, Staphylococcus aureus, Francisella tularensis, Coxiella burnetii, Borrelia burgdorferi, and Brucella spp.). From 673 birds we investigated fecal samples, from the remaining 33, blood samples. We detected Campylobacter 16S rDNA gene in 1.3% of birds, but none were of pathogenic Campylobacter jejuni and Campylobacter coli species. Escherichia coli 16S rDNA gene was found in 8.8% of the birds. Out of 34 birds that transported Yersinia enterocolitica strains (5.05%), only 1 carried a pathogenic isolate. Three birds (0.4%) were carriers of nonpathogenic Salmonella strains. Four avian samples (0.6%) were positive for Listeria monocytogenes and 1 (0.15%) was positive for Brucella spp. None of the birds tested carried the tick-borne pathogens C. burnetii or B. burgdorferi sensu lato. Antibiotic-resistant strains were detected, suggesting that migratory birds could be reservoirs and spreaders of bacterial pathogens as well as antibiotic resistance genes.

MEDYAN: Mechanochemical Simulations of Contraction and Polarity Alignment in Actomyosin Networks.

Active matter systems, and in particular the cell cytoskeleton, exhibit complex mechanochemical dynamics that are still not well understood. While prior computational models of cytoskeletal dynamics have lead to many conceptual insights, an important niche still needs to be filled with a high-resolution structural modeling framework, which includes a minimally-complete set of cytoskeletal chemistries, stochastically treats reaction and diffusion processes in three spatial dimensions, accurately and efficiently describes mechanical deformations of the filamentous network under stresses generated by molecular motors, and deeply couples mechanics and chemistry at high spatial resolution. To address this need, we propose a novel reactive coarse-grained force field, as well as a publicly available software package, named the Mechanochemical Dynamics of Active Networks (MEDYAN), for simulating active network evolution and dynamics (available at www.medyan.org). This model can be used to study the non-linear, far from equilibrium processes in active matter systems, in particular, comprised of interacting semi-flexible polymers embedded in a solution with complex reaction-diffusion processes. In this work, we applied MEDYAN to investigate a contractile actomyosin network consisting of actin filaments, alpha-actinin cross-linking proteins, and non-muscle myosin IIA mini-filaments. We found that these systems undergo a switch-like transition in simulations from a random network to ordered, bundled structures when cross-linker concentration is increased above a threshold value, inducing contraction driven by myosin II mini-filaments. Our simulations also show how myosin II mini-filaments, in tandem with cross-linkers, can produce a range of actin filament polarity distributions and alignment, which is crucially dependent on the rate of actin filament turnover and the actin filament's resulting super-diffusive behavior in the actomyosin-cross-linker system. We discuss the biological implications of these findings for the arc formation in lamellipodium-to-lamellum architectural remodeling. Lastly, our simulations produce force-dependent accumulation of myosin II, which is thought to be responsible for their mechanosensation ability, also spontaneously generating myosin II concentration gradients in the solution phase of the simulation volume.

Docking and Scoring with Target-Specific Pose Classifier Succeeds in Native-Like Pose Identification But Not Binding Affinity Prediction in the CSAR 2014 Benchmark Exercise.

The CSAR 2014 exercise provided an important benchmark for testing current approaches for pose identification and ligand ranking using three X-ray characterized proteins: Factor Xa (FXa), Spleen Tyrosine Kinase (SYK), and tRNA Methyltransferase (TRMD). In Phase 1 of the exercise, we employed Glide and MedusaDock docking software, both individually and in combination, with the special target-specific pose classifier trained to discriminate native-like from decoy poses. All approaches succeeded in the accurate detection of native and native-like poses. We then used Glide SP and MedusaScore scoring functions individually and in combination with the pose-scoring approach to predict relative binding affinities of the congeneric series of ligands in Phase 2 of the exercise. Similar to other participants in the CSAR 2014 exercise, we found that our models showed modest prediction accuracy. Quantitative structure-activity relationship (QSAR) models developed for the FXa ligands using available bioactivity data from ChEMBL showed relatively low prediction accuracy for the CSAR 2014 ligands of the same target. Interestingly, QSAR models built with CSAR data only yielded Spearman correlation coefficients as high as ρ = 0.69 for FXa and ρ = 0.79 for SYK based on 5-fold cross-validation. Virtual screening of the DUD library using the FXa structure was successful in discriminating between active compounds and decoys in spite of poor ranking accuracy of the underlying scoring functions. Our results suggest that two of the three common tasks associated with molecular docking, i.e., native-like pose identification and virtual screening, but not binding affinity prediction, could be accomplished successfully for the CSAR 2014 challenge data set.

The Impact of Collagen Fibril Polarity on Second Harmonic Generation Microscopy.

In this work, we report the implementation of interferometric second harmonic generation (SHG) microscopy with femtosecond pulses. As a proof of concept, we imaged the phase distribution of SHG signal from the complex collagen architecture of juvenile equine growth cartilage. The results are analyzed in respect to numerical simulations to extract the relative orientation of collagen fibrils within the tissue. Our results reveal large domains of constant phase together with regions of quasi-random phase, which are correlated to respectively high- and low-intensity regions in the standard SHG images. A comparison with polarization-resolved SHG highlights the crucial role of relative fibril polarity in determining the SHG signal intensity. Indeed, it appears that even a well-organized noncentrosymmetric structure emits low SHG signal intensity if it has no predominant local polarity. This work illustrates how the complex architecture of noncentrosymmetric scatterers at the nanoscale governs the coherent building of SHG signal within the focal volume and is a key advance toward a complete understanding of the structural origin of SHG signals from tissues.

Extra-slow-growing Tardiphaga strains isolated from nodules of Vavilovia formosa (Stev.) Fed.

Eleven extra-slow-growing strains were isolated from nodules of the relict legume Vavilovia formosa growing in North Ossetia (Caucasus) and Armenia. All isolates formed a single rrs cluster together with the type strain Tardiphaga robiniae LMG 26467(T), while the sequencing of the 16S-23S rDNA intergenic region (ITS) and housekeeping genes glnII, atpD, dnaK, gyrB, recA and rpoB divided them into three groups. North Ossetian isolates (in contrast to the Armenian ones) were clustered separately from the type strain LMG 26467(T). However, all isolates were classified as T. robiniae because the DNA-DNA relatedness between them and the type strain LMG 26467(T) was 69.6% minimum. Two symbiosis-related genes (nodM and nodT) were amplified in all isolated Tardiphaga strains. It was shown that the nodM gene phylogeny is similar to that of ITS and housekeeping genes. The presence of the other symbiosis-related genes in described Tardiphaga strains, which is recently described genus of rhizobia, as well as their ability to form nodules on any plants are under investigation.