Exact and efficient analytical calculation of the accessible surface areas and their gradients for macromolecules

Protecting Lyophilized Escherichia coli Adenylate Kinase

Drying protein-based drugs, usually via lyophilization, can facilitate storage at ambient temperature and improve accessibility but many proteins cannot withstand drying and must be formulated with protective additives called excipients. However, mechanisms of protection are poorly understood, precluding rational formulation design. To better understand dry proteins and their protection, we examine Escherichia coli adenylate kinase (AdK) lyophilized alone and with the additives trehalose, maltose, bovine serum albumin, cytosolic abundant heat soluble protein D, histidine, and arginine. We apply liquid-observed vapor exchange NMR to interrogate the residue-level structure in the presence and absence of additives. We pair these observations with differential scanning calorimetry data of lyophilized samples and AdK activity assays with and without heating. We show that the amino acids do not preserve the native structure as well as sugars or proteins and that after heating the most stable additives protect activity best.

Deciphering the etiology of undiagnosed ocular anomalies along with systemic alterations in pediatric patients through whole exome sequencing

Inherited and developmental eye diseases are quite diverse and numerous, and determining their genetic cause is challenging due to their high allelic and locus heterogeneity. New molecular approaches, such as whole exome sequencing (WES), have proven to be powerful molecular tools for addressing these cases. The present study used WES to identify the genetic etiology in ten unrelated Mexican pediatric patients with complex ocular anomalies and other systemic alterations of unknown etiology. The WES approach allowed us to identify five clinically relevant variants in the GZF1, NFIX, TRRAP, FGFR2 and PAX2 genes associated with Larsen, Malan, developmental delay with or without dysmorphic facies and autism, LADD1 and papillorenal syndromes. Mutations located in GZF1 and NFIX were classified as pathogenic, those in TRRAP and FGFR2 were classified as likely pathogenic variants, and those in PAX2 were classified as variants of unknown significance. Protein modeling of the two missense FGFR2 p.(Arg210Gln) and PAX2 p.(Met3Thr) variants showed that these changes could induce potential structural alterations in important functional regions of the proteins. Notably, four out of the five variants were not previously reported, except for the TRRAP gene. Consequently, WES enabled the identification of the genetic cause in 40% of the cases reported. All the syndromes reported herein are very rare, with phenotypes that may overlap with other genetic entities.

Drug deconjugation-assisted peptide mapping by LC-MS/MS to identify conjugation sites and quantify site occupancy for antibody-drug conjugates

Antibody-drug conjugates (ADCs) are a heterogeneous mixture of conjugated species with varied drug loadings. Depending on conjugation sites, linkers and drugs can exhibit different stability as influenced by the solvent-accessibility and local charge, resulting in different ADC efficacy, pharmacokinetics, and toxicity. Conjugation site analysis is critical for ADC structural characterization to assure product quality and consistency. It enables early conjugation studies at site-specific levels, confirms the absence of unexpected products to support conjugation process development, and aids in ensuring lot-to-lot consistency for comparability studies. Peptide mapping using liquid chromatography-tandem mass spectrometry is the industry standard method for analyzing conjugation sites. However, some concerns remain for this approach as the large and hydrophobic drug moieties often result in poor MS/MS fragmentation quality and impede the identification of conjugation sites. Additionally, the ionization discrepancy between conjugated and unconjugated peptides can lead to a relatively large bias for site occupancy calculation. In this work, we present a simple drug deconjugation-assisted peptide mapping method to identify and quantify the drug conjugation for ADCs with protease-cleavable linkers. Papain-based drug deconjugation was used to remove the highly hydrophobic drug moiety, which significantly improved the quantitation accuracy of conjugation level and the fragmentation quality. Sample preparation conditions were optimized to avoid introducing artificial modifications, allowing the tracking of initial sample status and subsequent changes of quality attributes during process development and stability assessment. This method was applied to analyze thermally-stressed ADC samples to monitor changes of site-specific conjugation levels, DAR, succinimide hydrolysis of the linker, and various PTMs. We believe this is an effective and straightforward tool for conjugation site analysis while simultaneously monitoring multiple quality attributes for ADCs with protease-cleavable linkers.

Three-Dimensional Interaction Homology: Deconstructing Residue-Residue and Residue-Lipid Interactions in Membrane Proteins

A method is described to deconstruct the network of hydropathic interactions within and between a protein's sidechain and its environment into residue-based three-dimensional maps. These maps encode favorable and unfavorable hydrophobic and polar interactions, in terms of spatial positions for optimal interactions, relative interaction strength, as well as character. In addition, these maps are backbone angle-dependent. After map calculation and clustering, a finite number of unique residue sidechain interaction maps exist for each backbone conformation, with the number related to the residue's size and interaction complexity. Structures for soluble proteins (~749,000 residues) and membrane proteins (~387,000 residues) were analyzed, with the latter group being subdivided into three subsets related to the residue's position in the membrane protein: soluble domain, core-facing transmembrane domain, and lipid-facing transmembrane domain. This work suggests that maps representing residue types and their backbone conformation can be reassembled to optimize the medium-to-high resolution details of a protein structure. In particular, the information encoded in maps constructed from the lipid-facing transmembrane residues appears to paint a clear picture of the protein-lipid interactions that are difficult to obtain experimentally.

Predicting inhibitor development using a random peptide phage-display library approach in the SIPPET cohort

ABSTRACT

Inhibitor development is the most severe complication of hemophilia A (HA) care and is associated with increased morbidity and mortality. This study aimed to use a novel immunoglobulin G epitope mapping method to explore the factor VIII (FVIII)-specific epitope profile in the SIPPET cohort population and to develop an epitope mapping-based inhibitor prediction model. The population consisted of 122 previously untreated patients with severe HA who were followed up for 50 days of exposure to FVIII or 3 years, whichever occurred first. Sampling was performed before FVIII treatment and at the end of the follow-up. The outcome was inhibitor development. The FVIII epitope repertoire was assessed by means of a novel random peptide phage-display assay. A least absolute shrinkage and selection operator (LASSO) regression model and a random forest model were fitted on posttreatment sample data and validated in pretreatment sample data. The predictive performance of these models was assessed by the C-statistic and a calibration plot. We identified 27 775 peptides putatively directed against FVIII, which were used as input for the statistical models. The C-statistic of the LASSO and random forest models were good at 0.78 (95% confidence interval [CI], 0.69-0.86) and 0.80 (95% CI, 0.72-0.89). Model calibration of both models was moderately good. Two statistical models, developed on data from a novel random peptide phage display assay, were used to predict inhibitor development before exposure to exogenous FVIII. These models can be used to set up diagnostic tests that predict the risk of inhibitor development before starting treatment with FVIII.

Computer-aided rational design strategy based on protein surface charge to improve the thermal stability of a novel esterase from Geobacillus jurassicus

OBJECTIVES

Although Geobacillus are significant thermophilic bacteria source, there are no reports of thermostable esterase gene in Geobacillus jurassicus or rational design strategies to increase the thermal stability of esterases.

RESULTS

Gene gju768 showed a highest similarity of 15.20% to esterases from Geobacillus sp. with detail enzymatic properties. Using a combination of Gibbs Unfolding Free Energy (∆∆G) calculator and the distance from the mutation site to the catalytic site (DsCα-Cα) to screen suitable mutation sites with elimination of negative surface charge, the mutants (D24N, E221Q, and E253Q) displayed stable mutants with higher thermal stability than the wild-type (WT). Mutant E253Q exhibited the best thermal stability, with a half-life (T1/2) at 65 °C of 32.4 min, which was 1.8-fold of the WT (17.9 min).

CONCLUSION

Cloning of gene gju768 and rational design based on surface charge engineering contributed to the identification of thermostable esterase from Geobacillus sp. and the exploration of evolutionary strategies for thermal stability.

Biochemical, structural, and computational analyses of two new clinically identified missense mutations of ALDH7A1

Aldehyde dehydrogenase 7A1 (ALDH7A1) catalyzes a step of lysine catabolism. Certain missense mutations in the ALDH7A1 gene cause pyridoxine dependent epilepsy (PDE), a rare autosomal neurometabolic disorder with recessive inheritance that affects almost 1:65,000 live births and is classically characterized by recurrent seizures from the neonatal period. We report a biochemical, structural, and computational study of two novel ALDH7A1 missense mutations that were identified in a child with rare recurrent seizures from the third month of life. The mutations affect two residues in the oligomer interfaces of ALDH7A1, Arg134 and Arg441 (Arg162 and Arg469 in the HGVS nomenclature). The corresponding enzyme variants R134S and R441C (p.Arg162Ser and p.Arg469Cys in the HGVS nomenclature) were expressed in Escherichia coli and purified. R134S and R441C have 10,000- and 50-fold lower catalytic efficiency than wild-type ALDH7A1, respectively. Sedimentation velocity analytical ultracentrifugation shows that R134S is defective in tetramerization, remaining locked in a dimeric state even in the presence of the tetramer-inducing coenzyme NAD+. Because the tetramer is the active form of ALDH7A1, the defect in oligomerization explains the very low catalytic activity of R134S. In contrast, R441C exhibits wild-type oligomerization behavior, and the 2.0 Å resolution crystal structure of R441C complexed with NAD+ revealed no obvious structural perturbations when compared to the wild-type enzyme structure. Molecular dynamics simulations suggest that the mutation of Arg441 to Cys may increase intersubunit ion pairs and alter the dynamics of the active site gate. Our biochemical, structural, and computational data on two novel clinical variants of ALDH7A1 add to the complexity of the molecular determinants underlying pyridoxine dependent epilepsy.

Ago2/CAV1 interaction potentiates metastasis via controlling Ago2 localization and miRNA action

Ago2 differentially regulates oncogenic and tumor-suppressive miRNAs in cancer cells. This discrepancy suggests a secondary event regulating Ago2/miRNA action in a context-dependent manner. We show here that a positive charge of Ago2 K212, that is preserved by SIR2-mediated Ago2 deacetylation in cancer cells, is responsible for the direct interaction between Ago2 and Caveolin-1 (CAV1). Through this interaction, CAV1 sequesters Ago2 on the plasma membranes and regulates miRNA-mediated translational repression in a compartment-dependent manner. Ago2/CAV1 interaction plays a role in miRNA-mediated mRNA suppression and in miRNA release via extracellular vesicles (EVs) from tumors into the circulation, which can be used as a biomarker of tumor progression. Increased Ago2/CAV1 interaction with tumor progression promotes aggressive cancer behaviors, including metastasis. Ago2/CAV1 interaction acts as a secondary event in miRNA-mediated suppression and increases the complexity of miRNA actions in cancer.

Protein painting for structural and binding site analysis via intracellular lysine reactivity profiling with o-phthalaldehyde

The three-dimensional structure and the molecular interaction of proteins determine their roles in many cellular processes. Chemical protein painting with protein mass spectrometry can identify changes in structural conformations and molecular interactions of proteins including their binding sites. Nevertheless, most current protein painting techniques identify protein targets and binding sites of drugs in vitro using a cell lysate or purified protein. Here, we tested 11 membrane-permeable lysine-reactive chemical probes for intracellular covalent labeling of endogenous proteins, which reveals ortho-phthalaldehyde (OPA) as the most reactive probe in the intracellular environment. An MS workflow and a new data analysis strategy termed RAPID (Reactive Amino acid Profiling by Inverse Detection) was developed to enhance detection sensitivity. RAPID with OPA successfully identified structural changes induced by the allosteric drug TEPP-46 on its target protein PKM2 and was applied to profile the conformation change of the proteome occurring in cells during thermal denaturation. The application of RAPID-OPA on cells treated with geldanamycin, selumetinib, and staurosporine successfully revealed their binding sites on target proteins. Thus, RAPID-OPA for cellular protein painting enables the identification of ligand-binding sites and detection of protein structural changes occurring in cells.

N4BP1 functions as a dimerization-dependent linear ubiquitin reader which regulates TNF signalling

Signalling through TNFR1 modulates proinflammatory gene transcription and programmed cell death, and its impairment causes autoimmune diseases and cancer. NEDD4-binding protein 1 (N4BP1) is a critical suppressor of proinflammatory cytokine production that acts as a regulator of innate immune signalling and inflammation. However, our current understanding about the molecular properties that enable N4BP1 to exert its suppressive potential remain limited. Here, we show that N4BP1 is a novel linear ubiquitin reader that negatively regulates NFκB signalling by its unique dimerization-dependent ubiquitin-binding module that we named LUBIN. Dimeric N4BP1 strategically positions two non-selective ubiquitin-binding domains to ensure preferential recognition of linear ubiquitin. Under proinflammatory conditions, N4BP1 is recruited to the nascent TNFR1 signalling complex, where it regulates duration of proinflammatory signalling in LUBIN-dependent manner. N4BP1 deficiency accelerates TNFα-induced cell death by increasing complex II assembly. Under proapoptotic conditions, caspase-8 mediates proteolytic processing of N4BP1, resulting in rapid degradation of N4BP1 by the 26 S proteasome, and acceleration of apoptosis. In summary, our findings demonstrate that N4BP1 dimerization creates a novel type of ubiquitin reader that selectively recognises linear ubiquitin which enables the timely and coordinated regulation of TNFR1-mediated inflammation and cell death.

Design, development, and validation of multi-epitope proteins for serological diagnosis of Zika virus infections and discrimination from dengue virus seropositivity

Zika virus (ZIKV), an arbovirus from the Flaviviridae family, is the causative agent of Zika fever, a mild and frequent oligosymptomatic disease in humans. Nonetheless, on rare occasions, ZIKV infection can be associated with Guillain-Barré Syndrome (GBS), and severe congenital complications, such as microcephaly. The oligosymptomatic disease, however, presents symptoms that are quite similar to those observed in infections caused by other frequent co-circulating arboviruses, including dengue virus (DENV). Moreover, the antigenic similarity between ZIKV and DENV, and even with other members of the Flaviviridae family, complicates serological testing due to the high cross-reactivity of antibodies. Here, we designed, produced in a prokaryotic expression system, and purified three multiepitope proteins (ZIKV-1, ZIKV-2, and ZIKV-3) for differential diagnosis of Zika. The proteins were evaluated as antigens in ELISA tests for the detection of anti-ZIKV IgG using ZIKV- and DENV-positive human sera. The recombinant proteins were able to bind and detect anti-ZIKV antibodies without cross-reactivity with DENV-positive sera and showed no reactivity with Chikungunya virus (CHIKV)- positive sera. ZIKV-1, ZIKV-2, and ZIKV-3 proteins presented 81.6%, 95%, and 66% sensitivity and 97%, 96%, and 84% specificity, respectively. Our results demonstrate the potential of the designed and expressed antigens in the development of specific diagnostic tests for the detection of IgG antibodies against ZIKV, especially in regions with the circulation of multiple arboviruses.

Viral nanoparticles: Current advances in design and development

Viral nanoparticles (VNPs) are self-assembling, adaptable delivery systems for vaccines and other therapeutic agents used in a variety of biomedical applications. The potential of viruses to invade and infect various hosts and cells renders them suitable as potential nanocarriers, possessing distinct functional characteristics, immunogenic properties, and improved biocompatibility and biodegradability. VNPs are frequently produced through precise genetic or chemical engineering, which involves adding diverse sequences or functional payloads to the capsid protein (CP). Several spherical and helical plant viruses, bacteriophages, and animal viruses are currently being used as VNPs, or non-infectious virus-like particles (VLPs). In addition to their broad use in cancer therapy, vaccine technology, diagnostics, and molecular imaging, VNPs have made important strides in the realms of tissue engineering, biosensing, and antimicrobial prophylaxis. They are also being used in energy storage cells due to their binding and piezoelectric properties. The large-scale production of VNPs for research, preclinical testing, and clinical use is fraught with difficulties, such as those relating to cost-effectiveness, scalability, and purity. Consequently, many plants- and microorganism-based platforms are being developed, and newer viruses are being explored. The goal of the current review is to provide an overview of these advances.

Traversing DNA-Protein Interactions Between Mesophilic and Thermophilic Bacteria: Implications from Their Cold Shock Response

Cold shock proteins (CSPs) are small, acidic proteins which contain a conserved nucleic acid-binding domain. These perform mRNA translation acting as "RNA chaperones" when triggered by low temperatures initiating their cold shock response. CSP- RNA interactions have been predominantly studied. Our focus will be CSP-DNA interaction examination, to analyse the diverse interaction patterns such as electrostatic, hydrogen and hydrophobic bonding in both thermophilic and mesophilic bacteria. The differences in the molecular mechanism of these contrasting bacterial proteins are studied. Computational techniques such as modelling, energy refinement, simulation and docking were operated to obtain data for comparative analysis. The thermostability factors which stabilise a thermophilic bacterium and their effect on their molecular regulation is investigated. Conformational deviation, atomic residual fluctuations, binding affinity, Electrostatic energy and Solvent Accessibility energy were determined during stimulation along with their conformational study. The study revealed that mesophilic bacteria E. coli CSP have higher binding affinity to DNA than thermophilic G. stearothermophilus. This was further evident by low conformation deviation and atomic fluctuations during simulation.

The unique salt bridge network in GlacPETase: a key to its stability

The extensive accumulation of polyethylene terephthalate (PET) has become a critical environmental issue. PET hydrolases can break down PET into its building blocks. Recently, we identified a glacial PET hydrolase GlacPETase sharing less than 31% amino acid identity with any known PET hydrolases. In this study, the crystal structure of GlacPETase was determined at 1.8 Å resolution, revealing unique structural features including a distinctive N-terminal disulfide bond and a specific salt bridge network. Site-directed mutagenesis demonstrated that the disruption of the N-terminal disulfide bond did not reduce GlacPETase's thermostability or its catalytic activity on PET. However, mutations in the salt bridges resulted in changes in melting temperature ranging from -8°C to +2°C and the activity on PET ranging from 17.5% to 145.5% compared to the wild type. Molecular dynamics simulations revealed that these salt bridges stabilized the GlacPETase's structure by maintaining their surrounding structure. Phylogenetic analysis indicated that GlacPETase represented a distinct branch within PET hydrolases-like proteins, with the salt bridges and disulfide bonds in this branch being relatively conserved. This research contributed to the improvement of our comprehension of the structural mechanisms that dictate the thermostability of PET hydrolases, highlighting the diverse characteristics and adaptability observed within PET hydrolases.IMPORTANCEThe pervasive problem of polyethylene terephthalate (PET) pollution in various terrestrial and marine environments is widely acknowledged and continues to escalate. PET hydrolases, such as GlacPETase in this study, offered a solution for breaking down PET. Its unique origin and less than 31% identity with any known PET hydrolases have driven us to resolve its structure. Here, we report the correlation between its unique structure and biochemical properties, focusing on an N-terminal disulfide bond and specific salt bridges. Through site-directed mutagenesis experiments and molecular dynamics simulations, the roles of the N-terminal disulfide bond and salt bridges were elucidated in GlacPETase. This research enhanced our understanding of the role of salt bridges in the thermostability of PET hydrolases, providing a valuable reference for the future engineering of PET hydrolases.

Residue coevolution and mutational landscape for OmpR and NarL response regulator subfamilies

DNA-binding response regulators (DBRRs) are a broad class of proteins that operate in tandem with their partner kinase proteins to form two-component signal transduction systems in bacteria. Typical DBRRs are composed of two domains where the conserved N-terminal domain accepts transduced signals and the evolutionarily diverse C-terminal domain binds to DNA. These domains are assumed to be functionally independent, and hence recombination of the two domains should yield novel DBRRs of arbitrary input/output response, which can be used as biosensors. This idea has been proved to be successful in some cases; yet, the error rate is not trivial. Improvement of the success rate of this technique requires a deeper understanding of the linker-domain and inter-domain residue interactions, which have not yet been thoroughly examined. Here, we studied residue coevolution of DBRRs of the two main subfamilies (OmpR and NarL) using large collections of bacterial amino acid sequences to extensively investigate the evolutionary signatures of linker-domain and inter-domain residue interactions. Coevolutionary analysis uncovered evolutionarily selected linker-domain and inter-domain residue interactions of known experimental structures, as well as previously unknown inter-domain residue interactions. We examined the possibility of these inter-domain residue interactions as contacts that stabilize an inactive conformation of the DBRR where DNA binding is inhibited for both subfamilies. The newly gained insights on linker-domain/inter-domain residue interactions and shared inactivation mechanisms improve the understanding of the functional mechanism of DBRRs, providing clues to efficiently create functional DBRR-based biosensors. Additionally, we show the feasibility of applying coevolutionary landscape models to predict the functionality of domain-swapped DBRR proteins. The presented result demonstrates that sequence information can be used to filter out bioengineered DBRR proteins that are predicted to be nonfunctional due to a high negative predictive value.

Regulatory sites of CaM-sensitive adenylyl cyclase AC8 revealed by cryo-EM and structural proteomics

Membrane adenylyl cyclase AC8 is regulated by G proteins and calmodulin (CaM), mediating the crosstalk between the cAMP pathway and Ca2+ signalling. Despite the importance of AC8 in physiology, the structural basis of its regulation by G proteins and CaM is not well defined. Here, we report the 3.5 Å resolution cryo-EM structure of the bovine AC8 bound to the stimulatory Gαs protein in the presence of Ca2+/CaM. The structure reveals the architecture of the ordered AC8 domains bound to Gαs and the small molecule activator forskolin. The extracellular surface of AC8 features a negatively charged pocket, a potential site for unknown interactors. Despite the well-resolved forskolin density, the captured state of AC8 does not favour tight nucleotide binding. The structural proteomics approaches, limited proteolysis and crosslinking mass spectrometry (LiP-MS and XL-MS), allowed us to identify the contact sites between AC8 and its regulators, CaM, Gαs, and Gβγ, as well as to infer the conformational changes induced by these interactions. Our results provide a framework for understanding the role of flexible regions in the mechanism of AC regulation.

Disaggregation of amyloid-beta fibrils via natural metabolites using long timescale replica exchange molecular dynamics simulation studies

CONTEXT

Amyloid fibrils are self-assembled fibrous protein aggregates that are associated with several presently incurable diseases such as Alzheimer's. disease that is characterized by the accumulation of amyloid fibrils in the brain, which leads to the formation of plaques and the death of brain cells. Disaggregation of amyloid fibrils is considered a promising approach to cure Alzheimer's disease. The mechanism of amyloid fibril formation is complex and not fully understood, making it difficult to develop drugs that can target the process. Diacetonamine and cystathionine are potential lead compounds to induce disaggregation of amyloid fibrils.

METHODS

In the current research, we have used long timescale molecular simulation studies and replica exchange molecular dynamics (REMD) for 1000 ns (1 μs) to examine the mechanisms by which natural metabolites can disaggregate amyloid-beta fibrils. Molecular docking was carried out using Glide and with prior protein minimization and ligand preparation. We focused on a screening a database of natural metabolites, as potential candidates for disaggregating amyloid fibrils. We used Desmond with OPLS 3e as a force field. MM-GBSA calculations were performed. Blood-brain barrier permeability, SASA, and radius of gyration parameters were calculated.

Viper Venom Phospholipase A2 Database: The Structural and Functional Anatomy of a Primary Toxin in Envenomation

Viper venom phospholipase A2 enzymes (vvPLA2s) and phospholipase A2-like (PLA2-like) proteins are two of the principal toxins in viper venom that are responsible for the severe myotoxic and neurotoxic effects caused by snakebite envenoming, among other pathologies. As snakebite envenoming is the deadliest neglected tropical disease, a complete understanding of these proteins' properties and their mechanisms of action is urgently needed. Therefore, we created a database comprising information on the holo-form, cofactor-bound 3D structure of 217 vvPLA2 and PLA2-like proteins in their physiologic environment, as well as 79 membrane-bound viper species from 24 genera, which we have made available to the scientific community to accelerate the development of new anti-snakebite drugs. In addition, the analysis of the sequenced, 3D structure of the database proteins reveals essential aspects of the anatomy of the proteins, their toxicity mechanisms, and the conserved binding site areas that may anchor universal interspecific inhibitors. Moreover, it pinpoints hypotheses for the molecular origin of the myotoxicity of the PLA2-like proteins. Altogether, this study provides an understanding of the diversity of these toxins and how they are conserved, and it indicates how to develop broad, interspecies, efficient small-molecule inhibitors to target the toxin's many mechanisms of action.

Insights into the pathogenesis of primary hyperoxaluria type I from the structural dynamics of alanine:glyoxylate aminotransferase variants

Primary hyperoxaluria type I (PH1) is caused by deficient alanine:glyoxylate aminotransferase (AGT) activity. PH1-causing mutations in AGT lead to protein mistargeting and aggregation. Here, we use hydrogen-deuterium exchange (HDX) to characterize the wild-type (WT), the LM (a polymorphism frequent in PH1 patients) and the LM G170R (the most common mutation in PH1) variants of AGT. We provide the first experimental analysis of AGT structural dynamics, showing that stability is heterogeneous in the native state and providing a blueprint for frustrated regions with potentially functional relevance. The LM and LM G170R variants only show local destabilization. Enzymatic transamination of the pyridoxal 5-phosphate cofactor bound to AGT hardly affects stability. Our study, thus, supports that AGT misfolding is not caused by dramatic effects on structural dynamics.

Discovery and Identification of a Novel Tag of HlyA60 for Protein Active Aggregate Formation in Escherichia coli

The strategy of active aggregation tag fusion expression with target proteins can solve the problems of restricted expression, inefficient purification, and laborious immobilization faced in the production of recombinant proteins in Escherichia coli. We localized a novel active aggregation peptide HlyA60 from the hemolysin A secretion system, which can effectively induce aggregate formation with satisfactory protein activities in E. coli after fusion expression with the protein of interest. Based on structural prediction and surface properties, the process of active aggregation of HlyA60 through electrostatic interactions and hydrophobic interactions was analyzed. To investigate the potential application of HlyA60 as an efficient aggregation tag, it was fused with acetyl xylan esterase and lipase A, separately. The resulting fusion proteins demonstrated active aggregation rates of 97.6 and 66.7%, respectively, leading to 1.9-fold and 1.7-fold increases in bacterial density at the end of fermentation. The AXE-HlyA60 fusion protein, which exhibited superior performance, was subjected to purification and immobilization. It was able to achieve column-free purification with an impressive 98.8% recovery and in situ immobilization; the immobilization enabled 30 cycles of reactions to take place with 85% residual activity maintained. Our findings provide a novel tool for efficiently producing recombinant proteins in E. coli.

Dissecting the Heterogeneous Glycan Profiles of Recombinant Coronavirus Spike Proteins with Individual Ion Mass Spectrometry

Surface-embedded glycoproteins, such as the spike protein trimers of coronaviruses MERS, SARS-CoV, and SARS-CoV-2, play a key role in viral function and are the target antigen for many vaccines. However, their significant glycan heterogeneity poses an analytical challenge. Here, we utilized individual ion mass spectrometry (I2MS), a multiplexed charge detection measurement with similarities to charge detection mass spectrometry (CDMS), in which a commercially available Orbitrap analyzer is used to directly produce mass profiles of these heterogeneous coronavirus spike protein trimers under native-like conditions. Analysis by I2MS shows that glycosylation contributes to the molecular mass of each protein trimer more significantly than expected by bottom-up techniques, highlighting the importance of obtaining complementary intact mass information when characterizing glycosylation of such heterogeneous proteins. Enzymatic dissection to remove sialic acid or N-linked glycans demonstrates that I2MS can be used to better understand the glycan profile from a native viewpoint. Deglycosylation of N-glycans followed by I2MS analysis indicates that the SARS-CoV-2 spike protein trimer contains glycans that are more difficult to remove than its MERS and SARS-CoV counterparts, and these differences are correlated with solvent accessibility. I2MS technology enables characterization of protein mass and intact glycan profile and is orthogonal to traditional mass analysis methods such as size exclusion chromatography-multiangle light scattering (SEC-MALS) and field flow fractionation-multiangle light scattering (FFF-MALS). An added advantage of I2MS is low sample use, requiring 100-fold less than other methodologies. This work highlights how I2MS technology can enable efficient development of vaccines and therapeutics for pharmaceutical development.

Chemoproteomic strategy identified p120-catenin glutathionylation regulates E-cadherin degradation and cell migration

Identification of cysteines with high oxidation susceptibility is important for understanding redox-mediated biological processes. In this report, we report a chemical proteomic strategy that finds cysteines with high susceptibility to S-glutathionylation. Our proteomic strategy, named clickable glutathione-based isotope-coded affinity tag (G-ICAT), identified 1,518 glutathionylated cysteines while determining their relative levels of glutathionylated and reduced forms upon adding hydrogen peroxide. Among identified cysteines, we demonstrated that CTNND1 (p120) C692 has high susceptibility to glutathionylation. Also, p120 wild type (WT), compared to C692S, induces its dissociation from E-cadherin under oxidative stress, such as glucose depletion. p120 and E-cadherin dissociation correlated with E-cadherin destabilization via its proteasomal degradation. Lastly, we showed that p120 WT, compared to C692S, increases migration and invasion of MCF7 cells under glucose depletion, supporting a model that p120 C692 glutathionylation increases cell migration and invasion by destabilization of E-cadherin, a core player in cell-cell adhesion.

Systematic Fe(II)-EDTA Method of Dose-Dependent Hydroxyl Radical Generation for Protein Oxidative Footprinting

Correlating the structure and dynamics of proteins with biological function is critical to understanding normal and dysfunctional cellular mechanisms. We describe a quantitative method of hydroxyl radical generation via Fe(II)-ethylenediaminetetraacetic acid (EDTA)-catalyzed Fenton chemistry that provides ready access to protein oxidative footprinting using equipment commonly found in research and process control laboratories. Robust and reproducible dose-dependent oxidation of protein samples is observed and quantitated by mass spectrometry with as fine a single residue resolution. An oxidation analysis of lysozyme provides a readily accessible benchmark for our method. The efficacy of our oxidation method is demonstrated by mapping the interface of a RAS-monobody complex, the surface of the NIST mAb, and the interface between PRC2 complex components. These studies are executed using standard laboratory tools and a few pennies of reagents; the mass spectrometry analysis can be streamlined to map the protein structure with single amino acid residue resolution.

The structure of a Lactobacillus helveticus chlorogenic acid esterase and the dynamics of its insertion domain provide insights into substrate binding

Chlorogenic acid esterases (ChlEs) are a useful class of enzymes that hydrolyze chlorogenic acid (CGA) into caffeic and quinic acids. ChlEs can break down CGA in foods to improve their sensory properties and release caffeic acid in the digestive system to improve the absorption of bioactive compounds. This work presents the structure, molecular dynamics, and biochemical characterization of a ChlE from Lactobacillus helveticus (Lh). Molecular dynamics simulations suggest that substrate access to the active site of LhChlE is modulated by two hairpin loops above the active site. Docking simulations and mutational analysis suggest that two residues within the loops, Gln145 and Lys164 , are important for CGA binding. Lys164 provides a slight substrate preference for CGA, whereas Gln145 is required for efficient turnover. This work is the first to examine the dynamics of a bacterial ChlE and provides insights on substrate binding preference and turnover in this type of enzyme.

Sequence basis for selectivity of ephrin-B2 ligand for Eph receptors and pathogenic henipavirus G glycoproteins

IMPORTANCE

Ephrin-B2 (EFNB2) is a ligand for six Eph receptors in humans and regulates multiple cell developmental and signaling processes. It also functions as the cell entry receptor for Nipah virus and Hendra virus, zoonotic viruses that can cause respiratory and/or neurological symptoms in humans with high mortality. Here, we investigate the sequence basis of EFNB2 specificity for binding the Nipah virus attachment G glycoprotein over Eph receptors. We then use this information to engineer EFNB2 as a soluble decoy receptor that specifically binds the attachment glycoproteins of the Nipah virus and other related henipaviruses to neutralize infection. These findings further mechanistic understanding of protein selectivity and may facilitate the development of diagnostics or therapeutics against henipavirus infection.

Periplasmic chitooligosaccharide-binding protein requires a three-domain organization for substrate translocation

Periplasmic solute-binding proteins (SBPs) specific for chitooligosaccharides, (GlcNAc)n (n = 2, 3, 4, 5 and 6), are involved in the uptake of chitinous nutrients and the negative control of chitin signal transduction in Vibrios. Most translocation processes by SBPs across the inner membrane have been explained thus far by two-domain open/closed mechanism. Here we propose three-domain mechanism of the (GlcNAc)n translocation based on experiments using a recombinant VcCBP, SBP specific for (GlcNAc)n from Vibrio cholerae. X-ray crystal structures of unliganded or (GlcNAc)3-liganded VcCBP solved at 1.2-1.6 Å revealed three distinct domains, the Upper1, Upper2 and Lower domains for this protein. Molecular dynamics simulation indicated that the motions of the three domains are independent and that in the (GlcNAc)3-liganded state the Upper2/Lower interface fluctuated more intensively, compared to the Upper1/Lower interface. The Upper1/Lower interface bound two GlcNAc residues tightly, while the Upper2/Lower interface appeared to loosen and release the bound sugar molecule. The three-domain mechanism proposed here was fully supported by binding data obtained by thermal unfolding experiments and ITC, and may be applicable to other translocation systems involving SBPs belonging to the same cluster.

The updated Structural Database of Allergenic Proteins (SDAP 2.0) provides 3D models for allergens and incorporated bioinformatics tools

Background

Allergenic proteins can cause IgE-mediated adverse reactions in sensitized individuals. Although the sequences of many allergenic proteins have been identified, bioinformatics data analysis with advanced computational methods and modeling is needed to identify the basis for IgE binding and cross-reactivity.

Objective

We aim to present the features and use of the updated Structural Database of Allergenic Proteins 2.0 (SDAP 2.0) webserver, a unique, publicly available resource to compare allergens using specially designed computational tools and new high-quality 3-D models for most known allergens.

Methods

Previously developed and novel software tools for identifying cross-reactive allergens using sequence and structure similarity are implemented in SDAP 2.0. A comprehensive set of high-quality 3-D models of most allergens was generated with the state-of-the-art AlphaFold 2 software. A graphics tool enables the interactive visualization of IgE epitopes on experimentally determined and modeled 3-D structures.

Results

A user can search for allergens similar to a given input sequence with the FASTA algorithm or the window-based World Health Organization/International Union of Immunological Societies (WHO/IUIS) guidelines on safety concerns of novel food products. Peptides similar to known IgE epitopes can be identified with the property distance tool and conformational epitopes by the Cross-React method. The updated database contains 1657 manually curated sequences including all allergens from the IUIS database, 334 experimentally determined X-ray or NMR structures, and 1565 3-D models. Each allergen/isoallergen is classified according to its protein family.

Conclusions

SDAP provides access to the steadily increasing information on allergenic structures and epitopes with integrated bioinformatics tools to identify and analyze their similarities. In addition to serving the research and regulatory community, it provides clinicians with tools to identify potential coallergies in a sensitive patient and can help companies to design hypoallergenic foods and immunotherapies.

Side-chain dynamics of the α1B -adrenergic receptor determined by NMR via methyl relaxation

G protein-coupled receptors (GPCRs) are medically important membrane proteins that sample inactive, intermediate, and active conformational states characterized by relatively slow interconversions (~μs-ms). On a faster timescale (~ps-ns), the conformational landscape of GPCRs is governed by the rapid dynamics of amino acid side chains. Such dynamics are essential for protein functions such as ligand recognition and allostery. Unfortunately, technical challenges have almost entirely precluded the study of side-chain dynamics for GPCRs. Here, we investigate the rapid side-chain dynamics of a thermostabilized α1B -adrenergic receptor (α1B -AR) as probed by methyl relaxation. We determined order parameters for Ile, Leu, and Val methyl groups in the presence of inverse agonists that bind orthosterically (prazosin, tamsulosin) or allosterically (conopeptide ρ-TIA). Despite the differences in the ligands, the receptor's overall side-chain dynamics are very similar, including those of the apo form. However, ρ-TIA increases the flexibility of Ile1764×56 and possibly of Ile2145×49 , adjacent to Pro2155×50 of the highly conserved P5×50 I3×40 F6×44 motif crucial for receptor activation, suggesting differences in the mechanisms for orthosteric and allosteric receptor inactivation. Overall, increased Ile side-chain rigidity was found for residues closer to the center of the membrane bilayer, correlating with denser packing and lower protein surface exposure. In contrast to two microbial membrane proteins, in α1B -AR Leu exhibited higher flexibility than Ile side chains on average, correlating with the presence of Leu in less densely packed areas and with higher protein-surface exposure than Ile. Our findings demonstrate the feasibility of studying receptor-wide side-chain dynamics in GPCRs to gain functional insights.

Computational analysis of carboxylesterase genes and proteins in non-pathogenic food bacterium Alicyclobacillus acidocaldarius: insights from proteogenomics

Over recent years, Alicyclobacillus acidocaldarius, a Gram-positive nonpathogenic rod-shaped thermo-acid-tolerant bacterium, has posed numerous challenges for the fruit juice industry. However, the bacterium's unique characteristics, particularly its nonpathogenic and thermophilic capabilities, offer significant opportunities for genetic exploration by biotechnologists. This study presents the computational proteogenomics report on the carboxylesterase (CE) enzyme in A. acidocaldarius, shedding light on structural and evolutional of CEs from this bacterium. Our analysis revealed that the average molecular weight of CEs in A. acidocaldarius was 41 kDa, with an isoelectric point around 5. The amino acid composition favored negative amino acids over positive ones. The aliphatic index and hydropathicity were approximately 88 and - 0.15, respectively. While the protein sequence showed no disulfide bonds in the CEs' structure, the presence of Cys amino acids was observed in the structure of CEs. Phylogenetic analysis presented more than 99% similarity between CEs, indicating their close evolutionary relationship. By applying homology modeling, the 3-dimensional structural models of the carboxylesterase were constructed, which with the help of structural conservation and solvent accessibility analysis highlighted key residues and regions responsible for enzyme stability and conformation. The specific patterns presented the total solvent accessibility of less than 25 (Å2) was in considerable position as well as Gly residues were noticeably have high accessibility to solvent in all structures. Ala was the more frequent amino acids in the conserved-SASA of carboxylesterases. Furthermore, unsupervised agglomerative hierarchical clustering based on solvent accessibility feature successfully clustered and even distinguished this enzyme from proteases from the same genome. These findings contribute to a deeper understanding of the nonpathogenic A. acidocaldarius carboxylesterase and its potential applications in biotechnology. Additionally, structural analysis of CEs would help to address potential solutions in fruit juice industry with utilization of computational structural biology.

SARS-CoV-2 antibodies recognize 23 distinct epitopic sites on the receptor binding domain

The COVID-19 pandemic and SARS-CoV-2 variants have dramatically illustrated the need for a better understanding of antigen (epitope)-antibody (paratope) interactions. To gain insight into the immunogenic characteristics of epitopic sites (ES), we systematically investigated the structures of 340 Abs and 83 nanobodies (Nbs) complexed with the Receptor Binding Domain (RBD) of the SARS-CoV-2 spike protein. We identified 23 distinct ES on the RBD surface and determined the frequencies of amino acid usage in the corresponding CDR paratopes. We describe a clustering method for analysis of ES similarities that reveals binding motifs of the paratopes and that provides insights for vaccine design and therapies for SARS-CoV-2, as well as a broader understanding of the structural basis of Ab-protein antigen (Ag) interactions.

Isolation, Structure Elucidation and in Vitro Anticancer Activity of Phytochemical Constituents of Goniothalamus wynaadensis Bedd. and Identification of α-Tubulin as a Putative Molecular Target by in Silico Study

The phytochemical analysis of ethyl acetate and methanol extract of Goniothalamus wynaadensis Bedd. leaves led to an isolation of eight (1-8) known molecules, among them seven (2-8) isolated for the first time from this species, which includes (+)-goniothalamin oxide (2), goniodiol-7-monoacetate (3), goniodiol-8-monoacetate (4), goniodiol (5), (+)-8-epi-9-deoxygoniopypyrone (6) etc. The phytochemical modification by acetylation of 3 and 4 gave goniodiol diacetate (9) with absolute configuration (6R, 7R, 8R) confirmed by single crystal X-ray diffraction. Compounds 3-9 were cytotoxic against breast, ovarian, prostate and colon cancer cell lines with IC50 <10 μM. Cell cycle analysis and Annexin-V assay on MDA-MB-231 cell using goniodiol-7-monoacetate (3) exhibited apoptotic response as well as necrotic response and showed cell proliferation arrest at G2/M phase. An in silico target identification for these molecules was carried out with an α-tubulin protein target by covalent docking. To gain an in-depth understanding and identify the stability of these protein-ligand complexes on thermodynamic energy levels, further assessment of the isolated molecules binding to the Cys-316 of α-tubulin was performed based on reaction energetic analysis via DFT studies which hinted the isolated molecules may be α-tubulin inhibitors similar to Pironetin. Molecular dynamics reiterated the observations.

Changes in Hemoglobin Properties in Complex with Glutathione and after Glutathionylation

Hemoglobin is the main protein of red blood cells that provides oxygen transport to all cells of the human body. The ability of hemoglobin to bind the main low-molecular-weight thiol of the cell glutathione, both covalently and noncovalently, is not only an important part of the antioxidant protection of red blood cells, but also affects its affinity for oxygen in both cases. In this study, the properties of oxyhemoglobin in complex with reduced glutathione (GSH) and properties of glutathionylated hemoglobin bound to glutathione via an SS bond were characterized. For this purpose, the methods of circular dichroism, Raman spectroscopy, infrared spectroscopy, tryptophan fluorescence, differential scanning fluorimetry, and molecular modeling were used. It was found that the glutathionylation of oxyhemoglobin caused changes in the secondary structure of the protein, reducing the alpha helicity, but did not affect the heme environment, tryptophan fluorescence, and the thermostability of the protein. In the noncovalent complex of oxyhemoglobin with reduced glutathione, the secondary structure of hemoglobin remained almost unchanged; however, changes in the heme environment and the microenvironment of tryptophans, as well as a decrease in the protein's thermal stability, were observed. Thus, the formation of a noncovalent complex of hemoglobin with glutathione makes a more significant effect on the tertiary and quaternary structure of hemoglobin than glutathionylation, which mainly affects the secondary structure of the protein. The obtained data are important for understanding the functioning of glutathionylated hemoglobin, which is a marker of oxidative stress, and hemoglobin in complex with GSH, which appears to deposit GSH and release it during deoxygenation to increase the antioxidant protection of cells.

Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium

Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis.

PKAD-2: New entries and expansion of functionalities of the database of experimentally measured pKa's of proteins

Almost all biological reactions are pH dependent and understanding the origin of pH dependence requires knowledge of the pKa's of ionizable groups. Here we report a new edition of PKAD, the PKAD-2, which is a database of experimentally measured pKa's of proteins, both wild type and mutant proteins. The new additions include 117 wild type and 54 mutant pKa values, resulting in total 1742 experimentally measured pKa's. The new edition of PKAD-2 includes 8 new wild type and 12 new mutant proteins, resulting in total of 220 proteins. This new edition incorporates a visual 3D image of the highlighted residue of interest within the corresponding protein or protein complex. Hydrogen bonds were identified, counted, and implemented as a search feature. Other new search features include the number of neighboring residues <4A from the heaviest atom of the side chain of a given amino acid. Here, we present PKAD-2 with the intention to continuously incorporate novel features and current data with the goal to be used as benchmark for computational methods.

Discovering functionally important sites in proteins

Proteins play important roles in biology, biotechnology and pharmacology, and missense variants are a common cause of disease. Discovering functionally important sites in proteins is a central but difficult problem because of the lack of large, systematic data sets. Sequence conservation can highlight residues that are functionally important but is often convoluted with a signal for preserving structural stability. We here present a machine learning method to predict functional sites by combining statistical models for protein sequences with biophysical models of stability. We train the model using multiplexed experimental data on variant effects and validate it broadly. We show how the model can be used to discover active sites, as well as regulatory and binding sites. We illustrate the utility of the model by prospective prediction and subsequent experimental validation on the functional consequences of missense variants in HPRT1 which may cause Lesch-Nyhan syndrome, and pinpoint the molecular mechanisms by which they cause disease.

Dimerization of European Robin Cryptochrome 4a

Homo-dimer formation is important for the function of many proteins. Although dimeric forms of cryptochromes (Cry) have been found by crystallography and were recently observed in vitro for European robin Cry4a, little is known about the dimerization of avian Crys and the role it could play in the mechanism of magnetic sensing in migratory birds. Here, we present a combined experimental and computational investigation of the dimerization of robin Cry4a resulting from covalent and non-covalent interactions. Experimental studies using native mass spectrometry, mass spectrometric analysis of disulfide bonds, chemical cross-linking, and photometric measurements show that disulfide-linked dimers are routinely formed, that their formation is promoted by exposure to blue light, and that the most likely cysteines are C317 and C412. Computational modeling and molecular dynamics simulations were used to generate and assess a number of possible dimer structures. The relevance of these findings to the proposed role of Cry4a in avian magnetoreception is discussed.

Multisite Labeling of Proteins Using the Ligand-Directed Reactivity of Triggerable Michael Acceptors

Targeted modification of endogenous proteins without genetic manipulation of protein expression machinery has a range of applications from chemical biology to drug discovery. Despite being demonstrated to be effective in various applications, target-specific protein labeling using ligand-directed strategies is limited by stringent amino acid selectivity. Here, we present highly reactive ligand-directed triggerable Michael acceptors (LD-TMAcs) that feature rapid protein labeling. Unlike previous approaches, the unique reactivity of LD-TMAcs enables multiple modifications on a single target protein, effectively mapping the ligand binding site. This capability is attributed to the tunable reactivity of TMAcs that enable the labeling of several amino acid functionalities via a binding-induced increase in local concentration while remaining fully dormant in the absence of protein binding. We demonstrate the target selectivity of these molecules in cell lysates using carbonic anhydrase as the model protein. Furthermore, we demonstrate the utility of this method by selectively labeling membrane-bound carbonic anhydrase XII in live cells. We envision that the unique features of LD-TMAcs will find use in target identification, investigation of binding/allosteric sites, and studying membrane proteins.

Arsenic binding to human metallothionein-3

Arsenic poisoning is of great concern with respect to its neurological toxicity, which is especially significant for young children. Human exposure to arsenic occurs worldwide from contaminated drinking water. In human physiology, one response to toxic metals is through coordination with the metallochaperone metallothionein (MT). Central nervous system expression of MT isoform 3 (MT3) is thought to be neuroprotective. We report for the first time on the metalation pathways of As3+ binding to apo-MT3 under physiological conditions, yielding the absolute binding constants (log Kn, n = 1-6) for each sequential As3+ binding event: 10.20, 10.02, 9.79, 9.48, 9.06, and 8.31 M-1. We report on the rate of the reaction of As3+ with apo-MT3 at pH 3.5 with rate constants (kn, n = 1-6) determined for each sequential As3+ binding event: 116.9, 101.2, 85.6, 64.0, 43.9, and 21.0 M-1 s-1. We further characterize the As3+ binding pathway to fully metalated Zn7MT3 and partially metalated Zn-MT3. As3+ binds rapidly with high binding constants under physiological conditions in a noncooperative manner, but is unable to replace the Zn2+ in fully-metalated Zn-MT3. As3+ binding to partially metalated Zn-MT3 takes place with a rearrangement of the Zn-binding profile. Our work shows that As 3+ rapidly and efficiently binds to both apo-MT3 and partially metalated Zn-MT3 at physiological pH.

Deep Mutational Scanning of an Oxygen-Independent Fluorescent Protein CreiLOV for Comprehensive Profiling of Mutational and Epistatic Effects

Oxygen-independent, flavin mononucleotide-based fluorescent proteins (FbFPs) are promising alternatives to green fluorescent protein in anaerobic contexts. Deep mutational scanning performs systematic profiling of protein sequence-function relationships but has not been applied to FbFPs. Focusing on CreiLOV from Chlamydomonas reinhardtii, we created and analyzed two comprehensive mutant collections: (1) single-residue, site-saturation mutagenesis libraries covering all 118 residues; and (2) a full combinatorial metagenesis library among 20 mutations at 15 residues, where mutation and residue selection was based on single-site mutagenesis results. Notably, the second type of library is indispensable to study higher-order epistasis but underrepresented in the literature. Using optimized FACS-seq assays, 2,185 (>92.5%) out of 2,360 possible single-site mutants and 165,428 (>89.7%) out of 184,320 possible combinatorial mutants were reliably assigned with fitness values. We constructed statistical and machine-learning models to analyze the CreiLOV data set, enabling accurate fitness prediction of higher-order mutants using lower-order mutagenesis data. In addition, we successfully isolated CreiLOV variants with improved fluorescence quantum yield and thermostability. This work provides new empirical data and design rules to engineer combinatorial protein variants.

Computational analysis of missense variant CYP4F2*3 (V433M) in association with human CYP4F2 dysfunction: a functional and structural impact

BACKGROUND

Cytochrome P450 4F2 (CYP4F2) enzyme is a member of the CYP4 family responsible for the metabolism of fatty acids, therapeutic drugs, and signaling molecules such as arachidonic acid, tocopherols, and vitamin K. Several reports have demonstrated that the missense variant CYP4F2*3 (V433M) causes decreased activity of CYP4F2 and inter-individual variations in warfarin dose in different ethnic groups. However, the molecular pathogenicity mechanism of missense V433M in CYP4F2 at the atomic level has not yet been completely elucidated.

METHODS AND RESULTS

In the current study, we evaluated the effect of the V433M substitution on CYP4F2 using 14 different bioinformatics tools. Further molecular dynamics (MD) simulations were performed to assess the impact of the V433M mutation on the CYP4F2 protein structure, stability, and dynamics. In addition, molecular docking was used to illustrate the effect of V433M on its interaction with vitamin K1. Based on our results, the CYP4F2*3 variant was a damaging amino acid substitution with a destabilizing nature. The simulation results showed that missense V433M affects the dynamics and stability of CYP4F2 by reducing its compactness and stability, which means that it tends to change the overall structural conformation and flexibility of CYP4F2. The docking results showed that the CYP4F2*3 variant decreased the binding affinity between vitamin K1 and CYP4F2, which reduced the activity of CYP4F2*3 compared to native CYP4F2.

CONCLUSIONS

This study determined the molecular pathogenicity mechanism of the CYP4F2*3 variant on the human CYP4F2 protein and provided new information for understanding the structure-function relationship of CYP4F2 and other CYP4 enzymes. These findings will aid in the development of effective drugs and treatment options.

Next Generation of Ovarian Cancer Detection Using Aptamers

Ovarian cancer is among the seven most common types of cancer in women, being the most fatal gynecological tumor, due to the difficulty of detection in early stages. Aptamers are important tools to improve tumor diagnosis through the recognition of specific molecules produced by tumors. Here, aptamers and their potential targets in ovarian cancer cells were analyzed by in silico approaches. Specific aptamers were selected by the Cell-SELEX method using Caov-3 and OvCar-3 cells. The five most frequent aptamers obtained from the last round of selection were computationally modeled. The potential targets for those aptamers in cells were proposed by analyzing proteomic data available for the Caov-3 and OvCar-3 cell lines. Overexpressed proteins for each cell were characterized as to their three-dimensional model, cell location, and electrostatic potential. As a result, four specific aptamers for ovarian tumors were selected: AptaC2, AptaC4, AptaO1, and AptaO2. Potential targets were identified for each aptamer through Molecular Docking, and the best complexes were AptaC2-FXYD3, AptaC4-ALPP, AptaO1-TSPAN15, and AptaO2-TSPAN15. In addition, AptaC2 and AptaO1 could detect different stages and subtypes of ovarian cancer tissue samples. The application of this technology makes it possible to propose new molecular biomarkers for the differential diagnosis of epithelial ovarian cancer.

Plant Defense Elicitation by the Hydrophobin Cerato-Ulmin and Correlation with Its Structural Features

Cerato-ulmin (CU) is a 75-amino-acid-long protein that belongs to the hydrophobin family. It self-assembles at hydrophobic-hydrophilic interfaces, forming films that reverse the wettability properties of the bound surface: a capability that may confer selective advantages to the fungus in colonizing and infecting elm trees. Here, we show for the first time that CU can elicit a defense reaction (induction of phytoalexin synthesis and ROS production) in non-host plants (Arabidopsis) and exerts its eliciting capacity more efficiently when in its soluble monomeric form. We identified two hydrophobic clusters on the protein's loops endowed with dynamical and physical properties compatible with the possibility of reversibly interconverting between a disordered conformation and a β-strand-rich conformation when interacting with hydrophilic or hydrophobic surfaces. We propose that the plasticity of those loops may be part of the molecular mechanism that governs the protein defense elicitation capability.

Cryo-EM structure of hnRNPDL-2 fibrils, a functional amyloid associated with limb-girdle muscular dystrophy D3

hnRNPDL is a ribonucleoprotein (RNP) involved in transcription and RNA-processing that hosts missense mutations causing limb-girdle muscular dystrophy D3 (LGMD D3). Mammalian-specific alternative splicing (AS) renders three natural isoforms, hnRNPDL-2 being predominant in humans. We present the cryo-electron microscopy structure of full-length hnRNPDL-2 amyloid fibrils, which are stable, non-toxic, and bind nucleic acids. The high-resolution amyloid core consists of a single Gly/Tyr-rich and highly hydrophilic filament containing internal water channels. The RNA binding domains are located as a solenoidal coat around the core. The architecture and activity of hnRNPDL-2 fibrils are reminiscent of functional amyloids, our results suggesting that LGMD D3 might be a loss-of-function disease associated with impaired fibrillation. Strikingly, the fibril core matches exon 6, absent in the soluble hnRNPDL-3 isoform. This provides structural evidence for AS controlling hnRNPDL assembly by precisely including/skipping an amyloid exon, a mechanism that holds the potential to generate functional diversity in RNPs.

Apo-metallothionein-3 cooperatively forms tightly compact structures under physiological conditions

Metallothioneins (MTs) are essential mammalian metal chaperones. MT isoform 1 (MT1) is expressed in the kidneys and isoform 3 (MT3) is expressed in nervous tissue. For MTs, the solution-based NMR structure was determined for metal-bound MT1 and MT2, and only one X-ray diffraction structure on a crystallized mixed metal-bound MT2 has been reported. The structure of solution-based metalated MT3 is partially known using NMR methods; however, little is known about the fluxional de novo apo-MT3 because the structure cannot be determined by traditional methods. Here, we used cysteine modification coupled with electrospray ionization mass spectrometry, denaturing reactions with guanidinium chloride, stopped-flow methods measuring cysteine modification and metalation, and ion mobility mass spectrometry to reveal that apo-MT3 adopts a compact structure under physiological conditions and an extended structure under denaturing conditions, with no intermediates. Compared with apo-MT1, we found that this compact apo-MT3 binds to a cysteine modifier more cooperatively at equilibrium and 0.5 times the rate, providing quantitative evidence that many of the 20 cysteines of apo-MT3 are less accessible than those of apo-MT1. In addition, this compact apo-MT3 can be identified as a distinct population using ion mobility mass spectrometry. Furthermore, proposed structural models can be calculated using molecular dynamics methods. Collectively, these findings provide support for MT3 acting as a noninducible regulator of the nervous system compared with MT1 as an inducible scavenger of trace metals and toxic metals in the kidneys.

Structural insights on the selective interaction of the histidine-rich piscidin antimicrobial peptide Of-Pis1 with membranes

Of-Pis1 is a potent piscidin antimicrobial peptide (AMP), recently isolated from rock bream (Oplegnathus fasciatus). This rich in histidines and glycines 24-amino acid peptide displays high and broad antimicrobial activity and no significant hemolytic toxicity against human erythrocytes, suggesting low toxicity. To better understand the mechanism of action of Of-Pis1 and its potential selectivity, using NMR and CD spectroscopies, we studied the interaction with eukaryotic and procaryotic membranes and membrane models. Anionic sodium dodecyl sulfate (SDS) and lipopolysaccharide (LPS) micelles were used to mimic procaryotic membranes, while zwitterionic dodecyl phosphocholine (DPC) was used as eukaryotic membrane surrogate. In an aqueous environment, Of-Pis1 adopts a flexible random coil conformation. In DPC and SDS instead, the N-terminal region of Of-Pis1 forms an amphipathic α-helix with the non-polar face in close contact with the micelles. Slower solvent exchange and higher pK_as of the histidine residues in SDS than in DPC suggest that Of-Pis1 interacts more tightly with SDS. Of-Pis1 also binds tightly and structurally perturbs LPS micelles. Of-Pis1 interacts with both Escherichia coli and mammalian cell membranes, but only in the presence of Escherichia coli membranes it populates the helical conformation. Furthermore, ligand-based NMR experiments support a tighter and more specific interaction with bacterial than with eukaryotic membranes. Overall, these data clearly show the selective interaction of this broadly active AMP with bacterial over eukaryotic membranes. The conformational information is discussed in terms of Of-Pis1 amino acid sequence and composition to provide insights useful to design more potent and selective AMPs.

Investigating the functional role of a buried interchain aromatic cluster in Escherichia coli GrpE dimer

Aromatic clusters in the core of proteins are often involved in imparting structural stability to proteins. However, their functional importance is not always clear. In this study, we investigate the thermosensing role of a phenylalanine cluster present in the GrpE homodimer. GrpE, which acts as a nucleotide exchange factor for the molecular chaperone DnaK, is well known for its thermosensing activity resulting from temperature-dependent structural changes that allow control of chaperone function. Using mutational analysis, we show that an interchain phenylalanine cluster in a four-helix bundle of the GrpE homodimer assists in the thermosensing ability of the co-chaperone. Substitution of aromatic residues with hydrophobic ones in the core of the four-helix bundle reduces the thermal stability of the bundle and that of a connected coiled-coil domain, which impacts thermosensing. Cell growth assays and SEM images of the mutants show filamentous growth of Escherichia coli cells at 42°C, which corroborates with the defect in thermosensing. Our work suggests that the interchain edge-to-face aromatic cluster is important for the propagation of the structural signal from the coiled-coil domain to the four-helical bundle of GrpE, thus facilitating GrpE-mediated thermosensing in bacteria.

Multiple polarity kinases inhibit phase separation of F-BAR protein Cdc15 and antagonize cytokinetic ring assembly in fission yeast

The F-BAR protein Cdc15 is essential for cytokinesis in Schizosaccharomyces pombe and plays a key role in attaching the cytokinetic ring (CR) to the plasma membrane (PM). Cdc15's abilities to bind to the membrane and oligomerize via its F-BAR domain are inhibited by phosphorylation of its intrinsically disordered region (IDR). Multiple cell polarity kinases regulate Cdc15 IDR phosphostate, and of these the DYRK kinase Pom1 phosphorylation sites on Cdc15 have been shown in vivo to prevent CR formation at cell tips. Here, we compared the ability of Pom1 to control Cdc15 phosphostate and cortical localization to that of other Cdc15 kinases: Kin1, Pck1, and Shk1. We identified distinct but overlapping cohorts of Cdc15 phosphorylation sites targeted by each kinase, and the number of sites correlated with each kinases' abilities to influence Cdc15 PM localization. Coarse-grained simulations predicted that cumulative IDR phosphorylation moves the IDRs of a dimer apart and toward the F-BAR tips. Further, simulations indicated that the overall negative charge of phosphorylation masks positively charged amino acids necessary for F-BAR oligomerization and membrane interaction. Finally, simulations suggested that dephosphorylated Cdc15 undergoes phase separation driven by IDR interactions. Indeed, dephosphorylated but not phosphorylated Cdc15 undergoes liquid-liquid phase separation to form droplets in vitro that recruit Cdc15 binding partners. In cells, Cdc15 phosphomutants also formed PM-bound condensates that recruit other CR components. Together, we propose that a threshold of Cdc15 phosphorylation by assorted kinases prevents Cdc15 condensation on the PM and antagonizes CR assembly.

A novel missense variant c.71G > T (p.Gly24Val) of the CRYBA4 gene contributes to autosomal-dominant congenital cataract in a Chinese family

PURPOSE

To investigate the potential genetic defects in a five-generation Chinese family with autosomal dominant congenital cataract (ADCC).

METHODS

Whole exome sequencing was performed to search the variants in the candidate genes associated with congenital cataract. Sanger sequencing was used to validate the variants and examine their co-segregation in the patients and their relatives. The potential effect of the variants was analyzed using several bioinformatic methods and further examined through Western blotting and co-immunoprecipitation.

RESULTS

A missense variant c. 71 G > T (p. Gly24Val) in the CRYBA4 gene, a known ADCC candidate gene, was identified to be heterozygously present in the patients and co-segregate with cataract in the family. The mutation was absent in all of the searched databases, including our in-house exome sequences of 10,000 Chinese. The alignments of the amino acid sequences of CRYBA4 in a variety of species revealed that the amino acid residue Gly24 was evolutionarily highly conserved, and the in silico analysis predicted that the missense mutation of Gly24Val was damaging for the protein structure and function of CRYBA4. Then, the in vitro expression analysis further revealed that the Gly24Val mutation in CRYBA4 inhibited its binding with CRYBB1. The impaired interaction of β-crystallin proteins may affect their water-solubility and contribute to the formation of precipitates in lens fiber cells.

CONCLUSION

We identified a novel missense variant in the CRYBA4 gene as a pathogenic mutation of ADCC in a Chinese family. Our finding expanded the CRYBA4 variation spectrum associated with congenital cataracts.

Small Molecular Antimicrobial Ligands of YspD are Potential Therapeutic Agents Against Yersinia enterocolitica Infection

YspD is a hydrophilic translocator forming the platform for assemblage of functional translocon. Exposure to the extra-cellular milieu makes YspD a potential therapeutic target. DoGSiteScorer predicted best druggable pocket (P0) within YspD, encompassing predominantly the C-terminal helical bundles and the long helices-9 & 5. COACH metaserver also identified ligand binding residues within the aforementioned druggable pocket mapping to helix-9. Amino acids of helix-9 are involved in oligomerization of YspD. Interaction of helix-9 and parts of C-terminal of YspD with hydrophobic translocator protein (YspB), is essential for translocation of bacterial effectors to initiate an infection. Helices-9 & 5 form an intramolecular coiled-coil structure, required for protein-protein interaction. Targeting intramolecular coiled-coil and parts of C-terminal would be important for functional inactivation of YspD. Solvent exposed surface in YspD, particularly in P0, enhances its accessibility to ligands. Nine small molecular inhibitors of TIIISS were identified and retrieved from ZINC15 database (drug-library) as putative drug candidates. Molecular docking of potential ligands with P0 was done using SwissDock server and Achilles Blind Docking server. Considering the "Significance" threshold of binding score and region of interaction, Salicylidene Acyl Hydrazide derivatives (INP0400) and Phenoxyacetamide derivative (MBX1641) were found to bind effectively with YspD. These potential ligands interact with functional domains of YspD including parts of C-terminal and the intramolecular coiled-coil, which may affect the oligomerization of YspD and disrupt the interaction of YspD with YspB, inhibiting formation of functional translocon. The identified small molecular antimicrobial ligands of YspD could be tested in vivo to attenuate Y. enterocolitica infection by deregulation of Ysa-Ysp TIIISS.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40011-022-01443-2.