Vibrational analysis of a solvated green fluorescent protein chromophore

Resonance Raman (RR) spectra of green fluorescent protein (GFP) model chromophores in solution have been simulated with the CASSCF/MM methodology. Although several reports on vibrational analysis of GFP model chromophores have been recently published, the RR spectra were simulated for the first time in explicit solution with the inclusion of the counterion, as these effects are crucial for unambiguously reproducing the vibrational band assignment in the anionic form of the GFP chromophore. This strategy allows for a one-to-one correspondence of the calculated vibrational modes to the observed RR bands, concerning both the location and intensity pattern. In addition, these simulations were complemented with total energy distribution calculations to aid in the unambiguous assignment of the measured spectra. The current study helps to clarify some of the previous RR bands assignments as well as producing some new assignment for the anionic form of GFP chromophore. The explicit solvent simulations and PCM-based calculations are compared to the measured spectra, and these results demonstrate that explicit solvent simulations provide better agreement with experiment, both in terms of vibrational frequencies and intensity distribution. Figure a Correlation of explicit hydration calculations (CASSCF/6-31G*/MM) for the HBI model chromophore and experimental RR data [21]; slope = 0.982, intercept = 27.210 and regression coefficient = 0.997. b Correlation of implicit PCM calculations (CASSCF/6-31G*) for the HBI model chromophore and experimental RR data [21], slope = 1.017, intercept = −48.838 and regression coefficient = 0.984     

Vibrational circular dichroism of tetraphenylporphyrin in peptide complexes? A computational study

The Raman and absorption spectra of tetraphenylporphyrin (TPP) were calculated and compared to experiment. The computation was based on the harmonic molecular force field and electric tensors obtained ab initio at the BPW91/6-31G* level. Good agreement was found between experimental and calculated frequencies and intensities. In order to estimate whether induced optical activity in chiral complexes interferes with the signal of peptide vibrations, the vibrational circular dichroism (VCD) spectra of TPP were simulated. The magnetic field perturbation theory (MFP) and the gauge-invariant atomic orbitals (GIAO) were used for the simulation. Such spectra were compared to theoretical VCD intensities of a model tripeptide as well to experimental spectra of a complex of the peptide and tetrakis(p-sulfonatophenyl)porphyrin (TSPP). No significant contribution to VCD signal from the TPP residue was found in experimental spectra. Thus, possible peptide conformational changes occurring during the complexation can be monitored directly in the amide I frequency region.     

Weight average approaches for predicting dynamical properties of biomolecules

Recent advances in atomistic molecular dynamics (MD) simulations of biomolecules allow us to explore their conformational spaces widely, observing large-scale conformational fluctuations or transitions between distinct structures. To reproduce or refine experimental data using MD simulations, structure ensembles, which are characterized by multiple structures and their statistical weights on the rugged free-energy landscapes, are often used. Here, we summarize weight average approaches for various experimental measurements. Weight average approaches are now applied to hybrid quantum mechanics/molecular mechanics MD simulations to predict fast vibrational motions in a protein with a high accuracy for better understanding of molecular functions from atomic structures.     

Computer simulations of the flexibility of a series of synthetic cyclic peptide analogues

The flexibility of a series of cyclic peptides derived from the epitope of a snake toxin is investigated using computer simulation techniques. Molecular dynamics (MD) simulations and vibrational analyses are performed on chemically constrained peptides modeled on the parent protein loop. In the 50 ps MD simulations, large variations in the atomic fluctuations are seen between the peptides, and can be related to the nature of the chemical constraints present in the molecules. Normal mode analyses are performed on energy-minimized configurations derived from the dynamics trajectories. The atomic fluctuations calculated from the normal modes are about 30% of those of the molecular dynamics for the more flexible peptides and 70% for the more constrained molecules. The calculated differences in flexibility between the molecules are much less significant in the harmonic approximation. © 1993 John Wiley & Sons, Inc.     

Developing polarized protein-specific charges for protein dynamics: MD free energy calculation of pKa shifts for Asp26/Asp20 in thioredoxin

Ab initio quantum mechanical calculation of protein in solution is carried out to generate polarized protein-specific charge(s) (PPC) for molecular dynamics (MD) stimulation of protein. The quantum calculation of protein is made possible by developing a fragment-based quantum chemistry approach in combination with the implicit continuum solvent model. The computed electron density of protein is utilized to derive PPCs that represent the polarized electrostatic state of protein near the native structure. These PPCs are atom-centered like those in the standard force fields and are thus computationally attractive for molecular dynamics simulation of protein. Extensive MD simulations have been carried out to investigate the effect of electronic polarization on the structure and dynamics of thioredoxin. Our study shows that the dynamics of thioredoxin is stabilized by electronic polarization through explicit comparison between MD results using PPC and AMBER charges. In particular, MD free-energy calculation using PPCs accurately reproduced the experimental value of pK_a shift for ionizable residue Asp²⁶ buried inside thioredoxin, whereas previous calculations using standard force fields overestimated pK_a shift by twice as much. Accurate prediction of pK_a shifts by rigorous MD free energy simulation for ionizable residues buried inside protein has been a significant challenge in computational biology for decades. This study presented strong evidence that electronic polarization of protein plays an important role in protein dynamics.     

Incoherent neutron scattering of copper azurin: a comparison with molecular dynamics simulation results

The low-frequency dynamics of copper azurin has been studied at different temperatures for a dry and deuterium hydrated sample by incoherent neutron scattering and the experimental results have been compared with molecular dynamics (MD) simulations carried out in the same temperature range. Experimental Debye-Waller factors are consistent with a dynamical transition at approximately 200 K which appears partially suppressed in the dry sample. Inelastic and quasielastic scattering indicate that hydration water modulates both vibrational and diffusive motions. The low-temperature experimental dynamical structure factor of the hydrated protein shows an excess of inelastic scattering peaking at about 3 meV and whose position is slightly shifted downwards in the dry sample. Such an excess is reminiscent of the “boson peak” observed in glass-like materials. This vibrational peak is quite well reproduced by MD simulations, although at a lower energy. The experimental quasielastic scattering of the two samples at 300 K shows a two-step relaxation behaviour with similar characteristic times, while the corresponding intensities differ only by a scale factor. Also, MD simulations confirm the two-step diffusive trend, but the slow process seems to be characterized by a decay faster than the experimental one. Comparison with incoherent neutron scattering studies carried out on proteins having different structure indicates that globular proteins display common elastic, quasielastic and inelastic features, with an almost similar hydration dependence, irrespective of their secondary and tertiary structure. Received: 12 October 1998 / Revised version: 19 February 1999 / Accepted: 1 March 1999     

Interactions of the Anticancer Drug Tamoxifen with Lipid Membranes

Interactions of the hydrophobic anticancer drug tamoxifen (TAM) with lipid model membranes were studied using calcein-encapsulated vesicle leakage, attenuated total reflection Fourier transform infrared (FTIR) spectroscopy, small-angle neutron scattering (SANS), atomic force microscopy (AFM) based force spectroscopy, and all-atom molecular dynamics (MD) simulations. The addition of TAM enhances membrane permeability, inducing calcein to translocate from the interior to the exterior of lipid vesicles. A large decrease in the FTIR absorption band’s magnitude was observed in the hydrocarbon chain region, suggesting suppressed bond vibrational dynamics. Bilayer thickening was determined from SANS data. Force spectroscopy measurements indicate that the lipid bilayer area compressibility modulus K_A is increased by a large amount after the incorporation of TAM. MD simulations show that TAM decreases the lipid area and increases chain order parameters. Moreover, orientational and positional analyses show that TAM exhibits a highly dynamic conformation within the lipid bilayer. Our detailed experimental and computational studies of TAM interacting with model lipid membranes shed new light on membrane modulation by TAM.     

Femtosecond formation dynamics of primary photoproducts of visual pigment rhodopsin

The coherent 11-cis-retinal photoisomerization dynamics in bovine rhodopsin was studied by femtosecond time-resolved laser absorption spectroscopy at 30-fs resolution. Femtosecond pulses of 500, 535, and 560 nm wavelength were used for rhodopsin excitation to produce different initial Franck-Condon states and relevant distinct values of the vibrational energy of the molecule in its electron excited state. Time evolution of the photoinduced rhodopsin absorption spectra was monitored after femtosecond excitation in the spectral range of 400–720 nm. Oscillations of the time-resolved absorption signals of rhodopsin photoproducts represented by photorhodopsin₅₇₀ with vibrationally-excited all-trans-retinal and rhodopsin₄₉₈ in its initial state with vibrationally-excited 11-cis-retinal were studied. These oscillations reflect the dynamics of coherent vibrational wave-packets in the ground state of photoproducts. Fourier analysis of these oscillatory components has revealed frequencies, amplitudes, and initial phases of different vibrational modes, along which the motion of wave-packets of both photoproducts occurs. The main vibrational modes established are 62, 160 cm⁻¹ and 44, 142 cm⁻¹ for photorhodopsin₅₇₀ and for rhodopsin₄₉₈, respectively. These vibrational modes are directly involved in the coherent reaction under the study, and their amplitudes in the power spectrum obtained through the Fourier transform of the kinetic curves depend on the excitation wavelength of rhodopsin.     

Vibrational circular dichroism and IR spectral analysis as a test of theoretical conformational modeling for a cyclic hexapeptide

A model cyclohexapeptide, cyclo-(Phe-(D)Pro-Gly-Arg-Gly-Asp) was synthesized and its IR and VCD spectra were used as a test of density functional theory (DFT) level predictions of spectral intensities for a peptide with a nonrepeating but partially constricted conformation. Peptide structure and flexibility was estimated by molecular dynamics (MD) simulations and the spectra were simulated using full quantum mechanical (QM) approaches for the complete peptide and for simplified models with truncated side chains. After simulated annealing, the backbone conformation of the ring structure is relatively stable, consisting of a normal beta-turn and a tight loop (no H-bond) which does not vary over short trajectories. Only in quite long MD runs at high temperatures do other conformations appear. MD simulations were carried out for the cyclic peptide in water and in TFE, which match experimental solvents, as well as with and without protonation of the Asp carboxyl group. DFT spectral simulations were made using the annealed structure and were extended to include basis set variation, to determine an optimal computational approach, and solvent simulation with a polarized continuum model (PCM). Stepwise full DFT simulation of spectra was done for various sequences with the same backbone geometry but based on (1) solely Gly residues, (2) Ala substitution except Gly and Pro, and (3) complete sequences with side chains. Additionally, a selection of structures was used to compute IR and VCD spectra with the optimal method to determine structural variation effects. The side chains, especially the Asp-COOH and Arg-NH(2) transitions, had an impact on the computed amide frequencies, IR intensities and VCD pattern. Since experimentally these groups would have little chirality, due to conformational variation, they do not impact the observed VCD spectra. Correcting for frequency shifts, the Ala model for the cyclopeptide gives the clearest representation of the amide VCD. The experimental sign pattern for the amide I' band in D(2)O and also the sharper, more intense amide I VCD band in TFE was seen to some degree in one conformer with Type II' turns, but the data favor a mix of structures.     

166 Millisecond-long molecular dynamics simulations of proteins on a special-purpose machine

Molecular dynamics (MD) simulation has long been recognized as a potentially powerful tool for understanding the structural, dynamic, and functional characteristics of proteins at an atomic level of detail. Many biologically important phenomena, however, occur over timescales that have previously fallen far outside the reach of MD technology. We have constructed a specialized, massively parallel machine, called Anton, that is capable of performing all-atom simulations of proteins in an explicitly represented solvent environment at a speed roughly two orders of magnitude beyond that of the previous state of the art. Using novel algorithms developed within our lab, the machine has now simulated the behavior of a number of proteins for periods as long as two milliseconds – approximately 200 times the length of the longest such simulation previously published. Such simulations have allowed us to observe and analyze key characteristics of the dynamics of proteins (including central elements of the process of protein folding) that were previously inaccessible to both computational and experimental study.     

Thermal stability of single‐domain antibodies estimated by molecular dynamics simulations

Single‐domain antibodies (sdAbs) function like regular antibodies, however, consist of only one domain. Because of their low molecular weight, sdAbs have advantages with respect to production and delivery to their targets and for applications such as antibody drugs and biosensors. Thus, sdAbs with high thermal stability are required. In this work, we chose seven sdAbs, which have a wide range of melting temperature (T_m) values and known structures. We applied molecular dynamics (MD) simulations to estimate their relative stability and compared them with the experimental data. High‐temperature MD simulations at 400 K and 500 K were executed with simulations at 300 K as a control. The fraction of native atomic contacts, Q, measured for the 400 K simulations showed a fairly good correlation with the T_m values. Interestingly, when the residues were classified by their hydrophobicity and size, the Q values of hydrophilic residues exhibited an even better correlation, suggesting that stabilization is correlated with favorable interactions of hydrophilic residues. Measuring the Q value on a per‐residue level enabled us to identify residues that contribute significantly to the instability and thus demonstrating how our analysis can be used in a mutant case study.     

Low-frequency vibrational modes and infrared absorbance of red, blue and green opsin

Vibrational excitations of low-frequency collective modes are essential for functionally important conformational transitions in proteins. We carried out an analysis of the low-frequency modes in the G protein coupled receptors (GPCR) family of cone opsins based on both normal-mode analysis and molecular dynamics (MD) simulations. Power spectra obtained by MD can be compared directly with normal modes. In agreement with existing experimental evidence related to transmembrane proteins, cone opsins have functionally important transitions that correspond to approximately 950 modes and are found below 80 cm⁻¹. This is in contrast to bacteriorhodopsin and rhodopsin, where the important low-frequency transition modes are below 50 cm⁻¹. We find that the density of states (DOS) profile of blue opsin in a solvent (e.g. water) has increased populations in the very lowest frequency modes (<15 cm⁻¹); this is indicative of the increased thermostability of blue opsin. From our work we found that, although light absorption behaves differently in blue, green and red opsins, their low-frequency vibrational motions are similar. The similarities and differences in the domain motions of blue, red and green opsins are discussed for several representative modes. In addition, the influence of the presence of a solvent is reported and compared with vacuum spectra. We thus demonstrate that terahertz spectroscopy of low-frequency modes might be relevant for identifying those vibrational degrees of freedom that correlate to known conformational changes in opsins. An erratum to this article can be found at     

CDSpecTech: A single software suite for multiple chiroptical spectroscopic analyses

The program CDSpecTech was developed to facilitate the analysis of chiroptical spectra, which include the following: vibrational circular dichroism (VCD) and corresponding vibrational absorption (VA) spectra; vibrational Raman optical activity (VROA) and corresponding vibrational Raman spectra; electronic circular dichroism (ECD) and corresponding electronic absorption (EA) spectra. In addition, the program allows for generating optical rotatory dispersion (ORD) as the Kramers–Kronig transform of ECD spectra. The simulation of theoretical spectra from transition strengths can be achieved using different bandshape profiles. The experimental and simulated theoretical spectra can be visually compared by displaying them together. A unique feature of CDSpecTech is performing spectral analysis using the ratio spectra; i.e., the dimensionless dissymmetry factor (DF) spectrum, which is the ratio of CD to absorption spectra, and the dimensionless circular intensity difference (CID) spectrum, which is the ratio of VROA to vibrational Raman spectra. The quantitative agreement between experimental and simulated theoretical spectra can also be assessed from the numerical similarity overlap between them. Two different similarity overlap methods are available. The program uses a graphical user interface which allows for ease of use and facilitates the analysis. All these features make CDSpecTech a valuable tool for the analysis of chiroptical spectra. The program is freely available on the World Wide Web.     

Electronic Circular Dichroism of [16]Helicene With Simplified TD‐DFT: Beyond the Single Structure Approach

The electronic circular dichroism (ECD) spectrum of the recently synthesized [16]helicene and a derivative comprising two triisopropylsilyloxy protection groups was computed by means of the very efficient simplified time‐dependent density functional theory (sTD‐DFT) approach. Different from many previous ECD studies of helicenes, nonequilibrium structure effects were accounted for by computing ECD spectra on "snapshots" obtained from a molecular dynamics (MD) simulation including solvent molecules. The trajectories are based on a molecule specific classical potential as obtained from the recently developed quantum chemically derived force field (QMDFF) scheme. The reduced computational cost in the MD simulation due to the use of the QMDFF (compared to ab‐initio MD) as well as the sTD‐DFT approach make realistic spectral simulations feasible for these compounds that comprise more than 100 atoms. While the ECD spectra of [16]helicene and its derivative computed vertically on the respective gas phase, equilibrium geometries show noticeable differences, these are “washed” out when nonequilibrium structures are taken into account. The computed spectra with two recommended density functionals (ωB97X and BHLYP) and extended basis sets compare very well with the experimental one. In addition we provide an estimate for the missing absolute intensities of the latter. The approach presented here could also be used in future studies to capture nonequilibrium effects, but also to systematically average ECD spectra over different conformations in more flexible molecules. Chirality 28:365–369, 2016. © 2016 Wiley Periodicals, Inc.     

Diffraction-based density restraints for membrane and membrane-peptide molecular dynamics simulations

We have recently shown that current molecular dynamics (MD) atomic force fields are not yet able to produce lipid bilayer structures that agree with experimentally-determined structures within experimental errors. Because of the many advantages offered by experimentally validated simulations, we have developed a novel restraint method for membrane MD simulations that uses experimental diffraction data. The restraints, introduced into the MD force field, act upon specified groups of atoms to restrain their mean positions and widths to values determined experimentally. The method was first tested using a simple liquid argon system, and then applied to a neat dioleoylphosphatidylcholine (DOPC) bilayer at 66% relative humidity and to the same bilayer containing the peptide melittin. Application of experiment-based restraints to the transbilayer double-bond and water distributions of neat DOPC bilayers led to distributions that agreed with the experimental values. Based upon the experimental structure, the restraints improved the simulated structure in some regions while introducing larger differences in others, as might be expected from imperfect force fields. For the DOPC-melittin system, the experimental transbilayer distribution of melittin was used as a restraint. The addition of the peptide caused perturbations of the simulated bilayer structure, but which were larger than observed experimentally. The melittin distribution of the simulation could be fit accurately to a Gaussian with parameters close to the observed ones, indicating that the restraints can be used to produce an ensemble of membrane-bound peptide conformations that are consistent with experiments. Such ensembles pave the way for understanding peptide-bilayer interactions at the atomic level.     

Exploring the energy landscape of protein folding using replica-exchange and conventional molecular dynamics simulations

Two independent replica-exchange molecular dynamics (REMD) simulations with an explicit water model were performed of the Trp-cage mini-protein. In the first REMD simulation, the replicas started from the native conformation, while in the second they started from a nonnative conformation. Initially, the first simulation yielded results qualitatively similar to those of two previously published REMD simulations: the protein appeared to be over-stabilized, with the predicted melting temperature 50-150K higher than the experimental value of 315K. However, as the first REMD simulation progressed, the protein unfolded at all temperatures. In our second REMD simulation, which starts from a nonnative conformation, there was no evidence of significant folding. Transitions from the unfolded to the folded state did not occur on the timescale of these simulations, despite the expected improvement in sampling of REMD over conventional molecular dynamics (MD) simulations. The combined 1.42 micros of simulation time was insufficient for REMD simulations with different starting structures to converge. Conventional MD simulations at a range of temperatures were also performed. In contrast to REMD, the conventional MD simulations provide an estimate of Tm in good agreement with experiment. Furthermore, the conventional MD is a fraction of the cost of REMD and continuous, realistic pathways of the unfolding process at atomic resolution are obtained.     

Molecular dynamics in cytochrome c oxidase Mössbauer spectra deconvolution

In this work low temperature molecular dynamics simulations of cytochrome c oxidase are used to predict an experimentally observable, namely Mössbauer spectra width. Predicted lineshapes are used to model Lorentzian doublets, with which published cytochrome c oxidase Mössbauer spectra were simulated. Molecular dynamics imposed constraints to spectral lineshapes permit to obtain useful information, like the presence of multiple chemical species in the binuclear center of cytochrome c oxidase. Moreover, a benchmark of quality for molecular dynamic simulations can be obtained. Despite the overwhelming importance of dynamics in electron–proton transfer systems, limited work has been devoted to unravel how much realistic are molecular dynamics simulations results. In this work, molecular dynamics based predictions are found to be in good agreement with published experimental spectra, showing that we can confidently rely on actual simulations. Molecular dynamics based deconvolution of Mössbauer spectra will lead to a renewed interest for application of this approach in bioenergetics.     

Conformation and dynamics of the three-helix bundle UBA domain of p62 from experiment and simulation

The ubiquitin associated domain of p62 is a small three-helix bundle of approximately 50 residues that mediates the recognition of polyubiquitin chains and ubiquitylated substrates. The solution structure of a 52 residue construct containing this domain has been characterized using heteronuclear nuclear magnetic resonance (NMR) methods. The resulting ensemble of NMR-derived structures was used in molecular dynamics (MD) simulations to investigate the equilibrium conformation and dynamics of this domain. NOE and (15)N relaxation data have been used to validate the structural ensemble produced by the MD simulations and show a good correlation for residues in regions of secondary structure. A similar approach was taken using an ensemble of structures from the MD simulations to calculate electronic circular dichroism (CD) and IR spectra from first principles with an encouraging correlation with the experimental CD and IR data.     

How does it really move? Recent progress in the investigation of protein nanosecond dynamics by NMR and simulation

Nuclear magnetic resonance (NMR) spin relaxation experiments currently probe molecular motions on timescales from picoseconds to nanoseconds. The detailed interpretation of these motions in atomic detail benefits from complementarity with the results from molecular dynamics (MD) simulations. In this mini-review, we describe the recent developments in experimental techniques to study the backbone dynamics from ¹⁵N relaxation and side-chain dynamics from ¹³C relaxation, discuss the different analysis approaches from model-free to dynamics detectors, and highlight the many ways that NMR relaxation experiments and MD simulations can be used together to improve the interpretation and gain insights into protein dynamics.     

Theoretical studies of the structure and molecular dynamics of a peptide crystal

Energy minimizations and molecular dynamics simulations have been performed on the cyclic peptide cyclo-(Ala-Pro-D-Phe)2 in both the isolated and crystal states. The results of these calculations have been analyzed, both to investigate our ability to reproduce experimental data (structure and vibrational and NMR spectra) and to investigate the effects of environment on the energy, structure, and dynamics of peptides. Comparison of the minimized and time-averaged crystal systems with the experimental peptide structure shows that the calculations have closely reproduced the experimental structure. Molecular dynamics of the isolated molecule has led to a new conformation, which is approximately equal to 8.5 kcal/mol more stable than the conformation that exists in the crystal, the latter conformation being stabilized by intermolecular (packing) forces. This illustrates the considerable effect that environment can have on the conformation of peptides. The crystal environment has also been shown to significantly reduce the dynamic conformational fluctuations seen for the isolated molecule. The behavior of the peptide during the isolated simulation also supports the experimental NMR observation of a symmetric structure that differs from the asymmetric, instantaneous structures which characterize the molecule during the dynamics. Calculations of vibrational frequencies of the peptide in the crystal and isolated states show the expected shifts in bond-stretching frequencies due to intermolecular interactions. Finally, we have calculated NMR coupling constants from the dynamics simulation of the isolated peptide and have compared these with the experimental values. This has led to a possible reinterpretation of the experimental data.