Hybrid Approaches to Structural Characterization of Conformational Ensembles of Complex Macromolecular Systems Combining NMR Residual Dipolar Couplings and Solution X-ray Scattering

Elucidation of the Mechanisms of Inter-domain Coupling in the Monomeric State of Enzyme I by High-pressure NMR

The catalytic cycle of Enzyme I (EI), a phosphotransferase enzyme responsible for converting phosphoenolpyruvate (PEP) into pyruvate, is characterized by a series of local and global conformational rearrangements. This multistep process includes a monomer-to-dimer transition, followed by an open-to-closed rearrangement of the dimeric complex upon PEP binding. In the present study, we investigate the thermodynamics of EI dimerization using a range of high-pressure solution NMR techniques complemented by SAXS experiments. 1H-15N TROSY and 1H-13C methyl TROSY NMR spectra combined with 15N relaxation measurements revealed that a native-like engineered variant of full-length EI fully dissociates into stable monomeric state above 1.5 kbar. Conformational ensembles of EI monomeric state were generated via a recently developed protocol combining coarse-grained molecular simulations with experimental backbone residual dipolar coupling measurements. Analysis of the structural ensembles provided detailed insights into the molecular mechanisms driving formation of the catalytically competent dimeric state, and reveals that each step of EI catalytical cycle is associated with a significant reduction in either inter- or intra-domain conformational entropy. Altogether, this study completes a large body work conducted by our group on EI and establishes a comprehensive structural and dynamical description of the catalytic cycle of this prototypical multidomain, oligomeric enzyme.

Solution Structure Ensembles of the Open and Closed Forms of the ∼130 kDa Enzyme I via AlphaFold Modeling, Coarse Grained Simulations, and NMR

Large-scale interdomain rearrangements are essential to protein function, governing the activity of large enzymes and molecular machineries. Yet, obtaining an atomic-resolution understanding of how the relative domain positioning is affected by external stimuli is a hard task in modern structural biology. Here, we show that combining structural modeling by AlphaFold2 with coarse-grained molecular dynamics simulations and NMR residual dipolar coupling data is sufficient to characterize the spatial domain organization of bacterial enzyme I (EI), a ∼130 kDa multidomain oligomeric protein that undergoes large-scale conformational changes during its catalytic cycle. In particular, we solve conformational ensembles for EI at two different experimental temperatures and demonstrate that a lower temperature favors sampling of the catalytically competent closed state of the enzyme. These results suggest a role for conformational entropy in the activation of EI and demonstrate the ability of our protocol to detect and characterize the effect of external stimuli (such as mutations, ligand binding, and post-translational modifications) on the interdomain organization of multidomain proteins. We expect the ensemble refinement protocol described here to be easily transferrable to the investigation of the structure and dynamics of other uncharted multidomain systems and have assembled a Google Colab page (https://potoyangroup.github.io/Seq2Ensemble/) to facilitate implementation of the presented methodology elsewhere.

NMR spectroscopy, excited states and relevance to problems in cell biology - transient pre-nucleation tetramerization of huntingtin and insights into Huntington's disease

Solution nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for analyzing three-dimensional structure and dynamics of macromolecules at atomic resolution. Recent advances have exploited the unique properties of NMR in exchanging systems to detect, characterize and visualize excited sparsely populated states of biological macromolecules and their complexes, which are only transient. These states are invisible to conventional biophysical techniques, and play a key role in many processes, including molecular recognition, protein folding, enzyme catalysis, assembly and fibril formation. All the NMR techniques make use of exchange between sparsely populated NMR-invisible and highly populated NMR-visible states to transfer a magnetization property from the invisible state to the visible one where it can be easily detected and quantified. There are three classes of NMR experiments that rely on differences in distance, chemical shift or transverse relaxation (molecular mass) between the NMR-visible and -invisible species. Here, I illustrate the application of these methods to unravel the complex mechanism of sub-millisecond pre-nucleation oligomerization of the N-terminal region of huntingtin, encoded by exon-1 of the huntingtin gene, where CAG expansion leads to Huntington's disease, a fatal autosomal-dominant neurodegenerative condition. I also discuss how inhibition of tetramerization blocks the much slower (by many orders of magnitude) process of fibril formation.

Large Chaperone Complexes Through the Lens of Nuclear Magnetic Resonance Spectroscopy

Molecular chaperones are the guardians of the proteome inside the cell. Chaperones recognize and bind unfolded or misfolded substrates, thereby preventing further aggregation; promoting correct protein folding; and, in some instances, even disaggregating already formed aggregates. Chaperones perform their function by means of an array of weak protein-protein interactions that take place over a wide range of timescales and are therefore invisible to structural techniques dependent upon the availability of highly homogeneous samples. Nuclear magnetic resonance (NMR) spectroscopy, however, is ideally suited to study dynamic, rapidly interconverting conformational states and protein-protein interactions in solution, even if these involve a high-molecular-weight component. In this review, we give a brief overview of the principles used by chaperones to bind their client proteins and describe NMR methods that have emerged as valuable tools to probe chaperone-substrate and chaperone-chaperone interactions. We then focus on a few systems for which the application of these methods has greatly increased our understanding of the mechanisms underlying chaperone functions. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Protein of a thousand faces: The tumor-suppressive and oncogenic responses of p53

The p53 protein is a pleiotropic regulator working as a tumor suppressor and as an oncogene. Depending on the cellular insult and the mutational status, p53 may trigger opposing activities such as cell death or survival, senescence and cell cycle arrest or proliferative signals, antioxidant or prooxidant activation, glycolysis, or oxidative phosphorylation, among others. By augmenting or repressing specific target genes or directly interacting with cellular partners, p53 accomplishes a particular set of activities. The mechanism in which p53 is activated depends on increased stability through post-translational modifications (PTMs) and the formation of higher-order structures (HOS). The intricate cell death and metabolic p53 response are reviewed in light of gaining stability via PTM and HOS formation in health and disease.

Energetics and J-coupling constants for Ala, Gly, and Val peptides demonstrated using ABEEM polarizable force field in vacuo and an aqueous solution

The development of an atom-bond electronegativity equalisation method at the σπ-level (ABEEM) polarisable force field (PFF) for peptides is presented. ABEEM PFF utilises a fluctuating charge model to explicitly describe...

Advances in NMR Spectroscopy of Weakly Aligned Biomolecular Systems

The measurement and application of residual dipolar couplings (RDCs) in solution NMR studies of biological macromolecules has become well established over the past quarter of a century. Numerous methods for generating the requisite anisotropic orientational molecular distribution have been demonstrated, each with its specific strengths and weaknesses. In parallel, an enormous number of pulse schemes have been introduced to measure the many different types of RDCs, ranging from the most widely measured backbone amide ¹⁵N-¹H RDCs, to ¹H-¹H RDCs and couplings between low-γ nuclei. Applications of RDCs range from structure validation and refinement to the determination of relative domain orientations, the measurement of backbone and domain motions, and de novo structure determination. Nevertheless, it appears that the power of the RDC methodology remains underutilized. This review aims to highlight the practical aspects of sample preparation and RDC measurement while describing some of the most straightforward applications that take advantage of the exceptionally precise information contained in such data. Some emphasis will be placed on more recent developments that enable the accurate measurement of RDCs in larger systems, which is key to the ongoing shift in focus of biological NMR spectroscopy from structure determination toward gaining improved understanding of how molecular flexibility drives protein function.

N 6-methyladenosine binding induces a metal-centered rearrangement that activates the human RNA demethylase Alkbh5

Alkbh5 catalyzes demethylation of the N ⁶-methyladenosine (m⁶A), an epigenetic mark that controls several physiological processes including carcinogenesis and stem cell differentiation. The activity of Alkbh5 comprises two coupled reactions. The first reaction involves decarboxylation of α-ketoglutarate (αKG) and formation of a Fe⁴⁺═O species. This oxyferryl intermediate oxidizes the m⁶A to reestablish the canonical base. Despite coupling between the two reactions being required for the correct Alkbh5 functioning, the mechanisms linking dioxygen activation to m⁶A binding are not fully understood. Here, we use solution NMR to investigate the structure and dynamics of apo and holo Alkbh5. We show that binding of m⁶A to Alkbh5 induces a metal-centered rearrangement of αKG that increases the exposed area of the metal, making it available for binding O₂ Our study reveals the molecular mechanisms underlying activation of Alkbh5, therefore opening new perspectives for the design of novel strategies to control gene expression and cancer progression.

Molecular insights on CALX-CBD12 interdomain dynamics from MD simulations, RDCs, and SAXS

Na⁺/Ca²⁺ exchangers (NCXs) are secondary active transporters that couple the translocation of Na⁺ with the transport of Ca²⁺ in the opposite direction. The exchanger is an essential Ca²⁺ extrusion mechanism in excitable cells. It consists of a transmembrane domain and a large intracellular loop that contains two Ca²⁺-binding domains, CBD1 and CBD2. The two CBDs are adjacent to each other and form a two-domain Ca²⁺ sensor called CBD12. Binding of intracellular Ca²⁺ to CBD12 activates the NCX but inhibits the NCX of Drosophila, CALX. NMR spectroscopy and SAXS studies showed that CALX and NCX CBD12 constructs display significant interdomain flexibility in the apo state but assume rigid interdomain arrangements in the Ca²⁺-bound state. However, detailed structure information on CBD12 in the apo state is missing. Structural characterization of proteins formed by two or more domains connected by flexible linkers is notoriously challenging and requires the combination of orthogonal information from multiple sources. As an attempt to characterize the conformational ensemble of CALX-CBD12 in the apo state, we applied molecular dynamics (MD) simulations, NMR (¹H-¹⁵N residual dipolar couplings), and small-angle x-ray scattering (SAXS) data in a combined strategy to select an ensemble of conformations in agreement with the experimental data. This joint approach demonstrated that CALX-CBD12 preferentially samples closed conformations, whereas the wide-open interdomain arrangement characteristic of the Ca²⁺-bound state is less frequently sampled. These results are consistent with the view that Ca²⁺ binding shifts the CBD12 conformational ensemble toward extended conformers, which could be a key step in the NCXs' allosteric regulation mechanism. This strategy, combining MD with NMR and SAXS, provides a powerful approach to select ensembles of conformations that could be applied to other flexible multidomain systems.

A Single Point Mutation Controls the Rate of Interconversion Between the g + and g - Rotamers of the Histidine 189 χ2 Angle That Activates Bacterial Enzyme I for Catalysis

Enzyme I (EI) of the bacterial phosphotransferase system (PTS) is a master regulator of bacterial metabolism and a promising target for development of a new class of broad-spectrum antibiotics. The catalytic activity of EI is mediated by several intradomain, interdomain, and intersubunit conformational equilibria. Therefore, in addition to its relevance as a drug target, EI is also a good model for investigating the dynamics/function relationship in multidomain, oligomeric proteins. Here, we use solution NMR and protein design to investigate how the conformational dynamics occurring within the N-terminal domain (EIN) affect the activity of EI. We show that the rotameric g⁺-to-g⁻ transition of the active site residue His¹⁸⁹ χ2 angle is decoupled from the state A-to-state B transition that describes a ∼90° rigid-body rearrangement of the EIN subdomains upon transition of the full-length enzyme to its catalytically competent closed form. In addition, we engineered EIN constructs with modulated conformational dynamics by hybridizing EIN from mesophilic and thermophilic species, and used these chimeras to assess the effect of increased or decreased active site flexibility on the enzymatic activity of EI. Our results indicate that the rate of the autophosphorylation reaction catalyzed by EI is independent from the kinetics of the g⁺-to-g⁻ rotameric transition that exposes the phosphorylation site on EIN to the incoming phosphoryl group. In addition, our work provides an example of how engineering of hybrid mesophilic/thermophilic chimeras can assist investigations of the dynamics/function relationship in proteins, therefore opening new possibilities in biophysics.

Structure elucidation of the elusive Enzyme I monomer reveals the molecular mechanisms linking oligomerization and enzymatic activity

Enzyme I (EI) is a phosphotransferase enzyme responsible for converting phosphoenolpyruvate (PEP) into pyruvate. This reaction initiates a five-step phosphorylation cascade in the bacterial phosphotransferase (PTS) transduction pathway. Under physiological conditions, EI exists in an equilibrium between a functional dimer and an inactive monomer. The monomer-dimer equilibrium is a crucial factor regulating EI activity and the phosphorylation state of the overall PTS. Experimental studies of EI's monomeric state have yet been hampered by the dimer's high thermodynamic stability, which prevents its characterization by standard structural techniques. In this study, we modified the dimerization domain of EI (EIC) by mutating three amino acids involved in the formation of intersubunit salt bridges. The engineered variant forms an active dimer in solution that can bind and hydrolyze PEP. Using hydrostatic pressure as an additional perturbation, we were then able to study the complete dissociation of the variant from 1 bar to 2.5 kbar in the absence and the presence of EI natural ligands. Backbone residual dipolar couplings collected under high-pressure conditions allowed us to determine the conformational ensemble of the isolated EIC monomeric state in solution. Our calculations reveal that three catalytic loops near the dimerization interface become unstructured upon monomerization, preventing the monomeric enzyme from binding its natural substrate. This study provides an atomic-level characterization of EI's monomeric state and highlights the role of the catalytic loops as allosteric connectors controlling both the activity and oligomerization of the enzyme.

Weak binding to the A2RE RNA rigidifies hnRNPA2 RRMs and reduces liquid-liquid phase separation and aggregation

hnRNPA2 is a major component of mRNA transport granules in oligodendrocytes and neurons. However, the structural details of how hnRNPA2 binds the A2 recognition element (A2RE) and if this sequence stimulates granule formation by enhancing phase separation of hnRNPA2 has not yet been studied. Using solution NMR and biophysical studies, we find that each of the two individual RRMs retain the domain structure observed in complex with RNA but are not rigidly confined (i.e. they move independently) in solution in the absence of RNA. hnRNPA2 RRMs bind the minimal rA2RE11 weakly but at least, and most likely, two hnRNPA2 molecules are able to simultaneously bind the longer 21mer myelin basic protein A2RE. Upon binding of the RNA, NMR chemical shift deviations are observed in both RRMs, suggesting both play a role in binding the A2RE11. Interestingly, addition of short A2RE RNAs or longer RNAs containing this sequence completely prevents in vitro phase separation of full-length hnRNPA2 and aggregation of the disease-associated mutants. These findings suggest that RRM interactions with specific recognition sequences alone do not account for nucleating granule formation, consistent with models where multivalent protein:RNA and protein:protein contacts form across many sites in granule proteins and long RNA transcripts.

A Time-Dependent Quantum Approach to Allostery and a Comparison With Light-Harvesting in Photosynthetic Phenomenon

The allosteric effect is one of the most important processes in regulating the function of proteins, and the elucidation of this phenomenon plays a significant role in understanding emergent behaviors in biological regulation. In this process, a perturbation, generated by a ligand in a part of the macromolecule (the allosteric site), moves along this system and reaches a specific (active) site, dozens of Ångströms away, with a great efficiency. The dynamics of this perturbation in the macromolecule can model precisely the allosteric process. In this article, we will be studying the general characteristics of allostery, using a time-dependent quantum approach to obtain rules that apply to this kind of process. Considering the perturbation as a wave that moves within the molecular system, we will characterize the allosteric process with three of the properties of this wave in the active site: (1) t_a, the characteristic time for reaching that site, (2) A_a, the amplitude of the wave in this site, and (3) B_a, its corresponding spectral broadening. These three parameters, together with the process mechanism and the perturbation efficiency in the process, can describe the phenomenon. One of the main purposes of this paper is to link the parameters t_a, A_a, and B_a and the perturbation efficiency to the characteristics of the system. There is another fundamental process for life that has some characteristics similar to allostery: the light-harvesting (LH) process in photosynthesis. Here, as in allostery, two distant macromolecular sites are involved—two sites dozens of Ångströms away. In both processes, it is particularly important that the perturbation is distributed efficiently without dissipating in the infinite degrees of freedom within the macromolecule. The importance of considering quantum effects in the LH process is well documented in literature, and the quantum coherences are experimentally proven by time-dependent spectroscopic techniques. Given the existing similarities between these two processes in macromolecules, in this work, we suggest using Quantum Mechanics (QM) to study allostery.

NMR Methods for Structural Characterization of Protein-Protein Complexes

Protein-protein interactions and the complexes thus formed are critical elements in a wide variety of cellular events that require an atomic-level description to understand them in detail. Such complexes typically constitute challenging systems to characterize and drive the development of innovative biophysical methods. NMR spectroscopy techniques can be applied to extract atomic resolution information on the binding interfaces, intermolecular affinity, and binding-induced conformational changes in protein-protein complexes formed in solution, in the cell membrane, and in large macromolecular assemblies. Here we discuss experimental techniques for the characterization of protein-protein complexes in both solution NMR and solid-state NMR spectroscopy. The approaches include solvent paramagnetic relaxation enhancement and chemical shift perturbations (CSPs) for the identification of binding interfaces, and the application of intermolecular nuclear Overhauser effect spectroscopy and residual dipolar couplings to obtain structural constraints of protein-protein complexes in solution. Complementary methods in solid-state NMR are described, with emphasis on the versatility provided by heteronuclear dipolar recoupling to extract intermolecular constraints in differentially labeled protein complexes. The methods described are of particular relevance to the analysis of membrane proteins, such as those involved in signal transduction pathways, since they can potentially be characterized by both solution and solid-state NMR techniques, and thus outline key developments in this frontier of structural biology.

Identifying dynamic, partially occupied residues using anomalous scattering

Structural studies of partially occupied, heterogeneous protein systems using crystallography are difficult. Here, methods are presented for detecting these states in crystals.

Although often presented as taking single ‘snapshots’ of the conformation of a protein, X-ray crystallography provides an averaged structure over time and space within the crystal. The important but difficult task of characterizing structural ensembles in crystals is typically limited to small conformational changes, such as multiple side-chain conformations. A crystallographic method was recently introduced that utilizes residual electron and anomalous density (READ) to characterize structural ensembles encompassing large-scale structural changes. Key to this method is an ability to accurately measure anomalous signals and distinguish them from noise or other anomalous scatterers. This report presents an optimized data-collection and analysis strategy for partially occupied iodine anomalous signals. Using the long-wavelength-optimized beamline I23 at Diamond Light Source, the ability to accurately distinguish the positions of anomalous scatterers with occupancies as low as ∼12% is demonstrated. The number and positions of these anomalous scatterers are consistent with previous biophysical, kinetic and structural data that suggest that the protein Im7 binds to the chaperone Spy in multiple partially occupied conformations. Finally, READ selections demonstrate that re-measured data using the new protocols are consistent with the previously characterized structural ensemble of the chaperone Spy with its client Im7. This study shows that a long-wavelength beamline results in easily validated anomalous signals that are strong enough to be used to detect and characterize highly disordered sections of crystal structures.

Tuning a timing device that regulates lateral root development in rice

Peptidyl Prolyl Isomerases (PPIases) accelerate cis-trans isomerization of prolyl peptide bonds. In rice, the PPIase LRT2 is essential for lateral root initiation. LRT2 displays in vitro isomerization of a highly conserved W-P peptide bond (¹⁰⁴W-P¹⁰⁵) in the natural substrate OsIAA11. OsIAA11 is a transcription repressor that, in response to the plant hormone auxin, is targeted to ubiquitin-mediated proteasomal degradation via specific recognition of the cis isomer of its ¹⁰⁴W-P¹⁰⁵ peptide bond. OsIAA11 controls transcription of specific genes, including its own, that are required for lateral root development. This auxin-responsive negative feedback circuit governs patterning and development of lateral roots along the primary root. The ability to tune LRT2 activity via mutagenesis is crucial for understanding and modeling the role of this bimodal switch in the auxin circuit and lateral root development. We present characterization of the thermal stability and isomerization rates of several LRT2 mutants acting on the OsIAA11 substrate. The thermally stable mutants display activities lower than that of wild-type (WT) LRT2. These include binding diminished but catalytically active P125K, binding incompetent W128A, and binding capable but catalytically incompetent H133Q mutations. Additionally, LRT2 homologs hCypA from human, TaCypA from Triticum aestivum (wheat) and PPIB from E. coli were shown to have 110, 50 and 60% of WT LRT2 activity on the OsIAA11 substrate. These studies identify several thermally stable LRT2 mutants with altered activities that will be useful for establishing relationships between cis-trans isomerization, auxin circuit dynamics, and lateral root development in rice.

View from Nuclear Magnetic Resonance Spectroscopy

Nuclear magnetic resonance (NMR) spectroscopy is one of the three major approaches for determining the structures of biological macromolecules. Historically, NMR spectroscopy was number two after X-ray crystallography in the rate of depositions to the Protein Data Bank (PDB). However, electron cryomicroscopy (CryoEM) recently surpassed NMR in this regard. NMR frequently is used in conjunction with X-ray or CryoEM in structure determinations. NMR has advantages over the other structural approaches in studies of conformational dynamics and interconverting conformational states of proteins and nucleic acids in solution. NMR spectroscopy, itself, can be considered as collection of hybrid methods in that structure determinations rely on the results of several separate magnetic resonance experiments that measure connectivities of magnetic-resonance-active nuclei through covalent bonds or through space or determine relative orientations of magnetic dipoles. NMR results frequently are combined with data from small-angle X-ray scattering or chemical crosslinking in developing structural models. NMR spectroscopy and CryoEM are particularly synergistic in that neither requires crystallization.

Applications of in Silico Methods for Design and Development of Drugs Targeting Protein-Protein Interactions

BACKGROUND

Identification of Protein-Protein Interactions (PPIs) is a major challenge in modern molecular biology and biochemistry research, due to the unquestionable role of proteins in cells, biological process and pathological states. Over the past decade, the PPIs have evolved from being considered a highly challenging field of research to being investigated and examined as targets for pharmacological intervention.

OBJECTIVE

Comprehension of protein interactions is crucial to known how proteins come together to build signalling pathways, to carry out their functions, or to cause diseases, when deregulated. Multiplicity and great amount of PPIs structures offer a huge number of new and potential targets for the treatment of different diseases.

METHODS

Computational techniques are becoming predominant in PPIs studies for their effectiveness, flexibility, accuracy and cost. As a matter of fact, there are effective in silico approaches which are able to identify PPIs and PPI site. Such methods for computational target prediction have been developed through molecular descriptors and data-mining procedures.

RESULTS

In this review, we present different types of interactions between protein-protein and the application of in silico methods for design and development of drugs targeting PPIs. We described computational approaches for the identification of possible targets on protein surface and to detect of stimulator/ inhibitor molecules.

CONCLUSION

A deeper study of the most recent bioinformatics methodologies for PPIs studies is vital for a better understanding of protein complexes and for discover new potential PPI modulators in therapeutic intervention.

Active Site Breathing of Human Alkbh5 Revealed by Solution NMR and Accelerated Molecular Dynamics

AlkB homolog 5 (Alkbh5) is one of nine members of the AlkB family, which are nonheme Fe²⁺/α-ketoglutarate-dependent dioxygenases that catalyze the oxidative demethylation of modified nucleotides and amino acids. Alkbh5 is highly selective for the N⁶-methyladenosine modification, an epigenetic mark that has spawned significant biological and pharmacological interest because of its involvement in important physiological processes, such as carcinogenesis and stem cell differentiation. Herein, we investigate the structure and dynamics of human Alkbh5 in solution. By using ¹⁵N and ¹³C_methyl relaxation dispersion and ¹⁵N-R₁ and R_1ρ NMR experiments, we show that the active site of apo Alkbh5 experiences conformational dynamics on multiple timescales. Consistent with this observation, backbone amide residual dipolar couplings measured for Alkbh5 in phage pf1 are inconsistent with the static crystal structure of the enzyme. We developed a simple approach that combines residual dipolar coupling data and accelerated molecular dynamics simulations to calculate a conformational ensemble of Alkbh5 that is fully consistent with the experimental NMR data. Our structural model reveals that Alkbh5 is more disordered in solution than what is observed in the crystal state and undergoes breathing motions that expand the active site and allow access to α-ketoglutarate. Disordered-to-ordered conformational changes induced by sequential substrate/cofactor binding events have been often invoked to interpret biochemical data on the activity and specificity of AlkB proteins. The structural ensemble reported in this work provides the first atomic-resolution model of an AlkB protein in its disordered conformational state to our knowledge.

Intermolecular hydrogen bonding H···Cl in crystal structure of palladium(II)-bis(diaminocarbene) complex

The reaction of bis(isocyanide)palladium complex cis-[PdCl2(CNXyl)2] (Xyl=2,6-Me2C6H3) with excess of 4,5-dichlorobenzene-1,2-amine in a C2H4Cl2/MeOH mixture affords monocationic bis(diaminocarbene) complex cis-[PdClC{(NHXyl)=NHC6H2Cl2 NH2}{C(NHXyl)=NHC6H2Cl2NH2}]Cl (3) in moderate yield (42%). Complex 3 exists in the solid phase in the H-bonded dimeric associate of two single charged organometallic cations and two chloride anions according to X-ray diffraction data. The Hirshfeld surface analysis for the X-ray structure of 3 reveals that the crystal packing is determined primarily by intermolecular contacts H–Cl, H–H, and H–C. The intermolecular hydrogen bonds N–H···Cl and C–H···Cl in the H-bonded dimeric associate of 3 were studied by DFT calculations and topological analysis of the electron density distribution within the framework of QTAIM method, and estimated energies of these supramolecular contacts vary from 1.6 to 9.1 kcal/mol. Such non-covalent bonding means that complex 3 is an anionic receptor for the chloride anions.

Probing the Atomic Structure of Transient Protein Contacts by Paramagnetic Relaxation Enhancement Solution NMR

Important biological processes, including enzyme catalysis, signaling, and protein folding, proceed through lowly populated (<5%) states that elude structural characterization by conventional techniques. Here, we describe the steps required for visualization of these sparsely populated conformations and transient protein-protein interactions using paramagnetic relaxation enhancement solution NMR. We describe experimental design, sample preparation, data acquisition and processing, and the basics of data analysis of structural ensembles.

Metadynamic metainference: Convergence towards force field independent structural ensembles of a disordered peptide

Metadynamic metainference has been recently introduced as a theoretical framework to determine structural ensembles by combining and weighting their noise multiple sources of experimental data with molecular mechanics force fields and metadynamics simulations. Here we build upon these initial developments to further extend and streamline the computational approach. We also show that metadynamic metainference can actually determine a structural ensemble for a disordered peptide that is essentially independent from the employed force field. We further show that it is possible to use a very computationally efficient implicit solvent force field in the place of very expensive state-of-the-art explicit solvent ones without a significant loss in accuracy.

Inferring properties of disordered chains from FRET transfer efficiencies

Förster resonance energy transfer (FRET) is a powerful tool for elucidating both structural and dynamic properties of unfolded or disordered biomolecules, especially in single-molecule experiments. However, the key observables, namely, the mean transfer efficiency and fluorescence lifetimes of the donor and acceptor chromophores, are averaged over a broad distribution of donor-acceptor distances. The inferred average properties of the ensemble therefore depend on the form of the model distribution chosen to describe the distance, as has been widely recognized. In addition, while the distribution for one type of polymer model may be appropriate for a chain under a given set of physico-chemical conditions, it may not be suitable for the same chain in a different environment so that even an apparently consistent application of the same model over all conditions may distort the apparent changes in chain dimensions with variation of temperature or solution composition. Here, we present an alternative and straightforward approach to determining ensemble properties from FRET data, in which the polymer scaling exponent is allowed to vary with solution conditions. In its simplest form, it requires either the mean FRET efficiency or fluorescence lifetime information. In order to test the accuracy of the method, we have utilized both synthetic FRET data from implicit and explicit solvent simulations for 30 different protein sequences, and experimental single-molecule FRET data for an intrinsically disordered and a denatured protein. In all cases, we find that the inferred radii of gyration are within 10% of the true values, thus providing higher accuracy than simpler polymer models. In addition, the scaling exponents obtained by our procedure are in good agreement with those determined directly from the molecular ensemble. Our approach can in principle be generalized to treating other ensemble-averaged functions of intramolecular distances from experimental data.

Hybrid Methods for Modeling Protein Structures Using Molecular Dynamics Simulations and Small-Angle X-Ray Scattering Data

Small-angle X-ray scattering (SAXS) is an efficient experimental tool to measure the overall shape of macromolecular structures in solution. However, due to the low resolution of SAXS data, high-resolution data obtained from X-ray crystallography or NMR and computational methods such as molecular dynamics (MD) simulations are complementary to SAXS data for understanding protein functions based on their structures at atomic resolution. Because MD simulations provide a physicochemically proper structural ensemble for flexible proteins in solution and a precise description of solvent effects, the hybrid analysis of SAXS and MD simulations is a promising method to estimate reasonable solution structures and structural ensembles in solution. Here, we review typical and useful in silico methods for modeling three dimensional protein structures, calculating theoretical SAXS profiles, and analyzing ensemble structures consistent with experimental SAXS profiles. We also review two examples of the hybrid analysis, termed MD-SAXS method in which MD simulations are carried out without any knowledge of experimental SAXS data, and the experimental SAXS data are used only to assess the consistency of the solution model from MD simulations with those observed in experiments. One example is an investigation of the intrinsic dynamics of EcoO109I using the computational method to obtain a theoretical profile from the trajectory of an MD simulation. The other example is a structural investigation of the vitamin D receptor ligand-binding domain using snapshots generated by MD simulations and assessment of the snapshots by experimental SAXS data.

Selecting Conformational Ensembles Using Residual Electron and Anomalous Density (READ)

Heterogeneous and dynamic biomolecular complexes play a central role in many cellular processes but are poorly understood due to experimental challenges in characterizing their structural ensembles. To address these difficulties, we developed a hybrid methodology that combines X-ray crystallography with ensemble selections typically used in NMR studies to determine structural ensembles of heterogeneous biomolecular complexes. The method, termed READ, for residual electron and anomalous density, enables the visualization of heterogeneous conformational ensembles of complexes within crystals. Here we present a detailed protocol for performing the ensemble selections to construct READ ensembles. From a diverse pool of binding poses, a selection scheme is used to determine a subset of conformations that maximizes agreement with the X-ray data. Overall, READ is a general approach for obtaining a high-resolution view of dynamic protein-protein complexes.

Dynamic Descriptions of Highly Flexible Molecules from NMR Dipolar Couplings: Physical Basis and Limitations

Biomolecules that control physiological function by changing their conformation play key roles in biology and remain poorly characterized. NMR dipolar couplings (DCs) depend intrinsically on both molecular shape and structural fluctuations, thereby providing the enticing prospect of tracking these conformational changes at atomic detail. Although this dual dependence has until now severely complicated analysis of DCs from highly dynamic systems, general approaches have recently been proposed that simplify interpretation of experimental DCs, by entirely eliminating molecular alignment from the analysis. Using simple and intuitive simulation of target ensembles, we investigate the impact of such approaches on the resulting descriptions of the conformational energy landscape. We find that ensemble descriptions of highly flexible systems derived from DCs without explicit consideration of the alignment properties of the constituent conformations can be compromised and inaccurate, despite exhibiting high correlation with experimental measurement.

Hybrid Methods in Iron-Sulfur Cluster Biogenesis

Hybrid methods, which combine and integrate several biochemical and biophysical techniques, have rapidly caught up in the last twenty years to provide a way to obtain a fuller description of proteins and molecular complexes with sizes and complexity otherwise not easily affordable. Here, we review the use of a robust hybrid methodology based on a mixture of NMR, SAXS, site directed mutagenesis and molecular docking which we have developed to determine the structure of weakly interacting molecular complexes. We applied this technique to gain insights into the structure of complexes formed amongst proteins involved in the molecular machine, which produces the essential iron-sulfur cluster prosthetic groups. Our results were validated both by X-ray structures and by other groups who adopted the same approach. We discuss the advantages and the limitations of our methodology and propose new avenues, which could improve it.

Protein Data Bank (PDB): The Single Global Macromolecular Structure Archive

The Protein Data Bank (PDB)--the single global repository of experimentally determined 3D structures of biological macromolecules and their complexes--was established in 1971, becoming the first open-access digital resource in the biological sciences. The PDB archive currently houses ~130,000 entries (May 2017). It is managed by the Worldwide Protein Data Bank organization (wwPDB; wwpdb.org), which includes the RCSB Protein Data Bank (RCSB PDB; rcsb.org), the Protein Data Bank Japan (PDBj; pdbj.org), the Protein Data Bank in Europe (PDBe; pdbe.org), and BioMagResBank (BMRB; www.bmrb.wisc.edu). The four wwPDB partners operate a unified global software system that enforces community-agreed data standards and supports data Deposition, Biocuration, and Validation of ~11,000 new PDB entries annually (deposit.wwpdb.org). The RCSB PDB currently acts as the archive keeper, ensuring disaster recovery of PDB data and coordinating weekly updates. wwPDB partners disseminate the same archival data from multiple FTP sites, while operating complementary websites that provide their own views of PDB data with selected value-added information and links to related data resources. At present, the PDB archives experimental data, associated metadata, and 3D-atomic level structural models derived from three well-established methods: crystallography, nuclear magnetic resonance spectroscopy (NMR), and electron microscopy (3DEM). wwPDB partners are working closely with experts in related experimental areas (small-angle scattering, chemical cross-linking/mass spectrometry, Forster energy resonance transfer or FRET, etc.) to establish a federation of data resources that will support sustainable archiving and validation of 3D structural models and experimental data derived from integrative or hybrid methods.

Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.