PubMed Health⌕ Search

Biomedical subjects

Wenjun Zheng

Publications and source records attributed to Wenjun Zheng.

15 recordsLinked to original sources

Experimental and theoretical study of the coordination of 1,2,4-triazolato, tetrazolato, and pentazolato ligands to the [K(18-crown-6)]+ fragment.

Treatment of 3,5-diisopropyltriazole, 3,5-diphenyltriazole, 3,5-di-3-pyridyltriazole, phenyltetrazole, pyrrolidinyltetrazole, or tert-butyltetrazole with equimolar quantities of potassium hydride and 18-crown-6 in tetrahydrofuran at ambient temperature led to slow hydrogen evolution and formation of (3,5-diisopropyl-1,2,4-triazolato)(18-crown-6)potassium (88%), (3,5-diphenyl-1,2,4-triazolato)(tetrahydrofuran)(18-crown-6)potassium (87%), (3,5-di-3-pyridyl-1,2,4-triazolato)(18-crown-6)potassium (81%), (phenyltetrazolato)(18-crown-6)potassium (94%), (pyrrolidinyltetrazolato)(18-crown-6)potassium (90%), and (tert-butyltetrazolato)(18-crown-6)potassium (94%) as colorless crystalline solids. (1,2,4-Triazolato)(18-crown-6)potassium was isolated as a hemi-hydrate in 81% yield upon treatment of 1,2,4-triazole with potassium metal in tetrahydrofuran. The X-ray crystal structures of these new complexes were determined, and the solid-state structures consist of the nitrogen heterocycles bonded to the (18-crown-6)potassium cationic fragments with eta2-bonding interactions. In addition, (3,5-diphenyl-1,2,4-triazolato)(tetrahydrofuran)(18-crown-6)potassium has one coordinated tetrahydrofuran ligand on the same face as the 3,5-diphenyl-1,2,4-triazolato ligand, while (3,5-di-3-pyridyl-1,2,4-triazolato)(18-crown-6)potassium forms a polymeric solid through coordination of the distal 3-pyridyl nitrogen atoms to the potassium ion on the face opposite the 1,2,4-triazolato ligand. The solid-state structures of the new complexes show variable asymmetry in the potassium-nitrogen distances within the eta2-interactions and also show variable bending of the heterocyclic C2N3 and CN4 cores toward the best plane of the 18-crown-6 ligand oxygen atoms. Molecular orbital and natural bond order calculations were carried out at the B3LYP/6-311G(d,p) level of theory on the model complex, (phenyltetrazolato)(18-crown-6)potassium, and demonstrate that the asymmetric potassium-nitrogen distances and bending of the CN4 core toward the 18-crown-6 ligand are due to hydrogen bond-like interactions between filled nitrogen-based orbitals and carbon-hydrogen sigma orbitals on the 18-crown-6 ligands. Calculations carried out on the model pentazolato complex (pentazolato)(18-crown-6)potassium predict a structure in which the pentazolato ligand N5 core is bent by 45 degrees toward the best plane of the 18-crown-6 oxygen atoms. Such bending is induced by the formation of intramolecular nitrogen-hydrogen-carbon hydrogen bonds. Examination of the solid-state structures of the new complexes reveals many intramolecular and intermolecular nitrogen-hydrogen distances of < or =3.0 A which support the presence of nitrogen-hydrogen-carbon hydrogen bonds.

Journal Article↗

Low-frequency normal modes that describe allosteric transitions in biological nanomachines are robust to sequence variations.

By representing the high-resolution crystal structures of a number of enzymes using the elastic network model, it has been shown that only a few low-frequency normal modes are needed to describe the large-scale domain movements that are triggered by ligand binding. Here we explore a link between the nearly invariant nature of the modes that describe functional dynamics at the mesoscopic level and the large evolutionary sequence variations at the residue level. By using a structural perturbation method (SPM), which probes the residue-specific response to perturbations (or mutations), we identify a sparse network of strongly conserved residues that transmit allosteric signals in three structurally unrelated biological nanomachines, namely, DNA polymerase, myosin motor, and the Escherichia coli chaperonin. Based on the response of every mode to perturbations, which are generated by interchanging specific sequence pairs in a multiple sequence alignment, we show that the functionally relevant low-frequency modes are most robust to sequence variations. Our work shows that robustness of dynamical modes at the mesoscopic level is encoded in the structure through a sparse network of residues that transmit allosteric signals.

Allosteric Regulation↗

Modeling protein conformational changes by iterative fitting of distance constraints using reoriented normal modes.

Recently we have developed a normal-modes-based algorithm that predicts the direction of protein conformational changes given the initial state crystal structure together with a small number of pairwise distance constraints for the end state. Here we significantly extend this method to accurately model both the direction and amplitude of protein conformational changes. The new protocol implements a multisteps search in the conformational space that is driven by iteratively minimizing the error of fitting the given distance constraints and simultaneously enforcing the restraint of low elastic energy. At each step, an incremental structural displacement is computed as a linear combination of the lowest 10 normal modes derived from an elastic network model, whose eigenvectors are reorientated to correct for the distortions caused by the structural displacements in the previous steps. We test this method on a list of 16 pairs of protein structures for which relatively large conformational changes are observed (root mean square deviation >3 angstroms), using up to 10 pairwise distance constraints selected by a fluctuation analysis of the initial state structures. This method has achieved a near-optimal performance in almost all cases, and in many cases the final structural models lie within root mean square deviation of 1 approximately 2 angstroms from the native end state structures.

Computer Simulation↗

Probing the local dynamics of nucleotide-binding pocket coupled to the global dynamics: myosin versus kinesin.

Based on the elastic network model, we develop a new analysis for protein complexes, which probes the local dynamics of a subsystem that is elastically coupled to a fluctuating environment. This method is applied to a comparative dynamical analysis of the nucleotide-binding pocket of two motor proteins-myosins and kinesins. In myosins, the observed structural changes in the nucleotide-pocket from the transition state to the rigorlike state are dominated by the lowest normal mode that involves significant movements in both switch I and switch II; in kinesins, the measured conformational changes in the nucleotide-pocket are also dominated by the lowest mode, which, however, only involves large movement in switch I. We then compute the global structural changes induced by the nucleotide-pocket deformations as described by the dominant pocket-mode, which yield encouraging results: in myosins, multiple hinge motions involving the opening/closing of the cleft between the upper and lower 50 -kDa subdomains and the swinging movement of the converter are induced, which are dominated by precisely the same global mode that has been recently identified by us as important to the dynamical correlations among the nucleotide-pocket, the actin-binding site, and the converter; in kinesins, the induced global conformational changes are well described by a highly collective global mode which hints for a dynamical pathway spanning from the nucleotide-pocket to the neck-linker via the H6 helix.

Actins↗

Fold recognition aided by constraints from small angle X-ray scattering data.

We performed a systematic exploration of the use of structural information derived from small angle X-ray scattering (SAXS) measurements to improve fold recognition. SAXS data provide the Fourier transform of the histogram of atomic pair distances (pair distribution function) for a given protein and hence can serve as a structural constraint on methods used to determine the native conformational fold of the protein. Here we used it to construct a similarity-based fitness score with which to evaluate candidate structures generated by a threading procedure. In order to combine the SAXS scores with the standard energy scores and other 1D profile-based scores used in threading, we made use both of a linear regression method and of a neural network-based technique to obtain optimal combined fitness scores and applied them to the ranking of candidate structures. Our results show that the use of SAXS data with gapless threading significantly improves the performance of fold recognition.

Computational Biology↗

Normal-modes-based prediction of protein conformational changes guided by distance constraints.

Based on the elastic network model, we develop a novel method that predicts the conformational change of a protein complex given its initial-state crystal structure together with a small set of pairwise distance constraints for the end state. The predicted conformational change, which is a linear combination of multiple low-frequency normal modes that are solved from the elastic network model, is computed as a response displacement induced by a perturbation to the system Hamiltonian that incorporates the given distance constraints. For a list of test cases, we find that the computed response displacement overlaps significantly with the measured conformational changes, when only a handful of pairwise constraints are used (</=10). The performance of this method is also shown to be robust against different choices of pairwise distance constraints and errors in their values. This method, if supplied with the experimentally derived distance constraints (for example, from NMR or other spectroscopic measurements), can be applied to the analysis of protein conformational changes toward transient states.

Algorithms↗

Identification of dynamical correlations within the myosin motor domain by the normal mode analysis of an elastic network model.

In order to systematically analyze functionally relevant dynamical correlations within macromolecular complexes, we have developed computational methods based on the normal mode analysis of an elastic network model. First, we define two types of dynamical correlations (fluctuation-based and density-based), which are computed by summing up contributions from all low-frequency normal modes up to a given cutoff. Then we use them to select dynamically important "hinge residues" whose elastic distortion affects the fluctuations of a large number of residues. Second, in order to clarify long-range dynamical correlations, we decompose the dynamical correlations to individual normal modes to identify the most relevant modes. We have applied these methods to the analysis of the motor domain of Dictyostelium myosin and have obtained the following three interesting results that shed light on its mechanism of force generation: first, we find the hinge residues are distributed over several key inter-subdomain joints (including the nucleotide-binding pocket, the relay helix, the SH1 helix, the strut between the upper 50 kDa and the lower 50 kDa subdomains), which is consistent with their hypothesized roles in modulating functionally relevant inter-subdomain conformational changes; second, a single mode 7 (for structure 1VOM) is found to dominate the fluctuation-based correlations between the converter/strut and the nucleotide-binding pocket, revealing a surprising simplicity for their intriguing roles in the force generation mechanism; finally, we find a negative density-based correlation between the strut and the nucleotide-binding pocket, which is consistent with the hypothesized signaling pathway that links the actin-binding site's opening/closing with the nucleotide-binding pocket's closing/opening.

Animals↗

Network of dynamically important residues in the open/closed transition in polymerases is strongly conserved.

The open/closed transition in polymerases is a crucial event in DNA replication and transcription. We hypothesize that the residues that transmit the signal for the open/closed transition are also strongly conserved. To identify the dynamically relevant residues, we use an elastic network model of polymerases and probe the residue-specific response to a local perturbation. In a variety of DNA/RNA polymerases, a network of residues spanning the fingers and palm domains is involved in the open/closed transition. The similarity in the network of residues responsible for large-scale domain movements supports the notion of a common induced-fit mechanism in the polymerase families for the formation of a closed ternary complex. Multiple sequence alignment shows that many of these residues are also strongly conserved. Residues with the largest sensitivity to local perturbations include those that are not so obviously involved in the polymerase catalysis. Our results suggest that mutations of the mechanical "hot spots" can compromise the efficiency of the enzyme.

Animals↗

A comparative study of motor-protein motions by using a simple elastic-network model.

In this work, we report on a study of the structure-function relationships for three families of motor proteins, including kinesins, myosins, and F1-ATPases, by using a version of the simple elastic-network model of large-scale protein motions originally proposed by Tirion [Tirion, M. (1996) Phys. Rev. Lett. 77, 1905-1908]. We find a surprising dichotomy between kinesins and the other motor proteins (myosins and F1-ATPase). For the latter, there exist one or two dominant lowest-frequency modes (one for myosin, two for F1-ATPase) obtained from normal-mode analysis of the elastic-network model, which overlap remarkably well with the measured conformational changes derived from pairs of solved crystal structures in different states. Furthermore, we find that the computed global conformational changes induced by the measured deformation of the nucleotide-binding pocket also overlap well with the measured conformational changes, which is consistent with the "nucleotide-binding-induced power-stroke" scenario. In contrast, for kinesins, this simplicity breaks down. Multiple modes are needed to generate the measured conformational changes, and the computed displacements induced by deforming the nucleotide-binding pocket also overlap poorly with the measured conformational changes, and are insufficient to explain the large-scale motion of the relay helix and the linker region. This finding may suggest the presence of two different mechanisms for myosins and kinesins, despite their strong evolutionary ties and structural similarities.

Adenosine Triphosphate↗

Syntheses, structures, and surface aromaticity of the new carbaalane [(AlH)(6)(AlNMe(3))(2)(CCH(2)R)(6)] (R = Ph, CH(2)SiMe(3)) and a stepwise functionalization of the inner and outer sphere of the cluster.

The reaction of the acetylene RC triple bond CH (R = Ph, CH(2)SiMe(3)) with an excess of AlH(3).NMe(3) in boiling toluene leads to the carbaalane [(AlH)(6)(AlNMe(3))(2)(CCH(2)R)(6)] (R = Ph 1, CH(2)SiMe(3) 2) in good yield. Treatment of 2 with BCl(3) under varying conditions gives the chlorinated products [(AlCl)(6)(AlNMe(3))(2)(CCH(2)CH(2)SiMe(3))(6)] 3 and [(AlCl)(6)(AlNMe(3))(2)(CCH(2)CH(2)SiMe(2)Cl)(6)] 4, respectively. The latter clearly demonstrates that the cluster can be stepwise functionalized within the inner and outer sphere. The X-ray single-crystal structures of 1, 2, and 4 have been determined. All compounds have in common that the central core consists of a cluster having eight aluminum and six carbon atoms. The bonding properties in this cluster are described as a new manifestation of three-dimensional surface aromaticity. Each Al(4)C fragment of the cube is formed by four bonds with three electron pairs, thus leading to a strong delocalization of the electrons. A phenomenological modeling using a three-dimensional Hückel scheme with fitted parameters to reproduce the energies from ab initio calculations revealed that the orbital scheme localized at one Al(4)C fragment possesses an orbital sextet with a large HOMO-LUMO gap. This is in line with the criteria of aromaticity. The idea of aromaticity was sustained also by qualitative valence bond reasons enumerating the different resonance structures by means of graph theoretical methods.

Journal Article↗

Protein structure prediction constrained by solution X-ray scattering data and structural homology identification.

Here we perform a systematic exploration of the use of distance constraints derived from small angle X-ray scattering (SAXS) measurements to filter candidate protein structures for the purpose of protein structure prediction. This is an intrinsically more complex task than that of applying distance constraints derived from NMR data where the identity of the pair of amino acid residues subject to a given distance constraint is known. SAXS, on the other hand, yields a histogram of pair distances (pair distribution function), but the identities of the pairs contributing to a given bin of the histogram are not known. Our study is based on an extension of the Levitt-Hinds coarse grained approach to ab initio protein structure prediction to generate a candidate set of C(alpha) backbones. In spite of the lack of specific residue information inherent in the SAXS data, our study shows that the implementation of a SAXS filter is capable of effectively purifying the set of native structure candidates and thus provides a substantial improvement in the reliability of protein structure prediction. We test the quality of our predicted C(alpha) backbones by doing structural homology searches against the Dali domain library, and find that the results are very encouraging. In spite of the lack of local structural details and limited modeling accuracy at the C(alpha) backbone level, we find that useful information about fold classification can be extracted from this procedure. This approach thus provides a way to use a SAXS data based structure prediction algorithm to generate potential structural homologies in cases where lack of sequence homology prevents identification of candidate folds for a given protein. Thus our approach has the potential to help in determination of the biological function of a protein based on structural homology instead of sequence homology.

Computational Biology↗