Comprehensive Whole-Genome Sequencing and Reporting of Drug Resistance Profiles on Clinical Cases of Mycobacterium tuberculosis in New York State.

Whole-genome sequencing (WGS) is a newer alternative for tuberculosis (TB) diagnostics and is capable of providing rapid drug resistance profiles while performing species identification and capturing the data necessary for genotyping. Our laboratory developed and validated a comprehensive and sensitive WGS assay to characterize Mycobacterium tuberculosis and other M. tuberculosis complex (MTBC) strains, composed of a novel DNA extraction, optimized library preparation, paired-end WGS, and an in-house-developed bioinformatics pipeline. This new assay was assessed using 608 MTBC isolates, with 146 isolates during the validation portion of this study and 462 samples received prospectively. In February 2016, this assay was implemented to test all clinical cases of MTBC in New York State, including isolates and early positive Bactec mycobacterial growth indicator tube (MGIT) 960 cultures from primary specimens. Since the inception of the assay, we have assessed the accuracy of identification of MTBC strains to the species level, concordance with culture-based drug susceptibility testing (DST), and turnaround time. Species identification by WGS was determined to be 99% accurate. Concordance between drug resistance profiles generated by WGS and culture-based DST methods was 96% for eight drugs, with an average resistance-predictive value of 93% and susceptible-predictive value of 96%. This single comprehensive WGS assay has replaced seven molecular assays and has resulted in resistance profiles being reported to physicians an average of 9 days sooner than with culture-based DST for first-line drugs and 32 days sooner for second-line drugs.

Clinical Sensitivity of Cystic Fibrosis Mutation Panels in a Diverse Population.

Infants are screened for cystic fibrosis (CF) in New York State (NYS) using an IRT-DNA algorithm. The purpose of this study was to validate and assess clinical validity of the US FDA-cleared Illumina MiSeqDx CF 139-Variant Assay (139-VA) in the diverse NYS CF population. The study included 439 infants with CF identified via newborn screening (NBS) from 2002 to 2012. All had been screened using the Abbott Molecular CF Genotyping Assay or the Hologic InPlex CF Molecular Test. All with CF and zero or one mutation were tested using the 139-VA. DNA extracted from dried blood spots was reliably and accurately genotyped using the 139-VA. Sixty-three additional mutations were identified. Clinical sensitivity of three panels ranged from 76.2% (23 mutations recommended for screening by ACMG/ACOG) to 79.7% (current NYS 39-mutation InPlex panel), up to 86.0% for the 139-VA. For all, sensitivity was highest in Whites and lowest in the Black population. Although the sample size was small, there was a nearly 20% increase in sensitivity for the Black CF population using the 139-VA (68.2%) over the ACMG/ACOG and InPlex panels (both 50.0%). Overall, the 139-VA is more sensitive than other commercially available panels, and could be considered for NBS, clinical, or research laboratories conducting CF screening.

Restricted domain mobility in the Candida albicans Ess1 prolyl isomerase.

Ess1 is a peptidyl prolyl cis/trans isomerase that is required for virulence of the pathogenic fungi Candida albicans and Cryptococcus neoformans. The enzyme isomerizes the phospho-Ser-Pro linkages in the C-terminal domain of RNA polymerase II. Its human homolog, Pin1, has been implicated in a wide range of human diseases, including cancer and Alzheimer's disease. Crystallographic and NMR studies have demonstrated that the sequence linking the catalytic isomerase domain and the substrate binding WW domain of Pin1 is unstructured and that the two domains are only loosely associated in the absence of the substrate. In contrast, the crystal structure of C. albicans Ess1 revealed a highly ordered linker that contains a three turn alpha-helix and extensive association between the two tightly juxtaposed domains. In part to address the concern that the marked differences in the domain interactions for the human and fungal structures might reflect crystal lattice effects, NMR chemical shift analysis and 15N relaxation measurements have been employed to confirm that the linker of the fungal protein is highly ordered in solution. With the exception of two loops within the active site of the isomerase domain, the local backbone geometry observed in the crystal structure appears to be well preserved throughout the protein chain. The marked differences in interdomain interactions and linker flexibility between the human and fungal enzymes provide a structural basis for therapeutic targeting of the fungal enzymes.

Intron-encoded homing endonuclease I-TevI also functions as a transcriptional autorepressor.

Customary binding sites of intron-encoded homing endonucleases lie within cognate intronless alleles, at the so-called homing sites. Here, we describe a novel, high-affinity binding site for I-TevI endonuclease, encoded within the group I td intron of phage T4. This site is an operator that overlaps the T4 late promoter, which drives I-TevI expression from within the td intron. I-TevI binds the operator and homing sites with equal affinity, and functions as a transcriptional autorepressor. Distinct sequence and spacing requirements of the catalytic domain result in reduced cleavage activity on operator DNA. Crystallographic studies showed that the overall interactions of the DNA-binding domain with the operator and homing sites are similar, but have some different hydrogen-bonding contacts. We present a model in which the flexibility in protein-DNA interactions allows I-TevI to bind variant intronless alleles to promote intron mobility while facilitating its function in autorepression, and thereby persistence in its host.

Crystallographic and mutational studies of Mycobacterium tuberculosis recA mini-inteins suggest a pivotal role for a highly conserved aspartate residue.

The 440 amino acid Mtu recA intein consists of independent protein-splicing and endonuclease domains. Previously, removal of the central endonuclease domain of the intein, and selection for function, generated a 168 residue mini-intein, DeltaI-SM, that had splicing activity similar to that of the full-length, wild-type protein. A D422G mutation (DeltaI-CM) increased C-terminal cleavage activity. Using the DeltaI-SM mini-intein structure (presented here) as a guide, we previously generated a highly active 139 residue mini-intein, DeltaDeltaI(hh)-SM, by replacing 36 amino acid residues in the residual endonuclease loop with a seven-residue beta-turn from the autoprocessing domain of Hedgehog protein. The three-dimensional structures of DeltaI-SM, DeltaDeltaI(hh)-SM, and two variants, DeltaDeltaI(hh)-CM and DeltaDeltaI(hh), have been determined to evaluate the effects of the minimization on intein integrity and to investigate the structural and functional consequences of the D422G mutation. These structural studies show that Asp422 is capable of interacting with both the N and C termini. These interactions are lacking in the CM variant, but are replaced by contacts with water molecules. Accordingly, additional mutagenesis of residue 422, combined with mutations that isolate N-terminal and C-terminal cleavage, showed that the side-chain of Asp422 plays a role in both N and C-terminal cleavage, thereby suggesting that this highly conserved residue regulates the balance between the two reactions.

Minimization and stabilization of the Mycobacterium tuberculosis recA intein.

Many naturally occurring inteins consist of two functionally independent domains, a protein-splicing domain and an endonuclease domain. In a previous study, a 168 amino acid residue mini-intein was generated by removal of the central endonuclease domain of the 440 residue Mycobacterium tuberculosis (Mtu) recA intein. In addition, directed evolution experiments identified a mutation, V67L, that improved the activity of the mini-intein significantly. A recent crystal structure shows that the loop connecting two beta-strands from the N-terminal and C-terminal intein subdomains of the mini-intein is disordered. The goals of the present study were to generate smaller mini-intein derivatives and to understand the basis for reversal of the splicing defect by the V67L mutation. Guided by the structural information, we generated a number of derivatives 135 to 152 residues in length, with V67 or L67. All of the new minimal inteins are functional in splicing. In vivo selection experiments for function showed that by removal of the loop region, 137 residues may be the lower limit for full protein-splicing activity. In addition, the activation effect of the V67L mutation was observed to be universal for mini-inteins longer than 137 residues. Structural and functional analyses indicate that the role of the mutation is in stabilization of the mini-intein core.

Family 18 chitinase-oligosaccharide substrate interaction: subsite preference and anomer selectivity of Serratia marcescens chitinase A.

The sizes and anomers of the products formed during the hydrolysis of chitin oligosaccharides by the Family 18 chitinase A (ChiA) from Serratia marcescens were analysed by hydrophilic interaction chromatography using a novel approach in which reactions were performed at 0 degrees C to stabilize the anomer conformations of the initial products. Crystallographic studies of the enzyme, having the structure of the complex of the ChiA E315L (Glu315-->Leu) mutant with a hexasaccharide, show that the oligosaccharide occupies subsites -4 to +2 in the substrate-binding cleft, consistent with the processing of beta-chitin by the release of disaccharide at the reducing end. Products of the hydrolysis of hexa- and penta-saccharides by wild-type ChiA, as well as by two mutants of the residues Trp275 and Phe396 important in binding the substrate at the +1 and +2 sites, show that the substrates only occupy sites -2 to +2 and that additional N -acetyl-D-glucosamines extend beyond the substrate-binding cleft at the reducing end. The subsites -3 and -4 are not used in this four-site binding mode. The explanation for these results is found in the high importance of individual binding sites for the processing of short oligosaccharides compared with the cumulative recognition and processive hydrolysis mechanism used to digest natural beta-chitin.

Green-lighting green fluorescent protein: faster and more efficient folding by eliminating a cis-trans peptide isomerization event.

Wild-type green fluorescent protein (GFP) folds on a time scale of minutes. The slow step in folding is a cis-trans peptide bond isomerization. The only conserved cis-peptide bond in the native GFP structure, at P89, was remodeled by the insertion of two residues, followed by iterative energy minimization and side chain design. The engineered GFP was synthesized and found to fold faster and more efficiently than its template protein, recovering 50% more of its fluorescence upon refolding. The slow phase of folding is faster and smaller in amplitude, and hysteresis in refolding has been eliminated. The elimination of a previously reported kinetically trapped state in refolding suggests that X-P89 is trans in the trapped state. A 2.55 Å resolution crystal structure revealed that the new variant contains only trans-peptide bonds, as designed. This is the first instance of a computationally remodeled fluorescent protein that folds faster and more efficiently than wild type.

Structure of catalytically competent intein caught in a redox trap with functional and evolutionary implications.

Here we describe self-splicing proteins, called inteins, that function as redox-responsive switches in bacteria. Redox regulation was achieved by engineering a disulfide bond between the intein's catalytic cysteine and a cysteine in the flanking 'extein' sequence. This interaction was validated by an X-ray structure, which includes a transient splice junction. A natural analog of the designed system was identified in Pyrococcus abyssi, suggesting an unprecedented form of adaptive, post-translational regulation.

Selection and structure of hyperactive inteins: peripheral changes relayed to the catalytic center.

Inteins are phylogenetically diverse self-splicing proteins that are of great functional, evolutionary, biotechnological, and medical interest. To address the relationship between intein structure and function, particularly with respect to regulating the splicing reaction, and to groom inteins for application, we developed a phage display system to extend current in vivo selection for enhanced intein function to selection in vitro. We thereby isolated inteins that can function under excursions in temperature, pH, and denaturing environment. Remarkably, most mutations mapped to the surface of the intein, remote from the active site. We chose two mutants with enhanced splicing activity for crystallography, one of which was also subjected to NMR analysis. These studies define a "ripple effect", whereby mutations in peripheral non-catalytic residues can cause subtle allosteric changes in the active-site environment in a way that facilitates intein activity. Altered salt-bridge formation and chemical shift changes of the mutant inteins provide a molecular rationale for their phenotypes. These fundamental insights will advance the utility of inteins in chemical biology, biotechnology, and medicine.

1H, 13C, and 15N NMR assignments of an engineered intein based on Mycobacterium tuberculosis RecA.

The backbone and side chain resonance assignments of an engineered intein based on Mycobacterium tuberculosis RecA have been determined based on triple-resonance experiments with the uniformly [(13)C,(15)N]-labeled protein.

Mechanism for intein C-terminal cleavage: a proposal from quantum mechanical calculations.

Inteins are autocatalytic protein cleavage and splicing elements. A cysteine to alanine mutation at the N-terminal of inteins inhibits splicing and isolates the C-terminal cleavage reaction. Experiments indicate an enhanced C-terminal cleavage reaction rate upon decreasing the solution pH for the cleavage mutant, which cannot be explained by the existing mechanistic framework. We use intein crystal structure data and the information about conserved amino acids to perform semiempirical PM3 calculations followed by high-level density functional theory calculations in both gas phase and implicit solvent environments. Based on these calculations, we propose a detailed "low pH" mechanism for intein C-terminal cleavage. Water plays an important role in the proposed reaction mechanism, acting as an acid as well as a base. The protonation of the scissile peptide bond nitrogen by a hydronium ion is an important first step in the reaction. That step is followed by the attack of the C-terminal asparagine side chain on its carbonyl carbon, causing succinimide formation and simultaneous peptide bond cleavage. The computed reaction energy barrier in the gas phase is approximately 33 kcal/mol and reduces to approximately 25 kcal/mol in solution, close to the 21 kcal/mol experimentally observed at pH 6.0. This mechanism is consistent with the observed increase in C-terminal cleavage activity at low pH for the cleavage mutant of the Mycobacterium tuberculosis RecA mini-intein.

The structure of the Candida albicans Ess1 prolyl isomerase reveals a well-ordered linker that restricts domain mobility.

Ess1 is a peptidyl-prolyl cis/trans isomerase (PPIase) that binds to the carboxy-terminal domain (CTD) of RNA polymerase II. Ess1 is thought to function by inducing conformational changes in the CTD that control the assembly of cofactor complexes on the transcription unit. Ess1 (also called Pin1) is highly conserved throughout the eukaryotic kingdom and is required for growth in some species, including the human fungal pathogen Candida albicans. Here we report the crystal structure of the C. albicansEss1 protein, determined at 1.6 A resolution. The structure reveals two domains, the WW and the isomerase domain, that have conformations essentially identical to those of human Pin1. However, the linker region that joins the two domains is quite different. In human Pin1, this linker is short and flexible, and part of it is unstructured. In contrast, the fungal Ess1 linker is highly ordered and contains a long alpha-helix. This structure results in a rigid juxtaposition of the WW and isomerase domains, in an orientation that is distinct from that observed in Pin1, and that eliminates a hydrophobic pocket between the domains that was implicated as the main substrate recognition site. These differences suggest distinct modes of interaction with long substrate molecules, such as the CTD of RNA polymerase II. We also show that C. albicans ess1(-)() mutants are attenuated for in vivo survival in mice. Together, these results suggest that CaEss1 might constitute a useful antifungal drug target, and that structural differences between the fungal and human enzymes could be exploited for drug design.

Three-dimensional structure of the R115E mutant of T4-bacteriophage 2'-deoxycytidylate deaminase.

2'-Deoxycytidylate deaminase (dCD) converts deoxycytidine 5'-monophosphate (dCMP) to deoxyuridine 5'-monophosphate and is a major supplier of the substrate for thymidylate synthase, an important enzyme in DNA synthesis and a major target for cancer chemotherapy. Wild-type dCD is allosterically regulated by the end products of its metabolic pathway, deoxycytidine 5'-triphosphate and deoxythymidine 5'-triphosphate, which act as an activator and an inhibitor, respectively. The first crystal structure of a dCD, in the form of the R115E mutant of the T4-bacteriophage enzyme complexed with the active site inhibitor pyrimidin-2-one deoxyribotide, has been determined at 2.2 A resolution. This mutant of dCD is active, even in the absence of the allosteric regulators. The molecular topology of dCD is related to that of cytidine deaminase (CDA) but with modifications for formation of the binding site for the phosphate group of dCMP. The enzyme has a zinc ion-based mechanism that is similar to that of CDA. A second zinc ion that is present in bacteriophage dCD, but absent in mammalian dCD and CDA, is important for the structural integrity of the enzyme and for the binding of the phosphate group of the substrate or inhibitor. Although the R115E mutant of dCD is a dimer in solution, it crystallizes as a hexamer, mimicking the natural state of the wild-type enzyme. Residues 112 and 115, which are known to be important for the binding of the allosteric regulators, are found in a pocket that is at the intersubunit interfaces in the hexamer but distant from the substrate-binding site. The substrate-binding site is composed of residues from a single protein molecule and is sequestered in a deep groove. This groove is located at the outer surface of the hexamer but ends at the subunit interface that also includes residue 115. It is proposed that the absence of subunit interactions at this interface in the dimeric R115E mutant renders the substrate-binding site accessible. In contrast, for the wild-type enzyme, binding of dCTP induces an allosteric effect that affects the subunit interactions and results in an increase in the accessibility of the binding site.

Catalytic domain structure and hypothesis for function of GIY-YIG intron endonuclease I-TevI.

I-TevI, a member of the GIY-YIG family of homing endonucleases, consists of an N-terminal catalytic domain and a C-terminal DNA-binding domain joined by a flexible linker. The GIY-YIG motif is in the N-terminal domain of I-TevI, which corresponds to a phylogenetically widespread catalytic cartridge that is often associated with mobile genetic elements. The crystal structure of the catalytic domain of I-TevI, the first of any GIY-YIG endonuclease, reveals a novel alpha/beta-fold with a central three-stranded antiparallel beta-sheet flanked by three helices. The most conserved and putative catalytic residues are located on a shallow, concave surface and include a metal coordination site. Similarities in the three-dimensional arrangement of the catalytically important residues and the cation-binding site with those of the His-Cys box endonuclease I-PpoI suggest the possibility of mechanistic relationships among these different families of homing endonucleases despite completely different folds.

Three-dimensional structure of human gamma -glutamyl hydrolase. A class I glatamine amidotransferase adapted for a complex substate.

gamma-Glutamyl hydrolase catalyzes the cleavage of the gamma-glutamyl chain of folylpoly-gamma-glutamyl substrates and is a central enzyme in folyl and antifolyl poly-gamma-glutamate metabolism. The crystal structure of human gamma-glutamyl hydrolase, determined at 1.6-A resolution, reveals that the protein is a homodimer. The overall structure of human gamma-glutamyl hydrolase contains 11 alpha-helices and 14 beta-strands, with a fold in which a central eight-stranded beta-sheet is sandwiched by three and five alpha-helices on each side. The topology is very similar to that of the class I glutamine amidotransferase domains, with the only major differences consisting of extensions in four loops and at the C terminus. These insertions are important for defining the substrate binding cleft and/or the dimer interface. Two sequence motifs are found in common between human gamma-glutamyl hydrolase and the class I glutamine amidotransferase family and include the catalytically essential residues, Cys-110 and His-220. These residues are located in the center of a large l-shaped cleft that is closed at one end and open at the other. Several conserved residues, including Glu-114, His-171, Gln-218, and Lys-223, may be important for substrate binding. Modeling of a methotrexate thioester intermediate, based on the corresponding complex of the glutamate thioester intermediate of Escherichia coli carbamoyl-phosphate synthetase, indicates that the substrate binds in an orientation with the pteroyl group toward the open end of the cleft.

Zinc finger as distance determinant in the flexible linker of intron endonuclease I-TevI.

I-TevI, the phage T4 td intron-encoded endonuclease, recognizes a lengthy DNA target and initiates intron mobility by introducing a double-strand break in the homing site. The enzyme uses both sequence and distance determinants to cleave the DNA 23-25 bp upstream of the intron insertion site. I-TevI consists of an N-terminal catalytic domain and a C-terminal DNA-binding domain separated by a long, flexible linker. The DNA-binding domain consists of three subdomains: a zinc finger, a minor-groove binding alpha-helix, and a helix-turn-helix. In this study, a mutational analysis was undertaken to assess the roles of these subdomains in substrate binding and cleavage. Surprisingly, the zinc finger is not required for DNA binding or catalysis. Rather, the zinc finger is a component of the linker and directs the catalytic domain to cleave the homing site at a fixed distance from the intron insertion site. When the cleavage site (CS) is shifted outside a given range, wild-type I-TevI defaults to the fixed distance, whereas zinc-finger mutants have lost the distance determinant and search out the displaced cleavage sequences. Although counterintuitive, a protein containing a 19-aa deletion of the zinc finger can extend further than can wild-type I-TevI to cleave a distant CS sequence, and a Cys-to-Ala mutant of the ligands for zinc, nominally a longer protein, can retract to cleave at a closer CS sequence. Models are presented for the novel function of the zinc finger, as a molecular constraint, whereby intramolecular protein-protein interactions position the catalytic domain by "catalytic clamp" and/or "linker-organizer" mechanisms.

Study of the structural dynamics of the E coli 70S ribosome using real-space refinement.

Cryo-EM density maps showing the 70S ribosome of E. coli in two different functional states related by a ratchet-like motion were analyzed using real-space refinement. Comparison of the two resulting atomic models shows that the ribosome changes from a compact structure to a looser one, coupled with the rearrangement of many of the proteins. Furthermore, in contrast to the unchanged inter-subunit bridges formed wholly by RNA, the bridges involving proteins undergo large conformational changes following the ratchet-like motion, suggesting an important role of ribosomal proteins in facilitating the dynamics of translation.