Tryptophan side chain conformers monitored by NMR and time-resolved fluorescence spectroscopies.

We have inserted a tryptophan (F77W) in the core of the regulatory domain of cardiac troponin C (cNTnC), and previously determined the structure of this mutant with and without the cosolvent trifluoroethanol (TFE). Interestingly, the orientations of the indole side chain of the Trp are in opposite directions in the two structures (Julien et al., Protein Sci 2009; 18:1165-1174). Fluorescence decay experiments for single Trp-containing proteins often show several lifetimes, which have been interpreted as reflecting conformational heterogeneity of the Trp side chain resulting from different rotamers. To test this interpretation, we monitored the effect of TFE on wild type, F77W and F77W-V82A calcium-saturated cNTnC using 2D (13)C-HSQC NMR and time-correlated single photon counting fluorescence spectroscopies. The time dependence of the Trp fluorescence decay was fit with three lifetimes. Addition of TFE caused a gradual, but limited decrease of the lifetimes due to dynamic quenching. For F77W cNTnC, the amplitude fractions of the lifetimes also changed upon addition of TFE-the long lifetime increased from 13 to 29%, while the middle lifetime decreased from 63 to 50% and the short lifetime remained relatively unchanged. For F77W-V82A cNTnC, comparable NMR changes are observed, confirming the switch in rotamer conformation, but only much smaller changes in fluorescence decay parameters were detected. These data indicate that the balance between the rotamer states can be changed without changing the lifetime amplitude fractions appreciably, and suggest that the environment(s) of the indole ring, responsible for the different lifetimes, can result from factors other than the intrinsic rotamer state of the tryptophan.

Tryptophan conformations associated with partial unfolding in ribonuclease T1.

The origin of the biexponential fluorescence decay of Trp in ribonuclease T1 under mildly destabilizing conditions, such as increased pH or temperature, or the presence of detergent, is still not understood. We have performed two extended replica-exchange molecular dynamics simulations to obtain a detailed representation of the native state at two protonation states corresponding to a high and low pH. At high pH, the appearance of partially unfolded states is evident. We found that this pH-induced destabilization originates from increased global repulsion as well as reduced local favorable electrostatic interactions and reduced H-bonding strength of His(27), His(40), and His(92). At high pH, alternative tryptophan rotamers appear and are linked to a distorted environment of the tryptophan, which also acts as a separate source of ground-state heterogeneity. The total population of these alternative conformations agrees well with the amplitude of the experimentally observed secondary fluorescence lifetime.

How do rotameric conformations influence the time-resolved fluorescence of tryptophan in proteins? A perspective based on molecular modeling and quantum chemistry.

We discuss the dynamics of tryptophan rotamers in the context of the non-exponential fluorescence decay in proteins. The central question is: how does the ground-state conformational heterogeneity influence the time evolution of tryptophan fluorescence? This problem is examined here from the theoretical perspective. Three methods at different levels of theory, and with different scopes and computational requirements are reviewed. The Dead-end elimination method is limited to side-chain dynamics and provides an efficient way to detect the stable tryptophan rotamers in a protein. Its application to the study of heterogeneous emission characteristics is illustrated. Molecular dynamics is aimed at the full phase space of the macromolecule in solution, but must rely on classical force fields and laws of evolution. We examine to what extent the molecular mechanics paradigm yields sufficiently accurate thermodynamic results, and what are the possible kinetic implications. Finally Quantum Chemistry is the only theoretical method that allows a direct assessment of the excited states. It is necessarily restricted to small molecular systems, and thus must be used in a hybrid combination with classical methods and electrostatic models. So far understanding of the emitting state has greatly progressed as a result of these calculations, but the actual treatment of the photophysical decay processes at the quantum level has not yet really started.

Exploration of the activation pathway of Deltaalpha-Chymotrypsin with molecular dynamics simulations and correlation with kinetic experiments.

Correlating the experimentally observed kinetics of protein conformational changes with theoretical predictions is a formidable and challenging task, due to the multitude of degrees of freedom (>5,000) in a protein and the huge gap between the timescale of the kinetic event of interest (ms) and the typical timescale of computer simulations (ns). In this study we show that using the targeted molecular dynamics (TMD) method it is possible to simulate conformational changes of the ms time range and to correlate multiple simulations of single pathways with ensemble experiments on both the structural and energetic basis. As a model system we chose to study the conformational change of rat-Deltaalpha-chymotrypsin from its inactive to its active conformation. This activation process has been analyzed previously by experimental and theoretical methods, i.e. fluorescence stopped-flow spectroscopy (FSF), molecular dynamics (MD) and TMD. Inspired by the results of these studies on the wild type (WT) enzyme, several mutants were constructed to alter the conformational pathway and studied by FSF measurements. In the present work WT and mutant N18G were subjected to multiple MD and subsequent TMD simulations. We report the existence of two main activation pathways, a feature of chymotrypsin activation that has been hitherto unknown. A method to correlate the energetics of the different pathways calculated by TMD and the kinetic parameters observed by experimental methods such as FSF is presented. Our work is relevant for experimental single molecule studies of enzymes in general.

Impact on replicative fitness of the G48E substitution in the protease of HIV-1: an in vitro and in silico evaluation.

We observed an unusual glycine-to-glutamate substitution at protease (PR) residue position 48 (G48E) in an African patient infected with a subtype A1 HIV-1 strain failing a saquinavir-containing regimen. Phenotypic analysis of protease inhibitor (PI) susceptibility showed that the G48E site-directed mutant, when introduced into an NL4-3 HIV-1 PR backbone, was slightly resistant to SQV (2-fold when compared with the wild-type virus). In addition, the G48E and G48E/V82A site-directed mutants were associated with a decrease in fitness, whereas a reversion to the wild type at position 48 was observed in vitro. Growth competition experiments using a novel growth competition assay based on enhanced green fluorescent protein- or Discosoma spp. red fluorescent protein-expressing viruses showed that the replicative fitness of the G48E virus was reduced to 55% compared with the parental NL4-3 virus. Synthesizing all possible site-directed mutants found in the patient strain is too time-consuming; therefore, a molecular dynamics (MD) simulation approach was used to understand why this mutation survived despite its fitness cost. These simulations documented that the G48E mutant interacted with PI resistance mutations (M46I, I54V, Q58E, and L63P) and with natural polymorphisms specific to subtype A1 (E35D, M36I, and R57K) that were present in the patient's virus. We hypothesize that the polymorphisms contained in the PR flap regions of the patient's virus may compensate for the presence of G48E, possibly by restoring the flexibility of the PR flaps. In summary, our results demonstrate that the G48E substitution, when introduced in the context of an HIV-1 subtype B strain, is highly unstable and gives rise to viruses with a poor replicative fitness in vitro. We also showed that when confronted with too many mutations to evaluate in vitro, MD simulations are helpful to draft hypotheses on how polymorphisms can interact with resistance mutations to stabilize their potential fitness cost.

An unusual red-edge excitation and time-dependent Stokes shift in the single tryptophan mutant protein DD-carboxypeptidase from Streptomyces: the role of dynamics and tryptophan rotamers.

The fluorescence emission of the single tryptophan (W233) of the mutant protein DD-carboxypeptidase from streptomyces is characterized by a red-edge excitation shift (REES), i.e., the phenomenon that the wavelength of maximum emission depends on the excitation wavelength. This phenomenon is an indication for a strongly reduced dynamic environment of the single tryptophan, which has a very low accessibility to the solvent. The REES shows, however, an unusual temperature and time dependence. This, together with the fluorescence lifetime analysis, showing three resolvable lifetimes, can be explained by the presence of three rotameric states that can be identified using the Dead-End Elimination method. The three individual lifetimes increase with increasing emission wavelength, indicating the presence of restricted protein dynamics within the rotameric states. This is confirmed by time-resolved anisotropy measurements that show dynamics within the rotamers but not among the rotamers. The global picture is that of a protein with a single buried tryptophan showing strongly restricted dynamics within three distinct rotameric states with different emission spectra and an anisotropic environment.

An alternate sucrose binding mode in the E203Q Arabidopsis invertase mutant: an X-ray crystallography and docking study.

In the present study, we report on the X-ray crystallographic structure of a GH32 invertase mutant, (i.e., the Arabidopsis thaliana cell-wall invertase 1-E203Q, AtcwINV1-mutant) in complex with sucrose. This structure was solved to reveal the features of sugar binding in the catalytic pocket. However, as demonstrated by the X-ray structure the sugar binding and the catalytic pocket arrangement is significantly altered as compared with what was expected based on previous X-ray structures on GH-J clan enzymes. We performed a series of docking and molecular dynamics simulations on various derivatives of AtcwINV1 to reveal the reasons behind this modified sugar binding. Our results demonstrate that the E203Q mutation introduced into the catalytic pocket triggers conformational changes that alter the wild type substrate binding. In addition, this study also reveals the putative productive sucrose binding modus in the wild type enzyme.