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.     

Two conformations of crystalline adenylate kinase.

Pig muscle adenylate kinase (EC2.7.4.3) can exist in three crystal forms, which are interconvertible. For crystal form A the enzyme structure is known in atomic detail. We report the X-ray diffraction analysis of crystal form B at 4.7 Å resolution and a comparison with the A form. During the transition from A to B the packing arrangement of the molecules changes slightly. Moreover, the individual molecule undergoes an appreciable conformational change: by displacing a chain segment of seven residues and two adjacent α-helices a hydrophobic pocket is opened deep in the cleft near the centre of the molecule. Concomitantly the β-pleated sheet is enlarged by about four hydrogen bonds in the B form. Several lines of evidence indicate that the observed conformational change is an intrinsic property of the molecule and is not induced by crystal packing forces.     

Simulating the minimum core for hydrophobic collapse in globular proteins.

To investigate the nature of hydrophobic collapse considered to be the driving force in protein folding, we have simulated aqueous solutions of two model hydrophobic solutes, methane and isobutylene. Using a novel methodology for determining contacts, we can precisely follow hydrophobic aggregation as it proceeds through three stages: dispersed, transition, and collapsed. Theoretical modeling of the cluster formation observed by simulation indicates that this aggregation is cooperative and that the simulations favor the formation of a single cluster midway through the transition stage. This defines a minimum solute hydrophobic core volume. We compare this with protein hydrophobic core volumes determined from solved crystal structures. Our analysis shows that the solute core volume roughly estimates the minimum core size required for independent hydrophobic stabilization of a protein and defines a limiting concentration of nonpolar residues that can cause hydrophobic collapse. These results suggest that the physical forces driving aggregation of hydrophobic molecules in water is indeed responsible for protein folding.     

Structural alteration of Escherichia coli Hsp31 by thermal unfolding increases chaperone activity

Escherichia coli Hsp31, encoded by hchA, is a heat-inducible molecular chaperone. We found that Hsp31 undergoes a conformational change via temperature-induced unfolding, generating a high molecular weight (HMW) form with enhanced chaperone activity. Although it has previously been reported that some subunits of the Hsp31 crystal structure show structural heterogeneity with increased hydrophobic surfaces, Hsp31 basically forms a dimer. We found that a C-terminal deletion (CΔ19) of Hsp31 exhibited structurally and functionally similar characteristics to that of the HMW form. Both the CΔ19 and HMW forms achieved a structure with considerably more β-sheets and less α-helices than the native dimeric form, exposing a portion of its hydrophobic surfaces. The structural alterations were determined from its spectral changes in circular dichroism, intrinsic fluorescence of tryptophan residues, and fluorescence of bis-ANS binding to a hydrophobic surface. Interestingly, during thermal transition, the dimeric Hsp31 undergoes a conformational change to the HMW species via the CΔ19 structure, as monitored with near-UV CD spectrum, implying that the CΔ19 resembles an intermediate state between the dimer and the HMW form. From these results, we propose that Hsp31 transforms itself into a fully functional chaperone by altering its tertiary and quaternary structures.     

The sequence TGAAKAVALVL from glyceraldehyde-3-phosphate dehydrogenase displays structural ambivalence and interconverts between alpha-helical and beta-hairpin conformations mediated by collapsed conformational states.

The peptide TGAAKAVALVL from glyceraldehyde-3-phosphate dehydrogenase adopts a helical conformation in the crystal structure and is a site for two hydrated helical segments, which are thought to be helical folding intermediates. Overlapping sequences of four to five residues from the peptide, sample both helical and strand conformations in known protein structures, which are dissimilar to glyceraldehyde-3-phosphate dehydrogenase suggesting that the peptide may have a structural ambivalence. Molecular dynamics simulations of the peptide sequence performed for a total simulation time of 1.2 micros, starting from the various initial conformations using GROMOS96 force field under NVT conditions, show that the peptide samples a large number of conformational forms with transitions from alpha-helix to beta-hairpin and vice versa. The peptide, therefore, displays a structural ambivalence. The mechanism from alpha-helix to beta-hairpin transition and vice versa reveals that the compact bends and turns conformational forms mediate such conformational transitions. These compact structures including helices and hairpins have similar hydrophobic radius of gyration (Rgh) values suggesting that similar hydrophobic interactions govern these conformational forms. The distribution of conformational energies is Gaussian with helix sampling lowest energy followed by the hairpins and coil. The lowest potential energy of the full helix may enable the peptide to take up helical conformation in the crystal structure of the glyceraldehyde-3-phosphate dehydrogenase, even though the peptide has a preference for hairpin too. The relevance of folding and unfolding events observed in our simulations to hydrophobic collapse model of protein folding are discussed.     

Molecular dynamics simulations of hydrophobic collapse of ubiquitin.

Nine nonnative conformations of ubiquitin, generated during two different thermal denaturation trajectories, were simulated under nearly native conditions (62 degrees C). The simulations included all protein and solvent atoms explicitly, and simulation times ranged from 1-2.4 ns. The starting structures had alpha-carbon root-mean-square deviations (RMSDs) from the crystal structure of 4-12 A and radii of gyration as high as 1.3 times that of the native state. In all but one case, the protein collapsed when the temperature was lowered and sampled conformations as compact as those reached in a control simulation beginning from the crystal structure. In contrast, the protein did not collapse when simulated in a 60% methanol:water mixture. The behavior of the protein depended on the starting structure: during simulation of the most native-like starting structures (<5 A RMSD to the crystal structure) the RMSD decreased, the number of native hydrogen bonds increased, and the secondary and tertiary structure increased. Intermediate starting structures (5-10 A RMSD) collapsed to the radius of gyration of the control simulation, hydrophobic residues were preferentially buried, and the protein acquired some native contacts. However, the protein did not refold. The least native starting structures (10-12 A RMSD) did not collapse as completely as the more native-like structures; instead, they experienced large fluctuations in radius of gyration and went through cycles of expansion and collapse, with improved burial of hydrophobic residues in successive collapsed states.     

On the Stability of Parainfluenza Virus 5 F Proteins

The crystal structure of the F protein (prefusion form) of the paramyxovirus parainfluenza virus 5 (PIV5) WR isolate was determined. We investigated the basis by which point mutations affect fusion in PIV5 isolates W3A and WR, which differ by two residues in the F ectodomain. The P22 stabilizing site acts through a local conformational change and a hydrophobic pocket interaction, whereas the S443 destabilizing site appears sensitive to both conformational effects and amino acid charge/polarity changes.     

Crystal structures of CheY from Thermotoga maritima do not support conventional explanations for the structural basis of enhanced thermostability.

The crystal structure of CheY protein from Thermotoga maritima has been determined in four crystal forms with and without Mg++ bound, at up to 1.9 A resolution. Structural comparisons with CheY from Escherichia coli shows substantial similarity in their folds, with some concerted changes propagating away from the active site that suggest how phosphorylated CheY, a signal transduction protein in bacterial chemotaxis, is recognized by its targets. A highly conserved segment of the protein (the "y-turn loop," residues 55-61), previously suggested to be a rigid recognition determinant, is for the first time seen in two alternative conformations in the different crystal structures. Although CheY from Thermotoga has much higher thermal stability than its mesophilic counterparts, comparison of structural features previously proposed to enhance thermostability such as hydrogen bonds, ion pairs, compactness, and hydrophobic surface burial would not suggest it to be so.     

Conserved water mediated recognition and the dynamics of active site Cys 331 and Tyr 411 in hydrated structure of human IMPDH‐II

Inosine monophosphate dehydrogenase (IMPDH) of human is involved in GMP biosynthesis pathway, increased level of IMPDH‐II (an isoform of enzyme) activity have found in leukemic and sarcoma cells. Modeling and extensive molecular dynamics simulation (15 ns) studies of IMPDH‐II (1B3O PDB structure) have indicated the intricate involvement of four conserved water molecules (W 1, W 2, W 3, and W 4) in the conformational transition or the mobilities of “flap” (residues 400–450) and “loop” (residues 325–342) regions in enzyme. The stabilization of active site residues Asn 303, Gly 324, Ser 329, Cys 331, Asp 364, and Tyr 411 through variable H‐bonding coordination from the conserved water molecular center seems interesting in the uninhibited hydrated form of human IMPDH‐II structures. This conformational transition or the flexibility of mobile regions, water molecular recognition to active site residues Cys 331 and Tyr 411, and the presence of a hydrophilic cavity ～540 ų (enclaved by the loop and flap region) near the C‐terminal surface of this enzyme may explore a rational hope toward the water mimic inhibitor or anticancer agent design for human. Copyright © 2010 John Wiley & Sons, Ltd.     

The active-inactive transition of human thymidylate synthase: targeted molecular dynamics simulations

Human thymidylate synthase (hTS) is an established anticancer target. It catalyses the production of 2'-deoxythymidine-5'-monophosphate, an essential building block for DNA synthesis. Because of the development of cellular drug resistance against current hTS inhibitors, alternative inhibition strategies are needed. hTS exists in two forms, active and inactive, defined by the conformation of the active-site (AS) loop, which carries the catalytic cysteine, C195. To investigate the mechanism of activation and inactivation, targeted molecular dynamics (TMD) simulations of the transitions between active and inactive states of hTS were performed. Analysis of changes in the dihedral angles in the AS loop during different TMD simulations revealed complex conformational transitions. Despite hTS being a homodimeric enzyme and the conformational transition significantly involving the dimer interface, the transition occurs in an asymmetric, sequential manner via an ensemble of pathways. In addition to C195, which reoriented during the simulations, other key residues in the rotation of the AS loop included W182 and R185. The interactions of the cognate bulky W182 residues at the dimer interface hindered the simultaneous twist of the AS loops in the hTS dimer. Interactions of R185, which is unique for hTS, with ligands at different allosteric sites affected the activation transition. In addition to providing insights into the activation/inactivation mechanism of hTS and how conformational transitions can occur in homodimeric proteins, our observations suggest that blocking of AS loop rotation by ligands binding in the large cavity between the loops could be one way to stabilize inactive hTS and inhibit the enzyme.     

The molecular basis of the temperature- and pH-induced conformational transitions in elastin-based peptides

Elastin undergoes an inverse temperature transition and collapses at high temperatures in both simulation and experiment. We investigated a pH-dependent modification of this transition by simulating a glutamic acid (Glu)-substituted elastin at varying pHs and temperatures. The Glu-substituted peptide collapsed at higher temperature than the unsubstituted elastin when Glu was charged. The charge effects could be reversed by neutralization of the Glu carboxyl groups at low pH, and in that case the peptide collapsed at a lower temperature. The collapse was accompanied by the formation of beta-turns and short distorted beta-sheets. Formation of contacts between hydrophobic side chains drives the collapse at high temperature, but interactions between water and polar groups (Glu and main chain) can attenuate this effect at high pH. The overall competition and balance of the polar and nonpolar groups determined the conformational states of the peptide. Water hydration contributed to the conformational transition, and the peptide and its hydration shell must be considered. Structurally, waters near polar residues mainly formed hydrogen bonds with the protein atoms, while waters around the hydrophobic side chains tended to be parallel to the peptide groups to maximize water-water interactions.     

The Conformational Transition Pathway of ATP Binding Cassette Transporter MsbA Revealed by Atomistic Simulations

ATP binding cassette transporters are integral membrane proteins that use the energy released from ATP hydrolysis at the two nucleotide binding domains (NBDs) to translocate a wide variety of substrates through a channel at the two transmembrane domains (TMDs) across the cell membranes. MsbA from Gram-negative bacteria is a lipid and multidrug resistance ATP binding cassette exporter that can undergo large scale conformational changes between the outward-facing and the inward-facing conformations revealed by crystal structures in different states. Here, we use targeted molecular dynamics simulation methods to explore the atomic details of the conformational transition from the outward-facing to the inward-facing states of MsbA. The molecular dynamics trajectories revealed a clear spatiotemporal order of the conformational movements. The disruption of the nucleotide binding sites at the NBD dimer interface is the very first event that initiates the following conformational changes, verifying the assumption that the conformational conversion is triggered by ATP hydrolysis. The conserved x-loops of the NBDs were identified to participate in the interaction network that stabilizes the cytoplasmic tetrahelix bundle of the TMDs and play an important role in mediating the cross-talk between the NBD and TMD. The movement of the NBD dimer is transmitted through x-loops to break the tetrahelix bundle, inducing the packing rearrangements of the transmembrane helices at the cytoplasmic side and the periplasmic side sequentially. The packing rearrangement within each periplasmic wing of TMD that results in exposure of the substrate binding sites occurred at the end stage of the trajectory, preventing the wrong timing of the binding site accessibility.     

Functional implications of structural differences between variants A and B of bovine beta-lactoglobulin

The structure of the trigonal crystal form of bovine beta-lactoglobulin variant B at pH 7.1 has been determined by X-ray diffraction methods at a resolution of 2.22 A and refined to values for R and Rfree of 0.239 and 0.286, respectively. By comparison with the structure of the trigonal crystal form of bovine beta-lactoglobulin variant A at pH 7.1, which was determined previously [Qin BY et al., 1998, Biochemistry 37:14014-14023], the structural consequences of the sequence differences D64G and V118A of variants A and B, respectively, have been investigated. Only minor differences in the core calyx structure occur. In the vicinity of the mutation site D64G on loop CD (residues 61-67), there are small changes in main-chain conformation, whereas the substitution V118A on beta-strand H is unaccompanied by changes in the surrounding structure, thereby creating a void volume and weakened hydrophobic interactions with a consequent loss of thermal stability relative to variant A. A conformational difference is found for the loop EF, implicated in the pH-dependent conformational change known as the Tanford transition, but it is not clear whether this reflects differences intrinsic to the variants in solution or differences in crystallization.     

Functional interactions at the interface between voltage-sensing and pore domains in the Shaker K(v) channel

Voltage-activated potassium (K(v)) channels contain a central pore domain that is partially surrounded by four voltage-sensing domains. Recent X-ray structures suggest that the two domains lack extensive protein-protein contacts within presumed transmembrane regions, but whether this is the case for functional channels embedded in lipid membranes remains to be tested. We investigated domain interactions in the Shaker K(v) channel by systematically mutating the pore domain and assessing tolerance by examining channel maturation, S4 gating charge movement, and channel opening. When mapped onto the X-ray structure of the K(v)1.2 channel the large number of permissive mutations support the notion of relatively independent domains, consistent with crystallographic studies. Inspection of the maps also identifies portions of the interface where residues are sensitive to mutation, an external cluster where mutations hinder voltage sensor activation, and an internal cluster where domain interactions between S4 and S5 helices from adjacent subunits appear crucial for the concerted opening transition.     

Integrated insights from simulation, experiment, and mutational analysis yield new details of LacI function

Protein structural change underlies many signal transduction processes. Although end-state structures are known for various allosteric proteins, intermediates are difficult to observe. Recently, targeted molecular dynamics simulation (TMD) was used to examine the conformational transition and predict relevant intermediates for wild-type lactose repressor (LacI). A catalog of involved residues suggests that the transition of this homodimer is asymmetric and that K84 is a prominent participant in the dynamic N-subdomain interface. Previous experiments indicated that hydrophobic substitutions at position 84 engender slowed, biphasic inducer binding kinetics, which might reflect the same phenomena observed in TMD. Here, we report biochemical confirmation that DNA and inducer binding remain allosterically linked in K84A and K84L, albeit with a differential smaller than that found in wild-type LacI. Other features of these mutant proteins are consistent with an allosteric conformational shift that approximates that of the wild type. As a consequence, these repressors can be utilized to explore an unanswered question about LacI function: How many inducers (one or two per dimer) are required to diminish operator affinity? The biphasic natures of the K84L and K84A inducer association rates allow direct correlation between the two distinct inducer binding events and operator release. Indeed, the kinetics of operator release for the K84A and K84L closely parallel those for the second inducer binding event. Together with implications from previous equilibrium results for wild-type and mutant proteins, these kinetic data demonstrate that binding of two inducers per dimeric DNA binding unit is required to release the operator in these variant LacI proteins.     

Conformational pathway for the kissing complex-->extended dimer transition of the SL1 stem-loop from genomic HIV-1 RNA as monitored by targeted molecular dynamics techniques

HIV-1 retroviral genomic RNA dimerization is initiated by loop-loop interactions between the SL1 stem-loops of two identical RNA molecules. The SL1-SL1 unstable resulting kissing complex (KC) then refolds irreversibly into a more stable complex called extended dimer (ED). Although the structures of both types of complex have been determined, very little is known about the conformational pathway corresponding to the transition, owing to the difficulty of observing experimentally intermediate conformations. In this study, we applied targeted molecular dynamics simulation techniques (TMD) to the phosphorus atoms for monitoring this pathway for the backbone, and a two-step strategy was adopted. In a first step, called TMD(-1), the dimer structure was constrained to progressively move away from KC without indicating the direction, until the RMSD from KC reaches 36A. A total of 20 TMD(-1) simulations were performed under different initial conditions and different simulation parameters. For RMSD ranging between 0 and 22A, the whole set of TMD(-1) simulations follows a similar pathway, then divergences are observed. None of the simulations leads to the ED structure. At RMSD=22A, the dimers look like two parallel Us, still linked by the initial loop-loop interaction, but the strands of the stems (the arms of the Us) are positioned in such a manner that they can form intramolecular as well as intermolecular Watson-Crick base-pairs. This family of structure is called UU. In a second step (TMD simulations), 18 structures were picked up along the pathways generated with TMD(-1) and were constrained to move toward ED by decreasing progressively their RMSD from ED. We found that only structures from the UU family are able to easily reach ED-like conformations of the backbones without exhibiting a large constraint energy.     

Crystal structure of maltooligosyltrehalose trehalohydrolase from Deinococcus radiodurans in complex with disaccharides

Trehalose (alpha-D-glucopyranosyl-1,1-alpha-D-glucopyranose) is a non-reducing diglucoside found in various organisms that serves as a carbohydrate reserve and as an agent that protects against a variety of physical and chemical stresses. Deinococcus radiodurans possesses an alternative biosynthesis pathway for the synthesis of trehalose from maltooligosaccharides. This reaction is mediated by two enzymes: maltooligosyltrehalose synthase (MTSase) and maltooligosyltrehalose trehalohydrolase (MTHase). Here, we present the 1.1A resolution crystal structure of MTHase. It consists of three major domains: two beta-sheet domains and a conserved glycosidase (beta/alpha)8 barrel catalytic domain. Three subdomains consisting of short insertions were identified within the catalytic domain. Subsequently, structures of MTHase in complex with maltose and trehalose were obtained at 1.2 A and 1.5 A resolution, respectively. These structures reveal the importance of the three inserted subdomains in providing the key residues required for substrate recognition. Trehalose is recognised specifically in the +1 and +2 binding subsites by an extensive hydrogen-bonding network and a strong hydrophobic stacking interaction in between two aromatic residues. Moreover, upon binding to maltose, which mimics the substrate sugar chain, a major concerted conformational change traps the sugar chain in the active site. The presence of magnesium in the active site of the MTHase-maltose complex suggests that MTHase activity may be regulated by divalent cations.     

Structural Insight into the Transmembrane Domain and the Juxtamembrane Region of the Erythropoietin Receptor in Micelles

Erythropoietin receptor (EpoR) dimerization is an important step in erythrocyte formation. Its transmembrane domain (TMD) and juxtamembrane (JM) region are essential for signal transduction across the membrane. A construct compassing residues S212–P259 and containing the TMD and JM region of the human EpoR was purified and reconstituted in detergent micelles. The solution structure of the construct was determined in dodecylphosphocholine (DPC) micelles by solution NMR spectroscopy. Structural and dynamic studies demonstrated that the TMD and JM region are an α-helix in DPC micelles, whereas residues S212–D224 at the N-terminus of the construct are not structured. The JM region is a helix that contains a hydrophobic patch formed by conserved hydrophobic residues (L253, I257, and W258). Nuclear Overhauser effect analysis, fluorescence spectroscopy, and paramagnetic relaxation enhancement experiments suggested that the JM region is exposed to the solvent. The structures of the TMD and JM region of the mouse EpoR were similar to those of the human EpoR.     

New insight into the pH-dependent conformational changes in bovine beta-lactoglobulin from Raman optical activity.

We have studied the conformation of beta-lactoglobulin in aqueous solution at room temperature over the pH range approximately 2.0-9.0 using vibrational Raman optical activity (ROA). The ROA spectra clearly show that the basic up and down beta-barrel core is preserved over the entire pH range, in agreement with other studies. However, from the shift of a sharp positive ROA band at approximately 1268 to approximately 1294 cm(-1) on going from pH values below that of the Tanford transition, which is centered at pH approximately 7.5, to values above, the Tanford transition appears to be associated with changes in the local conformations of residues in loop sequences possibly corresponding to a migration into the alpha-helical region of the Ramachandran surface from a nearby region. These changes may be related to those detected in X-ray crystal structures which revealed that the Tanford transition is associated with conformational changes in loops which form a doorway to the interior of the protein. The results illustrate how the ability of ROA to detect loop and turn structure separately from secondary structure is useful for studying conformational plasticity in proteins.     

Detailed conformational dynamics of juxtamembrane region and activation loop in c-Kit kinase activation process

The stem cell factor receptor (c-Kit) plays critical roles in initiating cell growth and proliferation. Its kinase functional abnormality has been thought to associate with several human cancers. The regulation of c-Kit kinase activity is achieved by phosphorylation on the residues Tyr568 and Tyr570 within juxtamembrane region (JMR) and subsequent structural transition of JMR and activation loop (A-loop). However, the detailed conformational dynamics of JMR and A-loop are far from clear, especially whether their conformational changes are coupled or not during the kinase activation transition. In this investigation, the complete conformational transition pathway was determined using a series of nanosecond conventional molecular dynamics (MD) and targeted molecular dynamics (TMD) simulations in explicit water systems. The results of the MD simulations show that the phosphorylation of residues Tyr568 and Tyr570 within JMR induces the detachment of JMR from the kinase C-lobe and increases the fluctuation in the structure of JMR, thus appearing to initiate the kinase activation process. During the course of the TMD simulation, which characterizes the conformational transition of c-Kit from autoinhibitory to activated state, the JMR undergoes a rapid departure from the allosteric binding site and drifts into solvent, followed by the conformational flip of A-loop from inactive (fold) state to active (extended) state. A change in the orientation of helix alphaC in response to the motion of JMR and A-loop has also been observed. The computational results presented here indicate that the dissociation of JMR from the kinase domain is prerequisite to c-Kit activation, which is consistent with previous experiments.