BAsE-Seq: a method for obtaining long viral haplotypes from short sequence reads

Epidemiologists aim to inform the design of public health interventions with evidence on the evolution, emergence and spread of infectious diseases. Sequencing of pathogen genomes, together with date, location, clinical manifestation and other relevant data about sample origins, can contribute to describing nearly every aspect of transmission dynamics, including local transmission and global spread. The analyses of these data have implications for all levels of clinical and public health practice, from institutional infection control to policies for surveillance, prevention and treatment. This review highlights the range of epidemiological questions that can be addressed from the combination of genome sequence and traditional ‘line lists’ (tables of epidemiological data where each line includes demographic and clinical features of infected individuals). We identify opportunities for these data to inform interventions that reduce disease incidence and prevalence. By considering current limitations of, and challenges to, interpreting these data, we aim to outline a research agenda to accelerate the genomics-driven transformation in public health microbiology.     

Studies on the DIDS-binding Site of Monocarboxylate Transporter 1 Suggest a Homology Model of the Open Conformation and a Plausible Translocation Cycle

Site-directed mutagenesis of MCT1 was performed on exofacial lysines Lys³⁸, Lys⁴⁵, Lys²⁸², and Lys⁴¹³. K38Q-MCT1 and K38R-MCT1 were inactive when expressed at the plasma membrane of Xenopus laevis oocytes, whereas K45R/K282R/K413R-MCT1 and K45Q/K282Q/K413Q-MCT1 were active. The former exhibited normal reversible and irreversible inhibition by DIDS, whereas the latter showed less reversible and no irreversible inhibition. K45Q/K413Q-MCT1 retained some irreversible inhibition, whereas K45Q/K282Q-MCT1 and K282Q/K413Q-MCT1 did not. These data suggest that the two DIDS SO₃⁻ groups interact with positively charged Lys²⁸² together with Lys⁴⁵ and/or Lys⁴¹³. This positions one DIDS isothiocyanate group close to Lys³⁸, leading to its covalent modification and irreversible inhibition. Additional mutagenesis revealed that DIDS cross-links MCT1 to its ancillary protein embigin using either Lys³⁸ or Lys²⁹⁰ of MCT1 and Lys¹⁶⁰ or Lys¹⁶⁴ of embigin. We have modeled a possible structure for the outward facing (open) conformation of MCT1 by employing modest rotations of the C-terminal domain of the inner facing conformation modeled previously. The resulting model structure has a DIDS-binding site consistent with experimental data and locates Lys³⁸ in a hydrophobic environment at the bottom of a substrate-binding channel. Our model suggests a translocation cycle in which Lys³⁸ accepts a proton before binding lactate. Both the lactate and proton are then passed through the channel via Asp³⁰²⁻ and Asp³⁰⁶⁺, an ion pair already identified as important for transport and located adjacent to Phe³⁶⁰, which controls channel selectivity. The cross-linking data have also been used to model a structure of MCT1 bound to embigin that is consistent with published data.Monocarboxylate transporter 1 (MCT1)³ is a member of the monocarboxylate transporter family (SLC16) of which there are 14 known members encoded by both the human and mouse genomes (1). All of the members of this family are thought to have 12 transmembrane alpha helices (TMs) with a large loop between TMs 6 and 7 and the C and N termini facing the cytosol (2, 3). The only members of the MCT family that have been shown to catalyze transport of monocarboxylates such as l-lactate across the plasma membrane are isoforms 1–4 (4–8). This transport is proton-linked and leads to the net uptake or release of lactic acid from cells, which is critical for metabolic pathways such as anaerobic glycolysis, gluconeogenesis, and lactate oxidation (9). MCT8 is a high affinity thyroid hormone transporter (10), whereas MCT10 (TAT1) is an aromatic amino acid transporter (11). The other members of the MCT family remain to be characterized.MCT1 is the most widely distributed member of the MCT family and was first identified as the lactate transporter present in red blood cells where its kinetics and substrate and inhibitor specificity were investigated in detail (9, 11, 12). These studies revealed that MCT1 can be inhibited by stilbene disulfonate derivatives such as DIDS and 4,4′-dibenzamido-stilbene-2,2′-disulfonate (DBDS). DIDS was shown to exhibit a rapid reversible inhibition of transport that was competitive with respect to l-lactate. This is followed by a slowly developing irreversible inhibition that is not exhibited by DBDS and is thought to be caused by one of the isothiocyanate groups of DIDS attacking a lysine residue on MCT1 (13–15). Prolonged incubation with DIDS also led to a fraction of the MCT1 becoming cross-linked to a 70-kDa glycoprotein that was identified as embigin, also known as gp70 (16). Embigin has a short intracellular C terminus, a single TM sequence containing a glutamic acid residue, and a large extracellular N terminus containing two immunoglobulin domains (17, 18). Subsequent studies revealed that either embigin or, more frequently, the homologous protein basigin (also known as CD147) is required as a chaperone to take MCT1 to the membrane (19) where the two proteins must remain associated for transport activity to be maintained (20, 21).Expression of MCTs 1, 2, and 4 in Xenopus laevis oocytes has enabled their further characterization and the effects of site-directed mutagenesis to be investigated (4, 5, 7, 8, 22–24). Such studies, together with homology modeling have enabled us to propose a three-dimensional structure of MCT1 based around the published structure of the Escherichia coli glycerol-3-phosphate transporter (Protein Data Bank 1PW4) (24). This model can account for the effects of mutating a range of amino acids, including some that disrupt the interaction with basigin, and has led to the proposal that the single TM of basigin or embigin lies between TMs 3 and 6 of MCT1. The model also reveals exofacial lysines that are present in MCT1 that might be responsible for the irreversible inhibition of MCT1 by DIDS and the cross-linking of MCT1 to embigin. In rat MCT1 these residues are Lys³⁸, Lys⁴⁵, Lys²⁸², Lys²⁸⁴, Lys²⁹⁰, and Lys⁴¹³. In this paper, we use site-directed mutagenesis of these lysine residues to identify which of them are involved in DIDS binding to MCT1. In addition we use site-directed mutagenesis of embigin to demonstrate that Lys¹⁶⁰ and Lys¹⁶⁴ are involved in its cross-linking to MCT1. Our new data allow us to propose a modified structural model of MCT1 in its outward facing conformation that binds DIDS. This model is consistent with the site-directed mutagenesis data and also suggests a mechanism for the translocation cycle that involves Lys³⁸ as well as Asp³⁰² and Arg³⁰⁶ that have already been identified as important for transport (23, 24). We have also been able to model a structure of MCT1 bound to embigin that is consistent with published data.     

The Ability of GAP1IP4BP To Function as a Rap1 GTPase-Activating Protein (GAP) Requires Its Ras GAP-Related Domain and an Arginine Finger Rather than an Asparagine Thumb

GAP1^IP4BP is a member of the GAP1 family of Ras GTPase-activating proteins (GAPs) that includes GAP1^m, CAPRI, and RASAL. Composed of a central Ras GAP-related domain (RasGRD), surrounded by amino-terminal C2 domains and a carboxy-terminal PH/Btk domain, these proteins, with the notable exception of GAP1^m, possess an unexpected arginine finger-dependent GAP activity on the Ras-related protein Rap1 (S. Kupzig, D. Deaconescu, D. Bouyoucef, S. A. Walker, Q. Liu, C. L. Polte, O. Daumke, T. Ishizaki, P. J. Lockyer, A. Wittinghofer, and P. J. Cullen, J. Biol. Chem. 281:9891-9900, 2006). Here, we have examined the mechanism through which GAP1^IP4BP can function as a Rap1 GAP. We show that deletion of domains on either side of the RasGRD, while not affecting Ras GAP activity, do dramatically perturb Rap1 GAP activity. By utilizing GAP1^IP4BP/GAP1^m chimeras, we establish that although the C2 and PH/Btk domains are required to stabilize the RasGRD, it is this domain which contains the catalytic machinery required for Rap1 GAP activity. Finally, a key residue in Rap1-specific GAPs is a catalytic asparagine, the so-called asparagine thumb. By generating a molecular model describing the predicted Rap1-binding site in the RasGRD of GAP1^IP4BP, we show that mutagenesis of individual asparagine or glutamine residues that lie in close proximity to the predicted binding site has no detectable effect on the in vivo Rap1 GAP activity of GAP1^IP4BP. In contrast, we present evidence consistent with a model in which the RasGRD of GAP1^IP4BP functions to stabilize the switch II region of Rap1, allowing stabilization of the transition state during GTP hydrolysis initiated by the arginine finger.The Ras-like family of small GTPases are ubiquitously expressed, evolutionarily conserved proteins that, by undergoing conformational changes in response to the alternate binding of GDP and GTP, function as binary switches (28, 31, 35). The GDP-bound “off” state and the GTP-bound “on” state recognize distinct effector proteins, thereby allowing the regulation of a variety of downstream signaling events (28, 31, 35). While Ras is the best-known and best-studied Ras-like GTPase, Rap1 has recently attracted considerable attention (reviewed in reference 20).Rap1 was originally identified through its ability, when overexpressed, to reverse the phenotype of K-Ras-transformed NIH 3T3 cells (19). As Ras and Rap1 have very similar effector regions, the ability of Rap1 to reverse the transformed phenotype appeared to arise through an ability to compete with K-Ras effectors. For example, Rap1 binds the Ras effector Raf1 but this does not lead to its activation (11). This is consistent with a simple model in which Rap1 functions as a Ras antagonist (6, 37). However, recent work has challenged this view. Increasing evidence points to Rap1 interacting with its own panel of effectors through which it controls cell-cell adhesion and cell-matrix interactions (reviewed in reference 20).Like that of other GTPases, the activation of Ras and Rap1 is regulated through guanine nucleotide exchange factors, which control activation by stimulating the exchange of GDP for GTP. Inactivation is driven by GTPase-activating proteins (GAPs). These enhance the intrinsic GTPase activity of Ras and Rap1, thereby leading to GTP hydrolysis. A wide variety of guanine nucleotide exchange factors and GAPs specific for these GTPases have been identified (14). Through the arrangement of different modular domains, these proteins are regulated following the activation of cell surface receptors. This occurs either through direct association with the activated receptor or indirectly through second messengers (4, 5, 14, 41).Mammalian proteins capable of functioning as Ras GAPs include NF1 (3, 27, 40), p120^GAP (38), the semaphorin 4D receptor plexin-B1 (29), and members of the GAP1 (reviewed in reference 41) and SynGAP (DAB2IP, nGAP, and SynGAP) families (10, 18, 39). These function as Ras GAPs by supplying a catalytic arginine residue—the arginine finger—into the active site of Ras. This stabilizes the transition state of the GTPase reaction, increasing the reaction rate by more than 1,000-fold (1, 33, 34).Rap1 GAPs include Rap GAPs I and II, the SPA-1 family (SPA-1, SPAR, SPAL, and E6TP1), and tuberin (16, 17, 26, 32). Unlike Ras, Rap1 does not possess the catalytic glutamine residue that is critical for GTP hydrolysis in Ras. This fundamental difference means that the mechanisms by which Ras and Rap1 GAPs function are distinct. Rap1 GAPs do not employ a catalytic arginine residue (8, 9); instead, they provide a catalytic asparagine—the asparagine thumb—to stimulate GTP hydrolysis (15). Here the asparagine carboxamide side chain has a function similar to that of the glutamine residue in Ras, stabilizing the position of the nucleophilic water and γ-phosphate in the transition complex (15, 36).Given such distinct catalytic mechanisms, surprisingly, some Ras GAPs, while having no detectable sequence homology with any Rap1 GAPs, are capable of stimulating the GTPase activity of Rap1. The first protein found to display such dual activities was GAP1^IP4BP (13) (also known as RASA3, GAPIII, and R-Ras GAP). This is a member of the GAP1 family, which also comprises GAP1^m, CAPRI, and RASAL (2, 23-25). These proteins are characterized by a domain architecture comprising amino-terminal tandem C2 domains, a highly conserved central Ras GAP-related domain (RasGRD), and a carboxy-terminal pleckstrin homology (PH) domain that is associated with a Bruton''s tyrosine kinase (Btk) motif (41). Consistent with the presence of the RasGRD, all proteins display Ras GAP activity, although each is differentially regulated following receptor stimulation (41). With the notable exception of GAP1^m, all GAP1 proteins also possess efficient Rap1 GAP activity (22). Such dual specificity is not restricted solely to GAP1 proteins. Recently, C2 domain-containing SynGAP—a neuronal Ras GAP—has also been shown to display Rap1 GAP activity (21), an activity that appears to require, alongside the RasGRD, the presence of a single C2 domain (30).Here we have examined the mechanism behind the dual Ras and Rap1 GAP activities of GAP1^IP4BP. Through the generation of a series of GAP1^IP4BP/GAP1^m chimeras, we have established that while the C2 domains of GAP1^IP4BP are required to stabilize the RasGRD, these domains do not supply catalytic residues required for Rap1 GAP activity. Rather, the Rap1 GAP catalytic machinery appears to reside solely within the RasGRD. By the site-directed mutagenesis of selected asparagine and glutamine residues within this domain—selected following the generation of a predicted molecular model of the GAP1^IP4BP RasGRD-Ras(Rap1) complex—we establish that the ability of GAP1^IP4BP to function as a Rap1 GAP does not occur via a mechanism that utilizes a classic asparagine thumb. Rather, we suggest that the GAP1^IP4BP RasGRD functions to stabilize the switch II region of Rap1 in a manner that allows a catalytic arginine finger from GAP1^IP4BP to drive the hydrolysis of GTP.     

Allometry between leg and body length of insects: lack of support for the size–grain hypothesis

• 1 The size–grain hypothesis ( Kaspari & Weiser, 1999 ) states that (1) as organisms decrease in size, they perceive their environment as being more rugose; (2) long legs allow organisms to step over obstacles but hinder them from entering small gaps; and (3) as the size of an organism decreases, the benefits of long legs begin to be outweighed by the costs of construction. Natural selection should therefore favour proportionally longer legs in larger organisms, thereby leading to a positive allometry between leg and body length (scaling exponent b > 1).
• 2 Here we compare the scaling exponent of leg‐to‐body length relationships among insects that walk, walk and fly, and predominantly fly. We measured the lengths of the hind tibia, hind femur, and body length of each species.
• 3 The taxa varied considerably in the scaling exponent b. In seven out of ten groups (Formicidae, Isoptera, Carabidae, Pentatomidae, Apidae, Lepidoptera, Odonata adult), b was significantly greater than one. However, there was no gradual decrease in b from walking to walking/flying to flying insects.
• 4 The results of the present study provide no support for the size–grain hypothesis. We propose that leg length is not only affected by the rugosity of the environment, but also by (1) functional adaptations, (2) phylogeny, (3) lifestyle, (4) the type of insect development (hemimetabolism or holometabolism), and (5) constraints of gas exchange.

     

Crystal Structure of the LasA Virulence Factor from Pseudomonas aeruginosa: Substrate Specificity and Mechanism of M23 Metallopeptidases

Pseudomonas aeruginosa is an opportunist Gram-negative bacterial pathogen responsible for a wide range of infections in immunocompromized individuals and is a leading cause of mortality in cystic fibrosis patients. A number of secreted virulence factors, including various proteolytic enzymes, contribute to the establishment and maintenance of Pseudomonas infection. One such is LasA, an M23 metallopeptidase related to autolytic glycylglycine endopeptidases such as Staphylococcus aureus lysostaphin and LytM, and to DD-endopeptidases involved in entry of bacteriophage to host bacteria. LasA is implicated in a range of processes related to Pseudomonas virulence, including stimulating ectodomain shedding of the cell surface heparan sulphate proteoglycan syndecan-1 and elastin degradation in connective tissue. Here we present crystal structures of active LasA as a complex with tartrate and in the uncomplexed form. While the overall fold resembles that of the other M23 family members, the LasA active site is less constricted and utilizes a different set of metal ligands. The active site of uncomplexed LasA contains a five-coordinate zinc ion with trigonal bipyramidal geometry and two metal-bound water molecules. Using these structures as a starting point, we propose a model for substrate binding by LasA that explains its activity against a wider range of substrates than those used by related lytic enzymes, and offer a catalytic mechanism for M23 metallopeptidases consistent with available structural and mutagenesis data. Our results highlight how LasA is a structurally distinct member of this endopeptidase family, consistent with its activity against a wider range of substrates and with its multiple roles in Pseudomonas virulence.     

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).     

A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors

Changeux et al. (Changeux et al. C. R. Biol. 343:33–39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4β2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.     

PARAMETERS USED FOR RECOGNITION OF DISTRESS CALLS IN TWO SPECIES: LARUS ARGENTATUS AND STURNUS VULGARIS

ABSTRACT

Distress calls are signals effective over a long distance. They are well known to evoke interspecific reactions. We suggest that the interspecificity phenomenon results from the use of similar laws of decoding by the species concerned. These laws must take into account the transmission channel which always has a great influence on long-range communication. We tested our hypothesis by broadcasting simplified synthetic calls to two species of birds: the herring gull and the starling. The various calls differed in terms of frequency modulation (FM). Two main conclusions emerged from this series of tests:

1. The parameters used for recognition are not sophisticated: a simple slope applied to a carrier frequency that corresponds to the acoustic shape of a distress call is sufficient to confer a distress meaning to the signal. The basic rules are the same for the gull and the starling, with differences only in the acceptance level of the species.

2. The system of recognition is based upon parameters not altered by the environment: the birds make use of the slow frequency modulations (FMI). In contrast, the fast frequency sweeps (FMII) which are modified during propagation do not seem to be utilised. The use of these characteristics of distress calls for recognition allows interspecificity and maximum efficiency for propagation over long distances.

Antibiotic resistance and pathogenicity factors in <Emphasis Type="Italic">Staphylococcus aureus</Emphasis> isolated from mastitic Sahiwal cattle

Methicillin-resistant Staphylococcus aureus (MRSA) poses a serious problem in dairy animals suffering from mastitis. In the present study, the distribution of mastitic MRSA and antibiotic resistance was studied in 107 strains of S. aureus isolated from milk samples from 195 infected udders. The characterizations pathogenic factors (adhesin and toxin genes) and antibiotic susceptibility of isolates were carried out using gene amplification and disc diffusion assays, respectively. A high prevalence of MRSA was observed in the tested isolates (13.1%). The isolates were also highly resistant to antibiotics, i.e. 36.4% were resistant to streptomycin, 33.6% to oxytetracycline, 29.9% to gentamicin and 26.2% each to chloramphenicol, pristinomycin and ciprofloxacin. A significant variation in the expression of pathogenic factors (Ig, coa and clf) was observed in these isolates. The overall distribution of adhesin genes ebp, fib, bbp, fnbB, cap5, cap8, map and cna in the isolates was found to be 69.1, 67.2, 6.5, 20.5, 60.7, 26.1, 81.3 and 8.4%, respectively. The presence of fib, fnbB, bbp and map genes was considerably greater in MRSA than in methicillin-susceptible S. aureus (MSSA) isolates. The proportions of toxin genes, namely, hlb, seb, sec, sed, seg and sei, in the isolates were found to be 94.3, 0.9, 8.4, 0.9, 10.2 and 49.5%, respectively. The proportions of agr genes I, II, III and IV were found to be 39.2, 27.1, 21.5 and 12.1%, respectively. A few isolates showed similar antibiotic-resistance patterns, which could be due to identical strains or the dissemination of the same strains among animals. These findings can be utilized in mastitis treatment programmes and antimicrobials strategies in organized herds.