Computational Insights into SARS-CoV-2 Main Protease Mutations and Nirmatrelvir Efficacy: The Effects of P132H and P132H-A173V.

Nirmatrelvir, a pivotal component of the oral antiviral Paxlovid for COVID-19, targets the SARS-CoV-2 main protease (Mpro) as a covalent inhibitor. Here, we employed combined computational methods to explore how the prevalent Omicron variant mutation P132H, alone and in combination with A173V (P132H-A173V), affects nirmatrelvir's efficacy. Our findings suggest that P132H enhances the noncovalent binding affinity of Mpro for nirmatrelvir, whereas P132H-A173V diminishes it. Although both mutants catalyze the rate-limiting step more efficiently than the wild-type (WT) Mpro, P132H slows the overall rate of covalent bond formation, whereas P132H-A173V accelerates it. Comprehensive analysis of noncovalent and covalent contributions to the overall binding free energy of the covalent complex suggests that P132H likely enhances Mpro sensitivity to nirmatrelvir, while P132H-A173V may confer resistance. Per-residue decompositions of the binding and activation free energies pinpoint key residues that significantly affect the binding affinity and reaction rates, revealing how the mutations modulate these effects. The mutation-induced conformational perturbations alter drug-protein local contact intensities and the electrostatic preorganization of the protein, affecting noncovalent binding affinity and the stability of key reaction states, respectively. Our findings inform the mechanisms of nirmatrelvir resistance and sensitivity, facilitating improved drug design and the detection of resistant strains.

Identification of and Mechanistic Insights into SARS-CoV-2 Main Protease Non-Covalent Inhibitors: An In-Silico Study.

The indispensable role of the SARS-CoV-2 main protease (Mpro) in the viral replication cycle and its dissimilarity to human proteases make Mpro a promising drug target. In order to identify the non-covalent Mpro inhibitors, we performed a comprehensive study using a combined computational strategy. We first screened the ZINC purchasable compound database using the pharmacophore model generated from the reference crystal structure of Mpro complexed with the inhibitor ML188. The hit compounds were then filtered by molecular docking and predicted parameters of drug-likeness and pharmacokinetics. The final molecular dynamics (MD) simulations identified three effective candidate inhibitors (ECIs) capable of maintaining binding within the substrate-binding cavity of Mpro. We further performed comparative analyses of the reference and effective complexes in terms of dynamics, thermodynamics, binding free energy (BFE), and interaction energies and modes. The results reveal that, when compared to the inter-molecular electrostatic forces/interactions, the inter-molecular van der Waals (vdW) forces/interactions are far more important in maintaining the association and determining the high affinity. Given the un-favorable effects of the inter-molecular electrostatic interactions-association destabilization by the competitive hydrogen bond (HB) interactions and the reduced binding affinity arising from the un-compensable increase in the electrostatic desolvation penalty-we suggest that enhancing the inter-molecular vdW interactions while avoiding introducing the deeply buried HBs may be a promising strategy in future inhibitor optimization.

Variances and covariances of linear summary statistics of segregating sites.

Each mutation in a population sample of DNA sequences can be classified by the number of sequences that inherit the mutant nucleotide, the resulting frequencies are known as mutations of different sizes or site frequency spectrum. Many summary statistics can be defined as a linear function of these frequencies. A flexible class of such linear summary statistics is explored analytically in this paper which include several well-known quantities, such as the number of segregating sizes and the mean number of nucleotide differences between two sequences. Some asymptotic variances and covariances are obtained while the analytical formulas for the variances and covariances of nine such linear summary statistics are derived, most of which are unknown to date. This study not only provides some theoretical foundations for exploring linear summary statistics, but also provides some newlinear summary statistics that may be utilized for analyzing sample polymorphism. Furthermore it is showed that a newly developed linear summary statistics has a smaller variance almost uniformly than Watterson's estimator, and that a class of linear summary statistics given too heavy weights on mutations of smaller sizes result in asymptotically non-zero variance.

Amitosis as a strategy of cell division-Insight from the proliferation of Tetrahymena thermophila macronuclei.

Cell division is a necessity of life which can be either mitotic or amitotic. While both are fundamental, amitosis is sometimes considered a relic of little importance in biology. Nevertheless, eukaryotes often have polyploid cells, including cancer cells, which may divide amitotically. To understand how amitosis ensures the completion of cell division, we turn to the macronuclei of ciliates. The grand scheme governing the proliferation of the macronuclei of ciliate cells, which involves chromosomal replication and amitosis, is currently unknown, which is crucial for developing population genetics model of ciliate populations. Using a novel model that encompasses a wide range of mechanisms together with experimental data of the composition of mating types at different stages derived from a single karyonide of Tetrahymena thermophila, we show that the chromosomal replication of the macronucleus has a strong head-start effect, with only about five copies of chromosomes replicated at a time and persistent reuse of the chromosomes involved in the early replication. Furthermore the fission of a fully grown macronucleus is non-random with regard to chromosome composition, with a strong tendency to push chromosomes and their replications to the same daughter cell.

Exploring the Cold-Adaptation Mechanism of Serine Hydroxymethyltransferase by Comparative Molecular Dynamics Simulations.

Cold-adapted enzymes feature a lower thermostability and higher catalytic activity compared to their warm-active homologues, which are considered as a consequence of increased flexibility of their molecular structures. The complexity of the (thermo)stability-flexibility-activity relationship makes it difficult to define the strategies and formulate a general theory for enzyme cold adaptation. Here, the psychrophilic serine hydroxymethyltransferase (pSHMT) from Psychromonas ingrahamii and its mesophilic counterpart, mSHMT from Escherichia coli, were subjected to µs-scale multiple-replica molecular dynamics (MD) simulations to explore the cold-adaptation mechanism of the dimeric SHMT. The comparative analyses of MD trajectories reveal that pSHMT exhibits larger structural fluctuations and inter-monomer positional movements, a higher global flexibility, and considerably enhanced local flexibility involving the surface loops and active sites. The largest-amplitude motion mode of pSHMT describes the trends of inter-monomer dissociation and enlargement of the active-site cavity, whereas that of mSHMT characterizes the opposite trends. Based on the comparison of the calculated structural parameters and constructed free energy landscapes (FELs) between the two enzymes, we discuss in-depth the physicochemical principles underlying the stability-flexibility-activity relationships and conclude that (i) pSHMT adopts the global-flexibility mechanism to adapt to the cold environment and, (ii) optimizing the protein-solvent interactions and loosening the inter-monomer association are the main strategies for pSHMT to enhance its flexibility.

Molecular dynamics simulations reveal distinct differences in conformational dynamics and thermodynamics between the unliganded and CD4-bound states of HIV-1 gp120.

The entry of human immunodeficiency virus type I (HIV-1) into host cells is initiated by binding to the cell-surface receptor CD4, which induces a conformational transition of the envelope (Env) glycoprotein gp120 from the closed, unliganded state to the open, CD4-bound state. Despite many available structures in these two states, detailed aspects on the dynamics and thermodynamics of gp120 remain elusive. Here, we performed microsecond-scale (µs-scale) multiple-replica molecular dynamics (MD) simulations to explore the differences in the conformational dynamics, protein motions, and thermodynamics between the unliganded and CD4-bound/complexed forms of gp120. Comparative analyses of MD trajectories reveal that CD4 binding promotes the structural deviations/changes and conformational flexibility, loosens the structural packing, and complicates the molecular motions of gp120. Comparison of the constructed free energy landscapes (FELs) reveals that the CD4-complexed gp120 has more conformational substates, larger conformational entropy, and lower thermostability than the unliganded form. Therefore, the unliganded conformation represents a structurally and energetically stable "ground state" for the full-length gp120. The observed great increase in the mobility of V1/V2 and V3 along with their more versatile movement directions in the CD4-bound gp120 compared to the unliganded form suggests that their orientations with respect to each other and to the structural core determine the differences in the conformational dynamics and thermodynamics between the two gp120 forms. The results presented here provide a basis by which to better understand the functional and immunological properties of gp120 and, furthermore, to deploy appropriate strategies for the development of anti-HIV-1 drugs or vaccines.

Large numbers of vertebrates began rapid population decline in the late 19th century.

Accelerated losses of biodiversity are a hallmark of the current era. Large declines of population size have been widely observed and currently 22,176 species are threatened by extinction. The time at which a threatened species began rapid population decline (RPD) and the rate of RPD provide important clues about the driving forces of population decline and anticipated extinction time. However, these parameters remain unknown for the vast majority of threatened species. Here we analyzed the genetic diversity data of nuclear and mitochondrial loci of 2,764 vertebrate species and found that the mean genetic diversity is lower in threatened species than in related nonthreatened species. Our coalescence-based modeling suggests that in many threatened species the RPD began ∼123 y ago (a 95% confidence interval of 20-260 y). This estimated date coincides with widespread industrialization and a profound change in global living ecosystems over the past two centuries. On average the population size declined by ∼25% every 10 y in a threatened species, and the population size was reduced to ∼5% of its ancestral size. Moreover, the ancestral size of threatened species was, on average, ∼22% smaller than that of nonthreatened species. Because the time period of RPD is short, the cumulative effect of RPD on genetic diversity is still not strong, so that the smaller ancestral size of threatened species may be the major cause of their reduced genetic diversity; RPD explains 24.1-37.5% of the difference in genetic diversity between threatened and nonthreatened species.

Common variants in SIRT1 and human longevity in a Chinese population.

BACKGROUND: The silent information regulator SIR2/SIRT1gene has been demonstrated as regulating lifespan in many model organisms, including yeast, worms, fruit flies and rodents. SIRT1, the human homolog of SIR2, is considered a candidate gene as a modifier of human life expectancy. METHODS: In the current study we included 616 long-lived individuals and 846 matched younger controls to investigate associations between 8 common single nucleotide polymorphisms (SNPs) (i.e., rs12778366, rs3758391, rs3740051, rs33957861, rs7896005, rs12413112, rs11599176 and rs4746720) in the SIRT1 gene and human longevity. RESULTS: The 8 SNPs had strong linkage disequilibrium (LD) and were in an LD block, which was characterized by 4 common haplotypes that capture 99.3% of the genetic variation present within it. We found no evidence for statistically significant associations between the tested SIRT1 SNPs and longevity at the allele, genotype or haplotype levels. CONCLUSIONS: Current findings show that several common variants in SIRT1 are not associated with human longevity.

Association of common variants in TOMM40/APOE/APOC1 region with human longevity in a Chinese population.

Apolipoprotein E (APOE), translocase of outer mitochondrial membrane 40 homolog (TOMM40) and apolipoprotein C-I (APOC1) may extend lifespan by marked delay or escape from age-related diseases. This study aimed to elucidate the association of human longevity with genetic variations in TOMM40/APOE/APOC1 region in a Chinese population. Ten tag single-nucleotide polymorphisms (SNPs) in the TOMM40/APOE/APOC1 region were successfully genotyped in 616 unrelated long-lived individuals and 846 younger controls. Of the 10 SNPs, rs7254892 in 5' upstream of TOMM40 showed significant association with human longevity (G/A-A/A vs G/G: odds ratio (OR)=1.59, 95% confidence interval (CI)=1.20-2.09, P=0.0011, Bonferroni corrected P (Pc)=0.033). The haplotype analysis suggested that individuals carrying the haplotype A-A-A-A-T-A-T-G-C-A (rs7254892-rs157580-rs2075649-rs2075650-rs157582-rs8106922-rs1160985-rs405697-rs439401-rs445925) tended to have longer lifespan than those carrying the most common haplotype G-G-A-A-C-A-C-A-T-G (OR=1.59, 95% CI=1.19-2.12, P=0.0018, Pc=0.0216). These findings indicated that variants in TOMM40/APOE/APOC1 region might be associated with human longevity. Further studies are needed to identify the causal genetic variants influencing human longevity.

Lack of association between polymorphisms in the SIRT6 gene and longevity in a Chinese population.

Sirtuin 6 (SIRT6) has recently been demonstrated to play an important role in the regulation of longevity in mammals. We therefore aimed to determine whether common variations in the SIRT6 gene are associated with human longevity. Five tag single nucleotide polymorphisms (SNPs) across the SIRT6 gene and its 5 kb up-/downstream region, including rs350852, rs350844, rs352493, rs4807546 and rs3760905, have been successfully determined in 616 unrelated Chinese long-lived individuals (LLIs) (mean age: 102.4 ± 2.3 years) and 846 younger controls (mean age: 48.9 ± 10.6 years) from Hainan Island, China. The allele and genotype frequencies of the five SNPs showed no statistically significant difference between the LLIs and controls (all P > 0.05). The five SNPs were in strong linkage disequilibrium and defined seven common haplotypes. Likewise, no association between these haplotypes and longevity was observed (all P > 0.05). The present study reveals that common genetic variations in the SIRT6 gene are not associated with human longevity.

Age-dependent transition from cell-level to population-level control in murine intestinal homeostasis revealed by coalescence analysis.

In multi-cellular organisms, tissue homeostasis is maintained by an exquisite balance between stem cell proliferation and differentiation. This equilibrium can be achieved either at the single cell level (a.k.a. cell asymmetry), where stem cells follow strict asymmetric divisions, or the population level (a.k.a. population asymmetry), where gains and losses in individual stem cell lineages are randomly distributed, but the net effect is homeostasis. In the mature mouse intestinal crypt, previous evidence has revealed a pattern of population asymmetry through predominantly symmetric divisions of stem cells. In this work, using population genetic theory together with previously published crypt single-cell data obtained at different mouse life stages, we reveal a strikingly dynamic pattern of stem cell homeostatic control. We find that single-cell asymmetric divisions are gradually replaced by stochastic population-level asymmetry as the mouse matures to adulthood. This lifelong process has important developmental and evolutionary implications in understanding how adult tissues maintain their homeostasis integrating the trade-off between intrinsic and extrinsic regulations.

Population-based resequencing of ANGPTL4 uncovers variations that reduce triglycerides and increase HDL.

Resequencing genes provides the opportunity to assess the full spectrum of variants that influence complex traits. Here we report the first application of resequencing to a large population (n = 3,551) to examine the role of the adipokine ANGPTL4 in lipid metabolism. Nonsynonymous variants in ANGPTL4 were more prevalent in individuals with triglyceride levels in the lowest quartile than in individuals with levels in the highest quartile (P = 0.016). One variant (E40K), present in approximately 3% of European Americans, was associated with significantly lower plasma levels of triglyceride and higher levels of high-density lipoprotein cholesterol in European Americans from the Atherosclerosis Risk in Communities Study and in Danes from the Copenhagen City Heart Study. The ratio of nonsynonymous to synonymous variants was higher in European Americans than in African Americans (4:1 versus 1.3:1), suggesting population-specific relaxation of purifying selection. Thus, resequencing of ANGPTL4 in a multiethnic population allowed analysis of the phenotypic effects of both rare and common variants while taking advantage of genetic variation arising from ethnic differences in population history.

Effect of the Solvent Temperatures on Dynamics of Serine Protease Proteinase K.

To obtain detailed information about the effect of the solvent temperatures on protein dynamics, multiple long molecular dynamics (MD) simulations of serine protease proteinase K with the solute and solvent coupled to different temperatures (either 300 or 180 K) have been performed. Comparative analyses demonstrate that the internal flexibility and mobility of proteinase K are strongly dependent on the solvent temperatures but weakly on the protein temperatures. The constructed free energy landscapes (FELs) at the high solvent temperatures exhibit a more rugged surface, broader spanning range, and higher minimum free energy level than do those at the low solvent temperatures. Comparison between the dynamic hydrogen bond (HB) numbers reveals that the high solvent temperatures intensify the competitive HB interactions between water molecules and protein surface atoms, and this in turn exacerbates the competitive HB interactions between protein internal atoms, thus enhancing the conformational flexibility and facilitating the collective motions of the protein. A refined FEL model was proposed to explain the role of the solvent mobility in facilitating the cascade amplification of microscopic motions of atoms and atomic groups into the global collective motions of the protein.

Persistent latent tuberculosis reactivation risk in United States immigrants.

RATIONALE: Current guidelines limit latent tuberculosis infection (LTBI) evaluation to persons in the United States less than or equal to 5 years based on the assumption that high TB rates among recent entrants are attributable to high LTBI reactivation risk, which declines over time. We hypothesized that high postarrival TB rates may instead be caused by imported active TB. OBJECTIVES: Estimate reactivation and imported TB in an immigrant cohort. METHODS: We linked preimmigration records from a cohort of California-bound Filipino immigrants during 2001-2010 with subsequent TB reports. TB was likely LTBI reactivation if the immigrant had no evidence of active TB at preimmigration examination, likely imported if preimmigration radiograph was abnormal and TB was reported less than or equal to 6 months after arrival, and likely reactivation of inactive TB if radiograph was abnormal but TB was reported more than 6 months after arrival. MEASUREMENTS AND MAIN RESULTS: Among 123,114 immigrants, 793 TB cases were reported. Within 1 year of preimmigration examination, 85% of TB was imported; 6 and 9% were reactivation of LTBI and inactive TB, respectively. Conversely, during Years 2-9 after U.S. entry, 76 and 24% were reactivation of LTBI and inactive TB, respectively. The rate of LTBI reactivation (32 per 100,000) did not decline during Years 1-9. CONCLUSIONS: High postarrival TB rates were caused by detection of imported TB through active postarrival surveillance. Among immigrants without active TB at baseline, reported TB did not decline over 9 years, indicating sustained high risk of LTBI reactivation. Revised guidelines should support LTBI screening and treatment more than 5 years after U.S. arrival.

Assessment of partial-mouth periodontal examination protocols for periodontitis surveillance.

OBJECTIVE: To evaluate bias associated with nine identified partial-mouth periodontal examination (PMPE) protocols in estimating periodontitis prevalence using the periodontitis case definition given by the Centers of Disease Control and Prevention and American Academy of Periodontology (CDC/AAP). MATERIAL AND METHODS: Prevalence from full-mouth examination was determined in a sample of 3667 adults ≥30 years old from the National Health and Nutrition Examination Survey (NHANES) 2009-2010. Prevalence, absolute bias, relative bias, sensitivity and inflation factor were derived for these protocols according to the CDC/AAP definition and half-reduced CDC/AAP definition as ≤50% of sites were measured. RESULTS: Bias in moderate and severe periodontitis prevalence ranged between 11.1-52.5% and 27.1-76.3% for full-mouth mesiobuccal-distolingual protocol and half-mouth mesiobuccal protocol respectively; according to the CDC/AAP definition. With half-reduced CDC/AAP definition, half-mouth four sites protocol provided small absolute bias (3.2%) and relative bias (9.3%) for the estimates of moderate periodontitis prevalence; corresponding biases for severe periodontitis were -1.2% and -10.2%. CONCLUSION: Periodontitis prevalence can be estimated with limited bias when a half-mouth four sites protocol and a half-reduced CDC/AAP case definition are used in combination.

Highly variable recessive lethal or nearly lethal mutation rates during germ-line development of male Drosophila melanogaster.

Each cell of higher organism adults is derived from a fertilized egg through a series of divisions, during which mutations can occur. Both the rate and timing of mutations can have profound impacts on both the individual and the population, because mutations that occur at early cell divisions will affect more tissues and are more likely to be transferred to the next generation. Using large-scale multigeneration screening experiments for recessive lethal or nearly lethal mutations of Drosophila melanogaster and recently developed statistical analysis, we show for male D. melanogaster that (i) mutation rates (for recessive lethal or nearly lethal) are highly variable during germ cell development; (ii) first cell cleavage has the highest mutation rate, which drops substantially in the second cleavage or the next few cleavages; (iii) the intermediate stages, after a few cleavages to right before spermatogenesis, have at least an order of magnitude smaller mutation rate; and (iv) spermatogenesis also harbors a fairly high mutation rate. Because germ-line lineage shares some (early) cell divisions with somatic cell lineage, the first conclusion is readily extended to a somatic cell lineage. It is conceivable that the first conclusion is true for most (if not all) higher organisms, whereas the other three conclusions are widely applicable, although the extent may differ from species to species. Therefore, conclusions or analyses that are based on equal mutation rates during development should be taken with caution. Furthermore, the statistical approach developed can be adopted for studying other organisms, including the human germ-line or somatic mutational patterns.

Estimating DNA polymorphism from next generation sequencing data with high error rate by dual sequencing applications.

BACKGROUND: As the error rate is high and the distribution of errors across sites is non-uniform in next generation sequencing (NGS) data, it has been a challenge to estimate DNA polymorphism (θ) accurately from NGS data. RESULTS: By computer simulations, we compare the two methods of data acquisition - sequencing each diploid individual separately and sequencing the pooled sample. Under the current NGS error rate, sequencing each individual separately offers little advantage unless the coverage per individual is high (>20X). We hence propose a new method for estimating θ from pooled samples that have been subjected to two separate rounds of DNA sequencing. Since errors from the two sequencing applications are usually non-overlapping, it is possible to separate low frequency polymorphisms from sequencing errors. Simulation results show that the dual applications method is reliable even when the error rate is high and θ is low. CONCLUSIONS: In studies of natural populations where the sequencing coverage is usually modest (~2X per individual), the dual applications method on pooled samples should be a reasonable choice.

Estimating population genetic parameters and comparing model goodness-of-fit using DNA sequences with error.

It is known that sequencing error can bias estimation of evolutionary or population genetic parameters. This problem is more prominent in deep resequencing studies because of their large sample size n, and a higher probability of error at each nucleotide site. We propose a new method based on the composite likelihood of the observed SNP configurations to infer population mutation rate theta = 4N(e)micro, population exponential growth rate R, and error rate epsilon, simultaneously. Using simulation, we show the combined effects of the parameters, theta, n, epsilon, and R on the accuracy of parameter estimation. We compared our maximum composite likelihood estimator (MCLE) of theta with other theta estimators that take into account the error. The results show the MCLE performs well when the sample size is large or the error rate is high. Using parametric bootstrap, composite likelihood can also be used as a statistic for testing the model goodness-of-fit of the observed DNA sequences. The MCLE method is applied to sequence data on the ANGPTL4 gene in 1832 African American and 1045 European American individuals.

Significantly fewer protein functional changing variants for lipid metabolism in Africans than in Europeans.

BACKGROUND: The disorders in metabolism of energy substances are usually related to some diseases, such as obesity, diabetes and cancer, etc. However, the genetic background for these disorders has not been well understood. In this study, we explored the genetic risk differences among human populations in metabolism (catabolism and biosynthesis) of energy substances, including lipids, carbohydrates and amino acids. RESULTS: Two genotype datasets (Hapmap and 1000 Genome) were used for this study. The genetic risks of protein functional changing variants (PFCVs) on genes involved in lipid, carbohydrate and amino acid metabolism were calculated using two genetic risk indices: the total number of PFCVs (Num) and the total possibly harmful score of PFCVs (R). Observations in these two genotype datasets consistently showed that Africans had lower genetic risk in lipid metabolism (both catabolic and biosynthetic processes) compared to Europeans. However this relationship was not observed in carbohydrate and amino acid metabolism. CONCLUSIONS: Our results suggested that Africans had higher efficiency of utilizing lipids as energy substances than Europeans. In other words, lipids might be more preferred as energy substances in Africans than in Europeans.

Assessing periodontitis in populations: a systematic review of the validity of partial-mouth examination protocols.

OBJECTIVE: To estimate bias associated with partial-mouth periodontal examination (PMPE) protocols regarding estimates of prevalence, severity and extent of clinical attachment loss (CAL), pocket depth (PD) and gingival recession (REC). MATERIAL AND METHODS: A search was made for articles published in English, from 1946 to 2012, which compared PMPE versus full-mouth periodontal examination protocols for CAL or PD ≥ 4 mm or REC ≥3 mm thresholds. PMPE protocols were evaluated for sensitivity of estimates of periodontitis prevalence, relative biases for severity and extent estimates. RESULTS: A review of the literature identified 12 studies which reported 32 PMPE protocols. Three PMPE protocols which had sensitivities ≥85% and relative biases ≤0.05 in absolute values for severity and extent estimates were as follows: (1) half-mouth six-sites, (2) diagonal quadrants six-sites and (3) full-mouth mesiobuccal-midbuccal-distobuccal (MB-B-DB). Two other PMPE protocols (full-mouth and half-mouth mesiobuccal-midbuccal-distolingual) performed well for prevalence and severity of periodontitis; however, their performance in estimates of extent was unknown. CONCLUSIONS: Among the 32 PMPE protocols listed, the half-mouth six-sites and full-mouth MB-B-DB protocols had the highest sensitivities for prevalence estimates and lowest relative biases for severity and extent estimates.