Large Supramolecular Complexes Research Articles

Assessing the domain-based local pair natural orbital (DLPNO) approximation for non-covalent interactions in sizable supramolecular complexes.

The titular domain-based local pair natural orbital (DLPNO) approximation is the most widely used method for extending correlated wave function models to large molecular systems, yet its fidelity for intermolecular interaction energies in large supramolecular complexes has not been thoroughly vetted. Non-covalent interactions are sensitive to tails of the electron density and involve nonlocal dispersion that is discarded or approximated if the screening of pair natural orbitals (PNOs) is too aggressive. Meanwhile, the accuracy of the DLPNO approximation is known to deteriorate as molecular size increases. Here, we test the DLPNO approximation at the level of second-order Møller-Plesset perturbation theory (MP2) and coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)] for a variety of large supramolecular complexes. DLPNO-MP2 interaction energies are within 3% of canonical values for small dimers with ≲10 heavy atoms, but for larger systems, the DLPNO approximation is often quite poor unless the results are extrapolated to the canonical limit where the threshold for discarding PNOs is taken to zero. Counterpoise correction proves to be essential in reducing errors with respect to canonical results. For a sequence of nanoscale graphene dimers up to (C96H24)2, extrapolated DLPNO-MP2 interaction energies agree with canonical values to within 1%, independent of system size, provided that the basis set does not contain diffuse functions; these cause the DLPNO approximation to behave erratically, such that results cannot be extrapolated in a meaningful way. DLPNO-CCSD(T) calculations are typically performed using looser PNO thresholds as compared to DLPNO-MP2, but this significantly impacts accuracy for large supramolecular complexes. Standard DLPNO-CCSD(T) settings afford errors of 2-6 kcal/mol for dimers involving coronene (C24H12) and circumcoronene (C54H18), even at the DLPNO-CCSD(T1) level.

Native architecture of a human GBP1 defense complex for cell-autonomous immunity to infection.

All living organisms deploy cell-autonomous defenses to combat infection. In plants and animals, large supramolecular complexes often activate immune proteins for protection. In this work, we resolved the native structure of a massive host-defense complex that polymerizes 30,000 guanylate-binding proteins (GBPs) over the surface of gram-negative bacteria inside human cells. Construction of this giant nanomachine took several minutes and remained stable for hours, required guanosine triphosphate hydrolysis, and recruited four GBPs plus caspase-4 and Gasdermin D as a cytokine and cell death immune signaling platform. Cryo-electron tomography suggests that GBP1 can adopt an extended conformation for bacterial membrane insertion to establish this platform, triggering lipopolysaccharide release that activated coassembled caspase-4. Our "open conformer" model provides a dynamic view into how the human GBP1 defense complex mobilizes innate immunity to infection.

Utility of methyl side chain probes for solution NMR studies of large proteins

Selective isotopic labeling of methyl side chain groups in proteins and other biomolecules, combined with advances in perdeuteration, new pulse sequences, and high field spectrometers with cryogenic probes, has revolutionized the field of solution nuclear magnetic resonance (NMR) spectroscopy by enabling characterization of macromolecular systems with molecular weights above 1 MDa in their native aqueous environment. This tutorial provides a basic overview for how modern NMR spectroscopists can utilize methyl side chain probes to study their system of interest. The advantages and limitations of methyl side chain probes are discussed. In addition, the preparation of selectively 13C-methyl labeled recombinant protein samples, strategies for manual and automated methyl NMR resonance assignment, and the application of methyl probes for characterization of dynamics and conformational exchange are discussed. A sneak preview for ways in which methyl probes are expected to continue to advance the field of biomolecular NMR towards new horizons in solution studies of large supramolecular complexes is also presented.

SNARE Modulators and SNARE Mimetic Peptides.

The soluble N-ethylmaleimide-sensitive factor (NSF) attachment protein (SNAP) receptor (SNARE) proteins play a central role in most forms of intracellular membrane trafficking, a key process that allows for membrane and biocargo shuffling between multiple compartments within the cell and extracellular environment. The structural organization of SNARE proteins is relatively simple, with several intrinsically disordered and folded elements (e.g., SNARE motif, N-terminal domain, transmembrane region) that interact with other SNAREs, SNARE-regulating proteins and biological membranes. In this review, we discuss recent advances in the development of functional peptides that can modify SNARE-binding interfaces and modulate SNARE function. The ability of the relatively short SNARE motif to assemble spontaneously into stable coiled coil tetrahelical bundles has inspired the development of reduced SNARE-mimetic systems that use peptides for biological membrane fusion and for making large supramolecular protein complexes. We evaluate two such systems, based on peptide-nucleic acids (PNAs) and coiled coil peptides. We also review how the self-assembly of SNARE motifs can be exploited to drive on-demand assembly of complex re-engineered polypeptides.

Protein dynamics developments for the large scale and cryoEM: case study of ProDy 2.0.

Cryo-electron microscopy (cryoEM) has become a well established technique with the potential to produce structures of large and dynamic supramolecular complexes that are not amenable to traditional approaches for studying structure and dynamics. The size and low resolution of such molecular systems often make structural modelling and molecular dynamics simulations challenging and computationally expensive. This, together with the growing wealth ofstructural data arising from cryoEM and other structural biology methods, hasdriven a trend in the computational biophysics community towards the development of new pipelines for analysing global dynamics using coarse-grained models and methods. At the centre of this trend has been a return to elastic network models, normal mode analysis (NMA) and ensemble analyses such as principal component analysis, and the growth of hybrid simulation methodologies that make use of them. Here, this field is reviewed with a focus on ProDy, the Python application programming interface for protein dynamics, which has been developed over the last decade. Two key developments in this area are highlighted: (i) ensemble NMA towards extracting and comparing the signature dynamics of homologous structures, aided by the recent SignDy pipeline, and (ii) pseudoatom fitting for more efficient global dynamics analyses of large and low-resolution supramolecular assemblies from cryoEM, revisited in the CryoDy pipeline. It is believed that such a renewal and extension of old models and methods in new pipelines will be critical for driving the field forward into the next cryoEM revolution.

Beyond-dipole van der Waals contributions within the many-body dispersion framework

By introducing a suitable range-separation of the Coulomb coupling in analogy to Ambrosetti et al (2014 J. Chem. Phys. 140 18A508), here we extend the many-body dispersion approach to include beyond-dipole van der Waals (vdW) interactions at a full many-body level, in combination with semi-local density functional theory. A reciprocal-space implementation is further introduced in order to efficiently treat periodic systems. Consistent reliability is found from molecular dimers to large supramolecular complexes and two-dimensional systems. The large weight of both many-body effects and multipolar terms illustrates how a correct description of vdW forces in large-scale systems requires full account of both contributions, beyond standard pairwise dipolar approaches.

Interaction of HMGB-domain proteins from differentiated cells with DNA modified by cisplatin

Using UV absorption spectroscopy and circular dichroism, we studied the interaction of DNA-cisplatin complexes with nonhistone chromosomal proteins HMGB1 and HMGB2. It was shown that the presence of platinum ions affected the ability of HMGB2 to form large supramolecular complexes with DNA. The role of C-terminal domain of the HMGB1/2 proteins is discussed.

In cellulo FRET-FLIM and single molecule tracking reveal the supra-molecular organization of the pyoverdine bio-synthetic enzymes in Pseudomonas aeruginosa.

The bio-synthesis of pyoverdine (PVD) in Pseudomonas aeruginosa involves multiple enzymatic steps including the action of non-ribosomal peptide synthetases (NRPSs). One hallmark of NRPS is their ability to make usage of non-proteinogenic amino-acids synthesized by co-expressed accessory enzymes. It is generally proposed that different enzymes of a secondary metabolic pathway assemble into large supra-molecular complexes. However, evidence for the assembly of sequential enzymes in the cellular context is sparse. Here, we used in cellulo single-molecule tracking and Förster resonance energy transfer measured by fluorescence lifetime microscopy (FRET-FLIM) to explore the spatial partitioning of the ornithine hydroxylase PvdA and its interactions with NRPS. We found PvdA was mostly diffusing bound to large complexes in the cytoplasm with a small exchangeable trapped fraction. FRET-FLIM clearly showed that PvdA is physically interacting with PvdJ, PvdI, PvdL, and PvdD, the four NRPS involved in the PVD pathway in Pseudomonas aeruginosa PAO1. The binding modes of PvdA were strikingly different according to the NRPS it is interacting with, suggesting that PvdA binding sites have co-evolved with the enzymatic active sites of NRPS. Our data provide evidence for strongly organized multi-enzymatic complexes responsible for the bio-synthesis of PVD and illustrate how binding sites have evolved to finely control the co-localization of sequential enzymes and promote metabolic pathway efficiency.

Supramolecular Complexes in Cell Death and Inflammation and Their Regulation by Autophagy.

Signaling activation is a tightly regulated process involving myriad posttranslational modifications such as phosphorylation/dephosphorylation, ubiquitylation/deubiquitylation, proteolytical cleavage events as well as translocation of proteins to new compartments within the cell. In addition to each of these events potentially regulating individual proteins, the assembly of very large supramolecular complexes has emerged as a common theme in signal transduction and is now known to regulate many signaling events. This is particularly evident in pathways regulating both inflammation and cell death/survival. Regulation of the assembly and silencing of these complexes plays important roles in immune signaling and inflammation and the fate of cells to either die or survive. Here we will give a summary of some of the better studied supramolecular complexes involved in inflammation and cell death, particularly with a focus on diseases caused by their autoactivation and the role autophagy either plays or may be playing in their regulation.

Polymorphisms of a Collagen-Like Adhesin Contributes to Legionella pneumophila Adhesion, Biofilm Formation Capacity and Clinical Prevalence.

Legionellosis is a severe respiratory illness caused by the inhalation of aerosolized water droplets contaminated with the opportunistic pathogen Legionella pneumophila. The ability of L. pneumophila to produce biofilms has been associated with its capacity to colonize and persist in human-made water reservoirs and distribution systems, which are the source of legionellosis outbreaks. Nevertheless, the factors that mediate L. pneumophila biofilm formation are largely unknown. In previous studies we reported that the adhesin Legionella collagen-like protein (Lcl), is required for auto-aggregation, attachment to multiple surfaces and the formation of biofilms. Lcl structure contains three distinguishable regions: An N-terminal region with a predicted signal sequence, a central region containing tandem collagen-like repeats (R-domain) and a C-terminal region (C-domain) with no significant homology to other known proteins. Lcl R-domain encodes tandem repeats of the collagenous tripeptide Gly-Xaa-Yaa (GXY), a motif that is key for the molecular organization of mammalian collagen and mediates the binding of collagenous proteins to different cellular and environmental ligands. Interestingly, Lcl is polymorphic in the number of GXY tandem repeats. In this study, we combined diverse biochemical, genetic, and cellular approaches to determine the role of Lcl domains and GXY repeats polymorphisms on the structural and functional properties of Lcl, as well as on bacterial attachment, aggregation and biofilm formation. Our results indicate that the R-domain is key for assembling Lcl collagenous triple-helices and has a more preponderate role over the C-domain in Lcl adhesin binding properties. We show that Lcl molecules oligomerize to form large supramolecular complexes to which both, R and C-domains are required. Furthermore, we found that the number of GXY tandem repeats encoded in Lcl R-domain correlates positively with the binding capabilities of Lcl and with the attachment and biofilm production capacity of L. pneumophila strains. Accordingly, the number of GXY tandem repeats in Lcl influences the clinical prevalence of L. pneumophila strains. Therefore, the number of Lcl tandem repeats could be considered as a potential predictor for virulence in L. pneumophila isolates.

Ubiquitin-Independent Degradation of Proteins in Proteasomes

Proteasomes are large supramolecular protein complexes present in all prokaryotic and eukaryotic cells, where they perform targeted degradation of intracellular proteins. Until recently, it was generally accepted that prior to proteolytic degradation in proteasomes the proteins had to be targeted by ubiquitination: ATP-dependent attachment of (typically four sequential) residues of the low-molecular protein, ubiquitin, which involves the ubiquitin-activating enzyme, ubiquitin-conjugating enzyme, and ubiquitin ligase. Cytoplasmic and nucleoplasmic proteins labeled in this way are then digested in 26S proteasomes. However, it becomes increasingly clear that using this route the cell eliminates only a part of unwanted proteins. Many proteins can be cleaved by the 20S proteasome in an ATP-independent manner and without previous ubiquitination. Ubiquitin-independent degradation of proteins in proteasomes is a relatively new area of studies of the role of the ubiquitin-proteasome system. However, recent data obtained in this direction already correct existing concepts about proteasomal degradation of proteins and its regulation. Ubiquitin-independent proteasome degradation needs the main structural precondition in proteins: the presence of unstructured regions in the amino acid sequences that provide interaction with the proteasome. Taking into consideration that in humans almost half of all genes encode proteins that contain a certain proportion of intrinsically disordered regions, it appears that the list of proteins undergoing ubiquitin-independent degradation will demonstrate a further increase. Since 26S proteasomes account for only 30% of the total proteasome content in mammalian cells, most of the proteasomes exist in the form of 20S complexes. The latter suggests that ubiquitin-independent proteolysis performed by the 20S proteasome is a natural process of removing damaged proteins from the cell and maintaining a constant level of intrinsically disordered proteins. In this case, the functional overload of proteasomes in aging and/or other types of pathological processes, if it is not accompanied by triggering more radical mechanisms for the elimination of damaged proteins, organelles, and whole cells, has the most serious consequences for the whole organism.

Structures and properties of large supramolecular coordination complexes predicted with the generalized energy-based fragmentation method.

The generalized energy-based fragmentation (GEBF) method has been extended to facilitate ab initio calculations of large supramolecular coordination complexes. For metal-containing coordination complexes, a special fragmentation scheme is proposed for GEBF calculations, in which coordinate bonds between metal ions and ligands are kept intact, and only single covalent bonds in organic ligands are cut into fragments. A simple strategy is exploited for the determination of the ground-state spin multiplicity of each metal ion so that the total spin of all metal-containing subsystems is assigned automatically. With this fragmentation scheme, the GEBF method is demonstrated to provide reliable energies, optimized geometry, nuclear magnetic resonance (NMR) properties and the infrared spectrum for a medium-sized supramolecular coordination complex, which are very consistent with those from the full-system quantum chemistry calculations. The GEBF method is then applied to two large supramolecular coordination complexes to illustrate its capability. For the trimetallic coordination complex Fe2Zn2(RuL2)2 (with 618 atoms), the calculated 1H chemical shifts from GEBF calculations with the B97-2 functional can account well for the experimental NMR spectrum. For the cage-guest complex Pd4L8(BF4-)3, the computed infrared spectrum obtained with the GEBF-M06-2X method can help assign the experimental peaks to the corresponding vibrational motions. The GEBF method combining with advanced electronic structure methods is expected to be a useful tool to understand and interpret structural and spectroscopic information of various supramolecular coordination complexes.

Ubiquitin-independent protein degradation in proteasomes

Proteasomes are large supramolecular protein complexes present in all prokaryotic and eukaryotic cells, where they perform targeted degradation of intracellular proteins. Until recently, it was generally accepted that prior proteolytic degradation in proteasomes the proteins had to be targeted by ubiquitination: the ATP-dependent addition of (typically four sequential) residues of the low-molecular ubiquitin protein, involving the ubiquitin-activating enzyme, ubiquitin-conjugating enzyme and ubiquitin ligase. The cytoplasm and nucleoplasm proteins labeled in this way are then digested in 26S proteasomes. However, in recent years it has become increasingly clear that using this route the cell eliminates only a part of unwanted proteins. Many proteins can be cleaved by the 20S proteasome in an ATP-independent manner and without previous ubiquitination. Ubiquitin-independent protein degradation in proteasomes is a relatively new area of studies of the role of the ubiquitin-proteasome system. However, recent data obtained in this direction already correct existing concepts about proteasomal degradation of proteins and its regulation. Ubiquitin-independent proteasome degradation needs the main structural precondition in proteins: the presence of unstructured regions in the amino acid sequences that provide interaction with the proteasome. Taking into consideration that in humans almost half of all genes encode proteins that contain a certain proportion of intrinsically disordered regions, it appears that the list of proteins undergoing ubiquitin-independent degradation will demonstrate further increase. Since 26S of proteasomes account for only 30% of the total proteasome content in mammalian cells, most of the proteasomes exist in the form of 20S complexes. The latter suggests that ubiquitin-independent proteolysis performed by the 20S proteasome is a natural process of removing damaged proteins from the cell and maintaining a constant level of intrinsically disordered proteins. In this case, the functional overload of proteasomes in aging and/or other types of pathological processes, if it is not accompanied by triggering more radical mechanisms for the elimination of damaged proteins, organelles and whole cells, has the most serious consequences for the whole organism.

DLPNO-CCSD(T) scaled methods for the accurate treatment of large supramolecular complexes.

In this work, we present scaled variants of the DLPNO-CCSD(T) method, dubbed as (LS)DLPNO-CCSD(T) and (NS)DLPNO-CCSD(T), to obtain accurate interaction energies in supramolecular complexes governed by noncovalent interactions. The novel scaled schemes are based on the linear combination of the DLPNO-CCSD(T) correlation energies calculated with the standard (LoosePNO and NormalPNO) and modified (Loose2PNO and Normal2PNO) DLPNO-CCSD(T) accuracy levels. The scaled DLPNO-CCSD(T) variants provide nearly TightPNO accuracy, which is essential for the quantification of weak noncovalent interactions, with a noticeable saving in computational cost. Importantly, the accuracy of the proposed schemes is preserved irrespective of the nature and strength of the supramolecular interaction. The (LS)DLPNO-CCSD(T) and (NS)DLPNO-CCSD(T) protocols have been used to study in depth the role of the CH-π versus π-π interactions in the supramolecular complex formed by the electron-donor truxene-tetrathiafulvalene (truxTTF) and the electron-acceptor hemifullerene (C30 H12 ). (NS)DLPNO-CCSD(T)/CBS calculations clearly reveal the higher stability of staggered (dominated by CH-π interactions) versus bowl-in-bowl (dominated by π-π interactions) arrangements in the truxTTF•C30 H12 heterodimer. Hemifullerene and similar carbon-based buckybowls are therefore expected to self-assemble with donor compounds in a richer way other than the typical concave-convex π-π arrangement found in fullerene-based aggregates. © 2017 Wiley Periodicals, Inc.

Nanoscale \u03c0\u2013\u03c0 stacked molecules are bound by collective charge fluctuations

Non-covalent π−π interactions are central to chemical and biological processes, yet the full understanding of their origin that would unite the simplicity of empirical approaches with the accuracy of quantum calculations is still missing. Here we employ a quantum-mechanical Hamiltonian model for van der Waals interactions, to demonstrate that intermolecular electron correlation in large supramolecular complexes at equilibrium distances is appropriately described by collective charge fluctuations. We visualize these fluctuations and provide connections both to orbital-based approaches to electron correlation, as well as to the simple London pairwise picture. The reported binding energies of ten supramolecular complexes obtained from the quantum-mechanical fluctuation model joined with density functional calculations are within 5% of the reference energies calculated with the diffusion quantum Monte-Carlo method. Our analysis suggests that π−π stacking in supramolecular complexes can be characterized by strong contributions to the binding energy from delocalized, collective charge fluctuations—in contrast to complexes with other types of bonding.

ChemInform Abstract: Small Atomic Orbital Basis Set First‐Principles Quantum Chemical Methods for Large Molecular and Periodic Systems: A Critical Analysis of Error Sources

AbstractReview: 123 refs.

Small Atomic Orbital Basis Set First-Principles Quantum Chemical Methods for Large Molecular and Periodic Systems: A Critical Analysis of Error Sources.

In quantum chemical computations the combination of Hartree–Fock or a density functional theory (DFT) approximation with relatively small atomic orbital basis sets of double‐zeta quality is still widely used, for example, in the popular B3LYP/6‐31G* approach. In this Review, we critically analyze the two main sources of error in such computations, that is, the basis set superposition error on the one hand and the missing London dispersion interactions on the other. We review various strategies to correct those errors and present exemplary calculations on mainly noncovalently bound systems of widely varying size. Energies and geometries of small dimers, large supramolecular complexes, and molecular crystals are covered. We conclude that it is not justified to rely on fortunate error compensation, as the main inconsistencies can be cured by modern correction schemes which clearly outperform the plain mean‐field methods.

Novel approaches to examine the regulation of voltage-gated calcium channels in the heart.

The cardiac L-type Ca(2+) channel plays a key role in cardiac excitation-contraction coupling, action potential duration, and gene expression. Abnormalities in CaV1.2 function, including increased long-opening-mode gating and blunted adrenergic responsiveness, are associated with heart failure and hypertrophy. The increased activation of CaV1.2, in turn, triggers Ca(2+) -responsive signaling pathways, which contribute to the pathogenesis of heart failure and hypertrophy. CaV1.2 in heart is associated with large supramolecular complexes that impact on channel trafficking, localization, turnover, and function. Much of the prevailing dogma relating to mechanisms underlying CaV1.2 trafficking and modulation is derived from studies using recombinant channels reconstituted in heterologous expression systems. However, recent results using knock-in mice indicate that several long-standing "facts" about CaV1.2 regulation derived from heterologous expression studies are not replicated in native heart, emphasizing the critical need for mechanistic studies in the context of actual cardiomyocytes. In this review, we discuss the use of the use of knockin and knockout mice as well as new tools, including doxycycline-induced expression of informative α1C mutants within cardiomyocytes to probe adrenergic regulation of CaV1.2.

Accurate Treatment of Large Supramolecular Complexes by Double-Hybrid Density Functionals Coupled with Nonlocal van der Waals Corrections.

In this work, we present a thorough assessment of the performance of some representative double-hybrid density functionals (revPBE0-DH-NL and B2PLYP-NL) as well as their parent hybrid and GGA counterparts, in combination with the most modern version of the nonlocal (NL) van der Waals correction to describe very large weakly interacting molecular systems dominated by noncovalent interactions. Prior to the assessment, an accurate and homogeneous set of reference interaction energies was computed for the supramolecular complexes constituting the L7 and S12L data sets by using the novel, precise, and efficient DLPNO-CCSD(T) method at the complete basis set limit (CBS). The correction of the basis set superposition error and the inclusion of the deformation energies (for the S12L set) have been crucial for obtaining precise DLPNO-CCSD(T)/CBS interaction energies. Among the density functionals evaluated, the double-hybrid revPBE0-DH-NL and B2PLYP-NL with the three-body dispersion correction provide remarkably accurate association energies very close to the chemical accuracy. Overall, the NL van der Waals approach combined with proper density functionals can be seen as an accurate and affordable computational tool for the modeling of large weakly bonded supramolecular systems.

DNA-Directed Assembly of Capture Tools for Constitutional Studies of Large Protein Complexes.

Large supramolecular protein complexes, such as the molecular machinery involved in gene regulation, cell signaling, or cell division, are key in all fundamental processes of life. Detailed elucidation of structure and dynamics of such complexes can be achieved by reverse-engineering parts of the complexes in order to probe their interactions with distinctive binding partners in vitro. The exploitation of DNA nanostructures to mimic partially assembled supramolecular protein complexes in which the presence and state of two or more proteins are decisive for binding of additional building blocks is reported here. To this end, four-way DNA Holliday junction motifs bearing a fluorescein and a biotin tag, for tracking and affinity capture, respectively, are site-specifically functionalized with centromeric protein (CENP) C and CENP-T. The latter serves as baits for binding of the so-called KMN component, thereby mimicking early stages of the assembly of kinetochores, structures that mediate and control the attachment of microtubules to chromosomes in the spindle apparatus. Results from pull-down experiments are consistent with the hypothesis that CENP-C and CENP-T may bind cooperatively to the KMN network.