Liling Zhao

Novel Strategies for Drug Discovery Based on Intrinsically Disordered Proteins (IDPs)

Intrinsically disordered proteins (IDPs) are proteins that usually do not adopt well-defined native structures when isolated in solution under physiological conditions. Numerous IDPs have close relationships with human diseases such as tumor, Parkinson disease, Alzheimer disease, diabetes, and so on. These disease-associated IDPs commonly play principal roles in the disease-associated protein-protein interaction networks. Most of them in the disease datasets have more interactants and hence the size of the disease-associated IDPs interaction network is simultaneously increased. For example, the tumor suppressor protein p53 is an intrinsically disordered protein and also a hub protein in the p53 interaction network; α-synuclein, an intrinsically disordered protein involved in Parkinson diseases, is also a hub of the protein network. The disease-associated IDPs may provide potential targets for drugs modulating protein-protein interaction networks. Therefore, novel strategies for drug discovery based on IDPs are in the ascendant. It is dependent on the features of IDPs to develop the novel strategies. It is found out that IDPs have unique structural features such as high flexibility and random coil-like conformations which enable them to participate in both the “one to many” and “many to one” interaction. Accordingly, in order to promote novel strategies for drug discovery, it is essential that more and more features of IDPs are revealed by experimental and computing methods.

Effects of Different Force Fields and Temperatures on the Structural Character of Abeta (12–28) Peptide in Aqueous Solution

The aim of this work is to investigate the effects of different force fields and temperatures on the structural character of Aβ (12–28) peptide in aqueous solution. Moreover, the structural character of Aβ (12–28) peptide is compared with other amyloid peptides (such as H1 and α-syn12 peptide). The two independent temperature replica exchange molecular dynamics (T-REMD) simulations were completed by using two different models (OPLS-AA/TIP4P and GROMOS 43A1/SPC). We compared the models by analyzing the distributions of backbone dihedral angles, the secondary structure propensity, the free energy surface and the formation of β-hairpin. The results show that the mostly populated conformation state is random coil for both models. The population of β-hairpin is below 8 percent for both models. However, the peptide modeled by GROMOS 43A1 form β-hairpin with turn located at residues F19-E22, while the peptide modeled by OPLS-AA form β-hairpin with turn located at residues L17-F20.

Computational Studies of a Mechanism for Binding and Drug Resistance in the Wild Type and Four Mutations of HIV-1 Protease with a GRL-0519 Inhibitor.

Drug resistance of mutations in HIV-1 protease (PR) is the most severe challenge to the long-term efficacy of HIV-1 PR inhibitor in highly active antiretroviral therapy. To elucidate the molecular mechanism of drug resistance associated with mutations (D30N, I50V, I54M, and V82A) and inhibitor (GRL-0519) complexes, we have performed five molecular dynamics (MD) simulations and calculated the binding free energies using the molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) method. The ranking of calculated binding free energies is in accordance with the experimental data. The free energy spectra of each residue and inhibitor interaction for all complexes show a similar binding model. Analysis based on the MD trajectories and contribution of each residues show that groups R2 and R3 mainly contribute van der Waals energies, while groups R1 and R4 contribute electrostatic interaction by hydrogen bonds. The drug resistance of D30N can be attributed to the decline in binding affinity of residues 28 and 29. The size of Val50 is smaller than Ile50 causes the residue to move, especially in chain A. The stable hydrophobic core, including the side chain of Ile54 in the wild type (WT) complex, became unstable in I54M because the side chain of Met54 is flexible with two alternative conformations. The binding affinity of Ala82 in V82A decreases relative to Val82 in WT. The present study could provide important guidance for the design of a potent new drug resisting the mutation inhibitors.

Bias-Exchange Metadynamics Simulation of Membrane Permeation of 20 Amino Acids

Thermodynamics of the permeation of amino acids from water to lipid bilayers is an important first step for understanding the mechanism of cell-permeating peptides and the thermodynamics of membrane protein structure and stability. In this work, we employed bias-exchange metadynamics simulations to simulate the membrane permeation of all 20 amino acids from water to the center of a dipalmitoylphosphatidylcholine (DPPC) membrane (consists of 256 lipids) by using both directional and torsion angles for conformational sampling. The overall accuracy for the free energy profiles obtained is supported by significant correlation coefficients (correlation coefficient at 0.5–0.6) between our results and previous experimental or computational studies. The free energy profiles indicated that (1) polar amino acids have larger free energy barriers than nonpolar amino acids; (2) negatively charged amino acids are the most difficult to enter into the membrane; and (3) conformational transitions for many amino acids during membrane crossing is the key for reduced free energy barriers. These results represent the first set of simulated free energy profiles of membrane crossing for all 20 amino acids.

Molecular Dynamics Simulations of Human Antimicrobial Peptide LL-37 in Model POPC and POPG Lipid Bilayers

Cathelicidins are a large family of cationic antimicrobial peptides (AMPs) found in mammals with broad spectrum antimicrobial activity. LL-37 is the sole amphipathic α-helical AMP from human Cathelicidins family. In addition to its bactericidal capability, LL-37 has antiviral, anti-tumor, and immunoregulatory activity. Despite many experimental studies, its molecular mechanism of action is not yet fully understood. Here, we performed three independent molecular dynamics simulations (600 ns or more) of a LL-37 peptide in the presence of 256 lipid bilayers with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) mimicking bacterial and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) mimicking mammalian membranes. We found that LL-37 can be quickly absorbed onto the POPG bilayer without loss of its helical conformation in the core region and with the helix lying in parallel to the bilayer. The POPG bilayer was deformed. In contrast, LL-37 is slower in reaching the POPC surface and loss much of its helical conformation during the interaction with the bilayer. LL-37 only partially entered the POPC bilayer without significant deformation of the membrane. The observed difference for different bilayers is largely due to the fact that LL-37 is positively charged, POPG is negatively charged, and POPC is neutral. Our simulation results demonstrated the initial stage of disruption of the bacterial membrane by LL-37 in atomic details. Comparison to experimental results on LL-37 and simulation studies in other systems was made.

Fast-Folding Pathways of the Thrombin-Binding Aptamer G-Quadruplex Revealed by a Markov State Model

G-quadruplex structures participate in many important cellular processes. For a better understanding of their functions, knowledge of the mechanism by which they fold into the functional native structures is necessary. In this work, we studied the folding process of the thrombin-binding aptamer G-quadruplex. Enabled by a computational paradigm that couples an advanced sampling method and a Markov state model, four folding intermediates were identified, including an antiparallel G-hairpin, two G-triplex structures, and a double-hairpin conformation. Likewise, a misfolded structure with a nonnative distribution of syn/anti guanines was also observed. Based on these states, a transition path analysis revealed three fast-folding pathways, along which the thrombin-binding aptamer would fold to the native state directly, with no evidence of potential nonnative competing conformations. The results also showed that the TGT-loop plays an important role in the folding process. The findings of this research may provide general insight about the folding of other G-quadruplex structures.

Effect of pH on the Aggregation of α-syn12 Dimer in Explicit Water by Replica-Exchange Molecular Dynamics Simulation

The dimeric structure of the N-terminal 12 residues drives the interaction of α-synuclein protein with membranes. Moreover, experimental studies indicated that the aggregation of α-synuclein is faster at low pH than neutral pH. Nevertheless, the effects of different pH on the structural characteristics of the α-syn12 dimer remain poorly understood. We performed 500 ns temperature replica exchange molecular dynamics (T-REMD) simulations of two α-syn12 peptides in explicit solvent. The free energy surfaces contain ten highly populated regions at physiological pH, while there are only three highly populated regions contained at acidic pH. The anti-parallel β-sheet conformations were found as the lowest free energy state. Additionally, these states are nearly flat with a very small barrier which indicates that these states can easily transit between themselves. The dimer undergoes a disorder to order transition from physiological pH to acidic pH and the α-syn12 dimer at acidic pH involves a faster dimerization process. Further, the Lys6–Asp2 contact may prevent the dimerization.

The Effect of C-Terminal Helix on the Stability of FF Domain Studied by Molecular Dynamics Simulation

To investigate the effect of C-terminal helix on the stability of the FF domain, we studied the native domain FF3-71 from human HYPA/FBP11 and the truncated version FF3-60 with C-terminal helix being deleted by molecular dynamics simulations with GROMACS package and GROMOS 43A1 force field. The results indicated that the structures of truncated version FF3-60 were evident different from those of native partner FF3-71. Compared with FF3-71, the FF3-60 lost some native contacts and exhibited some similar structural characters to those of intermediate state. The C-terminal helix played a major role in stabilizing the FF3-71 domain. To a certain degree, the FF domain had a tendency to form an intermediate state without the C-terminal helix. In our knowledge, this was the first study to examine the role of C-terminal helix of FF domain in detail by molecular dynamics simulations, which was useful to understand the three-state folding mechanism of the small FF domain.