Persistent SARS-CoV-2 infections may act as viral reservoirs that could seed future outbreaks1-5, give rise to highly divergent lineages6-8 and contribute to cases with post-acute COVID-19 sequelae (long COVID)9,10. However, the population prevalence of persistent infections, their viral load kinetics and evolutionary dynamics over the course of infections remain largely unknown. Here, using viral sequence data collected as part of a national infection survey, we identified 381 individuals with SARS-CoV-2 RNA at high titre persisting for at least 30 days, of which 54 had viral RNA persisting at least 60 days. We refer to these as \'persistent infections\' as available evidence suggests that they represent ongoing viral replication, although the persistence of non-replicating RNA cannot be ruled out in all. Individuals with persistent infection had more than 50% higher odds of self-reporting long COVID than individuals with non-persistent infection. We estimate that 0.1-0.5% of infections may become persistent with typically rebounding high viral loads and last for at least 60 days. In some individuals, we identified many viral amino acid substitutions, indicating periods of strong positive selection, whereas others had no consensus change in the sequences for prolonged periods, consistent with weak selection. Substitutions included mutations that are lineage defining for SARS-CoV-2 variants, at target sites for monoclonal antibodies and/or are commonly found in immunocompromised people11-14. This work has profound implications for understanding and characterizing SARS-CoV-2 infection, epidemiology and evolution.

Estimating parameters of amino acid substitution models is a crucial task in bioinformatics. The maximum likelihood (ML) approach has been proposed to estimate amino acid substitution models from large datasets. The quality of newly estimated models is normally assessed by comparing with the existing models in building ML trees. Two important questions remained are the correlation of the estimated models with the true models and the required size of the training datasets to estimate reliable models. In this article, we performed a simulation study to answer these two questions based on simulated data. We simulated genome datasets with different numbers of genes/alignments based on predefined models (called true models) and predefined trees (called true trees). The simulated datasets were used to estimate amino acid substitution model using the ML estimation methods. Our experiments showed that models estimated by the ML methods from simulated datasets with more than 100 genes have high correlations with the true models. The estimated models performed well in building ML trees in comparison with the true models. The results suggest that amino acid substitution models estimated by the ML methods from large genome datasets are a reliable tool for analyzing amino acid sequences.

The nucleocapsid component of SARS-CoV2 is involved in the viral genome packaging. GammaP.1(Brazil) and the 20 C-US(USA) variants had a high frequency of the P80R and P67S mutations respectively in the RNA-binding domain of the nucleocapsid. Since RNA-binding domain participates in the electrostatic interactions with the viral genome, the study of the effects of proline substitutions on the flexibility of the protein will be meaningful. It evinced that the trajectory of the wildtype and mutants was stable during the simulation and exhibited distinct changes in the flexibility of the protein. Moreover, the beta-hairpin loop region of the protein structures exhibited high amplitude fluctuations and dominant motions. Additionally, modulations were detected in the drug binding site. Besides, the extent of correlation and anti-correlation motions involving the protruding region, helix, and the other RNA binding sites differed between the wildtype and mutants. The secondary structure analysis disclosed the variation in the occurrence pattern of the secondary structure elements between the proteins. Protein-ssRNA interaction analysis was also done to detect the amino acid contacts with ssRNA. R44, R59, and Y61 residues of the wildtype and P80R mutant exhibited different duration contacts with the ssRNA. It was also noticed that R44, R59, and Y61 of the wildtype and P80R formed hydrogen bonds with the ssRNA. However in P67S, residues T43, R44, R45, R40, R59, and R41 displayed contacts and formed hydrogen bonds with ssRNA. Binding free energy was also calculated and was lowest for P67S than wildtype andP80R. Thus, proline substitutions influence the structure of the RNA-binding domain and may modulate viral genome packaging besides the host-immune response.Communicated by Ramaswamy H. Sarma.

Steroid 11β-hydroxylase deficiency (11β-OHD) is a rare autosomal recessive disorder caused by pathogenic variants of CYP11B1 gene. This study aimed to perform molecular analysis of a Chinese 11β-OHD series and in vitro functional study of twenty CYP11B1 missense variants. Twelve Chinese patients with clinical diagnosis of 11β-OHD were included in the study to analyze their molecular etiology. Genomic DNA of patients was extracted to be sequenced all coding exons and intronic flanking sequences of CYP11B1. Fourteen missense variants found in 12 patients mentioned above along with 6 missense variants previously reported by our team were evaluated functionally. Amino acid substitutions were analyzed with computational program to determine their effects on the three-dimensional structure of CYP11B1 protein. Clinical characteristics and hormone levels at baseline of the 18 patients carrying 18 missense variants aforementioned were recorded to perform genotype-phenotype correlation. A total of 21 rare variants including 9 novel and 12 recurrent ones were identified in 12 patients, out of which 17 were missense, 2 were nonsense, 1 was a splice site variant, and 1 was a deletion-insertion variant. Results of in vitro functional study revealed that 3 out of 20 missense mutants (p.Leu3Pro, p.Gly267Ser, and p.Ala367Ser) had partial enzyme activity and the other 17 had little enzymatic activity. The impairment degree of enzymatic activity in vitro functional study was also reflected in the severity degree of interaction change between the wild-type/mutant-type amino acid and its adjacent amino acids in three-dimensional model. In conclusion, the addition of 9 novel variants expands the spectrum of CYP11B1 pathogenic variants. Our results demonstrate that twenty CYP11B1 variants lead to impaired 11β-hydroxylase activity in vitro. Visualizing these variants in the three-dimensional model structure of CYP11B1 protein can provide a plausible explanation for the results measured in vitro.

UNASSIGNED: Here we report two rare α-globin chain variants in two unrelated families: Hb Val de Marne [α133(H16) Ser > Arg (AGC > CGC); HBA2: c.400A > C] and Hb Dongguan [α52(E6) Ser > Cys (TCT > TGT); HBA1: c.158C > G]. Notably, HBA2: c.400A > C is an unreported new variant in the third exon of the α2 gene, and simple heterozygous unstable Hb Dongguan haematological characteristics are proposed for the first time.
UNASSIGNED: Hb analysis was performed by using capillary electrophoresis (CE). Twenty-three common mutations were detected using a suspension array system. Mutations were identified by DNA sequencing.
UNASSIGNED: The CE results showed an abnormal peak with incomplete separation from Hb A at zone 8 in two members of Family 1. DNA sequencing confirmed the presence of Hb Val de Marne [α133(H16) Ser > Arg (AGC > CGC); HBA2: c.400A > C]. Five members of Family 2 exhibited an abnormal peak at zone 11, and DNA sequencing confirmed the presence of Hb Dongguan [α52(E6) Ser > Cys (TCT > TGT); HBA1: c.158C > G].
UNASSIGNED: The discovery of HBA2: C.400A > C expands the existing spectrum of α-globin variants. The carriers of simple heterozygous Hb Dongguan generally do not have obvious clinical symptoms. The information in this study will help clinicians understand the screening, molecular diagnosis and clinical significance of Hb variants.

Coronavirus disease 2019 (COVID-19) caused appalling conditions over the globe, which is currently faced by the entire human population. One of the primary reasons behind the uncontrollable situation is the lack of specific therapeutics. In such conditions, drug repurposing of available drugs (viz. Chloroquine, Lopinavir, etc.) has been proposed, but various clinical and preclinical investigations indicated the toxicity and adverse side effects of these drugs. This study explores the inhibition potency of phytochemicals from Tinospora cordifolia (Giloy) against SARS CoV-2 drugable targets (spike glycoprotein and Mpro proteins) using molecular docking and MD simulation studies. ADMET, virtual screening, MD simulation, postsimulation analysis (RMSD, RMSF, Rg, SASA, PCA, FES) and MM-PBSA calculations were carried out to predict the inhibition efficacy of the phytochemicals against SARS CoV-2 targets. Tinospora compounds showed better binding affinity than the corresponding reference. Their binding affinity ranges from -9.63 to -5.68 kcal/mole with spike protein and -10.27 to -7.25 kcal/mole with main protease. Further 100 ns exhaustive simulation studies and MM-PBSA calculations supported favorable and stable binding of them. This work identifies Nine Tinospora compounds as potential inhibitors. Among those, 7-desacetoxy-6,7-dehydrogedunin was found to inhibit both spike (7NEG) and Mpro (7MGS and 6LU7) proteins, and Columbin was found to inhibit selected spike targets (7NEG and 7NX7). In all the analyses, these compounds performed well and confirms the stable binding. Hence the identified compounds, advocated as potential inhibitors can be taken for further in vitro and in vivo experimental validation to determine their anti-SARS-CoV-2 potential.Communicated by Ramaswamy H. Sarma.

Respiratory syncytial virus (RSV) is the leading cause of lower respiratory tract infection in children and infants. To date, there is no effective vaccine available against RSV. Heparan sulfate is a type of glycosaminoglycan that aids in the attachment of the RSV to the host cell membrane via the G protein. In the present study, the effect of amino acid substitution on the structure and stability of the ectodomain G protein was studied. Further, it was investigated whether mutation (K117A) in the CX3C motif of G protein alters the binding with heparan sulfate. The point mutation significantly affects the conformational stability of the G protein. The mutant protein showed a low binding affinity with heparan sulfate as compared to the wild-type G protein, as determined by fluorescence quenching, isothermal titration calorimetry (ITC), and molecular docking studies. The low binding affinity and decreased stability suggested that this mutation may play an important role in prevention of attachment of virion to the host cell receptors. Collectively, this investigation suggests that mutation in the CX3C motif of G protein may likely improve the efficacy and safety of the RSV vaccine.

The COVID-19 pandemic continues to threaten healthcare systems worldwide due to the limited access to vaccines, suboptimal treatment options, and the continuous emergence of new and more transmissible SARS-CoV-2 variants. Reverse-genetics studies of viral genes and mutations have proven highly valuable in advancing basic virus research, leading to the development of therapeutics. We developed a functional and highly versatile full-length SARS-CoV-2 infectious system by cloning the sequence of a COVID-19 associated virus isolate (DK-AHH1) into a bacterial artificial chromosome (BAC). Viruses recovered after RNA-transfection of in vitro transcripts into Vero E6 cells showed growth kinetics and remdesivir susceptibility similar to the DK-AHH1 virus isolate. Insertion of reporter genes, green fluorescent protein, and nanoluciferase into the ORF7 genomic region led to high levels of reporter activity, which facilitated high throughput treatment experiments. We found that putative coronavirus remdesivir resistance-associated substitutions F480L and V570L-and naturally found polymorphisms A97V, P323L, and N491S, all in nsp12-did not decrease SARS-CoV-2 susceptibility to remdesivir. A nanoluciferase reporter clone with deletion of spike (S), envelope (E), and membrane (M) proteins exhibited high levels of transient replication, was inhibited by remdesivir, and therefore could function as an efficient non-infectious subgenomic replicon system. The developed SARS-CoV-2 reverse-genetics systems, including recombinants to modify infectious viruses and non-infectious subgenomic replicons with autonomous genomic RNA replication, will permit high-throughput cell culture studies-providing fundamental understanding of basic biology of this coronavirus. We have proven the utility of the systems in rapidly introducing mutations in nsp12 and studying their effect on the efficacy of remdesivir, which is used worldwide for the treatment of COVID-19. Our system provides a platform to effectively test the antiviral activity of drugs and the phenotype of SARS-CoV-2 mutants.

Penicillium echinulatum 2HH is an ascomycete well known for its production of cellulolytic enzymes. Understanding lignocellulolytic and sugar uptake systems is essential to obtain efficient fungi strains for the production of bioethanol. In this study we performed a genome-wide functional annotation of carbohydrate-active enzymes and sugar transporters involved in the lignocellulolytic system of P. echinulatum 2HH and S1M29 strains (wildtype and mutant, respectively) and eleven related fungi. Additionally, signal peptide and orthology prediction were carried out. We encountered a diverse assortment of cellulolytic enzymes in P. echinulatum, especially in terms of β-glucosidases and endoglucanases. Other enzymes required for the breakdown of cellulosic biomass were also found, including cellobiohydrolases, lytic cellulose monooxygenases and cellobiose dehydrogenases. The S1M29 mutant, which is known to produce an increased cellulase activity, and the 2HH wild type strain of P. echinulatum did not show significant differences between their enzymatic repertoire. Nevertheless, we unveiled an amino acid substitution for a predicted intracellular β-glucosidase of the mutant, which might contribute to hyperexpression of cellulases through a cellodextrin induction pathway. Most of the P. echinulatum enzymes presented orthologs in P. oxalicum 114-2, supporting the presence of highly similar cellulolytic mechanisms and a close phylogenetic relationship between these fungi. A phylogenetic analysis of intracellular β-glucosidases and sugar transporters allowed us to identify several proteins potentially involved in the accumulation of intracellular cellodextrins. These may prove valuable targets in the genetic engineering of P. echinulatum focused on industrial cellulases production. Our study marks an important step in characterizing and understanding the molecular mechanisms employed by P. echinulatum in the enzymatic hydrolysis of lignocellulosic biomass.

The influenza virus is an important respiratory pathogen that causes many incidences of diseases and even death each year. One of the primary factors of this virus is the Neuraminidase surface protein, which causes the virus to leave the host cell and spread to new target cells. The main antiviral medication for influenza is designed as a protein inhibitor ligand that prevents further spread of the disease, and eventually relieves the emerged symptoms. The effectiveness of such inhibitory drugs is highly associated with their binding affinity. In this paper, the binding affinity of an herbal ligand of Capsaicin bound to Neuraminidase of the influenza virus is investigated using steered molecular dynamics (SMD) simulation. Since mutations of the virus directly impact the binding affinity of the inhibitory drugs, different mutations were generated by using Mutagenesis module. The rapid spread of infection during the avian influenza A/H5N1 epidemic has raised concerns about far more dangerous consequences if the virus becomes resistant to current drugs. Currently, oseltamivir (Tamiflu), zanamivir (Relenza), pramivir (Rapivab), and laninamivir (Inavir) are increasingly used to treat the flu. However, with the rapid evolution of the virus, some drug-resistant strains are emerging. Therefore, it is very important to seek alternative therapies and identify the roots of drug resistance. Obtained results demonstrated a reduced binding affinity for the applied mutations. This reduction in binding affinity will cause the virus mutation to become resistant to the drug, which will spread the disease and make it more difficult to treat. From a molecular prospect, this decrease in binding affinity is due to the loss of a number of effective bonds between the ligand and the receptor, which occurs with mutations of the wild-type (WT) species. The results of the present study can be used in the rational design of novel drugs that are compatible with specific mutations.