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.

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.

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.

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.

Experiments have shown that cholesterol influences the membrane permeability of small molecules, amino acids, and cell-penetrating peptides. However, their exact translocation mechanisms under the influence of cholesterol remain poorly understood. Given the practical importance of cell-penetrating peptides and the existence of varied cholesterol contents in different cell types, it is necessary to examine the permeation of amino acids in cholesterol-containing membranes at atomic level of details. Here, bias-exchange metadynamics simulations were employed to investigate the molecular mechanism of the membrane permeation of two amino acids Arg and Trp important for cell-penetrating peptides in the presence of different concentrations of cholesterol. We found that the free energy barrier of Arg+ (the protonated form) permeation increased linearly as the cholesterol concentration increased, whereas the barrier of Trp permeation had a rapid increase from 0 mol. % to 20 mol. % cholesterol-containing membranes and nearly unchanged from 20 mol. % to 40 mol. % cholesterol-containing membranes. Arg0 becomes slightly more stable than Arg+ at the center of the dipalmitoylphosphatidylcholine (DPPC) membrane with 40 mol. % cholesterol concentrations. As a result, Arg+ has a similar permeability as Trp at 0 mol. % and 20 mol. % cholesterol, but a significantly lower permeability than Trp at 40 mol. % cholesterol. This difference is caused by the gradual reduction of water defects for Arg+ as the cholesterol concentration increases but lack of water defects for Trp in cholesterol-containing membranes. Strong but different orientation dependence between Arg+ and Trp permeations is observed. These results provide an improved microscopic understanding of amino-acid permeation through cholesterol-containing DPPC membrane systems.

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.

Background: The studies on protein folding/unfolding indicate that the native state topology is an important determinant of protein folding mechanism. The folding/unfolding behaviors of proteins which have similar topologies have been studied under Cartesian space and the results indicate that some proteins share the similar folding/unfolding characters.

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.