A New Method for Determining Structure Ensemble: Application to a RNA Binding Di-Domain Protein

Conformational Dynamics of the Periplasmic Chaperone SurA

The periplasmic protein SurA is the primary chaperone involved in the biogenesis of bacterial outer membrane proteins and is a potential antibacterial drug target. The three-dimensional structure of SurA can be divided into three parts, a core module formed by the N- and C-terminal regions and two peptidyl-prolyl isomerase (PPIase) domains, P1 and P2. Despite the determination of the structures of several SurA-peptide complexes, the functional mechanism of this chaperone remains elusive and the roles of the two PPIase domains are yet unclear. Herein, we characterize the conformational dynamics of SurA by using solution nuclear magnetic resonance and single-molecule fluorescence resonance energy transfer methods. We demonstrate a "closed-to-open" structural transition of the P1 domain that is correlated with both chaperone activity and peptide binding and show that the flexible P2 domain can also occupy conformations that closely contact the NC core module. Our results offer a structural basis for the counteracting roles of the two PPIase domains in regulating the SurA chaperone activity.

Conformational states of Zika virus non-structural protein 3 determined by molecular dynamics simulations with small-angle X-Ray scattering data

Zika virus (ZIKV) has become a great public health emergency. Its non-structural protein 3 (NS3) is a key enzyme in viral replication and has been considered as a potential therapeutic target. A conformational characterization of ZIKV NS3 is critical for a comprehensive understanding of its molecular interactions and functions. However, the high conformational flexibility of solution NS3 obstacles the structural characterization of NS3 solely from the experimental observable that averages over its heterogeneous conformations. Here, we employed replica exchange with solute tempering (REST) method to simulate the di-domain protein ZIKV NS3. Three independent MD simulations identified a conserved conformational ensemble of NS3, consisting of a major conformational state and several minor states from compact to loose conformations. The major state agrees well with the scattering profile from small-angle X-ray scattering (SAXS) experiments. Moreover, the simulated ensemble is supported by a direct data-fitting result that requires both short- and long-range structural contacts to recover the experimental data. We discussed the interplay between simulation and experiment in ensemble construction of flexible biomolecules and shed light on the physically derived conformational ensembles.

A Method for Determining Structure Ensemble of Large Disordered Protein: Application to a Mechanosensing Protein

Structure characterization of intrinsically disordered proteins (IDPs) remains a key obstacle in understanding their functional mechanisms. Due to the highly dynamic feature of IDPs, structure ensembles instead of static unique structures are often derived from experimental data. Several state-of-the-art computational methods have been developed to select an optimal ensemble from a pregenerated structure pool, but they suffer from low efficiency for large IDPs. Here we present a matching pursuit genetic algorithm (MPGA) for structure ensemble determination, which takes advantages from both matching pursuit (MP) to reduce the search space and genetic algorithm (GA) to reduce the restriction on constraint types. The MPGA method is validated using a reference ensemble with predefined structures. In comparison with the conventional GA, MPGA takes much less computational time for large IDPs. The utility of the method is demonstrated by application to structure ensemble determination of a mechanosensing protein domain with 306 amino acids. The structure ensemble determined reveals that the N-terminal region 1-240 is more compact than the C-terminal region 240-306. The unique structural feature explains why only a small portion of YXXP tyrosine residues can be phosphorylated easily by kinases in the absence of extension force and why the phosphorylation is force-dependent.

Investigating energy-based pool structure selection in the structure ensemble modeling with experimental distance constraints: The example from a multidomain protein Pub1

The structural variations of multidomain proteins with flexible parts mediate many biological processes, and a structure ensemble can be determined by selecting a weighted combination of representative structures from a simulated structure pool, producing the best fit to experimental constraints such as interatomic distance. In this study, a hybrid structure-based and physics-based atomistic force field with an efficient sampling strategy is adopted to simulate a model di-domain protein against experimental paramagnetic relaxation enhancement (PRE) data that correspond to distance constraints. The molecular dynamics simulations produce a wide range of conformations depicted on a protein energy landscape. Subsequently, a conformational ensemble recovered with low-energy structures and the minimum-size restraint is identified in good agreement with experimental PRE rates, and the result is also supported by chemical shift perturbations and small-angle X-ray scattering data. It is illustrated that the regularizations of energy and ensemble-size prevent an arbitrary interpretation of protein conformations. Moreover, energy is found to serve as a critical control to refine the structure pool and prevent data overfitting, because the absence of energy regularization exposes ensemble construction to the noise from high-energy structures and causes a more ambiguous representation of protein conformations. Finally, we perform structure-ensemble optimizations with a topology-based structure pool, to enhance the understanding on the ensemble results from different sources of pool candidates.

Visualizing single-stranded nucleic acids in solution

Single-stranded nucleic acids (ssNAs) are ubiquitous in many key cellular functions. Their flexibility limits both the number of high-resolution structures available, leaving only a small number of protein-ssNA crystal structures, while forcing solution investigations to report ensemble averages. A description of the conformational distributions of ssNAs is essential to more fully characterize biologically relevant interactions. We combine small angle X-ray scattering (SAXS) with ensemble-optimization methods (EOM) to dynamically build and refine sets of ssNA structures. By constructing candidate chains in representative dinucleotide steps and refining the models against SAXS data, a broad array of structures can be obtained to match varying solution conditions and strand sequences. In addition to the distribution of large scale structural parameters, this approach reveals, for the first time, intricate details of the phosphate backbone and underlying strand conformations. Such information on unperturbed strands will critically inform a detailed understanding of an array of problems including protein-ssNA binding, RNA folding and the polymer nature of NAs. In addition, this scheme, which couples EOM selection with an iteratively refining pool to give confidence in the underlying structures, is likely extendable to the study of other flexible systems.

Coexistence of multiple minor states of fatty acid binding protein and their functional relevance

Proteins are dynamic over a wide range of timescales, but determining the number of distinct dynamic processes and identifying functionally relevant dynamics are still challenging. Here we present the study on human intestinal fatty acid binding protein (hIFABP) using a novel analysis of ¹⁵N relaxation dispersion (RD) and chemical shift saturation transfer (CEST) experiments. Through combined analysis of the two types of experiments, we found that hIFABP exists in a four-state equilibrium in which three minor states interconvert directly with the major state. According to conversion rates from the major “closed” state to minor states, these minor states are irrelevant to the function of fatty acid transport. Based on chemical shifts of the minor states which could not be determined from RD data alone but were extracted from a combined analysis of RD and CEST data, we found that all the minor states are native-like. This conclusion is further supported by hydrogen-deuterium exchange experiments. Direct conversions between the native state and native-like intermediate states may suggest parallel multitrack unfolding/folding pathways of hIFABP. Moreover, hydrogen-deuterium exchange data indicate the existence of another locally unfolded minor state that is relevant to the fatty acid entry process.