Comparative Protein Structure Modeling Using MODELLER

Multiple Functions of Compatible Solute Ectoine and Strategies for Constructing Overproducers for Biobased Production

Ectoine and its derivative 5-hydroxyectoine are compatible solutes initially found in the hyperhalophilic bacterium Ectothiorhodospira halochloris, which inhabits the desert in Egypt. The habitat of ectoine producers implies the primary function of ectoine as a cytoprotectant against harsh conditions such as high salinity, drought, and high radiation. More extensive and in-depth studies have revealed the multiple functions of ectoine in its native producer bacterial cells and other types of cells and its biomolecular components (such as proteins and DNA) as a general protective agent. Its chemical properties as a bio-based amino acid derivative make it attractive for basic scientific research and related industries, such as the food/agricultural industry, cosmetic manufacturing, biologics, and therapeutic agent preparation. This article first discusses the functions and applications of ectoine and 5-hydroxyectoine. Subsequently, more emphasis was placed on advances in bio-based ectoine and/or 5-hydroxyectoine production. Strategies for developing more robust cell factories for highly efficient ectoine and/or 5-hydroxyectoine production are further discussed. We hope this review will provide a valuable reference for studies on the bio-based production of ectoine and 5-hydroxyectoine.

High-throughput prioritization of target proteins for development of new antileishmanial compounds

Leishmaniasis, a vector-borne disease, is caused by the infection of Leishmania spp., obligate intracellular protozoan parasites. Presently, human vaccines are unavailable, and the primary treatment relies heavily on systemic drugs, often presenting with suboptimal formulations and substantial toxicity, making new drugs a high priority for LMIC countries burdened by the disease, but a low priority in the agenda of most pharmaceutical companies due to unattractive profit margins. New ways to accelerate the discovery of new, or the repositioning of existing drugs, are needed. To address this challenge, our study aimed to identify potential protein targets shared among clinically-relevant Leishmania species. We employed a subtractive proteomics and comparative genomics approach, integrating high-throughput multi-omics data to classify these targets based on different druggability metrics. This effort resulted in the ranking of 6502 ortholog groups of protein targets across 14 pathogenic Leishmania species. Among the top 20 highly ranked groups, metabolic processes known to be attractive drug targets, including the ubiquitination pathway, aminoacyl-tRNA synthetases, and purine synthesis, were rediscovered. Additionally, we unveiled novel promising targets such as the nicotinate phosphoribosyltransferase enzyme and dihydrolipoamide succinyltransferases. These groups exhibited appealing druggability features, including less than 40% sequence identity to the human host proteome, predicted essentiality, structural classification as highly druggable or druggable, and expression levels above the 50th percentile in the amastigote form. The resources presented in this work also represent a comprehensive collection of integrated data regarding trypanosomatid biology.

Deafness DFNB128 Associated with a Recessive Variant of Human MAP3K1 Recapitulates Hearing Loss of Map3k1-Deficient Mice

Deafness in vertebrates is associated with variants of hundreds of genes. Yet, many mutant genes causing rare forms of deafness remain to be discovered. A consanguineous Pakistani family segregating nonsyndromic deafness in two sibships were studied using microarrays and exome sequencing. A 1.2 Mb locus (DFNB128) on chromosome 5q11.2 encompassing six genes was identified. In one of the two sibships of this family, a novel homozygous recessive variant NM_005921.2:c.4460G>A p.(Arg1487His) in the kinase domain of MAP3K1 co-segregated with nonsyndromic deafness. There are two previously reported Map3k1-kinase-deficient mouse models that are associated with recessively inherited syndromic deafness. MAP3K1 phosphorylates serine and threonine and functions in a signaling pathway where pathogenic variants of HGF, MET, and GAB1 were previously reported to be associated with human deafness DFNB39, DFNB97, and DFNB26, respectively. Our single-cell transcriptome data of mouse cochlea mRNA show expression of Map3k1 and its signaling partners in several inner ear cell types suggesting a requirement of wild-type MAP3K1 for normal hearing. In contrast to dominant variants of MAP3K1 associated with Disorders of Sex Development 46,XY sex-reversal, our computational modeling of the recessive substitution p.(Arg1487His) predicts a subtle structural alteration in MAP3K1, consistent with the limited phenotype of nonsyndromic deafness.

Exploring the Activation Mechanism of the GPR183 Receptor

The G protein-coupled receptors (GPCRs) play a pivotal role in numerous biological processes as crucial cell membrane receptors. However, the dynamic mechanisms underlying the activation of GPR183, a specific GPCR, remain largely elusive. To address this, we employed computational simulation techniques to elucidate the activation process and key events associated with GPR183, including conformational changes from inactive to active state, binding interactions with the Gi protein complex, and GDP release. Our findings demonstrate that the association between GPR183 and the Gi protein involves the formation of receptor-specific conformations, the gradual proximity of the Gi protein to the binding pocket, and fine adjustments of the protein conformation, ultimately leading to a stable GPR183-Gi complex characterized by a high energy barrier. The presence of Gi protein partially promotes GPR183 activation, which is consistent with the observation of GPCR constitutive activity test experiments, thus illustrating the reliability of our calculations. Moreover, our study suggests the existence of a stable partially activated state preceding complete activation, providing novel avenues for future investigations. In addition, the relevance of GPR183 for various diseases, such as colitis, the response of eosinophils to Mycobacterium tuberculosis infection, antiviral properties, and pulmonary inflammation, has been emphasized, underscoring its therapeutic potential. Consequently, understanding the activation process of GPR183 through molecular dynamic simulations offers valuable kinetic insights that can aid in the development of targeted therapies.

Binding characteristics of the doxepin E/Z-isomers to the histamine H1 receptor revealed by receptor-bound ligand analysis and molecular dynamics study

Doxepin is an antihistamine and tricyclic antidepressant that binds to the histamine H1 receptor (H1R) with high affinity. Doxepin is an 85:15 mixture of the E- and Z-isomers. The Z-isomer is well known to be more effective than the E-isomer, whereas based on the crystal structure of the H1R/doxepin complex, the hydroxyl group of Thr1123.37 is close enough to form a hydrogen bond with the oxygen atom of the E-isomer. The detailed binding characteristics and reasons for the differences remain unclear. In this study, we analyzed doxepin isomers bound to the receptor following extraction from a purified H1R protein complexed with doxepin. The ratio of the E- and Z-isomers bound to wild-type (WT) H1R was 55:45, indicating that the Z-isomer was bound to WT H1R with an approximately 5.2-fold higher affinity than the E-isomer. For the T1123.37V mutant, the E/Z ratio was 89:11, indicating that both isomers have similar affinities. Free energy calculations using molecular dynamics (MD) simulations also reproduced the experimental results of the relative binding free energy differences between the isomers for WT and T1123.37V. Furthermore, MD simulations revealed that the hydroxyl group of T1123.37 did not form hydrogen bonds with the E-isomer, but with the adjacent residues in the binding pocket. Analysis of the receptor-bound doxepin and MD simulations suggested that the hydroxyl group of T1123.37 contributes to the formation of a chemical environment in the binding pocket, which is slightly more favorable for the Z-isomer without hydrogen bonding with doxepin.

Molecular dynamics simulations of lipid-protein interactions in SLC4 proteins

The SLC4 family of secondary bicarbonate transporters is responsible for the transport of HCO3-, CO32-, Cl-, Na+, K+, NH3, and H+, which are necessary for regulation of pH and ion homeostasis. They are widely expressed in numerous tissues throughout the body and function in different cell types with different membrane properties. Potential lipid roles in SLC4 function have been reported in experimental studies, focusing mostly on two members of the family: AE1 (Cl-/HCO3- exchanger) and NBCe1 (Na+-CO32-cotransporter). Previous computational studies of the outward-facing state of AE1 with model lipid membranes revealed enhanced protein-lipid interactions between cholesterol (CHOL) and phosphatidylinositol bisphosphate (PIP2). However, the protein-lipid interactions in other members of the family and other conformation states are still poorly understood and this precludes the detailed studies of a potential regulatory role for lipids in the SLC4 family. In this work, we performed coarse-grained and atomistic molecular dynamics simulations on three members of the SLC4 family with different transport modes: AE1, NBCe1, and NDCBE (an Na+-CO32-/Cl- exchanger), in model HEK293 membranes consisting of CHOL, PIP2, phosphatidylcholine, phosphatidylethanolamine, phosphatidylserine, and sphingomyelin. The recently resolved inward-facing state of AE1 was also included in the simulations. Lipid-protein contact analysis of the simulated trajectories was performed with the ProLint server, which provides a multitude of visualization tools for illustration of areas of enhanced lipid-protein contact and identification of putative lipid binding sites within the protein matrix. We observed enrichment of CHOL and PIP2 around all proteins with subtle differences in their distribution depending on the protein type and conformation state. Putative binding sites were identified for CHOL, PIP2, phosphatidylcholine, and sphingomyelin in the three studied proteins, and their potential roles in the SLC4 transport function, conformational transition, and protein dimerization are discussed.

Exploring the structural assembly of rice ADP-glucose pyrophosphorylase subunits using MD simulation

ADP-glucose pyrophosphorylase plays a pivotal role as an allosteric enzyme, essential for starch biosynthesis in plants. The higher plant AGPase comparises of a pair of large and a pair of small subunits to form a heterotetrameric complex. Growing evidence indicates that each subunit plays a distinct role in regulating the underlying mechanism of starch biosynthesis. In the rice genome, there are four large subunit genes (OsL1-L4) and three small subunit genes (OsS1, OsS2a, and OsS2b). While the structural assembly of cytosolic rice AGPase subunits (OsL2:OsS2b) has been elucidated, there is currently no such documented research available for plastidial rice AGPases (OsL1:OsS1). In this study, we employed protein modeling and MD simulation approaches to gain insights into the structural association of plastidial rice AGPase subunits. Our results demonstrate that the heterotetrameric association of OsL1:OsS1 is very similar to that of cytosolic OsL2:OsS2b and potato AGPase heterotetramer (StLS:StSS). Moreover, the yeast-two-hybrid results on OsL1:OsS1, which resemble StLS:StSS, suggest a differential protein assembly for OsL2:OsS2b. Thus, the regulatory and catalytic mechanisms for plastidial AGPases (OsL1:OsS1) could be different in rice culm and developing endosperm compared to those of OsL2:OsS2b, which are predominantly found in rice endosperm.

Recombinant cellobiose dehydrogenase from Thermothelomyces thermophilus: Its functional characterization and applicability in cellobionic acid production

The fungus Thermothelomyces thermophilus is a thermotolerant microorganism that has been explored as a reservoir for enzymes (hydrolytic enzymes and oxidoreductases). The functional analysis of a recombinant cellobiose dehydrogenase (MtCDHB) from T. thermophilus demonstrated a thermophilic behavior, an optimal pH in alkaline conditions for inter-domain electron transfer, and catalytic activity on cellooligosaccharides with different degree of polymerization. Its applicability was evaluated to the sustainable production of cellobionic acid (CBA), a potential pharmaceutical and cosmetic ingredient rarely commercialized. Dissolving pulp was used as a disaccharide source for MtCDHB. Initially, recombinant exoglucanases (MtCBHI and MtCBHII) from T. thermophilus hydrolyzed the dissolving pulp, resulting in 87% cellobiose yield, which was subsequently converted into CBA by MtCDHB, achieving a 66% CBA yield after 24 h. These findings highlight the potential of MtCDHB as a novel approach to obtaining CBA through the bioconversion of a plant-based source.

Mutations of voltage-gated sodium channel contribute to pyrethroid resistance in Panonychus citri

Insecticide resistance in Panonychus citri is a major obstacle to mite control in citrus orchards. Pyrethroid insecticides are continually used to control mites in China, although resistance to pyrethroids has evolved in some populations. Here, the resistance to the pyrethroid fenpropathrin was investigated and 7 out of 8 field-collected populations of P. citri exhibited a high level of resistance, ranging from 171-fold to 15 391-fold higher than the susceptible (SS) comparison strain. Three voltage-gated sodium channel (VGSC) mutations were identified in the tested populations: L1031V, F1747L, and F1751I. Amplicon sequencing was used to evaluate the frequency of these mutations in the 19 field populations. L1031V and F1751I were present in all populations at frequencies of 11.6%-82.1% and 0.5%-31.8%, respectively, whereas the F1747L mutation was only present in 12 populations from Chongqing, Sichuan, Guangxi, and Yunnan provinces. Introduction of these mutations singly or in combination into transgenic flies significantly increased their resistance to fenpropathrin and these flies also exhibited reduced mortality after exposure to the pyrethroids permethrin and β-cypermethrin. Panonychus citri VGSC homology modeling and ligand docking indicate that F1747 and F1751 form direct binding contacts with pyrethroids, which are lost with mutation, whereas L1031 mutation may diminish pyrethroid effects through an allosteric mechanism. Overall, the results provide molecular markers for monitoring pest resistance to pyrethroids and offer new insights into the basis of pyrethroid actions on sodium channels.

High-throughput system for the thermostability analysis of proteins

Thermal stability of proteins is a primary metric for evaluating their physical properties. Although researchers attempted to predict it using machine learning frameworks, their performance has been dependent on the quality and quantity of published data. This is due to the technical limitation that thermodynamic characterization of protein denaturation by fluorescence or calorimetry in a high-throughput manner has been challenging. Obtaining a melting curve that derives solely from the target protein requires laborious purification, making it far from practical to prepare a hundred or more samples in a single workflow. Here, we aimed to overcome this throughput limitation by leveraging the high protein secretion efficacy of Brevibacillus and consecutive treatment with plate-scale purification methodologies. By handling the entire process of expression, purification, and analysis on a per-plate basis, we enabled the direct observation of protein denaturation in 384 samples within 4 days. To demonstrate a practical application of the system, we conducted a comprehensive analysis of 186 single mutants of a single-chain variable fragment of nivolumab, harvesting the melting temperature (Tm) ranging from -9.3 up to +10.8°C compared to the wild-type sequence. Our findings will allow for data-driven stabilization in protein design and streamlining the rational approaches.

L-serine treatment in patients with GRIN-related encephalopathy: a phase 2A, non-randomized study

GRIN-related disorders are rare developmental encephalopathies with variable manifestations and limited therapeutic options. Here, we present the first non-randomized, open-label, single-arm trial (NCT04646447) designed to evaluate the tolerability and efficacy of L-serine in children with GRIN genetic variants leading to loss-of-function. In this phase 2A trial, patients aged 2-18 years with GRIN loss-of-function pathogenic variants received L-serine for 52 weeks. Primary end points included safety and efficacy by measuring changes in the Vineland Adaptive Behavior Scales, Bayley Scales, age-appropriate Wechsler Scales, Gross Motor Function-88, Sleep Disturbance Scale for Children, Pediatric Quality of Life Inventory, Child Behavior Checklist and the Caregiver-Teacher Report Form following 12 months of treatment. Secondary outcomes included seizure frequency and intensity reduction and EEG improvement. Assessments were performed 3 months and 1 day before starting treatment and 1, 3, 6 and 12 months after beginning the supplement. Twenty-four participants were enrolled (13 males/11 females, mean age 9.8 years, SD 4.8), 23 of whom completed the study. Patients had GRIN2B, GRIN1 and GRIN2A variants (12, 6 and 5 cases, respectively). Their clinical phenotypes showed 91% had intellectual disability (61% severe), 83% had behavioural problems, 78% had movement disorders and 58% had epilepsy. Based on the Vineland Adaptive Behavior Composite standard scores, nine children were classified as mildly impaired (cut-off score > 55), whereas 14 were assigned to the clinically severe group. An improvement was detected in the Daily Living Skills domain (P = 0035) from the Vineland Scales within the mild group. Expressive (P = 0.005), Personal (P = 0.003), Community (P = 0.009), Interpersonal (P = 0.005) and Fine Motor (P = 0.031) subdomains improved for the whole cohort, although improvement was mostly found in the mild group. The Growth Scale Values in the Cognitive subdomain of the Bayley-III Scale showed a significant improvement in the severe group (P = 0.016), with a mean increase of 21.6 points. L-serine treatment was associated with significant improvement in the median Gross Motor Function-88 total score (P = 0.002) and the mean Pediatric Quality of Life total score (P = 0.00068), regardless of severity. L-serine normalized the EEG pattern in five children and the frequency of seizures in one clinically affected child. One patient discontinued treatment due to irritability and insomnia. The trial provides evidence that L-serine is a safe treatment for children with GRIN loss-of-function variants, having the potential to improve adaptive behaviour, motor function and quality of life, with a better response to the treatment in mild phenotypes.

Joint genotypic and phenotypic outcome modeling improves base editing variant effect quantification

CRISPR base editing screens enable analysis of disease-associated variants at scale; however, variable efficiency and precision confounds the assessment of variant-induced phenotypes. Here, we provide an integrated experimental and computational pipeline that improves estimation of variant effects in base editing screens. We use a reporter construct to measure guide RNA (gRNA) editing outcomes alongside their phenotypic consequences and introduce base editor screen analysis with activity normalization (BEAN), a Bayesian network that uses per-guide editing outcomes provided by the reporter and target site chromatin accessibility to estimate variant impacts. BEAN outperforms existing tools in variant effect quantification. We use BEAN to pinpoint common regulatory variants that alter low-density lipoprotein (LDL) uptake, implicating previously unreported genes. Additionally, through saturation base editing of LDLR, we accurately quantify missense variant pathogenicity that is consistent with measurements in UK Biobank patients and identify underlying structural mechanisms. This work provides a widely applicable approach to improve the power of base editing screens for disease-associated variant characterization.

Computation-Guided Rational Design of Cysteine-Less Protein Variants in Engineered hCGL

The engineered human cystathionine-γ-lyase (hCGL) resulting in enhanced activity toward both cysteine and cystine unveils a potential robust antitumor activity. However, the presence of cysteine residues has the potential to induce oligomerization or incorrect disulfide bonding, which may decrease the bioavailability of biopharmaceuticals. Through a meticulous design process targeting the cysteine residues within engineered hCGL, a set of potential beneficial mutants were obtained by virtual screening employing Rosetta and ABACUS. Experimental measurements have revealed that most of the mutants showed increased activity toward both substrates l-Cys and CSSC. Furthermore, mutants C109V and C229D demonstrated Tm value increases of 8.2 and 1.8 °C, respectively. After an 80 min incubation at 60 °C, mutant C229D still maintained high residual activity. Unexpectedly, mutant C109V, displaying activity approximately 2-fold higher than the activity of wild type (WT) for both substrates, showed disappointing instability in plasma, which suggests that computational design still requires further consideration. Analysis of their structure and molecular dynamics (MD) simulation revealed the impact of hydrophobic interaction, hydrogen bonds, and near-attack conformation (NAC) stability on activity and stability. This study acquired information about mutants that exhibit enhanced activity or thermal resistance and serve as valuable guidance for subsequent specific cysteine modifications.

Uncovering the molecular interactions underlying MBD2 and MBD3 phase separation

Chromatin organization controls DNA's accessibility to regulatory factors to influence gene expression. Heterochromatin, or transcriptionally silent chromatin enriched in methylated DNA and methylated histone tails, self-assembles through multivalent interactions with its associated proteins into a condensed, but dynamic state. Liquid-liquid phase separation (LLPS) of key heterochromatin regulators, such as heterochromatin protein 1 (HP1), plays an essential role in heterochromatin assembly and function. Methyl-CpG-binding protein 2 (MeCP2), the most studied member of the methyl-CpG-binding domain (MBD) family of proteins, has been recently shown to undergo LLPS in the absence and presence of methylated DNA. These studies provide a new mechanistic framework for understanding the role of methylated DNA and its readers in heterochromatin formation. However, the details of the molecular interactions by which other MBD family members undergo LLPS to mediate genome organization and transcriptional regulation are not fully understood. Here, we focus on two MBD proteins, MBD2 and MBD3, that have distinct but interdependent roles in gene regulation. Using an integrated computational and experimental approach, we uncover the homotypic and heterotypic interactions governing MBD2 and MBD3 phase separation and DNA's influence on this process. We show that despite sharing the highest sequence identity and structural homology among all the MBD protein family members, MBD2 and MBD3 exhibit differing residue patterns resulting in distinct phase separation mechanisms. Understanding the molecular underpinnings of MBD protein condensation offers insights into the higher-order, LLPS-mediated organization of heterochromatin.

Potential Involvement of the South American Lungfish Intelectin-2 in Innate-Associated Immune Modulation

Intelectins belong to a family of lectins with specific and transitory carbohydrate interaction capabilities. These interactions are related to the activity of agglutinating pathogens, as intelectins play a significant role in immunity. Despite the prominent immune defense function of intelectins, limited information about its structural characteristics and carbohydrate interaction properties is available. This study investigated an intelectin transcript identified in RNA-seq data obtained from the South American lungfish (Lepidosiren paradoxa), namely LpITLN2-B. The structural analyses predicted LpITLN2-B to be a homo-trimeric globular protein with the fibrinogen-like functional domain (FReD), exhibiting a molecular mass of 57 kDa. The quaternary structure is subdivided into three monomers, A, B, and C, and each domain comprises 11 β-sheets: an anti-parallel β-sheet, a β-hairpin, and a disordered β-sheet structure. Molecular docking demonstrates a significant interaction with disaccharides rather than monosaccharides. The preferential interaction with disaccharides highlights the potential interaction with pathogen molecules, such as LPS and Poly(I:C). The hemagglutination assay inhibited lectins activity, especially maltose and sucrose, highlighting lectin activity in L. paradoxa samples. Overall, our results show the potential relevance of LpITLN2-B in L. paradoxa immune defense against pathogens.

Comprehensive Characterization of LAT1 Cholesterol-Binding Sites

The human L-type amino acid transporter 1 (LAT1; SLC7A5), is an amino acid exchanger protein, primarily found in the blood-brain barrier, placenta, and testis, where it plays a key role in amino acid homeostasis. Cholesterol is an essential lipid that has been highlighted to play a role in regulating the activity of membrane transporters, such as LAT1, yet little is known about the molecular mechanisms driving this phenomenon. Here we perform a comprehensive computational analysis to investigate cholesterol's role in LAT1 structure and function, focusing on four cholesterol-binding sites (CHOL1-4) identified in a recent LAT1-apo inward-open conformation cryo-EM structure. Through a series of independent molecular dynamics (MD) simulations, molecular docking, MM/GBSA free energy calculations, and other analysis tools, we explored the interactions between LAT1 and cholesterol. Our findings suggest that CHOL3 forms the most stable and favorable interactions with LAT1. Principal component analysis (PCA) and center of mass (COM) distance assessments show that CHOL3 binding stabilizes the inward-open state of LAT1 by preserving the spatial arrangement of the hash and bundle domains. Additionally, we propose an alternative cholesterol-binding site for originally assigned CHOL1. Overall, this study improves the understanding of cholesterol's modulatory effect on LAT1 and proposes candidate sites for the discovery of future allosteric ligands with rational design.

Structural basis for antiepileptic drugs and botulinum neurotoxin recognition of SV2A

More than one percent of people have epilepsy worldwide. Levetiracetam (LEV) is a successful new-generation antiepileptic drug (AED), and its derivative, brivaracetam (BRV), shows improved efficacy. Synaptic vesicle glycoprotein 2a (SV2A), a putative membrane transporter in the synaptic vesicles (SVs), has been identified as a target of LEV and BRV. SV2A also serves as a receptor for botulinum neurotoxin (BoNT), which is the most toxic protein and has paradoxically emerged as a potent reagent for therapeutic and cosmetic applications. Nevertheless, no structural analysis on AEDs and BoNT recognition by full-length SV2A has been available. Here we describe the cryo-electron microscopy structures of the full-length SV2A in complex with the BoNT receptor-binding domain, BoNT/A2 HC, and either LEV or BRV. The large fourth luminal domain of SV2A binds to BoNT/A2 HC through protein-protein and protein-glycan interactions. LEV and BRV occupy the putative substrate-binding site in an outward-open conformation. A propyl group in BRV creates additional contacts with SV2A, explaining its higher binding affinity than that of LEV, which was further supported by label-free spectral shift assay. Numerous LEV derivatives have been developed as AEDs and positron emission tomography (PET) tracers for neuroimaging. Our work provides a structural framework for AEDs and BoNT recognition of SV2A and a blueprint for the rational design of additional AEDs and PET tracers.

Lethal perturbation of an Escherichia coli regulatory network is triggered by a restriction-modification system's regulator and can be mitigated by excision of the cryptic prophage Rac

Bacterial gene regulatory networks orchestrate responses to environmental challenges. Horizontal gene transfer can bring in genes with regulatory potential, such as new transcription factors (TFs), and this can disrupt existing networks. Serious regulatory perturbations may even result in cell death. Here, we show the impact on Escherichia coli of importing a promiscuous TF that has adventitious transcriptional effects within the cryptic Rac prophage. A cascade of regulatory network perturbations occurred on a global level. The TF, a C regulatory protein, normally controls a Type II restriction-modification system, but in E. coli K-12 interferes with expression of the RacR repressor gene, resulting in de-repression of the normally-silent Rac ydaT gene. YdaT is a prophage-encoded TF with pleiotropic effects on E. coli physiology. In turn, YdaT alters expression of a variety of bacterial regulons normally controlled by the RcsA TF, resulting in deficient lipopolysaccharide biosynthesis and cell division. At the same time, insufficient RacR repressor results in Rac DNA excision, halting Rac gene expression due to loss of the replication-defective Rac prophage. Overall, Rac induction appears to counteract the lethal toxicity of YdaT. We show here that E. coli rewires its regulatory network, so as to minimize the adverse regulatory effects of the imported C TF. This complex set of interactions may reflect the ability of bacteria to protect themselves by having robust mechanisms to maintain their regulatory networks, and/or suggest that regulatory C proteins from mobile operons are under selection to manipulate their host's regulatory networks for their own benefit.

Computational Mutagenesis of Antibody Fragments: Disentangling Side Chains from ΔΔG Predictions

The development of highly potent antibodies and antibody fragments as binding agents holds significant implications in fields such as biosensing and biotherapeutics. Their binding strength is intricately linked to the arrangement and composition of residues at the binding interface. Computational techniques offer a robust means to predict the three-dimensional structure of these complexes and to assess the affinity changes resulting from mutations. Given the interdependence of structure and affinity prediction, our objective here is to disentangle their roles. We aim to evaluate independently six side-chain reconstruction methods and ten binding affinity estimation techniques. This evaluation was pivotal in predicting affinity alterations due to single mutations, a key step in computational affinity maturation protocols. Our analysis focuses on a data set comprising 27 distinct antibody/hen egg white lysozyme complexes, each with crystal structures and experimentally determined binding affinities. Using six different side-chain reconstruction methods, we transformed each structure into its corresponding mutant via in silico single-point mutations. Subsequently, these structures undergo minimization and molecular dynamics simulation. We therefore estimate ΔΔG values based on the original crystal structure, its energy-minimized form, and the ensuing molecular dynamics trajectories. Our research underscores the critical importance of selecting reliable side-chain reconstruction methods and conducting thorough molecular dynamics simulations to accurately predict the impact of mutations. In summary, our study demonstrates that the integration of conformational sampling and scoring is a potent approach to precisely characterizing mutation processes in single-point mutagenesis protocols and crucial for computational antibody design.

ZP2 cleavage blocks polyspermy by modulating the architecture of the egg coat

Following the fertilization of an egg by a single sperm, the egg coat or zona pellucida (ZP) hardens and polyspermy is irreversibly blocked. These events are associated with the cleavage of the N-terminal region (NTR) of glycoprotein ZP2, a major subunit of ZP filaments. ZP2 processing is thought to inactivate sperm binding to the ZP, but its molecular consequences and connection with ZP hardening are unknown. Biochemical and structural studies show that cleavage of ZP2 triggers its oligomerization. Moreover, the structure of a native vertebrate egg coat filament, combined with AlphaFold predictions of human ZP polymers, reveals that two protofilaments consisting of type I (ZP3) and type II (ZP1/ZP2/ZP4) components interlock into a left-handed double helix from which the NTRs of type II subunits protrude. Together, these data suggest that oligomerization of cleaved ZP2 NTRs extensively cross-links ZP filaments, rigidifying the egg coat and making it physically impenetrable to sperm.

Phage proteins target and co-opt host ribosomes immediately upon infection

Bacteriophages must seize control of the host gene expression machinery to replicate. To bypass bacterial anti-phage defence systems, this host takeover occurs immediately upon infection. A general understanding of phage mechanisms for immediate targeting of host transcription and translation processes is lacking. Here we introduce an integrative high-throughput approach to uncover phage-encoded proteins that target the gene expression machinery of Pseudomonas aeruginosa immediately upon infection with the jumbo phage ΦKZ. By integrating biochemical, genetic and structural analyses, we identify an abundant and conserved phage factor ΦKZ014 that targets the large ribosomal subunit by binding the 5S ribosomal RNA, and rapidly promotes replication in several clinical isolates. ΦKZ014 is among the earliest ΦKZ proteins expressed after infection and remains bound to ribosomes during the entire translation cycle. Our study provides a strategy to decipher molecular components of phage-mediated host takeover and argues that phage genomes represent an untapped discovery space for proteins that modulate the host gene expression machinery.

Predicting Mutational Effects on Ca2+-Activated Chloride Conduction of TMEM16A Based on a Simulation Study

The dysfunction and defects of ion channels are associated with many human diseases, especially for loss-of-function mutations in ion channels such as cystic fibrosis transmembrane conductance regulator mutations in cystic fibrosis. Understanding ion channels is of great current importance for both medical and fundamental purposes. Such an understanding should include the ability to predict mutational effects and describe functional and mechanistic effects. In this work, we introduce an approach to predict mutational effects based on kinetic information (including reaction barriers and transition state locations) obtained by studying the working mechanism of target proteins. Specifically, we take the Ca2+-activated chloride channel TMEM16A as an example and utilize the computational biology model to predict the mutational effects of key residues. Encouragingly, we verified our predictions through electrophysiological experiments, demonstrating a 94% prediction accuracy regarding mutational directions. The mutational strength assessed by Pearson's correlation coefficient is -0.80 between our calculations and the experimental results. These findings suggest that the proposed methodology is reliable and can provide valuable guidance for revealing functional mechanisms and identifying key residues of the TMEM16A channel. The proposed approach can be extended to a broad scope of biophysical systems.

Mobile barrier mechanisms for Na+-coupled symport in an MFS sugar transporter

While many 3D structures of cation-coupled transporters have been determined, the mechanistic details governing the obligatory coupling and functional regulations still remain elusive. The bacterial melibiose transporter (MelB) is a prototype of major facilitator superfamily transporters. With a conformation-selective nanobody, we determined a low-sugar affinity inward-facing Na+-bound cryoEM structure. The available outward-facing sugar-bound structures showed that the N- and C-terminal residues of the inner barrier contribute to the sugar selectivity. The inward-open conformation shows that the sugar selectivity pocket is also broken when the inner barrier is broken. Isothermal titration calorimetry measurements revealed that this inward-facing conformation trapped by this nanobody exhibited a greatly decreased sugar-binding affinity, suggesting the mechanisms for substrate intracellular release and accumulation. While the inner/outer barrier shift directly regulates the sugar-binding affinity, it has little or no effect on the cation binding, which is supported by molecular dynamics simulations. Furthermore, the hydron/deuterium exchange mass spectrometry analyses allowed us to identify dynamic regions; some regions are involved in the functionally important inner barrier-specific salt-bridge network, which indicates their critical roles in the barrier switching mechanisms for transport. These complementary results provided structural and dynamic insights into the mobile barrier mechanism for cation-coupled symport.

Altered Protein Dynamics and a More Reactive Catalytic Cysteine in a Neurodegeneration-associated UCHL1 Mutant

A mutant of ubiquitin C-terminal hydrolase L1 (UCHL1) detected in early-onset neurodegenerative patients, UCHL1R178Q, showed higher catalytic activity than wild-type UCHL1 (UCHL1WT). Lying within the active-site pocket, the arginine is part of an interaction network that holds the catalytic histidine in an inactive arrangement. However, the structural basis and mechanism of enzymatic activation upon glutamine substitution was not understood. We combined X-ray crystallography, protein nuclear magnetic resonance (NMR) analysis, enzyme kinetics, covalent inhibition analysis, and biophysical measurements to delineate activating factors in the mutant. While the crystal structure of UCHL1R178Q showed nearly the same arrangement of the catalytic residues and active-site pocket, the mutation caused extensive alteration in the chemical environment and dynamics of more than 30 residues, some as far as 15 Å away from the site of mutation. Significant broadening of backbone amide resonances in the HSQC spectra indicates considerable backbone dynamics changes in several residues, in agreement with solution small-angle X-ray scattering (SAXS) analyses which indicate an overall increase in protein flexibility. Enzyme kinetics show the activation is due to a kcat effect despite a slightly weakened substrate affinity. In line with this, the mutant shows a higher second-order rate constant (kinact/Ki) in a reaction with a substrate-derived irreversible inhibitor, Ub-VME, compared to the wild-type enzyme, an observation indicative of a more reactive catalytic cysteine in the mutant. Together, the observations underscore structural plasticity as a factor contributing to enzyme kinetic behavior which can be modulated through mutational effects.

Identification of Residues Potentially Involved in Optical Shifts in the Water-Soluble Chlorophyll a-Binding Protein through Molecular Dynamics Simulations

Reversible light and thermally induced spectral shifts are universally observed in a wide variety of pigment-protein complexes at temperatures ranging from cryogenic to ambient. In this paper, we employed large-scale molecular dynamics (MD) simulations of a prototypical pigment-protein complex to better understand these shifts at a molecular scale. Although multiple mechanisms have been proposed over the years, no verification of these proposals via MD simulations has thus far been performed; our work represents the first step in this direction. From simulations of the water-soluble chlorophyll-binding protein complex, we determined that rearrangements of long hydrogen bonds were unlikely to be the origin of the multiwell landscape features necessary to explain observed spectral shifts. We also assessed small motions of amino acid residues and identified side chain rotations of some of these residues as likely candidates for the origin of relevant multiwell landscape features. The protein free-energy landscapes associated with side chain rotations feature energy barriers of around 1100-1600 cm-1, in agreement with optical spectroscopy results, with the most promising residue type associated with experimental signatures being serine, which possesses a symmetric triple-well landscape and moment of inertia of a relevant magnitude.

DockOpt: A Tool for Automatic Optimization of Docking Models

Molecular docking is a widely used technique for leveraging protein structure for ligand discovery, but it remains difficult to utilize due to limitations that have not been adequately addressed. Despite some progress toward automation, docking still requires expert guidance, hindering its adoption by a broader range of investigators. To make docking more accessible, we developed a new utility called DockOpt, which automates the creation, evaluation, and optimization of docking models prior to their deployment in large-scale prospective screens. DockOpt outperforms our previous automated pipeline across all 43 targets in the DUDE-Z benchmark data set, and the generated models for 84% of targets demonstrate sufficient enrichment to warrant their use in prospective screens, with normalized LogAUC values of at least 15%. DockOpt is available as part of the Python package Pydock3 included in the UCSF DOCK 3.8 distribution, which is available for free to academic researchers at https://dock.compbio.ucsf.edu and free for everyone upon registration at https://tldr.docking.org.

Surface Charge Can Modulate Phase Separation of Multidomain Proteins

Phase separation has emerged as an important mechanism explaining the formation of certain biomolecular condensates. Biological phase separation is often driven by the multivalent interactions of modular protein domains. Beyond valency, the physical features of folded domains that promote phase separation are poorly understood. We used a model system─the small ubiquitin modifier (SUMO) and its peptide ligand, the SUMO interaction motif (SIM)─to examine how domain surface charge influences multivalency-driven phase separation. Phase separation of polySUMO and polySIM was altered by pH via a change in the protonation state of SUMO surface histidines. These effects were recapitulated by histidine mutations, which modulated SUMO solubility and polySUMO-polySIM phase separation in parallel and were quantitatively explained by atomistic modeling of weak interactions among proteins in the system. Thus, surface charge can tune the phase separation of multivalent proteins, suggesting a means of controlling phase separation biologically, evolutionarily, and therapeutically.

Computational Design of Peptide Assemblies

With the ongoing development of peptide self-assembling materials, there is growing interest in exploring novel functional peptide sequences. From short peptides to long polypeptides, as the functionality increases, the sequence space is also expanding exponentially. Consequently, attempting to explore all functional sequences comprehensively through experience and experiments alone has become impractical. By utilizing computational methods, especially artificial intelligence enhanced molecular dynamics (MD) simulation and de novo peptide design, there has been a significant expansion in the exploration of sequence space. Through these methods, a variety of supramolecular functional materials, including fibers, two-dimensional arrays, nanocages, etc., have been designed by meticulously controlling the inter- and intramolecular interactions. In this review, we first provide a brief overview of the current main computational methods and then focus on the computational design methods for various self-assembled peptide materials. Additionally, we introduce some representative protein self-assemblies to offer guidance for the design of self-assembling peptides.

Identification and prioritisation of potential vaccine candidates using subtractive proteomics and designing of a multi-epitope vaccine against Wuchereria bancrofti

This study employed subtractive proteomics and immunoinformatics to analyze the Wuchereria bancrofti proteome and identify potential therapeutic targets, with a focus on designing a vaccine against the parasite species. A comprehensive bioinformatics analysis of the parasite's proteome identified 51 probable therapeutic targets, among which "Kunitz/bovine pancreatic trypsin inhibitor domain-containing protein" was identified as the most promising vaccine candidate. The candidate protein was used to design a multi-epitope vaccine, incorporating B-cell and T-cell epitopes identified through various tools. The vaccine construct underwent extensive analysis of its antigenic, physical, and chemical features, including the determination of secondary and tertiary structures. Docking and molecular dynamics simulations were performed with HLA alleles, Toll-like receptor 4 (TLR4), and TLR3 to assess its potential to elicit the human immune response. Immune simulation analysis confirmed the predicted vaccine's strong binding affinity with immunoglobulins, indicating its potential efficacy in generating an immune response. However, experimental validation and testing of this multi-epitope vaccine construct would be needed to assess its potential against W. bancrofti and even for a broader range of lymphatic filarial infections given the similarities between W. bancrofti and Brugia.

Uncovering the clinical relevance of unclassified variants in DNA repair genes: a focus on BRCA negative Tunisian cancer families

Introduction: Recent advances in sequencing technologies have significantly increased our capability to acquire large amounts of genetic data. However, the clinical relevance of the generated data continues to be challenging particularly with the identification of Variants of Uncertain Significance (VUSs) whose pathogenicity remains unclear. In the current report, we aim to evaluate the clinical relevance and the pathogenicity of VUSs in DNA repair genes among Tunisian breast cancer families. Methods: A total of 67 unsolved breast cancer cases have been investigated. The pathogenicity of VUSs identified within 26 DNA repair genes was assessed using different in silico prediction tools including SIFT, PolyPhen2, Align-GVGD and VarSEAK. Effects on the 3D structure were evaluated using the stability predictor DynaMut and molecular dynamics simulation with NAMD. Family segregation analysis was also performed. Results: Among a total of 37 VUSs identified, 11 variants are likely deleterious affecting ATM, BLM, CHEK2, ERCC3, FANCC, FANCG, MSH2, PMS2 and RAD50 genes. The BLM variant, c.3254dupT, is novel and seems to be associated with increased risk of breast, endometrial and colon cancer. Moreover, c.6115G>A in ATM and c.592+3A>T in CHEK2 were of keen interest identified in families with multiple breast cancer cases and their familial cosegregation with disease has been also confirmed. In addition, functional in silico analyses revealed that the ATM variant may lead to protein immobilization and rigidification thus decreasing its activity. We have also shown that FANCC and FANCG variants may lead to protein destabilization and alteration of the structure compactness which may affect FANCC and FANCG protein activity. Conclusion: Our findings revealed that VUSs in DNA repair genes might be associated with increased cancer risk and highlight the need for variant reclassification for better disease management. This will help to improve the genetic diagnosis and therapeutic strategies of cancer patients not only in Tunisia but also in neighboring countries.

Elucidating the binding mechanism of SARS-CoV-2 NSP6-TBK1 and structure-based designing of phytocompounds inhibitors for instigating the host immune response

SARS-CoV-2, also referred to as severe acute respiratory syndrome coronavirus 2, is the virus responsible for causing COVID-19, an infectious disease that emerged in Wuhan, China, in December 2019. Among its crucial functions, NSP6 plays a vital role in evading the human immune system by directly interacting with a receptor called TANK-binding kinase (TBK1), leading to the suppression of IFNβ production. Consequently, in the present study we used the structural and biophysical approaches to analyze the effect of newly emerged mutations on the binding of NSP6 and TBK1. Among the identified mutations, four (F35G, L37F, L125F, and I162T) were found to significantly destabilize the structure of NSP6. Furthermore, the molecular docking analysis highlighted that the mutant NSP6 displayed its highest binding affinity with TBK1, exhibiting docking scores of -1436.2 for the wildtype and -1723.2, -1788.6, -1510.2, and -1551.7 for the F35G, L37F, L125F, and I162T mutants, respectively. This suggests the potential for an enhanced immune system evasion capability of NSP6. Particularly, the F35G mutation exhibited the strongest binding affinity, supported by a calculated binding free energy of -172.19 kcal/mol. To disrupt the binding between NSP6 and TBK1, we conducted virtual drug screening to develop a novel inhibitor derived from natural products. From this screening, we identified the top 5 hit compounds as the most promising candidates with a docking score of -6.59 kcal/mol, -6.52 kcal/mol, -6.32 kcal/mol, -6.22 kcal/mol, and -6.21 kcal/mol. The molecular dynamic simulation of top 3 hits further verified the dynamic stability of drugs-NSP6 complexes. In conclusion, this study provides valuable insight into the higher infectivity of the SARS-CoV-2 new variants and a strong rationale for the development of novel drugs against NSP6.

Remote homology identification of the Drosophila melanogaster ortholog of the RNA Polymerase I subunit Rpa34/POLR1G

Highly conserved orthologous proteins are easily identified by sequence homology alone, whereas poorly conserved orthologs require additional structural information to be identified. All Drosophila orthologs of RNA polymerase I, II, and III subunits-except one-have been identified by sequence homology. Here, we identified CG11076 as the missing Rpa34/POLR1G ortholog in Drosophila . Remote homology detection and secondary structure analysis showed that CG11076 is predicted to have high structural conservation with Rpa34/POLR1G, and phylogenetic analysis demonstrated that these proteins are closely related. Our work underscores the importance of utilizing both sequence and structure to identify highly divergent orthologous proteins in different species.

Molecular insights into atypical modes of β-arrestin interaction with seven transmembrane receptors

β-arrestins (βarrs) are multifunctional proteins involved in signaling and regulation of seven transmembrane receptors (7TMRs), and their interaction is driven primarily by agonist-induced receptor activation and phosphorylation. Here, we present seven cryo-electron microscopy structures of βarrs either in the basal state, activated by the muscarinic receptor subtype 2 (M2R) through its third intracellular loop, or activated by the βarr-biased decoy D6 receptor (D6R). Combined with biochemical, cellular, and biophysical experiments, these structural snapshots allow the visualization of atypical engagement of βarrs with 7TMRs and also reveal a structural transition in the carboxyl terminus of βarr2 from a β strand to an α helix upon activation by D6R. Our study provides previously unanticipated molecular insights into the structural and functional diversity encoded in 7TMR-βarr complexes with direct implications for exploring novel therapeutic avenues.

Conformational transitions of the HIV-1 Gag polyprotein upon multimerization and gRNA binding

During the HIV-1 assembly process, the Gag polyprotein multimerizes at the producer cell plasma membrane, resulting in the formation of spherical immature virus particles. Gag-genomic RNA (gRNA) interactions play a crucial role in the multimerization process, which is yet to be fully understood. We performed large-scale all-atom molecular dynamics simulations of membrane-bound full-length Gag dimer, hexamer, and 18-mer. The inter-domain dynamic correlation of Gag, quantified by the heterogeneous elastic network model applied to the simulated trajectories, is observed to be altered by implicit gRNA binding, as well as the multimerization state of the Gag. The lateral dynamics of our simulated membrane-bound Gag proteins, with and without gRNA binding, agree with prior experimental data and help to validate our simulation models and methods. The gRNA binding is observed to affect mainly the SP1 domain of the 18-mer and the matrix-capsid linker domain of the hexamer. In the absence of gRNA binding, the independent dynamical motion of the nucleocapsid domain results in a collapsed state of the dimeric Gag. Unlike stable SP1 helices in the six-helix bundle, without IP6 binding, the SP1 domain undergoes a spontaneous helix-to-coil transition in the dimeric Gag. Together, our findings reveal conformational switches of Gag at different stages of the multimerization process and predict that the gRNA binding reinforces an efficient binding surface of Gag for multimerization, and also regulates the dynamic organization of the local membrane region itself.

Computational design of scFv anti-receptor binding domain of SARS-CoV-2 spike protein based on antibody S230 anti-SARS-CoV-1

Developing therapeutics such as neutralizing antibodies targeting the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein is essential to halt the Covid-19 infection. However, antibody production is expensive and relatively inaccessible to many low-income countries. Therefore, a more efficient and smaller antibody fragment, such as a single-chain variable fragment (scFv), derived from a known neutralizing antibody structure, is of interest due to the lower cost of recombinant protein production and the ability to tailor scFvs against circulating viruses. In this study, we used computational design to create an scFv based on the structure of a known neutralizing antibody, S230, for SARS-CoV-1. By analyzing the interaction of S230 with the RBD of both SARS-CoV-1 and SARS-CoV-2, five mutations were introduced to improve the binding of the scFv to the RBD of SARS-CoV-2. These mutations were Ser32Thr, Trp99Val, Asn57Val, Lys65Glu, and Tyr106Ile. Molecular dynamics simulations were used to evaluate the stability and affinity of the designed scFv. Our results showed that the designed scFv improved binding to the RBD of SARS-CoV-2 compared to the original S230, as indicated by principal component analysis, distance analysis, and MM/GBSA interaction energy. Furthermore, a positive result in a spot test lateral flow assay of the expressed scFv against the RBD indicated that the mutations did not alter the protein's structure. The designed scFv showed a negative result when tested against human serum albumin as a negative control, indicating reasonable specificity. We hope that this study will be useful in designing a specific and low-cost therapeutic agent, particularly during early outbreaks when information on neutralizing antibodies is limited.Communicated by Ramaswamy H. Sarma.

Unique Substrate Recognition and Sodium-Substrate Binding Stoichiometry in a Bacterial Serotonin Transporter, TuriSERT

All resolved high-resolution structures of the transporters in the neurotransmitter sodium symporter (NSS) family revealed that the NSS members share common structural and mechanistic features for substrate and ion binding and transport. However, a recently reported bacterial orthologue of the human serotonin transporter (hSERT), TuriSERT, possesses a structural characteristic specific for amino acid substrate binding but does transport a biogenic amine. The unique structural feature of TuriSERT requires a novel configuration for coordinating its substrate and ions. In the present study, we characterized TuriSERT expressed in Escherichia coli cells with a fluorescent substrate by biochemical, structural, and pharmacological approaches. Substrate transport by TuriSERT requires Na+ but not Cl-. Replacement of Asp262 by asparagine renders TuriSERT Cl--dependent. Substitutions of the corresponding Na1 residues did not alter Na+ dependence on substrate transport, whereas the mutation of a Na2 site residue led to a loss of transport activity, suggesting that Na+ binds only to the Na2 site in TuriSERT. In addition, substitutions of several residues essential for recognizing 5-hydroxytryptamine (5-HT) in hSERT had little effect on 5-HT displacement potency in transport assay for TuriSERT. In contrast, mutations of the residues that are proposed to coordinate with 5-HT in our docking model dramatically reduced 5-HT displacement. Furthermore, our results indicated that all tested antidepressants showed a weak inhibitory effect on TuriSERT. The present study demonstrated the existence of a unique substrate binding site and 1:1 stoichiometry of sodium-substrate binding in TuriSERT, a novel structural finding for the NSS transporters.

Innate extracellular Hsp70 inflammatory properties are mediated by the interaction of Siglec-E and LOX-1 receptors

Innate immune responses to cell damage-associated molecular patterns induce a controlled degree of inflammation, ideally avoiding the promotion of intense unwanted inflammatory adverse events. When released by damaged cells, Hsp70 can stimulate different responses that range from immune activation to immune suppression. The effects of Hsp70 are mediated through innate receptors expressed primarily by myeloid cells, such as dendritic cells (DCs). The regulatory innate receptors that bind to extracellular mouse Hsp70 (mHsp70) are not fully characterized, and neither are their potential interactions with activating innate receptors. Here, we describe that extracellular mHsp70 interacts with a receptor complex formed by inhibitory Siglec-E and activating LOX-1 on DCs. We also find that this interaction takes place within lipid microdomains, and Siglec-E acts as a negative regulator of LOX-1-mediated innate activation upon mHsp70 or oxidized LDL binding. Thus, HSP70 can both bind to and modulate the interaction of inhibitory and activating innate receptors on the cell surface. These findings add another dimension of regulatory mechanism to how self-molecules contribute to dampening of exacerbated inflammatory responses.

Mobile barrier mechanisms for Na+-coupled symport in an MFS sugar transporter

While many 3D structures of cation-coupled transporters have been determined, the mechanistic details governing the obligatory coupling and functional regulations still remain elusive. The bacterial melibiose transporter (MelB) is a prototype of the Na+-coupled major facilitator superfamily transporters. With a conformational nanobody (Nb), we determined a low-sugar affinity inward-facing Na+-bound cryoEM structure. Collectively with the available outward-facing sugar-bound structures, both the outer and inner barriers were localized. The N- and C-terminal residues of the inner barrier contribute to the sugar selectivity pocket. When the inner barrier is broken as shown in the inward-open conformation, the sugar selectivity pocket is also broken. The binding assays by isothermal titration calorimetry revealed that this inward-facing conformation trapped by the conformation-selective Nb exhibited a greatly decreased sugar-binding affinity, suggesting the mechanisms for the substrate intracellular release and accumulation. While the inner/outer barrier shift directly regulates the sugar-binding affinity, it has little or no effect on the cation binding, which is also supported by molecular dynamics simulations. Furthermore, the use of this Nb in combination with the hydron/deuterium exchange mass spectrometry allowed us to identify dynamic regions; some regions are involved in the functionally important inner barrier-specific salt-bridge network, which indicates their critical roles in the barrier switching mechanisms for transport. These complementary results provided structural and dynamic insights into the mobile barrier mechanism for cation-coupled symport.

Understanding Cysteine Reactivity in Protein Environments with Electric Fields

The role cysteine residues play in proteins is mediated by their protonation state, whereby the thiolate form of the side chain is highly reactive while the thiol form is more inert. However, the pKa of cysteine residues is hard to predict as it can differ widely from its reference value in solution, an effect that is accentuated by local effects in the heterogeneous protein environment. Here, we present a new approach to the prediction of cysteine reactivity based on electric field calculations at the thiol/thiolate group. We validated our approach by predicting the protonation state of cysteine residues in different protein environments (in the active site, at the protein surface, and buried within the protein interior), including Cys-25 in papaya protease omega, which was proven problematic for the more traditional constant pH molecular dynamics (MD) technique. We predict pKa shifts consistent with experimental observations, and the decomposition of the electric fields into contributions from molecular fragments provides a direct handle to rationalize local pH and pKa effects in proteins without introducing parameters other than those of the force field used for MD simulations.

Pyrrole-Based Enaminones as Building Blocks for the Synthesis of Indolizines and Pyrrolo[1,2-a]pyrazines Showing Potent Antifungal Activity

As a new approach, pyrrolo[1,2-a]pyrazines were synthesized through the cyclization of 2-formylpyrrole-based enaminones in the presence of ammonium acetate. The enaminones were prepared with a straightforward method, reacting the corresponding alkyl 2-(2-formyl-1H-pyrrol-1-yl)acetates, 2-(2-formyl-1H-pyrrol-1-yl)acetonitrile, and 2-(2-formyl-1H-pyrrol-1-yl)acetophenones with DMFDMA. Analogous enaminones elaborated from alkyl (E)-3-(1H-pyrrol-2-yl)acrylates were treated with a Lewis acid to afford indolizines. The antifungal activity of the series of substituted pyrroles, pyrrole-based enaminones, pyrrolo[1,2-a]pyrazines, and indolizines was evaluated on six Candida spp., including two multidrug-resistant ones. Compared to the reference drugs, most test compounds produced a more robust antifungal effect. Docking analysis suggests that the inhibition of yeast growth was probably mediated by the interaction of the compounds with the catalytic site of HMGR of the Candida species.

Structure-Based Virtual Screening of Helicobacter pylori SecA Inhibitors

The human bacterial pathogen Helicobacter pylori causes a range of gastric diseases. The killing rate of Helicobacter pylori is declining year by year because of high antibiotics resistance. It is urgent to develop new target and novel anti- Helicobacter pylori drugs. As an "energy pump" for bacterial cells, SecA is essential for bacterial growth and drives bacterial protein transmembrane transport, moreover SecA is absent in mammals, all of which nominate SecA as an attractive antimicrobial target. Here, we provided a structure-based virtual screening method to screen the 3D-diversity natural-product-like screening library against SecA for novel anti- Helicobacter pylori inhibitors with novel scaffolds. In this study, homology modeling was used to construct the three-dimensional structure of Helicobacter pylori SecA. Two rounds of molecular docking were then used to find new small-molecule inhibitors of SecA, identifying six lead candidates that maintained key interactions with the binding pocket. After that, molecular dynamics simulations were used to explore more accurate ligand-receptor binding modes in states close to natural conditions. Encouragingly, all six compounds were relatively stable during the simulation. Apart from that the binding free energy calculation based on MM/PBSA demonstrated favorable results of < -13.642 kcal/mol. Finally, ADME-T analysis indicated that these compounds were also sufficiently druggable. All six compounds can be well combined with the crystal structure, which further facilitate the development of SecA inhibitors and lead compounds against Helicobacter pylori.

Structural insights into the mechanism of GTP initiation of microtubule assembly

In eukaryotes, the dynamic assembly of microtubules (MT) plays an important role in numerous cellular processes. The underlying mechanism of GTP triggering MT assembly is still unknown. Here, we present cryo-EM structures of tubulin heterodimer at their GTP- and GDP-bound states, intermediate assembly states of GTP-tubulin, and final assembly stages of MT. Both GTP- and GDP-tubulin heterodimers adopt similar curved conformations with subtle flexibility differences. In head-to-tail oligomers of tubulin heterodimers, the inter-dimer interface of GDP-tubulin exhibits greater flexibility, particularly in tangential bending. Cryo-EM of the intermediate assembly states reveals two types of tubulin lateral contacts, "Tube-bond" and "MT-bond". Further, molecular dynamics (MD) simulations show that GTP triggers lateral contact formation in MT assembly in multiple sequential steps, gradually straightening the curved tubulin heterodimers. Therefore, we propose a flexible model of GTP-initiated MT assembly, including the formation of longitudinal and lateral contacts, to explain the nucleation and assembly of MT.

Joint genotypic and phenotypic outcome modeling improves base editing variant effect quantification

CRISPR base editing screens are powerful tools for studying disease-associated variants at scale. However, the efficiency and precision of base editing perturbations vary, confounding the assessment of variant-induced phenotypic effects. Here, we provide an integrated pipeline that improves the estimation of variant impact in base editing screens. We perform high-throughput ABE8e-SpRY base editing screens with an integrated reporter construct to measure the editing efficiency and outcomes of each gRNA alongside their phenotypic consequences. We introduce BEAN, a Bayesian network that accounts for per-guide editing outcomes and target site chromatin accessibility to estimate variant impacts. We show this pipeline attains superior performance compared to existing tools in variant classification and effect size quantification. We use BEAN to pinpoint common variants that alter LDL uptake, implicating novel genes. Additionally, through saturation base editing of LDLR, we enable accurate quantitative prediction of the effects of missense variants on LDL-C levels, which aligns with measurements in UK Biobank individuals, and identify structural mechanisms underlying variant pathogenicity. This work provides a widely applicable approach to improve the power of base editor screens for disease-associated variant characterization.

In Silico Analysis of Pyeongwi-San Involved in Inflammatory Bowel Disease Treatment Using Network Pharmacology, Molecular Docking, and Molecular Dynamics

BACKGOUND

Pyeongwi-san (PWS) is a widely used formula for treating digestive disorders in Korea and China. Inflammatory bowel disease (IBD) is characterized by progressive inflammation of the gastrointestinal tract. Emerging evidence supports the protective effect of PWS against IBD, but specific mechanisms are still elusive.

METHODS

Active compounds of PWS were screened from the medicinal materials and chemical compounds in Northeast Asian traditional medicine (TM-MC) in the consideration of drug-likeness and oral bioavailability. Target candidates of active compounds were predicted using the ChEMBL database. IBD-related targets were obtained from the GeneCards and DisGeNET databases. The network of composition-targets-disease was constructed. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were analyzed. Molecular docking was used to simulate the binding affinity of active compounds on target proteins and molecular dynamics was used to validate the molecular docking result.

RESULTS

A total of 26 core target proteins of PWS were related to IBD. Enrichment analysis suggested that PWS is highly associated with tumor necrosis factor signaling pathway, apoptosis, and the collapse of tight junctions. Moreover, molecular docking and molecular dynamics simulation proposed β-eudesmol and (3R,6R,7S)-1,10-bisaboladien-3-ol to ameliorate IBD through the binding to TNF and MMP9, respectively.

CONCLUSION

Present in silico analysis revealed potential pathways and insight of PWS to regulate IBD. These results imply that the therapeutic effect of PWS might be achieved via an inhibitory effect.

In silico identification of modulators of J domain protein-Hsp70 interactions in Plasmodium falciparum: a drug repurposing strategy against malaria

Plasmodium falciparum is a unicellular, intracellular protozoan parasite, and the causative agent of malaria in humans, a deadly vector borne infectious disease. A key phase of malaria pathology, is the invasion of human erythrocytes, resulting in drastic remodeling by exported parasite proteins, including molecular chaperones and co-chaperones. The survival of the parasite within the human host is mediated by P. falciparum heat shock protein 70s (PfHsp70s) and J domain proteins (PfJDPs), functioning as chaperones-co-chaperones partnerships. Two complexes have been shown to be important for survival and pathology of the malaria parasite: PfHsp70-x-PFE0055c (exported); and PfHsp70-2-PfSec63 (endoplasmic reticulum). Virtual screening was conducted on the drug repurposing library, the Pandemic Response Box, to identify small-molecules that could specifically disrupt these chaperone complexes. Five top ranked compounds possessing preferential binding affinity for the malarial chaperone system compared to the human system, were identified; three top PfHsp70-PfJDP binders, MBX 1641, zoliflodacin and itraconazole; and two top J domain binders, ezetimibe and a benzo-diazepinone. These compounds were validated by repeat molecular dockings and molecular dynamics simulation, resulting in all the compounds, except for MBX 1461, being confirmed to bind preferentially to the malarial chaperone system. A detailed contact analysis of the PfHsp70-PfJDP binders identified two different types of modulators, those that potentially inhibit complex formation (MBX 1461), and those that potentially stabilize the complex (zoliflodacin and itraconazole). These data suggested that zoliflodacin and itraconazole are potential novel modulators specific to the malarial system. A detailed contact analysis of the J domain binders (ezetimibe and the benzo-diazepinone), revealed that they bound with not only greater affinity but also a better pose to the malarial J domain compared to that of the human system. These data suggested that ezetimibe and the benzo-diazepinone are potential specific inhibitors of the malarial chaperone system. Both itraconazole and ezetimibe are FDA-approved drugs, possess anti-malarial activity and have recently been repurposed for the treatment of cancer. This is the first time that such drug-like compounds have been identified as potential modulators of PfHsp70-PfJDP complexes, and they represent novel candidates for validation and development into anti-malarial drugs.

Expression of the readthrough transcript CiDRE in alveolar macrophages boosts SARS-CoV-2 susceptibility and promotes COVID-19 severity

Lung infection during severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) via the angiotensin-I-converting enzyme 2 (ACE2) receptor induces a cytokine storm. However, the precise mechanisms involved in severe COVID-19 pneumonia are unknown. Here, we showed that interleukin-10 (IL-10) induced the expression of ACE2 in normal alveolar macrophages, causing them to become vectors for SARS-CoV-2. The inhibition of this system in hamster models attenuated SARS-CoV-2 pathogenicity. Genome-wide association and quantitative trait locus analyses identified a IFNAR2-IL10RB readthrough transcript, COVID-19 infectivity-enhancing dual receptor (CiDRE), which was highly expressed in patients harboring COVID-19 risk variants at the IFNAR2 locus. We showed that CiDRE exerted synergistic effects via the IL-10-ACE2 axis in alveolar macrophages and functioned as a decoy receptor for type I interferons. Collectively, our data show that high IL-10 and CiDRE expression are potential risk factors for severe COVID-19. Thus, IL-10R and CiDRE inhibitors might be useful COVID-19 therapies.

Molecular dynamics of mismatch detection-How MutS uses indirect readout to find errors in DNA

The mismatch repair protein MutS safeguards genomic integrity by finding and initiating repair of basepairing errors in DNA. Single-molecule studies show MutS diffusing on DNA, presumably scanning for mispaired/unpaired bases, and crystal structures show a characteristic "mismatch-recognition" complex with DNA enclosed within MutS and kinked at the site of error. But how MutS goes from scanning thousands of Watson-Crick basepairs to recognizing rare mismatches remains unanswered, largely because atomic-resolution data on the search process are lacking. Here, 10 μs all-atom molecular dynamics simulations of Thermus aquaticus MutS bound to homoduplex DNA and T-bulge DNA illuminate the structural dynamics underlying the search mechanism. MutS-DNA interactions constitute a multistep mechanism to check DNA over two helical turns for its 1) shape, through contacts with the sugar-phosphate backbone, 2) conformational flexibility, through bending/unbending engineered by large-scale motions of the clamp domain, and 3) local deformability, through basepair destabilizing contacts. Thus, MutS can localize a potential target by indirect readout due to lower energetic costs of bending mismatched DNA and identify a site that distorts easily due to weaker base stacking and pairing as a mismatch. The MutS signature Phe-X-Glu motif can then lock in the mismatch-recognition complex to initiate repair.

M, Katona G. Survivin prevents the polycomb repressor complex 2 from methylating histone 3 lysine 27

This study investigates the role of survivin in epigenetic control of gene transcription through interaction with the polycomb repressive complex 2 (PRC2). PRC2 is responsible for silencing gene expression by trimethylating lysine 27 on histone 3. We observed differential expression of PRC2 subunits in CD4+ T cells with varying levels of survivin expression, and ChIP-seq results indicated that survivin colocalizes with PRC2 along DNA. Inhibition of survivin resulted in a significant increase in H3K27 trimethylation, implying that survivin prevents PRC2 from functioning. Peptide microarray showed that survivin interacts with peptides from PRC2 subunits, and machine learning revealed that amino acid composition contains relevant information for predicting survivin interaction. NMR and BLI experiments supported the interaction of survivin with PRC2 subunit EZH2. Finally, protein-protein docking revealed that the survivin-EZH2 interaction interface overlaps with catalytic residues of EZH2, potentially inhibiting its H3K27 methylation activity. These findings suggest that survivin inhibits PRC2 function.

Histone H3 lysine 27 crotonylation mediates gene transcriptional repression in chromatin

Histone lysine acylation, including acetylation and crotonylation, plays a pivotal role in gene transcription in health and diseases. However, our understanding of histone lysine acylation has been limited to gene transcriptional activation. Here, we report that histone H3 lysine 27 crotonylation (H3K27cr) directs gene transcriptional repression rather than activation. Specifically, H3K27cr in chromatin is selectively recognized by the YEATS domain of GAS41 in complex with SIN3A-HDAC1 co-repressors. Proto-oncogenic transcription factor MYC recruits GAS41/SIN3A-HDAC1 complex to repress genes in chromatin, including cell-cycle inhibitor p21. GAS41 knockout or H3K27cr-binding depletion results in p21 de-repression, cell-cycle arrest, and tumor growth inhibition in mice, explaining a causal relationship between GAS41 and MYC gene amplification and p21 downregulation in colorectal cancer. Our study suggests that H3K27 crotonylation signifies a previously unrecognized, distinct chromatin state for gene transcriptional repression in contrast to H3K27 trimethylation for transcriptional silencing and H3K27 acetylation for transcriptional activation.

Coarse-grained molecular dynamics simulations of lipid-protein interactions in SLC4 proteins

The SLC4 family of secondary bicarbonate transporters is responsible for the transport of HCO 3 -, CO 3 2- , Cl - , Na + , K + , NH 3 and H + necessary for regulation of pH and ion homeostasis. They are widely expressed in numerous tissues throughout the body and function in different cell types with different membrane properties. Potential lipid roles in SLC4 function have been reported in experimental studies, focusing mostly on two members of the family: AE1 (Cl - /HCO 3 - exchanger) and NBCe1 (Na + -CO 3 2- cotransporter). Previous computational studies of the outward facing (OF) state of AE1 with model lipid membranes revealed enhanced protein-lipid interactions between cholesterol (CHOL) and phosphatidylinositol bisphosphate (PIP2). However, the protein-lipid interactions in other members of the family and other conformation states are still poorly understood and this precludes the detailed studies of a potential regulatory role for lipids in the SLC4 family. In this work, we performed multiple 50 µs coarse-grained molecular dynamics simulations on three members of the SLC4 family with different transport modes: AE1, NBCe1 and NDCBE (a Na + -CO 3 2- /Cl - exchanger), in model HEK293 membranes consisting of CHOL, PIP2, phosphatidylcholine (POPC), phosphatidylethanolamine (POPE), phosphatidylserine (POPS), and sphingomyelin (POSM). The recently resolved inward-facing (IF) state of AE1 was also included in the simulations. Lipid-protein contact analysis of the simulated trajectories was performed with the ProLint server, which provides a multitude of visualization tools for illustration of areas of enhanced lipid-protein contact and identification of putative lipid binding sites within the protein matrix. We observed enrichment of CHOL and PIP2 around all proteins with subtle differences in their distribution depending on the protein type and conformation state. Putative binding sites were identified for CHOL, PIP2, POPC, and POSM in the three studied proteins and their potential roles in the SLC4 transport function, conformational transition and protein dimerization were discussed.

Statement of significance

The SLC4 protein family is involved in critical physiological processes like pH and blood pressure regulation and maintenance of ion homeostasis. Its members can be found in various tissues. A number of studies suggest possible lipid regulation of the SLC4 function. However, the protein-lipid interactions in the SLC4 family are still poorly understood. Here we make use of long coarse-grained molecular dynamics simulations to assess the protein-lipid interactions in three SLC4 proteins with different transport modes, AE1, NBCe1, and NDCBE. We identify putative lipid binding sites for several lipid types of potential mechanistic importance, discuss them in the framework of the known experimental data and provide a necessary basis for further studies on lipid regulation of SLC4 function.