Outcomes of the EMDataResource cryo-EM Ligand Modeling Challenge.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

RIBFIND2: Identifying rigid bodies in protein and nucleic acid structures.

Molecular structures are often fitted into cryo-EM maps by flexible fitting. When this requires large conformational changes, identifying rigid bodies can help optimize the model-map fit. Tools for identifying rigid bodies in protein structures exist, however an equivalent for nucleic acid structures is lacking. With the increase in cryo-EM maps containing RNA and progress in RNA structure prediction, there is a need for such tools. We previously developed RIBFIND, a program for clustering protein secondary structures into rigid bodies. In RIBFIND2, this approach is extended to nucleic acid structures. RIBFIND2 can identify biologically relevant rigid bodies in important groups of complex RNA structures, capturing a wide range of dynamics, including large rigid-body movements. The usefulness of RIBFIND2-assigned rigid bodies in cryo-EM model refinement was demonstrated on three examples, with two conformations each: Group II Intron complexed IEP, Internal Ribosome Entry Site and the Processome, using cryo-EM maps at 2.7-5 Å resolution. A hierarchical refinement approach, performed on progressively smaller sets of RIBFIND2 rigid bodies, was clearly shown to have an advantage over classical all-atom refinement. RIBFIND2 is available via a web server with structure visualization and as a standalone tool.

Cryo-EM model validation recommendations based on outcomes of the 2019 EMDataResource challenge.

This paper describes outcomes of the 2019 Cryo-EM Model Challenge. The goals were to (1) assess the quality of models that can be produced from cryogenic electron microscopy (cryo-EM) maps using current modeling software, (2) evaluate reproducibility of modeling results from different software developers and users and (3) compare performance of current metrics used for model evaluation, particularly Fit-to-Map metrics, with focus on near-atomic resolution. Our findings demonstrate the relatively high accuracy and reproducibility of cryo-EM models derived by 13 participating teams from four benchmark maps, including three forming a resolution series (1.8 to 3.1 Å). The results permit specific recommendations to be made about validating near-atomic cryo-EM structures both in the context of individual experiments and structure data archives such as the Protein Data Bank. We recommend the adoption of multiple scoring parameters to provide full and objective annotation and assessment of the model, reflective of the observed cryo-EM map density.

Frequency and phenotypic spectrum of KMT2B dystonia in childhood: A single-center cohort study.

BACKGROUND: Childhood-onset dystonia is often genetically determined. Recently, KMT2B variants have been recognized as an important cause of childhood-onset dystonia. OBJECTIVE: To define the frequency of KMT2B mutations in a cohort of dystonic patients aged <18 years at onset, the associated clinical and radiological phenotype, and the natural history of disease. METHODS: Whole-exome sequencing or customized gene panels were used to screen a cohort of 65 patients who had previously tested negative for all other known dystonia-associated genes. RESULTS: We identified 14 patients (21.5%) carrying KMT2B variants, of which 1 was classified as a variant of unknown significance. We also identified 2 additional patients carrying pathogenic mutations in GNAO1 and ATM. Overall, we established a definitive genetic diagnosis in 23% of cases. We observed a spectrum of clinical manifestations in KMT2B variant carriers, ranging from generalized dystonia to short stature or intellectual disability alone, even within the same family. In 78.5% of cases, dystonia involved the lower limbs at onset, with later caudocranial generalization. Eight patients underwent pallidal DBS with a median decrease of Burke-Fahn-Marsden Dystonia Rating Scale-Motor score of 38.5% in the long term. We also report on 4 asymptomatic carriers, suggesting that some KMT2B mutations may be associated with incomplete disease penetrance. CONCLUSIONS: KMT2B mutations are frequent in childhood-onset dystonia and cause a complex neurodevelopmental syndrome, often featuring growth retardation and intellectual disability as additional phenotypic features. A dramatic and long-lasting response to DBS is characteristic of DYT-KMT2B dystonia. © 2019 International Parkinson and Movement Disorder Society.

Use of evolutionary information in the fitting of atomic level protein models in low resolution cryo-EM map of a protein assembly improves the accuracy of the fitting.

Protein-protein interface residues, especially those at the core of the interface, exhibit higher conservation than residues in solvent exposed regions. Here, we explore the ability of this differential conservation to evaluate fittings of atomic models in low-resolution cryo-EM maps and select models from the ensemble of solutions that are often proposed by different model fitting techniques. As a prelude, using a non-redundant and high-resolution structural dataset involving 125 permanent and 95 transient complexes, we confirm that core interface residues are conserved significantly better than nearby non-interface residues and this result is used in the cryo-EM map analysis. From the analysis of inter-component interfaces in a set of fitted models associated with low-resolution cryo-EM maps of ribosomes, chaperones and proteasomes we note that a few poorly conserved residues occur at interfaces. Interestingly a few conserved residues are not in the interface, though they are close to the interface. These observations raise the potential requirement of refitting the models in the cryo-EM maps. We show that sampling an ensemble of models and selection of models with high residue conservation at the interface and in good agreement with the density helps in improving the accuracy of the fit. This study indicates that evolutionary information can serve as an additional input to improve and validate fitting of atomic models in cryo-EM density maps.

Genome sequencing of herb Tulsi (Ocimum tenuiflorum) unravels key genes behind its strong medicinal properties.

BACKGROUND: Krishna Tulsi, a member of Lamiaceae family, is a herb well known for its spiritual, religious and medicinal importance in India. The common name of this plant is 'Tulsi' (or 'Tulasi' or 'Thulasi') and is considered sacred by Hindus. We present the draft genome of Ocimum tenuiflurum L (subtype Krishna Tulsi) in this report. The paired-end and mate-pair sequence libraries were generated for the whole genome sequenced with the Illumina Hiseq 1000, resulting in an assembled genome of 374 Mb, with a genome coverage of 61 % (612 Mb estimated genome size). We have also studied transcriptomes (RNA-Seq) of two subtypes of O. tenuiflorum, Krishna and Rama Tulsi and report the relative expression of genes in both the varieties. RESULTS: The pathways leading to the production of medicinally-important specialized metabolites have been studied in detail, in relation to similar pathways in Arabidopsis thaliana and other plants. Expression levels of anthocyanin biosynthesis-related genes in leaf samples of Krishna Tulsi were observed to be relatively high, explaining the purple colouration of Krishna Tulsi leaves. The expression of six important genes identified from genome data were validated by performing q-RT-PCR in different tissues of five different species, which shows the high extent of urosolic acid-producing genes in young leaves of the Rama subtype. In addition, the presence of eugenol and ursolic acid, implied as potential drugs in the cure of many diseases including cancer was confirmed using mass spectrometry. CONCLUSIONS: The availability of the whole genome of O.tenuiflorum and our sequence analysis suggests that small amino acid changes at the functional sites of genes involved in metabolite synthesis pathways confer special medicinal properties to this herb.

Inhibition of the interaction between NS3 protease and HCV IRES with a small peptide: a novel therapeutic strategy.

Recently, we have demonstrated that the protease domain of NS3 alone can bind specifically to hepatitis C virus (HCV) internal ribosome entry site (IRES) near the initiator AUG, dislodges human La protein and inhibits translation in favor of viral RNA replication. Here, by using a computational approach, the contact points of the protease on the HCV IRES were putatively mapped. A 30-mer NS3 peptide was designed from the predicted RNA-binding region that retained RNA-binding ability and also inhibited IRES-mediated translation. This peptide was truncated to 15 mer and this also demonstrated ability to inhibit HCV RNA-directed translation as well as replication. More importantly, its activity was tested in an in vivo mouse model by encapsulating the peptide in Sendai virus virosomes followed by intravenous delivery. The study demonstrates for the first time that the HCV NS3-IRES RNA interaction can be selectively inhibited using a small peptide and reports a strategy to deliver the peptide into the liver.

Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Mosaic RAS/MAPK variants cause sporadic vascular malformations which respond to targeted therapy.

BACKGROUND: Sporadic vascular malformations (VMs) are complex congenital anomalies of blood vessels that lead to stroke, life-threatening bleeds, disfigurement, overgrowth, and/or pain. Therapeutic options are severely limited, and multidisciplinary management remains challenging, particularly for high-flow arteriovenous malformations (AVM). METHODS: To investigate the pathogenesis of sporadic intracranial and extracranial VMs in 160 children in which known genetic causes had been excluded, we sequenced DNA from affected tissue and optimized analysis for detection of low mutant allele frequency. RESULTS: We discovered multiple mosaic-activating variants in 4 genes of the RAS/MAPK pathway, KRAS, NRAS, BRAF, and MAP2K1, a pathway commonly activated in cancer and responsible for the germline RAS-opathies. These variants were more frequent in high-flow than low-flow VMs. In vitro characterization and 2 transgenic zebrafish AVM models that recapitulated the human phenotype validated the pathogenesis of the mutant alleles. Importantly, treatment of AVM-BRAF mutant zebrafish with the BRAF inhibitor vemurafinib restored blood flow in AVM. CONCLUSION: Our findings uncover a major cause of sporadic VMs of different clinical types and thereby offer the potential of personalized medical treatment by repurposing existing licensed cancer therapies. FUNDING: This work was funded or supported by grants from the AVM Butterfly Charity, the Wellcome Trust (UK), the Medical Research Council (UK), the UK National Institute for Health Research, the L'Oreal-Melanoma Research Alliance, the European Research Council, and the National Human Genome Research Institute (US).

Mutations in the histone methyltransferase gene KMT2B cause complex early-onset dystonia.

Histone lysine methylation, mediated by mixed-lineage leukemia (MLL) proteins, is now known to be critical in the regulation of gene expression, genomic stability, cell cycle and nuclear architecture. Despite MLL proteins being postulated as essential for normal development, little is known about the specific functions of the different MLL lysine methyltransferases. Here we report heterozygous variants in the gene KMT2B (also known as MLL4) in 27 unrelated individuals with a complex progressive childhood-onset dystonia, often associated with a typical facial appearance and characteristic brain magnetic resonance imaging findings. Over time, the majority of affected individuals developed prominent cervical, cranial and laryngeal dystonia. Marked clinical benefit, including the restoration of independent ambulation in some cases, was observed following deep brain stimulation (DBS). These findings highlight a clinically recognizable and potentially treatable form of genetic dystonia, demonstrating the crucial role of KMT2B in the physiological control of voluntary movement.

Protein flexibility in the light of structural alphabets.

Protein structures are valuable tools to understand protein function. Nonetheless, proteins are often considered as rigid macromolecules while their structures exhibit specific flexibility, which is essential to complete their functions. Analyses of protein structures and dynamics are often performed with a simplified three-state description, i.e., the classical secondary structures. More precise and complete description of protein backbone conformation can be obtained using libraries of small protein fragments that are able to approximate every part of protein structures. These libraries, called structural alphabets (SAs), have been widely used in structure analysis field, from definition of ligand binding sites to superimposition of protein structures. SAs are also well suited to analyze the dynamics of protein structures. Here, we review innovative approaches that investigate protein flexibility based on SAs description. Coupled to various sources of experimental data (e.g., B-factor) and computational methodology (e.g., Molecular Dynamic simulation), SAs turn out to be powerful tools to analyze protein dynamics, e.g., to examine allosteric mechanisms in large set of structures in complexes, to identify order/disorder transition. SAs were also shown to be quite efficient to predict protein flexibility from amino-acid sequence. Finally, in this review, we exemplify the interest of SAs for studying flexibility with different cases of proteins implicated in pathologies and diseases.

3PFDB+: improved search protocol and update for the identification of representatives of protein sequence domain families.

Protein domain families are usually classified on the basis of similarity of amino acid sequences. Selection of a single representative sequence for each family provides targets for structure determination or modeling and also enables fast sequence searches to associate new members to a family. Such a selection could be challenging since some of these domain families exhibit huge variation depending on the number of members in the family, the average family sequence length or the extent of sequence divergence within a family. We had earlier created 3PFDB database as a repository of best representative sequences, selected from each PFAM domain family on the basis of high coverage. In this study, we have improved the database using more efficient strategies for the initial generation of sequence profiles and implement two independent methods, FASSM and HMMER, for identifying family members. HMMER employs a global sequence similarity search, while FASSM relies on motif identification and matching. This improved and updated database, 3PFDB+ generated in this study, provides representative sequences and profiles for PFAM families, with 13 519 family representatives having more than 90% family coverage. The representative sequence is also highlighted in a two-dimensional plot, which reflects the relative divergence between family members. Representatives belonging to small families with short sequences are mainly associated with low coverage. The set of sequences not recognized by the family representative profiles, highlight several potential false or weak family associations in PFAM. Partial domains and fragments dominate such cases, along with sequences that are highly diverged or different from other family members. Some of these outliers were also predicted to have different secondary structure contents, which reflect different putative structure or functional roles for these domain sequences. Database URL: http://caps.ncbs.res.in/3pfdbplus/.