Christian A. Hanke

Quantitative FRET studies and integrative modeling unravel the structure and dynamics of biomolecular systems

Förster Resonance Energy Transfer (FRET) combined with single-molecule spectroscopy probes macromolecular structure and dynamics and identifies coexisting conformational states. We review recent methodological developments in integrative structural modeling by satisfying spatial restraints on networks of FRET pairs (hybrid-FRET). We discuss procedures to incorporate prior structural knowledge and to obtain optimal distance networks. Finally, a workflow for hybrid-FRET is presented that automates integrative structural modeling and experiment planning to put hybrid-FRET on rails. To test this workflow, we simulate realistic single-molecule experiments and resolve three protein conformers, exchanging at 30μs and 10ms, with accuracies of 1-3Å RMSD versus the target structure. Incorporation of data from other spectroscopies and imaging is also discussed.

Rigidity theory for biomolecules: concepts, software, and applications

The mechanical heterogeneity of biomolecular structures is intimately linked to their diverse biological functions. Applying rigidity theory to biomolecules identifies this heterogeneous composition of flexible and rigid regions, which can aid in the understanding of biomolecular stability and long‐ranged information transfer through biomolecules, and yield valuable information for rational drug design and protein engineering. We review fundamental concepts in rigidity theory, ways to represent biomolecules as constraint networks, and methodological and algorithmic developments for analyzing such networks and linking the results to biomolecular function. Software packages for performing rigidity analyses on biomolecules in an efficient, automated way are described, as are rigidity analyses on biomolecules including the ribosome, viruses, or transmembrane proteins. The analyses address questions of allosteric mechanisms, mutation effects on (thermo‐)stability, protein (un‐)folding, and coarse‐graining of biomolecules. We advocate that the application of rigidity theory to biomolecules has matured in such a way that it could be broadly applied as a computational biophysical method to scrutinize biomolecular function from a structure‐based point of view and to complement approaches focused on biomolecular dynamics. We discuss possibilities to improve constraint network representations and to perform large‐scale and prospective studies. WIREs Comput Mol Sci 2017, 7:e1311. doi: 10.1002/wcms.1311

Precision and accuracy of single-molecule FRET measurements—a multi-laboratory benchmark study

Single-molecule Förster resonance energy transfer (smFRET) is increasingly being used to determine distances, structures, and dynamics of biomolecules in vitro and in vivo. However, generalized protocols and FRET standards to ensure the reproducibility and accuracy of measurements of FRET efficiencies are currently lacking. Here we report the results of a comparative blind study in which 20 labs determined the FRET efficiencies (E) of several dye-labeled DNA duplexes. Using a unified, straightforward method, we obtained FRET efficiencies with s.d. between ±0.02 and ±0.05. We suggest experimental and computational procedures for converting FRET efficiencies into accurate distances, and discuss potential uncertainties in the experiment and the modeling. Our quantitative assessment of the reproducibility of intensity-based smFRET measurements and a unified correction procedure represents an important step toward the validation of distance networks, with the ultimate aim of achieving reliable structural models of biomolecular systems by smFRET-based hybrid methods.A multi-laboratory study finds that single-molecule FRET is a reproducible and reliable approach for determining accurate distances in dye-labeled DNA duplexes.

Force Field Dependence of Riboswitch Dynamics

Riboswitches are noncoding regulatory elements that control gene expression in response to the presence of metabolites, which bind to the aptamer domain. Metabolite binding appears to occur through a combination of conformational selection and induced fit mechanism. This demands to characterize the structural dynamics of the apo state of aptamer domains. In principle, molecular dynamics (MD) simulations can give insights at the atomistic level into the dynamics of the aptamer domain. However, it is unclear to what extent contemporary force fields can bias such insights. Here, we show that the Amber force field ff99 yields the best agreement with detailed experimental observations on differences in the structural dynamics of wild type and mutant aptamer domains of the guanine-sensing riboswitch (Gsw), including a pronounced influence of Mg2+. In contrast, applying ff99 with parmbsc0 and parmχOL modifications (denoted ff10) results in strongly damped motions and overly stable tertiary loop-loop interactions. These results are based on 58 MD simulations with an aggregate simulation time>11 μs, careful modeling of Mg2+ ions, and thorough statistical testing. Our results suggest that the moderate stabilization of the χ-anti region in ff10 can have an unwanted damping effect on functionally relevant structural dynamics of marginally stable RNA systems. This suggestion is supported by crystal structure analyses of Gsw aptamer domains that reveal χ torsions with high-anti values in the most mobile regions. We expect that future RNA force field development will benefit from considering marginally stable RNA systems and optimization toward good representations of dynamics in addition to structural characteristics.

Ligand-mediated and tertiary interactions cooperatively stabilize the P1 region in the guanine-sensing riboswitch

Riboswitches are genetic regulatory elements that control gene expression depending on ligand binding. The guanine-sensing riboswitch (Gsw) binds ligands at a three-way junction formed by paired regions P1, P2, and P3. Loops L2 and L3 cap the P2 and P3 helices and form tertiary interactions. Part of P1 belongs to the switching sequence dictating the fate of the mRNA. Previous studies revealed an intricate relationship between ligand binding and presence of the tertiary interactions, and between ligand binding and influence on the P1 region. However, no information is available on the interplay among these three main regions in Gsw. Here we show that stabilization of the L2-L3 region by tertiary interactions, and the ligand binding site by ligand binding, cooperatively influences the structural stability of terminal base pairs in the P1 region in the presence of Mg2+ ions. The results are based on molecular dynamics simulations with an aggregate simulation time of ~10 μs across multiple systems of the unbound state of the Gsw aptamer and a G37A/C61U mutant, and rigidity analyses. The results could explain why the three-way junction is a central structural element also in other riboswitches and how the cooperative effect could become contextual with respect to intracellular Mg2+ concentration. The results suggest that the transmission of allosteric information to P1 can be entropy-dominated.

Publisher Correction: Precision and accuracy of single-molecule FRET measurements—a multi-laboratory benchmark study

This paper was originally published under standard Springer Nature copyright. As of the date of this correction, the Analysis is available online as an open-access paper with a CC-BY license. No other part of the paper has been changed.

Tertiary interactions in the unbound guanine-sensing riboswitch focus functional conformational variability on the binding site

Riboswitches are genetic regulatory elements mainly found in bacteria, which regulate gene expression based on the availability of a ligand. Purine-sensing riboswitches, including the guanine-sensing riboswitch (Gsw), possess tertiary interactions connecting the L2 and L3 loops. These interactions are important for ligand binding to the aptamer. However, atomic-level structural knowledge about the unbound state and how the tertiary interactions influence the conformational heterogeneity of the aptamer is still scarce. We performed replica exchange molecular dynamics simulations of the aptamer domain of wild-type Gsw and a G37A/C61U mutant, which exhibits destabilized tertiary interactions, at different Mg2+ concentrations with an aggregate simulation time of ∼16 μs, and subsequently obtained free-energy landscapes. Our data provide evidence that suggests that the unbound state of wild-type Gsw is conformationally rather homogeneous from a global viewpoint, yet the ligand binding site shows functionally necessary mobility required for ligand binding. For the mutant, the data suggest a heterogeneous ensemble, in particular without Mg2+. Hence, the tertiary interactions focus functional conformational variability on the binding site region of wild-type Gsw. Our data allow speculating that already the weakening of the tertiary interactions by two hydrogen bonds shifts the kinetics of folding from downhill folding without traps or intermediate states for wild-type Gsw to a folding including intermediates and misfolded structures for the mutant. A slowed-down folding of the aptamer might favor a decision during transcriptional regulation for the off-path, even if the ligand binds.