Biomedical subjects

Buyong Ma

Publications and source records attributed to Buyong Ma.

At least 19 recordsLinked to original sources

PubMed2006

Consensus features in amyloid fibrils: sheet-sheet recognition via a (polar or nonpolar) zipper structure.

Amyloid fibrils characterized as highly intractable thread-like species are associated with many neurodegenerative diseases. Although neither the mechanism of amyloid formation nor the origin of amyloid toxicity is currently completely understood, the detailed three-dimensional atomic structures of the yeast protein Sup35 and Abeta amyloid protein determined by recent experiments provide the first and important step towards the comprehension of the pathogenesis and aggregation mechanisms of amyloid diseases. By analyzing these two amyloid peptides which have available crystal structures and other amyloid sequences with proposed structures using computational simulations, we delineate three common features in amyloid organizations and amyloid structures. These could contribute to an improved understanding of the molecular mechanism of amyloid formation, the nature of the aggregation driving forces that stabilize these structures and the development of potential therapeutic agents against amyloid diseases.

Amino Acid Motifs↗

PubMed2006

Simulations as analytical tools to understand protein aggregation and predict amyloid conformation.

Computational tools are increasingly being applied to solve the protein aggregation problem, providing insight into amyloid structures and aggregation mechanisms. The paradigm of Abeta amyloid structure elucidation provides an example of an innovative experimental design and endeavor, echoing the computational testing of possible molecular associations, all reflected in the current Ma-Nussinov-Tycko model of the Abeta amyloid. Simulations have shown that dimer formation can lock some misfolded conformations, and catalyze the shift of the equilibrium away from the native state. In most cases, a stable amyloid seed requires at least two-layered beta-sheets with properly registered side-chains. Under kinetic control, the final protein aggregations are the outcome of maximizing the van der Waals interactions between side chains and backbone hydrogen bonds.

Amyloid↗

PubMed2006

Interaction of protegrin-1 with lipid bilayers: membrane thinning effect.

Protegrins (PG) are important in defending host tissues, preventing infection via an attack on the membrane surface of invading microorganisms. Protegrins have powerful antibiotic abilities, but the molecular-level mechanisms underlying the interactions of their beta-sheet motifs with the membrane are not known. Protegrin-1 (PG-1) is composed of 18 amino acids with a high content of basic residues and two disulfide bonds. Here we focused on the stability of PG-1 at the amphipathic interface in lipid bilayers and on the details of the peptide-membrane interactions. We simulated all-atom models of the PG-1 monomer with explicit water and lipid bilayers composed of both homogeneous POPC (palmitoyl-oleyl-phosphatidylcholine) lipids and a mixture of POPC/POPG (palmitoyl-oleyl-phosphatidylglycerol) (4:1) lipids. We observed that local thinning of the lipid bilayers mediated by the peptide is enhanced in the lipid bilayer containing POPG, consistent with experimental results of selective membrane targeting. The beta-hairpin motif of PG-1 is conserved in both lipid settings, whereas it is highly bent in aqueous solution. The conformational dynamics of PG-1, especially the highly charged beta-hairpin turn region, are found to be mostly responsible for disturbing the membrane. Even though the eventual membrane disruption requires PG-1 oligomers, our simulations clearly show the first step of the monomeric effects. The thinning effects in the bilayer should relate to pore/channel formation in the lipid bilayer and thus be responsible for further defects in the membrane caused by oligomer.

Antimicrobial Cationic Peptides↗

PubMed2006

Structural stability and dynamics of an amyloid-forming peptide GNNQQNY from the yeast prion sup-35.

A seven amino acid yeast prion sup-35 fragment (GNNQQNY) forms amyloid fibrils. The availability of its detailed atomic oligomeric structure makes it a good model for studying the early stage of aggregation. Here we perform long all-atom explicit solvent molecular simulations of various sizes and arrangements of oligomer seeds of the wild-type and its mutants to study its stability and dynamics. Previous studies have suggested that the early stage rate-limiting step of oligomer formation occurs in high-order oligomers. Our simulations show that with the increase in the number of strands even from a dimer to a trimer, oligomer stability increases dramatically. This suggests that the minimal nucleus seed for GNNQQNY fibril formation could be small and is likely three or four peptides, in agreement with experiment, and that higher-order oligomers do not dissociate quickly since they have small diffusion coefficients and thus slow kinetics. Further, for the hydrophilic polar GNNQQNY, there are no hydrogen bonds and no hydrophobic interactions between adjacent beta-sheets. Simulations suggest that within the sheet, the driving forces to associate and stabilize are interstrand backbone-backbone and side chain-side chain hydrogen bonds, whereas between the sheets, shape-complementary by the dry polar steric zipper via the side chains of Asn-2, Gln-4, and Asn-6 holds the sheets together, as proposed in an earlier study. Since the polar side chains of Asn-2, Gln-4, and Asn-6 act as a hook to bind two neighboring sheets together, these geometric restraints reduce the conformational search for the correct side chain packing to a two-dimensional problem of intersheet side chain interactions. Mutant simulations show that substitution of Asn-2, Gln-4, or Asn-6 by Ala would disrupt this steric zipper, leading to unstable oligomers.

Amino Acid Sequence↗

PubMed2006

Comparison of the human and worm p53 structures suggests a way for enhancing stability.

Maintaining the native conformation is essential for the proper function of tumor suppressor protein p53. However, p53 is a low-stability protein that can easily lose its function upon structural perturbations such as those resulting from missense mutations, leading to the development of cancer. Therefore, it is important to develop strategies to design stable p53 which still maintains its normal function. Here, we compare the stabilities of the human and worm p53 core domains using molecular dynamics simulations. We find that the worm p53 is significantly more stable than the human form. Detailed analysis of the structural fluctuations shows that the stability difference lies in the peripheral structural motifs that contrast in their structural features and flexibility. The most dramatic difference in stability originates from loop L1, from the turn between helix H1 and beta-strand S5, and from the turn that connects beta-strands S7 and S8. Structural analysis shows significant differences for these motifs between the two proteins. Loop L1 lacks secondary structure, and the turns between helix H1 and strand S5 and between strands S7 and S8 are much longer in the human form p53. On the basis of these differences, we designed a mutant by shortening the turn between strands S7 and S8 to enhance the stability. Surprisingly, this mutant was very stable when probed by molecular dynamics simulations. In addition, the stabilization was not localized in the turn region. Loop L1 was also significantly stabilized. Our results show that stabilizing peripheral structural motifs can greatly enhance the stability of the p53 core domain and therefore is likely to be a viable alternative in the development of stable p53. In addition, loop- or turn-related mutants with different stabilities may also be used to probe the relationship between function, a particular structural motif, and its flexibility.

Animals↗

PubMed2006

The stability of monomeric intermediates controls amyloid formation: Abeta25-35 and its N27Q mutant.

The structure and stabilities of the intermediates affect protein folding as well as misfolding and amyloid formation. By applying Kramer's theory of barrier crossing and a Morse-function-like energy landscape, we show that intermediates with medium stability dramatically increase the rate of amyloid formation; on the other hand, very stable and very unstable intermediates sharply decrease amyloid formation. Remarkably, extensive molecular dynamics simulations and conformational energy landscape analysis of Abeta25-35 and its N27Q mutant corroborate the mathematical description. Both experimental and current simulation results indicate that the core of the amyloid structure of Abeta25-35 formed from residues 28-35. A single mutation of N27Q of Abeta25-35 makes the Abeta25-35 N27Q amyloid-free. Energy landscape calculations show that Abeta25-35 has extended intermediates with medium stability that are prone to form amyloids, whereas the extended intermediates for Abeta25-35 N27Q split into stable and very unstable species that are not disposed to form amyloids. The results explain the contribution of both alpha-helical and beta-strand intermediates to amyloid formation. The results also indicate that the structure and stability of the intermediates, as well as of the native folded and the amyloid states can be targeted in drug design. One conceivable approach is to stabilize the intermediates to deter amyloid formation.

Amino Acid Substitution↗

PubMed2006

"Similarity trap" in protein-protein interactions could be carcinogenic: simulations of p53 core domain complexed with 53BP1 and BRCA1 BRCT domains.

Similar binding sites often imply similar protein-protein interactions and similar functions; however, similar binding sites may also constitute traps for nonfunctional associations. How are similar sites distinguished to prevent misassociations? BRCT domains from breast cancer-susceptibility gene product BRCA1 and protein 53BP1 have similar structures yet different binding behaviors with p53 core domain. 53BP1-BRCT domain forms a stable complex with p53. In contrast, BRCA1-p53 interaction is weak or other mechanisms operate. To delineate the difference, we designed 13 BRCA1-BRCT mutants and computationally investigated the structural and stability changes compared to the experimental p53-53BP1 structure. Interestingly, of the 13, the 2 mutations that are cancerous and involve nonconserved residues are those that enforced p53 core domain binding with BRCA1-BRCT in a way similar to p53-53BP1 binding. Hence, falling into the "similarity trap" may disrupt normal BRCA1 and p53 functions. Our results illustrate how this trap is avoided in the native state.

Amino Acid Sequence↗

PubMed2005

CD4 binding partially locks the bridging sheet in gp120 but leaves the beta2/3 strands flexible.

The structure of the free form HIV gp120, critical for therapeutic agent development, is unavailable due to its high flexibility. Previous thermodynamic data, structural analysis and simulation results have suggested a large conformational change in the core domain upon CD4 binding. The bridging sheet, which consists of four beta-strands with beta20/21 nestling against the inner/outer domains and beta2/3 facing outward, more exposed to the solvent, was proposed to be unfolded in the native state. In order to test this proposition and to characterize the native conformations, we performed potential mean force (PMF) molecular dynamics (MD) simulations on the CD4-bound crystal structure. We pushed the bridging sheet away from the inner and outer domain to explore the accessible conformational space for the bridging sheet. In addition, we performed conventional MD simulations on structures with the bridging sheet partially unfolded to investigate the stability of the association between the inner and outer domains. Based on the free energy profiles, we find that the whole bridging sheet is unlikely to unfold without other concurrent conformational changes. On the other hand, the partial bridging sheet, beta strands 2/3, can switch its conformation from the folded to the unfolded state. Furthermore, relaxation of conformation with partially unfolded bridging sheet through MD simulations leads to a conformation with beta strands 20/21 quickly re-anchoring against the inner and outer domains. Such a conformation, although lacking some of the hydrophobic interactions present in the CD4-bound structure, displayed high stability as further indicated by other restrained MD simulations. The relevance of this conformation to the free form structure and the pathway for conformational change from the free form to the CD4-bound structure is discussed in detail in light of the available unliganded SIV gp120 crystal structure.

Amino Acid Motifs↗

PubMed2005

Comparison of the protein-protein interfaces in the p53-DNA crystal structures: towards elucidation of the biological interface.

p53, the tumor suppressor protein, functions as a dimer of dimers. However, how the tetramer binds to the DNA is still an open question. In the crystal structure, three copies of the p53 monomers (containing chains A, B, and C) were crystallized with the DNA-consensus element. Although the structure provides crucial data on the p53-DNA contacts, the active oligomeric state is unclear because the two dimeric (A-B and B-C) interfaces present in the crystal cannot both exist in the tetramer. Here, we address the question of which of these two dimeric interfaces may be more biologically relevant. We analyze the sequence and structural properties of the p53-p53 dimeric interfaces and carry out extensive molecular dynamics simulations of the crystal structures of the human and mouse p53 dimers. We find that the A-B interface residues are more conserved than those of the B-C. Molecular dynamics simulations show that the A-B interface can provide a stable DNA-binding motif in the dimeric state, unlike B-C. Our results indicate that the interface between chains A-B in the p53-DNA complex constitutes a better candidate for a stable biological interface, whereas the B-C interface is more likely to be due to crystal packing. Thus, they have significant implications toward our understanding of DNA binding by p53 as well as p53-mediated interactions with other proteins.

Amino Acid Motifs↗

PubMed2005

In the quest for stable rescuing mutants of p53: computational mutagenesis of flexible loop L1.

p53 is a protein with marginal stability. Its transcriptional functions are often inactivated by single missense mutations, shown to be associated with half of all human cancers. Here, we aim to design stable functional p53 mutants. We target loop L1, one of the most mobile structural motifs in the p53 core domain (p53C). Specifically, we selected Ser116 in the middle of loop L1 and mutated it to 14 other amino acids. All resulting mutants were subjected to molecular dynamics simulations, revealing a wide spectrum of stabilities. Among these, mutant S116M displayed a remarkable stability, with a structural deviation comparable to that of the experimental quadruple mutant M133L/V203A/N239Y/N268D that is thermodynamically more stable than that of the wild type by 2.6 kcal/mol. Structural analysis showed that the high stability of the S116M mutant was indeed due to the preservation of the p53C loop L1 conformation and the reduction of mobility in that region. The differential stabilities conferred by the single mutations are rationalized based on the geometries and chemical properties of the side chains introduced into this site. Linearity (i.e., nonbranched), moderate size, and balanced hydrophobic and hydrophilic properties of the side chain are crucial to the stabilizing effect of the residue substitutions.

Amino Acid Substitution↗

PubMed2005

Protein-protein interactions: organization, cooperativity and mapping in a bottom-up Systems Biology approach.

Understanding and ultimately predicting protein associations is immensely important for functional genomics and drug design. Here, we propose that binding sites have preferred organizations. First, the hot spots cluster within densely packed 'hot regions'. Within these regions, they form networks of interactions. Thus, hot spots located within a hot region contribute cooperatively to the stability of the complex. However, the contributions of separate, independent hot regions are additive. Moreover, hot spots are often already pre-organized in the unbound (free) protein states. Describing a binding site through independent local hot regions has implications for binding site definition, design and parametrization for prediction. The compactness and cooperativity emphasize the similarity between binding and folding. This proposition is grounded in computation and experiment. It explains why summation of the interactions may over-estimate the stability of the complex. Furthermore, statistically, charge-charge coupling of the hot spots is disfavored. However, since within the highly packed regions the solvent is screened, the electrostatic contributions are strengthened. Thus, we propose a new description of protein binding sites: a site consists of (one or a few) self-contained cooperative regions. Since the residue hot spots are those conserved by evolution, proteins binding multiple partners at the same sites are expected to use all or some combination of these regions.

Binding Sites↗

PubMed2005

The contribution of the Trp/Met/Phe residues to physical interactions of p53 with cellular proteins.

Dynamic molecular interaction networks underlie biological phenomena. Among the many genes which are involved, p53 plays a central role in networks controlling cellular life and death. It not only operates as a tumor suppressor, but also helps regulate hundreds of genes in response to various types of stress. To accomplish these functions as a guardian of the genome, p53 interacts extensively with both nucleic acids and proteins. This paper examines the physical interfaces of the p53 protein with cellular proteins. Previously, in the analysis of the structures of protein-protein complexes, we have observed that amino acids Trp, Met and Phe are important for protein-protein interactions in general. Here we show that these residues are critical for the many functions of p53. Several clusters of the Trp/Met/Phe residues are involved in the p53 protein-protein interactions. Phe19/Trp23 in the TA1 region extensively binds to the transcriptional factors and the MDM2 protein. Trp53/Phe54 in the TA2 region is crucial for transactivation and DNA replication. Met243 in the core domain interacts with 53BP1, 53BP2 and Rad 51 proteins. Met384/Phe385 in the C-terminal region interacts with the S100B protein and the Bromodomain of the CBP protein. Thus, these residues may assist in elucidating the p53 interactions when structural data are not available.

Biophysics↗

PubMed2005

From structure to function: methods and applications.

The rapid increase in experimental data along with recent progress in computational methods has brought modern biology a step closer toward solving one of the most challenging problems: prediction of protein function. Comprehension of protein function at its most basic level requires understanding of molecular interactions. Currently, it is becoming universally accepted that the scale of the accumulated data for analysis and for prediction necessitate highly efficient computational tools with appropriate application capabilities. The review presents the up-to-date advances in computational methods for structural pattern discovery and for prediction of molecular associations. We focus on their applications toward a range of biological problems and highlight the advantages of the combination of these methods and their integration with biological experiments. We provide examples, synergistically merging structural modeling, rigid and flexible structural alignment and detection of conserved structural patterns and docking (rigid and flexible with hinge-bending movements). We hope the review will lead to a broader utilization of computational methods, and their cross-fertilization with experiment.

Algorithms↗

PubMed2004

A comparative study of amyloid fibril formation by residues 15-19 of the human calcitonin hormone: a single beta-sheet model with a small hydrophobic core.

Experimentally, the human calcitonin hormone (hCT) can form highly stable amyloid protofibrils. Further, a peptide consisting of hCT residues 15-19, DFNKF, was shown to create highly ordered fibrils, similar to those formed by the entire hormone sequence. However, there are limited experimental data regarding the detailed 3D arrangement of either of these fibrils. We have modeled the DFNKF protofibril, using molecular dynamics simulations. We tested the stabilities of single sheet and of various multi sheet models. Remarkably, our most ordered and stable model consists of a parallel-stranded, single beta-sheet with a relatively insignificant hydrophobic core. We investigate the chemical and physical interactions responsible for the high structural organization of this single beta-sheet amyloid fibril. We observe that the most important chemical interactions contributing to the stability of the DFNKF organization are electrostatic, specifically between the Lys and the C terminus, between the Asp and N terminus, and a hydrogen bond network between the Asn side-chains of adjacent strands. Additionally, we observe hydrophobic and aromatic pi stacking interactions. We further simulated truncated filaments, FNKF and DFNK. Our tetra-peptide mutant simulations assume models similar to the penta-peptide. Experimentally, the FNKF does not create fibrils while DFNK does, albeit short and less ordered than DFNKF. In the simulations, the FNKF system was less stable than the DFNK and DFNKF. DFNK also lost many of its original interactions becoming less organized, however, many contacts were maintained. Thus, our results emphasize the role played by specific amino acid interactions. To further study specific interactions, we have mutated the penta-peptide, simulating DANKF, DFNKA and EFNKF. Here we describe the model, its relationship to experiment and its implications to amyloid organization.

Amino Acid Sequence↗

PubMed2004

How similar are protein folding and protein binding nuclei? Examination of vibrational motions of energy hot spots and conserved residues.

The underlying physico-chemical principles of the interactions between domains in protein folding are similar to those between protein molecules in binding. Here we show that conserved residues and experimental hot spots at intermolecular binding interfaces overlap residues that vibrate with high frequencies. Similarly, conserved residues and hot spots are found in protein cores and are also observed to vibrate with high frequencies. In both cases, these residues contribute significantly to the stability. Hence, these observations validate the proposition that binding and folding are similar processes. In both packing plays a critical role, rationalizing the residue conservation and the experimental alanine scanning hot spots. We further show that high-frequency vibrating residues distinguish between protein binding sites and the remainder of the protein surface.

Amino Acids↗

PubMed2004

Hot regions in protein--protein interactions: the organization and contribution of structurally conserved hot spot residues.

Structurally conserved residues at protein-protein interfaces correlate with the experimental alanine-scanning hot spots. Here, we investigate the organization of these conserved, computational hot spots and their contribution to the stability of protein associations. We find that computational hot spots are not homogeneously distributed along the protein interfaces; rather they are clustered within locally tightly packed regions. Within the dense clusters, they form a network of interactions and consequently their contributions to the stability of the complex are cooperative; however the contributions of independent clusters are additive. This suggests that the binding free energy is not a simple summation of the single hot spot residue contributions. As expected, around the hot spot residues we observe moderately conserved residues, further highlighting the crucial role of the conserved interactions in the local densely packed environment. The conserved occurrence of these organizations suggests that they are advantageous for protein-protein associations. Interestingly, the total number of hydrogen bonds and salt bridges contributed by hot spots is as expected. Thus, H-bond forming residues may use a "hot spot for water exclusion" mechanism. Since conserved residues are located within highly packed regions, water molecules are easily removed upon binding, strengthening electrostatic contributions of charge-charge interactions. Hence, the picture that emerges is that protein-protein associations are optimized locally, with the clustered, networked, highly packed structurally conserved residues contributing dominantly and cooperatively to the stability of the complex. When addressing the crucial question of "what are the preferred ways of proteins to associate", these findings point toward a critical involvement of hot regions in protein-protein interactions.

Binding Sites↗

PubMed2004

Protein-protein interactions: hot spots and structurally conserved residues often locate in complemented pockets that pre-organized in the unbound states: implications for docking.

Energetic hot spots account for a significant portion of the total binding free energy and correlate with structurally conserved interface residues. Here, we map experimentally determined hot spots and structurally conserved residues to investigate their geometrical organization. Unfilled pockets are pockets that remain unfilled after protein-protein complexation, while complemented pockets are pockets that disappear upon binding, representing tightly fit regions. We find that structurally conserved residues and energetic hot spots are strongly favored to be located in complemented pockets, and are disfavored in unfilled pockets. For the three available protein-protein complexes with complemented pockets where both members of the complex were alanine-scanned, 62% of all hot spots (DeltaDeltaG>2kcal/mol) are within these pockets, and 60% of the residues in the complemented pockets are hot spots. 93% of all red-hot residues (DeltaDeltaG>/=4kcal/mol) either protrude into or are located in complemented pockets. The occurrence of hot spots and conserved residues in complemented pockets highlights the role of local tight packing in protein associations, and rationalizes their energetic contribution and conservation. Complemented pockets and their corresponding protruding residues emerge among the most important geometric features in protein-protein interactions. By screening the solvent, this organization shields backbone hydrogen bonds and charge-charge interactions. Complemented pockets often pre-exist binding. For 18 protein-protein complexes with complemented pockets whose unbound structures are available, in 16 the pockets are identified to pre-exist in the unbound structures. The root-mean-squared deviations of the atoms lining the pockets between the bound and unbound states is as small as 0.9A, suggesting that such pockets constitute features of the populated native state that may be used in docking.

Models, Molecular↗

PubMed2004

Is allostery an intrinsic property of all dynamic proteins?

Allostery involves coupling of conformational changes between two widely separated binding sites. The common view holds that allosteric proteins are symmetric oligomers, with each subunit existing in "at least" two conformational states with a different affinity for ligands. Recent observations such as the allosteric behavior of myoglobin, a classical example of a nonallosteric protein, call into question the existing allosteric dogma. Here we argue that all (nonfibrous) proteins are potentially allosteric. Allostery is a consequence of re-distributions of protein conformational ensembles. In a nonallosteric protein, the binding site shape may not show a concerted second-site change and enzyme kinetics may not reflect an allosteric transition. Nevertheless, appropriate ligands, point mutations, or external conditions may facilitate a population shift, leading a presumably nonallosteric protein to behave allosterically. In principle, practically any potential drug binding to the protein surface can alter the conformational redistribution. The question is its effectiveness in the redistribution of the ensemble, affecting the protein binding sites and its function. Here, we review experimental observations validating this view of protein allostery.

Allosteric Regulation↗