Prediction of hydrodynamic and other solution properties of rigid proteins from atomic- and residue-level models

Diffusion of proteins in crowded solutions studied by docking-based modeling

The diffusion of proteins is significantly affected by macromolecular crowding. Molecular simulations accounting for protein interactions at atomic resolution are useful for characterizing the diffusion patterns in crowded environments. We present a comprehensive analysis of protein diffusion under different crowding conditions based on our recent docking-based approach simulating an intracellular crowded environment by sampling the intermolecular energy landscape using the Markov Chain Monte Carlo protocol. The procedure was extensively benchmarked, and the results are in very good agreement with the available experimental and theoretical data. The translational and rotational diffusion rates were determined for different types of proteins under crowding conditions in a broad range of concentrations. A protein system representing most abundant protein types in the E. coli cytoplasm was simulated, as well as large systems of other proteins of varying sizes in heterogeneous and self-crowding solutions. Dynamics of individual proteins was analyzed as a function of concentration and different diffusion rates in homogeneous and heterogeneous crowding. Smaller proteins diffused faster in heterogeneous crowding of larger molecules, compared to their diffusion in the self-crowded solution. Larger proteins displayed the opposite behavior, diffusing faster in the self-crowded solution. The results show the predictive power of our structure-based simulation approach for long timescales of cell-size systems at atomic resolution.

The ribosome lowers the entropic penalty of protein folding

Most proteins fold during biosynthesis on the ribosome1, and co-translational folding energetics, pathways and outcomes of many proteins have been found to differ considerably from those in refolding studies2-10. The origin of this folding modulation by the ribosome has remained unknown. Here we have determined atomistic structures of the unfolded state of a model protein on and off the ribosome, which reveal that the ribosome structurally expands the unfolded nascent chain and increases its solvation, resulting in its entropic destabilization relative to the peptide chain in isolation. Quantitative 19F NMR experiments confirm that this destabilization reduces the entropic penalty of folding by up to 30 kcal mol-1 and promotes formation of partially folded intermediates on the ribosome, an observation that extends to other protein domains and is obligate for some proteins to acquire their active conformation. The thermodynamic effects also contribute to the ribosome protecting the nascent chain from mutation-induced unfolding, which suggests a crucial role of the ribosome in supporting protein evolution. By correlating nascent chain structure and dynamics to their folding energetics and post-translational outcomes, our findings establish the physical basis of the distinct thermodynamics of co-translational protein folding.

Structure-driven development of a biomimetic rare earth artificial metalloprotein

The 2011 discovery of the first rare earth-dependent enzyme in methylotrophic Methylobacterium extorquens AM1 prompted intensive research toward understanding the unique chemistry at play in these systems. This enzyme, an alcohol dehydrogenase (ADH), features a La3+ ion closely associated with redox-active coenzyme pyrroloquinoline quinone (PQQ) and is structurally homologous to the Ca2+-dependent ADH from the same organism. AM1 also produces a periplasmic PQQ-binding protein, PqqT, which we have now structurally characterized to 1.46-Å resolution by X-ray diffraction. This crystal structure reveals a Lys residue hydrogen-bonded to PQQ at the site analogously occupied by a Lewis acidic cation in ADH. Accordingly, we prepared K142A- and K142D-PqqT variants to assess the relevance of this site toward metal binding. Isothermal titration calorimetry experiments and titrations monitored by UV-Vis absorption and emission spectroscopies support that K142D-PqqT binds tightly (Kd = 0.6 ± 0.2 μM) to La3+ in the presence of bound PQQ and produces spectral signatures consistent with those of ADH enzymes. These spectral signatures are not observed for WT- or K142A-variants or upon addition of Ca2+ to PQQ ⸦ K142D-PqqT. Addition of benzyl alcohol to La3+-bound PQQ ⸦ K142D-PqqT (but not Ca2+-bound PQQ ⸦ K142D-PqqT, or La3+-bound PQQ ⸦ WT-PqqT) produces spectroscopic changes associated with PQQ reduction, and chemical trapping experiments reveal the production of benzaldehyde, supporting ADH activity. By creating a metal binding site that mimics native ADH enzymes, we present a rare earth-dependent artificial metalloenzyme primed for future mechanistic, biocatalytic, and biosensing applications.

Rapid and reversible dissolution of biomolecular condensates using light-controlled recruitment of a solubility tag

Biomolecular condensates are broadly implicated in both normal cellular regulation and disease. Consequently, several chemical biology and optogenetic approaches have been developed to induce phase separation of a protein of interest. However, few tools are available to perform the converse function - dissolving a condensate of interest on demand. Such a tool would aid in testing whether the condensate plays specific functional roles. Here we show that light-gated recruitment of a solubilizing domain, maltose-binding protein (MBP), results in rapid and controlled dissolution of condensates formed from proteins of interest. Our optogenetic MBP-based dissolution strategy (OptoMBP) is rapid, reversible, and can be spatially controlled with subcellular precision. We also provide a proof-of-principle application of OptoMBP by disrupting condensation of the oncogenic fusion protein FUS-CHOP and reverting FUS-CHOP driven transcriptional changes. We envision that the OptoMBP system could be broadly useful for disrupting constitutive protein condensates to probe their biological functions.

Dynamics of side chains in poly(quinoxaline-2,3-diyl)s studied via quasielastic neutron scattering

The side chain dynamics of poly(quinoxaline-2,3-diyl)s (PQXs) are expected to influence their conformation. To investigate these dynamics experimentally, quasielastic neutron scattering (QENS) was performed for PQXs in deuterated tetrahydrofuran (THF-d8) and deuterated 1,1,2-trichloroethane/THF (1,1,2-TCE-d3/THF-d8), in which they formed right-handed and left-handed helical structures, respectively. The mean-square displacement of the PQX side chains in 1,1,2-TCE-d3/THF-d8 was lower than that in THF-d8. Furthermore, QENS complementary studies and molecular dynamics simulations unraveled a coupling between the main-chain and side chain dynamics of PQXs, suggesting the possibility of controlling the main-chain helical chirality through the dynamics of chiral side chains. These insights present a novel strategy for the design of synthetic helical macromolecules with precise chirality control.

Analysis of homodimer formation in 12-oxophytodienoate reductase 3 in solutio and crystallo challenges the physiological role of the dimer

12-oxophytodienoate reductase 3 (OPR3) is a key enzyme in the biosynthesis of jasmonoyl-L-isoleucine, the receptor-active form of jasmonic acid and crucial signaling molecule in plant defense. OPR3 was initially crystallized as a self-inhibitory dimer, implying that homodimerization regulates enzymatic activity in response to biotic and abiotic stresses. Since a sulfate ion is bound to Y364, mimicking a phosphorylated tyrosine, it was suggested that dimer formation might be controlled by reversible phosphorylation of Y364 in vivo. To investigate OPR3 homodimerization and its potential physiological role in more detail, we performed analytical gel filtration and dynamic light scattering on wild-type OPR3 and three variants (R283D, R283E, and Y364P). The experiments revealed a rapid and highly sensitive monomer-dimer equilibrium for all OPR3 constructs. We crystallized all constructs with and without sulfate to examine its effect on the dimerization process and whether reversible phosphorylation of Y364 triggers homodimerization in vivo. All OPR3 constructs crystallized in their monomeric and dimeric forms independent of the presence of sulfate. Even variant Y364P, lacking the putative phosphorylation site, was crystallized as a self-inhibitory homodimer, indicating that Y364 is not required for dimerization. Generally, the homodimer is relatively weak, and our results raise doubts about its physiological role in regulating jasmonate biosynthesis.

E41K mutation activates Bruton's tyrosine kinase by stabilizing an inositol hexakisphosphate-dependent invisible dimer

Bruton's tyrosine kinase (BTK) regulates diverse cellular signaling of the innate and adaptive immune system in response to microbial pathogens. Downregulation or constitutive activation of BTK is reported in patients with autoimmune diseases or various B-cell leukemias. BTK is a multidomain protein tyrosine kinase that adopts an Src-like autoinhibited conformation maintained by the interaction between the kinase and PH-TH domains. The PH-TH domain plays a central role in regulating BTK function. BTK is activated by binding to PIP3 at the plasma membrane upon stimulation by the B-cell receptor (BCR). The PIP3 binding allows dimerization of the PH-TH domain and subsequent transphosphorylation of the activation loop. Alternatively, a recent study shows that the multivalent T-cell-independent (TI) antigen induces BCR response by activating BTK independent of PIP3 binding. It was proposed that a transiently stable IP6-dependent PH-TH dimer may activate BTK during BCR activation by the TI antigens. However, no IP6-dependent PH-TH dimer has been identified yet. Here, we investigated a constitutively active PH-TH mutant (E41K) to determine if the elusive IP6-dependent PH-TH dimer exists. We showed that the constitutively active E41K mutation activates BTK by stabilizing the IP6-dependent PH-TH dimer. We observed that a downregulating mutation in the PH-TH domain (R28H) linked to X-linked agammaglobulinemia impairs BTK activation at the membrane and in the cytosol by preventing PH-TH dimerization. We conclude that the IP6 dynamically remodels the BTK active fraction between the membrane and the cytoplasm. Stimulating with IP6 increases the cytosolic fraction of the activated BTK.

Characterization and structural basis for the brightness of mCLIFY: a novel monomeric and circularly permuted bright yellow fluorescent protein

We present mCLIFY: a monomeric, bright, yellow, and long-lived fluorescent protein (FP) created by circular permutation of YPet, the brightest yellow FP from Aequorea Victoria for use in cellular and in vitro single molecule studies. mCLIFY retains the enhanced photophysical properties of YPET as a monomer at concentrations ≤ 40 μM. In contrast, we determined that YPet has a dimerization dissociation constant (K D 1-2) of 3.4 μM. Dimerization of YPet can cause homo-FRET, which underlies quantitative errors due to dimerization and homo-FRET. We determined the atomic structure of mCLIFY at 1.57 Å resolution and used its similarity with Venus for guided chromophore-targeted substitution studies to provide insights into its enhanced photophysical properties. The mutation V58L within the chromophore pocket improved quantum yield and extinction coefficient, making mCLIFY ~30% brighter than Venus. The extensive characterization of the photophysical and structural properties of YPet and mCLIFY presented here allowed us to reveal the basis of their long lifetimes and enhanced brightness and the basis of YPet's dimerization.

Strategy for modeling higher-order G-quadruplex structures recalcitrant to NMR determination

Guanine-rich nucleic acids can form intramolecularly folded four-stranded structures known as G-quadruplexes (G4s). Traditionally, G4 research has focused on short, highly modified DNA or RNA sequences that form well-defined homogeneous compact structures. However, the existence of longer sequences with multiple G4 repeats, from proto-oncogene promoters to telomeres, suggests the potential for more complex higher-order structures with multiple G4 units that might offer selective drug-targeting sites for therapeutic development. These larger structures present significant challenges for structural characterization by traditional high-resolution methods like multi-dimensional NMR and X-ray crystallography due to their molecular complexity. To address this current challenge, we have developed an integrated structural biology (ISB) platform, combining experimental and computational methods to determine self-consistent molecular models of higher-order G4s (xG4s). Here we outline our ISB method using two recent examples from our lab, an extended c-Myc promoter and long human telomere G4 repeats, that highlights the utility and generality of our approach to characterizing biologically relevant xG4s.

Location Is Everything: Influence of His-Tag Fusion Site on Properties of Adenylosuccinate Synthetase from Helicobacter pylori

The requirement for fast and dependable protein purification methods is constant, either for functional studies of natural proteins or for the production of biotechnological protein products. The original procedure has to be formulated for each individual protein, and this demanding task was significantly simplified by the introduction of affinity tags. Helicobacter pylori adenylosuccinate synthetase (AdSS) is present in solution in a dynamic equilibrium of monomers and biologically active homodimers. The addition of the His6-tag on the C-terminus (C-His-AdSS) was proven to have a negligible effect on the characteristics of this enzyme. This paper shows that the same enzyme with the His6-tag fused on its N-terminus (N-His-AdSS) has a high tendency to precipitate. Circular dichroism and X-ray diffraction studies do not detect any structural change that could explain this propensity. However, the dynamic light scattering, differential scanning fluorimetry, and analytical ultracentrifugation measurements indicate that the monomer of this construct is prone to aggregation, which shifts the equilibrium towards the insoluble precipitant. In agreement, enzyme kinetics measurements showed reduced enzyme activity, but preserved affinity for the substrates, in comparison with the wild-type and C-His-AdSS. The presented results reinforce the notion that testing the influence of the tag on protein properties should not be overlooked.

Quantitative analysis of the slow exchange process by 19F NMR in the presence of scalar and dipolar couplings: applications to the ribose 2'-19F probe in nucleic acids

Solution NMR spectroscopy is a particularly powerful technique for characterizing the functional dynamics of biomolecules, which is typically achieved through the quantitative characterization of chemical exchange processes via the measurement of spin relaxation rates. In addition to the conventional nuclei such as 15N and 13C, which are abundant in biomolecules, fluorine-19 (19F) has recently garnered attention and is being widely used as a site-specific spin probe. While 19F offers the advantages of high sensitivity and low background, it can be susceptible to artifacts in quantitative relaxation analyses due to a multitude of dipolar and scalar coupling interactions with nearby 1H spins. In this study, we focused on the ribose 2'-19F spin probe in nucleic acids and investigated the effects of 1H-19F spin interactions on the quantitative characterization of slow exchange processes on the millisecond time scale. We demonstrated that the 1H-19F dipolar coupling can significantly affect the interpretation of 19F chemical exchange saturation transfer (CEST) experiments when 1H decoupling is applied, while the 1H-19F interactions have a lesser impact on Carr-Purcell-Meiboom-Gill relaxation dispersion applications. We also proposed a modified CEST scheme to alleviate these artifacts along with experimental verifications on self-complementary RNA systems. The theoretical framework presented in this study can be widely applied to various 19F spin systems where 1H-19F interactions are operative, further expanding the utility of 19F relaxation-based NMR experiments.

Structural characterization of PHOX2B and its DNA interaction shed light on the molecular basis of the +7Ala variant pathogenicity in CCHS

An expansion of poly-alanine up to +13 residues in the C-terminus of the transcription factor PHOX2B underlies the onset of congenital central hypoventilation syndrome (CCHS). Recent studies demonstrated that the alanine tract expansion influences PHOX2B folding and activity. Therefore, structural information on PHOX2B is an important target for obtaining clues to elucidate the insurgence of the alanine expansion-related syndrome and also for defining a viable therapy. Here we report by NMR spectroscopy the structural characterization of the homeodomain (HD) of PHOX2B and HD + C-terminus PHOX2B protein, free and in the presence of the target DNA. The obtained structural data are then exploited to obtain a structural model of the PHOX2B-DNA interaction. In addition, the variant +7Ala, responsible for one of the most frequent forms of the syndrome, was analysed, showing different conformational proprieties in solution and a strong propensity to aggregation. Our data suggest that the elongated poly-alanine tract would be related to disease onset through a loss-of-function mechanism. Overall, this study paves the way for the future rational design of therapeutic drugs, suggesting as a possible therapeutic route the use of specific anti-aggregating molecules capable of preventing variant aggregation and possibly restoring the DNA-binding activity of PHOX2B.

Stability and deformation of biomolecular condensates under the action of shear flow

Biomolecular condensates play a key role in cytoplasmic compartmentalization and cell functioning. Despite extensive research on the physico-chemical, thermodynamic, or crowding aspects of the formation and stabilization of the condensates, one less studied feature is the role of external perturbative fluid flow. In fact, in living cells, shear stress may arise from streaming or active transport processes. Here, we investigate how biomolecular condensates are deformed under different types of shear flows. We first model Couette flow perturbations via two-way coupling between the condensate dynamics and fluid flow by deploying Lattice Boltzmann Molecular Dynamics. We then show that a simplified approach where the shear flow acts as a static perturbation (one-way coupling) reproduces the main features of the condensate deformation and dynamics as a function of the shear rate. With this approach, which can be easily implemented in molecular dynamics simulations, we analyze the behavior of biomolecular condensates described through residue-based coarse-grained models, including intrinsically disordered proteins and protein/RNA mixtures. At lower shear rates, the fluid triggers the deformation of the condensate (spherical to oblated object), while at higher shear rates, it becomes extremely deformed (oblated or elongated object). At very high shear rates, the condensates are fragmented. We also compare how condensates of different sizes and composition respond to shear perturbation, and how their internal structure is altered by external flow. Finally, we consider the Poiseuille flow that realistically models the behavior in microfluidic devices in order to suggest potential experimental designs for investigating fluid perturbations in vitro.

Functional diversity of YbbN/CnoX proteins: Insights from a comparative analysis of three thioredoxin-like oxidoreductases from Pseudomonas aeruginosa, Xylella fastidiosa and Escherichia coli

YbbN/CnoX are proteins that display a Thioredoxin (Trx) domain linked to a tetratricopeptide domain. YbbN from Escherichia coli (EcYbbN) displays a co-chaperone (holdase) activity that is induced by HOCl. Here, we compared EcYbbN with YbbN proteins from Xylella fastidiosa (XfYbbN) and from Pseudomonas aeruginosa (PaYbbN). EcYbbN presents a redox active Cys residue at Trx domain (Cys63), 24 residues away from SQHC motif (SQHC[N24]C) that can form mixed disulfides with target proteins. In contrast, XfYbbN and PaYbbN present two Cys residues in the CXXC (CAPC) motif, while only PaYbbN shows the Cys residue equivalent to Cys63 of EcYbbN. Our phylogenetic analysis revealed that most of the YbbN proteins are in the bacteria domain of life and that their members can be divided into four groups according to the conserved Cys residues. EcYbbN (SQHC[N24]C), XfYbbN (CAPC[N24]V) and PaYbbN (CAPC[N24]C) are representatives of three sub-families. In contrast to EcYbbN, both XfYbbN and PaYbbN: (1) reduced an artificial disulfide (DTNB) and (2) supported the peroxidase activity of Peroxiredoxin Q from X. fastidiosa, suggesting that these proteins might function similarly to the canonical Trx enzymes. Indeed, XfYbbN was reduced by XfTrx reductase with a high catalytic efficiency (kcat/Km = 1.27 x 107 M-1 s-1), similar to the canonical XfTrx (XfTsnC). Furthermore, EcYbbN and XfYbbN, but not PaYbbN displayed HOCl-induced holdase activity. Remarkably, EcYbbN gained disulfide reductase activity while lost the HOCl-activated chaperone function, when the SQHC was replaced by CQHC. In contrast, the XfYbbN CAPA mutant lost the disulfide reductase activity, while kept its HOCl-induced chaperone function. In all cases, the induction of the holdase activity was accompanied by YbbN oligomerization. Finally, we showed that deletion of ybbN gene did not render in P. aeruginosa more sensitive stressful treatments. Therefore, YbbN/CnoX proteins display distinct properties, depending on the presence of the three conserved Cys residues.

NMR characterization of RNA binding property of the DEAD-box RNA helicase DDX3X and its implications for helicase activity

The DEAD-box RNA helicase (DDX) plays a central role in many aspects of RNA metabolism by remodeling the defined structure of RNA molecules. While a number of structural studies have revealed the atomistic details of the interaction between DDX and RNA ligands, the molecular mechanism of how this molecule unwinds a structured RNA into an unstructured single-stranded RNA (ssRNA) has largely remained elusive. This is due to challenges in structurally characterizing the unwinding intermediate state and the lack of thermodynamic details underlying this process. In this study, we use solution nuclear magnetic resonance (NMR) spectroscopy to characterize the interaction of human DDX3X, a member of the DDX family, with various RNA ligands. Our results show that the inherent binding affinity of DDX3X for ssRNA is significantly higher than that for structured RNA elements. This preferential binding, accompanied by the formation of a domain-closed conformation in complex with ssRNA, effectively stabilizes the denatured ssRNA state and thus underlies the unwinding activity of DDX3X. Our results provide a thermodynamic and structural basis for the DDX function, whereby DDX can recognize and remodel a distinct set of structured RNAs to participate in a wide range of physiological processes.

Molecular and thermodynamic determinants of self-assembly and hetero-oligomerization in the enterobacterial thermo-osmo-regulatory protein H-NS

Environmentally regulated gene expression is critical for bacterial survival under stress conditions, including extremes in temperature, osmolarity and nutrient availability. Here, we dissect the thermo- and osmo-responsory behavior of the transcriptional repressor H-NS, an archetypal nucleoid-condensing sensory protein, ubiquitous in enterobacteria that infect the mammalian gut. Through experiments and thermodynamic modeling, we show that H-NS exhibits osmolarity, temperature and concentration dependent self-association, with a highly polydisperse native ensemble dominated by monomers, dimers, tetramers and octamers. The relative population of these oligomeric states is determined by an interplay between dimerization and higher-order oligomerization, which in turn drives a competition between weak homo- versus hetero-oligomerization of protein-protein and protein-DNA complexes. A phosphomimetic mutation, Y61E, fully eliminates higher-order self-assembly and preserves only dimerization while weakening DNA binding, highlighting that oligomerization is a prerequisite for strong DNA binding. We further demonstrate the presence of long-distance thermodynamic connectivity between dimerization and oligomerization sites on H-NS which influences the binding of the co-repressor Cnu, and switches the DNA binding mode of the hetero-oligomeric H-NS:Cnu complex. Our work thus uncovers important organizational principles in H-NS including a multi-layered thermodynamic control, and provides a molecular framework broadly applicable to other thermo-osmo sensory proteins that employ similar mechanisms to regulate gene expression.

Flexible structural arrangement and DNA-binding properties of protein p6 from Bacillus subtillis phage φ29

The genome-organizing protein p6 of Bacillus subtilis bacteriophage φ29 plays an essential role in viral development by activating the initiation of DNA replication and participating in the early-to-late transcriptional switch. These activities require the formation of a nucleoprotein complex in which the DNA adopts a right-handed superhelix wrapping around a multimeric p6 scaffold, restraining positive supercoiling and compacting the viral genome. Due to the absence of homologous structures, prior attempts to unveil p6's structural architecture failed. Here, we employed AlphaFold2 to engineer rational p6 constructs yielding crystals for three-dimensional structure determination. Our findings reveal a novel fold adopted by p6 that sheds light on its self-association mechanism and its interaction with DNA. By means of protein-DNA docking and molecular dynamic simulations, we have generated a comprehensive structural model for the nucleoprotein complex that consistently aligns with its established biochemical and thermodynamic parameters. Besides, through analytical ultracentrifugation, we have confirmed the hydrodynamic properties of the nucleocomplex, further validating in solution our proposed model. Importantly, the disclosed structure not only provides a highly accurate explanation for previously experimental data accumulated over decades, but also enhances our holistic understanding of the structural and functional attributes of protein p6 during φ29 infection.

Programming interchangeable and reversible heterooligomeric protein self-assembly using a bifunctional ligand

Protein design for self-assembly allows us to explore the emergence of protein-protein interfaces through various chemical interactions. Heterooligomers, unlike homooligomers, inherently offer a comprehensive range of structural and functional variations. Besides, the macromolecular repertoire and their applications would significantly expand if protein components could be easily interchangeable. This study demonstrates that a rationally designed bifunctional linker containing an enzyme inhibitor and maleimide can guide the formation of diverse protein heterooligomers in an easily applicable and exchangeable manner without extensive sequence optimizations. As proof of concept, we selected four structurally and functionally unrelated proteins, carbonic anhydrase, aldolase, acetyltransferase, and encapsulin, as building block proteins. The combinations of two proteins with the bifunctional linker yielded four two-component heterooligomers with discrete sizes, shapes, and enzyme activities. Besides, the overall size and formation kinetics of the heterooligomers alter upon adding metal chelators, acidic buffer components, and reducing agents, showing the reversibility and tunability in the protein self-assembly. Given that the functional groups of both the linker and protein components are readily interchangeable, our work broadens the scope of protein-assembled architectures and their potential applications as functional biomaterials.

Biophysical and structural characterization of a multifunctional viral genome packaging motor

The large dsDNA viruses replicate their DNA as concatemers consisting of multiple covalently linked genomes. Genome packaging is catalyzed by a terminase enzyme that excises individual genomes from concatemers and packages them into preassembled procapsids. These disparate tasks are catalyzed by terminase alternating between two distinct states-a stable nuclease that excises individual genomes and a dynamic motor that translocates DNA into the procapsid. It was proposed that bacteriophage λ terminase assembles as an anti-parallel dimer-of-dimers nuclease complex at the packaging initiation site. In contrast, all characterized packaging motors are composed of five terminase subunits bound to the procapsid in a parallel orientation. Here, we describe biophysical and structural characterization of the λ holoenzyme complex assembled in solution. Analytical ultracentrifugation, small angle X-ray scattering, and native mass spectrometry indicate that 5 subunits assemble a cone-shaped terminase complex. Classification of cryoEM images reveals starfish-like rings with skewed pentameric symmetry and one special subunit. We propose a model wherein nuclease domains of two subunits alternate between a dimeric head-to-head arrangement for genome maturation and a fully parallel arrangement during genome packaging. Given that genome packaging is strongly conserved in both prokaryotic and eukaryotic viruses, the results have broad biological implications.

Rapid and reversible dissolution of biomolecular condensates using light-controlled recruitment of a solubility tag

Biomolecular condensates are broadly implicated in both normal cellular regulation and disease. Consequently, several chemical biology and optogenetic approaches have been developed to induce phase separation of a protein of interest. However, few tools are available to perform the converse function-dissolving a condensate of interest on demand. Such a tool would aid in testing whether the condensate plays specific functional roles, a major question in cell biology and drug development. Here we report an optogenetic approach to selectively dissolve a condensate of interest in a reversible and spatially controlled manner. We show that light-gated recruitment of maltose-binding protein (MBP), a commonly used solubilizing domain in protein purification, results in rapid and controlled dissolution of condensates formed from proteins of interest. Our optogenetic MBP-based dissolution strategy (OptoMBP) is rapid, reversible, and can be spatially controlled with subcellular precision. We also provide a proof-of-principle application of OptoMBP, showing that disrupting condensation of the oncogenic fusion protein FUS-CHOP results in reversion of FUS-CHOP driven transcriptional changes. We envision that the OptoMBP system could be broadly useful for disrupting constitutive protein condensates to probe their biological functions.

Ton motor conformational switch and peptidoglycan role in bacterial nutrient uptake

Active nutrient uptake is fundamental for survival and pathogenicity of Gram-negative bacteria, which operate a multi-protein Ton system to transport essential nutrients like metals and vitamins. This system harnesses the proton motive force at the inner membrane to energize the import through the outer membrane, but the mechanism of energy transfer remains enigmatic. Here, we study the periplasmic domain of ExbD, a crucial component of the proton channel of the Ton system. We show that this domain is a dynamic dimer switching between two conformations representing the proton channel's open and closed states. By in vivo phenotypic assays we demonstrate that this conformational switch is essential for the nutrient uptake by bacteria. The open state of ExbD triggers a disorder to order transition of TonB, enabling TonB to supply energy to the nutrient transporter. We also reveal the anchoring role of the peptidoglycan layer in this mechanism. Herein, we propose a mechanistic model for the Ton system, emphasizing ExbD duality and the pivotal catalytic role of peptidoglycan. Sequence analysis suggests that this mechanism is conserved in other systems energizing gliding motility and membrane integrity. Our study fills important gaps in understanding bacterial motor mechanism and proposes novel antibacterial strategies.

Acetylation regulates the oligomerization state and activity of RNase J, the Helicobacter pylori major ribonuclease

In the gastric pathogen Helicobacter pylori, post-transcriptional regulation relies strongly on the activity of the essential ribonuclease RNase J. Here, we elucidated the crystal and cryo-EM structures of RNase J and determined that it assembles into dimers and tetramers in vitro. We found that RNase J extracted from H. pylori is acetylated on multiple lysine residues. Alanine substitution of several of these residues impacts on H. pylori morphology, and thus on RNase J function in vivo. Mutations of Lysine 649 modulates RNase J oligomerization in vitro, which in turn influences ribonuclease activity in vitro. Our structural analyses of RNase J reveal loops that gate access to the active site and rationalizes how acetylation state of K649 can influence activity. We propose acetylation as a regulatory level controlling the activity of RNase J and its potential cooperation with other enzymes of RNA metabolism in H. pylori.

The Onset of Molecule-Spanning Dynamics in Heat Shock Protein Hsp90

Protein dynamics have been investigated on a wide range of time scales. Nano- and picosecond dynamics have been assigned to local fluctuations, while slower dynamics have been attributed to larger conformational changes. However, it is largely unknown how fast (local) fluctuations can lead to slow global (allosteric) changes. Here, fast molecule-spanning dynamics on the 100 to 200 ns time scale in the heat shock protein 90 (Hsp90) are shown. Global real-space movements are assigned to dynamic modes on this time scale, which is possible by a combination of single-molecule fluorescence, quasi-elastic neutron scattering and all-atom molecular dynamics (MD) simulations. The time scale of these dynamic modes depends on the conformational state of the Hsp90 dimer. In addition, the dynamic modes are affected to various degrees by Sba1, a co-chaperone of Hsp90, depending on the location within Hsp90, which is in very good agreement with MD simulations. Altogether, this data is best described by fast molecule-spanning dynamics, which precede larger conformational changes in Hsp90 and might be the molecular basis for allostery. This integrative approach provides comprehensive insights into molecule-spanning dynamics on the nanosecond time scale for a multi-domain protein.

Targeting NOX2 with Bivalent Small-Molecule p47phox-p22phox Inhibitors

Nicotinamide adenine dinucleotide phosphate oxidase isoform 2 (NOX2) is an enzymatic complex whose function is the regulated generation of reactive oxygen species (ROS). NOX2 activity is central to redox signaling events and antibacterial response, but excessive ROS production by NOX2 leads to oxidative stress and inflammation in a range of diseases. The protein-protein interaction between the NOX2 subunits p47phox and p22phox is essential for NOX2 activation, thus p47phox is a potential drug target. Previously, we identified 2-aminoquinoline as a fragment hit toward p47phoxSH3A-B and converted it to a bivalent small-molecule p47phox-p22phox inhibitor (Ki = 20 μM). Here, we systematically optimized the bivalent compounds by exploring linker types and positioning as well as substituents on the 2-aminoquinoline part and characterized the bivalent binding mode with biophysical methods. We identified several compounds with submicromolar binding affinities and cellular activity and thereby demonstrated that p47phox can be targeted by potent small molecules.

Analytical Ultracentrifugation Detects Quaternary Rearrangements and Antibody-Induced Conformational Selection of the SARS-CoV-2 Spike Trimer

Analytical ultracentrifugation (AUC) analysis shows that the SARS-CoV-2 trimeric Spike (S) protein adopts different quaternary conformations in solution. The relative abundance of the "open" and "close" conformations is temperature-dependent, and samples with different storage temperature history have different open/close distributions. Neutralizing antibodies (NAbs) targeting the S receptor binding domain (RBD) do not alter the conformer populations; by contrast, a NAb targeting a cryptic conformational epitope skews the Spike trimer toward an open conformation. The results highlight AUC, which is typically applied for molecular mass determination of biomolecules as a powerful tool for detecting functionally relevant quaternary protein conformations.

The solution structure of the unbound IgG Fc receptor CD64 resembles its crystal structure: Implications for function

FcγRI (CD64) is the only high-affinity Fcγ receptor found on monocytes, macrophages, eosinophils, neutrophils and dendritic cells. It binds immunoglobulin G (IgG) antibody-antigen complexes at its Fc region to trigger key immune responses. CD64 contains three immunoglobulin-fold extracellular domains (D1, D2 and D3) and a membrane-spanning region. Despite the importance of CD64, no solution structure for this is known to date. To investigate this, we used analytical ultracentrifugation, small-angle X-ray scattering, and atomistic modelling. Analytical ultracentrifugation revealed that CD64 was monomeric with a sedimentation coefficient s020,w of 2.53 S, together with some dimer. Small-angle X-ray scattering showed that its radius of gyration RG was 3.3-3.4 nm and increased at higher concentrations to indicate low dimerization. Monte Carlo modelling implemented in the SASSIE-web package generated 279,162 physically-realistic trial CD64 structures. From these, the scattering best-fit models at the lowest measured concentrations that minimised dimers revealed that the D1, D2 and D3 domains were structurally similar to those seen in three CD64 crystal structures, but showed previously unreported flexibility between D1, D2 and D3. Despite the limitations of the scattering data, the superimposition of the CD64 solution structures onto crystal structures of the IgG Fc-CD64 complex showed that the CD64 domains do not sterically clash with the IgG Fc region, i.e. the solution structure of CD64 was sufficiently compact to allow IgG to bind to its high-affinity Fcγ receptor. This improved understanding may result in novel approaches to inhibit CD64 function, and opens the way for the solution study of the full-length CD64-IgG complex.

Biophysical Determinants for the Viscosity of Concentrated Monoclonal Antibody Solutions

Monoclonal antibodies (mAbs) are particularly relevant for therapeutics due to their high specificity and versatility, and mAb-based drugs are hence used to treat numerous diseases. The increased patient compliance of self-administration motivates the formulation of products for subcutaneous (SC) administration. The associated challenge is to formulate highly concentrated antibody solutions to achieve a significant therapeutic effect, while limiting their viscosity and preserving their physicochemical stability. Protein-protein interactions (PPIs) are in fact the root cause of several potential problems concerning the stability, manufacturability, and delivery of a drug product. The understanding of macroscopic viscosity requires an in-depth knowledge on protein diffusion, PPIs, and self-association/aggregation. Here, we study the self-diffusion of different mAbs of the IgG1 subtype in aqueous solution as a function of the concentration and temperature by quasi-elastic neutron scattering (QENS). QENS allows us to probe the short-time self-diffusion of the molecules and therefore to determine the hydrodynamic mAb cluster size and to gain information on the internal mAb dynamics. Small-angle neutron scattering (SANS) is jointly employed to probe structural details and to understand the nature and intensity of PPIs. Complementary information is provided by molecular dynamics (MD) simulations and viscometry, thus obtaining a comprehensive picture of mAb diffusion.

Ton Motor Conformational Switch and Peptidoglycan Role in Bacterial Nutrient Uptake

Active nutrient uptake is fundamental for survival and pathogenicity of Gram-negative bacteria, which operate a multi-protein Ton system to transport essential nutrients like metals and vitamins. This system harnesses the proton motive force at the inner membrane to energize the import through the outer membrane, but the mechanism of energy transfer remains enigmatic. Here, we study the periplasmic domain of ExbD, a crucial component of the proton channel of the Ton system. We show that this domain is a dynamic dimer switching between two conformations representing the proton channel's open and closed states. By in vivo phenotypic assays we demonstrate that this conformational switch is essential for the nutrient uptake by bacteria. The open state of ExbD triggers a disorder to order transition of TonB, enabling TonB to supply energy to the nutrient transporter. We also reveal the anchoring role of the peptidoglycan layer in this mechanism. Herein, we propose a mechanistic model for the Ton system, emphasizing ExbD duality and the pivotal catalytic role of peptidoglycan. Sequence analysis suggests that this mechanism is conserved in other systems energizing gliding motility and membrane integrity. Our study fills important gaps in understanding bacterial motor mechanism and proposes novel antibacterial strategies.

Fast prediction of antibiotic permeability through membrane channels using Brownian dynamics

The efficient permeation across the Gram-negative bacterial membrane is an important step in the overall process of antibacterial action of a molecule and the one that has posed a significant hurdle on the way toward approved antibiotics. Predicting the permeability for a large library of molecules and assessing the effect of different molecular transformations on permeation rates of a given molecule is critical to the development of effective antibiotics. We present a computational approach for obtaining estimates of molecular permeability through a porin channel in a matter of hours using a Brownian dynamics approach. The fast sampling using a temperature acceleration scheme enables the approximate estimation of permeability using the inhomogeneous solubility diffusion model. Although the method is a significant approximation to similar all-atom approaches tested previously, we show that the present approach predicts permeabilities that correlate fairly well with the respective experimental permeation rates from liposome swelling experiments and accumulation rates from antibiotic accumulation assays, and is significantly, i.e., about 14 times, faster compared with a previously reported approach. The possible applications of the scheme in high-throughput screening for fast permeators are discussed.

Anomalous Oligomerization Behavior of E. coli Aquaporin Z in Detergent and in Nanodiscs

Aquaporins are tetrameric integral membrane proteins that act as water channels, and can also permeabilize membranes to other solutes. The monomer appears to be the functional form despite all aquaporins being organized as tetramers, which therefore must provide a clear functional advantage. In addition to this quaternary organization, some aquaporins can act as adhesion molecules in membrane junctions, when tetramers located in opposing membranes interact via their extracellular domains. These stacked forms have been observed in a range of aquaporins, whether using lipidic membrane environments, in electron crystallography, or using detergent micelles, in single-particle cryo-electron microscopy (cryo-EM). In the latter technique, structural studies can be performed when the aquaporin is reconstituted into nanodiscs of lipids that are surrounded by a protein scaffold. During attempts to study E. coli Aquaporin Z (AqpZ), we have found that in some conditions these nanodiscs tend to form filaments that appear to be either thicker head-to-tail or thinner side-to-side stacks of nanodiscs. Nanodisc oligomerization was observed using orthogonal analytical techniques analytical ultra-centrifugation and mass photometry, although the nature of the oligomers (head-to-tail or side-to-side) could not be determined. Using the latter technique, the AqpZ tetramer itself formed oligomers of increasing size when solubilized only in detergent, which is consistent with multiple stacking of AqpZ tetramers. We observed images consistent with both of these filaments in negative staining EM conditions, but only thicker filaments in cryo-EM conditions. We hypothesize that the apparent nanodisc side-to-side arrangement that can only be visualized in negative staining conditions is related to artifacts due to the sample preparation. Filaments of any kind were not observed in EM when nanodiscs did not contain AqpZ, or after addition of detergent into the nanodisc cryo-EM preparation, at concentrations that did not disrupt nanodisc formation. To our knowledge, these filaments have not been observed in nanodiscs preparations of other membrane proteins. AqpZ, like other aquaporins has a charge asymmetry between the cytoplasmic (more positive) and the extracellular sides, which may explain the likely head-to-tail stacking observed, both in nanodisc preparations and also in detergent micelles.

Optimized OPEP Force Field for Simulation of Crowded Protein Solutions

Macromolecular crowding has profound effects on the mobility of proteins, with strong implications on the rates of intracellular processes. To describe the dynamics of crowded environments, detailed molecular models are needed, capturing the structures and interactions arising in the crowded system. In this work, we present OPEPv7, which is a coarse-grained force field at amino-acid resolution, suited for rigid-body simulations of the structure and dynamics of crowded solutions formed by globular proteins. Using the OPEP protein model as a starting point, we have refined the intermolecular interactions to match the experimentally observed dynamical slowdown caused by crowding. The resulting force field successfully reproduces the diffusion slowdown in homogeneous and heterogeneous protein solutions at different crowding conditions. Coupled with the lattice Boltzmann technique, it allows the study of dynamical phenomena in protein assemblies and opens the way for the in silico rheology of protein solutions.

Varying molecular interactions explain aspects of crowder-dependent enzyme function of a viral protease

Biochemical processes in cells, including enzyme-catalyzed reactions, occur in crowded conditions with various background macromolecules occupying up to 40% of cytoplasm's volume. Viral enzymes in the host cell also encounter such crowded conditions as they often function at the endoplasmic reticulum membranes. We focus on an enzyme encoded by the hepatitis C virus, the NS3/4A protease, which is crucial for viral replication. We have previously found experimentally that synthetic crowders, polyethylene glycol (PEG) and branched polysucrose (Ficoll), differently affect the kinetic parameters of peptide hydrolysis catalyzed by NS3/4A. To gain understanding of the reasons for such behavior, we perform atomistic molecular dynamics simulations of NS3/4A in the presence of either PEG or Ficoll crowders and with and without the peptide substrates. We find that both crowder types make nanosecond long contacts with the protease and slow down its diffusion. However, they also affect the enzyme structural dynamics; crowders induce functionally relevant helical structures in the disordered parts of the protease cofactor, NS4A, with the PEG effect being more pronounced. Overall, PEG interactions with NS3/4A are slightly stronger but Ficoll forms more hydrogen bonds with NS3. The crowders also interact with substrates; we find that the substrate diffusion is reduced much more in the presence of PEG than Ficoll. However, contrary to NS3, the substrate interacts more strongly with Ficoll than with PEG crowders, with the substrate diffusion being similar to crowder diffusion. Importantly, crowders also affect the substrate-enzyme interactions. We observe that both PEG and Ficoll enhance the presence of substrates near the active site, especially near catalytic H57 but Ficoll crowders increase substrate binding more than PEG molecules.

Self-assembly of dengue virus empty capsid-like particles in solution

Nucleocapsid (NC) assembly is an essential step of the virus replication cycle. It ensures genome protection and transmission among hosts. Flaviviruses are human viruses for which envelope structure is well known, whereas no information on NC organization is available. Here we designed a dengue virus capsid protein (DENVC) mutant in which a highly positive spot conferred by arginine 85 in α4-helix was replaced by a cysteine residue, simultaneously removing the positive charge and restricting the intermolecular motion through the formation of a disulfide cross-link. We showed that the mutant self-assembles into capsid-like particles (CLP) in solution without nucleic acids. Using biophysical techniques, we investigated capsid assembly thermodynamics, showing that an efficient assembly is related to an increased DENVC stability due to α4/α4' motion restriction. To our knowledge, this is the first time that flaviviruses' empty capsid assembly is obtained in solution, revealing the R85C mutant as a powerful tool to understand the NC assembly mechanism.

Therapeutic phosphorodiamidate morpholino oligonucleotides: Physical properties, solution structures, and folding thermodynamics

Elucidating the structure-function relationships for therapeutic RNA mimicking phosphorodiamidate morpholino oligonucleotides (PMOs) is challenging due to the lack of information about their structures. While PMOs have been approved by the US Food and Drug Administration for treatment of Duchenne muscular dystrophy, no structural information on these unique, charge-neutral, and stable molecules is available. We performed circular dichroism and solution viscosity measurements combined with molecular dynamics simulations and machine learning to resolve solution structures of 22-mer, 25-mer, and 30-mer length PMOs. The PMO conformational dynamics are defined by the competition between non-polar nucleobases and uncharged phosphorodiamidate groups for shielding from solvent exposure. PMO molecules form non-canonical, partially helical, stable folded structures with a small 1.4- to 1.7-nm radius of gyration, low count of three to six base pairs and six to nine base stacks, characterized by -34 to -51 kcal/mol free energy, -57 to -103 kcal/mol enthalpy, and -23 to -53 kcal/mol entropy for folding. The 4.5- to 6.2-cm³/g intrinsic viscosity and Huggins constant of 4.5-9.9 are indicative of extended and aggregating systems. The results obtained highlight the importance of the conformational ensemble view of PMO solution structures, thermodynamic stability of their non-canonical structures, and concentration-dependent viscosity properties. These principles form a paradigm to understand the structure-properties-function relationship for therapeutic PMOs to advance the design of new RNA-mimic-based drugs.

Unfolding and Aggregation Pathways of Variable Domains from Immunoglobulin Light Chains

Light chain amyloidosis is the most common form of systemic amyloidosis. This disease is caused by the formation and deposition of amyloid fibers made from immunoglobulin light chains. Environmental conditions such as pH and temperature can affect protein structure and induce the development of these fibers. Several studies have shed light on the native state, stability, dynamics, and final amyloid state of these proteins; however, the initiation process and the fibril formation pathway remain poorly understood structurally and kinetically. To study this, we analyzed the unfolding and aggregation process of the 6aJL2 protein under acidic conditions, with temperature changes, and upon mutation, using biophysical and computational techniques. Our results suggest that the differences in amyloidogenicity displayed by 6aJL2 under these conditions are caused by traversing different aggregation pathways, including unfolded intermediates and the formation of oligomers.

The dynamic nature of netrin-1 and the structural basis for glycosaminoglycan fragment-induced filament formation

Netrin-1 is a bifunctional chemotropic guidance cue that plays key roles in diverse cellular processes including axon pathfinding, cell migration, adhesion, differentiation, and survival. Here, we present a molecular understanding of netrin-1 mediated interactions with glycosaminoglycan chains of diverse heparan sulfate proteoglycans (HSPGs) and short heparin oligosaccharides. Whereas interactions with HSPGs act as platform to co-localise netrin-1 close to the cell surface, heparin oligosaccharides have a significant impact on the highly dynamic behaviour of netrin-1. Remarkably, the monomer-dimer equilibrium of netrin-1 in solution is abolished in the presence of heparin oligosaccharides and replaced with highly hierarchical and distinct super assemblies leading to unique, yet unknown netrin-1 filament formation. In our integrated approach we provide a molecular mechanism for the filament assembly which opens fresh paths towards a molecular understanding of netrin-1 functions.

Emergence of an Auxin Sensing Domain in Plant-Associated Bacteria

Bacteria have evolved a sophisticated array of signal transduction systems that allow them to adapt their physiology and metabolism to changing environmental conditions. Typically, these systems recognize signals through dedicated ligand binding domains (LBDs) to ultimately trigger a diversity of physiological responses. Nonetheless, an increasing number of reports reveal that signal transduction receptors also bind antagonists to inhibit responses mediated by agonists. The mechanisms by which antagonists block the downstream signaling cascade remain largely unknown. To advance our knowledge in this field, we used the LysR-type transcriptional regulator AdmX as a model. AdmX activates the expression of an antibiotic biosynthetic cluster in the rhizobacterium Serratia plymuthica. AdmX specifically recognizes the auxin phytohormone indole-3-acetic acid (IAA) and its biosynthetic intermediate indole-3-pyruvic acid (IPA) as signals. However, only IAA, but not IPA, was shown to regulate antibiotic production in S. plymuthica. Here, we report the high-resolution structures of the LBD of AdmX in complex with IAA and IPA. We found that IAA and IPA compete for binding to AdmX. Although IAA and IPA binding does not alter the oligomeric state of AdmX, IPA binding causes a higher degree of compactness in the protein structure. Molecular dynamics simulations revealed significant differences in the binding modes of IAA and IPA by AdmX, and the inspection of the three-dimensional structures evidenced differential agonist- and antagonist-mediated structural changes. Key residues for auxin binding were identified and an auxin recognition motif defined. Phylogenetic clustering supports the recent evolutionary emergence of this motif specifically in plant-associated enterobacteria. IMPORTANCE Although antagonists were found to bind different bacterial signal transduction receptors, we are still at the early stages of understanding the molecular details by which these molecules exert their inhibitory effects. Here, we provide insight into the structural changes resulting from the binding of an agonist and an antagonist to a sensor protein. Our data indicate that agonist and antagonist recognition is characterized by small conformational differences in the LBDs that can be efficiently transmitted to the output domain to modulate the final response. LBDs are subject to strong selective pressures and are rapidly evolving domains. An increasing number of reports support the idea that environmental factors drive the evolution of sensor domains. Given the recent evolutionary history of AdmX homologs, as well as their narrow phyletic distribution within plant-associated bacteria, our results are in accordance with a plant-mediated evolutionary process that resulted in the emergence of receptor proteins that specifically sense auxin phytohormones.

In-solution structure and oligomerization of human histone deacetylase 6 - an integrative approach

Human histone deacetylase 6 (HDAC6) is a structurally unique, multidomain protein implicated in a variety of physiological processes including cytoskeletal remodelling and the maintenance of cellular homeostasis. Our current understanding of the HDAC6 structure is limited to isolated domains, and a holistic picture of the full-length protein structure, including possible domain interactions, is missing. Here, we used an integrative structural biology approach to build a solution model of HDAC6 by combining experimental data from several orthogonal biophysical techniques complemented by molecular modelling. We show that HDAC6 is best described as a mosaic of folded and intrinsically disordered domains that in-solution adopts an ensemble of conformations without any stable interactions between structured domains. Furthermore, HDAC6 forms dimers/higher oligomers in a concentration-dependent manner, and its oligomerization is mediated via the positively charged N-terminal microtubule-binding domain. Our findings provide the first insights into the structure of full-length human HDAC6 and can be used as a basis for further research into structure function and physiological studies of this unique deacetylase.

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.

Assessment of models for calculating the hydrodynamic radius of intrinsically disordered proteins

Diffusion measurements by pulsed-field gradient NMR and fluorescence correlation spectroscopy can be used to probe the hydrodynamic radius of proteins, which contains information about the overall dimension of a protein in solution. The comparison of this value with structural models of intrinsically disordered proteins is nonetheless impaired by the uncertainty of the accuracy of the methods for computing the hydrodynamic radius from atomic coordinates. To tackle this issue, we here build conformational ensembles of 11 intrinsically disordered proteins that we ensure are in agreement with measurements of compaction by small-angle x-ray scattering. We then use these ensembles to identify the forward model that more closely fits the radii derived from pulsed-field gradient NMR diffusion experiments. Of the models we examined, we find that the Kirkwood-Riseman equation provides the best description of the hydrodynamic radius probed by pulsed-field gradient NMR experiments. While some minor discrepancies remain, our results enable better use of measurements of the hydrodynamic radius in integrative modeling and for force field benchmarking and parameterization.

Less is more: A simple methyl-TROSY based pulse scheme offers improved sensitivity in applications to high molecular weight complexes

The HMQC pulse sequence and variants thereof have been exploited in studies of high molecular weight protein complexes, taking advantage of the fact that fast and slow relaxing magnetization components are sequestered along two distinct magnetization transfer pathways. Despite the simplicity of the HMQC scheme an even shorter version can be designed, based on elimination of the terminal refocusing period, as a further means of increasing signal. Here we present such an experiment, and show that significant sensitivity gains, in some cases by factors of two or more, are realized in studies of proteins varying in molecular masses from 100 kDa to 1 MDa.

Predicting molecular properties of α-synuclein using force fields for intrinsically disordered proteins

Independent force field validation is an essential practice to keep track of developments and for performing meaningful Molecular Dynamics simulations. In this work, atomistic force fields for intrinsically disordered proteins (IDP) are tested by simulating the archetypical IDP α-synuclein in solution for 2.5 μs. Four combinations of protein and water force fields were tested: ff19SB/OPC, ff19SB/TIP4P-D, ff03CMAP/TIP4P-D, and a99SB-disp/TIP4P-disp, with four independent repeat simulations for each combination. We compare our simulations to the results of a 73 μs simulation using the a99SB-disp/TIP4P-disp combination, provided by D. E. Shaw Research. From the trajectories, we predict a range of experimental observations of α-synuclein and compare them to literature data. This includes protein radius of gyration and hydration, intramolecular distances, NMR chemical shifts, and ³ J-couplings. Both ff19SB/TIP4P-D and a99SB-disp/TIP4P-disp produce extended conformational ensembles of α-synuclein that agree well with experimental radius of gyration and intramolecular distances while a99SB-disp/TIP4P-disp reproduces a balanced α-synuclein secondary structure content. It was found that ff19SB/OPC and ff03CMAP/TIP4P-D produce overly compact conformational ensembles and show discrepancies in the secondary structure content compared to the experimental data.

Observing reorientation dynamics with Time-Resolved fluorescence and molecular dynamics in varying periodic boundary conditions

This work presents a combined study of time-resolved fluorescence spectroscopy and all-atom molecular dynamics simulation to investigate periodic boundary conditions' and water models' influence on the orientation dynamics and translational and rotational diffusion of peptides in solution. We have characterized the effects of solvent box size and water model choice on the dynamics of two peptide systems, NATA and WK5. Computationally, translational, and rotational diffusion and internal fluctuations are investigated through all-atom molecular dynamics simulation with two water models and different box sizes. These results are compared with time-resolved fluorescence anisotropy decay (FAD) measurements. The associated time constant and orientation dynamics from FAD measurement along the ¹L_b axis provided baseline data to validate molecular dynamics simulation. The modeling results show that diffusion rates vary roughly in inverse proportion to water model viscosity, as one would expect. Corrections for finite box size are significant for translational diffusion and insignificant for rotational diffusion. This study also finds that internal dynamics described by autocorrelation functions and kinetic network models are relatively insensitive to both box size and water model properties. Our observation suggests that different peptide properties respond differently to a change in simulation conditions.Communicated by Ramaswamy H. Sarma.

A solution structure analysis reveals a bent collagen triple helix in the complement activation recognition molecule mannan-binding lectin

Collagen triple helices are critical in the function of mannan-binding lectin (MBL), an oligomeric recognition molecule in complement activation. The MBL collagen regions form complexes with the serine proteases MASP-1 and MASP-2 in order to activate complement, and mutations lead to common immunodeficiencies. To evaluate their structure-function properties, we studied the solution structures of four MBL-like collagen peptides. The thermal stability of the MBL collagen region was much reduced by the presence of a GQG interruption in the typical (X-Y-Gly)n repeat compared to controls. Experimental solution structural data were collected using analytical ultracentrifugation and small angle X-ray and neutron scattering. As controls, we included two standard Pro-Hyp-Gly collagen peptides (POG)_10-13, as well as three more peptides with diverse (X-Y-Gly)n sequences that represented other collagen features. These data were quantitatively compared with atomistic linear collagen models derived from crystal structures and 12,000 conformations obtained from molecular dynamics simulations. All four MBL peptides were bent to varying degrees up to 85^o in the best-fit molecular dynamics models. The best-fit benchmark peptides (POG)_n were more linear but exhibited a degree of conformational flexibility. The remaining three peptides showed mostly linear solution structures. In conclusion, the collagen helix is not strictly linear, the degree of flexibility in the triple helix depends on its sequence, and the triple helix with the GQG interruption showed a pronounced bend. The bend in MBL GQG peptides resembles the bend in the collagen of complement C1q and may be key for lectin pathway activation.

Exploring Highly Conserved Regions of SARS-CoV-2 Spike S2 Subunit as Targets for Fusion Inhibition Using Chimeric Proteins

Since the beginning of the COVID-19 pandemic, considerable efforts have been made to develop protective vaccines against SARS-CoV-2 infection. However, immunity tends to decline within a few months, and new virus variants are emerging with increased transmissibility and capacity to evade natural or vaccine-acquired immunity. Therefore, new robust strategies are needed to combat SARS-CoV-2 infection. The viral spike composed of S1 and S2 subunits mediates viral attachment and membrane fusion to infect the host cell. In this process, interaction between the highly conserved heptad repeat 1 and 2 regions (HR1 and HR2) of S2 is crucial and for this reason; these regions are promising targets to fight SARS-CoV-2. Here, we describe the design and characterization of chimeric proteins that structurally imitate the S2 HR1 region in a trimeric coiled-coil conformation. We biophysically characterized the proteins and determined their capacity to bind the HR2 region, as well as their inhibitory activity of SARS-CoV-2 infection in vitro. HR1 mimetic proteins showed conformational heterogeneity and a propensity to form oligomers. Moreover, their structure is composed of subdomains with varied stability. Interestingly, the full HR1 proteins showed high affinity for HR2-derived peptides and SARS-CoV-2 inhibitory activity, whereas smaller proteins mimicking HR1 subdomains had a decreased affinity for their complementary HR2 region and did not inhibit the virus. The results provide insight into effective strategies to create mimetic proteins with broad inhibitory activity and therapeutic potential against SARS-CoV-2.

Novel chimeric proteins mimicking SARS-CoV-2 spike epitopes with broad inhibitory activity

SARS-CoV-2 spike (S) protein mediates virus attachment to the cells and fusion between viral and cell membranes. Membrane fusion is driven by mutual interaction between the highly conserved heptad-repeat regions 1 and 2 (HR1 and HR2) of the S2 subunit of the spike. For this reason, these S2 regions are interesting therapeutic targets for COVID-19. Although HR1 and HR2 have been described as transiently exposed during the fusion process, no significant antibody responses against these S2 regions have been reported. Here we designed chimeric proteins that imitate highly stable HR1 helical trimers and strongly bind to HR2. The proteins have broad inhibitory activity against WT B.1 and BA.1 viruses. Sera from COVID-19 convalescent donors showed significant levels of reactive antibodies (IgG and IgA) against the HR1 mimetic proteins, whereas these antibody responses were absent in sera from uninfected donors. Moreover, both inhibitory activity and antigenicity of the proteins correlate positively with their structural stability but not with the number of amino acid changes in their HR1 sequences, indicating a conformational and conserved nature of the involved epitopes. Our results reveal previously undetected spike epitopes that may guide the design of new robust COVID-19 vaccines and therapies.

Tuning gold-based surface functionalization for streptavidin detection: A combined simulative and experimental study

A rationally designed gold-functionalized surface capable of capturing a target protein is presented using the biotin–streptavidin pair as a proof-of-concept. We carried out multiscale simulations to shed light on the binding mechanism of streptavidin on four differently biotinylated surfaces. Brownian Dynamics simulations were used to reveal the preferred initial orientation of streptavidin over the surfaces, whereas classical molecular dynamics was used to refine the binding poses and to investigate the fundamental forces involved in binding, and the binding kinetics. We assessed the binding events and the stability of the streptavidin attachment through a quartz crystal microbalance with dissipation monitoring (QCM-D). The sensing element comprises of biotinylated polyethylene glycol chains grafted on the sensor’s gold surface via thiol-Au chemistry. Finally, we compared the results from experiments and simulations. We found that the confined biotin moieties can specifically capture streptavidin from the liquid phase and provide guidelines on how to exploit the microscopic parameters obtained from simulations to guide the design of further biosensors with enhanced sensitivity.

G-quadruplex DNA: A Longer Story

G-quadruplexes (G4s) are distinctive four-stranded DNA or RNA structures found within cells that are thought to play functional roles in gene regulation and transcription, translation, recombination, and DNA damage/repair. While G4 structures can be uni-, bi-, or tetramolecular with respect to strands, folded unimolecular conformations are most significant in vivo. Unimolecular G4 can potentially form in sequences with runs of guanines interspersed with what will become loops in the folded structure: 5'G_xL_yG_xL_yG_xL_yG_x, where x is typically 2-4 and y is highly variable. Such sequences are highly conserved and specifically located in genomes. In the folded structure, guanines from each run combine to form planar tetrads with four hydrogen-bonded guanine bases; these tetrads stack on one another to produce four strand segments aligned in specific parallel or antiparallel orientations, connected by the loop sequences. Three types of loops (lateral, diagonal, or "propeller") have been identified. The stacked tetrads form a central cavity that features strong coordination sites for monovalent cations that stabilize the G4 structure, with potassium or sodium preferred. A single monomeric G4 typically forms from a sequence containing roughly 20-30 nucleotides. Such short sequences have been the primary focus of X-ray crystallographic or NMR studies that have produced high-resolution structures of a variety of monomeric G4 conformations. These structures are often used as the basis for drug design efforts to modulate G4 function.We believe that the focus on monomeric G4 structures formed by such short sequences is perhaps myopic. Such short sequences for structural studies are often arbitrarily selected and removed from their native genomic sequence context, and then are often changed from their native sequences by base substitutions or deletions intended to optimize the formation of a homogeneous G4 conformation. We believe instead that G-quadruplexes prefer company and that in a longer natural sequence context multiple adjacent G4 units can form to combine into more complex multimeric G4 structures with richer topographies than simple monomeric forms. Bioinformatic searches of the human genome show that longer sequences with the potential for forming multiple G4 units are common. Telomeric DNA, for example, has a single-stranded overhang of hundreds of nucleotides with the requisite repetitive sequence with the potential for formation of multiple G4s. Numerous extended promoter sequences have similar potentials for multimeric G4 formation. X-ray crystallography and NMR methods are challenged by these longer sequences (>30 nt), so other tools are needed to explore the possible multimeric G4 landscape. We have implemented an integrated structural biology approach to address this challenge. This approach integrates experimental biophysical results with atomic-level molecular modeling and molecular dynamics simulations that provide quantitatively testable model structures. In every long sequence we have studied so far, we found that multimeric G4 structures readily form, with a surprising diversity of structures dependent on the exact native sequence used. In some cases, stable hairpin duplexes form along with G4 units to provide an even richer landscape. This Account provides an overview of our approach and recent progress and provides a new perspective on the G-quadruplex folding landscape.

Transient Prenyltransferase-Cyclase Association in Fusicoccadiene Synthase, an Assembly-Line Terpene Synthase

Fusicoccadiene synthase from the fungus Phomopsis amygdali (PaFS) is an assembly-line terpene synthase that catalyzes the first two steps in the biosynthesis of Fusiccocin A, a diterpene glycoside. The C-terminal prenyltransferase domain of PaFS catalyzes the condensation of one molecule of C5 dimethylallyl diphosphate and three molecules of C5 isopentenyl diphosphate to form C20 geranylgeranyl diphosphate, which then transits to the cyclase domain for cyclization to form fusicoccadiene. Previous structural studies of PaFS using electron microscopy (EM) revealed a central octameric prenyltransferase core with eight cyclase domains tethered in random distal positions through flexible 70-residue linkers. However, proximal prenyltransferase-cyclase configurations could be captured by covalent cross-linking and observed by cryo-EM and mass spectrometry. Here, we use cryo-EM to show that proximally configured prenyltransferase-cyclase complexes are observable even in the absence of covalent cross-linking; moreover, such complexes can involve multiple cyclase domains. A conserved basic patch on the prenyltransferase domain comprises the primary touchpoint with the cyclase domain. These results support a model for transient prenyltransferase-cyclase association in which the cyclase domains of PaFS are in facile equilibrium between proximal associated and random distal positions relative to the central prenyltransferase octamer. The results of biophysical measurements using small-angle X-ray scattering, analytical ultracentrifugation, dynamic light scattering, and size-exclusion chromatography in-line with multi-angle light scattering are consistent with this model. This model accordingly provides a framework for understanding substrate transit between the prenyltransferase and cyclase domains as well as the cooperativity observed for geranylgeranyl diphosphate cyclization.