Mass-weighted molecular dynamics simulation of cyclic polypeptides.

A modified molecular dynamics (MD) method in which atomic masses are weighted was developed previously for studying the conformational flexibility of neuroregulating tetrapeptide Phe-Met-Arg-Phe-amide (FMRF-amide). The method has now been applied to longer and constrained molecules, namely a disulfide-linked cyclic hexapeptide, c[CYFQNC], and its linear and "pseudo-cyclic" analogues. The sampling of dehedral conformational space of teh linear hexapeptide in mass-weighted MD simulations was found to be improved significantly over conventional MD simulations, as in the case of the shorter FMRF-amide molecule studied previously. In the cyclic hexapeptide, the internal constraint of the molecule due to the intramolecular disulfide bond (hence the absence of free terminals in the molecule) does not adversely affect the significant improvement of conformational sampling in mass-weighted MD simulations over normal MD simulations. The pseudo-cyclic polypeptide is identical to the linear CYFQNC molecule in amino acid sequence (i.e., side chains of the cysteine residues are reduced), but the positions of its two terminal heavy atoms were held fixed in space such that the molecule has a nearly cyclic conformation. For this molecule, the mass-weighted MD simulation generated a wide range of polypeptide backbone conformations covering the internal dihedral degrees of freedom; moreover, the physical space of the pseudo-cyclic structure was also sampled in a complete revolution of the entire molecular fragment about the two fixed termini during the simulation. These characteristics suggest that mass-weighted MD can also be an extremely useful method for conformational analyses of constrained molecules and, in particular, for modeling loops on protein surfaces.     

Mass-weighted molecular dynamics simulation of the protein-ligand complex of rhizopuspepsin and inhibitor.

The mass-weighted molecular dynamics simulation method was developed previously for sampling the multidimensional conformational space of linear and cyclic polypeptides and studying their conformational flexibility. Herein results from molecular dynamics simulations of the protein-ligand complex of the aspartyl protease rhizopuspepsin and a polypeptide inhibitor are reported. The dihedral conformational space sampling for the linear peptide inhibitor in situ was found to be increased in the mass-weighted simulation as in other molecular systems previously studied. More significantly, the physical space of the enzyme binding pocket was also sampled efficiently in the simulations and multiple binding sites were identified for the inhibitor. These results suggest that it may be possible now to study, by computer simulations, the putative initial enzyme-inhibitor complex suggested experimentally from the time-dependent kinetics of enzyme inhibition by slow-binding inhibitors (Morrison, J. F., and C. T. Walsh. 1988. Adv. Enzymol. 61:201), and/or conformational substates in protein-ligand complexes suggested in the study of reassociation dynamics of myoglobin and carbon monoxide following photolysis (Austin, R. H., K. W. Beeson, L. Eisenstein, H. Frauenfelder, and I. C. Gunsalus. 1975. Biochemistry. 14:5355). Moreover, the intermediate binding steps and the molecular flexibility of the inhibitor shown in the MWMD simulation may have crucial roles in the ligand binding process.     

Lipid bilayer polypeptide interactions studied by molecular dynamics simulation

A model membrane with a polypeptide alpha-helix inserted has been simulated by molecular dynamics at a temperature well above the gel/liquid crystalline phase transition temperature. Order parameters of the lipids and other equilibrium and dynamic quantities have been calculated. Three systems, polyglycine constrained into an alphahelical configuration, glycophorin with similarly conformationally constrained backbone and finally glycophorin free to change its backbone conformation, have been studied. In all cases there was an ordering of the chains close to the helix. This effect was, however, much smaller for glycophorin with its rather bulky side chains than for polyglycine. The dynamics of the lipids were affected by the neighbouring helix, not drastically however. Lateral diffusion and reorientational time correlations of lipids close to the helix were slower than for the bulk ones, but not more than two or three times. Thus, we did not find any evidence of bound or frozen boundary lipids.     

Modeling ganglioside headgroups by conformational analysis and molecular dynamics

The conformations and dynamics of gangliosides GM1, GM2, 6-GM2 and GM4 have been studied by computational means, and the results compared to NMR data. Unconstrained conformational searches were run using the AMBER* force field augmented by MNDO derived parameters for the Neu5Ac anomeric torsion, the GB/SA water solvation model, and the MC/EM alogorithm; extended (10–12[emsp4 ]ns) dynamic simulations in GB/SA water were performed with the MC/SD protocol, and the stored structures were minimized. The overall mobility of the Neu5Ac2,3Gal linkage and the position of its minimum energy conformation have been shown to depend mainly on the presence or the absence of a GalNAc residue at the adjacent position. The best quantitative agreement with the available NOE data was achieved after minimization of the structures stored during the MC/SD dynamic runs. The latter protocol appears to reproduce satisfactorily the available experimental data, and can be used with confidence to build three-dimensional models of ganglioside headgroups.     

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.     

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.     

Exploring the conformational dynamics and membrane interactions of PorB from C. glutamicum: a multi-scale molecular dynamics simulation study

Members of the gram-positive mycolata bacteria have unusual cell envelopes which help them to avoid the immune system and the effects of most antibiotics, whilst rendering them permeable to solutes of importance in industrial bioconversion. It is therefore of interest to understand the molecular mechanisms for this selective permeability. PorB is an unusual porin from the outer membrane (OM) of Corynebacterium glutamicum. It has been proposed as an atypical α-helical, symmetrical homo-pentameric architecture, with an unusual distribution of polar amino acids on its surface. The proposed structure is too short to traverse a typical phospholipid bilayer, in contrast with the β-barrel porins of Gram-negative bacteria. Nevertheless, it has been shown to form small anion-selective channels in membranes typical of Escherichia coli. To further understand its function, we have performed ~400 ns of all-atom and ~270 μs of coarse-grained simulations of PorB in a range of membrane mimetic and phospholipid milieus. Our results suggest that PorB can undergo spontaneous conformational rearrangements that allow it to adapt to its local lipid environment. We speculate that the increased flexibility of this α-helical porin in comparison with rigid β-barrels may be an adaptation for the heterogeneous mycolic OM, and explains its demonstrated ability to form measurable pores with phospholipid membranes.     

LeuT conformational sampling utilizing accelerated molecular dynamics and principal component analysis

Monoamine transporters (MATs) function by coupling ion gradients to the transport of dopamine, norepinephrine, or serotonin. Despite their importance in regulating neurotransmission, the exact conformational mechanism by which MATs function remains elusive. To this end, we have performed seven 250 ns accelerated molecular dynamics simulations of the leucine transporter, a model for neurotransmitter MATs. By varying the presence of binding-pocket leucine substrate and sodium ions, we have sampled plausible conformational states representative of the substrate transport cycle. The resulting trajectories were analyzed using principal component analysis of transmembrane helices 1b and 6a. This analysis revealed seven unique structures: two of the obtained conformations are similar to the currently published crystallographic structures, one conformation is similar to a proposed open inward structure, and four conformations represent novel structures of potential importance to the transport cycle. Further analysis reveals that the presence of binding-pocket sodium ions is necessary to stabilize the locked-occluded and open-inward conformations.     

Docking and molecular dynamics simulation studies on glycation-induced conformational changes of human paraoxonase 1

Human paraoxonase 1 (huPON1) is a calcium-dependent esterase responsible for hydrolysis of a wide variety of substrates including organophosphates, esters, lactones, and paraoxon. Although its natural substrate is unknown, the action of PON as an antioxidant is well documented. Because recent reports have suggested glycation may induce reduced PON activity in diabetes, we investigated the structural features of huPON1 and its glycated mutant by template-based modeling, docking, and molecular dynamics (MD) simulations. Our results corroborated the importance of the His115–His134 dyad in both the lactonase and paraoxonase activity of huPON1. Structural alterations in the glycated model reflected weak interactions between the docked substrate and the active site cleft. We also used MD simulation to gain insight into glycation-induced conformational changes of huPON1 and the implication of this on depleted enzymatic activity. The catalytic calcium found on the surface interacts with the side chain oxygen of residues, including Asn224, Asn270, Asn168, Asp269, and Glu53, and this interaction with the respective residues undergoes minor displacement on glycation. The root-mean-square fluctuation had high motional flexibility in the non-glycated model whereas the conformation of the glycated structure was comparatively stable. Our findings emphasize the consequence of glycation-induced alterations and their effect on overall enzymatic 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.     

Path of nascent polypeptide in exit tunnel revealed by molecular dynamics simulation of ribosome

Molecular dynamics simulations were carried out on Thermus thermophilus 70S ribosome with and without a nascent polypeptide inside the exit tunnel. Modeling of the polypeptide in the tunnel revealed two possible paths: one over Arg⁹² of L22 and one under (from the viewpoint of 50S on top of 30S). A strong interaction between L4 and Arg⁹² was observed without the polypeptide and when it passed over Arg⁹². However, when the polypeptide passed under, Arg⁹² repositioned to interact with Ade²⁰⁵⁹ of 23S rRNA. Using steered molecular dynamics the polypeptide could be pulled through the L4-L22 constriction when situated under Arg⁹², but did not move when over. These results suggest that the tunnel is closed by the Arg⁹²-L4 interaction before elongation of the polypeptide and the tunnel leads the entering polypeptide from the peptidyl transferase center to the passage under Arg⁹², causing Arg⁹² to switch to an open position. It is possible, therefore, that Arg⁹² plays the role of a gate, opening and closing the tunnel at L4-L22. There is some disagreement over whether the tunnel is dynamic or rigid. At least within the timescale of our simulations conformational analysis showed that global motions mainly involve relative movement of the 50S and 30S subunits and seem not to affect the conformation of the tunnel.     

Large-scale conformational changes of Trypanosoma cruzi proline racemase predicted by accelerated molecular dynamics simulation

Chagas' disease, caused by the protozoan parasite Trypanosoma cruzi (T. cruzi), is a life-threatening illness affecting 11-18 million people. Currently available treatments are limited, with unacceptable efficacy and safety profiles. Recent studies have revealed an essential T. cruzi proline racemase enzyme (TcPR) as an attractive candidate for improved chemotherapeutic intervention. Conformational changes associated with substrate binding to TcPR are believed to expose critical residues that elicit a host mitogenic B-cell response, a process contributing to parasite persistence and immune system evasion. Characterization of the conformational states of TcPR requires access to long-time-scale motions that are currently inaccessible by standard molecular dynamics simulations. Here we describe advanced accelerated molecular dynamics that extend the effective simulation time and capture large-scale motions of functional relevance. Conservation and fragment mapping analyses identified potential conformational epitopes located in the vicinity of newly identified transient binding pockets. The newly identified open TcPR conformations revealed by this study along with knowledge of the closed to open interconversion mechanism advances our understanding of TcPR function. The results and the strategy adopted in this work constitute an important step toward the rationalization of the molecular basis behind the mitogenic B-cell response of TcPR and provide new insights for future structure-based drug discovery.     

Probing the ATP-induced conformational flexibility of the PcrA helicase protein using molecular dynamics simulation

Helicases are enzymes that unwind double-stranded DNA (dsDNA) into its single-stranded components. It is important to understand the binding and unbinding of ATP from the active sites of helicases, as this knowledge can be used to elucidate the functionality of helicases during the unwinding of dsDNA. In this work, we investigated the unbinding of ATP and its effect on the active-site residues of the helicase PcrA using molecular dynamic simulations. To mimic the unbinding process of ATP from the active site of the helicase, we simulated the application of an external force that pulls ATP from the active site and computed the free-energy change during this process. We estimated an energy cost of ~85 kJ/mol for the transformation of the helicase from the ATP-bound state (1QHH) to the ATP-free state (1PJR). Unbinding led to conformational changes in the residues of the protein at the active site. Some of the residues at the ATP-binding site were significantly reoriented when the ATP was pulled. We observed a clear competition between reorientation of the residues and energy stabilization by hydrogen bonds between the ATP and active-site residues. We also checked the flexibility of the PcrA protein using a principal component analysis of domain motion. We found that the ATP-free state of the helicase is more flexible than the ATP-bound state.     

Crystal structure,conformational analysis,and molecular dynamics of tetra-O-methyl-(+) -catechin

The structure of tetra-O-methyl- (+) -catechin has been determined in the crystalline state. Two independent molecules, denoted structure A and structure B, exist in the unit cell. Crystals are triclinic, space group P1, a = 4.8125(2) Å, b = 12.9148(8) Å, c = 13.8862(11) Å, α = 86.962(6) °, β = 89.120(5)°, γ = 88.044(5)°, Z = 2, D_c = 1.336 g cm^?3, R = 0.033 for 6830 observations. The heterocyclic rings of the crystal structures are compared to previous results for 8-bromotetra-O-methyl-(+)-catechin, penta-O-acetyl-(+)-catechin, and (?) -epicatechin. One of the two molecules has a heterocyclic ring conformation similar to that observed previously for (?)-epicatechin, and the other has a heterocyclic ring conformation similar to one predicted earlier in a theoretical analysis of dimers of (+)-catechin and (?) -epicatechin. Both structure A and structure B in the crystal have heterocyclic ring conformations that place the dimethoxyphenyl substituent at C(2) in the equatorial position. However, this heterocyclic ring conformation does not explain the proton nmr coupling constant measured in solution. Molecular dynamics simulations show an equatorial ? axial interconversion of the heterocyclic ring, which can explain the nmr results. © 1993 John Wiley & Sons, Inc.     

Diffusional behavior and guest conformational analysis of hexadecane-1,16-diol and hexadecane in urea crystal model via molecular dynamics simulation approach

Diffusion at the atomic or molecular level is a source of many physical, chemical, and biological processes taking place in plentiful materials. This work is an endeavor toward investigating the diffusional behavior of two different type of guests, hexadecane-1,16-diol and hexadecane enclathration in urea tunnel architecture, whereby the correlation of the diffusion mechanism with the guest’s structural and conformational properties is explored. To carry out this study, molecular dynamics simulation approach is adopted. It is found that hexadecane-1,16-diol exhibit slower diffusion with an average diffusion coefficient value \( \sim 1.83\times {10}^{-10} \), where hexadecane diffuse more rapidly with an average diffusion coefficient value \( \sim 2.58\times {10}^{-9} \). It is also observed that the structural properties influence the guest’s travel distance and torsion angle distribution of the trans and gauche conformational proportion. Furthermore, the observed high energy barrier accounted for hexadecane-1,16-diol and low energy barrier for hexadecane along urea tunnel systems was analyzed. The comparison of our obtained results are in close agreement with the available experimental measurements, i.e., gauche proportion properties between two different guest molecules correlate well with Raman spectroscopy investigation on α,ω-dihalogenoalkane/urea inclusion compounds. Our calculations also successfully endorse the structure-property relation between the two systems.     

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.     

Electrostatic interaction-mediated conformational changes of adipocyte fatty acid binding protein probed by molecular dynamics simulation

Abstract

Adipocyte fatty acid binding protein (A-FABP) is a potential drug target for treatment of diabetes, obesity and atherosclerosis. Molecular dynamics (MD) simulations, principal component (PC) analysis and binding free energy calculations were combined to probe effect of electrostatic interactions of residues R78, R106 and R126 with inhibitors ZGB, ZGC and IBP on structural stability of three inhibitor/A-FABP complexes. The results indicate that mutation R126A produces significant influence on polar interactions of three inhibitors with A-FABP and these interactions are main force for driving the conformational change of A-FABP. Analyses on hydrogen bond interactions show that the decrease in hydrogen bonding interactions of residues R126 and Y128 with three inhibitors and the increase in that of K58 with inhibitors ZGC and IBP in the R126A mutated systems mostly regulate the conformational changes of A-FABP. This work shows that R126A can generate a significant perturbation on structural stability of A-FABP, which implies that R126 is of significance in inhibitor bindings. We expect that this study can provide a theoretical guidance for design of potent inhibitors targeting A-FABP.

Communicated by Ramaswamy H. Sarma     

Molecular dynamics simulation of transmembrane polypeptide orientational fluctuations

The orientation and motion of a model lysine-terminated transmembrane polypeptide were investigated by molecular dynamics simulation. Recent 2H NMR studies of synthetic polypeptides with deuterated alanine side chains suggest that such transmembrane polypeptides undergo fast, axially symmetric reorientation about the bilayer normal but have a preferred average azimuthal orientation about the helix axis. In this work, interactions that might contribute to this behavior were investigated in a simulated system consisting of 64 molecules of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and one alpha-helical polypeptide with the sequence acetyl-KK-(LA)11-KK-amide. In one simulation, initiated with the peptide oriented along the bilayer normal, the system was allowed to evolve for 8.5 ns at 1 atm of pressure and a temperature of 55 degrees C. A second simulation was initiated with the peptide orientation chosen to match a set of experimentally observed alanine methyl deuteron quadrupole splittings and allowed to proceed for 10 ns. Simulated alanine methyl group orientations were found to be inequivalent, a result that is consistent with 2H NMR observations of specifically labeled polypeptides in POPC bilayers. Helix tilt varied substantially over the durations of both simulations. In the first simulation, the peptide tended toward an orientation about the helix axis similar to that suggested by experiment. In the second simulation, orientation about the helix axis tended to return to this value after an excursion. These results provide some insight into how interactions at the bilayer surface can constrain reorientation about the helix axis while accommodating large changes in helix tilt.