Dispersed polaron simulations of electron transfer in photosynthetic reaction centers.

A microscopic method for simulating quantum mechanical, nuclear tunneling effects in biological electron transfer reactions is presented and applied to several electron transfer steps in photosynthetic bacterial reaction centers. In this "dispersed polaron" method the fluctuations of the protein and the electron carriers are projected as effective normal modes onto an appropriate reaction coordinate and used to evaluate the quantum mechanical rate constant. The simulations, based on the crystallographic structure of the reaction center from Rhodopseudomonas viridis, focus on electron transfer from a bacteriopheophytin to a quinone and the subsequent back-reaction. The rates of both of these reactions are almost independent of temperature or even increase with decreasing temperature. The simulations reproduce this unusual temperature dependence in a qualitative way, without the use of adjustable parameters for the protein's Franck-Condon factors. The observed dependence of the back-reaction on the free energy of the reaction also is reproduced, including the special behavior in the "inverted region."

Electrostatic control of charge separation in bacterial photosynthesis.

Electrostatic interaction energies of the electron carriers with their surroundings in a photosynthetic bacterial reaction center are calculated. The calculations are based on the detailed crystal structure of reaction centers from Rhodopseu-domonas viridis, and use an iterative, self-consistent procedure to evaluate the effects of induced dipoles in the protein and the surrounding membrane. To obtain the free energies of radical-pair states, the calculated electrostatic interaction energies are combined with the experimentally measured midpoint redox potentials of the electron carriers and of bacteriochlorophyll (BChl) and bacteriopheophytin (BPh) in vitro. The P+HL- radical-pair, in which an electron has moved from the primary electron donor (P) to a BPh on the 'L' side of the reaction center (HL), is found to lie approx. 2.0 kcal/mol below the lowest excited singlet state (P*), when the radical-pair is formed in the static crystallographic structure. The reorganization energy for the subsequent relaxation of P+HL- is calculated to be 5.0 kcal/mol, so that the relaxed radical-pair lies about 7 kcal/mol below P*. The unrelaxed P+BL- radical-pair, in which the electron acceptor is the accessory BChl located between P and HL, appears to be essentially isoenergetic with P*.P+BM-, in which an electron moves to the BChl on the 'M' side, is calculated to lie about 5.5 kcal/mol above P*. These results have an estimated error range of +/- 2.5 kcal/mol. They are shown to be relatively insensitive to various details of the model, including the charge distribution in P+, the atomic charges used for the amino acid residues, the boundaries of the structural region that is considered microscopically and the treatments of the histidyl ligands of P and of potentially ionizable amino acids. The calculated free energies are consistent with rapid electron transfer from P* to HL by way of BL, and with a much slower electron transfer to the pigments on the M side. Tyrosine M208 appears to play a particularly important role in lowering the energy of P+BL-. Electrostatic interactions with the protein favor localization of the positive charge of P+ on PM, one of the two BChl molecules that make up the electron donor.

[Genetic relationship among four Yunnan populations revealed by sequence of mtDNA D-loop].

mtDNA D-loop noncoding region 16048-16569 and the following 1-41 (563 bp) in 99 individuals of four Yunnan ethnic minorities (Dai, Wa, Lahu and Tibetan) were sequenced and then a phylogenetic tree was reconstructed by Neighbor-Joining method. These 99 mtDNA lineages were classified into 3 genotype groups in the tree. All lineages with 9 bp deletion in the COII/tRAN(Lys) intergenic region were clustered in group I, some individuals of Dai, Lahu, Wa and only 2 Tibetan individuals clustered in group II, individuals of all four populations were included in group III. A phylogenetic tree of the four populations was constructed by NJ method on the basis of estimate of net genetic distance among them. Our results showed, the genetic distance among Dai, Wa and Lahu is very close, but far from Tibetan, their genetic distance is similar to their geographic distance. Although both as descendants of ancient Di-Qiang tribe in history and speaking similar language, Lahu and Tibetan are not closely related. This result indicates that there are different origins of these two populations.

[Genetic relationships among six Chinese populations revealed by analysis of 30 autosomal STRs].

30 autosomal STRs of 6 Chinese populations (Bai, Naxi, Tu, Sala, Han in Shandong, She) were amplified by multiplex PCRs using fluerescein-labelled primers. Shriver's Dsw was estimated on the basis of the results of the genescanning and genotyping after running unnatural PAGE of the PCRs' products on ABI 377 sequencer. Phylogenetic trees were constructed by using Neighbor-Joining and UPGMA method based on Dsw, and then the genetic relationships among them were analyzed referring to some relative informations. Our results indicated that the genetic distance between Sala and Tu is near, 0.033. But the distances between Sala and other four populations are far, over 0.12; Tu is close to Naxi and Shandong Han, and the distances are 0.038, 0.063 respectively; The distance between Bai and Han is the nearest, 0.007, but there is a distance, 0.075, between Bai and Naxi, and a far distance, 0.112, between Bai and Tu; The distance between Naxi and Han is 0.100 and the distances between She and other 5 populations are all over 0.12. In both of the NJ and UPGMA phylogenetic trees, Naxi, Tu and Sala is one cluster and Bai and Han is another cluster. She is a single branch. These results, basically consistent with 6 Chinese populations' geographic distribution and histories, can provide some genetic information to comprehensively study their origin, migration, formation and development with their historical records and archaeological evidence.

Reorganization energy of the initial electron-transfer step in photosynthetic bacterial reaction centers.

The reorganization energy (lambda) for electron transfer from the primary electron donor (P*) to the adjacent bacteriochlorophyll (B) in photosynthetic bacterial reaction centers is explored by molecular-dynamics simulations. Relatively long (40 ps) molecular-dynamics trajectories are used, rather than free energy perturbation techniques. When the surroundings of the reaction center are modeled as a membrane, lambda for P* B --> P+ B- is found to be approximately 1.6 kcal/mol. The results are not sensitive to the treatment of the protein's ionizable groups, but surrounding the reaction center with water gives higher values of lambda (approximately 6.5 kcal/mol). In light of the evidence that P+ B- lies slightly below P* in energy, the small lambda obtained with the membrane model is consistent with the speed and temperature independence of photochemical charge separation. The calculated reorganization energy is smaller than would be expected if the molecular dynamics trajectories had sampled the full conformational space of the system. Because the system does not relax completely on the time scale of electron transfer, the lambda obtained here probably is more pertinent than the larger value that would be obtained for a fully equilibrated system.

Microscopic simulation of quantum dynamics and nuclear tunneling in bacterial reaction centers.

The "dispersed polaron" version of the semiclassical trajectory approach is used to evaluate the quantum mechanical nuclear tunneling effects in the charge recombination reaction, P(+)Q(-)→PQ, in photosynthetic bacterial reaction centers, The cclculations are based on the crystallographic structure of reaction centers from Rhodopseudomonas viridis. They succeed in capturing the temperature dependence of the rate constant without using adjustable parameters. This provides the first example of a microscopic simulation of quantum mechanical nuclear tunneling in a biological system.

Calculations of antibody-antigen interactions: microscopic and semi-microscopic evaluation of the free energies of binding of phosphorylcholine analogs to McPC603.

The study of antibody-antigen interactions should greatly benefit from the development of quantitative models for the evaluation of binding free energies in proteins. The present work addresses this challenge by considering the test case of the binding free energies of phosphorylcholine analogs to the murine myeloma protein McPC603. This includes the evaluation of the differential binding energy as well as the absolute binding energies and their corresponding electrostatic contributions. Four different approaches are examined: the Protein Dipoles Langevin Dipoles (PDLD) method, the semi-microscopic PDLD (PDLD/S) method, a free energy perturbation (FEP) method based on an adiabatic charging procedure and a linear response approximation that accelerates the FEP calculation. The PDLD electrostatic calculations are augmented by estimates of the relevant hydrophobic and steric contributions. The determination of the hydrophobic energy involves an approach which considers the modification of the effective surface area of the solute by local field effects. The steric contributions are analyzed in terms of the corresponding reorganization energies. This treatment, which considers the protein as a harmonic system, views the steric forces as the restoring forces for the electrostatic interactions. The FEP method is found to give unreliable results with regular cut-off radii and starts to give quantitative results only in very expensive treatment with very large cut-off radii. The PDLD and PDLD/S methods are much faster than the FEP approach and give reasonable results for both the relative and absolute binding energies. The speed and simplicity of the PDLD/S method make it an effective strategy for interactive docking studies and indeed such an option is incorporated in the program MOLARIS. A component analysis of the different energy contributions of the FEP treatment and a similar PDLD analysis indicate that electrostatic effects provide the largest contribution to the differential binding energy, while the hydrophobic and steric contributions are much smaller. This finding lends further support to the idea that electrostatic interactions play a major role in determining the antigen specificity of McPC603.

Examining methods for calculations of binding free energies: LRA, LIE, PDLD-LRA, and PDLD/S-LRA calculations of ligands binding to an HIV protease.

Several strategies for evaluation of the protein-ligand binding free energies are examined. Particular emphasis is placed on the Linear Response Approximation (LRA) (Lee et. al., Prot Eng 1992;5:215-228) and the Linear Interaction Energy (LIE) method (Aqvist et. al., Prot Eng 1994;7:385-391). The performance of the Protein Dipoles Langevin Dipoles (PDLD) method and its semi-microscopic version (the PDLD/S method) is also considered. The examination is done by using these methods in the evaluating of the binding free energies of neutral C2-symmetric cyclic urea-based molecules to Human Immunodeficiency Virus (HIV) protease. Our starting point is the introduction of a thermodynamic cycle that decomposes the total binding free energy to electrostatic and non-electrostatic contributions. This cycle is closely related to the cycle introduced in our original LRA study (Lee et. al., Prot Eng 1992;5:215-228). The electrostatic contribution is evaluated within the LRA formulation by averaging the protein-ligand (and/or solvent-ligand) electrostatic energy over trajectories that are propagated on the potentials of both the polar and non-polar (where all residual charges are set to zero) states of the ligand. This average involves a scaling factor of 0.5 for the contributions from each state and this factor is being used in both the LRA and LIE methods. The difference is, however, that the LIE method neglects the contribution from trajectories over the potential of the non-polar state. This approximation is entirely valid in studies of ligands in water but not necessarily in active sites of proteins. It is found in the present case that the contribution from the non-polar states to the protein-ligand binding energy is rather small. Nevertheless, it is clearly expected that this term is not negligible in cases where the protein provides preorganized environment to stabilize the residual charges of the ligand. This contribution can be particularly important in cases of charged ligands. The analysis of the non-electrostatic term is much more complex. It is concluded that within the LRA method one has to complete the relevant thermodynamic cycle by evaluating the binding free energy of the "non-polar" ligand, l;, where all the residual charges are set to zero. It is shown that the LIE term, which involves the scaling of the van der Waals interaction by a constant beta (usually in the order of 0.15 to 0.25), corresponds to this part of the cycle. In order to elucidate the nature of this non-electrostatic term and the origin of the scaling constant beta, it is important to evaluate explicitly the different contributions to the binding energy of the non-polar ligand, DeltaG(bind,l;). Since this cannot be done at present (for relatively large ligands) by rigorous free energy perturbation approaches, we evaluate DeltaG(bind,l;) by the PDLD approach, augmented by microscopic calculations of the change in configurational entropy upon binding. This evaluation takes into account the van der Waals, hydrophobic, water penetration and entropic contributions, which are the most important free energy contributions that make up the total DeltaG(bind,l;). The sum of these contributions is scaled by a factor straight theta and it is argued that obtaining a quantitative balance between these contributions should result in straight theta = 1. By doing so we should have a reliable estimate of the value of the LIE beta and a way to understand its origin. The present approach gives straight theta values between 0.5 and 0.73, depending on the approximation used. This is encouraging but still not satisfying. Nevertheless, one might be able to use our PDLD approach to estimate the change of the LIE straight theta between different protein active sites. It is pointed out that the LIE method is quite similar to our original approach where the electrostatic term was evaluated by the LRA method and the non-electrostatic term by the PDLD method (with its vdw, solvation,

The minimal polymorphism of class II E alpha chains is not due to the functional neutrality of mutations.

Given the extensive allelic amino acid sequence polymorphism present in the first domain of A alpha, A beta, and E beta chains and its profound effects on class II function, the minimal polymorphism in the mouse E alpha chain (and in its human homologue DR alpha) is paradox. Two possible explanations for the lack of polymorphism in E alpha are: (1) the E alpha chain plays such a uniquely critical structural/functional role in antigen presentation, T-cell activation, repertoire selection, and/or pairing with E beta or other proteins for expression that it cannot vary, and mutations are selected against; (2) the E alpha chain plays a less significant role than the outer domains of other major histocompatibility complex (MHC) proteins in determining the interactions with processed peptides or with T-cell receptor (TCR), so there is no selective pressure to maintain new mutations. To explore this question we compared the ability of transfectants expressing wild type (wt) E alpha E beta d and mutant E alpha wt E beta d proteins to present peptides and bacterial superantigens to T-cell hybridomas. Mutations at the E alpha amino acid positions 31, 52, and 65&66, to residues that represent allelic alternatives in A alpha chains, significantly reduced activation of peptide-specific T hybridomas, and mutations at 71 sometimes enhanced T-cell stimulation. None of the E alpha mutations reduced, and some enhanced, superantigen stimulation of T-cell hybridomas. These results argue against the hypothesis that E alpha chains are minimally polymorphic because mutations in E alpha are functionally neutral.

How important are entropic contributions to enzyme catalysis?

The idea that enzymes accelerate their reactions by entropic effects has played a major role in many prominent proposals about the origin of enzyme catalysis. This idea implies that the binding to an enzyme active site freezes the motion of the reacting fragments and eliminates their entropic contributions, (delta S(cat)(double dagger))', to the activation energy. It is also implied that the binding entropy is equal to the activation entropy, (delta S(w)(double dagger))', of the corresponding solution reaction. It is, however, difficult to examine this idea by experimental approaches. The present paper defines the entropic proposal in a rigorous way and develops a computer simulation approach that determines (delta S(double dagger))'. This approach allows us to evaluate the differences between (delta S(double dagger))' of an enzymatic reaction and of the corresponding reference reaction in solution. Our approach is used in a study of the entropic contribution to the catalytic reaction of subtilisin. It is found that this contribution is much smaller than previously thought. This result is due to the following: (i) Many of the motions that are free in the reactants state of the reference solution reaction are also free at the transition state. (ii) The binding to the enzyme does not completely freeze the motion of the reacting fragments so that (delta S(double dagger))' in the enzymes is not zero. (iii) The binding entropy is not necessarily equal to (delta S(w)(double dagger))'.

Mitochondrial DNA polymorphisms in Yunnan nationalities in China.

Nucleotide sequences of the D-loop region of human mitochondrial DNA from four Yunnan nationalities, Dai, Wa, Lahu, and Tibetan, were analyzed. Based on a comparison of 563-bp sequences in 99 people, 66 different sequence types were observed. Of these, 64 were unique to their respective populations, whereas only 2 types were shared between the Lahu and Wa nationalities. The D-loop sequence variation and phylogenetic analysis suggested that the 99 mtDNA lineages were classified into eight clusters in the phylogenetic tree. All lineages that had a 9-bp deletion in the COII/tRNALYs intergenic region appeared in one cluster in the D-loop tree, suggesting a single event of the deletion in the Yunnan nationalities studied. Genetic distances, based on net nucleotide diversities between populations including Han Chinese and mainland Japanese, revealed that the Dai, Wa, Lahu, and Han Chinese are closely related to each other, while Tibetan and mainland Japanese formed a single cluster. The bootstrap probability of separation between the Dai-Wa-Lahu-Chinese clade and the Tibetan-Japanese clade was 99%, indicating that there are at least two different origins among minority groups in Yunnan province. Although the genetic distance between Tibetan and Japanese within the clade is rather long, the results may shed light on the origins of mainland Japanese.

Association of Fcgamma receptor IIb polymorphism with susceptibility to systemic lupus erythematosus in Chinese: a common susceptibility gene in the Asian populations.

The association of Fcgamma receptor (FcgammaR) polymorphisms with systemic lupus erythematosus (SLE) has been demonstrated in various populations; however, the results have been inconsistent. We recently identified a single-nucleotide polymorphism encoding a non-synonymous substitution, Ile232Thr (I232T), of FCGR2B and its association with SLE in Japanese and in Thais. Multiple functional FcgammaR genes with polymorphisms (FCGR2A, FCGR2B, FCGR3A, and FCGR3B) cluster in 1q23, and some of them are in linkage disequilibrium (LD). To differentiate contributions from multiple-linked loci, comparison of different populations may provide useful information. In this study, we analyzed the above four FCGR polymorphisms of the Chinese patients and controls for the association with SLE. FCGR2A-H131R, FCGR2B-I232T, FCGR3A-F176V, and FCGR3B genotypes were determined in 167 Chinese patients with SLE and 129 healthy controls. Association was examined using case-control analysis. Allele frequencies of FCGR2B-232T and FCGR3A-176F were significantly increased in SLE [odds ratio (OR) = 1.67 and OR = 1.41, respectively]. Interestingly, while these alleles had a tendency of positive LD in the controls, FCGR2B-232T was in positive association with FCGR3A-176V in SLE, suggesting that these two alleles were associated with SLE in an independent manner. Comparison between SLE with and without nephritis indicated significant association of FCGR2B-232T with nephritis (OR = 2.65). When the present results were combined with our previous data on the Japanese and the Thais using meta-analytic methods, highly significant and independent association was observed for FCGR2B and FCGR3A genotypes. These results strongly suggested that FCGR2B is a common susceptibility factor to SLE in the Asians.

Genetic relationship of populations in China.

Despite the fact that the continuity of morphology of fossil specimens of modern humans found in China has repeatedly challenged the Out-of-Africa hypothesis, Chinese populations are underrepresented in genetic studies. Genetic profiles of 28 populations sampled in China supported the distinction between southern and northern populations, while the latter are biphyletic. Linguistic boundaries are often transgressed across language families studied, reflecting substantial gene flow between populations. Nevertheless, genetic evidence does not support an independent origin of Homo sapiens in China. The phylogeny also suggested that it is more likely that ancestors of the populations currently residing in East Asia entered from Southeast Asia.