Efficient generation of feasible pathways for protein conformational transitions

We develop a computationally efficient method to simulate the transition of a protein between two conformations. Our method is based on a coarse-grained elastic network model in which distances between spatially proximal amino acids are interpolated between the values specified by the two end conformations. The computational speed of this method depends strongly on the choice of cutoff distance used to define interactions as measured by the density of entries of the constant linking/contact matrix. To circumvent this problem we introduce the concept of using a cutoff based on a maximum number of nearest neighbors. This generates linking matrices that are both sparse and uniform, hence allowing for efficient computations that are independent of the arbitrariness of cutoff distance choices. Simulation results demonstrate that the method developed here reliably generates feasible intermediate conformations, because our method observes steric constraints and produces monotonic changes in virtual bond and torsion angles. Applications are readily made to large proteins, and we demonstrate our method on lactate dehydrogenase, citrate synthase, and lactoferrin. We also illustrate how this framework can be used to complement experimental techniques that partially observe protein motions.     

Amaranthus caudatus lectin with polyproline II fold: conformational and functional transitions and molecular dynamics

Polyproline II (PPII) fold, a peculiar structural element was detected in the Amaranthus caudatus seed lectin (ACL) based on far UV circular dichroism spectrum, conformational transitions of the lectin, and a distinct isodichroic point in thermal denaturation. It was confirmed using PolyprOnline database to estimate the percentage of amino acids contributing to PPII fold and showed the values as 13.5 and 13.9% for PROSS and XTLSSTR, respectively. Investigations of the functional and conformational transitions of ACL during thermal-, pH-, and guanidine hydrochloride (GdnHCl)-induced denaturation were carried out using biochemical and biophysical techniques and molecular dynamics (MD) simulations approach. The lectin got aggregated at 60°C with instantaneous structural alterations. The aggregation-prone regions in ACL were predicted using online servers viz. AGGRESCAN, AmylPred, FoldAmyloid, and Waltz that were represented by Visual Molecular Dynamics tools. Nine conserved regions were identified by these softwares as being ‘hot-spots’ for aggregation. MD simulation studies of the lectin at 60°C revealed increase in radius of gyration. The loss of PPII fold in 2.0 M GdnHCl was reversible. The partially unfolded intermediate of ACL with diminished PPII fold formed at pH 1.0 was stable up to 90°C. The polyproline II fold has been rarely detected in lectins, ACL being the second after the potato lectin.     

Targeted molecular dynamics simulation studies of binding and conformational changes in E. coli MurD

Enzymes involved in the biosynthesis of bacterial peptidoglycan, an essential cell wall polymer unique to prokaryotic cells, represent a highly interesting target for antibacterial drug design. Structural studies of E. coli MurD, a three-domain ATP hydrolysis driven muramyl ligase revealed two inactive open conformations of the enzyme with a distinct C-terminal domain position. It was hypothesized that the rigid body rotation of this domain brings the enzyme to its closed active conformation, a structure, which was also determined experimentally. Targeted molecular dynamics 1 ns-length simulations were performed in order to examine the substrate binding process and gain insight into structural changes in the enzyme that occur during the conformational transitions into the active conformation. The key interactions essential for the conformational transitions and substrate binding were identified. The results of such studies provide an important step toward more powerful exploitation of experimental protein structures in structure-based inhibitor design.     

5 Nanosecond molecular dynamics and NMR study of conformational transitions in the sialyl-Lewis X antigen

The range of internal motions of the sialyl Lewis-X (SLe^x) tetraseccharide(NeuNAc     

Vasopressin conformational fluctuations: a molecular dynamics study

Molecular dynamics simulations have been used to study the conformational fluctuations of the oligopeptide hormone vasopressin. Starting coordinates for these simulations were built upon the crystal structure of pressinoic acid, the cyclic ring moiety of vasopressin, recently determined by x-ray diffraction. Coordinates for the additional tripeptide “tail” of vasopressin were selected by arbitrary positioning of this segment using interactive computer graphics. Two such starting configurations were minimized to relax strains, and long dynamics simulations (20 and 40 ps) in vacuo were then conducted following extensive heating and equilibration sequences (36 ps). In these studies, vasopressin was found to undergo few substantial conformational changes at 300 K on the time scale simulated, in contrast to the results of a shorter previous simulation, but comparable structural transitions were observed during the equilibration periods. The pressinoic acid structure was found to be a reasonably stable possible conformation for vasopressin in vacuum on this time scale.     

A soft-modeling approach to interpret thermodynamic and conformational transitions of polynucleotides

Multivariate outputs from the experimental monitoring of biochemical processes are usually difficult to interpret applying methods based on a priori chemical models. Curve resolution methods are model-free procedures, generally known as soft-modeling methods, which obtain the concentration profiles and instrumental responses of each individual species involved in a multivariate monitored process without making any kind of external assumption. Of the curve resolution methods available, the alternating least squares (ALS) is proposed here because of its ability to operate on one or on several matrices. Furthermore, ALS allows the introduction of information related to the internal data structure and to the general features of the concentration profiles and instrumental responses through the input of suitable constraints in the iterative resolution procedure. The ALS potential is tested on several data sets coming from the multivariate spectrometric monitoring of polyuridylic (polyU), polycytidylic (polyC), and polyadenylic (polyA) protonation equilibria in dioxane/water 30% (v/v). Information concerning the evolution of the concentration profiles and the spectra of each individual species involved in the acid-base equilibria, the presence and pattern of polyelectrolyte effects, and the presence of conformational transitions associated or not with the proton uptake process is presented.     

The conformational transition pathways of ATP-binding cassette transporter BtuCD revealed by targeted molecular dynamics simulation

BtuCD is a member of the ATP-binding cassette transporters in Escherichia coli that imports vitamin B(12) into the cell by utilizing the energy of ATP hydrolysis. Crystal structures of BtuCD and its homologous protein HI1470/1 in various conformational states support the "alternating access" mechanism which proposes the conformational transitions of the substrate translocation pathway at transmembrane domain (TMD) between the outward-facing and inward-facing states. The conformational transition at TMD is assumed to couple with the movement of the cytoplasmic nucleotide-binding domains (NBDs) driven by ATP hydrolysis/binding. In this study, we performed targeted molecular dynamics (MD) simulations to explore the atomic details of the conformational transitions of BtuCD importer. The outward-facing to inward-facing (O→I) transition was found to be initiated by the conformational movement of NBDs. The subsequent reorientation of the substrate translocation pathway at TMD began with the closing of the periplasmic gate, followed by the opening of the cytoplamic gate in the last stage of the conformational transition due to the extensive hydrophobic interactions at this region, consistent with the functional requirement of unidirectional transport of the substrates. The reverse inward-facing to outward-facing (I→O) transition was found to exhibit intrinsic diversity of the conformational transition pathways and significant structural asymmetry, suggesting that the asymmetric crystal structure of BtuCD-F is an intermediate state in this process.     

Targeted molecular dynamics simulation studies of calcium binding and conformational change in the C-terminal half of gelsolin

Gelsolin consists of six related domains (G1-G6) and the C-terminal half (G4-G6) acts as a calcium sensor during the activation of the whole molecule, a process that involves large domain movements. In this study, we used targeted molecular dynamics simulations to elucidate the conformational transitions of G4-G6 at an atomic level. Domains G4 and G6 are initially ruptured, followed by a rotation of G6 by approximately 90 degrees , which is the dominant conformational change. During this period, local conformational changes occur at the G4 and G5 calcium-binding sites, facilitating large changes in interdomain distances. Alterations in the binding affinities of the calcium ions in these three domains appear to be related to local conformational changes at their binding sites. Analysis of the relative stabilities of the G4-G6-bound calcium ions suggests that they bind first to G6, then to G4, and finally to G5.     

Targeted gene correction: a new strategy for molecular medicine

Advances, over the past 20 years, in the genetic manipulation of mammalian cells form the scientific basis of gene therapy. A number of strategies are presently being used to replace or augment a dysfunctional gene with a correct copy of itself. Now, a novel approach to correct the dysfunctional gene in the chromosome is being developed. Data obtained from biochemical, cell-based and animal studies suggest that the era of gene repair is dawning. It is now conceivable that inherited and non-inherited disorders might be treated with a small molecular tool designed to fix the mutation directly. Here, the conceptualization of the technique and its barriers to success are discussed.     

Moving beyond biotic homogenization: searching for new insights into vegetation dynamics

Biotic homogenization has been predicted to occur in cities across the world. However, the empirical evidence has been less than convincing. Lososová et al. explore the middle ground between these two points of view in this issue of Journal of Vegetation Science. They take a more sophisticated approach, linking homogenization to bigger questions of vegetation assembly in urban environments.     

Experimental validation of molecular dynamics simulations of lipid bilayers: a new approach

A novel protocol has been developed for comparing the structural properties of lipid bilayers determined by simulation with those determined by diffraction experiments, which makes it possible to test critically the ability of molecular dynamics simulations to reproduce experimental data. This model-independent method consists of analyzing data from molecular dynamics bilayer simulations in the same way as experimental data by determining the structure factors of the system and, via Fourier reconstruction, the overall transbilayer scattering-density profiles. Multi-nanosecond molecular dynamics simulations of a dioleoylphosphatidylcholine bilayer at 66% RH (5.4 waters/lipid) were performed in the constant pressure and temperature ensemble using the united-atom GROMACS and the all-atom CHARMM22/27 force fields with the GROMACS and NAMD software packages, respectively. The quality of the simulated bilayer structures was evaluated by comparing simulation with experimental results for bilayer thickness, area/lipid, individual molecular-component distributions, continuous and discrete structure factors, and overall scattering-density profiles. Neither the GROMACS nor the CHARMM22/27 simulations reproduced experimental data within experimental error. The widths of the simulated terminal methyl distributions showed a particularly strong disagreement with the experimentally observed distributions. A comparison of the older CHARMM22 with the newer CHARMM27 force fields shows that significant progress is being made in the development of atomic force fields for describing lipid bilayer systems empirically.     

pH‐dependent conformational dynamics of beta‐secretase 1: A molecular dynamics study

Beta‐secretase 1 (BACE‐1) is an aspartyl protease implicated in the overproduction of β‐amyloid fibrils responsible for Alzheimer disease. The process of β‐amyloid genesis is known to be pH dependent, with an activity peak between solution pH of 3.5 and 5.5. We have studied the pH‐dependent dynamics of BACE‐1 to better understand the pH dependent mechanism. We have implemented support for graphics processor unit (GPU) accelerated constant pH molecular dynamics within the AMBER molecular dynamics software package and employed this to determine the relative population of different aspartyl dyad protonation states in the pH range of greatest β‐amyloid production, followed by conventional molecular dynamics to explore the differences among the various aspartyl dyad protonation states. We observed a difference in dynamics between double‐protonated, mono‐protonated, and double‐deprotonated states over the known pH range of higher activity. These differences include Tyr 71‐aspartyl dyad proximity and active water lifetime. This work indicates that Tyr 71 stabilizes catalytic water in the aspartyl dyad active site, enabling BACE‐1 activity.     

Investigation of the conformational dynamics of the A2A adenosine receptor by molecular dynamics simulation

The structural behavior of the ligand-free form of adenosine receptor A_2A in an explicit membrane-mimicking environment was investigated by molecular dynamics (MD) simulation. Principal components analysis was applied to the series of MD snapshots and to a collection of X-ray structures of the A_2A receptor. The resulting charts revealed a correlation in the dynamic behavior of the receptor observed in the MD trajectories and in the experimental dataset. The most pronounced structural dynamics in the A_2A receptor were observed in the intracellular part: TM 5 and 6 with the connecting loop, just as generally recognized in crystallographic studies and attributed to receptor activation. There are grounds for supposing that this pattern of intramolecular motions ensues directly from the spatial architecture (fold) of the A_2A receptor.     

Iterative cluster-NMA: A tool for generating conformational transitions in proteins

Computational models provide insight into the structure-function relationship in proteins. These approaches, especially those based on normal mode analysis, can identify the accessible motion space around a given equilibrium structure. The large magnitude, collective motions identified by these methods are often well aligned with the general direction of the expected conformational transitions. However, these motions cannot realistically be extrapolated beyond the local neighborhood of the starting conformation. In this article, the iterative cluster-NMA (icNMA) method is presented for traversing the energy landscape from a starting conformation to a desired goal conformation. This is accomplished by allowing the evolving geometry of the intermediate structures to define the local accessible motion space, and thus produce an appropriate displacement. Following the derivation of the icNMA method, a set of sample simulations are performed to probe the robustness of the model. A detailed analysis of beta1,4-galactosyltransferase-T1 is also given, to highlight many of the capabilities of icNMA. Remarkably, during the transition, a helix is seen to be extended by an additional turn, emphasizing a new unknown role for secondary structures to absorb slack during transitions. The transition pathway for adenylate kinase, which has been frequently studied in the literature, is also discussed.     

Photoinduced conformational dynamics of a photoswitchable peptide: a nonequilibrium molecular dynamics simulation study

Employing nonequilibrium molecular dynamics simulations, a comprehensive computational study of the photoinduced conformational dynamics of a photoswitchable bicyclic azobenzene octapeptide is presented. The calculation of time-dependent probability distributions along various global and local reaction coordinates reveals that the conformational rearrangement of the peptide is rather complex and occurs on at least four timescales: 1) After photoexcitation, the azobenzene unit of the molecule undergoes nonadiabatic photoisomerization within 0.2 ps. 2) On the picosecond timescale, the cooling (13 ps) and the stretching (14 ps) of the photoexcited peptide is observed. 3) Most reaction coordinates exhibit a 50-100 ps component reflecting a fast conformational rearrangement. 4) The 500-1000 ps component observed in the simulation accounts for the slow diffusion-controlled conformational equilibration of the system. The simulation of the photoinduced molecular processes is in remarkable agreement with time-resolved optical and infrared experiments, although the calculated cooling as well as the initial conformational rearrangements of the peptide appear to be somewhat too slow. Based on an ab initio parameterized vibrational Hamiltonian, the time-dependent amide I frequency shift is calculated. Both intramolecular and solvent-induced contributions to the frequency shift were found to change by < or = 2 cm(-1), in reasonable agreement with experiment. The potential of transient infrared spectra to characterize the conformational dynamics of peptides is discussed in some detail.     

Molecular dynamics simulations of forced conformational transitions in 1,6-linked polysaccharides

Recent atomic force microscopy stretching measurements of single polysaccharide molecules suggest that their elasticity is governed by force-induced conformational transitions of the pyranose ring. However, the mechanism of these transitions and the mechanics of the pyranose ring are not fully understood. Here we use steered molecular dynamics simulations of the stretching process to unravel the mechanism of forced conformational transitions in 1,6 linked polysaccharides. In contrast to most sugars, 1,6 linked polysaccharides have an extra bond in their inter-residue linkage, C5-C6, around which restricted rotations occur and this additional degree of freedom increases the mechanical complexity of these polymers. By comparing the computational results with the atomic force microscopy data we determine that forced rotations around the C5-C6 bond have a significant and different impact on the elasticity of alpha- and beta-linked polysaccharides. Beta-linkages of a polysaccharide pustulan force the rotation around the C5-C6 bonds and produce a Hookean-like elasticity but do not affect the conformation of the pyranose rings. However, alpha-linkages of dextran induce compound conformational transitions that include simultaneous rotations around the C5-C6 bonds and chair-boat transitions of the pyranose rings. These previously not-recognized transitions are responsible for the characteristic plateau in the force-extension relationship of dextran.     

A new method to monitor the rate of conformational transitions in RNA.

Many RNAs need Mg2+to produce stable tertiary structures. Here we describe a simple method to measure the rate and activation parameters of tertiary structure unfolding that exploits this Mg2+dependence. Our approach is based on mixing an RNA solution with excess EDTA in a stopped-flow instrument equipped with an absorbance detector, under conditions of temperature and ionic strength where, after chelation of Mg2+, tertiary structure unfolds. We have demonstrated the utility of this method by studying phenylalanine-specific transfer RNA from yeast (tRNAPhe) because the unfolding rates and the corresponding activation parameters have been determined previously and provide a benchmark for our technique. We find that within error, our stopped-flow method reproduces both the rate and activation enthalpy for tertiary unfolding of yeast tRNAPhe measured previously by temperature-jump relaxation kinetics. Since many different RNAs require divalent magnesium for tertiary structure stabilization, this technique should be applicable to study the folding of other RNAs.     

A generalized conformational energy function of DNA derived from molecular dynamics simulations

Proteins recognize DNA sequences by two different mechanisms. The first is direct readout, in which recognition is mediated by direct interactions between the protein and the DNA bases. The second is indirect readout, which is caused by the dependence of conformation and the deformability of the DNA structure on the sequence. Various energy functions have been proposed to evaluate the contribution of indirect readout to the free-energy changes in complex formations. We developed a new generalized energy function to estimate the dependence of the deformability of DNA on the sequence. This function was derived from molecular dynamics simulations previously conducted on B-DNA dodecamers, each of which had one possible tetramer sequence embedded at its center. By taking the logarithm of the probability distribution function (PDF) for the base-step parameters of the central base-pair step of the tetramer, its ability to distinguish the native sequence from random ones was superior to that with the previous method that approximated the energy function in harmonic form. From a comparison of the energy profiles calculated with these two methods, we found that the harmonic approximation caused significant errors in the conformational energies of the tetramers that adopted multiple stable conformations.