Rounding up: Engineering 12-membered rings from the cyclic 11-mer TRAP

The protein TRAP (trp RNA binding attenuation protein) forms a highly thermostable ring-shaped 11-mer. By linking in tandem two, three, or four DNA sequences encoding TRAP monomers, we have engineered new rings that consist of 12 TRAP subunits and bind 12 ligand molecules. The hydrogen bonding pattern and buried surface area within and between subunits are essentially identical between the 11-mer and 12-mer crystal structures. Why do the artificial proteins choose to make single 12-mer rings? The 12-mer rings are highly sterically strained by their peptide linkers and far from thermostable. That proteins choose to adopt a strained conformation of few subunits rather than an unstrained one with 11 subunits demonstrates the importance of entropic factors in controlling protein-protein interactions in general.     

Interaction of the trp RNA-binding attenuation protein (TRAP) with anti-TRAP

The trp RNA-binding attenuation protein (TRAP) negatively regulates expression of the tryptophan biosynthesis genes of Bacillus subtilis. In the presence of tryptophan, TRAP is activated to bind to the 5'-leader region of the trp mRNA resulting in termination prior to the structural genes. In addition, accumulation of uncharged tRNA(Trp) induces synthesis of anti-TRAP (AT), which binds to TRAP and inhibits its function. Both of these proteins consist of oligomers of identical subunits. Here, we characterize the self-association of each of these proteins and the TRAP-AT interaction in free solution using equilibrium and velocity analytical ultracentrifugation. TRAP exists as a stable 11-mer in the absence and in the presence of tryptophan. Tryptophan binding induces a conformational change in TRAP. AT exists in a reversible equilibrium between trimer and dodecamer with an equilibrium constant of approximately 3 x 10(14)M(-3). About 20% of the trimer is incompetent to form dodecamer. The AT equilibrium is slow on the time-scale of the velocity experiment. Formation of TRAP-AT complexes occurs only in the presence of tryptophan. A complex containing one TRAP 11-mer and one AT 12-mer forms with high affinity. At higher ratios of TRAP:AT complexes containing two TRAP 11-mers and one AT 12-mer are detected. A model for the structure of the complex is proposed.     

Mechanisms of allosteric gene regulation by NMR quantification of microsecond-millisecond protein dynamics

The trp RNA-binding attenuation protein (TRAP) is a paradigmatic allosteric protein that regulates the tryptophan biosynthetic genes associated with the trp operon in bacilli. The ring-shaped 11-mer TRAP is activated for recognition of a specific trp-mRNA target by binding up to 11 tryptophan molecules. To characterize the mechanisms of tryptophan-induced TRAP activation, we have performed methyl relaxation dispersion (MRD) nuclear magnetic resonance (NMR) experiments that probe the time-dependent structure of TRAP in the microsecond-to-millisecond "chemical exchange" time window. We find significant side chain flexibility localized to the RNA and tryptophan binding sites of the apo protein and that these dynamics are dramatically reduced upon ligand binding. Analysis of the MRD NMR data provides insights into the structural nature of transiently populated conformations sampled in solution by apo TRAP. The MRD data are inconsistent with global two-state exchange, indicating that conformational sampling in apo TRAP is asynchronous. These findings imply a temporally heterogeneous population of structures that are incompatible with RNA binding and substantiate the study of TRAP as a paradigm for probing and understanding essential dynamics in allosteric, regulatory proteins.     

Computer simulations of the flexibility of a series of synthetic cyclic peptide analogues

The flexibility of a series of cyclic peptides derived from the epitope of a snake toxin is investigated using computer simulation techniques. Molecular dynamics (MD) simulations and vibrational analyses are performed on chemically constrained peptides modeled on the parent protein loop. In the 50 ps MD simulations, large variations in the atomic fluctuations are seen between the peptides, and can be related to the nature of the chemical constraints present in the molecules. Normal mode analyses are performed on energy-minimized configurations derived from the dynamics trajectories. The atomic fluctuations calculated from the normal modes are about 30% of those of the molecular dynamics for the more flexible peptides and 70% for the more constrained molecules. The calculated differences in flexibility between the molecules are much less significant in the harmonic approximation. © 1993 John Wiley & Sons, Inc.     

Trp RNA-Binding Attenuation Protein: Modifying Symmetry and Stability of a Circular Oligomer

Background

Subunit number is amongst the most important structural parameters that determine size, symmetry and geometry of a circular protein oligomer. The L-tryptophan biosynthesis regulator, TRAP, present in several Bacilli, is a good model system for investigating determinants of the oligomeric state. A short segment of C-terminal residues defines whether TRAP forms an 11-mer or 12-mer assembly. To understand which oligomeric state is more stable, we examine the stability of several wild type and mutant TRAP proteins.

Methodology/Principal Findings

Among the wild type B. stearothermophilus, B. halodurans and B. subtilis TRAP, we find that the former is the most stable whilst the latter is the least. Thermal stability of all TRAP is shown to increase with L-tryptophan concentration. We also find that mutant TRAP molecules that are truncated at the C-terminus - and hence induced to form 12-mers, distinct from their 11-mer wild type counterparts - have increased melting temperatures. We show that the same effect can be achieved by a point mutation S72N at a subunit interface, which leads to exclusion of C-terminal residues from the interface. Our findings are supported by dye-based scanning fluorimetry, CD spectroscopy, and by crystal structure and mass spectrometry analysis of the B. subtilis S72N TRAP.

Conclusions/Significance

We conclude that the oligomeric state of a circular protein can be changed by introducing a point mutation at a subunit interface. Exclusion (or deletion) of the C-terminus from the subunit interface has a major impact on properties of TRAP oligomers, making them more stable, and we argue that the cause of these changes is the altered oligomeric state. The more stable TRAP oligomers could be used in potential applications of TRAP in bionanotechnology.     

Characterization of a trp RNA-binding attenuation protein (TRAP) mutant with tryptophan independent RNA binding activity

TRAP (trp RNA-binding attenuation protein) is an 11 subunit RNA-binding protein that regulates expression of genes involved in tryptophan metabolism (trp) in Bacillus subtilis in response to changes in intracellular tryptophan concentration. When activated by binding up to 11 tryptophan residues, TRAP binds to the mRNAs of several trp genes and down-regulates their expression. Recently, a TRAP mutant was found that binds RNA in the absence of tryptophan. In this mutant protein, Thr30, which is part of the tryptophan-binding site, is replaced with Val (T30V). We have compared the RNA-binding properties of T30V and wild-type (WT) TRAP, as well as of a series of hetero-11-mers containing mixtures of WT and T30V TRAP subunits. The most significant difference between the interaction of T30V and WT TRAP with RNA is that the affinity of T30V TRAP is more dependent on ionic strength. Analysis of the hetero-11-mers allowed us to examine how subunits interact within an 11-mer with regard to binding to tryptophan or RNA. Our data suggest that individual subunits retain properties similar to those observed when they are in homo-11-mers and that individual G/UAG triplets within the RNA can bind to TRAP differently.     

Domain motions in bacteriophage T4 lysozyme: A comparison between molecular dynamics and crystallographic data

A comparison of a series of extended molecular dynamics (MD) simulations of bacteriophage T4 lysozyme in solvent with X-ray data is presented. Essential dynamics analyses were used to derive collective fluctuations from both the simulated trajectories and a distribution of crystallographic conformations. In both cases the main collective fluctuations describe domain motions. The protein consists of an N- and C-terminal domain connected by a long helix. The analysis of the distribution of crystallographic conformations reveals that the N-terminal helix rotates together with either of these two domains. The main domain fluctuation describes a closure mode of the two domains in which the N-terminal helix rotates concertedly with the C-terminal domain, while the domain fluctuation with second largest amplitude corresponds to a twisting mode of the two domains, with the N-terminal helix rotating concertedly with the N-terminal domain. For the closure mode, the difference in hinge-bending angle between the most open and most closed X-ray structure along this mode is 49 degrees. In the MD simulation that shows the largest fluctuation along this mode, a rotation of 45 degrees was observed. Although the twisting mode has much less freedom than the closure mode in the distribution of crystallographic conformations, experimental results suggest that it might be functionally important. Interestingly, the twisting mode is sampled more extensively in all MD simulations than it is in the distribution of X-ray conformations. Proteins 31:116–127, 1998. © 1998 Wiley-Liss, Inc.     

A kinetic model for the internal motions of proteins: diffusion between multiple harmonic wells.

The dynamics of collective protein motions derived from Molecular Dynamics simulations have been studied for two small model proteins: initiation factor I and the B1 domain of Protein G. First, we compared the structural fluctuations, obtained by local harmonic approximations in different energy minima, with the ones revealed by large scale molecular dynamics (MD) simulations. It was found that a limited set of harmonic wells can be used to approximate the configurational fluctuations of these proteins, although any single harmonic approximation cannot properly describe their dynamics. Subsequently, the kinetics of the main (essential) collective protein motions were characterized. A dual-diffusion behavior was observed in which a fast type of diffusion switches to a much slower type in a typical time of about 1-3 ps. From these results, the large backbone conformational fluctuations of a protein may be considered as "hopping" between multiple harmonic wells on a basically flat free energy surface.     

Simulating nanoscale functional motions of biomolecules

We are describing efficient dynamics simulation methods for the characterization of functional motion of biomolecules on the nanometer scale. Multivariate statistical methods are widely used to extract and enhance functional collective motions from molecular dynamics (MD) simulations. A dimension reduction in MD is often realized through a principal component analysis (PCA) or a singular value decomposition (SVD) of the trajectory. Normal mode analysis (NMA) is a related collective coordinate space approach, which involves the decomposition of the motion into vibration modes based on an elastic model. Using the myosin motor protein as an example we describe a hybrid technique termed amplified collective motions (ACM) that enhances sampling of conformational space through a combination of normal modes with atomic level MD. Unfortunately, the forced orthogonalization of modes in collective coordinate space leads to complex dependencies that are not necessarily consistent with the symmetry of biological macromolecules and assemblies. In many biological molecules, such as HIV-1 protease, reflective or rotational symmetries are present that are broken using standard orthogonal basis functions. We present a method to compute the plane of reflective symmetry or the axis of rotational symmetry from the trajectory frames. Moreover, we develop an SVD that best approximates the given trajectory while respecting the symmetry. Finally, we describe a local feature analysis (LFA) to construct a topographic representation of functional dynamics in terms of local features. The LFA representations are low-dimensional, and provide a reduced basis set for collective motions, but unlike global collective modes they are sparsely distributed and spatially localized. This yields a more reliable assignment of essential dynamics modes across different MD time windows.     

Creating hetero-11-mers composed of wild-type and mutant subunits to study RNA binding to TRAP.

TRAP (trp RNA-binding attenuation protein) is an RNA-binding protein that regulates expression of the tryptophan biosynthetic genes in Bacillus subtilis by binding to RNA targets that contain multiple GAG and UAG repeats. TRAP is composed of 11 identical subunits arranged symmetrically in a ring. The secondary structure of the protein consists entirely of antiparallel beta-sheets, beta-turns, and loops. We show here that the TRAP 11-mer can be reversibly denatured into unfolded monomers by guanidine hydrochloride. Removing the denaturant allows the protein to spontaneously renature into fully functional 11-mers. Based on this finding, we developed a subunit mixing method to hybridize wild-type and mutant subunits into heteromeric 11-mers by denaturation followed by subunit mixing renaturation. This method allows the study of subunit cooperativity in protein-ligand interaction such as RNA binding. Our data further support and extend the previously proposed two-step model for RNA binding to TRAP by showing that the initiation of binding requires at least one fully active subunit in the protein combined with one fully functional repeat in the RNA. The initiation complex tethers the RNA on the protein, thus allowing cooperative interaction with the remainder of the repeats.     

A statistical approach to the interpretation of molecular dynamics simulations of calmodulin equilibrium dynamics

A sample of 35 independent molecular dynamics (MD) simulations of calmodulin (CaM) equilibrium dynamics was prepared from different but equally plausible initial conditions (20 simulations of the wild-type protein and 15 simulations of the D129N mutant). CaM's radius of gyration and backbone mean-square fluctuations were analyzed for the effect of the D129N mutation, and simulations were compared with experiments. Statistical tests were employed for quantitative comparisons at the desired error level. The computational model predicted statistically significant compaction of CaM relative to the crystal structure, consistent with the results of small-angle X-ray scattering (SAXS) experiments. This effect was not observed in several previously reported studies of (Ca2+)(4)-CaM, which relied on a single MD run. In contrast to radius of gyration, backbone mean-square fluctuations showed a distinctly non-normal and positively skewed distribution for nearly all residues. Furthermore, the D129N mutation affected the backbone dynamics in a complex manner and reduced the mobility of Glu123, Met124, Ile125, Arg126, and Glu127 located in the adjacent alpha-helix G. The implications of these observations for the comparisons of MD simulations with experiments are discussed. The proposed approach may be useful in studies of protein equilibrium dynamics where MD simulations fall short of properly sampling the conformational space, and when the comparison with experiments is affected by the reproducibility of the computational model.     

Minimum sequence requirements for the binding of paromomycin to the rRNA decoding site A

We have recently introduced a computational methodology that combines molecular dynamics (MD) simulations, free-energy calculations, and in vitro binding assays to predict the minimum RNA structural requirements for selective, high-affinity RNA binding to small-molecule ligands. Here, we show that this methodology can be applied to the conformationally flexible aminoglycoside antibiotic paromomycin. A RNA consisting of an 11-mer:10-mer duplex that contains one 16S ribosome RNA decoding A-site bound to paromomycin was simulated for 4 ns. The methodology predicts that the 11-mer:10-mer duplex binds to paromomycin with high affinity, whereas smaller RNA duplexes lose complex stability and the ability to bind paromomycin. The predicted high-affinity binding to paromomycin of the 11-mer:10-mer duplex was confirmed experimentally (EC(50) = 0.28 microM), as well as the inability of smaller complexes to bind. Our simulations show good agreement with experiment for dynamic and structural properties of the isolated A-site, including hydrogen-bonding networks and RNA structural rearrangements upon ligand binding. The results suggest that MD simulations can supplement in vitro methods as a tool for predicting minimum RNA-binding motifs for both small, rigid ligands, and large, flexible ligands when structural information is available.     

Molecular dynamics simulations highlight mobile regions in proteins: A novel suggestion for converting a murine V(H) domain into a more tractable species

The V(H) region of the murine antibody 1F7 has been identified as a single-domain chorismate mutase, but a tendency to denature and aggregate has hampered its biochemical characterization. Standard mutagenesis approaches targeting antibody chain dimerization areas have been exhausted. We describe a new approach to the problem, where we use molecular dynamics (MD) simulations to find the differences between the untractable protein and the known soluble V(H) domain from a llama antibody. MD simulations of proteins yield information on the relative stability and fluctuations of parts of the proteins. By comparing simulation results of two related proteins their differences in stability and fluctuations can be analyzed and may suggest mutations aimed at (de)stabilization of one of the two proteins. For the mouse versus llama simulations, this approach highlights an untried area in the protein which shows increased fluctuations. The replacement of this eight-residue segment with the corresponding llama sequence gave a chimeric mutant which shows significantly decreased fluctuations. We see this as a general scheme to generate suggestions for mutagenesis experiments, not only obviously generalizable to other immunoglobulin domains, but to other protein systems as well.     

Collective motions in proteins: a covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations.

A method is described for identifying collective motions in proteins from molecular dynamics trajectories or normal mode simulations. The method makes use of the covariances of atomic positional fluctuations. It is illustrated by an analysis of the bovine pancreatic trypsin inhibitor. Comparison of the covariance and cross-correlation matrices shows that the relative motions have many similar features in the different simulations. Many regions of the protein, especially regions of secondary structure, move in a correlated manner. Anharmonic effects, which are included in the molecular dynamics simulations but not in the normal analysis, are of some importance in determining the larger scale collective motions, but not the more local fluctuations. Comparisons of molecular dynamics simulations in the present and absence of solvent indicate that the environment is of significance for the long-range motions.     

Molecular modelling of the mass density of single proteins

Using molecular dynamics (MD) simulations, the density of single proteins and its temperature dependence was modelled starting from the experimentally determined protein structure and a generic, transferable force field, without the need of prior parameterization. Although all proteins consist of the same 20 amino acids, their density in aqueous solution varies up to 10% and the thermal expansion coefficient up to twofold. To model the protein density, systematic MD simulations were carried out for 10 proteins with a broad range of densities (1.32-1.43 g/cm(3)) and molecular weights (7-97 kDa). The simulated densities deviated by less than 1.4% from their experimental values that were available for four proteins. Further analyses of protein density showed that it can be essentially described as a consequence of amino acid composition. For five proteins, the density was simulated at different temperatures. The simulated thermal expansion coefficients ranged between 4.3 and 7.1 × 10(-4) K(-1) and were similar to the experimentally determined values of ribonuclease-A and lysozyme (deviations of 2.4 and 14.6%, respectively). Further analyses indicated that the thermal expansion coefficient is linked to the temperature dependence of atomic fluctuations: proteins with a high thermal expansion coefficient show a low increase in flexibility at increasing temperature. A low increase in atomic fluctuations with temperature has been previously described as a possible mechanism of thermostability. Thus, a high thermal expansion coefficient might contribute to protein thermostability.     

Anisotropic atomic motions in high-resolution protein crystallography molecular dynamics simulations

Molecular dynamics (MD) simulations using empirical force fields are popular for the study of proteins. In this work, we compare anisotropic atomic fluctuations in nanosecond-timescale MD simulations with those observed in an ultra-high-resolution crystal structure of crambin. In order to make our comparisons, we have developed a compact graphical technique for assessing agreement between spatial atomic distributions determined by MD simulations and observed anisotropic temperature factors.     

Isolating a trimer intermediate in the self-assembly of E2 protein cage

Understanding the self-assembly mechanism of caged proteins provides clues to develop their potential applications in nanotechnology, such as a nanoscale drug delivery system. The E2 protein from Bacillus stearothermophilus , with a virus-like caged structure, has drawn much attention for its potential application as a nanocapsule. To investigate its self-assembly process from subunits to a spherical protein cage, we truncate the C-terminus of the E2 subunit. The redesigned protein subunit shows dynamic transition between monomer and trimer, but not the integrate 60-mer. The results indicate the role of the trimer as the intermediate and building block during the self-assembly of the E2 protein cage. In combination with the molecular dynamics simulations results, we conclude that the C-terminus modulates the self-assembly of the E2 protein cage from trimer to 60-mer. This investigation elucidates the role of the intersubunit interactions in engineering other functionalities in other caged structure proteins.     

Molecular dynamics simulation of proteins under high pressure: Structure,function and thermodynamics

BackgroundMolecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.Scope of reviewFirst, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.Major conclusionsRecent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.General significanceMD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.     

REACH coarse-grained biomolecular simulation: transferability between different protein structural classes

Coarse graining of protein interactions provides a means of simulating large biological systems. The REACH (Realistic Extension Algorithm via Covariance Hessian) coarse-graining method, in which the force constants of a residue-scale elastic network model are calculated from the variance-covariance matrix obtained from atomistic molecular dynamics (MD) simulation, involves direct mapping between scales without the need for iterative optimization. Here, the transferability of the REACH force field is examined between protein molecules of different structural classes. As test cases, myoglobin (all α), plastocyanin (all β), and dihydrofolate reductase (α/β) are taken. The force constants derived are found to be closely similar in all three proteins. An MD version of REACH is presented, and low-temperature coarse-grained (CG) REACH MD simulations of the three proteins are compared with atomistic MD results. The mean-square fluctuations of the atomistic MD are well reproduced by the CGMD. Model functions for the CG interactions, derived by averaging over the three proteins, are also shown to produce fluctuations in good agreement with the atomistic MD. The results indicate that, similarly to the use of atomistic force fields, it is now possible to use a single, generic REACH force field for all protein studies, without having first to derive parameters from atomistic MD simulation for each individual system studied. The REACH method is thus likely to be a reliable way of determining spatiotemporal motion of a variety of proteins without the need for expensive computation of long atomistic MD simulations.     

Principal component and normal mode analysis of proteins; a quantitative comparison using the GroEL subunit

Principal component analysis (PCA) and normal mode analysis (NMA) have emerged as two invaluable tools for studying conformational changes in proteins. To compare these approaches for studying protein dynamics, we have used a subunit of the GroEL chaperone, whose dynamics is well characterized. We first show that both PCA on trajectories from molecular dynamics (MD) simulations and NMA reveal a general dynamical behavior in agreement with what has previously been described for GroEL. We thus compare the reproducibility of PCA on independent MD runs and subsequently investigate the influence of the length of the MD simulations. We show that there is a relatively poor one-to-one correspondence between eigenvectors obtained from two independent runs and conclude that caution should be taken when analyzing principal components individually. We also observe that increasing the simulation length does not improve the agreement with the experimental structural difference. In fact, relatively short MD simulations are sufficient for this purpose. We observe a rapid convergence of the eigenvectors (after ca. 6 ns). Although there is not always a clear one-to-one correspondence, there is a qualitatively good agreement between the movements described by the first five modes obtained with the three different approaches; PCA, all-atoms NMA, and coarse-grained NMA. It is particularly interesting to relate this to the computational cost of the three methods. The results we obtain on the GroEL subunit contribute to the generalization of robust and reproducible strategies for the study of protein dynamics, using either NMA or PCA of trajectories from MD simulations.