Measurement of Mesoscale Conformational Dynamics of Freely Diffusing Molecules with Tracking FCS

Few techniques are suited to probe the structure and dynamics of molecular complexes at the mesoscale level ($～$100–1000 nm). We have developed a single-molecule technique that uses tracking fluorescence correlation spectroscopy (tFCS) to probe the conformation and dynamics of mesoscale molecular assemblies. tFCS measures the distance fluctuations between two fluorescently labeled sites within an untethered, freely diffusing biomolecule. To achieve subdiffraction spatial resolution, we developed a feedback scheme that allows us to maintain the molecule at an optimal position within the laser intensity gradient for fluorescence correlation spectroscopy. We characterized tFCS spatial sensitivity by measuring the Brownian end-to-end dynamics of DNA molecules as short as 1000 bp. We demonstrate that tFCS detects changes in the compaction of reconstituted nucleosome arrays and can assay transient protein-mediated interactions between distant sites in an individual DNA molecule. Our measurements highlight the applicability of tFCS to a wide variety of biochemical processes involving mesoscale conformational dynamics.     

Lipid Fingerprints and Cofactor Dynamics of Light-Harvesting Complex II in Different Membranes

Plant light-harvesting complex II (LHCII) is the key antenna complex for plant photosynthesis. We present coarse-grained molecular dynamics simulations of monomeric and trimeric LHCII in a realistic thylakoid membrane environment based on the Martini force field. The coarse-grained protein model has been optimized with respect to atomistic reference simulations. Our simulations provide detailed insights in the thylakoid lipid fingerprint of LHCII which compares well with experimental data from membrane protein purification. Comparing the monomer and trimeric LHCII reveals a stabilizing effect of trimerization on the chromophores as well as the protein. Moreover, the average chromophore distance shortens in the trimer leading to stronger excitonic couplings. When changing the native thylakoid environment to a model membrane the protein flexibility remains constant, whereas the chromophore flexibility is reduced. Overall, the presented LHCII model lays the foundation to investigate the μs dynamics of this key antenna protein of plants.     

Intrinsically Disordered Protein Ntr2 Modulates the Spliceosomal RNA Helicase Brr2

Precursor messenger RNA splicing is mediated by the spliceosome, a large and dynamic molecular machine composed of five small nuclear RNAs and numerous proteins. Many spliceosomal proteins are predicted to be intrinsically disordered or contain large disordered regions, but experimental validation of these predictions is scarce, and the precise functions of these proteins are often unclear. Here, we show via circular dichroism spectroscopy, dynamic light scattering, and NMR spectroscopy that the yeast spliceosomal disassembly factor Ntr2 is largely intrinsically disordered. Peptide SPOT analyses, analytical size-exclusion chromatography, and surface plasmon resonance measurements revealed that Ntr2 uses an N-terminal region to bind the C-terminal helicase unit of the Brr2 RNA helicase, an enzyme involved in spliceosome activation and implicated in splicing catalysis and spliceosome disassembly. NMR analyses suggested that Ntr2 does not adopt a tertiary structure and likely remains disordered upon complex formation. RNA binding and unwinding studies showed that Ntr2 downregulates Brr2 helicase activity in vitro by modulating the fraction of helicase molecules productively bound to the RNA substrate. Our data clarify the nature of a physical link between Brr2 and Ntr2, and point to the possibility of a functional Ntr2-Brr2 interplay during splicing.     

Simultaneous Determination of Protein Structure and Dynamics Using Cryo-Electron Microscopy

Cryo-electron microscopy is rapidly emerging as a powerful technique to determine the structures of complex macromolecular systems elusive to other techniques. Because many of these systems are highly dynamical, characterizing their movements is also a crucial step to unravel their biological functions. To achieve this goal, we report an integrative modeling approach to simultaneously determine structure and dynamics of macromolecular systems from cryo-electron microscopy density maps. By quantifying the level of noise in the data and dealing with their ensemble-averaged nature, this approach enables the integration of multiple sources of information to model ensembles of structures and infer their populations. We illustrate the method by characterizing structure and dynamics of the integral membrane receptor STRA6, thus providing insights into the mechanisms by which it interacts with retinol binding protein and translocates retinol across the membrane.     

Structure of a Therapeutic Full-Length Anti-NPRA IgG4 Antibody: Dissecting Conformational Diversity

We report the x-ray crystal structure of intact, full-length human immunoglobulin (IgG4) at 1.8 Å resolution. The data for IgG4 (S228P), an antibody targeting the natriuretic peptide receptor A, show a previously unrecognized type of Fab-Fc orientation with a distorted λ-shape in which one Fab-arm is oriented toward the Fc portion. Detailed structural analysis by x-ray crystallography and molecular simulations suggest that this is one of several conformations coexisting in a dynamic equilibrium state. These results were confirmed by small angle x-ray scattering in solution. Furthermore, electron microscopy supported these findings by preserving molecule classes of different conformations. This study fosters our understanding of IgG4 in particular and our appreciation of antibody flexibility in general. Moreover, we give insights into potential biological implications, specifically for the interaction of human anti-natriuretic peptide receptor A IgG4 with the neonatal Fc receptor, Fcγ receptors, and complement-activating C1q by considering conformational flexibility.     

Neurons with Complex Karyotypes Are Rare in Aged Human Neocortex

1. Download high-res image (220KB)
2. Download full-size image

     

Integrating NMR,SAXS, and Atomistic Simulations: Structure and Dynamics of a Two-Domain Protein

Multidomain proteins with two or more independently folded functional domains are prevalent in nature. Whereas most multidomain proteins are linked linearly in sequence, roughly one-tenth possess domain insertions where a guest domain is implanted into a loop of a host domain, such that the two domains are connected by a pair of interdomain linkers. Here, we characterized the influence of the interdomain linkers on the structure and dynamics of a domain-insertion protein in which the guest LysM domain is inserted into a central loop of the host CVNH domain. Expanding upon our previous crystallographic and NMR studies, we applied SAXS in combination with NMR paramagnetic relaxation enhancement to construct a structural model of the overall two-domain system. Although the two domains have no fixed relative orientation, certain orientations were found to be preferred over others. We also assessed the accuracies of molecular mechanics force fields in modeling the structure and dynamics of tethered multidomain proteins by integrating our experimental results with microsecond-scale atomistic molecular dynamics simulations. In particular, our evaluation of two different combinations of the latest force fields and water models revealed that both combinations accurately reproduce certain structural and dynamical properties, but are inaccurate for others. Overall, our study illustrates the value of integrating experimental NMR and SAXS studies with long timescale atomistic simulations for characterizing structural ensembles of flexibly linked multidomain systems.     

Structure of the Human Core Centromeric Nucleosome Complex

Download : Download video (54MB)     

ESR Resolves the C Terminus Structure of the Ligand-free Human Glutathione S-Transferase A1-1

Nitroxide- and Cu²⁺-based electron spin resonance (ESR) are combined to provide insight into the conformational states of the functionally important α-helix of the human glutathione S-transferase A1. Distance measurements on various spin-labeled dimeric human glutathione S-transferase A1-1 all result in bimodal distance distributions, indicating that the C-terminus exists in two distinct conformations in solution, one of which closely matches that found in the crystal structure of the ligand-bound enzyme. These measurements permit the generation of a model of the unliganded conformation. Room temperature ESR indicates that the second conformation has high mobility, potentially enabling the enzyme’s high degree of substrate promiscuity. This model is then validated using computational modeling and further Cu²⁺-based ESR distance measurements. Cu²⁺-based ESR also provides evidence that the secondary structure of the second conformation is of helical nature. Addition of S-hexyl glutathione results in a shift in relative populations, favoring the state that is similar to the previously known structure of the ligand-bound enzyme.     

Topological Constraints and Their Conformational Entropic Penalties on RNA Folds

Functional RNAs can fold into intricate structures using a number of different secondary and tertiary structural motifs. Many factors contribute to the overall free energy of the target fold. This study aims at quantifying the entropic costs coming from the loss of conformational freedom when the sugar-phosphate backbone is subjected to constraints imposed by secondary and tertiary contacts. Motivated by insights from topology theory, we design a diagrammatic scheme to represent different types of RNA structures so that constraints associated with a folded structure may be segregated into mutually independent subsets, enabling the total conformational entropy loss to be easily calculated as a sum of independent terms. We used high-throughput Monte Carlo simulations to simulate large ensembles of single-stranded RNA sequences in solution to validate the assumptions behind our diagrammatic scheme, examining the entropic costs for hairpin initiation and formation of many multiway junctions. Our diagrammatic scheme aids in the factorization of secondary/tertiary constraints into distinct topological classes and facilitates the discovery of interrelationships among multiple constraints on RNA folds. This perspective, which to our knowledge is novel, leads to useful insights into the inner workings of some functional RNA sequences, demonstrating how they might operate by transforming their structures among different topological classes.     

Rotamer Dynamics: Analysis of Rotamers in Molecular Dynamics Simulations of Proteins

Given by χ torsional angles, rotamers describe the side-chain conformations of amino acid residues in a protein based on the rotational isomers (hence the word rotamer). Constructed rotamer libraries, based on either protein crystal structures or dynamics studies, are the tools for classifying rotamers (torsional angles) in a way that reflect their frequency in nature. Rotamer libraries are routinely used in structure modeling and evaluation. In this perspective article, we would like to encourage researchers to apply rotamer analyses beyond their traditional use. Molecular dynamics (MD) of proteins highlight the in silico behavior of molecules in solution and thus can identify favorable side-chain conformations. In this article, we used simple computational tools to study rotamer dynamics (RD) in MD simulations. First, we isolated each frame in the MD trajectories in separate Protein Data Bank files via the cpptraj module in AMBER. Then, we extracted torsional angles via the Bio3D module in R language. The classification of torsional angles was also done in R according to the penultimate rotamer library. RD analysis is useful for various applications such as protein folding, study of rotamer-rotamer relationship in protein-protein interaction, real-time correlation between secondary structures and rotamers, study of flexibility of side chains in binding site for molecular docking preparations, use of RD as guide in functional analysis and study of structural changes caused by mutations, providing parameters for improving coarse-grained MD accuracy and speed, and many others. Major challenges facing RD to emerge as a new scientific field involve the validation of results via easy, inexpensive wet-lab methods. This realm is yet to be explored.     

Probing Conformational Changes during the Gating Cycle of a Potassium Channel in Lipid Bilayers

Ion conduction across the cellular membrane requires the simultaneous opening of activation and inactivation gates of the K⁺ channel pore. The bacterial KcsA channel has served as a powerful system for dissecting the structural changes that are related to four major functional states associated with K⁺ gating. Yet, the direct observation of the full gating cycle of KcsA has remained structurally elusive, and crystal structures mimicking these gating events require mutations in or stabilization of functionally relevant channel segments. Here, we found that changes in lipid composition strongly increased the KcsA open probability. This enabled us to probe all four major gating states in native-like membranes by combining electrophysiological and solid-state NMR experiments. In contrast to previous crystallographic views, we found that the selectivity filter and turret region, coupled to the surrounding bilayer, were actively involved in channel gating. The increase in overall steady-state open probability was accompanied by a reduction in activation-gate opening, underscoring the important role of the surrounding lipid bilayer in the delicate conformational coupling of the inactivation and activation gates.     

Thermodynamics of Conformational Transitions in a Disordered Protein Backbone Model

Conformational entropy is expected to contribute significantly to the thermodynamics of structural transitions in intrinsically disordered proteins or regions in response to protein/ligand binding, posttranslational modifications, and environmental changes. We calculated the backbone (dihedral) conformational entropy of oligoglycine $\left({\text{Gly}}_N\right)$, a protein backbone mimic and model intrinsically disordered region, as a function of chain length ($N=3$, 4, 5, 10, and 15) from simulations using three different approaches. The backbone conformational entropy scales linearly with chain length with a slope consistent with the entropy of folding of well-structured proteins. The entropic contributions of second-order dihedral correlations are predominantly through intraresidue ?-ψ pairs, suggesting that oligoglycine may be thermodynamically modeled as a system of independent glycine residues. We find the backbone conformational entropy to be largely independent of global structural parameters, like the end-to-end distance and radius of gyration. We introduce a framework referred to herein as “ensemble confinement” to estimate the loss (gain) of conformational free energy and its entropic component when individual residues are constrained to (released from) particular regions of the ?-ψ map. Quantitatively, we show that our protein backbone model resists ordering/folding with a significant, unfavorable ensemble confinement free energy because of the loss of a substantial portion of the absolute backbone entropy. Proteins can couple this free-energy reservoir to distal binding events as a regulatory mechanism to promote or suppress binding.     

Molecular Dynamics of the Association of L-Selectin and FERM Regulated by PIP2

Phosphatidylinositol 4,5-bisphosphate (PIP2) acts as a signaling lipid, mediating membrane trafficking and recruitment of proteins to membranes. A key example is the PIP2-dependent regulation of the adhesion of L-selectin to the cytoskeleton adaptors of the N-terminal subdomain of ezrin-radixin-moesin (FERM). The molecular details of the mediating behavior of multivalent anionic PIP2 lipids in this process, however, remain unclear. Here, we use coarse-grained molecular dynamics simulation to explore the mechanistic details of PIP2 in the transformation, translocation, and association of the FERM/L-selectin complex. We compare membranes of different compositions and find that anionic phospholipids are necessary for both FERM and the cytoplasmic domain of L-selectin to absorb on the membrane surface. The subsequent formation of the FERM/L-selectin complex is strongly favored by the presence of PIP2, which clusters around both proteins and triggers a conformational transition in the cytoplasmic domain of L-selectin. We are able to quantify the effect of PIP2 on the association free energy of the complex by means of a potential of mean force. We conclude that PIP2 behaves as an adhesive agent to enhance the stability of the FERM/L-selectin complex and identify key residues involved. The molecular information revealed in this study highlights the specific role of membrane lipids such as PIP2 in protein translocation and potential signaling.     

Architectural Dynamics of CaMKII-Actin Networks

Calcium-calmodulin-dependent kinase II (CaMKII) has an important role in dendritic spine remodeling upon synaptic stimulation. Using fluorescence video microscopy and image analysis, we investigated the architectural dynamics of rhodamine-phalloidin stabilized filamentous actin (F-actin) networks cross-linked by CaMKII. We used automated image analysis to identify F-actin bundles and crossover junctions and developed a dimensionless metric to characterize network architecture. Similar networks were formed by three different CaMKII species with a 10-fold length difference in the linker region between the kinase domain and holoenzyme hub, implying linker length is not a primary determinant of F-actin cross-linking. Electron micrographs showed that at physiological molar ratios, single CaMKII holoenzymes cross-linked multiple F-actin filaments at random, whereas at higher CaMKII/F-actin ratios, filaments bundled. Light microscopy established that the random network architecture resisted macromolecular crowding with polyethylene glycol and blocked ATP-powered compaction by myosin-II miniature filaments. Importantly, the networks disassembled after the addition of calcium-calmodulin and were then spaced within 3 min into compacted foci by myosin motors or more slowly (30 min) aggregated by crowding. Single-molecule total internal reflection fluorescence microscopy showed CaMKII dissociation from surface-immobilized globular actin exhibited a monoexponential dwell-time distribution, whereas CaMKII bound to F-actin networks had a long-lived fraction, trapped at crossover junctions. Release of CaMKII from F-actin, triggered by calcium-calmodulin, was too rapid to measure with flow-cell exchange (<20 s). The residual bound fraction was reduced substantially upon addition of an N-methyl-D-aspartate receptor peptide analog but not ATP. These results provide mechanistic insights to CaMKII-actin interactions at the collective network and single-molecule level. Our findings argue that CaMKII-actin networks in dendritic spines maintain spine size against physical stress. Upon synaptic stimulation, CaMKII is disengaged by calcium-calmodulin, triggering network disassembly, expansion, and subsequent compaction by myosin motors with kinetics compatible with the times recorded for the poststimulus changes in spine volume.