Aspects of model building applied to the C-terminal domain of the L12 protein from chloroplast ribosomes: a molecular dynamics study

A 170 picosecond molecular dynamics trajectory has been calculated starting from a model-built structure of chloroplast CTF. Local conformational changes occur during the equilibration period. Thereafter, a dynamically stable structure is attained. The conformational changes involve a turn connecting two structural subdomains which has an amino acid insertion and several substitutions with respect to the E. coli sequence. Potential energy minimisation alone fails to detect such a change. The overall folding and atomic positional fluctuations are very similar to those found in MD simulations of the E. coli molecule. The combined use of computer graphics based model building and MD calculations has lead to a thermally stable putative structure for the chloroplast CTF.     

Locally accessible conformations of proteins: multiple molecular dynamics simulations of crambin.

Multiple molecular dynamics (MD) simulations of crambin with different initial atomic velocities are used to sample conformations in the vicinity of the native structure. Individual trajectories of length up to 5 ns sample only a fraction of the conformational distribution generated by ten independent 120 ps trajectories at 300 K. The backbone atom conformational space distribution is analyzed using principal components analysis (PCA). Four different major conformational regions are found. In general, a trajectory samples only one region and few transitions between the regions are observed. Consequently, the averages of structural and dynamic properties over the ten trajectories differ significantly from those obtained from individual trajectories. The nature of the conformational sampling has important consequences for the utilization of MD simulations for a wide range of problems, such as comparisons with X-ray or NMR data. The overall average structure is significantly closer to the X-ray structure than any of the individual trajectory average structures. The high frequency (less than 10 ps) atomic fluctuations from the ten trajectories tend to be similar, but the lower frequency (100 ps) motions are different. To improve conformational sampling in molecular dynamics simulations of proteins, as in nucleic acids, multiple trajectories with different initial conditions should be used rather than a single long trajectory.     

A comparison of the dynamic behavior of monomeric and dimeric insulin shows structural rearrangements in the active monomer

Molecular dynamics (MD) simulations (5-10ns in length) and normal mode analyses were performed for the monomer and dimer of native porcine insulin in aqueous solution; both starting structures were obtained from an insulin hexamer. Several simulations were done to confirm that the results obtained are meaningful. The insulin dimer is very stable during the simulation and remains very close to the starting X-ray structure; the RMS fluctuations calculated from the MD simulation agree with the experimental B-factors. Correlated motions were found within each of the two monomers; they can be explained by persistent non-bonded interactions and disulfide bridges. The correlated motions between residues B24 and B26 of the two monomers are due to non-bonded interactions between the side-chains and backbone atoms. For the isolated monomer in solution, the A chain and the helix of the B chain are found to be stable during 5ns and 10ns MD simulations. However, the N-terminal and the C-terminal parts of the B chain are very flexible. The C-terminal part of the B chain moves away from the X-ray conformation after 0.5-2.5ns and exposes the N-terminal residues of the A chain that are thought to be important for the binding of insulin to its receptor. Our results thus support the hypothesis that, when monomeric insulin is released from the hexamer (or the dimer in our study), the C-terminal end of the monomer (residues B25-B30) is rearranged to allow binding to the insulin receptor. The greater flexibility of the C-terminal part of the beta chain in the B24 (Phe-->Gly) mutant is in accord with the NMR results. The details of the backbone and side-chain motions are presented. The transition between the starting conformation and the more dynamic structure of the monomers is characterized by displacements of the backbone of Phe B25 and Tyr B26; of these, Phe B25 has been implicated in insulin activation.     

Inherent protein structural flexibility at the RNA-binding interface of L30e

The Saccharomyces cerevisiae ribosomal protein L30 autoregulates its own expression by binding to a purine-rich internal loop in its pre-mRNA and mRNA. NMR studies of L30 and its RNA complex showed that both the internal loop of the RNA as well as a region of the protein become substantially more ordered upon binding. A crystal structure of a maltose binding protein (MBP)-L30 fusion protein with two copies in the asymmetric unit has been determined. The flexible RNA-binding region in the L30 copies has two distinct conformations, one resembles the RNA bound form solved by NMR and the other is unique. Structure prediction algorithms also had difficulty accurately predicting this region, which is consistent with conformational flexibility seen in the NMR and X-ray crystallography studies. Inherent conformational flexibility may be a hallmark of regions involved in intermolecular interactions.     

Recruiting monomer for dimer formation: resolving the antagonistic mechanisms of novel immune check point inhibitors against Programmed Death Ligand-1 in cancer immunotherapy

The design of small molecule antagonists against Programmed Death Ligand-1 (PD-L1) has been the recent highlight of the immune checkpoint blockade therapy. This interventive approach has been potentiated by the development of BMS compounds; BMS-1001 and BMS-1166, which exert their therapeutic activities by inducing dimerisation of PD-L1; a molecular mechanism that has remained unclear. For the first time, we resolve the dynamical events that underlie the antagonistic mechanisms of BMS-1001 and BMS-1166 when bound to PD-L1 using an all-atom molecular dynamics (MD) simulations approach and free binding energy Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) calculations. Time-scale dynamical findings revealed that upon binding a PD-L1 monomer, the BMS-compounds gradually facilitated the ‘inbound’ motion of another PD-L1 monomer in the same conformational phase space up till dimer formation. Moreover, the non-liganded PD-L1 monomer exhibited the highest structural flexibility and atomistic motions relative to the BMS-liganded monomer as revealed by post-MD trajectory analyses using root mean square deviation (RMSD) and root mean square fluctuations (RMSF) parameters. Trajectory investigations into ligand motions also revealed that the BMS compounds exhibited mechanistic transitions from the monomeric binding site (monomer A) where they were initially bound, to the second monomeric site (monomer B) where they were strongly bound, followed by eventual high-affinity interactions at the tunnel-like binding cleft formed upon the dimerisation of both PD-L1 monomers. These findings present a model that describes the mechanism by which the BMS compounds induce PD-L1 dimerisation and could further enhance the design of highly selective and novel monomeric recruiters of PD-L1 in cancer immunotherapy.     

Insulin assembly damps conformational fluctuations: Raman analysis of amide I linewidths in native states and fibrils

The crystal structure of insulin has been investigated in a variety of dimeric and hexameric assemblies. Interest in dynamics has been stimulated by conformational variability among crystal forms and evidence suggesting that the functional monomer undergoes a conformational change on receptor binding. Here, we employ Raman spectroscopy and Raman microscopy to investigate well-defined oligomeric species: monomeric and dimeric analogs in solution, native T(6) and R(6) hexamers in solution and corresponding polycrystalline samples. Remarkably, linewidths of Raman bands associated with the polypeptide backbone (amide I) exhibit progressive narrowing with successive self-assembly. Whereas dimerization damps fluctuations at an intermolecular beta-sheet, deconvolution of the amide I band indicates that formation of hexamers stabilizes both helical and non-helical elements. Although the structure of a monomer in solution resembles a crystallographic protomer, its encagement in a native assembly damps main-chain fluctuations. Further narrowing of a beta-sheet-specific amide I band is observed on reorganization of insulin in a cross-beta fibril. Enhanced flexibility of the native insulin monomer is in accord with molecular dynamics simulations. Such conformational fluctuations may initiate formation of an amyloidogenic nucleus and enable induced fit on receptor binding.     

Insight into the activity of SARS main protease: Molecular dynamics study of dimeric and monomeric form of enzyme

The phenomenon that SARS coronavirus main protease (SARS M(pro)) dimer is the main functional form has been confirmed by experiment. However, because of the absence of structural information of the monomer, the reasons for this remain unknown. To investigate it, two molecular dynamics (MD) simulations in water for dimer and monomer models have been carried out, using the crystal structure of protomer A of the dimer as the starting structure for the monomer. During the MD simulation of dimer, three interest phenomena of protomer A have been observed: (i) the distance between NE2 of His41 and SG of Cys145 averages 3.72 A, which agrees well with the experimental observations made by X-ray crystallography; (ii) His163 and Glu166 form the "tooth" conformational properties, resulting in the specificity for glutamine at substrate P1 site; and (iii) the substrate-binding pocket formed by loop 140-146 and loop 184-197 is large enough to accommodate the substrate analog. However, during the MD simulation of the monomer complex, the three structural characteristics are all absent, which results directly in the inactivation of the monomer. Throughout the MD simulation of the dimer, the N-terminus of protomer B forms stable hydrogen bonds with Phe140 and Glu166, through which His163, Glu166, and loop 140-146 are kept active form. Furthermore, a water-bridge has been found between the N-terminus of protomer B and Gly170, which stabilizes His172 and avoids it moving toward Tyr161 to disrupt the H-bond between Tyr161 and His163, stabilizing the conformation of His163. The interactions between the N-terminus and another monomer maintain the activity of dimer.     

Brownian and essential dynamics studies of the HIV-1 integrase catalytic domain

The three-dimensional structure of the active site region of the enzyme HIV-1 integrase is not unambiguously known. This region includes a flexible peptide loop that cannot be well resolved in crystallographic determinations. Here we present two different computational approaches with different levels of resolution and on different time-scales to understand this flexibility and to analyze the dynamics of this part of the protein. We have used molecular dynamics simulations with an atomic model to simulate the region in a realistic and reliable way for 1 ns. It is found that parts of the loop wind up after 300 ps to extend an existing helix. This indicates that the helix is longer than in the earlier crystal structures that were used as basis for this study. Very recent crystal data confirms this finding, underlining the predictive value of accurate MD simulations. Essential dynamics analysis of the MD trajectory yields an anharmonic motion of this loop. We have supplemented the MD data with a much lower resolution Brownian dynamics simulation of 600 ns length. It provides ideas about the slow-motion dynamics of the loop. It is found that the loop explores a conformational space much larger than in the MD trajectory, leading to a "gating"-like motion with respect to the active site.     

Dynamic properties of the N-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two 3-dimensional domain-swapped dimers with different quaternary structures. One dimer is characterized by the swapping of the C-terminal region (C-Dimer) and presents a rather loose structure. The other dimer (N-Dimer) exhibits a very compact structure with exchange of the N-terminal helix. Here we report the results of a molecular dynamics/essential dynamics (MD/ED) study carried out on the N-Dimer. This investigation, which represents the first MD/ED analysis on a three-dimensional domain-swapped enzyme, provides information on the dynamic properties of the active site residues as well as on the global motions of the dimer subunits. In particular, the analysis of the flexibility of the active site residues agrees well with recent crystallographic and site-directed mutagenesis studies on monomeric RNase A, thus indicating that domain swapping does not affect the dynamics of the active sites. A slight but significant rearrangement of N-Dimer quaternary structure, favored by the formation of additional hydrogen bonds at subunit interface, has been observed during the MD simulation. The analysis of collective movements reveals that each subunit of the dimer retains the functional breathing motion observed for RNase A. Interestingly, the breathing motion of the two subunits is dynamically coupled, as they open and close in phase. These correlated motions indicate the presence of active site intercommunications in this dimer. On these bases, we propose a speculative mechanism that may explain negative cooperativity in systems preserving structural symmetry during the allosteric transitions.     

Conformational study of a nine residue fragment of the antigenic loop of foot-and-mouth disease virus.

The nine-residue peptide Ac-TASARGDLA-NHMe was selected as model peptide in order to understand the conformational features of the antigenic loop of foot-and-mouth disease virus (FMDV). A throughout exploration of the conformational space has been carried out by means of molecular dynamics (MD) and energy minimization. The calculations have been carried out using the AMBER force field. Solvent effects have been included by an effective dielectric constant of epsilon = 4r. The lowest energy conformation presents a secondary structure constituted by an alpha-helix at the N-terminal end followed by two gamma-turns in the central region. The rest of the accessible minima found present also a high tendency to form gamma-turns. Finally, a 100 ps MD trajectory calculation at 298 K suggest a stability of the secondary structure elements of the lowest energy conformation.     

Structure and atomic fluctuation patterns of potato carboxypeptidase a inhibitor protein

Molecular dynamics (MD) simulation methods were applied to the study of the structural and dynamic fluctuation properties of the potato carboxypeptidase A inhibitor protein (PCI) immersed in a bath of 1259 water molecules. A trajectory of 200 ps was generated at constant temperature and pressure. The crystallographic structure of PCI, as found in its complex with bovine carboxy-peptidase A (CPA), was used to seed the MD simulation. Analyses show that the structure of the PCI core is fairly rigid and stable in itself, and that little deformation is caused by the protein-protein interactions found in the PCI-CPA complex. The N-terminal tail fluctuates to approach the core structure and appears as a relatively disordered region. In contrast, the conformations of the C-terminal tail, which is involved in the inhibitory mechanism, fluctuates in the neighborhood of the X-ray structure in orientations which facilitate CPA binding. Comparison with the structural entries for PCI in water obtained from both 2D-NMR experiments and X-ray data shows that important features of the MD structural results fluctuates between the initial crystal values and those obtained from the NMR solution structure. This fluctuation is not uniform; minor regions move away from the X-ray conformation while they do not approach the NMR conformation. The results reported suggest that the trajectory is long enough to show a behavior that is consistent with the conformational space available to the protein in solution.Abbreviations CPA Carboxypeptidase - DG Distance Geometry - NMR Nuclear Magnetic Resonance - NIS Non Inertial Solvent - MD Molecular Dynamics - PBC Periodic Boundary Conditions - PCI Potato Carboxypeptidase Inhibitor - RMSD Root Mean Square Deviation - a.m.u. Atomic mass units Correspondence to: O. Tapia     

Structure, dynamics, and elasticity of free 16s rRNA helix 44 studied by molecular dynamics simulations

Molecular dynamics (MD) simulations were employed to investigate the structure, dynamics, and local base-pair step deformability of the free 16S ribosomal helix 44 from Thermus thermophilus and of a canonical A-RNA double helix. While helix 44 is bent in the crystal structure of the small ribosomal subunit, the simulated helix 44 is intrinsically straight. It shows, however, substantial instantaneous bends that are isotropic. The spontaneous motions seen in simulations achieve large degrees of bending seen in the X-ray structure and would be entirely sufficient to allow the dynamics of the upper part of helix 44 evidenced by cryo-electron microscopic studies. Analysis of local base-pair step deformability reveals a patch of flexible steps in the upper part of helix 44 and in the area proximal to the bulge bases, suggesting that the upper part of helix 44 has enhanced flexibility. The simulations identify two conformational substates of the second bulge area (bottom part of the helix) with distinct base pairing. In agreement with nuclear magnetic resonance (NMR) and X-ray studies, a flipped out conformational substate of conserved 1492A is seen in the first bulge area. Molecular dynamics (MD) simulations reveal a number of reversible alpha-gamma backbone flips that correspond to transitions between two known A-RNA backbone families. The flipped substates do not cumulate along the trajectory and lead to a modest transient reduction of helical twist with no significant influence on the overall geometry of the duplexes. Despite their considerable flexibility, the simulated structures are very stable with no indication of substantial force field inaccuracies.     

Molecular dynamics of sickle and normal hemoglobins

Molecular dynamics (MD) simulations have been carried out for 62.5 ps on crystal structures of deoxy sickle cell hemoglobin (HbS) and normal deoxy hemoglobin (HbA) using the CHARMM MD algorithm, with a time step of 0.001 ps. In the trajectory analysis of the 12.5–62.5 (50 ps) simulation, oscillations of the radius of gyration and solvent-accessible surface area were calculated. HbS exhibited a general contraction during the simulation, while HbA exhibited a nearly constant size. The average deviations of simulated structures from the starting structures were found to be 1.8 Å for HbA and 2.3 Å for HbS. The average rms amplitudes of atomic motions (atomic flexibility) were about 0.7 Å for HbA and about 1.0 Å for HbS. The amplitudes of backbone motion correlate well with temperature factors derived from x-ray crystallography. A comparison of flexibility between the α- and β-chains in both HbA and HbS indicates that the β-chains generally exhibited greater flexibility than the α-chains, and that the HbS β-chains exhibit greater flexibility in the N-terminal and D- and F-helix regions than do those of HbA. The average amplitude of backbone torsional oscillations was about 9°, a value comparable with that of other simulations, with enhanced torsional oscillation occurring primarily at the ends of helices or in loop regions between helices. Comparison of atomic flexibility and torsional oscillation results suggests that the increased β-chain flexibility results from relatively concerted motions of secondary structure elements. The increased flexibility may play an important role in HbS polymerization. Time course analysis of conformational energy of association, hydrogen bonding and hydrophobic bonding (as calculated from solvent accessibility) shows that all three of these factors contribute to the stability of subunit association for both hemoglobins. © 1993 John Wiley & Sons, Inc.     

Structure and fluctuations of bacteriophage T4 glutaredoxin modelled by molecular dynamics

A 120 ps molecular dynamics (MD) trajectory for bacteriophage T4 glutaredoxin was calculated including non-inertial solvent effects. The potential energy attains an equilibrated regime after the first 20 ps. The r.m.s. difference of all non-hydrogen atoms between X-ray and average MD structures for the regular secondary structure is 0.99A which shows that the MD simulation reproduces the essentials of the structure with high accuracy. Loop displacements are detected, shown by the larger full structure all non-hydrogen atom r.m.s. difference of 1.2A. The fluctuation pattern derived from MD agrees fairly well with that derived from X-ray isotropic temperature factors. The active site is a stable structural region in this MD modellization. Structural changes are put in context with the protein's function.     

Brownian and Essential Dynamics Studies of the HIV-1 Integrase Catalytic Domain

Abstract

The three-dimensional structure of the active site region of the enzyme HIV-1 integrase is not unambiguously known. This region includes a flexible peptide loop that cannot be well resolved in crystallographic determinations. Here we present two different computional approaches with different levels of resolution and on different time-scales to understand this flexibility and to analyze the dynamics of this part of the protein. We have used molecular dynamics simulations with an atomic model to simulate the region in a realistic and reliable way for 1 ns. It is found that parts of the loop wind up after 300 ps to extend an existing helix. This indicates that the helix is longer than in the earlier crystal structures that were used as basis for this study. Very recent crystal data confirms this finding, underlining the predictive value of accurate MD simulations. Essential dynamics analysis of the MD trajectory yields an anharmonic motion of this loop. We have supplemented the MD data with a much lower resolution Brownian dynamics simulation of 600 ns length. It provides ideas about the slow-motion dynamics of the loop. It is found that the loop explores a conformational space much larger than in the MD trajectory, leading to a “gating”-like motion with respect to the active site.     

Importance of folded monomer and extended antiparallel dimer structures as enkephalin active conformation. Molecular dynamics simulations of [Met5]enkephalin in water

Simulations of the molecular dynamics of the [Met5]enkephalin monomer and dimer structures in water have been carried out. The dynamic trajectories have been analyzed in terms of the distances between intra- or intermolecular polar atoms. The time-correlated conformational transitions of an extended monomer structure have been converged into a stationary state among the beta-bend folded forms. However, the dynamics simulation of an extended antiparallel dimer structure has shown no noticeable conformation change. These results imply that both the beta-bend monomer and the extended dimer structures exist together as the fundamental conformation of enkephalins.     

Molecular-dynamics simulations of [Met5]- and [D-Ala2,Met5]-enkephalins. Biological implication of monomeric folded and dimeric unfolded conformations.

To investigate the biologically active conformation of enkephalin, molecular-dynamics simulations were applied to [Met5]- and [D-Ala2,Met5]-enkephalins. The dynamic trajectory of monomeric extended [Met5]-enkephalin was analysed in terms of relative mobility between respective torsions of backbone chain. After 10 ps of the dynamics simulation, the conformational transition was converged into a stationary state among the beta-bend folded forms, where they are stabilized by several intramolecular hydrogen-bond formations. Similar conformational transition was also observed in the dynamics simulation of [D-Ala2,Met5]enkephalin, which is a more mu-receptor-specific peptide than [Met5]enkephalin. The geometrical correspondence between the monomeric enkephalin conformation in the stationary state and morphine molecule (a mu-specific rigid opiate) was surveyed by virtue of the triangular substructures generated by choosing three functional atoms in each molecule, and good resemblances were observed. On the other hand, the dynamics simulation of the antiparallel extended [Met5]enkephalin dimer showed a trajectory different from that of the monomeric one. Two intermolecular hydrogen bonds at Tyr1 (NH3+)...Met5(CO2-) end residues were held throughout the 100 ps simulation, the dimeric structure being consequently kept. The conformational transition of the backbone chains from the antiparallel extended form to the twisted one took place via an intermediate state. Many conformations revealed during the dynamics simulation showed that the relative orientations of each two Tyr1, Gly3, Phe4 and Met5 residues in the dimer are nearly related by a pseudo-C2-symmetry respectively, and both halves of the dimer structure could be further fitted to the monomeric folded enkephalin conformation. The monomeric and dimeric conformations of enkephalin at their stationary states are discussed in relation to the substrate-specificity for mu- and delta-opioid receptors.     

Conformational ensemble of an intrinsically flexible loop in mitochondrial import protein Tim21 studied by modeling and molecular dynamics simulations

BackgroundTim21, a subunit of a highly dynamic translocase of the inner mitochondrial membrane (TIM23) complex, translocates proteins by interacting with subunits in the translocase of the outer membrane (TOM) complex and Tim23 channel in the TIM23 complex. A loop segment in Tim21, which is in close proximity of the binding site of Tim23, has different conformations in X-ray, NMR and new crystal contact-free space (CCFS) structures. MD simulations can provide information on the structure and dynamics of the loop in solution.MethodsThe conformational ensemble of the loop was characterized using loop modeling and molecular dynamics (MD) simulations.ResultsMD simulations confirmed mobility of the loop. Multidimensional scaling and clustering were used to characterize the dynamic conformational ensemble of the loop. Free energy landscape showed that the CCFS crystal structure occupied a low energy region as compared to the conventional X-ray crystal structure. Analysis of crystal packing indicates that the CCFS provides larger conformational space for the motions of the loop.ConclusionsOur work reported the conformational ensemble of the loop in solution, which is in agreement with the structure obtained from CCFS approach. The combination of the experimental techniques and computational methods is beneficial for studying highly flexible regions of proteins.General significanceComputational methods, such as loop modeling and MD simulations, have proved to be useful for studying conformational flexibility of proteins. These methods in integration with experimental techniques such as CCFS has the potential to transform the studies on flexible regions of proteins.     

Domain motions of the Mip protein from Legionella pneumophila

The homodimeric 45.6 kDa (total mass) Mip protein, a virulence factor from Legionella pneumophila, was investigated with solution NMR spectroscopy and molecular dynamics (MD) simulations. Two Mip monomers are dimerized via an N-terminal helix bundle that is connected via a long alpha-helix to a C-terminal FKBP domain in each subunit. More than 85% of the amino acids were identified in triple-resonance NMR spectra. (15)N relaxation analysis showed a bimodal distribution of R(1)/R(2) values, with the lower ratio in the N-terminal domain. Relaxation dispersion measurements confirmed that these reduced ratios did not originate from conformational exchange. Thus, two different correlation times (tau(c)) can be deduced, reflecting partly uncoupled motions of both domains. Relaxation data of a Mip(77)(-)(213) monomer mutant were similar to those observed in the dimer, corroborating that the FKBP domain, including part of the connecting helix, behaves as one dynamic entity. MD simulations (18 ns) of the Mip dimer also yielded two different correlation times for the two domains and thus confirm the independence of the domain motions. Principal component analysis of the dihedral space covariance matrix calculated from the MD trajectory suggests a flexible region in the long connecting helix that acts as a hinge between the two domains. Such motion provides a possible explanation of how Mip can bind to complex molecular components of the extracellular matrix and mediate alveolar damage and bacterial spread in the lung.