Atomistic ensemble modeling and small-angle neutron scattering of intrinsically disordered protein complexes: applied to minichromosome maintenance protein

SAXS/MC studies of the mixed-folded protein Cdt1 reveal monomeric, folded over conformations

Cdt1 is a protein critical for DNA replication licensing and is well-established to be a binding partner of the minichromosome maintenance (MCM) complex. Cdt1 has also been demonstrated to have an emerging, "moonlighting" role at the kinetochore via direct binding to microtubules and to the Ndc80 complex. However, it is not known how the structure and conformations of Cdt1 could allow for these multiple, completely unique sets of protein complexes. And while there exist multiple robust methods to study entirely folded or entirely unfolded proteins, structure-function studies of combined, mixed folded/disordered proteins remain challenging. It this work, we employ multiple orthogonal biophysical and computational techniques to provide a detailed structural characterization of human Cdt1 92-546. DSF and DSCD show both folded winged helix (WH) domains of Cdt1 are relatively unstable. CD and NMR show the N-terminal and the linker regions are intrinsically disordered. Using DLS and SEC-MALS, we show that Cdt1 is polydisperse, monomeric at high concentrations, and without any apparent inter-molecular self-association. SEC-SAXS of the monomer in solution enabled computational modeling of the protein in silico. Using the program SASSIE, we performed rigid body Monte Carlo simulations to generate a conformational ensemble. Using experimental SAXS data, we filtered for conformations which did and did not fit our data. We observe that neither fully extended nor extremely compact Cdt1 conformations are consistent with our SAXS data. The best fit models have the N-terminal and linker regions extended into solution and the two folded domains close to each other in apparent "folded over" conformations. The best fit Cdt1 conformations are consistent with a function as a scaffold protein which may be sterically blocked without the presence of binding partners. Our studies also provide a template for combining experimental and computational biophysical techniques to study mixed-folded proteins.

Targeting Intrinsically Disordered Proteins through Dynamic Interactions

Intrinsically disordered proteins (IDPs) are over-represented in major disease pathways and have attracted significant interest in understanding if and how they may be targeted using small molecules for therapeutic purposes. While most existing studies have focused on extending the traditional structure-centric drug design strategies and emphasized exploring pre-existing structure features of IDPs for specific binding, several examples have also emerged to suggest that small molecules could achieve specificity in binding IDPs and affect their function through dynamic and transient interactions. These dynamic interactions can modulate the disordered conformational ensemble and often lead to modest compaction to shield functionally important interaction sites. Much work remains to be done on further elucidation of the molecular basis of the dynamic small molecule-IDP interaction and determining how it can be exploited for targeting IDPs in practice. These efforts will rely critically on an integrated experimental and computational framework for disordered protein ensemble characterization. In particular, exciting advances have been made in recent years in enhanced sampling techniques, Graphic Processing Unit (GPU)-computing, and protein force field optimization, which have now allowed rigorous physics-based atomistic simulations to generate reliable structure ensembles for nontrivial IDPs of modest sizes. Such de novo atomistic simulations will play crucial roles in exploring the exciting opportunity of targeting IDPs through dynamic interactions.

Unwinding 20 Years of the Archaeal Minichromosome Maintenance Helicase

Replicative DNA helicases are essential cellular enzymes that unwind duplex DNA in front of the replication fork during chromosomal DNA replication. Replicative helicases were discovered, beginning in the 1970s, in bacteria, bacteriophages, viruses, and eukarya, and, in the mid-1990s, in archaea. This year marks the 20th anniversary of the first report on the archaeal replicative helicase, the minichromosome maintenance (MCM) protein. This minireview summarizes 2 decades of work on the archaeal MCM.

Monte Carlo simulation algorithm for B-DNA

Understanding the structure-function relationship of biomolecules containing DNA has motivated experiments aimed at determining molecular structure using methods such as small-angle X-ray and neutron scattering (SAXS and SANS). SAXS and SANS are useful for determining macromolecular shape in solution, a process which benefits by using atomistic models that reproduce the scattering data. The variety of algorithms available for creating and modifying model DNA structures lack the ability to rapidly modify all-atom models to generate structure ensembles. This article describes a Monte Carlo algorithm for simulating DNA, not with the goal of predicting an equilibrium structure, but rather to generate an ensemble of plausible structures which can be filtered using experimental results to identify a sub-ensemble of conformations that reproduce the solution scattering of DNA macromolecules. The algorithm generates an ensemble of atomic structures through an iterative cycle in which B-DNA is represented using a wormlike bead-rod model, new configurations are generated by sampling bend and twist moves, then atomic detail is recovered by back mapping from the final coarse-grained configuration. Using this algorithm on commodity computing hardware, one can rapidly generate an ensemble of atomic level models, each model representing a physically realistic configuration that could be further studied using molecular dynamics. © 2016 Wiley Periodicals, Inc.

Self-assembly and soluble aggregate behavior of computationally designed coiled-coil peptide bundles

Coiled-coil peptides have proven useful in a range of materials applications ranging from the formation of well-defined fibrils to responsive hydrogels. The ability to design from first principles their oligomerization and subsequent higher order assembly offers their expanded use in producing new materials. Toward these ends, homo-tetrameric, antiparallel, coiled-coil, peptide bundles have been designed computationally, synthesized via solid-phase methods, and their solution behavior characterized. Two different bundle-forming peptides were designed and examined. Within the targeted coiled coil structure, both bundles contained the same hydrophobic core residues. However, different exterior residues on the two different designs yielded sequences with different distributions of charged residues and two different expected isoelectric points of pI 4.4 and pI 10.5. Both coiled-coil bundles were extremely stable with respect to temperature (Tm > 80 C) and remained soluble in solution even at high (millimolar) peptide concentrations. The coiled-coil tetramer was confirmed to be the dominant species in solution by analytical sedimentation studies and by small-angle neutron scattering, where the scattering form factor is well represented by a cylinder model with the dimensions of the targeted coiled coil. At high concentrations (5-15 mM), evidence of interbundle structure was observed via neutron scattering. At these concentrations, the synthetic bundles form soluble aggregates, and interbundle distances can be determined via a structure factor fit to scattering data. The data support the successful design of robust coiled-coil bundles. Despite their different sequences, each sequence forms loosely associated but soluble aggregates of the bundles, suggesting similar dissociated states for each. The behavior of the dispersed bundles is similar to that observed for natural proteins.

Preparing monodisperse macromolecular samples for successful biological small-angle X-ray and neutron-scattering experiments

Small-angle X-ray scattering (SAXS) and small-angle neutron scattering (SANS) are techniques used to extract structural parameters and determine the overall structures and shapes of biological macromolecules, complexes and assemblies in solution. The scattering intensities measured from a sample contain contributions from all atoms within the illuminated sample volume, including the solvent and buffer components, as well as the macromolecules of interest. To obtain structural information, it is essential to prepare an exactly matched solvent blank so that background scattering contributions can be accurately subtracted from the sample scattering to obtain the net scattering from the macromolecules in the sample. In addition, sample heterogeneity caused by contaminants, aggregates, mismatched solvents, radiation damage or other factors can severely influence and complicate data analysis, so it is essential that the samples be pure and monodisperse for the duration of the experiment. This protocol outlines the basic physics of SAXS and SANS, and it reveals how the underlying conceptual principles of the techniques ultimately 'translate' into practical laboratory guidance for the production of samples of sufficiently high quality for scattering experiments. The procedure describes how to prepare and characterize protein and nucleic acid samples for both SAXS and SANS using gel electrophoresis, size-exclusion chromatography (SEC) and light scattering. Also included are procedures that are specific to X-rays (in-line SEC-SAXS) and neutrons, specifically preparing samples for contrast matching or variation experiments and deuterium labeling of proteins.

Archaeal MCM Proteins as an Analog for the Eukaryotic Mcm2-7 Helicase to Reveal Essential Features of Structure and Function

In eukaryotes, the replicative helicase is the large multisubunit CMG complex consisting of the Mcm2–7 hexameric ring, Cdc45, and the tetrameric GINS complex. The Mcm2–7 ring assembles from six different, related proteins and forms the core of this complex. In archaea, a homologous MCM hexameric ring functions as the replicative helicase at the replication fork. Archaeal MCM proteins form thermostable homohexamers, facilitating their use as models of the eukaryotic Mcm2–7 helicase. Here we review archaeal MCM helicase structure and function and how the archaeal findings relate to the eukaryotic Mcm2–7 ring.

Deuterium Labeling Together with Contrast Variation Small-Angle Neutron Scattering Suggests How Skp Captures and Releases Unfolded Outer Membrane Proteins

In Gram-negative bacteria, the chaperone protein Skp forms specific and stable complexes with membrane proteins while they are transported across the periplasm to the outer membrane. The jellyfish-like architecture of Skp is similar to the eukaryotic and archaeal prefoldins and the mitochondrial Tim chaperones, that is the α-helical "tentacles" extend from a β-strand "body" to create an internal cavity. Contrast variation small-angle neutron scattering (SANS) experiments on Skp alone in solution and bound in two different complexes to unfolded outer membrane proteins (uOMPs), OmpA and OmpW, demonstrate that the helical tentacles of Skp bind their substrate in a clamp-like mechanism in a conformation similar to that previously observed in the apo crystal structure of Skp. Deuteration of the uOMP component combined with contrast variation analysis allowed the shapes of Skp and uOMP as well as the location of uOMP with respect to Skp to be determined in both complexes. This represents unique information that could not be obtained without deuterium labeling of the uOMPs. The data yield the first direct structural evidence that the α-helical Skp tentacles move closer together on binding its substrate and that the structure of Skp is different when binding different uOMPs. This work presents, by example, a tutorial on performing SANS experiments using both deuterium labeling and contrast variation, including SANS theory, sample preparation, data collection, sample quality validation, data analysis, and structure modeling.

Archaeal DNA replication

DNA replication is essential for all life forms. Although the process is fundamentally conserved in the three domains of life, bioinformatic, biochemical, structural, and genetic studies have demonstrated that the process and the proteins involved in archaeal DNA replication are more similar to those in eukaryal DNA replication than in bacterial DNA replication, but have some archaeal-specific features. The archaeal replication system, however, is not monolithic, and there are some differences in the replication process between different species. In this review, the current knowledge of the mechanisms governing DNA replication in Archaea is summarized. The general features of the replication process as well as some of the differences are discussed.

On the energy components governing molecular recognition in the framework of continuum approaches

Molecular recognition is a process that brings together several biological macromolecules to form a complex and one of the most important characteristics of the process is the binding free energy. Various approaches exist to model the binding free energy, provided the knowledge of the 3D structures of bound and unbound molecules. Among them, continuum approaches are quite appealing due to their computational efficiency while at the same time providing predictions with reasonable accuracy. Here we review recent developments in the field emphasizing on the importance of adopting adequate description of physical processes taking place upon the binding. In particular, we focus on the efforts aiming at capturing some of the atomistic details of the binding phenomena into the continuum framework. When possible, the energy components are reviewed independently of each other. However, it is pointed out that rigorous approaches should consider all energy contributions on the same footage. The two major schemes for utilizing the individual energy components to predict binding affinity are outlined as well.

Rapid and accurate calculation of small-angle scattering profiles using the golden ratio

Calculating the scattering intensity of anN-atom system is a numerically exhaustingO(N2) task. A simple approximation technique that scales linearly with the number of atoms is presented. Using an exact expression for the scattering intensityI(q) at a given wavevectorq, the rotationally averaged intensityI(q) is computed by evaluatingI(q) in several scattering directions. The orientations of theqvectors are taken from a quasi-uniform spherical grid generated by the golden ratio. Using various biomolecules as examples, this technique is compared with an established multipole expansion method. For a given level of speed, the technique is more accurate than the multipole expansion for anisotropically shaped molecules, while comparable in accuracy for globular shapes. The processing time scales sub-linearly inNwhen the atoms are identical and lie on a lattice. The procedure is easily implemented and should accelerate the analysis of small-angle scattering data.