Conformational transitions in protein-protein association: binding of fasciculin-2 to acetylcholinesterase

The neurotoxin fasciculin-2 (FAS2) is a picomolar inhibitor of synaptic acetylcholinesterase (AChE). The dynamics of binding between FAS2 and AChE is influenced by conformational fluctuations both before and after protein encounter. Submicrosecond molecular dynamics trajectories of apo forms of fasciculin, corresponding to different conformational substates, are reported here with reference to the conformational changes of loop I of this three-fingered toxin. This highly flexible loop exhibits an ensemble of conformations within each substate corresponding to its functions. The high energy barrier found between the two major substates leads to transitions that are slow on the timescale of the diffusional encounter of noninteracting FAS2 and AChE. The more stable of the two apo substates may not be the one observed in the complex with AChE. It seems likely that the more stable apo form binds rapidly to AChE and conformational readjustments then occur in the resulting encounter complex.     

A molecular dynamics analysis of protein structural elements

The relation between protein secondary structure and internal motions was examined by using molecular dynamics to calculate positional fluctuations of individual helix, beta-sheet, and loop structural elements in free and substrate-bound hen egg-white lysozyme. The time development of the fluctuations revealed a general correspondence between structure and dynamics; the fluctuations of the helices and beta-sheets converged within the 101 psec period of the simulation and were lower than average in magnitude, while the fluctuations of the loop regions were not converged and were mostly larger than average in magnitude. Notable exceptions to this pattern occurred in the substrate-bound simulation. A loop region (residues 101-107) of the active site cleft had significantly reduced motion due to interactions with the substrate. Moreover, part of a loop and a 3(10) helix (residues of 67-88) not in contact with the substrate showed a marked increase in fluctuations. That these differences in dynamics of free and substrate-bound lysozyme did not result simply from sampling errors was established by an analysis of the variations in the fluctuations of the two halves of the 101 psec simulation of free lysozyme. Concerted transitions of four to five mainchain phi and psi angles between dihedral wells were shown to be responsible for large coordinate shifts in the loops. These transitions displaced six or fewer residues and took place either abruptly, in 1 psec or less, or with a diffusive character over 5-10 psec. Displacements of rigid secondary structures involved longer timescale motions in bound lysozyme; a 0.5 A rms change in the position of a helix occurred over the 55 psec simulation period. This helix reorientation within the protein appears to be a response to substrate binding. There was little correlation between the solvent accessible surface area and the dynamics of the different structural elements.     

Convergence of molecular dynamics simulations of membrane proteins

The central question in evaluating almost any result from a molecular dynamics simulation is whether the calculation has converged. Unfortunately, assessing the ergodicity of a single trajectory is very difficult to do. In this work, we assess the sampling of molecular dynamics simulations of the membrane protein rhodopsin by comparing the results from 26 independent trajectories, each run for 100 ns. By examining principal components and cluster populations, we show that rhodopsin's fluctuations are not well described by 100 ns of dynamics, and that the sampling is not fully converged even for individual loops. The results serve as a reminder of the caution required when interpreting molecular dynamics simulations of macromolecules.     

Stability of a structural scaffold upon activity transfer: X-ray structure of a three fingers chimeric protein

Fasciculin 2 and toxin alpha proteins belong to the same structural family of three-fingered snake toxins. They act on different targets, but in each case the binding region involves residues from loops I and II. The superimposition of the two structures suggests that these functional regions correspond to structurally distinct zones. Loop I, half of loop II and the C-terminal residue of fasciculin 2 were therefore transferred into the toxin alpha. The inhibition constant of the resulting chimera is only 15-fold lower than that of fasciculin 2, and as expected the potency of binding to the toxin alpha target has been lost. In order to understand the structure-function relationship between the chimera and its "parent" molecules, we solved its structure by X-ray crystallography. The protein crystallized in space group P3(1)21 with a=b=58.5 A, and c=62.3 A. The crystal structure was solved by molecular replacement and refined to 2.1 A resolution. The structure belongs to the three-fingered snake toxin family with a core of four disulphide bridges from which emerge the three loops I, II and III. Superimposition of the chimera on fasciculin 2 or toxin alpha revealed an overall fold intermediate between those of the two parent molecules. The regions corresponding to toxin alpha and to fasciculin 2 retained their respective geometries. In addition, the chimera protein displayed a structural behaviour similar to that of fasciculin 2, i.e. dimerization in the crystal structure of fasciculin 2, and the geometry of the region that binds to acetylcholinesterase. In conclusion, this structure shows that the chimera retains the general structural characteristics of three-fingered toxins, and the structural specificity of the transferred function.     

A molecular dynamics simulation of bacteriophage T4 lysozyme.

An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. The simulation was carried out with all hydrogen atoms modeled explicitly, the inclusion of all 152 crystallographic waters and at a temperature of 300 K. Temporal analysis of the trajectory versus energy, hydrogen bond stability, r.m.s. deviation from the starting crystal structure and radius of gyration, demonstrates that the simulation was both stable and representative of the average experimental structure. Average structural properties were calculated from the enzyme trajectory and compared with the crystal structure. The mean value of the C alpha displacements of the average simulated structure from the X-ray structure was 1.1 +/- 0.1 A; differences of the backbone phi and psi angles between the average simulated structure and the crystal structure were also examined. Thermal-B factors were calculated from the simulation for heavy and backbone atoms and both were in good agreement with experimental values. Relationships between protein secondary structure elements and internal motions were studied by examining the positional fluctuations of individual helix, sheet and turn structures. The structural integrity in the secondary structure units was preserved throughout the simulation; however, the A helix did show some unusually high atomic fluctuations. The largest backbone atom r.m.s. fluctuations were found in non-secondary structure regions; similar results were observed for r.m.s. fluctuations of non-secondary structure phi and psi angles. In general, the calculated values of r.m.s. fluctuations were quite small for the secondary structure elements. In contrast, surface loops and turns exhibited much larger values, being able to sample larger regions of conformational space. The C alpha difference distance matrix and super-positioning analyses comparing the X-ray structure with the average dynamics structure suggest that a 'hinge-bending' motion occurs between the N- and C-terminal domains.     

Simulations of a protein crystal: explicit treatment of crystallization conditions links theory and experiment in the streptavidin-biotin complex

A 250 ns molecular dynamics simulation of the biotin-liganded streptavidin crystal lattice, including cryoprotectant molecules and crystallization salts, is compared to a 250 ns simulation of the lattice solvated with pure water. The simulation using detailed crystallization conditions preserves the initial X-ray structure better than the simulation using pure water, even though the protein molecules display comparable mobility in either simulation. Atomic fluctuations computed from the simulation with crystallization conditions closely reproduce fluctuations derived from experimental temperature factors (correlation coefficient of 0.88, omitting two N-terminal residues with very high experimental B-factors). In contrast, fluctuations calculated from the simulation with pure water were less accurate, particularly for two of the streptavidin loops exposed to solvent in the crystal lattice. Finally, we obtain good agreement between the water and cryoprotectant densities obtained from the simulated crystallization conditions and the electron density due to solvent molecules in the X-ray structure. Our results suggest that detailed lattice simulations with realistic crystallization conditions can be used to assess potential function parameters, validate simulation protocols, and obtain valuable insights that solution-phase simulations do not easily provide. We anticipate that this will prove to be a powerful strategy for molecular dynamics simulations of biomolecules.     

A comparison of structural and dynamic properties of different simulation methods applied to SH3.

The dynamic and static properties of molecular dynamics simulations using various methods for treating solvent were compared. The SH3 protein domain was chosen as a test case because of its small size and high surface-to-volume ratio. The simulations were analyzed in structural terms by examining crystal packing, distribution of polar residues, and conservation of secondary structure. In addition, the "essential dynamics" method was applied to compare each of the molecular dynamics trajectories with a full solvent simulation. This method proved to be a powerful tool for the comparison of large concerted atomic motions in SH3. It identified methods of simulation that yielded significantly different dynamic properties compared to the full solvent simulation. Simulating SH3 using the stochastic dynamics algorithm with a vacuum (reduced charge) force field produced properties close to those of the full solvent simulation. The application of a recently described solvation term did not improve the dynamic properties. The large concerted atomic motions in the full solvent simulation as revealed by the essential dynamics method were analyzed for possible biological implications. Two loops, which have been shown to be involved in ligand binding, were seen to move in concert to open and close the ligand-binding site.     

Molecular dynamics simulations of human glutathione transferase P1-1: conformational fluctuations of the apo-structure.

We have investigated by molecular dynamics simulations the conformational fluctuations of the monomer of human apo-glutathione transferase P1-1. After attainment of steady-state dynamics, the structural fluctuations involve mainly the protein segments that participate also in the holo-apo transition discussed in the accompanying article (Stella et al., 1999:37:1-9.). The most mobile region is the C-terminal segment of helix 2. In contrast, helices 1, 6, 7, and 8 constitute a relatively rigid protein core. An "essential dynamics" analysis of the simulation shows that the largest fluctuations involve specific regions of glutathione transferases. In such regions, atomic motions are correlated. Motions of helix 2 are accounted for by the second most prominent principal component, which reveals a fluctuation between two distinct conformations. The residues that constitute the H-site undergo a breathing motion, possibly relevant during the binding of hydrophobic cosubstrates. Based on our simulation, several experimental findings can be rationalized, including the viscosity-dependent reactivity of Cys 47 and Cys 101 as well as the selective proteolysis of the peptide bond between Lys 44 and Ala 45. We have also modeled the structural changes that lead to the formation of an intrachain disulfide bridge between cysteines 47 and 101 and to the inactivation of the enzyme. The resulting structure maintains essentially the native fold except for helix 2, which closes the G-site. Proteins 1999;37:10-19.     

Motions and structural variability within toxins: implication for their use as scaffolds for protein engineering

Animal toxins are small proteins built on the basis of a few disulfide bonded frameworks. Because of their high variability in sequence and biologic function, these proteins are now used as templates for protein engineering. Here we report the extensive characterization of the structure and dynamics of two toxin folds, the "three-finger" fold and the short alpha/beta scorpion fold found in snake and scorpion venoms, respectively. These two folds have a very different architecture; the short alpha/beta scorpion fold is highly compact, whereas the "three-finger" fold is a beta structure presenting large flexible loops. First, the crystal structure of the snake toxin alpha was solved at 1.8-A resolution. Then, long molecular dynamics simulations (10 ns) in water boxes of the snake toxin alpha and the scorpion charybdotoxin were performed, starting either from the crystal or the solution structure. For both proteins, the crystal structure is stabilized by more hydrogen bonds than the solution structure, and the trajectory starting from the X-ray structure is more stable than the trajectory started from the NMR structure. The trajectories started from the X-ray structure are in agreement with the experimental NMR and X-ray data about the protein dynamics. Both proteins exhibit fast motions with an amplitude correlated to their secondary structure. In contrast, slower motions are essentially only observed in toxin alpha. The regions submitted to rare motions during the simulations are those that exhibit millisecond time-scale motions. Lastly, the structural variations within each fold family are described. The localization and the amplitude of these variations suggest that the regions presenting large-scale motions should be those tolerant to large insertions or deletions.     

Assembly and stability of Salmonella enterica ser. Typhi TolC protein in POPE and DMPE

In this work we assessed the suitability of two different lipid membranes for the simulation of a TolC protein from Salmonella enterica serovar Typhi. The TolC protein family is found in many pathogenic Gram-negative bacteria including Vibrio cholera and Pseudomonas aeruginosa and acts as an outer membrane channel for expulsion of drug and toxin from the cell. In S. typhi, the causative agent for typhoid fever, the TolC outer membrane protein is an antigen for the pathogen. The lipid environment is an important modulator of membrane protein structure and function. We evaluated the conformation of the TolC protein in the presence of DMPE and POPE bilayers using molecular dynamics simulation. The S. typhi TolC protein exhibited similar conformational dynamics to TolC and its homologues. Conformational flexibility of the protein is seen in the C-terminal, extracellular loops, and α-helical region. Despite differences in the two lipids, significant similarities in the motion of the protein in POPE and DMPE were observed, including the rotational motion of the C-terminal residues and the partially open extracellular loops. However, analysis of the trajectories demonstrated effects of hydrophobic matching of the TolC protein in the membrane, particularly in the lengthening of the lipids and subtle movements of the protein’s β-barrel towards the lower leaflet in DMPE. The study exhibited the use of molecular dynamics simulation in revealing the differential effect of membrane proteins and lipids on each other. In this study, POPE is potentially a more suitable model for future simulation of the S. typhi TolC protein.     

Molecular dynamics of myoglobin at 298 degrees K. Results from a 300-ps computer simulation.

We have carried out a very long (300 ps) molecular dynamics simulation of the protein myoglobin. This trajectory is approximately three times longer than the longest previous molecular dynamics simulation of a protein, and ten times longer than protein simulations of comparable size (1,423 atoms in our model). Here we report results from this long simulation concerning the average structure, the mean square fluctuations of atoms about the average structure, and the nuclear magnetic resonance order parameters for various groups in myoglobin. The results demonstrate that the average coordinates change very slowly during the simulation. The relative atomic mobilities are well described by the simulation. For both the mean square atomic fluctuations and the order parameters, however, there are significant quantitative differences when values calculated using shorter portions of the trajectory are compared with results obtained for the entire 300-ps simulation. The implications of this result for obtaining converged properties from protein molecular dynamics simulations for comparison with experiment are discussed.     

Molecular dynamics simulations of human rhinovirus and an antiviral compound

The human rhinovirus 14 (HRV14) protomer, with or without the antiviral compound WIN 52084s, was simulated using molecular dynamics and rotational symmetry boundary conditions to model the effect of the entire icosahedral capsid. The protein asymmetrical unit, comprising four capsid proteins (VP1, VP2, VP3, and VP4) and two calcium ions, was solvated both on the exterior and the interior to fill the inside of the capsid. The stability of the simulations of this large system (~800 residues and 6,650 water molecules) is comparable to more conventional globular protein simulations. The influence of the antiviral compound on compressibility and positional fluctuations is reported. The compressibility, estimated from the density fluctuations in the region of the binding pocket, was found to be greater with WIN 52084s bound than without the drug, substantiating previous computations on reduced viral systems. An increase in compressibility correlates with an entropically more favorable system. In contrast to the increase in density fluctuations and compressibility, the positional fluctuations decreased dramatically for the external loops of VP1 and the N-terminus of VP3 when WIN 52084s is bound. Most of these VP1 and VP3 loops are found near the fivefold axis, a region whose mobility was not considered in reduced systems, but can be observed with this simulation of the full viral protomer. Altered loop flexibility is consistent with changes in proteolytic sensitivity observed experimentally. Moreover, decreased flexibility in these intraprotomeric loops is noteworthy since the externalization of VP4, part of VP1, and RNA during the uncoating process is thought to involve areas near the fivefold axis. Both the decrease in positional fluctuations at the fivefold axis and the increase in compressibility near the WIN pocket are discussed in relationship to the antiviral activity of stabilizing the virus against uncoating.     

Domain analysis of fatty acid synthase protein (NP_217040) from Mycobacterium tuberculosis H37Rv--a bioinformatics study

Different domains of fatty acid synthase (FAS) protein of Mycobacterium tuberculosis H37Rv, involved in mycolic acid synthesis were analyzed using various bioinformatics tools. Based on different database searches (CDD and Pfam), FAS protein of Mycobacterium tuberculosis was grouped into eight domains, five of which showed close similarity with pdb templates (1MLA, 1IQ6A, 2BMOA, and 1J3NA). Based on the PSI blast analysis, 3D structures of only five domains were predicted using MODELLER software, and loop modeling was done for only those regions that were predicted as loops by predict protein server. Compared to the original structure, the loop modeled structure showed a lower DOPE score value for FAS protein. The X-ray determined templates that were used for predicting the 3D structure suggest that, FAS protein has "Malonyl-coenzyme A-Hydratase-Nitrobenzene dioxygenase-3-oxoacyl-(acp) synthase" activity. Accuracy of the prediction of 3D structure of different domains of FAS protein was further validated by Ramachandran plot and PROCHECK (G-value).     

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.     

Importance of the CMAP correction to the CHARMM22 protein force field: dynamics of hen lysozyme

The recently developed CMAP correction to the CHARMM22 force field (C22) is evaluated from 25 ns molecular dynamics simulations on hen lysozyme. Substantial deviations from experimental backbone root mean-square fluctuations and N-H NMR order parameters obtained in the C22 trajectories (especially in the loops) are eliminated by the CMAP correction. Thus, the C22/CMAP force field yields improved dynamical and structural properties of proteins in molecular dynamics simulations.     

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.     

Defining the nature of thermal intermediate in 3 state folding proteins: apoflavodoxin, a study case

The early stages of the thermal unfolding of apoflavodoxin have been determined by using atomistic multi microsecond-scale molecular dynamics (MD) simulations complemented with a variety of experimental techniques. Results strongly suggest that the intermediate is reached very early in the thermal unfolding process and that it has the properties of an "activated" form of the native state, where thermal fluctuations in the loops break loop-loop contacts. The unrestrained loops gain then kinetic energy corrupting short secondary structure elements without corrupting the core of the protein. The MD-derived ensembles agree with experimental observables and draw a picture of the intermediate state inconsistent with a well-defined structure and characteristic of a typical partially disordered protein. Our results allow us to speculate that proteins with a well packed core connected by long loops might behave as partially disordered proteins under native conditions, or alternatively behave as three state folders. Small details in the sequence, easily tunable by evolution, can yield to one or the other type of proteins.     

Molecular modelling and dynamics of CA2 missense mutations causative to carbonic anhydrase 2 deficiency syndrome

Abstract

Carbonic anhydrase 2 (CA2) enzyme deficiency caused by CA2 gene mutations is an inherited disorder characterized by symptoms like osteopetrosis, renal tubular acidosis, and cerebral calcification. This study has collected the CA2 deficiency causal missense mutations and assessed their pathogenicity using diverse computational programs. The 3D protein models for all missense mutations were built, and analyzed for structural divergence, protein stability, and molecular dynamics properties. We found M-CAP as the most sensitive prediction method to measure the deleterious potential of CA2 missense mutations. Free energy dynamics of tertiary structure models of CA2 mutants with DUET, mCSM, and SDM based consensus methods predicted only 50% of the variants as destabilizing. Superimposition of native and mutant CA2 models revealed the minor structural fluctuations at the amino acid residue level but not at the whole protein structure level. Near native molecular dynamic simulation analysis indicated that CA2 causative missense variants result in residue level fluctuation pattern in the protein structure. This study expands the understanding of genotype-protein phenotype correlations underlying CA2 variant pathogenicity and presents a potential avenue for modifying the CA2 deficiency by targeting biophysical structural features of CA2 protein.

Communicated by Ramaswamy H. Sarma     

Molecular dynamics effects on protein electrostatics

Electrostatic calculations have been carried out on a number of structural conformers of tuna cytochrome c. Conformers were generated using molecular dynamics simulations with a range of solvent simulating, macroscopic dielectric formalisms, and one solvent model that explicitly included solvent water molecules. Structures generated using the lowest dielectric models were relatively tight, with side chains collapsed on the surface, while those from the higher dielectric models had more internal and external fluidity, with surface side chains exploring a fuller range of conformational space. The average structure generated with the explicitly solvated model corresponded most closely with the crystal structure. Individual pK values, overall titration curves, and electrostatic potential surfaces were calculated for average structures and structures along each simulation. Differences between structural conformers within each simulation give rise to substantial changes in calculated local electrostatic interactions, resulting in pK value fluctuations for individual sites in the protein that vary by 0.3-2.0 pK units from the calculated time average. These variations are due to the thermal side chain reorientations that produce fluctuations in charge site separations. Properties like overall titration curves and pH dependent stability are not as sensitive to side chain fluctuations within a simulation, but there are substantial effects between simulations due to marked differences in average side chain behavior. These findings underscore the importance of proper dielectric formalism in molecular dynamics simulations when used to generate alternate solution structures from a crystal structure, and suggest that conformers significantly removed from the average structure have altered electrostatic properties that may prove important in episodic protein properties such as catalysis.