Computational investigations of structural changes resulting from point mutations in a collagen-like peptide

The results of 0.5-1.0 ns molecular dynamics simulations of the collagen-like peptides [(POG)4(POA)(POG)4]3 and [(POG)9]3 (POG: proline-hydroxyproline-glycine) are presented. All simulations were performed using the AMBER-94 molecular mechanical force field with a shell of TIP3P waters surrounding the peptides. The initial geometries for the collagen-like peptides included an x-ray crystallographic structure, a computer-generated structure, a [(POG)9]3 structure modeled from the x-ray structure, and the x-ray structure with crystallographic waters replaced with a shell of modeled TIP3P waters. We examined the molecular dynamics peptide residue rms deviation fluctuations, dihedral angles, molecular and chain end-to-end distances, helical parameters, and peptide-peptide and peptide-solvent hydrogen-bonding patterns. Our molecular dynamics simulations of [(POG)4(POA)(POG)4]3 show average structures and internal coordinates similar to the x-ray crystallographic structure. Our results demonstrate that molecular dynamics can be used to reproduce the experimental structures of collagen-like peptides. We have demonstrated the feasibility of using the AMBER-94 molecular mechanical force field, which was parameterized to model nucleic acids and globular proteins, for fibril proteins. We provide a new interpretation of peptide-solvent hydrogen bonding and a peptide-peptide hydrogen bonding pattern not previously reported in x-ray studies. Last, we report on the differences; in particular with respect to main-chain dihedral angles and hydrogen bonding, between the native and mutant collagen-like peptides.     

3 Nsec molecular dynamics simulation of the protein ubiquitin and comparison with X-ray crystal and solution NMR structures.

Mainly due to computational limitations, past protein molecular dynamics simulations have rarely been extended to 300 psec; we are not aware of any published results beyond 350 psec. The present work compares a 3000 psec simulation of the protein ubiquitin with the available x-ray crystallographic and solution NMR structures. Aside from experimental structure availability, ubiquitin was studied because of its relatively small size (76 amino acids) and lack of disulfide bridges. An implicit solvent model was used except for explicit treatment of waters of crystallization. We found that the simulated average structure retains most of the character of the starting x-ray crystal structure. In two highly surface accessible regions, the simulation was not in agreement with the x-ray structure. In addition, there are six backbone-backbone hydrogen bonds that are in conflict between the solution NMR and x-ray crystallographic structures; two are bonds that the NMR does not locate, and four are ones that the two methods disagree upon the donor. Concerning these six backbone-backbone hydrogen bonds, the present simulation agrees with the solution NMR structure in five out-of-the six cases, in that if a hydrogen bond is present in the x-ray structure and not in the NMR structure, the bond breaks within 700 psec. Of the two hydrogen bonds that are found in the NMR structure and not in the x-ray structure, one forms at 1400 psec and the other forms rarely. The present results suggest that relatively long molecular dynamics simulations, that use protein x-ray crystal coordinates for the starting structure and a computationally efficient solvent representation, may be used to gain an understanding of conformational and dynamic differences between the solid-crystal and dilute-solution states.     

The open nucleotide pocket of the profilin/actin x-ray structure is unstable and closes in the absence of profilin

The open nucleotide pocket conformation of actin in the profilin:actinCaATP x-ray structure has been hypothesized to be a crucial intermediate for nucleotide exchange in the actin depolymerization/polymerization cycle. The requirement for ancillary modification of actin for crystallization leads to ambiguities in this interpretation, however. We have used molecular dynamics simulations to model the thermodynamic properties of the actin x-ray structure, outside the crystal lattice, in an aqueous environment with profilin removed. Our simulations show that the open-nucleotide-pocket, profilin-free structure is actually unstable, and closes. The coordination of actin to the nucleotide in the molecular-dynamics-derived closed structure is virtually identical to that in the closed profilin:actinSrATP x-ray structure. Thus, there is currently no thermodynamically stable structure representing the open-nucleotide-pocket state of actin.     

A computer model-built structure of the T. congolense trypanothione reductase in complex with NADP

A three-dimensional structure is engineered for the Trypanosoma congolense trypanothione reductase (TpR) using the sequence homology with glutathione reductase (GR) and lipoamide dehydrogenase, molecular graphics, energy optimization and molecular dynamics techniques. The model was extended to include the complex with the coenzyme nicotinamide adenine dinucleotide phosphate (NADP). The TpR-NADP structure is compared with X-ray data from the glutathione reductase complex with the reduced NADP (NADPH). A model of TpR-NADP including the trypanothione substrate is presented, and an electron-transfer mechanism is proposed.     

Active site of trypanothione reductase. A target for rational drug design.

The X-ray crystal structure of the enzyme trypanothione reductase, isolated from the trypanosomatid organism Crithidia fasciculata, has been solved by molecular replacement. The search model was the crystal structure of human glutathione reductase that shares approximately 40% sequence identity. The trypanosomal enzyme crystallizes in the tetragonal space group P4(1) with unit cell lengths of a = 128.9 A and c = 92.3 A. The asymmetric unit consists of a homodimer of approximate molecular mass 108 kDa. We present the structural detail of the active site as derived from the crystallographic model obtained at an intermediate stage of the analysis using diffraction data to 2.8 A resolution with an R-factor of 23.2%. This model has root-mean-square deviations from ideal geometry of 0.026 A for bond lengths and 4.7 degrees for bond angles. The trypanosomid enzyme assumes a similar biological function to glutathione reductase and, although similar in topology to human glutathione reductase, has an enlarged active site and a number of amino acid differences, steric and electrostatic, which allows it to process only the unique substrate trypanothione and not glutathione. This protein represents a prime target for chemotherapy of several debilitating tropical diseases caused by protozoan parasites belonging to the genera Trypanosoma and Leishmania. The structural differences between the parasite and host enzymes and their substrates thus provides a rational basis for the design of new drugs active against trypanosomes. In addition, our model explains the results of site-directed mutagenesis experiments, carried out on recombinant trypanothione reductase and glutathione reductases, designed by consideration of the crystal structure of human glutathione reductase.     

Probing protein structure by solvent perturbation of NMR spectra: the surface accessibility of bovine pancreatic trypsin inhibitor.

In the absence of specific interactions, the relative attenuation of protein NMR signals due to added stable free radicals such as TEMPOL should reflect the solvent accessibility of the molecular surface. The quantitative correlation between observed attenuation and surface accessibility was investigated with a model system, i.e., the small protein bovine pancreatic trypsin inhibitor. A detailed discussion is presented on the reliability and limits of the approach, and guidelines are provided for data acquisition, treatment, and interpretation. The NMR-derived accessibilities are compared with those obtained from x-ray diffraction and molecular dynamics data. Although the time-averaged accessibilities from molecular dynamics are ideally suited to fit the NMR data, better agreement was observed between the paramagnetic attenuations of the fingerprint cross-peaks of homonuclear proton spectra and the total NH and H alpha accessibilities calculated from x-ray coordinates, than from time-averaged molecular dynamics simulations. In addition, the solvent perturbation response appears to be a promising approach for detecting the thermal conformational evolution of secondary structure elements in proteins.     

Molecular dynamics computations and solid state nuclear magnetic resonance of the gramicidin cation channel.

This paper reports on a coupled approach to determining the structure of the gramicidin A ion channel, utilizing solid state nuclear magnetic resonance (NMR) of isotopically labeled gramicidin channels aligned parallel to the magnetic field direction, and molecular dynamics (MD). MD computations using an idealized right-handed beta-helix as a starting point produce a refined molecular structure that is in excellent agreement with atomic resolution solid state NMR data. The data provided by NMR and MD are complementary to each other. When applied in a coordinated manner they provide a powerful approach to structure determination in molecular systems not readily amenable to x-ray diffraction.     

Simulation-based fitting of protein-protein interaction potentials to SAXS experiments

We present a new method for computing interaction potentials of solvated proteins directly from small-angle x-ray scattering data. An ensemble of proteins is modeled by Monte Carlo or molecular dynamics simulation. The global x-ray scattering of the whole model ensemble is then computed at each snapshot of the simulation, and averaged to obtain the x-ray scattering intensity. Finally, the interaction potential parameters are adjusted by an optimization algorithm, and the procedure is iterated until the best agreement between simulation and experiment is obtained. This new approach obviates the need for approximations that must be made in simplified analytical models. We apply the method to lambda repressor fragment 6-85 and fyn-SH3. With the increased availability of fast computer clusters, Monte Carlo and molecular dynamics analysis using residue-level or even atomistic potentials may soon become feasible.     

Short-range order and collective dynamics of DMPC bilayers: a comparison between molecular dynamics simulations, X-ray, and neutron scattering experiments

We present an extensive comparison of short-range order and short wavelength dynamics of a hydrated phospholipid bilayer derived by molecular dynamics simulations, elastic x-ray, and inelastic neutron scattering experiments. The quantities that are compared between simulation and experiment include static and dynamic structure factors, reciprocal space mappings, and electron density profiles. We show that the simultaneous use of molecular dynamics and diffraction data can help to extract real space properties like the area per lipid and the lipid chain ordering from experimental data. In addition, we assert that the interchain distance can be computed to high accuracy from the interchain correlation peak of the structure factor. Moreover, it is found that the position of the interchain correlation peak is not affected by the area per lipid, while its correlation length decreases linearly with the area per lipid. This finding allows us to relate a property of the structure factor quantitatively to the area per lipid. Finally, the short wavelength dynamics obtained from the simulations and from inelastic neutron scattering are analyzed and compared. The conventional interpretation in terms of the three-effective-eigenmode model is found to be only partly suitable to describe the complex fluid dynamics of lipid chains.     

Stability and structure of oligomers of the Alzheimer peptide Abeta16-22: from the dimer to the 32-mer

Several neurodegenerative diseases such as Alzheimer's, Parkinson's, and Huntington's diseases are associated with amyloid fibrils formed by different polypeptides. We probe the structure and stability of oligomers of different sizes of the fragment Abeta(16-22) of the Alzheimer beta-amyloid peptide using atomic-detail molecular dynamics simulations with explicit solvent. We find that only large oligomers form a stable beta-sheet aggregate, the minimum nucleus size being of the order of 8-16 peptides. This effect is attributed to better hydrophobic contacts and a better shielding of backbone-backbone hydrogen bonds from the solvent in bigger assemblies. Moreover, the observed stability of beta-sheet aggregates with a different number of layers can be explained on the basis of their solvent-accessible surface area. Depending on the stacking interface between the sheets, we observe straight or twisted structures, which could be linked to the experimentally observed polymorphism of amyloid fibrils. To compare our 32-mer structure to experimental data, we calculate its x-ray diffraction pattern. Good agreement is found between experimentally and theoretically determined reflections, suggesting that our model indeed closely resembles the structures found in vitro.     

Homology modeling and molecular dynamics simulation studies of an inward rectifier potassium channel

A homology model has been generated for the pore-forming domain of Kir6.2, a component of an ATP-sensitive K channel, based on the x-ray structure of the bacterial channel KcsA. Analysis of the lipid-exposed and pore-lining surfaces of the model reveals them to be compatible with the known features of membrane proteins and Kir channels, respectively. The Kir6.2 homology model was used as the starting point for nanosecond-duration molecular dynamics simulations in a solvated phospholipid bilayer. The overall drift from the model structure was comparable to that seen for KcsA in previous similar simulations. Preliminary analysis of the interactions of the Kir6.2 channel model with K(+) ions and water molecules during these simulations suggests that concerted single-file motion of K(+) ions and water through the selectivity filter occurs. This is similar to such motion observed in simulations of KcsA. This suggests that a single-filing mechanism is conserved between different K channel structures and may be robust to changes in simulation details. Comparison of Kir6.2 and KcsA suggests some degree of flexibility in the filter, thus complicating models of ion selectivity based upon a rigid filter.     

The Structure of Human Apolipoprotein A-IV as Revealed by Stable Isotope-assisted Cross-linking,Molecular Dynamics,and Small Angle X-ray Scattering

Apolipoprotein (apo)A-IV plays important roles in dietary lipid and glucose metabolism, and knowledge of its structure is required to fully understand the molecular basis of these functions. However, typical of the entire class of exchangeable apolipoproteins, its dynamic nature and affinity for lipid has posed challenges to traditional high resolution structural approaches. We previously reported an x-ray crystal structure of a dimeric truncation mutant of apoA-IV, which showed a unique helix-swapping molecular interface. Unfortunately, the structures of the N and C termini that are important for lipid binding were not visualized. To build a more complete model, we used chemical cross-linking to derive distance constraints across the full-length protein. The approach was enhanced with stable isotope labeling to overcome ambiguities in determining molecular span of the cross-links given the remarkable similarities in the monomeric and dimeric apoA-IV structures. Using 51 distance constraints, we created a starting model for full-length monomeric apoA-IV and then subjected it to two modeling approaches: (i) molecular dynamics simulations and (ii) fitting to small angle x-ray scattering data. This resulted in the most detailed models yet for lipid-free monomeric or dimeric apoA-IV. Importantly, these models were of sufficient detail to direct the experimental identification of new functional residues that participate in a “clasp” mechanism to modulate apoA-IV lipid affinity. The isotope-assisted cross-linking approach should prove useful for further study of this family of apolipoproteins in both the lipid-free and -bound states.     

Hierarchical approach to predicting permeation in ion channels

A hierarchical computational strategy combining molecular modeling, electrostatics calculations, molecular dynamics, and Brownian dynamics simulations is developed and implemented to compute electrophysiologically measurable properties of the KcsA potassium channel. Models for a series of channels with different pore sizes are developed from the known x-ray structure, using insights into the gating conformational changes as suggested by a variety of published experiments. Information on the pH dependence of the channel gating is incorporated into the calculation of potential profiles for K(+) ions inside the channel, which are then combined with K(+) ion mobilities inside the channel, as computed by molecular dynamics simulations, to provide inputs into Brownian dynamics simulations for computing ion fluxes. The open model structure has a conductance of approximately 110 pS under symmetric 250 mM K(+) conditions, in reasonable agreement with experiments for the largest conducting substate. The dimensions of this channel are consistent with electrophysiologically determined size dependence of quaternary ammonium ion blocking from the intracellular end of this channel as well as with direct structural evidence that tetrabutylammonium ions can enter into the interior cavity of the channel. Realistic values of Ussing flux ratio exponents, distribution of ions within the channel, and shapes of the current-voltage and current-concentration curves are obtained. The Brownian dynamics calculations suggest passage of ions through the selectivity filter proceeds by a "knock-off" mechanism involving three ions, as has been previously inferred from functional and structural studies of barium ion blocking. These results suggest that the present calculations capture the essential nature of K(+) ion permeation in the KcsA channel and provide a proof-of-concept for the integrated microscopic/mesoscopic multitiered approach for predicting ion channel function from structure, which can be applied to other channel structures.     

Solution structure of the glycosylphosphatidylinositol membrane anchor glycan of Trypanosoma brucei variant surface glycoprotein

The average solution conformation of the glycosylphosphatidylinositol (GPI) membrane anchor of Trypanosoma brucei variant surface glycoprotein (VSG) has been determined by using a combination of two-dimensional 1H-1H NMR methods together with molecular orbital calculations and restrained molecular dynamics simulations. This allows the generation of a model to describe the orientation of the glycan with respect to the membrane. This shows that the glycan exists in an extended configuration along the plane of the membrane and spans an area of 600 A2, which is similar to the cross-sectional area of a monomeric N-terminal VSG domain. Taken together, these observations suggest a possible space-filling role for the GPI anchor that may maintain the integrity of the VSG coat. The potential importance of the GPI glycan as a chemotherapeutic target is discussed in light of these observations.     

Elbow flexibility and ligand-induced domain rearrangements in antibody Fab NC6.8: large effects of a small hapten

Four 700-ps molecular dynamics simulations were carried out to analyze the structural dynamics of the antigen-binding antibody fragment NC6.8, which is known to exhibit large structural changes upon complexation. The first simulation was started from the x-ray structure of the uncomplexed Fab and produced trajectory averages that closely match the crystallographic results. It allowed assessment of the flexibility of the Fab, revealing an elbow motion of the variable domains with respect to the constant domains. The second simulation was started from the uncomplexed x-ray structure after insertion of the ligand into the binding site. This perturbation resulted in a significantly altered trajectory, with quaternary structural changes corresponding in many aspects to the experimental differences between complexed and uncomplexed state. The observed trend toward a smaller elbow angle and a higher flexion of the H-chain could also be seen in the third simulation, which was started from the x-ray structure of the complex. The changes were revealed to be a clear consequence of the complexation with the ligand because in the fourth simulation (started from the experimental complex structure after removal of the hapten) the Fab remained close to its initial structure. Analyses of the quaternary structure and the binding site of Fab NC6.8 are presented for all four simulations, and possible interpretations are discussed.     

The use of NMR chemical shifts to analyse the MD trajectories: simulation of bovine pancreatic trypsin inhibitor dynamics in water as a test case for solvent influences.

In this paper the NMR secondary chemical shifts, that are estimated from a set of 3D-structures, are compared with the observed ones to appraise the behaviour of a known x-ray diffraction structure (of the bovine pancreatic trypsin inhibitor protein) when various molecular dynamics are applied. The results of a 200 ps molecular dynamics under various conditions are analysed and different ways to modify the molecular dynamics are considered. With the purpose of avoiding the time-consuming explicit representation of the solvent (water) molecules, an attempt was made to understand the role of the solvent and to develop an implicit representation, which may be refined. A simulation of hydrophobic effects in an aqueous environment is also proposed which seems to provide a better approximation of the observed solution structure of the protein.     

Molecular dynamics as a tool to detect protein foldability. A mutant of domain B1 of protein G with non-native secondary structure propensities

The usefulness of molecular dynamics to assess the structural integrity of mutants containing several mutations has been investigated. Our goal was to determine whether molecular dynamics would be able to discriminate mutants of a protein having a close-to-wild-type fold, from those that are not folded under the same conditions. We used as a model the B1 domain of protein G in which we replaced the unique central alpha-helix by the sequence of the second beta-hairpin, which has a strong intrinsic propensity to form this secondary structure in solution. In the resulting protein, one-third of the secondary structure has been replaced by a non-native one. Models of the mutants were built based on the three-dimensional structure of the wild-type GB1 domain. During 2 ns of molecular dynamics simulations on these models, mutants containing up to 10 mutations in the helix retained the native fold, while another mutant with an additional mutation unfolded. This result is in agreement with our circular dichroism and NMR experiments, which indicated that the former mutants fold into a structure similar to the wild-type, as opposed to the latter mutant which is partly unfolded. Additionally, a mutant containing six mutations scattered through the surface of the domain, and which is unfolded, was also detected by the simulation. This study suggests that molecular dynamics calculations could be performed on molecular models of mutants of a protein to evaluate their foldability, prior to a mutagenesis experiment.     

The presence of a histidine-aspartic acid pair in the active site of 2-hydroxyacid dehydrogenases. X-ray refinement of cytoplasmic malate dehydrogenase

The structure of cytoplasmic malate dehydrogenase has been partially refined by crystallographic least squares methods. Using x-ray phases based on the refined coordinates, analysis of the resultant electron density maps has led to a new model of cytoplasmic malate dehydrogenase and a tentative "x-ray sequence." The two crystallographically independent subunits comprising the dimeric enzyme are nearly identical in structure and are related to each other by roughly 2-fold rotational symmetry. The best fit of the molecular structure of cytoplasmic malate dehydrogenase to that of lactate dehydrogenase has been obtained by least squares methods. The active sites of these two enzymes contain similarly oriented His-Asp pairs linked by a hydrogen bond which may function as a proton relay system during catalysis. This pair could also provide an explanation for the relatively stronger binding by cytoplasmic malate dehydrogenase and lactate dehydrogenase of NADH versus NAD. Similar His-Asp pairs have been observed in the serine proteases, thermolysin, and phospholipase A2, and the His-Asp pair may play a similar functional role in all of these enzymes.     

Heme structure and orientation in single monolayers of cytochrome c on polar and nonpolar soft surfaces

Polarized x-ray absorption fine structure (XAFS) spectroscopy has been performed in fluorescence mode under total external reflection conditions on frozen hydrated single monolayers of yeast cytochrome c (YCC). The protein molecules were vectorially oriented within the monolayer by tethering their naturally occurring and unique surface cysteine residues to the sulfhydryl-endgroups at the surface of a mixed organic self-assembled monolayer, itself covalently attached to an ultrapure silicon wafer. The sulfhydryl-endgroups were isolated by dilution with either methyl- or hydroxyl-endgroups, producing macroscopically nonpolar or uncharged-polar soft surfaces, respectively. Independent information on the heme-plane orientation relative to the monolayer plane was obtained experimentally via optical linear dichroism. The polarized XAFS data have been analyzed both qualitatively and by a global mapping approach limited to systematically altering the various iron-ligand distances within a model for the local atomic environment of the heme prosthetic group, and comparing the theoretically generated XAFS spectra with those obtained experimentally. A similar analysis of unpolarized XAFS data from a frozen solution of YCC was performed using either the heme environment from the NMR solution or the x-ray crystallographic data for YCC as the model structure. All resulting iron-ligand distances were then used in molecular dynamics (MD) computer simulations of YCC in these three systems to investigate the possible effects of anisotropic ligand motions on the fits of the calculated to the experimental XAFS spectra.     

The crystal structure of trypanothione reductase from the human pathogen Trypanosoma cruzi at 2.3 A resolution.

Trypanothione reductase (TR) is an NADPH-dependent flavoprotein unique to protozoan parasites from the genera Trypanosoma and Leishmania and is an important target for the design of improved trypanocidal drugs. We present details of the structure of TR from the human pathogen Trypanosoma cruzi, the agent responsible for Chagas' disease or South American trypanosomiasis. The structure has been solved by molecular replacement, using as the starting model the structure of the enzyme from the nonpathogenic Crithidia fasciculata, and refined to an R-factor of 18.9% for 53,868 reflections with F > or = sigma F between 8.0 and 2.3 A resolution. The model comprises two subunits (968 residues), two FAD prosthetic groups, two maleate ions, and 419 water molecules. The accuracy and geometry of the enzyme model is improved with respect to the C. fasciculata enzyme model. The new structure is described and specific features of the enzyme involved in substrate interactions are compared with previous models of TR and related glutathione reductases from human and Escherichia coli. Structural differences at the edge of the active sites suggest an explanation for the differing specificities toward glutathionylspermidine disulfide.