Gouri S. Jas

The complete folding pathway of a protein from nanoseconds to microseconds

Combining experimental and simulation data to describe all of the structures and the pathways involved in folding a protein is problematical. Transition states can be mapped experimentally by φ values, but the denatured state is very difficult to analyse under conditions that favour folding. Also computer simulation at atomic resolution is currently limited to about a microsecond or less. Ultrafast-folding proteins fold and unfold on timescales accessible by both approaches, so here we study the folding pathway of the three-helix bundle protein Engrailed homeodomain. Experimentally, the protein collapses in a microsecond to give an intermediate with much native α-helical secondary structure, which is the major component of the denatured state under conditions that favour folding. A mutant protein shows this state to be compact and contain dynamic, native-like helices with unstructured side chains. In the transition state between this and the native state, the structure of the helices is nearly fully formed and their docking is in progress, approximating to a classical diffusion–collision model. Molecular dynamics simulations give rate constants and structural details highly consistent with experiment, thereby completing the description of folding at atomic resolution.

The Structure of DNA within Cationic Lipid/DNA Complexes

The structure of DNA within CLDCs used for gene delivery is controversial. Previous studies using CD have been interpreted to indicate that the DNA is converted from normal B to C form in complexes. This investigation reexamines this interpretation using CD of model complexes, FTIR as well as Raman spectroscopy and molecular dynamics simulations to address this issue. CD spectra of supercoiled plasmid DNA undergo a significant loss of rotational strength in the signal near 275 nm upon interaction with either the cationic lipid dimethyldioctadecylammonium bromide or 1,2-dioleoyltrimethylammonium propane. This loss of rotational strength is shown, however, by both FTIR and Raman spectroscopy to occur within the parameters of the B-type conformation. Contributions of absorption flattening and differential scattering to the CD spectra of complexes are unable to account for the observed spectra. Model studies of the CD of complexes prepared from synthetic oligonucleotides of varying length suggest that significant reductions in rotational strength can occur within short stretches of DNA. Furthermore, some alteration in the hydrogen bonding of bases within CLDCs is indicated in the FTIR and Raman spectroscopy results. In addition, alterations in base stacking interactions as well as hydrogen bonding are suggested by molecular dynamics simulations. A global interpretation of all of the data suggests the DNA component of CLDCs remains in a variant B form in which base/base interactions are perturbed.

Ab initio calculations of S1 excited state vibrational spectra of benzene, naphthalene and anthracene

Abstract Normal-mode calculations are presented for the lowest singlet excited states S 1 of benzene, naphthalene and anthracene. Optimized geometries and Cartesian harmonic force constants of the excited states are obtained from ab initio calculations at the CIS/6-31G, CIS/6-31G ∗ , CIS/6-311G and CIS/6-311G ∗ levels for benzene, and at CIS/6-31G and CIS/6-31G ∗ levels for naphthalene and anthracene. Normal-mode analysis is performed in internal coordinates, yielding vibrational frequencies and forms of normal modes for the parent molecules and their perdeuterated derivatives. The results are compared with corresponding properties of the ground state S 0 calculated at the Hartree-Fock level and with available experimental data. The overall changes in molecular geometry and vibrational spectra upon S 0 → S 1 excitation are small. For benzene we find excellent agreement of the entire calculated S 1 vibrational spectrum with experimental results. For naphthalene and anthracene the calculated vibrational frequencies are in good agreement with the limited available experimental data. Our calculations provide predictions of frequencies and types of normal modes for the complete S 1 vibrational spectra of these molecules, and suggest reassignments of several normal-mode frequencies compared to the limited experimental results that are currently available.

Influence of temperature and viscosity on anthracene rotational diffusion in organic solvents: Molecular dynamics simulations and fluorescence anisotropy study

Molecular dynamics simulations and fluorescence anisotropy decay measurements are used to investigate the rotational diffusion of anthracene in two organic solvents—cyclohexane and 2-propanol—at several temperatures. Molecular dynamics simulations of 1 ns length were performed for anthracene in cyclohexane (at 280, 296, and 310 K) and in 2-propanol (at 296 K). The calculated time constants for reorientation of the short in-plane axis were 7–9 and 11–16 ps at 296 K in cyclohexane and 2-propanol, respectively, in excellent agreement with corresponding fluorescence depolarization measurements of 8 and 14 ps. The measured rotational reorientation times and the calculated average rotational diffusion coefficients varied in accord with Debye–Stokes–Einstein theory. Their magnitudes were close to values predicted for an ellipsoid of shape and size equivalent to an anthracene molecule, and exhibited predictable variation with external conditions—increasing with temperature and decreasing with solvent viscosity. H...

Unassisted Transport of N-acetyl-L-tryptophanamide through Membrane: Experiment and Simulation of Kinetics

Cellular transport machinery, such as channels and pumps, is working against the background of unassisted material transport through membranes. The permeation of a blocked tryptophan through a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) membrane is investigated to probe unassisted or physical transport. The transport rate is measured experimentally and modeled computationally. The time scale measured by parallel artificial membrane permeation assay (PAMPA) experiments is ~8 h. Simulations with the milestoning algorithm suggest mean first passage time (MFPT) of ~4 h and the presence of a large barrier at the center of the bilayer. A similar calculation with the solubility-diffusion model yields a MFPT of ~15 min. This permeation rate is 9 orders of magnitude slower than the permeation rate of only a tryptophan side chain (computed by us and others). This difference suggests critical dependence of transport time on permeant size and hydrophilicity. Analysis of the simulation results suggests that the permeant partially preserves hydrogen bonding of the peptide backbone to water and lipid molecules even when it is moving closer to the bilayer center. As a consequence, defects of the membrane structure are developed to assist permeation.

Structure and Dynamics of Calcium-activated Calmodulin in Solution

Abstract Two 4-ns molecular dynamics simulations of calcium loaded calmodulin in solution have been performed, using both standard nonbonded cutoffs and Ewald summation to treat electrostatic interactions. Our simulation results are generally consistent with solution experimental studies of calmodulin structure and dynamics, including NMR, cross-linking, fluorescence and x-ray scattering. The most interesting result of the molecular dynamics simulations is the detection of large-scale structural fluctuations of calmodulin in solution. The globular N- and C-terminal domains tend to move approximately like rigid bodies, with fluctuations of interdomain distances within a 7 Å range and of interdomain angles by up to 60 deg. Essential dynamics analysis indicates that the three dominant types of motion involve bending of the central helix in two perpendicular planes and a twist in which the domains rotate in opposite directions around the central helix. In the more realistic Ewald trajectory the protein backbone remains mostly within a 2–3 Å root-mean-square distance from the crystal structure, the secondary structure within the domains is conserved and middle part of the central helix becomes disordered. The central helix itself exhibits limited fluctuations, with its bend angle exploring the 0–50° range and the end-to-end distance falling in 39–43 Å. The results of the two simulations were similar in many respects. However, the cutoff trajectory exhibited a larger deviation from the crystal, loss of several helical hydrogen bonds in the N-terminal domain and lack of structural disorder in the central helix.

Helix Formation in a Pentapeptide: Experiment and Force-field Dependent Dynamics

We used a combined approach of experiment and simulation to determine the helical population and folding pathway of a small helix forming blocked pentapeptide, Ac-(Ala)(5)-NH(2). Experimental structural characterization of this blocked peptide was carried out with far UV circular dichroism spectroscopy, FTIR, and NMR measurements. These measurements confirm the presence of the α-helical state in a buffer solution. Direct molecular dynamics and replica-exchange simulations of the pentapeptide were performed using several popular force fields with explicit solvent. The simulations yielded statistically reliable estimates of helix populations, melting curves, folding, and nucleation times. The distributions of conformer populations are used to measure folding cooperativity. Finally, a statistical analysis of the sample of helix-coil transition paths was performed. The details of the calculated helix populations, folding kinetics and pathways vary with the employed force field. Interestingly, the helix populations, folding, and unfolding times obtained from most of the studied force fields are in qualitative agreement with each other and with available experimental data, with the deviations corresponding to several kcal/mol in energy at 300 K. Most of the force fields also predict qualitatively similar transition paths, with unfolding initiated at the C-terminus. Accuracy of potential energy parameters, rather than conformational sampling may be the limiting factor in current molecular simulations.

Experiments and comprehensive simulations of the formation of a helical turn

We investigate the kinetics and thermodynamics of a helical turn formation in the peptide Ac-WAAAH-NH(2). NMR measurements indicate that this peptide has significant tendency to form a structure of a helical turn, while temperature dependent CD establishes the helix fraction at different temperatures. Molecular dynamics and milestoning simulations agree with experimental observables and suggest an atomically detailed picture for the turn formation. Using a network representation, two alternative mechanisms of folding are identified: (i) a direct co-operative mechanism from the unfolded to the folded state without intermediate formation of hydrogen bonds and (ii) an indirect mechanism with structural intermediates with two residues in a helical conformation. This picture is consistent with kinetic measurements that reveal two experimental time scales of sub-nanosecond and several nanoseconds.

Reorientations of Aromatic Amino Acids and Their Side Chain Models: Anisotropy Measurements and Molecular Dynamics Simulations

We present a study of reorientation dynamics of the three aromatic amino acids and their side chain models in aqueous solution. Experimentally, anisotropy decay measurements with picosecond time resolution were performed for blocked tryptophan, tyrosine, and phenyalanine and model compounds p-cresol and 3-methylindole. Computationally, rotational diffusion was modeled by molecular dynamics simulations for the three aromatic residues and their side chain models: benzene, toluene, phenol, p-cresol, indole, and 3-methylindole in explicit water. Our simulations used the CHARMM protein force field and associated TIP3P water model and tend to underestimate the rotational correlation times. However, the simulations yield several interesting qualitative insights into reorientational motions that complement the experimental measurements. The effects of substituent and temperature on reorientations of the parent compounds are well reproduced computationally. Additionally, simulations indicate strongly anisotropic reorientations for most of the studied compounds and a separation of time scales between conformational dynamics and rotational diffusion. Comparison with continuum hydrodynamic models suggests that we may consider that the blocked amino acids move under stick boundary conditions, while the dynamics for most of the model compounds falls between stick and slip conditions. Our systematic treatment of blocked amino acids, starting from the parent compounds (benzene, phenol, and indole) provides a baseline for understanding the anisotropy decay signals of more complicated peptide systems.

Free-energy simulations of the oxidation of C-terminal methionines in calmodulin

In the course of aging or under conditions of oxidative stress, methionine residues of calmodulin undergo oxidation, leading to loss of biological activity of the protein. We have performed free‐energy simulations of the effects of C‐terminal methionine side‐chain oxidation on the properties of calmodulin. The simulation results indicate that oxidation should have a destabilizing effect on all three protein functional states: calcium free, calcium loaded, and with both calcium and target peptide bound. Because the different states are destabilized by different amounts, this leads to a more complex pattern in the observable effects on protein thermal stability, calcium affinity, and binding of a target peptide. The influence of oxidation on the free energy of CaM unfolding is estimated by comparing the free‐energy cost of oxidizing a Met residue in a Gly‐Met‐Gly peptide and in the protein. The protein thermal stability of the oxidized forms is lowered by a moderate amount 1–3 kcal/mol, in qualitative agreement with experimental results of 0.3 kcal/mol. The calculated changes in affinity for calcium and for the target peptide show opposing trends. Oxidation at position 144 is predicted to enhance peptide binding and weaken calcium binding, whereas oxidation at 145 weakens peptide binding and enhances affinity for calcium. The lower affinity of Met 145‐oxidized calmodulin toward the target peptide correlates with experimentally observed lowering of calmodulin‐activated Ca‐ATPase activity when oxidized calmodulin from aged rat brains is used. Thus, our simulations suggest that Met 145 is the oxidation site in the C‐terminal fragment of calmodulin. The microscopic mechanism behind the calculated free energy changes appears to be a greater affinity for water of the oxidized Met side‐chain relative to normal Met. Structures with Met exposed to solvent had consistently lower free energies than those with buried Met sidechains. Proteins 2002;48:257–268.