Dihedral Restraints Research Articles

Restraint validation of biomolecular structures determined by NMR in the Protein Data Bank

Biomolecular structure analysis from experimental NMR studies generally relies on restraints derived from a combination of experimental and knowledge-based data. A challenge for the structural biology community has been a lack of standards for representing these restraints, preventing the establishment of uniform methods of model-vs-data structure validation against restraints and limiting interoperability between restraint-based structure modeling programs. The NEF and NMR-STAR formats provide a standardized approach for representing commonly used NMR restraints. Using these restraint formats, a standardized validation system for assessing structural models of biopolymers against restraints has been developed and implemented in the wwPDB OneDep data deposition-validation-biocuration system. The resulting wwPDB restraint violation report provides a model vs. data assessment of biomolecule structures determined using distance and dihedral restraints, with extensions to other restraint types currently being implemented. These tools are useful for assessing NMR models, as well as for assessing biomolecular structure predictions based on distance restraints.

Minimal Increments of Hydrophobic Collapse within the N-Terminus of the Neuropeptide Galanin.

The neuropeptide galanin has a 35-year history as an intriguing target in drug design owing to its implication as a potential anticonvulsant and neuronal trophic factor among many other therapeutically interesting functions including analgesia and mood alteration. In this study, we report the structural characterization of three synthetic fragments of the galanin N-terminus in buffered aqueous solution: hGal(2-12)KK, hGal(1-12)KK, and hGal(1-17)KK. High-field two-dimensional 1H-1H nuclear magnetic resonance (NMR) data were acquired for these fragments and used to derive distance restraints. We further utilized modified hydrogen bonding and dihedral restraints to reflect chemical shift patterns in the data, which revealed the signature of a weakly folded helix. Together, these sets of restraints were used to generate NMR structures of all three fragments, which depict a core of hydrophobic residues that cluster together regardless of the presence of a helical structure, and correspond to residues in the N-terminus of galanin that have been previously shown to be critical for binding its receptors. The helical structure only appears following the inclusion of Gly(1) in the sequence, and at longer sequence lengths, unlike many other peptides, the helix does not propagate. Rather, a few turns of poorly ordered helix appear to be a secondary consequence of clusters of hydrophobic sidechains that are conserved across all of the peptides in this study; the helices themselves appear ordered as a consequence of this clustering, and these clusters compare directly to those observed recently to make contacts between galanin and two of its receptor subtypes. Collapsed hydrophobic residues therefore organize and compose the functional core of human galanin and raise interesting questions about the nature of the conformational order in ligands that bind cell surface receptors.

Concurrent Identification and Characterization of Protein Structure and Continuous Internal Dynamics with REDCRAFT.

Internal dynamics of proteins can play a critical role in the biological function of some proteins. Several well documented instances have been reported such as MBP, DHFR, hTS, DGCR8, and NSP1 of the SARS-CoV family of viruses. Despite the importance of internal dynamics of proteins, there currently are very few approaches that allow for meaningful separation of internal dynamics from structural aspects using experimental data. Here we present a computational approach named REDCRAFT that allows for concurrent characterization of protein structure and dynamics. Here, we have subjected DHFR (PDB-ID 1RX2), a 159-residue protein, to a fictitious, mixed mode model of internal dynamics. In this simulation, DHFR was segmented into 7 regions where 4 of the fragments were fixed with respect to each other, two regions underwent rigid-body dynamics, and one region experienced uncorrelated and melting event. The two dynamical and rigid-body segments experienced an average orientational modification of 7° and 12° respectively. Observable RDC data for backbone C′-N, N-HN, and C′-HN were generated from 102 uniformly sampled frames that described the molecular trajectory. The structure calculation of DHFR with REDCRAFT by using traditional Ramachandran restraint produced a structure with 29 Å of structural difference measured over the backbone atoms (bb-rmsd) over the entire length of the protein and an average bb-rmsd of more than 4.7 Å over each of the dynamical fragments. The same exercise repeated with context-specific dihedral restraints generated by PDBMine produced a structure with bb-rmsd of 21 Å over the entire length of the protein but with bb-rmsd of less than 3 Å over each of the fragments. Finally, utilization of the Dynamic Profile generated by REDCRAFT allowed for the identification of different dynamical regions of the protein and the recovery of individual fragments with bb-rmsd of less than 1 Å. Following the recovery of the fragments, our assembly procedure of domains (larger segments consisting of multiple fragments with a common dynamical profile) correctly assembled the four fragments that are rigid with respect to each other, categorized the two domains that underwent rigid-body dynamics, and identified one dynamical region for which no conserved structure could be defined. In conclusion, our approach was successful in identifying the dynamical domains, recovery of structure where it is meaningful, and relative assembly of the domains when possible.

NMR refinement and peptide folding using the GROMACS software

Nuclear magnetic resonance spectroscopy is used routinely for studying the three-dimensional structures and dynamics of proteins and nucleic acids. Structure determination is usually done by adding restraints based upon NMR data to a classical energy function and performing restrained molecular simulations. Here we report on the implementation of a script to extract NMR restraints from a NMR-STAR file and export it to the GROMACS software. With this package it is possible to model distance restraints, dihedral restraints and orientation restraints. The output from the script is validated by performing simulations with and without restraints, including the ab initio refinement of one peptide.

A three-dimensional potential of mean force to improve backbone and sidechain hydrogen bond geometry in Xplor-NIH protein structure determination.

We introduce a new hydrogen bonding potential of mean force generated from high-quality crystal structures for use in Xplor-NIH structure calculations. This term applies to hydrogen bonds involving both backbone and sidechain atoms. When used in structure refinement calculations of 10 example protein systems with experimental distance, dihedral and residual dipolar coupling restraints, we demonstrate that the new term has superior performance to the previously developed hydrogen bonding potential of mean force used in Xplor-NIH.

Prediction of disulfide dihedral angles using chemical shifts.

Cystine residues result from the formation of disulfide bonds between pairs of cysteine residues. This cross linking of the backbone is essential for the structure and activity of peptides and proteins. The conformation of a cystine side chain can be described using five dihedral angles, χ1, χ2, χ3, χ2', and χ1', with cystines favouring certain combinations of these angles. 2D NMR spectroscopy is ideally suited for structure determination of disulfide-rich peptides, because of their small size and constrained nature. However, only limited information of the cystine side chain conformation can be determined by NMR spectroscopy, leading to ambiguity in the deduced 3D structures. Resolving accurate structures is important as disulfide-rich peptides have proven to be promising drug candidates in a number of fields, either as bioactive leads or scaffolds. Using a database of NMR chemical shifts combined with crystallographic structures, we have developed a method called DISH that uses support vector machines to predict the dihedral angles of cysteine side chains. It is able to successfully predict χ2 angles with 91% accuracy, and has improved performance over existing prediction methods for χ1 angles, with 87% accuracy. For 81% of cysteine residues, DISH successfully predicted both the χ1 and χ2 angles. By revisiting published solution structures of peptides determined using NMR spectroscopy, we assessed the impact of additional cystine dihedral restraints on the quality of 3D models. DISH improved the resolution and accuracy, highlighting the potential for improving the understanding of structure-activity relationships and rational development of peptide drugs.

An editor for the generation and customization of geometry restraints.

Chemical restraints for use in macromolecular structure refinement are produced by a variety of methods, including a number of programs that use chemical information to generate the required bond, angle, dihedral, chiral and planar restraints. These programs help to automate the process and therefore minimize the errors that could otherwise occur if it were performed manually. Furthermore, restraint-dictionary generation programs can incorporate chemical and other prior knowledge to provide reasonable choices of types and values. However, the use of restraints to define the geometry of a molecule is an approximation introduced with efficiency in mind. The representation of a bond as a parabolic function is a convenience and does not reflect the true variability in even the simplest of molecules. Another complicating factor is the interplay of the molecule with other parts of the macromolecular model. Finally, difficult situations arise from molecules with rare or unusual moieties that may not have their conformational space fully explored. These factors give rise to the need for an interactive editor for WYSIWYG interactions with the restraints and molecule. Restraints Editor, Especially Ligands (REEL) is a graphical user interface for simple and error-free editing along with additional features to provide greater control of the restraint dictionaries in macromolecular refinement.

Probing Selection Mechanism of the Most Favorable Conformation of a Dipeptide in Chaotropic and Kosmotropic Solution.

Chaotropes like urea and guanidinium chloride (GdmCl) tend to destabilize, and kosmotropes like proline tend to stabilize folded structures of peptides and proteins. Here, we combine fluorescence anisotropy decay measurements and molecular dynamics simulations to gain a microscopic understanding of the molecular mechanism for shifting conformational preferences in aqueous, GdmCl, urea, and proline solutions of a simple model dipeptide, N-acetyl-tryptophan-amide (NATA). Measured anisotropy decay of NATA as a function of temperature, pH, and cosolvent concentrations showed reorientations moderately slower in GdmCl and urea and substantially slower in proline compared to those of aqueous environment. A small change in pH significantly slows orientation time in water and GdmCl and less markedly in urea. Computationally, we use molecular dynamics with dihedral restraints to separately analyze the motions and interactions of the representative NATA conformers in the four different solvent environments. This novel analysis provides a dissection of the observed overall diffusion rates into contributions from individual dipeptide conformations. The variation of rotational diffusion rates with conformation are quite large. Population-weighted averaging or using properties of the major cluster reproduces the dynamical features of the full unrestrained dynamics. Additionally, we correlate the observable diffusion rates with microscopic features of conformer size, shape, and solvation. This analysis uncovered underlying differences in detailed atomistic behavior of the three cosolvents-urea, GdmCl, and proline. For both urea and the pure water system we find good agreement with hydrodynamic theory, with diffusion rates primarily correlated with conformer size and shape. In contrast, for GdmCl and proline solutions, the variation in conformer diffusion rates was mostly determined by specific interactions with the cosolvents. We also find preferences for different molecular shapes by the three cosolvents, with increased preferential solvation of smaller and more spherical conformers by urea and larger and more elongated conformers by GdmCl and proline. Additionally, our results provide a basis for a simple approximate model of the effects of pH lowering on dipeptide conformational equilibria. The translational diffusion rates of NATA are less sensitive to conformations, but variation with solvation strength is similar to rotational diffusion. Our results, combining experiment and simulation, show that we can identify the individual peptide conformers with definite microscopic properties of shape, size, and solvation, that are responsible for producing physical observables, such as translational and orientational diffusion in the complex solvent environments of denaturants and osmolytes.

Human furin Cys198 imposes dihedral and positional restraints on His194 for optimal Ser386-proton transfer

Inhibitors of human furin may represent the clinical remedy for very aggressive cancer, viral, and bacterial infections. Most of the currently available inhibitors are weak in terms of potency, drug-likeness, and furin specificity thereby necessitating the development of newer compounds especially mechanism-based inhibitors. Here, the roles of active site Cys198 (C198), His194 (H194), and Ser386 (S386) were investigated using computational-site-directed mutagenesis and molecular dynamics (MD) simulation. Data were obtained from six (6) biosystems: wildtype (C198/S386), furin-C198G (S386), furin-S386G (C198), and their peptide (nascent hydrolyzed peptide H2N-RTRR-CO2) bound complexes. The results strongly supported that in wildtype furin but not S386G and C198G mutants, following S386/scissile carbon attack (4.0 Å), the peptide retracted from the active site, representing peptide release post hydrolysis. Furthermore, in S386G mutant, C194 side chain thiolate ion may act as the nucleophile replacement but competing electron-rich centers (H194, H364) and energetically unattainable geometric strain on the peptide may constitute the limiting factors. In biosystems not complexed with peptide (representative of pre-attack state), C198 preferentially engaged H194 imidazole moiety via sulfur–π bond system causing a dihedral and positional restraints on the imidazole ring for ultimate alignment of its NE2 hydrogen atom with the side chain enolate oxygen of S364 required for optimal proton transfer. In summary, small-molecular-weight compounds with dual serine and cysteine protease inhibitory actions may represent a new class of potent and furin-selective compounds for future clinical applications.

Conformational Study on Pheromonotropin neuropeptide using NMR and Molecular Dynamics

Pheromonotropic neuropeptide, Pseudaletia pheromonotropin, is an 18 amino acid peptide with the sequence Lys-Leu-Ser-Tyr-Asp-Asp-Lys-Val-Phe-Glu-Asn-Val-Glu-Phe-Thr-Pro-Arg-Leu which is widely distributed in female moths of Bombyx mori. This peptide is structurally related to leucopyrokinin, an insect myotropic neuropeptide, responsible for melanization and reddish coloration hormone (MRCH) activity. A combined NMR and Molecular Dynamics (MD) methods were used to fathom the conformational behaviour in water which was used as a solvent. The structure was investigated using 1D-NMR and 2D-NMR (COSY, TOCSY and ROESY) experiments. The conformation was built by constrained MD simulations using distance and dihedral restraints from NMR data in the GROMACS simulation package. The peptide predominately adopts a β-sheet structure in water.

Physics Based Structure Refinement in Casp11 using Geometric Unfolding and Hierarchically Restrained Replica Exchange Molecular Dynamics

We have developed a new physics based approach to the protein refinement problem by mimicking the mechanism of chaperons that rehabilitate proteins. The template structure is unfolded by selectively pulling on different portions of the protein using the geometric based technique FRODA, and refolding the protein using hierarchically restrained replica exchange molecular dynamics. FRODA unfolding is used to create a diverse set of topologies for surveying near the native like structures from a template. The unfolding trajectories are then used to find energetic restraints to enforce contacts and dihedral restraints. An REMD simulation is performed for the entire ensemble using consensus and reservoir techniques, which allow multiple structural candidates to “swap” into the replica cascade at the highest temperature replica and the most favorable folds to propagate to the lowest temperature replica. The restraints are added in a hierarchical fashion where local contacts are restrained first followed by the addition of non-local restraints to narrow the conformational search toward the native state. After REMD structures are clustered, refined structures are selected based on the highest populated cluster, RMSD and DFIRE score.

Different inward and outward conduction mechanisms in NaVMs suggested by molecular dynamics simulations.

Rapid and selective ion transport is essential for the generation and regulation of electrical signaling pathways in living organisms. Here, we use molecular dynamics (MD) simulations with an applied membrane potential to investigate the ion flux of bacterial sodium channel NaVMs. 5.9 µs simulations with 500 mM NaCl suggest different mechanisms for inward and outward flux. The predicted inward conductance rate of ∼27±3 pS, agrees with experiment. The estimated outward conductance rate is 15±3 pS, which is considerably lower. Comparing inward and outward flux, the mean ion dwell time in the selectivity filter (SF) is prolonged from 13.5±0.6 ns to 20.1±1.1 ns. Analysis of the Na+ distribution revealed distinct patterns for influx and efflux events. In 32.0±5.9% of the simulation time, the E53 side chains adopted a flipped conformation during outward conduction, whereas this conformational change was rarely observed (2.7±0.5%) during influx. Further, simulations with dihedral restraints revealed that influx is less affected by the E53 conformational flexibility. In contrast, during outward conduction, our simulations indicate that the flipped E53 conformation provides direct coordination for Na+. The free energy profile (potential of mean force calculations) indicates that this conformational change lowers the putative barriers between sites SCEN and SHFS during outward conduction. We hypothesize that during an action potential, the increased Na+ outward transition propensities at depolarizing potentials might increase the probability of E53 conformational changes in the SF. Subsequently, this might be a first step towards initiating slow inactivation.

Structure and Activity of the Peptidyl-Prolyl Isomerase Domain from the Histone Chaperone Fpr4 toward Histone H3 Proline Isomerization

The FK506-binding protein (FKBP) family of peptidyl-prolyl isomerases (PPIases) is characterized by a common catalytic domain that binds to the inhibitors FK506 and rapamycin. As one of four FKBPs within the yeast Saccharomyces cerevisiae, Fpr4 has been described as a histone chaperone, and is in addition implicated in epigenetic function in part due to its mediation of cis-trans conversion of proline residues within histone tails. To better understand the molecular details of this activity, we have determined the solution structure of the Fpr4 C-terminal PPIase domain by using NMR spectroscopy. This canonical FKBP domain actively increases the rate of isomerization of three decapeptides derived from the N terminus of yeast histone H3, whereas maintaining intrinsic cis and trans populations. Observation of the uncatalyzed and Fpr4-catalyzed isomerization rates at equilibrium demonstrate Pro(16) and Pro(30) of histone H3 as the major proline targets of Fpr4, with little activity shown against Pro(38). This alternate ranking of the three target prolines, as compared with affinity determination or the classical chymotrypsin-based fluorescent assay, reveals the mechanistic importance of substrate residues C-terminal to the peptidyl-prolyl bond.

CING: an integrated residue-based structure validation program suite.

We present a suite of programs, named CING for Common Interface for NMR Structure Generation that provides for a residue-based, integrated validation of the structural NMR ensemble in conjunction with the experimental restraints and other input data. External validation programs and new internal validation routines compare the NMR-derived models with empirical data, measured chemical shifts, distance- and dihedral restraints and the results are visualized in a dynamic Web 2.0 report. A red–orange–green score is used for residues and restraints to direct the user to those critiques that warrant further investigation. Overall green scores below ~20 % accompanied by red scores over ~50 % are strongly indicative of poorly modelled structures. The publically accessible, secure iCing webserver (https://nmr.le.ac.uk) allows individual users to upload the NMR data and run a CING validation analysis.Electronic supplementary materialThe online version of this article (doi:10.1007/s10858-012-9669-7) contains supplementary material, which is available to authorized users.

TCRep 3D: An Automated In Silico Approach to Study the Structural Properties of TCR Repertoires

TCRep 3D is an automated systematic approach for TCR-peptide-MHC class I structure prediction, based on homology and ab initio modeling. It has been considerably generalized from former studies to be applicable to large repertoires of TCR. First, the location of the complementary determining regions of the target sequences are automatically identified by a sequence alignment strategy against a database of TCR Vα and Vβ chains. A structure-based alignment ensures automated identification of CDR3 loops. The CDR are then modeled in the environment of the complex, in an ab initio approach based on a simulated annealing protocol. During this step, dihedral restraints are applied to drive the CDR1 and CDR2 loops towards their canonical conformations, described by Al-Lazikani et. al. We developed a new automated algorithm that determines additional restraints to iteratively converge towards TCR conformations making frequent hydrogen bonds with the pMHC. We demonstrated that our approach outperforms popular scoring methods (Anolea, Dope and Modeller) in predicting relevant CDR conformations. Finally, this modeling approach has been successfully applied to experimentally determined sequences of TCR that recognize the NY-ESO-1 cancer testis antigen. This analysis revealed a mechanism of selection of TCR through the presence of a single conserved amino acid in all CDR3β sequences. The important structural modifications predicted in silico and the associated dramatic loss of experimental binding affinity upon mutation of this amino acid show the good correspondence between the predicted structures and their biological activities. To our knowledge, this is the first systematic approach that was developed for large TCR repertoire structural modeling.

Interaction with Polyglutamine-expanded Huntingtin Alters Cellular Distribution and RNA Processing of Huntingtin Yeast Two-hybrid Protein A (HYPA)

Huntington disease (HD) is an autosomal inherited disorder that causes the deterioration of brain cells. The polyglutamine (polyQ) expansion of huntingtin (Htt) is implicated in the pathogenesis of HD via interaction with an RNA splicing factor, Htt yeast two-hybrid protein A/forming-binding protein 11 (HYPA/FBP11). Besides the pathogenic polyQ expansion, Htt also contains a proline-rich region (PRR) located exactly in the C terminus to the polyQ tract. However, how the polyQ expansion influences the PRR-mediated protein interaction and how this abnormal interaction leads to the biological consequence remain elusive. Our NMR structural analysis indicates that the PRR motif of Htt cooperatively interacts with the tandem WW domains of HYPA through domain chaperoning effect of WW1 on WW2. The polyQ-expanded Htt sequesters HYPA to the cytosolic location and then significantly reduces the efficiency of pre-mRNA splicing. We propose that the toxic gain-of-function of the polyQ-expanded Htt that causes dysfunction of cellular RNA processing contributes to the pathogenesis of HD.

Solution Structure of the Guanine Nucleotide-binding STAS Domain of SLC26-related SulP Protein Rv1739c from Mycobacterium tuberculosis

The structure and intrinsic activities of conserved STAS domains of the ubiquitous SulP/SLC26 anion transporter superfamily have until recently remained unknown. Here we report the heteronuclear, multidimensional NMR spectroscopy solution structure of the STAS domain from the SulP/SLC26 putative anion transporter Rv1739c of Mycobacterium tuberculosis. The 0.87-Å root mean square deviation structure revealed a four-stranded β-sheet with five interspersed α-helices, resembling the anti-σ factor antagonist fold. Rv1739c STAS was shown to be a guanine nucleotide-binding protein, as revealed by nucleotide-dependent quench of intrinsic STAS fluorescence and photoaffinity labeling. NMR chemical shift perturbation analysis partnered with in silico docking calculations identified solvent-exposed STAS residues involved in nucleotide binding. Rv1739c STAS was not an in vitro substrate of mycobacterial kinases or anti-σ factors. These results demonstrate that Rv1739c STAS binds guanine nucleotides at physiological concentrations and undergoes a ligand-induced conformational change but, unlike anti-σ factor antagonists, may not mediate signals via phosphorylation.

High Resolution 1H-Detected Solid-State NMR Spectroscopy of Protein Aliphatic Resonances: Access to Tertiary Structure Information

Biological magic angle spinning (MAS) solid-state nuclear magnetic resonance spectroscopy has developed rapidly over the past two decades. For the structure determination of a protein by solid-state NMR, routinely (13)C,(13)C distance restraints as well as dihedral restraints are employed. In protonated samples, this is achieved by growing the bacterium on a medium which contains [1,3]-(13)C glycerol or [2]-(13)C glycerol to dilute the (13)C spin system. Labeling schemes, which rely on heteronuclei, are insensitive both for detection and in terms of quantification of distances, since they are relying on low-γ nuclei. Proton detection can in principle provide a gain in sensitivity by a factor of 8 and 31, compared to the (13)C or (15)N detected version of the experiment. We report here a new labeling scheme, which enables (1)H-detection of aliphatic resonances with high resolution in MAS solid-state NMR spectroscopy. We prepared microcrystals of the SH3 domain of chicken α-spectrin with 5% protonation at nonexchangeable sites and obtained line widths on the order of 25 Hz for aliphatic (1)H resonances. We show further that (13)C resolved 3D-(1)H,(1)H correlation experiments yield access to long-range proton-proton distances in the protein.

Simultaneous structure and dynamics of a membrane protein using REDCRAFT: Membrane-bound form of Pf1 coat protein

A strategy for simultaneous study of the structure and internal dynamics of a membrane protein is described using the REDCRAFT algorithm. The membrane-bound form of the Pf1 major coat protein (mbPf1) was used as an example. First, synthetic data is utilized to validate the simultaneous study of structure and dynamics with REDCRAFT using dihedral restraints and backbone N–H RDCs from two different alignments. Subsequently, the validated analysis is applied to experimental data and confirms that REDCRAFT produces meaningful structures from sparse RDC data. Furthermore, simulated data from a two-state jump motion is used to illustrate the necessity for simultaneous consideration of structure and dynamics. Disregarding internal dynamics during the course of structure determination is shown to produce an average-state that is not related to the two intermediate states. During the analysis of RDC data from the dynamic model, REDCRAFT appropriately identifies the region separating the static and dynamic domains of the protein. Finally, analysis of experimental data strongly suggests the existence of internal motion between the amphipathic and the transmembrane helices of the membrane-bound form of the protein. The ability to perform fragmented structure determination of each domain without a priori assumption of the order tensors allows an independent determination of the order tensors, which yields a more comprehensive description of protein structure and dynamics and is particularly relevant to the study of membrane proteins.

Structure of a DNA Polymerase α-Primase Domain That Docks on the SV40 Helicase and Activates the Viral Primosome

DNA polymerase alpha-primase (pol-prim) plays a central role in DNA replication in higher eukaryotes, initiating synthesis on both leading and lagging strand single-stranded DNA templates. Pol-prim consists of a primase heterodimer that synthesizes RNA primers, a DNA polymerase that extends them, and a fourth subunit, p68 (also termed B-subunit), that is thought to regulate the complex. Although significant knowledge about single-subunit primases of prokaryotes has accumulated, the functions and regulation of pol-prim remain poorly understood. In the SV40 replication model, the p68 subunit is required for primosome activity and binds directly to the hexameric viral helicase T antigen, suggesting a functional link between T antigen-p68 interaction and primosome activity. To explore this link, we first mapped the interacting regions of the two proteins and discovered a previously unrecognized N-terminal globular domain of p68 (p68N) that physically interacts with the T antigen helicase domain. NMR spectroscopy was used to determine the solution structure of p68N and map its interface with the T antigen helicase domain. Structure-guided mutagenesis of p68 residues in the interface diminished T antigen-p68 interaction, confirming the interaction site. SV40 primosome activity of corresponding pol-prim mutants decreased in proportion to the reduction in p68N-T antigen affinity, confirming that p68-T antigen interaction is vital for primosome function. A model is presented for how this interaction regulates SV40 primosome activity, and the implications of our findings are discussed in regard to the molecular mechanisms of eukaryotic DNA replication initiation.