The pH-dependent conformational states of kyotorphin: a constant-pH molecular dynamics study

An extensive conformational study of the analgesic dipeptide kyotorphin (L-Tyr-L-Arg) at different pH values was performed using a constant-pH molecular dynamics method. This dipeptide showed a remarkable pH-dependent conformational variety. The protonation of the N-terminal amine was identified as a key element in the transition between the more extended and the more packed conformational states, as monitored by the dihedral angle defined by the atoms 1Cbeta-1Calpha-2Calpha-2Cbeta. The principal-component analysis of kyotorphin identified two major conformational populations (the extended trans and the packed cis) together with conformations that occur exclusively at extreme pH values. Other, less stable conformations were also identified, which help us to understand the transitions between the predominant populations. The fitting of kyotorphin's conformational space to the structure of morphine resulted in a set of conformers that were able to fulfill most of the constraints for the mu-receptor. These results suggest that there may be strong similarities between the kyotorphin receptor and the structural family of opioid receptors.     

C-terminal actin-binding sites of smooth muscle caldesmon switch actin between conformational states

Caldesmon is a component of the thin filaments of smooth muscles where it is believed to play an essential role in regulating the thin filaments’ interaction with myosin and hence contractility. We studied the effects of caldesmon and two recombinant fragments CaDH1 (residues 506–793) and CaDH2 (residues 683–767) on the structure of actin–tropomyosin by making measurements of the fluorescence polarisation of probes specifically attached to actin. CaDH1, like the parent molecule caldesmon, is an inhibitor of actin–tropomyosin interaction with myosin whilst CaDH2 is an activator. The F-actin in permeabilised and myosin free rabbit skeletal muscle ‘ghost’ fibres was labelled by tetramethyl rhodamine-isothiocyanate (TRITC)–phalloidin or fluorescein-5′-isothiocyanate (FITC) at lysine 61. Fluorescence polarisation measurements were made and the parameters Φ_A, Φ_E, Θ_1/2 and N were calculated. Φ_A and Φ_E are angles between the fiber axis and the absorption and emission dipoles, respectively; Θ_1/2 is the angle between the F-actin filament axis and the fiber axis; N is the relative number of randomly oriented fluorophores. Actin–tropomyosin interaction with myosin subfragment-1 induced changes in the parameters of the polarised fluorescence that are typical of strong binding of myosin to actin and of the ‘on’ conformational state of actin. Caldesmon and CaDH1 (as well as troponin in the absence of Ca²⁺) diminished the effect of S-1, whereas CaDH2 (as well as troponin in the presence of Ca²⁺) enhanced the effect of S1. Thus the structural evidence correlates with biochemical evidence that C-terminal actin-binding sites of caldesmon can modulate the structural transition of actin monomers between ‘off’ (caldesmon and CaDH1) and ‘on’ (S-1 and CaDH2) states in a manner analogous to troponin.     

Mechanical response and conformational changes of alpha-actinin domains during unfolding: a molecular dynamics study

Alpha-actinin is a cytoskeleton-binding protein involved in the assembly and regulation of the actin filaments. In this work molecular dynamics method was applied to investigate the mechanical behaviour of the human skeletal muscle α-actinin. Five configurations were unfolded at an elongation speed of 0.1 nm/ps in order to investigate the conformational changes occurring during the extension process. Moreover, a sensitivity analysis at different velocities was performed for one of the R2–R3 spectrin-like repeat configuration extracted in order to evaluate the effect of the pulling speed on the mechanical behaviour of the molecule. Two different behaviours were recognized with respect to the pulling speed. In particular, at speed higher than 0.025 nm/ps a continuous rearrangement without evident force peaks was obtained, on the contrary at lower speed evident peaks in the range 500–750 pN were detected. R3 repeat resulted more stable than R2 during mechanical unfolding, due to the lower hydrophobic surface available to the solvent. The characterization of the R2–R3 units can be useful for the development of cytoskeleton network models based on stiffness values obtained by analyses performed at the molecular level.     

Lipid domains of mycobacteria studied with fluorescent molecular probes

The complex mycobacterial cell envelope is recognized as a critical factor in our failure to control tuberculosis, leprosy and other non-tuberculous pathogens. Although its composition has been extensively determined, many details regarding the organization of the envelope remain uncertain. This is particularly so for the non-covalently bound lipids, whose natural distribution may be disrupted by conventional biochemical or cytological techniques. In order to study the native organization of lipid domains in the mycobacterial envelope, we have applied a range of fluorescent lipophilic probes to live mycobacteria, including Mycobacterium smegmatis, Mycobacterium tuberculosis, Mycobacterium avium, Mycobacterium gadium and Mycobacterium aurum, and analysed the resultant signals by fluorescence microscopy and digital image processing. Five key features were observed: (i) the presence of both envelope and intracellular lipid domains; (ii) differential localization of probes into these domains influenced predominantly by their hydrophobicity, as modelled by their calculated octanol:water partition coefficients and by their amphiphilicities; (iii) uneven distribution of lipophilic material in the envelope; (iv) selective labelling of septal regions of the envelope; and (v) modification of labelling patterns by additional treatments such as fluorescence quenching antibodies, detergents and solvents. Using this last approach, a coherent cell envelope lipid domain was demonstrated outside the cytoplasmic membrane and, for the first time, the proposed covalently linked mycolyl-arabinogalactan-peptidoglycan macromolecular complex was imaged directly. The use of fluorescent probes and high-resolution fluorescence microscopy has enabled us to obtain a coherent view of distinct lipid domains in mycobacteria. Further application of this approach will facilitate understanding of the role of lipids in the physiology of these organisms.     

Nucleotide-mediated conformational changes of monomeric actin and Arp3 studied by molecular dynamics simulations

Members of the actin family of proteins exhibit different biochemical properties when ATP, ADP-P_i, ADP, or no nucleotide is bound. We used molecular dynamics simulations to study the effect of nucleotides on the behavior of actin and actin-related protein 3 (Arp3). In all of the actin simulations, the nucleotide cleft stayed closed, as in most crystal structures. ADP was much more mobile within the cleft than ATP, despite the fact that both nucleotides adopt identical conformations in actin crystal structures. The nucleotide cleft of Arp3 opened in most simulations with ATP, ADP, and no bound nucleotide. Deletion of a C-terminal region of Arp3 that extends beyond the conserved actin sequence reduced the tendency of the Arp3 cleft to open. When the Arp3 cleft opened, we observed multiple instances of partial release of the nucleotide. Cleft opening in Arp3 also allowed us to observe correlated movements of the phosphate clamp, cleft mouth, and barbed-end groove, providing a way for changes in the nucleotide state to be relayed to other parts of Arp3. The DNase binding loop of actin was highly flexible regardless of the nucleotide state. The conformation of Ser14/Thr14 in the P1 loop was sensitive to the presence of the γ-phosphate, but other changes observed in crystal structures were not correlated with the nucleotide state on nanosecond timescales. The divalent cation occupied three positions in the nucleotide cleft, one of which was not previously observed in actin or Arp2/3 complex structures. In sum, these simulations show that subtle differences in structures of actin family proteins have profound effects on their nucleotide-driven behavior.     

Influence of a fluorobenzene nucleobase analogue on the conformational flexibility of RNA studied by molecular dynamics simulations

Chemically modified bases are frequently used to stabilize nucleic acids, to study the driving forces for nucleic acid structure formation and to tune DNA and RNA hybridization conditions. In particular, fluorobenzene and fluorobenzimidazole base analogues can act as universal bases able to pair with any natural base and to stabilize RNA duplex formation. Although these base analogues are compatible with an A-form RNA geometry, little is known about the influence on the fine structure and conformational dynamics of RNA. In the present study, nano-second molecular dynamics (MD) simulations have been performed to characterize the dynamics of RNA duplexes containing a central 1'-deoxy-1'-(2,4-difluorophenyl)-beta-D-ribofuranose base pair or opposite to an adenine base. For comparison, RNA with a central uridine:adenine pair and a 1'-deoxy-1'-(phenyl)-beta-D-ribofuranose opposite to an adenine was also investigated. The MD simulations indicate a stable overall A-form geometry for the RNAs with base analogues. However, the presence of the base analogues caused a locally enhanced mobility of the central bases inducing mainly base pair shear and opening motions. No stable 'base-paired' geometry was found for the base analogue pair or the base analogue:adenine pairs, which explains in part the universal base character of these analogues. Instead, the conformational fluctuations of the base analogues lead to an enhanced accessibility of the bases in the major and minor grooves of the helix compared with a regular base pair.     

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.     

Calnexin discriminates between protein conformational states and functions as a molecular chaperone in vitro.

Although calnexin is thought to function as a molecular chaperone for glycoproteins, a prevalent view is that it cannot distinguish between protein conformational states, binding solely through its lectin site to monoglucosylated oligosaccharides. Using purified components in vitro, calnexin effectively prevented the aggregation not only of glycoproteins bearing monoglucosylated oligosaccharides but also proteins lacking N-glycans, an effect enhanced by ATP. It also suppressed the thermal denaturation of nonglycosylated proteins and enhanced their refolding in conjunction with other cellular components. Calnexin formed stable complexes with unfolded conformers of these proteins but not with the native molecules. Therefore, in addition to being a lectin, calnexin functions as a bona fide molecular chaperone capable of interacting with polypeptide segments of folding glycoproteins.     

Effects of conformational activation of integrin alpha 1I and alpha 2I domains on selective recognition of laminin and collagen subtypes

Collagen receptor integrins alpha 1 beta 1 and alpha 2 beta 1 can selectively recognize different collagen subtypes. Here we show that their alpha I domains can discriminate between laminin isoforms as well: alpha 1I and alpha 2I recognized laminin-111, -211 and -511, whereas their binding to laminin-411 was negligible. Residue Arg-218 in alpha1 was found to be instrumental in high-avidity binding. The gain-of-function mutation E318W makes the alpha 2I domain to adopt the "open" high-affinity conformation, while the wild-type alpha 2I domain favors the "closed" low-affinity conformation. The E318W mutation markedly increased alpha 2I domain binding to the laminins (-111, -211 and -511), leading us to propose that the activation state of the alpha 2 beta 1 integrin defines its role as a laminin receptor. However, neither wild-type nor alpha 2IE318W domain could bind to laminin-411. alpha 2IE318W also bound tighter to all collagens than alpha 2I wild-type, but it showed reduced ability to discriminate between collagens I, IV and IX. The corresponding mutation, E317A, in the alpha 1I domain transformed the domain into a high-avidity binder of collagens I and IV. Thus, our results indicate that conformational activation of integrin alpha 1I and alpha 2I domains leads to high-avidity binding to otherwise disfavored collagen subtypes.     

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.     

Multiple conformational states and gate opening of outer membrane protein TolC revealed by molecular dynamics simulations

Outer membrane protein TolC serves as an exit duct for exporting substances out of cell. The occluded periplasmic entrance of TolC is required to open for substrate transport, although the opening mechanism remains elusive. In this study, systematic molecular dynamics (MD) simulations for wild type TolC and six mutants were performed to explore the conformational dynamics of TolC. The periplasmic gate was shown to sample multiple conformational states with various degrees of gating opening. The gate opening was facilitated by all mutations except Y362F, which adopts an even more closed state than wild type TolC. The interprotomer salt‐bridge R367–D153 is turned out to be crucial for periplasmic gate opening. The mutations that disrupt the interactions at the periplasmic tip may affect the stability of the trimeric assembly of TolC. Structural asymmetry of the periplasmic gate was observed to be opening size dependent. Asymmetric conformations are found in moderately opening states, while the most and the least opening states are often more symmetric. Finally, it is shown that lowering pH can remarkably stabilize the closed state of the periplasmic gate. Proteins 2014; 82:2169–2179. © 2014 Wiley Periodicals, Inc.     

The structure of a membrane-spanning polypeptide studied by molecular dynamics

We have performed a molecular dynamics simulation of a 46-residue segment of glycophorin which includes the hydrophobic membrane-spanning region of this protein. The presence of a membrane and of water is taken into account in a continuum approximation which makes use of phenomenological hydrophobic energies. The initial -helical conformation and the membrane incorporation of the hydrophobic segment remain stable for the length of the simulation which is 100 ps. Moreover, when the hydrophobic segment is partially shifted out of the membrane, it moves back into the membrane. Superimposed on these deterministic effects one also observes thermal fluctuations in the form of bending and tilting of the membrane-spanning helix.     

The STATIS method: Characterization of conformational states of flexible molecules from molecular dynamics simulations in solution

STATIS, a data analysis method used when data can be expressed as matrices, seems particularly well suited to characterize the internal molecular motions and conformational states extracted from MD trajectories. We first outline this method and the “adapted STATIS” method. Applications are presented for 18-crown-6 (simulated for 2 nsec in acetonitrile solution) and for the (L30)₂Cu⁺ catenate (stimulated for 150 psec in chloroform). STATIS should be valuable for the classification of molecular conformations and simplified visualization of MD trajectories.     

The conformational behaviour of complexes of alpha-cyclodextrin with p-chlorophenol and p-hydroxybenzoic acid in water as studied by molecular dynamics simulations.

Molecular dynamics simulations were performed to obtain information about the conformational behaviour and stabilization of alpha-cyclodextrin (alpha CD) complexes in water. Simulations of p-chlorophenol and p-hydroxybenzoic acid in alpha CD showed that the complex is a very flexible system. The guest compound rotates inside the cavity and partly moves in and out. alpha CD continuously adapts its conformation to the orientation of the guest compound (or vice versa): the hexagon of the glycosidic oxygen atoms is stretched parallel to the plane of the aromatic ring of the guest compound during 80% of the simulation. This suggests that Van der Waals interactions play an important role in the stabilization of the complex. Each intramolecular hydrogen bond between neighbouring glucose units is formed during 30-80% of the simulation. Hydrogen bonds between alpha CD and the guest compound, on the other hand, are rarely formed. Thus, intermolecular hydrogen bonding seems to play a minor role in the stabilization of alpha CD complexes.     

The conformational dynamics of a metastable serpin studied by hydrogen exchange and mass spectrometry

Serpins are a class of protease inhibitors that initially fold to a metastable structure and subsequently undergo a large conformational change to a stable structure when they inhibit their target proteases. How serpins are able to achieve this remarkable conformational rearrangement is still not understood. To address the question of how the dynamic properties of the metastable form may facilitate the conformational change, hydrogen/deuterium exchange and mass spectrometry were employed to probe the conformational dynamics of the serpin human alpha(1)-antitrypsin (alpha(1)AT). It was found that the F helix, which in the crystal structure appears to physically block the conformational change, is highly dynamic in the metastable form. In particular, the C-terminal half of the F helix appears to spend a substantial fraction of time in a partially unfolded state. In contrast, beta-strands 3A and 5A, which must separate to accommodate insertion of the reactive center loop (RCL), are not conformationally flexible in the metastable state but are rigid and stable. The conformational lability required for loop insertion must therefore be triggered during the inhibition reaction. Beta-strand 1C, which anchors the distal end of the RCL and thus prevents transition to the so-called latent form, is also stable, consistent with the observation that alpha(1)AT does not spontaneously adopt the latent form. A surprising degree of flexibility is seen in beta-strand 6A, and it is speculated that this flexibility may deter the formation of edge-edge polymers.     

Unfolding of the cold shock protein studied with biased molecular dynamics

The cold shock protein from Bacillus caldolyticus is a small beta-barrel protein that folds in a two-state mechanism. For the native protein and for several mutants, a wealth of experimental data are available on stability and folding, so that it is an optimal system to study this process. We compare data from unfolding simulations (trajectories of 5 and up to 12 ns) obtained with a bias potential at room temperature and from unbiased thermal unfolding simulations with experimental data. The unfolding patterns derived from the trajectories starting from different native-like conformations and subject to different unfolding conditions agree. The transition state found in the simulations of unfolding is close to the native structure in agreement with experiment. Moreover, a lower value of the free energy barrier of unfolding was found for the mutant R3E than for the mutant E46A and the native protein, as indicated by experimental data. The first unfolding event involves the three-stranded beta-sheet whose decomposition corresponds to the transition state. In contrast to conclusions drawn from experiments, we found that the two-stranded beta-strand forms the most stable substructure, which decomposes very late in the unfolding process. However, assuming that this structure forms very early in the folding process, our findings would not contradict the experiments but require a different interpretation of them.     

A FRET-based sensor reveals large ATP hydrolysis-induced conformational changes and three distinct states of the molecular motor myosin

The molecular motor myosin is proposed to bind to actin and swing its light-chain binding region through a large angle to produce an approximately 10 nm step in motion coupled to changes in the nucleotide state at the active site. To date, however, direct dynamic measurements have largely failed to show changes of that magnitude. Here, we use a cysteine engineering approach to create a high resolution, FRET-based sensor that reports a large, approximately 70 degree nucleotide-dependent angle change of the light-chain binding region. The combination of steady-state and time-resolved fluorescence resonance energy transfer measurements unexpectedly reveals two distinct prestroke states. The measurements also show that bound Mg.ADP.Pi, and not bound Mg.ATP, induces the myosin to adopt the prestroke states.