3V: cavity, channel and cleft volume calculator and extractor

The structure of inactivated mature tick-borne encephalitis virus at 3.0 Å resolution

Tick-borne encephalitis virus (TBEV) causes a severe disease, tick-borne encephalitis (TBE), that has a substantial epidemiological importance for Northern Eurasia. Between 10,000 and 15,000 TBE cases are registered annually despite the availability of effective formaldehyde-inactivated full-virion vaccines due to insufficient vaccination coverage, as well as sporadic cases of vaccine breakthrough. The development of improved vaccines would benefit from the atomic resolution structure of the antigen. Here we report the refined single-particle cryo-electron microscopy (cryo-EM) structure of the inactivated mature TBEV vaccine strain Sofjin-Chumakov (Far-Eastern subtype) at a resolution of 3.0 Å. The increase of the resolution with respect to the previously published structures of TBEV strains Hypr and Kuutsalo-14 (European subtype) was reached due to improvement of the virus sample quality achieved by the optimized preparation methods. All the surface epitopes of TBEV were structurally conserved in the inactivated virions. ELISA studies with monoclonal antibodies supported the hypothesis of TBEV protein shell cross-linking upon inactivation with formaldehyde.

GolpHCat (TMEM87A), a unique voltage-dependent cation channel in Golgi apparatus, contributes to Golgi-pH maintenance and hippocampus-dependent memory

Impaired ion channels regulating Golgi pH lead to structural alterations in the Golgi apparatus, such as fragmentation, which is found, along with cognitive impairment, in Alzheimer's disease. However, the causal relationship between altered Golgi structure and cognitive impairment remains elusive due to the lack of understanding of ion channels in the Golgi apparatus of brain cells. Here, we identify that a transmembrane protein TMEM87A, renamed Golgi-pH-regulating cation channel (GolpHCat), expressed in astrocytes and neurons that contributes to hippocampus-dependent memory. We find that GolpHCat displays unique voltage-dependent currents, which is potently inhibited by gluconate. Additionally, we gain structural insights into the ion conduction through GolpHCat at the molecular level by determining three high-resolution cryogenic-electron microscopy structures of human GolpHCat. GolpHCat-knockout mice show fragmented Golgi morphology and altered protein glycosylation and functions in the hippocampus, leading to impaired spatial memory. These findings suggest a molecular target for Golgi-related diseases and cognitive impairment.

CageCavityCalc (C3): A Computational Tool for Calculating and Visualizing Cavities in Molecular Cages

Organic(porous) and metal-organic cages are promising biomimetic platforms with diverse applications spanning recognition, sensing, and catalysis. The key to the emergence of these functions is the presence of well-defined inner cavities capable of binding a wide range of guest molecules and modulating their properties. However, despite the myriad cage architectures currently available, the rational design of structurally diverse and functional cages with specific host-guest properties remains challenging. Efficiently predicting such properties is critical for accelerating the discovery of novel functional cages. Herein, we introduce CageCavityCalc (C3), a Python-based tool for calculating the cavity size of molecular cages. The code is available on GitHub at https://github.com/VicenteMartiCentelles/CageCavityCalc. C3 utilizes a novel algorithm that enables the rapid calculation of cavity sizes for a wide range of molecular structures and porous systems. Moreover, C3 facilitates easy visualization of the computed cavity size alongside hydrophobic and electrostatic potentials, providing insights into host-guest interactions within the cage. Furthermore, the calculated cavity can be visualized using widely available visualization software, such as PyMol, VMD, or ChimeraX. To enhance user accessibility, a PyMol plugin has been created, allowing nonspecialists to use this tool without requiring computer programming expertise. We anticipate that the deployment of this computational tool will significantly streamline cage cavity calculations, thereby accelerating the discovery of functional cages.

Shy is a Proteobacterial steroid hydratase which catalyzes steroid side chain degradation without requiring a catalytically inert partner domain

Shy (side chain hydratase) and Sal (side chain aldolase), are involved in successive reactions in the pathway of bile acid side chain catabolism in Proteobacteria. Untagged Shy co-purified with His-tagged Sal indicating that the two enzymes form a complex. Shy contains a MaoC and a DUF35 domain. When co-expressed with Sal, the DUF35 domain but not the MaoC domain of Shy was observed to co-purify with Sal, indicating Sal interacts with Shy through its DUF35 domain. The MaoC domain of Shy (ShyMaoC) remained catalytically viable and could hydrate cholyl-enoyl-CoA with similar catalytic efficiency as in the Shy-Sal complex. Sal expressed with the DUF35 domain of Shy (Sal-ShyDUF35) was similarly competent for the retroaldol cleavage of cholyl-3-OH-CoA. ShyMaoC showed a preference for C5 side chain bile acid substrates, exhibiting low activity towards C3 side chain substrates. The ShyMaoC structure was determined by X-ray crystallography, showing a hot dog fold with a short central helix surrounded by a twisted anti-parallel β-sheet. Modeling and mutagenesis studies suggest that the bile acid substrate occupies the large open cleft formed by the truncated central helix and repositioning of the active site housing. ShyMaoC therefore contains two substrate binding sites per homodimer, making it distinct from previously characterized MaoC steroid hydratases that are (pseudo)-heterodimers with one substrate binding site per dimer. The characterization of Shy provides insight into how MaoC family hydratases have adapted to accommodate large polycyclic substrates that can facilitate future engineering of these enzymes to produce novel steroid pharmaceuticals.

Cryo-EM structure of cadmium-bound human ABCB6

ATP-binding cassette transporter B6 (ABCB6), a protein essential for heme biosynthesis in mitochondria, also functions as a heavy metal efflux pump. Here, we present cryo-electron microscopy structures of human ABCB6 bound to a cadmium Cd(II) ion in the presence of antioxidant thiol peptides glutathione (GSH) and phytochelatin 2 (PC2) at resolutions of 3.2 and 3.1 Å, respectively. The overall folding of the two structures resembles the inward-facing apo state but with less separation between the two halves of the transporter. Two GSH molecules are symmetrically bound to the Cd(II) ion in a bent conformation, with the central cysteine protruding towards the metal. The N-terminal glutamate and C-terminal glycine of GSH do not directly interact with Cd(II) but contribute to neutralizing positive charges of the binding cavity by forming hydrogen bonds and van der Waals interactions with nearby residues. In the presence of PC2, Cd(II) binding to ABCB6 is similar to that observed with GSH, except that two cysteine residues of each PC2 molecule participate in Cd(II) coordination to form a tetrathiolate. Structural comparison of human ABCB6 and its homologous Atm-type transporters indicate that their distinct substrate specificity might be attributed to variations in the capping residues situated at the top of the substrate-binding cavity.

DrugMap: A quantitative pan-cancer analysis of cysteine ligandability

Cysteine-focused chemical proteomic platforms have accelerated the clinical development of covalent inhibitors for a wide range of targets in cancer. However, how different oncogenic contexts influence cysteine targeting remains unknown. To address this question, we have developed "DrugMap," an atlas of cysteine ligandability compiled across 416 cancer cell lines. We unexpectedly find that cysteine ligandability varies across cancer cell lines, and we attribute this to differences in cellular redox states, protein conformational changes, and genetic mutations. Leveraging these findings, we identify actionable cysteines in NF-κB1 and SOX10 and develop corresponding covalent ligands that block the activity of these transcription factors. We demonstrate that the NF-κB1 probe blocks DNA binding, whereas the SOX10 ligand increases SOX10-SOX10 interactions and disrupts melanoma transcriptional signaling. Our findings reveal heterogeneity in cysteine ligandability across cancers, pinpoint cell-intrinsic features driving cysteine targeting, and illustrate the use of covalent probes to disrupt oncogenic transcription-factor activity.

Spatial arrangement and density variations in the cell envelope of Deinococcus radiodurans

The cell envelope of the poly-extremophile bacterium Deinococcus radiodurans is renowned for its highly organized structure and unique functional characteristics. In this bacterium, a precise regularity characterizes not just the S-layer, but it also extends to the underlying cell envelope layers, resulting in a dense and tightly arranged configuration. This regularity is attributed to a minimum of three protein complexes located at the outer membrane level. Together, they constitute a recurring structural unit that extends across the cell envelope, effectively tiling the entirety of the cell body. Nevertheless, a comprehensive grasp of the vacant spaces within each layer and their functional roles remains limited. In this study, we delve into these aspects by integrating the state of the art with structural calculations. This approach provides crucial evidence supporting an evolutive pressure intricately linked to surface phenomena depending on the environmental conditions.

Identification of an inhibitory pocket in falcilysin provides a new avenue for malaria drug development

Identification of new druggable protein targets remains the key challenge in the current antimalarial development efforts. Here we used mass-spectrometry-based cellular thermal shift assay (MS-CETSA) to identify potential targets of several antimalarials and drug candidates. We found that falcilysin (FLN) is a common binding partner for several drug candidates such as MK-4815, MMV000848, and MMV665806 but also interacts with quinoline drugs such as chloroquine and mefloquine. Enzymatic assays showed that these compounds can inhibit FLN proteolytic activity. Their interaction with FLN was explored systematically by isothermal titration calorimetry and X-ray crystallography, revealing a shared hydrophobic pocket in the catalytic chamber of the enzyme. Characterization of transgenic cell lines with lowered FLN expression demonstrated statistically significant increases in susceptibility toward MK-4815, MMV000848, and several quinolines. Importantly, the hydrophobic pocket of FLN appears amenable to inhibition and the structures reported here can guide the development of novel drugs against malaria.

Studying the Effects of Dissolved Noble Gases and High Hydrostatic Pressure on the Spherical DOPC Bilayer Using Molecular Dynamic Simulations

Fine-grained molecular dynamics simulations have been conducted to depict lipid objects enclosed in water and interacting with a series of noble gases dissolved in the medium. The simple point-charge (SPC) water system, featuring a boundary composed of 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC) molecules, maintained stability throughout the simulation under standard conditions. This allowed for the accurate modeling of the effects of hydrostatic pressure at an ambient pressure of 25 bar. The chosen pressure references the 240 m depth of seawater: the horizon frequently used by commercial divers, who comprise the primary patient population of the neurological complication of inert gas narcosis and the consequences of high-pressure neurological syndrome. To quantify and validate the neurological effects of noble gases and discriminate them from high hydrostatic pressure, we reduced the dissolved gas molar concentration to 1.5%, three times smaller than what we previously tested for the planar bilayer (3.5%). The nucleation and growth of xenon, argon and neon nanobubbles proved consistent with the data from the planar bilayer simulations. On the other hand, hyperbaric helium induces only a residual distorting effect on the liposome, with no significant condensed gas fraction observed within the hydrophobic core. The bubbles were distributed over a large volume-both in the bulk solvent and in the lipid phase-thereby causing substantial membrane distortion. This finding serves as evidence of the validity of the multisite distortion hypothesis for the neurological effect of inert gases at high pressure.

Chaperone-assisted cryo-EM structure of P. aeruginosa PhuR reveals molecular basis for heme binding

Pathogenic bacteria, such as Pseudomonas aeruginosa, depend on scavenging heme for the acquisition of iron, an essential nutrient. The TonB-dependent transporter (TBDT) PhuR is the major heme uptake protein in P. aeruginosa clinical isolates. However, a comprehensive understanding of heme recognition and TBDT transport mechanisms, especially PhuR, remains limited. In this study, we employed single-particle cryogenic electron microscopy (cryo-EM) and a phage display-generated synthetic antibody (sAB) as a fiducial marker to enable the determination of a high-resolution (2.5 Å) structure of PhuR with a bound heme. Notably, the structure reveals iron coordination by Y529 on a conserved extracellular loop, shedding light on the role of tyrosine in heme binding. Biochemical assays and negative-stain EM demonstrated that the sAB specifically targets the heme-bound state of PhuR. These findings provide insights into PhuR's heme binding and offer a template for developing conformation-specific sABs against outer membrane proteins (OMPs) for structure-function investigations.

Small angle X-ray scattering analysis of thermophilic cytochrome P450 CYP119 and the effects of the N-terminal histidine tag

Combining size exclusion chromatography-small angle X-ray scattering (SEC-SAXS) and molecular dynamics (MD) analysis is a promising approach to investigate protein behavior in solution, particularly for understanding conformational changes due to substrate binding in cytochrome P450s (CYPs). This study investigates conformational changes in CYP119, a thermophilic CYP from Sulfolobus acidocaldarius that exhibits structural flexibility similar to mammalian CYPs. Although the crystal structure of ligand-free (open state) and ligand-bound (closed state) forms of CYP119 is known, the overall structure of the enzyme in solution has not been explored until now. It was found that theoretical scattering profiles from the crystal structures of CYP119 did not align with the SAXS data, but conformers from MD simulations, particularly starting from the open state (46 % of all frames), agreed well. Interestingly, a small percentage of closed-state conformers also fit the data (9 %), suggesting ligand-free CYP119 samples ligand-bound conformations. Ab initio SAXS models for N-His tagged CYP119 revealed a tail-like unfolded structure impacting protein flexibility, which was confirmed by in silico modeling. SEC-SAXS analysis of N-His CYP119 indicated pentameric structures in addition to monomers in solution, affecting the stability and activity of the enzyme. This study adds insights into the conformational dynamics of CYP119 in solution.

Structural basis of prostaglandin efflux by MRP4

Multidrug resistance protein 4 (MRP4) is a broadly expressed ATP-binding cassette transporter that is unique among the MRP subfamily for transporting prostanoids, a group of signaling molecules derived from unsaturated fatty acids. To better understand the basis of the substrate selectivity of MRP4, we used cryogenic-electron microscopy to determine six structures of nanodisc-reconstituted MRP4 at various stages throughout its transport cycle. Substrate-bound structures of MRP4 in complex with PGE1, PGE2 and the sulfonated-sterol DHEA-S reveal a common binding site that accommodates a diverse set of organic anions and suggest an allosteric mechanism for substrate-induced enhancement of MRP4 ATPase activity. Our structure of a catalytically compromised MRP4 mutant bound to ATP-Mg2+ is outward-occluded, a conformation previously unobserved in the MRP subfamily and consistent with an alternating-access transport mechanism. Our study provides insights into the endogenous function of this versatile efflux transporter and establishes a basis for MRP4-targeted drug design.

Application of Machine Learning in the Quantitative Analysis of the Surface Characteristics of Highly Abundant Cytoplasmic Proteins: Toward AI-Based Biomimetics

Proteins in the crowded environment of human cells have often been studied regarding nonspecific interactions, misfolding, and aggregation, which may cause cellular malfunction and disease. Specifically, proteins with high abundance are more susceptible to these issues due to the law of mass action. Therefore, the surfaces of highly abundant cytoplasmic (HAC) proteins directly exposed to the environment can exhibit specific physicochemical, structural, and geometrical characteristics that reduce nonspecific interactions and adapt to the environment. However, the quantitative relationships between the overall surface descriptors still need clarification. Here, we used machine learning to identify HAC proteins using hydrophobicity, charge, roughness, secondary structures, and B-factor from the protein surfaces and quantified the contribution of each descriptor. First, several supervised learning algorithms were compared to solve binary classification problems for the surfaces of HAC and extracellular proteins. Then, logistic regression was used for the feature importance analysis of descriptors considering model performance (80.2% accuracy and 87.6% AUC) and interpretability. The HAC proteins showed positive correlations with negatively and positively charged areas but negative correlations with hydrophobicity, the B-factor, the proportion of beta structures, roughness, and the proportion of disordered regions. Finally, the details of each descriptor could be explained concerning adaptative surface strategies of HAC proteins to regulate nonspecific interactions, protein folding, flexibility, stability, and adsorption. This study presented a novel approach using various surface descriptors to identify HAC proteins and provided quantitative design rules for the surfaces well-suited to human cellular crowded environments.

Discovery and structural characterization of a thermostable bacterial monoamine oxidase

Monoamine oxidases (MAOs) are pivotal regulators of neurotransmitters in mammals, while microbial MAOs have been shown to be valuable biocatalysts for enantioselective synthesis of pharmaceutical compounds or precursors thereof. To extend the knowledge of how MAOs function at the molecular level and in order to provide more biocatalytic tools, we set out to identify and study a robust bacterial variant: a MAO from the thermophile Thermoanaerobacterales bacterium (MAOTb ). MAOTb is highly thermostable with melting temperatures above 73 °C and is well expressed in Escherichia coli. Substrate screening revealed that the oxidase is most efficient with n-alkylamines with n-heptylamine being the best substrate. Presteady-state kinetic analysis shows that reduced MAOTb rapidly reacts with molecular oxygen, confirming that it is a bona fide oxidase. The crystal structure of MAOTb was resolved at 1.5 Å and showed an exceptionally high similarity with the two human MAOs, MAO A and MAO B. The active site of MAOTb resembles mostly the architecture of human MAO A, including the cysteinyl protein-FAD linkage. Yet, the bacterial MAO lacks a C-terminal extension found in human MAOs, which explains why it is expressed and purified as a soluble protein, while the mammalian counterparts are anchored to the membrane through an α-helix. MAOTb also displays a slightly different active site access tunnel, which may explain the specificity toward long aliphatic amines. Being an easy-to-express, thermostable enzyme, for which a high-resolution structure was elucidated, this bacterial MAO may develop into a valuable biocatalyst for synthetic chemistry or biosensing.

Charge neutralization and β-elimination cleavage mechanism of family 42 L-rhamnose-α-1,4-D-glucuronate lyase revealed using neutron crystallography

Gum arabic (GA) is widely used as an emulsion stabilizer and edible coating and consists of a complex carbohydrate moiety with a rhamnosyl-glucuronate group capping the non-reducing ends. Enzymes that can specifically cleave the glycosidic chains of GA and modify their properties are valuable for structural analysis and industrial application. Cryogenic X-ray crystal structure of GA-specific L-rhamnose-α-1,4-D-glucuronate lyase from Fusarium oxysporum (FoRham1), belonging to the polysaccharide lyase (PL) family 42, has been previously reported. To determine the specific reaction mechanism based on its hydrogen-containing enzyme structure, we performed joint X-ray/neutron crystallography of FoRham1. Large crystals were grown in the presence of L-rhamnose (a reaction product), and neutron and X-ray diffraction datasets were collected at room temperature at 1.80 and 1.25 Å resolutions, respectively. The active site contained L-rhamnose and acetate, the latter being a partial analog of glucuronate. Incomplete H/D exchange between Arg166 and acetate suggested that a strong salt-bridge interaction was maintained. Doubly deuterated His105 and deuterated Tyr150 supported the interaction between Arg166 and the acetate. The unique hydrogen-rich environment functions as a charge neutralizer for glucuronate and stabilizes the oxyanion intermediate. The NE2 atom of His85 was deprotonated and formed a hydrogen bond with the deuterated O1 hydroxy of L-rhamnose, indicating the function of His85 as the base/acid catalyst for bond cleavage via β-elimination. Asp83 functions as a pivot between the two catalytic histidine residues by bridging them. This His-His-Asp structural motif is conserved in the PL 24, 25, and 42 families.

Crystal structure and solution scattering of Geobacillus stearothermophilus S9 peptidase reveal structural adaptations for carboxypeptidase activity

Acylaminoacyl peptidases (AAPs) play a pivotal role in various pathological conditions and are recognized as potential therapeutic targets. AAPs exhibit a wide range of activities, such as acylated amino acid-dependent aminopeptidase, endopeptidase, and less studied carboxypeptidase activity. We have determined the crystal structure of an AAP from Geobacillus stearothermophilus (S9gs) at 2.0 Å resolution. Despite being annotated as an aminopeptidase in the NCBI database, our enzymatic characterization proved S9gs to be a carboxypeptidase. Solution-scattering studies showed that S9gs exists as a tetramer in solution, and crystal structure analysis revealed adaptations responsible for the carboxypeptidase activity of S9gs. The findings present a hypothesis for substrate selection, substrate entry, and product exit from the active site, enriching our understanding of this rare carboxypeptidase.

Dynamic PRC1-CBX8 stabilizes a porous structure of chromatin condensates

The compaction of chromatin is a prevalent paradigm in gene repression. Chromatin compaction is commonly thought to repress transcription by restricting chromatin accessibility. However, the spatial organisation and dynamics of chromatin compacted by gene-repressing factors are unknown. Using cryo-electron tomography, we solved the threedimensional structure of chromatin condensed by the Polycomb Repressive Complex 1 (PRC1) in a complex with CBX8. PRC1-condensed chromatin is porous and stabilised through multivalent dynamic interactions of PRC1 with chromatin. Mechanistically, positively charged residues on the internally disordered regions (IDRs) of CBX8 mask negative charges on the DNA to stabilize the condensed state of chromatin. Within condensates, PRC1 remains dynamic while maintaining a static chromatin structure. In differentiated mouse embryonic stem cells, CBX8-bound chromatin remains accessible. These findings challenge the idea of rigidly compacted polycomb domains and instead provides a mechanistic framework for dynamic and accessible PRC1-chromatin condensates.

Tunnel engineering of gas-converting enzymes for inhibitor retardation and substrate acceleration

Carbon monoxide dehydrogenase (CODH), formate dehydrogenase (FDH), hydrogenase (H2ase), and nitrogenase (N2ase) are crucial enzymatic catalysts that facilitate the conversion of industrially significant gases such as CO, CO2, H2, and N2. The tunnels in the gas-converting enzymes serve as conduits for these low molecular weight gases to access deeply buried catalytic sites. The identification of the substrate tunnels is imperative for comprehending the substrate selectivity mechanism underlying these gas-converting enzymes. This knowledge also holds substantial value for industrial applications, particularly in addressing the challenges associated with separation and utilization of byproduct gases. In this comprehensive review, we delve into the emerging field of tunnel engineering, presenting a range of approaches and analyses. Additionally, we propose methodologies for the systematic design of enzymes, with the ultimate goal of advancing protein engineering strategies.

Strain-Induced Heteromorphosis Multi-Cavity Cages: Tension-Driven Self-Expansion Strategy for Controllable Enhancement of Complexity in Supramolecular Assembly

Coordinative supramolecular cages with adjustable cavities have found extensive applications in various fields, but the cavity modification strategies for multi-functional structures are still challenging. Here, we present a tension-driven self-expansion strategy for construction of multi-cavity cages with high structural complexity. Under the regulation of strain-induced capping ligands, unprecedented heteromorphosis triple-cavity cages S2 /S4 were obtained based on a metallo-organic ligand (MOL) scaffold. The heteromorphosis cages exhibited significant higher cavity diversity than the homomorphous double-cavity cages S1 /S3 ; all of the cages were thoroughly characterized through various analytical techniques including (1D and 2D) NMR, ESI-MS, TWIM-MS, AFM, and SAXS analyses. Furthermore, the encapsulation of porphyrin in the cavities of these multi-cavity cages were investigated. This research opens up new possibilities for the architecture of heteromorphosis supramolecular cages via precisely controlled "scaffold-capping" assembly with preorganized ligands, which could have potential applications in the development of multifunctional structures with higher complexity.

Polyoxometalate-Organic Hybrid "Calixarenes" as Supramolecular Hosts

Calixarenes are among the most useful and versatile macrocycles in supramolecular chemistry. The one thing that has not changed in the 80 years since their discovery, despite numerous derivatizations, is their fully organic, covalent scaffolds. Here, we report a new type of organic-inorganic hybrid "calixarenes" constructed by means of coordination-driven assembly. Replacing acetate ligands on the {SiW10 Cr2 (OAc)2 } clusters with 5-hydroxyisophthalates allows these 95° inorganic building blocks to be linked into bowl-shaped, hybrid "calix[n]arenes" (n=3, 4). With a large concave cavity, the metal-organic calix[4]arene can accommodate nanometer-sized polyoxoanions in an entropically driven process. The development of hybrid variants of calixarenes is expected to expand the scope of their physicochemical properties, guest/substrate binding, and applications on multiple fronts.

Molecular Architecture of Salmonella Typhimurium Virus P22 Genome Ejection Machinery

Bacteriophage P22 is a prototypical member of the Podoviridae superfamily. Since its discovery in 1952, P22 has become a paradigm for phage transduction and a model for icosahedral viral capsid assembly. Here, we describe the complete architecture of the P22 tail apparatus (gp1, gp4, gp10, gp9, and gp26) and the potential location and organization of P22 ejection proteins (gp7, gp20, and gp16), determined using cryo-EM localized reconstruction, genetic knockouts, and biochemical analysis. We found that the tail apparatus exists in two equivalent conformations, rotated by ∼6° relative to the capsid. Portal protomers make unique contacts with coat subunits in both conformations, explaining the 12:5 symmetry mismatch. The tail assembles around the hexameric tail hub (gp10), which folds into an interrupted β-propeller characterized by an apical insertion domain. The tail hub connects proximally to the dodecameric portal protein and head-to-tail adapter (gp4), distally to the trimeric tail needle (gp26), and laterally to six trimeric tailspikes (gp9) that attach asymmetrically to gp10 insertion domain. Cryo-EM analysis of P22 mutants lacking the ejection proteins gp7 or gp20 and biochemical analysis of purified recombinant proteins suggest that gp7 and gp20 form a molecular complex associated with the tail apparatus via the portal protein barrel. We identified a putative signal transduction pathway from the tailspike to the tail needle, mediated by three flexible loops in the tail hub, that explains how lipopolysaccharide (LPS) is sufficient to trigger the ejection of the P22 DNA in vitro.

Mitochondrial Acyl Carrier Protein of Leishmania major Displays Features Distinct from the Canonical Type II ACP

Prokaryotes synthesize fatty acids using a type II synthesis pathway (FAS). In this process, the central player, i.e., the acyl carrier protein (ACP), sequesters the growing acyl chain in its internal hydrophobic cavity. As the acyl chain length increases, the cavity expands in size, which is reflected in the NMR chemical shift perturbations and crystal structures of the acyl-ACP intermediates. A few eukaryotic organelles, such as plastids and mitochondria, also harbor type II fatty acid synthesis machinery. Plastid FAS from spinach and Plasmodium falciparum has been characterized at the molecular level, but the mitochondrial pathway remains unexplored. Here, we report NMR studies of the mitochondrial acyl-acyl carrier protein intermediates of Leishmania major (acyl-LmACP). Our studies show that LmACP experiences remarkably small conformational changes upon acylation, with perturbations confined to helices II and III only. CastP determined that the cavity size of apo-LmACP (PDB entry 5ZWT) is less than that of Escherichia coli ACP (PDB 1T8K). Thus, the small chemical shift perturbations observed in the LmACP intermediates, coupled with CastP results, suggest an unusually small cavity when fully expanded. The faster rate of C8-LmACP chain hydrolysis compared to E. coli ACP (EcACP) also supports these convictions. Structure comparison of LmACP with other type II ACP disclosed unique differences in the helix I and loop I conformations, as well as several residues present there. Numerous hydrophobic residues in helix I and loop I (conserved in all mitochondrial ACPs) are substituted with hydrophilic residues in the bacterial/plastid type II ACP. For instance, Phe and leucine at positions 14 and 34 in LmACP are substituted with a hydrophilic residue and Ala in bacterial/plastid type II ACP. Mutation of Leu 34 to Ala (corresponding residue in EcACP) resulted in a complete loss of structure, underscoring its importance in maintaining the ACP fold. Thus, our NMR studies, combined with insights from the crystal structure, highlight several unique features of LmACP, distinct from the prokaryote and plastid type II ACP. Given the high sequence identity, the features might be conserved in all mitochondrial ACPs.

Self-Assembly of Cluster-Mediated 3D Catenanes with Size-Specific Recognition Behavior

Although interlocked three-dimensional molecules display unique properties associated with their spatial structures, their synthesis and study of their host-guest properties remain challenging. We report the formation of a novel [2]catenane, [Et4N]@[(Tp*WS3Cu3Cl)2(cis-bpype)3]2(OTf)5 ([Et4N][1](OTf)5), by self-assembly of the cluster node [Tp*WS3Cu3Cl]+ and the organic linker (Z)-1,2-diphenyl-1,2-bis(4-(pyridin-4-yl)phenyl)ethene (cis-bpype). Single-crystal X-ray and NMR analyses established that [1]4+ is formed by the interpenetration of two cluster-organic cages. Unique cation-in-cation host-guest complexes were observed with this catenane. The crystalline, empty catenane was formed by taking advantage of the electrostatic repulsion-induced weak binding of the host. Encapsulation experiments also reveal that the empty catenane can adaptively encapsulate cations such as [Et4N]+ and [Pr4N]+ in the cross cavity but is unable to encapsulate [Bu4N]+ and [Me4N]+, although the size of the latter is compatible with that of the cavity. Theoretical calculations and volume analysis allow to unravel the ingenious role of catenane structures and the interplay between electrostatic repulsion and attractive noncovalent interactions for size-specific recognition behavior in host-guest systems involving species with similar electric charges.

Structural basis of peptide secretion for Quorum sensing by ComA

Quorum sensing (QS) is a crucial regulatory mechanism controlling bacterial signalling and holds promise for novel therapies against antimicrobial resistance. In Gram-positive bacteria, such as Streptococcus pneumoniae, ComA is a conserved efflux pump responsible for the maturation and secretion of peptide signals, including the competence-stimulating peptide (CSP), yet its structure and function remain unclear. Here, we functionally characterize ComA as an ABC transporter with high ATP affinity and determined its cryo-EM structures in the presence or absence of CSP or nucleotides. Our findings reveal a network of strong electrostatic interactions unique to ComA at the intracellular gate, a putative binding pocket for two CSP molecules, and negatively charged residues facilitating CSP translocation. Mutations of these residues affect ComA's peptidase activity in-vitro and prevent CSP export in-vivo. We demonstrate that ATP-Mg2+ triggers the outward-facing conformation of ComA for CSP release, rather than ATP alone. Our study provides molecular insights into the QS signal peptide secretion, highlighting potential targets for QS-targeting drugs.

Dynamic geometry design of cyclic peptide architectures for RNA structure

Designing inhibitors for RNA is still challenging due to the bottleneck of maintaining the binding interaction of inhibitor-RNA accompanied by subtle RNA flexibility. Thus, the current approach usually needs to screen thousands of candidate inhibitors for binding. Here, we propose a dynamic geometry design approach to enrich the hits with only a tiny pool of designed geometrically compatible scaffold candidates. First, our method uses graph-based tree decomposition to explore the complementarity rigid binding cyclic peptide and design the amino acid side chain length and charge to fit the RNA pocket. Then, we perform an energy-based dynamical network algorithm to optimize the inhibitor-RNA hydrogen bonds. Dynamic geometry-guided design yields successful inhibitors with low micromolar binding affinity scaffolds and experimentally competes with the natural RNA chaperone. The results indicate that the dynamic geometry method yields higher efficiency and accuracy than traditional methods. The strategy could be further optimized to design the length and chirality by adopting nonstandard amino acids and facilitating RNA engineering for biological or medical applications.

DrugMap: A quantitative pan-cancer analysis of cysteine ligandability

Cysteine-focused chemical proteomic platforms have accelerated the clinical development of covalent inhibitors of a wide-range of targets in cancer. However, how different oncogenic contexts influence cysteine targeting remains unknown. To address this question, we have developed DrugMap , an atlas of cysteine ligandability compiled across 416 cancer cell lines. We unexpectedly find that cysteine ligandability varies across cancer cell lines, and we attribute this to differences in cellular redox states, protein conformational changes, and genetic mutations. Leveraging these findings, we identify actionable cysteines in NFκB1 and SOX10 and develop corresponding covalent ligands that block the activity of these transcription factors. We demonstrate that the NFκB1 probe blocks DNA binding, whereas the SOX10 ligand increases SOX10-SOX10 interactions and disrupts melanoma transcriptional signaling. Our findings reveal heterogeneity in cysteine ligandability across cancers, pinpoint cell-intrinsic features driving cysteine targeting, and illustrate the use of covalent probes to disrupt oncogenic transcription factor activity.

Specifying conformational heterogeneity of multi-domain proteins at atomic resolution

The conformational landscape of multi-domain proteins is inherently linked to their specific functions. This also holds for polyubiquitin chains that are assembled by two or more ubiquitin domains connected by a flexible linker thus showing a large interdomain mobility. However, molecular recognition and signal transduction are associated with particular conformational substates that are populated in solution. Here, we apply high-resolution NMR spectroscopy in combination with dual-scale MD simulations to explore the conformational space of K6-, K29-, and K33-linked diubiquitin molecules. The conformational ensembles are evaluated utilizing a paramagnetic cosolute reporting on solvent exposure plus a set of complementary NMR parameters. This approach unravels a conformational heterogeneity of diubiquitins and explains the diversity of structural models that have been determined for K6-, K29-, and K33-linked diubiquitins in free and ligand-bound states so far. We propose a general application of the approach developed here to demystify multi-domain proteins occurring in nature.

TMH Stab-pred: Predicting the stability of α-helical membrane proteins using sequence and structural features

The folding and stability of transmembrane proteins (TMPs) are governed by the insertion of secondary structural elements into the cell membrane followed by their assembly. Understanding the important features that dictate the stability of TMPs is important for elucidating their functions. In this work, we related sequence and structure-based parameters with free energy (ΔG0) of α-helical membrane proteins. Our results showed that the free energy transfer of hydrophobic peptides, relative contact order, total interaction energy, number of hydrogen bonds and lipid accessibility of transmembrane regions are important for stability. Further, we have developed multiple-regression models to predict the stability of α-helical membrane proteins using these features and our method can predict the stability with a correlation and mean absolute error (MAE) of 0.89 and 1.21 kcal/mol, respectively, on jack-knife test. The method was validated with a blind test set of three recently reported experimental ΔG0, which could predict the stability within an average MAE of 0.51 kcal/mol. Further, we developed a webserver for predicting the stability and it is freely available at (https://web.iitm.ac.in/bioinfo2/TMHS/). The importance of selected parameters and limitations are discussed.

Intermolecular hydrogen bonding in calix[5]arene derived cavitands regulates the molecular recognition of fullerenes

We present a readily available calix[5]arene derived cavitand receptor that is stabilized in the closed cone conformer through intermolecular hydrogen bonding with methanol molecules. The receptor features a highly spherical aromatic surface that binds C60 and C70 fullerenes effectively, and the binding event can be regulated allosterically by the addition of methanol.

Structural basis of peptidoglycan synthesis by E. coli RodA-PBP2 complex

Peptidoglycan (PG) is an essential structural component of the bacterial cell wall that is synthetized during cell division and elongation. PG forms an extracellular polymer crucial for cellular viability, the synthesis of which is the target of many antibiotics. PG assembly requires a glycosyltransferase (GT) to generate a glycan polymer using a Lipid II substrate, which is then crosslinked to the existing PG via a transpeptidase (TP) reaction. A Shape, Elongation, Division and Sporulation (SEDS) GT enzyme and a Class B Penicillin Binding Protein (PBP) form the core of the multi-protein complex required for PG assembly. Here we used single particle cryo-electron microscopy to determine the structure of a cell elongation-specific E. coli RodA-PBP2 complex. We combine this information with biochemical, genetic, spectroscopic, and computational analyses to identify the Lipid II binding sites and propose a mechanism for Lipid II polymerization. Our data suggest a hypothesis for the movement of the glycan strand from the Lipid II polymerization site of RodA towards the TP site of PBP2, functionally linking these two central enzymatic activities required for cell wall peptidoglycan biosynthesis.

Crystal Structure of the Catalytic Domain of a Botulinum Neurotoxin Homologue from Enterococcus faecium: Potential Insights into Substrate Recognition

Clostridium botulinum neurotoxins (BoNTs) are the most potent toxins known, causing the deadly disease botulism. They function through Zn2+-dependent endopeptidase cleavage of SNARE (soluble N-ethylmaleimide-sensitive factor attachment protein receptor) proteins, preventing vesicular fusion and subsequent neurotransmitter release from motor neurons. Several serotypes of BoNTs produced by Clostridium botulinum (BoNT/A-/G and/X) have been well-characterised over the years. However, a BoNT-like gene (homologue of BoNT) was recently identified in the non-clostridial species, Enterococcus faecium, which is the leading cause of hospital-acquired multi-drug resistant infections. Here, we report the crystal structure of the catalytic domain of a BoNT homologue from Enterococcus faecium (LC/En) at 2.0 Å resolution. Detailed structural analysis in comparison with the full-length BoNT/En AlphaFold2-predicted structure, LC/A (from BoNT/A), and LC/F (from BoNT/F) revealed putative subsites and exosites (including loops 1-5) involved in recognition of LC/En substrates. LC/En also appears to possess a conserved autoproteolytic cleavage site whose function is yet to be established.

Structural modeling of cytochrome P450 51 from a deep-sea fish points to a novel structural feature in other CYP51s

Cytochromes P450 (CYP), enzymes involved in the metabolism of endogenous and xenobiotic substrates, provide an excellent model system to study how membrane proteins with unique functions have catalytically adapted through evolution. Molecular adaptation of deep-sea proteins to high hydrostatic pressure remains poorly understood. Herein, we have characterized recombinant cytochrome P450 sterol 14α-demethylase (CYP51), an essential enzyme of cholesterol biosynthesis, from an abyssal fish species, Coryphaenoides armatus. C. armatus CYP51 was heterologously expressed in Escherichia coli following N-terminal truncation and purified to homogeneity. Recombinant C. armatus CYP51 bound its sterol substrate lanosterol giving a Type I binding spectra (KD 15 μM) and catalyzed lanosterol 14α-demethylation turnover at 5.8 nmol/min/nmol P450. C. armatus CYP51 also bound the azole antifungals ketoconazole (KD 0.12 μM) and propiconazole (KD 0.54 μM) as determined by Type II absorbance spectra. Comparison of C. armatus CYP51 primary sequence and modeled structures with other CYP51s identified amino acid substitutions that may confer an ability to function under pressures of the deep sea and revealed heretofore undescribed internal cavities in human and other non-deep sea CYP51s. The functional significance of these cavities is not known. PROLOGUE: This paper is dedicated in memory of Michael Waterman and Tsuneo Omura, who as good friends and colleagues enriched our lives. They continue to inspire us.

Chaperone-assisted cryo-EM structure of P. aeruginosa PhuR reveals molecular basis for heme uptake

Pathogenic bacteria, such as Pseudomonas aeruginosa, depend on scavenging heme for the acquisition of iron, an essential nutrient. The TonB-dependent transporter (TBDT) PhuR is the major heme uptake protein in P. aeruginosa clinical isolates. However, a comprehensive understanding of heme recognition and TBDT transport mechanisms, especially PhuR, remains limited. In this study, we employed single-particle cryogenic electron microscopy (cryo-EM) and a phage display-generated synthetic antibody (sAB) as a fiducial marker to enable the determination of a high-resolution (2.5 Å) structure of PhuR with a bound heme. Notably, the structure reveals iron coordination by Y529 on a conserved extracellular loop, shedding light on the role of tyrosine in heme binding. Biochemical assays and negative-stain EM demonstrated that the sAB specifically targets the heme-bound state of PhuR. These findings provide insights into PhuR's heme binding and offer a template for developing conformation-specific sABs against outer membrane proteins (OMPs) for structure-function investigations.

New Crystal Form of Human Neuropilin-1 b1 Fragment with Six Electrostatic Mutations Complexed with KDKPPR Peptide Ligand

Neuropilin 1 (NRP1), a cell-surface co-receptor of a number of growth factors and other signaling molecules, has long been the focus of attention due to its association with the development and the progression of several types of cancer. For example, the KDKPPR peptide has recently been combined with a photosensitizer and a contrast agent to bind NRP1 for the detection and treatment by photodynamic therapy of glioblastoma, an aggressive brain cancer. The main therapeutic target is a pocket of the fragment b1 of NRP1 (NRP1-b1), in which vascular endothelial growth factors (VEGFs) bind. In the crystal packing of native human NRP1-b1, the VEGF-binding site is obstructed by a crystallographic symmetry neighbor protein, which prevents the binding of ligands. Six charged amino acids located at the protein surface were mutated to allow the protein to form a new crystal packing. The structure of the mutated fragment b1 complexed with the KDKPPR peptide was determined by X-ray crystallography. The variant crystallized in a new crystal form with the VEGF-binding cleft exposed to the solvent and, as expected, filled by the C-terminal moiety of the peptide. The atomic interactions were analyzed using new approaches based on a multipolar electron density model. Among other things, these methods indicated the role played by Asp320 and Glu348 in the electrostatic steering of the ligand in its binding site. Molecular dynamics simulations were carried out to further analyze the peptide binding and motion of the wild-type and mutant proteins. The simulations revealed that specific loops interacting with the peptide exhibited mobility in both the unbound and bound forms.

Structural and photophysical characterization of the small ultra-red fluorescent protein

The small Ultra-Red Fluorescent Protein (smURFP) represents a new class of fluorescent protein with exceptional photostability and brightness derived from allophycocyanin in a previous directed evolution. Here, we report the smURFP crystal structure to better understand properties and enable further engineering of improved variants. We compare this structure to the structures of allophycocyanin and smURFP mutants to identify the structural origins of the molecular brightness. We then use a structure-guided approach to develop monomeric smURFP variants that fluoresce with phycocyanobilin but not biliverdin. Furthermore, we measure smURFP photophysical properties necessary for advanced imaging modalities, such as those relevant for two-photon, fluorescence lifetime, and single-molecule imaging. We observe that smURFP has the largest two-photon cross-section measured for a fluorescent protein, and that it produces more photons than organic dyes. Altogether, this study expands our understanding of the smURFP, which will inform future engineering toward optimal FPs compatible with whole organism studies.

The effects of biliverdin on pressure-induced unfolding of apomyoglobin: The specific role of Zn2+ ions

Apomyoglobin (apoMb), a model protein in biochemistry, exhibits a strong propensity to bind various ligands, which makes it a good candidate as a carrier of bioactive hydrophobic drugs. The stability of its hydrophobic pocket determines its potential as a carrier of bioactive compounds. High pressure (HP) is a potent tool for studying protein stability, revealing the specific role of hydrophobic cavities in unfolding. We probed the effects of biliverdin (BV) binding and its complex with Zn²⁺ ions on the structure and HP stability of apoMb. CD spectroscopy and SAXS measurements revealed that BV and BV-Zn²⁺ complexes make the apoMb structure more compact with higher α-helical content. We performed in-situ HP measurements of apoMb intrinsic fluorescence to demonstrate the ability of BV to stabilise apoMb structure at HP conditions. Furthermore, the presence of Zn²⁺ within the apoMb-BV complex significantly enhances the BV stabilisation effect. In-situ visible absorption study of BV chromophore confirmed the ability of Zn²⁺ to increase the stability of apoMb-BV complex under HP: the onset of complex dissociation is shifted by ~100 MPa in the presence of Zn²⁺. By combining HP-fluorescence and HP-visible absorption spectroscopy, our strategy highlights the crucial role of tetrapyrrole-metal complexes in stabilising apoMb hydrophobic pocket.

Dynamic Interconversion and Induced-Fit Guest Binding with Two Macrocycle-Based Coordination Cages

We report the syntheses and host-guest chemistry of two interconvertible coordination cages, Pd₂L₂ and Pd₁L₁, from a dynamic macrocycle ligand (L) and a cis-blocking (tmen)Pd(NO₃)₂ (tmen = tetramethylethylenediamine) unit (Pd). The water-soluble macrocyclic L, which can bind various polycyclic aromatic hydrocarbon (PAH) guests in its cis-conformation, was constructed via four pyridinium bonds between two 2,4,6-tri(4-pyridyl)-1,3,5-triazine [TPT] panels and two p-xylene bridges. We selectively formed each cage either by changing the reaction concentration/solvent/temperature or through induced-fit guest encapsulation, while direct assembly of L and Pd resulted in a mixture of Pd₂L₂ and Pd₁L₁ in equilibrium. X-ray structures of the free ligand and the host-guest complexes confirmed the induce-fit adaptive changes in the ligand's conformation and the cage's cavity. This work demonstrates a useful strategy for designing multistimuli-responsive supramolecular hosts by coordination self-assembly with macrocyclic ligands featuring rich conformational freedom.

Cavity Characterization in Supramolecular Cages

Confining molecular guests within artificial hosts has provided a major driving force in the rational design of supramolecular cages with tailored properties. Over the last 30 years, a set of design strategies have been developed that enabled the controlled synthesis of a myriad of cages. Recently, there has been a growing interest in involving in silico methods in this toolbox. Cavity shape and size are important parameters that can be easily accessed by inexpensive geometric algorithms. Although these algorithms are well developed for the detection of nonartificial cavities (e.g., enzymes), they are not routinely used for the rational design of supramolecular cages. In order to test the capabilities of this tool, we have evaluated the performance and characteristics of seven different cavity characterization software in the context of 22 analogues of well-known supramolecular cages. Among the tested software, KVFinder project and Fpocket proved to be the most software to characterize supramolecular cavities. With the results of this work, we aim to popularize this underused technique within the supramolecular community.

PDA-Pred: Predicting the binding affinity of protein-DNA complexes using machine learning techniques and structural features

Protein-DNA interactions play an important role in various biological processes such as gene expression, replication, and transcription. Understanding the important features that dictate the binding affinity of protein-DNA complexes and predicting their affinities is important for elucidating their recognition mechanisms. In this work, we have collected the experimental binding free energy (ΔG) for a set of 391 Protein-DNA complexes and derived several structure-based features such as interaction energy, contact potentials, volume and surface area of binding site residues, base step parameters of the DNA and contacts between different types of atoms. Our analysis on relationship between binding affinity and structural features revealed that the important factors mainly depend on the number of DNA strands as well as functional and structural classes of proteins. Specifically, binding site properties such as number of atom contacts between the DNA and protein, volume of protein binding sites and interaction-based features such as interaction energies and contact potentials are important to understand the binding affinity. Further, we developed multiple regression equations for predicting the binding affinity of protein-DNA complexes belonging to different structural and functional classes. Our method showed an average correlation and mean absolute error of 0.78 and 0.98 kcal/mol, respectively, between the experimental and predicted binding affinities on a jack-knife test. We have developed a webserver, PDA-PreD (Protein-DNA Binding affinity predictor), for predicting the affinity of protein-DNA complexes and it is freely available at https://web.iitm.ac.in/bioinfo2/pdapred/.

RPflex: A Coarse-Grained Network Model for RNA Pocket Flexibility Study

RNA regulates various biological processes, such as gene regulation, RNA splicing, and intracellular signal transduction. RNA’s conformational dynamics play crucial roles in performing its diverse functions. Thus, it is essential to explore the flexibility characteristics of RNA, especially pocket flexibility. Here, we propose a computational approach, RPflex, to analyze pocket flexibility using the coarse-grained network model. We first clustered 3154 pockets into 297 groups by similarity calculation based on the coarse-grained lattice model. Then, we introduced the flexibility score to quantify the flexibility by global pocket features. The results show strong correlations between the flexibility scores and root-mean-square fluctuation (RMSF) values, with Pearson correlation coefficients of 0.60, 0.76, and 0.53 in Testing Sets I–III. Considering both flexibility score and network calculations, the Pearson correlation coefficient was increased to 0.71 in flexible pockets on Testing Set IV. The network calculations reveal that the long-range interaction changes contributed most to flexibility. In addition, the hydrogen bonds in the base–base interactions greatly stabilize the RNA structure, while backbone interactions determine RNA folding. The computational analysis of pocket flexibility could facilitate RNA engineering for biological or medical applications.

Human VH-based chimeric antigen receptor T cells targeting glypican 3 eliminate tumors in preclinical models of HCC

BACKGROUND AND AIMS

Efficacy of chimeric antigen receptor (CAR) T cells for treating solid tumors, including HCC, remains a challenge. Nanobodies are emerging building blocks of CAR T cells due to their small size and high expression. Membrane proximal sites have been shown as attractive epitopes of CAR T cells. However, current CAR formats are not tailored toward nanobodies or targeting membrane distal epitopes.

APPROACH AND RESULTS

Using hYP7 Fv (membrane proximal) and HN3 VH nanobody (membrane distal) as GPC3 targeting elements, we sought to determine how hinges and transmembrane portions of varying structures and sizes affect CAR T-cell function. We generated multiple permutations of CAR T cells containing CD8, CD28, IgG4, and Fc domains. We show that engineered HN3 CAR T cells can be improved by 2 independent, synergistic changes in the hinge and transmembrane domains. The T cells expressing the HN3 CAR which contains the hinge region of IgG4 and the CD28 transmembrane domain (HN3-IgG4H-CD28TM) exhibited high cytotoxic activity and caused complete HCC tumor eradication in immunodeficient mice. HN3-IgG4H-CD28TM CAR T cells were enriched for cytotoxic-memory CD8+ T cells and NFAT signals, and reduced β catenin levels in HCC cells.

CONCLUSION

Our findings indicate that altering the hinge and transmembrane domains of a nanobody-based CAR targeting a distal GPC3 epitope, in contrast to a membrane proximal epitope, lead to robust T-cell signaling and induce swift and durable eradication of HCC tumors.

Laboratory Scale Production of Complex Proteins Using Charge Complimentary Nanoenvironments

Protein refolding is a crucial procedure in bacterial recombinant expression. Aggregation and misfolding are the two challenges that can affect the overall yield and specific activity of the folded proteins. We demonstrated the in vitro use of nanoscale "thermostable exoshells" (tES) to encapsulate, fold and release diverse protein substrates. With tES, the soluble yield, functional yield, and specific activity increased from 2-fold to >100-fold when compared to folding in its absence. On average, the soluble yield was determined to be 6.5 mg/100 mg of tES for a set of 12 diverse substrates evaluated. The electrostatic charge complementation between the tES interior and the protein substrate was considered as the primary determinant for functional folding. We thus describe a useful and simple method for in vitro folding that has been evaluated and implemented in our laboratory.

Heat-induced transitions of an empty minute virus of mice capsid in explicit water: all-atom MD simulation

The capsid-like structure of the virus-based protein nanoparticles (NPs) can serve as bionanomaterials, with applications in biomedicines and nanotechnology. Release of packaged material from these nanocontainers is associated with subtle conformational changes of the NP structure, which in vitro, is readily accomplished by heating. Characterizing the structural changes as a function of temperature may provide fresh insights into nanomaterial/antiviral strategies. Here, we have calculated heat induced changes in the properties of an empty minute virus of mice particle using large-scale ≈ 3.0 × 10⁶ all-atom molecular dynamics simulations. We focus on two heat induced structural changes of the NP, namely, dynamical transition (DT) and breathing transition (BT), both characterized by sudden and sharp change of measured parameters at temperatures, T_DT and T_BT, respectively. While DT is assessed by mean-square fluctuation of hydrogen atoms of the NP, BT is monitored through internal volume and permeation rate of water molecules through the NP. Both the transitions, resulting primarily from collective atomistic motion, are found to occur at temperatures widely separated from one another (T_BT>T_DT). The breathing motions, responsible for the translocation events of the packaged materials through the NP to kick off, are further probed by computing atomic resolution stresses from NVE simulations. Distribution of equilibrium atomistic stresses on the NP reveals a largely asymmetric nature and suggests structural breathing may actually represent large dynamic changes in the hotspot regions, far from the NP pores, which is in remarkable resemblance with recently conducted hydrogen-deuterium exchange coupled to mass spectrometry experiment. Communicated by Ramaswamy H. Sarma.

Bacterial crystalline cellulose secretion via a supramolecular BcsHD scaffold

Cellulose, the most abundant biopolymer on Earth, is not only the predominant constituent of plants but also a key extracellular polysaccharide in the biofilms of many bacterial species. Depending on the producers, chemical modifications, and three-dimensional assemblies, bacterial cellulose (BC) can present diverse degrees of crystallinity. Highly ordered, or crystalline, cellulose presents great economical relevance due to its ever-growing number of biotechnological applications. Even if some acetic acid bacteria have long been identified as BC superproducers, the molecular mechanisms determining the secretion of crystalline versus amorphous cellulose remain largely unknown. Here, we present structural and mechanistic insights into the role of the accessory subunits BcsH (CcpAx) and BcsD (CesD) that determine crystalline BC secretion in the Gluconacetobacter lineage. We show that oligomeric BcsH drives the assembly of BcsD into a supramolecular cytoskeletal scaffold that likely stabilizes the cellulose-extruding synthase nanoarrays through an unexpected inside-out mechanism for secretion system assembly.

Biophysical characterization of calcium-binding and modulatory-domain dynamics in a pentameric ligand-gated ion channel

Pentameric ligand-gated ion channels (pLGICs) perform electrochemical signal transduction in organisms ranging from bacteria to humans. Among the prokaryotic pLGICs, there is architectural diversity involving N-terminal domains (NTDs) not found in eukaryotic relatives, exemplified by the calcium-sensitive channel (DeCLIC) from a Desulfofustis deltaproteobacterium, which has an NTD in addition to the canonical pLGIC structure. Here, we have characterized the structure and dynamics of DeCLIC through cryoelectron microscopy (cryo-EM), small-angle neutron scattering (SANS), and molecular dynamics (MD) simulations. In the presence and absence of calcium, cryo-EM yielded structures with alternative conformations of the calcium-binding site. SANS profiles further revealed conformational diversity at room temperature beyond that observed in static structures, shown through MD to be largely attributable to rigid-body motions of the NTD relative to the protein core, with expanded and asymmetric conformations improving the fit of the SANS data. This work reveals the range of motion available to the DeCLIC NTD and calcium-binding site, expanding the conformational landscape of the pLGIC family. Further, these findings demonstrate the power of combining low-resolution scattering, high-resolution structural, and MD simulation data to elucidate interfacial interactions that are highly conserved in the pLGIC family.

Docking-based long timescale simulation of cell-size protein systems at atomic resolution

Advances in computational modeling have led to an increasing focus on larger biomolecular systems, up to the level of a cell. Protein interactions are a central component of cellular processes. Techniques for modeling protein interactions have been divided between two fields: protein docking (predicting the static structures of protein complexes) and molecular simulation (modeling the dynamics of protein association, for relatively short simulation times at atomic resolution). Our study combined the two approaches to reach very long simulation times. The study makes the model more adequate to the real cells, to explore cellular processes at atomic resolution to better understand molecular mechanisms of life, and to use this knowledge to improve our ability to treat diseases.

Computational methodologies are increasingly addressing modeling of the whole cell at the molecular level. Proteins and their interactions are the key component of cellular processes. Techniques for modeling protein interactions, thus far, have included protein docking and molecular simulation. The latter approaches account for the dynamics of the interactions but are relatively slow, if carried out at all-atom resolution, or are significantly coarse grained. Protein docking algorithms are far more efficient in sampling spatial coordinates. However, they do not account for the kinetics of the association (i.e., they do not involve the time coordinate). Our proof-of-concept study bridges the two modeling approaches, developing an approach that can reach unprecedented simulation timescales at all-atom resolution. The global intermolecular energy landscape of a large system of proteins was mapped by the pairwise fast Fourier transform docking and sampled in space and time by Monte Carlo simulations. The simulation protocol was parametrized on existing data and validated on a number of observations from experiments and molecular dynamics simulations. The simulation protocol performed consistently across very different systems of proteins at different protein concentrations. It recapitulated data on the previously observed protein diffusion rates and aggregation. The speed of calculation allows reaching second-long trajectories of protein systems that approach the size of the cells, at atomic resolution.

Crystal structure and biochemical analysis suggest that YjoB ATPase is a putative substrate-specific molecular chaperone

AAA+ ATPases are ubiquitous proteins associated with most cellular processes, including DNA unwinding and protein unfolding. Their functional and structural properties are typically determined by domains and motifs added to the conserved ATPases domain. Currently, the molecular function and structure of many ATPases remain elusive. Here, we report the crystal structure and biochemical analyses of YjoB, a Bacillus subtilis AAA+ protein. The crystal structure revealed that the YjoB hexamer forms a bucket hat-shaped structure with a porous chamber. Biochemical analyses showed that YjoB prevents the aggregation of vegetative catalase KatA and gluconeogenesis-specific glyceraldehyde-3 phosphate dehydrogenase GapB but not citrate synthase, a conventional substrate. Structural and biochemical analyses further showed that the internal chamber of YjoB is necessary for inhibition of substrate aggregation. Our results suggest that YjoB, conserved in the class Bacilli, is a potential molecular chaperone acting in the starvation/stationary phases of B. subtilis growth.

Structure of Mycobacterium tuberculosis Cya, an evolutionary ancestor of the mammalian membrane adenylyl cyclases

Mycobacterium tuberculosis adenylyl cyclase (AC) Rv1625c / Cya is an evolutionary ancestor of the mammalian membrane ACs and a model system for studies of their structure and function. Although the vital role of ACs in cellular signaling is well established, the function of their transmembrane (TM) regions remains unknown. Here we describe the cryo-EM structure of Cya bound to a stabilizing nanobody at 3.6 Å resolution. The TM helices 1-5 form a structurally conserved domain that facilitates the assembly of the helical and catalytic domains. The TM region contains discrete pockets accessible from the extracellular and cytosolic side of the membrane. Neutralization of the negatively charged extracellular pocket Ex1 destabilizes the cytosolic helical domain and reduces the catalytic activity of the enzyme. The TM domain acts as a functional component of Cya, guiding the assembly of the catalytic domain and providing the means for direct regulation of catalytic activity in response to extracellular ligands.

Interaction of the new inhibitor paxlovid (PF-07321332) and ivermectin with the monomer of the main protease SARS-CoV-2: A volumetric study based on molecular dynamics, elastic networks, classical thermodynamics and SPT

The COVID-19 pandemic has accelerated the study of drugs, most notably ivermectin and more recently Paxlovid (PF-07321332) which is in phase III clinical trials with experimental data showing covalent binding to the viral protease M^pro. Theoretical developments of catalytic site-directed docking support thermodynamically feasible non-covalent binding to M^pro. Here we show that Paxlovid binds non-covalently at regions other than the catalytic sites with energies stronger than reported and at the same binding site as the ivermectin B1a homologue, all through theoretical methodologies, including blind docking. We volumetrically characterize the non-covalent interaction of the ivermectin homologues (avermectins B1a and B1b) and Paxlovid with the mM^pro monomer, through molecular dynamics and scaled particle theory (SPT). Using the fluctuation-dissipation theorem (FDT), we estimated the electric dipole moment fluctuations at the surface of each of complex involved in this study, with similar trends to that observed in the interaction volume. Using fluctuations of the intrinsic volume and the number of flexible fragments of proteins using anisotropic and Gaussian elastic networks (ANM+GNM) suggests the complexes with ivermectin are more dynamic and flexible than the unbound monomer. In contrast, the binding of Paxlovid to mM^pro shows that the mM^pro-PF complex is the least structurally dynamic of all the species measured in this investigation. The results support a differential molecular mechanism of the ivermectin and PF homologues in the mM^pro monomer. Finally, the results showed that Paxlovid despite beingbound in different sites through covalent or non-covalent forms behaves similarly in terms of its structural flexibility and volumetric behaviour.