A mean field model of ligand-protein interactions: implications for the structural assessment of human immunodeficiency virus type 1 protease complexes and receptor-specific binding.

We propose a general mean field model of ligand-protein interactions to determine the thermodynamic equilibrium of a system at finite temperature. The method is employed in structural assessments of two human immuno-deficiency virus type 1 protease complexes where the gross effects of protein flexibility are incorporated by utilizing a data base of crystal structures. Analysis of the energy spectra for these complexes has revealed that structural and thermo-dynamic aspects of molecular recognition can be rationalized on the basis of the extent of frustration in the binding energy landscape. In particular, the relationship between receptor-specific binding of these ligands to human immunodeficiency virus type 1 protease and a minimal frustration principle is analyzed.

Multiple pathways for ultrafast transduction of light energy in the photosynthetic reaction center of Rhodobacter sphaeroides

A pathway of electron transfer is described that operates in the wild-type reaction center (RC) of the photosynthetic bacterium Rhodobacter sphaeroides. The pathway does not involve the excited state of the special pair dimer of bacteriochlorophylls (P*), but instead is driven by the excited state of the monomeric bacteriochlorophyll (BA*) present in the active branch of pigments along which electron transfer occurs. Pump-probe experiments were performed at 77 K on membrane-bound RCs by using different excitation wavelengths, to investigate the formation of the charge separated state P+HA−. In experiments in which P or BA was selectively excited at 880 nm or 796 nm, respectively, the formation of P+HA− was associated with similar time constants of 1.5 ps and 1.7 ps. However, the spectral changes associated with the two time constants are very different. Global analysis of the transient spectra shows that a mixture of P+BA− and P* is formed in parallel from BA* on a subpicosecond time scale. In contrast, excitation of the inactive branch monomeric bacteriochlorophyll (BB) and the high exciton component of P (P+) resulted in electron transfer only after relaxation to P*. The multiple pathways for primary electron transfer in the bacterial RC are discussed with regard to the mechanism of charge separation in the RC of photosystem II from higher plants.

Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

Probing the interaction between two single molecules: fluorescence resonance energy transfer between a single donor and a single acceptor.

We extend the sensitivity of fluorescence resonance energy transfer (FRET) to the single molecule level by measuring energy transfer between a single donor fluorophore and a single acceptor fluorophore. Near-field scanning optical microscopy (NSOM) is used to obtain simultaneous dual color images and emission spectra from donor and acceptor fluorophores linked by a short DNA molecule. Photodestruction dynamics of the donor or acceptor are used to determine the presence and efficiency of energy transfer. The classical equations used to measure energy transfer on ensembles of fluorophores are modified for single-molecule measurements. In contrast to ensemble measurements, dynamic events on a molecular scale are observable in single pair FRET measurements because they are not canceled out by random averaging. Monitoring conformational changes, such as rotations and distance changes on a nanometer scale, within single biological macromolecules, may be possible with single pair FRET.

The CuA center of cytochrome-c oxidase: electronic structure and spectra of models compared to the properties of CuA domains.

The electronic structure and spectrum of several models of the binuclear metal site in soluble CuA domains of cytochrome-c oxidase have been calculated by the use of an extended version of the complete neglect of differential overlap/spectroscopic method. The experimental spectra have two strong transitions of nearly equal intensity around 500 nm and a near-IR transition close to 800 nm. The model that best reproduces these features consists of a dimer of two blue (type 1) copper centers, in which each Cu atom replaces the missing imidazole on the other Cu atom. Thus, both Cu atoms have one cysteine sulfur atom and one imidazole nitrogen atom as ligands, and there are no bridging ligands but a direct Cu-Cu bond. According to the calculations, the two strong bands in the visible region originate from exciton coupling of the dipoles of the two copper monomers, and the near-IR band is a charge-transfer transition between the two Cu atoms. The known amino acid sequence has been used to construct a molecular model of the CuA site by the use of a template and energy minimization. In this model, the two ligand cysteine residues are in one turn of an alpha-helix, whereas one ligand histidine is in a loop following this helix and the other one is in a beta-strand.