RESUMEN
Coarse-grained models designed for intrinsically disordered proteins and regions (IDP/Rs) usually omit some bonded potentials (e.g., angular and dihedral potentials) as a conventional strategy to enhance backbone flexibility. However, a notable drawback of this approach is the generation of inaccurate backbone conformations. Here, we addressed this problem by introducing residue-specific angular, refined dihedral, and correction map (CMAP) potentials, derived based on the statistics from a customized coil database. These bonded potentials were integrated into the existing Mpipi model, resulting in a new model, denoted as the "Mpipi+" model. Results show that the Mpipi+ model can improve backbone conformations. More importantly, it can markedly improve the secondary structure propensity (SSP) based on the experimental chemical shift and, consequently, succeed in capturing transient secondary structures. Moreover, the Mpipi+ model preserves the liquid-liquid phase separation (LLPS) propensities of IDPs.
Asunto(s)
Proteínas Intrínsecamente Desordenadas , Proteínas Intrínsecamente Desordenadas/química , Estructura Secundaria de Proteína , Modelos Moleculares , Conformación ProteicaRESUMEN
Coronins play critical roles in actin network formation. The diverse functions of coronins are regulated by the structured N-terminal ß propeller and the C-terminal coiled coil (CC). However, less is known about a middle "unique region" (UR), which is an intrinsically disordered region (IDR). The UR/IDR is an evolutionarily conserved signature in the coronin family. By integrating biochemical and cell biology experiments, coarse-grained simulations, and protein engineering, we find that the IDR optimizes the biochemical activities of coronins in vivo and in vitro. The budding yeast coronin IDR plays essential roles in regulating Crn1 activity by fine-tuning CC oligomerization and maintaining Crn1 as a tetramer. The IDR-guided optimization of Crn1 oligomerization is critical for F-actin cross-linking and regulation of Arp2/3-mediated actin polymerization. The final oligomerization status and homogeneity of Crn1 are contributed by three examined factors: helix packing, the energy landscape of the CC, and the length and molecular grammar of the IDR.
Asunto(s)
Citoesqueleto de Actina , Actinas , Proteínas Intrínsecamente Desordenadas , Citoesqueleto de Actina/metabolismo , Actinas/metabolismo , Polimerizacion , Proteínas Intrínsecamente Desordenadas/metabolismo , Proteínas Intrínsecamente Desordenadas/fisiología , Proteínas de Microfilamentos/metabolismo , Proteínas de Microfilamentos/fisiología , Saccharomyces cerevisiae/genética , Humanos , Proteínas de Saccharomyces cerevisiae/metabolismo , Proteínas de Saccharomyces cerevisiae/fisiologíaRESUMEN
Actin nucleation is achieved by collaborative teamwork of actin nucleator factors (NFs) and nucleation-promoting factors (NPFs) into functional protein complexes. Selective inter- and intramolecular interactions between the nucleation complex constituents enable diverse modes of complex assembly in initiating actin polymerization on demand. Budding yeast has two formins, Bni1 and Bnr1, which are teamed up with different NPFs. However, the selective pairing between formin NFs and NPFs into the nucleation core for actin polymerization is not completely understood. By examining the functions and interactions of NPFs and NFs via biochemistry, genetics, and mathematical modeling approaches, we found that two NPFs, Aip5 and Bud6, showed joint teamwork effort with Bni1 and Bnr1, respectively, by interacting with the C-terminal intrinsically disordered region (IDR) of formin, in which two NPFs work together to promote formin-mediated actin nucleation. Although the C-terminal IDRs of Bni1 and Bnr1 are distinct in length, each formin IDR orchestrates the recruitment of Bud6 and Aip5 cooperatively by different positioning strategies to form a functional complex. Our study demonstrated the dynamic assembly of the actin nucleation complex by recruiting multiple partners in budding yeast, which may be a general feature for effective actin nucleation by formins.
Asunto(s)
Proteínas de Microfilamentos/metabolismo , Proteínas de Saccharomyces cerevisiae/metabolismo , Actinas/metabolismo , Secuencia de Aminoácidos , Proteínas del Citoesqueleto/metabolismo , Forminas , Proteínas de Microfilamentos/genética , Saccharomyces cerevisiae/metabolismo , Proteínas de Saccharomyces cerevisiae/genéticaRESUMEN
Candida albicans is a dimorphic fungus that converts from a yeast form to a hyphae form during infection. This switch requires the formation of actin cable to coordinate polarized cell growth. It's known that nucleation of this cable requires a multiprotein complex localized at the tip called the polarisome, but the mechanisms underpinning this process were unclear. Here, we found that C. albicans Aip5, a homolog of polarisome component ScAip5 in Saccharomyces cerevisiae that nucleates actin polymerization and synergizes with the formin ScBni1, regulates actin assembly and hyphae growth synergistically with other polarisome proteins Bni1, Bud6, and Spa2. The C terminus of Aip5 binds directly to G-actin, Bni1, and the C-terminal of Bud6, which form the core of the nucleation complex to polymerize F-actin. Based on insights from structural biology and molecular dynamic simulations, we propose a possible complex conformation of the actin nucleation core, which provides cooperative positioning and supports the synergistic actin nucleation activity of a tri-protein complex Bni1-Bud6-Aip5. Together with known interactions of Bni1 with Bud6 and Aip5 in S. cerevisiae, our findings unravel molecular mechanisms of C. albicans by which the tri-protein complex coordinates the actin nucleation in actin cable assembly and hyphal growth, which is likely a conserved mechanism in different filamentous fungi and yeast.
Asunto(s)
Actinas/metabolismo , Candida albicans/crecimiento & desarrollo , Candida albicans/metabolismo , Proteínas Fúngicas/metabolismo , PolimerizacionRESUMEN
We explored the stability of the dengue virus envelope (E) protein dimer since it is widely assumed that the E protein dimer is stabilized by drug ligands or antibodies in an acidic environment, neutralizing the virus's ability to fuse with human cells. During this process, a large conformational change of the E protein dimer is required. We performed Molecular Dynamics simulations to mimic the conformational change and stability of the dimer in neutral and acidic conditions with the well-tempered metadynamics method. Furthermore, as a few neutralizing antibodies discovered from dengue patients were reported, we used the same simulation method to examine the influence of a selected antibody on the dimer stability in both neutral and acidic conditions. We also investigated the antibody's influence on a point-mutated E protein that had been reported to interrupt the protein-antibody interaction and result in more than 95% loss of the antibody's binding ability. Our simulation results are highly consistent with the experimental conclusion that binding of the antibody to the E protein dimer neutralizes the virus, especially in a low pH condition, while the mutation of W101A or N153A significantly reduces the antibody's ability in stabilizing the E protein dimer. We demonstrate that well-tempered metadynamics can be used to accurately explore the antibody's interaction on large protein complexes such as the E protein dimer, and the computational approach in this work is promising in future antibody development.
Asunto(s)
Virus del Dengue/química , Multimerización de Proteína , Proteínas del Envoltorio Viral/química , Envoltura Viral/química , Anticuerpos Monoclonales/química , Anticuerpos Monoclonales/metabolismo , Anticuerpos Neutralizantes/química , Anticuerpos Neutralizantes/metabolismo , Anticuerpos Antivirales/química , Anticuerpos Antivirales/metabolismo , Humanos , Concentración de Iones de Hidrógeno , Ligandos , Simulación de Dinámica Molecular , Mutación Puntual , Estabilidad Proteica , Proteínas del Envoltorio Viral/genética , Proteínas del Envoltorio Viral/metabolismoRESUMEN
Studies on the Bin-Amphiphysin-Rvs (BAR) domain have advanced a fundamental understanding of how proteins deform membrane. We previously showed that a BAR domain in tandem with a Pleckstrin Homology (PH domain) underlies the assembly of ACAP1 (Arfgap with Coil-coil, Ankryin repeat, and PH domain I) into an unusual lattice structure that also uncovers a new paradigm for how a BAR protein deforms membrane. Here, we initially pursued computation-based refinement of the ACAP1 lattice to identify its critical protein contacts. Simulation studies then revealed how ACAP1, which dimerizes into a symmetrical structure in solution, is recruited asymmetrically to the membrane through dynamic behavior. We also pursued electron microscopy (EM)-based structural studies, which shed further insight into the dynamic nature of the ACAP1 lattice assembly. As ACAP1 is an unconventional BAR protein, our findings broaden the understanding of the mechanistic spectrum by which proteins assemble into higher-ordered structures to achieve membrane deformation.
Asunto(s)
Proteínas Activadoras de GTPasa/metabolismo , Proteínas de la Membrana/metabolismo , Membrana Celular/metabolismo , Dimerización , Proteínas Activadoras de GTPasa/química , Humanos , Dominios Homólogos a Pleckstrina , Conformación ProteicaRESUMEN
Zika virus (ZIKV) has become a great public health emergency. Its non-structural protein 3 (NS3) is a key enzyme in viral replication and has been considered as a potential therapeutic target. A conformational characterization of ZIKV NS3 is critical for a comprehensive understanding of its molecular interactions and functions. However, the high conformational flexibility of solution NS3 obstacles the structural characterization of NS3 solely from the experimental observable that averages over its heterogeneous conformations. Here, we employed replica exchange with solute tempering (REST) method to simulate the di-domain protein ZIKV NS3. Three independent MD simulations identified a conserved conformational ensemble of NS3, consisting of a major conformational state and several minor states from compact to loose conformations. The major state agrees well with the scattering profile from small-angle X-ray scattering (SAXS) experiments. Moreover, the simulated ensemble is supported by a direct data-fitting result that requires both short- and long-range structural contacts to recover the experimental data. We discussed the interplay between simulation and experiment in ensemble construction of flexible biomolecules and shed light on the physically derived conformational ensembles.
Asunto(s)
Simulación de Dinámica Molecular , Péptido Hidrolasas/química , Dispersión del Ángulo Pequeño , Proteínas Virales/química , Difracción de Rayos X , Dominios Proteicos , Serina EndopeptidasasRESUMEN
Structure characterization of intrinsically disordered proteins (IDPs) remains a key obstacle in understanding their functional mechanisms. Due to the highly dynamic feature of IDPs, structure ensembles instead of static unique structures are often derived from experimental data. Several state-of-the-art computational methods have been developed to select an optimal ensemble from a pregenerated structure pool, but they suffer from low efficiency for large IDPs. Here we present a matching pursuit genetic algorithm (MPGA) for structure ensemble determination, which takes advantages from both matching pursuit (MP) to reduce the search space and genetic algorithm (GA) to reduce the restriction on constraint types. The MPGA method is validated using a reference ensemble with predefined structures. In comparison with the conventional GA, MPGA takes much less computational time for large IDPs. The utility of the method is demonstrated by application to structure ensemble determination of a mechanosensing protein domain with 306 amino acids. The structure ensemble determined reveals that the N-terminal region 1-240 is more compact than the C-terminal region 240-306. The unique structural feature explains why only a small portion of YXXP tyrosine residues can be phosphorylated easily by kinases in the absence of extension force and why the phosphorylation is force-dependent.
Asunto(s)
Algoritmos , Proteínas Intrínsecamente Desordenadas/química , Proteínas Intrínsecamente Desordenadas/genética , Modelos Moleculares , Conformación ProteicaRESUMEN
Double-stranded RNA (dsRNA) structures form triplexes and RNA-protein complexes through binding to single-stranded RNA (ssRNA) regions and proteins, respectively, for diverse biological functions. Hence, targeting dsRNAs through major-groove triplex formation is a promising strategy for the development of chemical probes and potential therapeutics. Short (e.g., 6-10 mer) chemically-modified Peptide Nucleic Acids (PNAs) have been developed that bind to dsRNAs sequence specifically at physiological conditions. For example, a PNA incorporating a modified base thio-pseudoisocytosine (L) has an enhanced recognition of a G-C pair in an RNA duplex through major-groove L·G-C base triple formation at physiological pH, with reduced pH dependence as observed for C+·G-C base triple formation. Currently, an unmodified T base is often incorporated into PNAs to recognize a Watson-Crick A-U pair through major-groove T·A-U base triple formation. A substitution of the 5-methyl group in T by hydrogen and halogen atoms (F, Cl, Br, and I) causes a decrease of the pKa of N3 nitrogen atom, which may result in improved hydrogen bonding in addition to enhanced base stacking interactions. Here, we synthesized a series of PNAs incorporating uracil and halouracils, followed by binding studies by non-denaturing polyacrylamide gel electrophoresis, circular dichroism, and thermal melting. Our results suggest that replacing T with uracil and halouracils may enhance the recognition of an A-U pair by PNA·RNA2 triplex formation in a sequence-dependent manner, underscoring the importance of local stacking interactions. Incorporating bromouracils and chlorouracils into a PNA results in a significantly reduced pH dependence of triplex formation even for PNAs containing C bases, likely due to an upshift of the apparent pKa of N3 atoms of C bases. Thus, halogenation and other chemical modifications may be utilized to enhance hydrogen bonding of the adjacent base triples and thus triplex formation. Furthermore, our experimental and computational modelling data suggest that PNA·RNA2 triplexes may be stabilized by incorporating a BrUL step but not an LBrU step, in dsRNA-binding PNAs.
Asunto(s)
Emparejamiento Base/genética , Halógenos/química , Conformación de Ácido Nucleico , Ácidos Nucleicos de Péptidos/química , ARN Bicatenario/síntesis química , Uracilo/análogos & derivados , Uracilo/química , Bromouracilo/química , Línea Celular Tumoral , Biología Computacional/métodos , Simulación por Computador , Halogenación , Células HeLa , Humanos , Enlace de Hidrógeno , Secuencias Invertidas Repetidas/genética , MicroARNs/genética , Proteínas de Unión al ARN/químicaRESUMEN
The structural variations of multidomain proteins with flexible parts mediate many biological processes, and a structure ensemble can be determined by selecting a weighted combination of representative structures from a simulated structure pool, producing the best fit to experimental constraints such as interatomic distance. In this study, a hybrid structure-based and physics-based atomistic force field with an efficient sampling strategy is adopted to simulate a model di-domain protein against experimental paramagnetic relaxation enhancement (PRE) data that correspond to distance constraints. The molecular dynamics simulations produce a wide range of conformations depicted on a protein energy landscape. Subsequently, a conformational ensemble recovered with low-energy structures and the minimum-size restraint is identified in good agreement with experimental PRE rates, and the result is also supported by chemical shift perturbations and small-angle X-ray scattering data. It is illustrated that the regularizations of energy and ensemble-size prevent an arbitrary interpretation of protein conformations. Moreover, energy is found to serve as a critical control to refine the structure pool and prevent data overfitting, because the absence of energy regularization exposes ensemble construction to the noise from high-energy structures and causes a more ambiguous representation of protein conformations. Finally, we perform structure-ensemble optimizations with a topology-based structure pool, to enhance the understanding on the ensemble results from different sources of pool candidates.
Asunto(s)
Simulación de Dinámica Molecular , Proteínas de Unión a Poli(A)/química , Proteínas de Saccharomyces cerevisiae/química , Aminoácidos/química , Sitios de Unión , Espectroscopía de Resonancia por Spin del Electrón , Unión Proteica , Dominios Proteicos , Estructura Secundaria de Proteína , Saccharomyces cerevisiae , Relación Estructura-Actividad , TermodinámicaRESUMEN
Integrative structure modeling is an emerging method for structural determination of protein-protein complexes that are challenging for conventional structural techniques. Here, we provide a practical protocol for implementing our integrated iSPOT platform by integrating three different biophysical techniques: small-angle X-ray scattering (SAXS), hydroxyl radical footprinting, and computational docking simulations. Specifically, individual techniques are described from experimental and/or computational perspectives, and complementary structural information from these different techniques are integrated for accurate characterization of the structures of large protein-protein complexes.
Asunto(s)
Espectrometría de Masas/métodos , Impresión Molecular/métodos , Complejos Multiproteicos/ultraestructura , Proteínas/ultraestructura , Dispersión del Ángulo Pequeño , Humanos , Radical Hidroxilo/química , Simulación del Acoplamiento Molecular , Complejos Multiproteicos/química , Unión Proteica , Conformación Proteica , Proteínas/química , Difracción de Rayos XRESUMEN
An ensemble-modeling scheme incorporating coarse-grained simulations with experimental small-angle X-ray scattering (SAXS) data is applied to dengue virus 2 (DENV2) nonstructural protein 5 (NS5). NS5 serves a key role in viral replication through its two domains that are connected by a 10-residue polypeptide segment. A set of representative structures is generated from a simulated structure pool using SAXS data fitting by the non-negativity least squares (NNLS) or standard ensemble optimization method (EOM) based on a genetic algorithm (GA). It is found that a proper low-energy threshold of the structure pool is necessary to produce a conformational ensemble of two representative structures by both NNLS and GA that agrees well with the experimental SAXS profile. The stability of the constructed ensemble is validated also by molecular dynamics simulations with an all-atom force field. The constructed ensemble successfully revealed the domain-domain orientation and domain-contacting interface of DENV2 NS5. Using experimental data fitting and additional investigations with synthesized data, it is found that energy restraint on the conformational pool is necessary to avoid overinterpretation of experimental data by spurious conformational representations.
Asunto(s)
Virus del Dengue/metabolismo , Proteínas no Estructurales Virales/química , Modelos Químicos , Simulación de Dinámica Molecular , Conformación Proteica , Dispersión del Ángulo Pequeño , Difracción de Rayos XRESUMEN
ACAP1 proteins were previously reported to specifically bind PIP2-containing cell membranes and form well-structured protein lattices in order to conduct membrane tubulation. We carried out molecular dynamics simulations to characterize orientation of the PH domains with respect to the BAR domains inside the protein dimer. Followed by molecular dynamics simulations, we present a comprehensive orientation analysis of PH domain under different states including unbound and bound with lipids. We have examined two binding pockets on the PH domain and present PMF profiles of the two pockets to account for their preference to PIP2 lipids. Combining orientation analysis and studies of binding pockets, our simulations results reveal valuable molecular basis for protein-lipid interactions of ACAP1 proteins during membrane remodeling process.
Asunto(s)
Proteínas Activadoras de GTPasa/química , Membrana Dobles de Lípidos/química , Simulación de Dinámica Molecular , Humanos , Lípidos/química , Dominios Homólogos a PleckstrinaRESUMEN
Elastic network models (ENMs) based on simple harmonic potential energy function have been proven over the last decade to be reliable computational models for understanding the intrinsic dynamics of biomacromolecules. In the original ENMs, the spring constants for different contact pairs are assumed to be identical, while there are a number of recent developments to determine non-uniform spring constants from atomistic force fields or experimental information. In particular, the fluctuation matching approaches in coarse-grained modeling can be applied to build more realistic heterogeneous ENMs, using information from an atomistic force field or experimental B-factors. The same type of approaches is further implemented to parameterize heterogeneous structure-based models, which can be considered as a natural extension of ENMs in terms of the potential energy function. In this review, we give an overview of different fluctuation matching methods adopted for ENMs and structure-based models, including an improved formulation and algorithm based on the relative entropy scheme.
Asunto(s)
Elasticidad , Sustancias Macromoleculares/química , Modelos Moleculares , Sustancias Macromoleculares/metabolismoRESUMEN
Minus-one ribosomal frameshifting is a translational recoding mechanism widely utilized by many RNA viruses to generate accurate ratios of structural and catalytic proteins. An RNA pseudoknot structure located in the overlapping region of the gag and pro genes of Simian Retrovirus type 1 (SRV-1) stimulates frameshifting. However, the experimental characterization of SRV-1 pseudoknot (un)folding dynamics and the effect of the base triple formation is lacking. Here, we report the results of our single-molecule nanomanipulation using optical tweezers and theoretical simulation by steered molecular dynamics. Our results directly reveal that the energetic coupling between loop 2 and stem 1 via minor-groove base triple formation enhances the mechanical stability. The terminal base pair in stem 1 (directly in contact with a translating ribosome at the slippery site) also affects the mechanical stability of the pseudoknot. The -1 frameshifting efficiency is positively correlated with the cooperative one-step unfolding force and inversely correlated with the one-step mechanical unfolding rate at zero force. A significantly improved correlation was observed between -1 frameshifting efficiency and unfolding rate at forces of 15-35 pN, consistent with the fact that the ribosome is a force-generating molecular motor with helicase activity. No correlation was observed between thermal stability and -1 frameshifting efficiency.
Asunto(s)
Sistema de Lectura Ribosómico , Productos del Gen gag/química , ARN Viral/genética , Retrovirus de los Simios , Ribosomas/química , Electroforesis en Gel de Poliacrilamida , Cinética , Espectroscopía de Resonancia Magnética , Simulación de Dinámica Molecular , Mutación , Conformación de Ácido Nucleico , Oligonucleótidos/química , Pinzas Ópticas , Desnaturalización Proteica , Pliegue de Proteína , ARN Helicasas/química , ARN Mensajero/metabolismo , TermodinámicaRESUMEN
Structure modelling via small-angle X-ray scattering (SAXS) data generally requires intensive computations of scattering intensity from any given biomolecular structure, where the accurate evaluation of SAXS profiles using coarse-grained (CG) methods is vital to improve computational efficiency. To date, most CG SAXS computing methods have been based on a single-bead-per-residue approximation but have neglected structural correlations between amino acids. To improve the accuracy of scattering calculations, accurate CG form factors of amino acids are now derived using a rigorous optimization strategy, termed electron-density matching (EDM), to best fit electron-density distributions of protein structures. This EDM method is compared with and tested against other CG SAXS computing methods, and the resulting CG SAXS profiles from EDM agree better with all-atom theoretical SAXS data. By including the protein hydration shell represented by explicit CG water molecules and the correction of protein excluded volume, the developed CG form factors also reproduce the selected experimental SAXS profiles with very small deviations. Taken together, these EDM-derived CG form factors present an accurate and efficient computational approach for SAXS computing, especially when higher molecular details (represented by the q range of the SAXS data) become necessary for effective structure modelling.
RESUMEN
The single-mutation of genes associated with Alzheimer's disease (AD) increases the production of Aß peptides. An elevated concentration of Aß peptides is prone to aggregation into oligomers and further deposition as plaque. Aß plaques and neurofibrillary tangles are two hallmarks of AD. In this review, we provide a broad overview of the diverses sources that could lead to AD, which include genetic origins, Aß peptides and tau protein. We shall discuss on tau protein and tau accumulation, which result in neurofibrillary tangles. We detail the mechanisms of Aß aggregation, fibril formation and its polymorphism. We then show the possible links between Aß and tau pathology. Furthermore, we summarize the structural data of Aß and its precursor protein obtained via Nuclear Magnetic Resonance (NMR) or X-ray crystallography. At the end, we go through the C-terminal and N-terminal truncated Aß variants. We wish to draw reader's attention to two predominant and toxic Aß species, namely Aß4-42 and pyroglutamate amyloid-beta peptides, which have been neglected for more than a decade and may be crucial in Aß pathogenesis due to their dominant presence in the AD brain.
Asunto(s)
Enfermedad de Alzheimer/patología , Enfermedad de Alzheimer/tratamiento farmacológico , Enfermedad de Alzheimer/genética , Péptidos beta-Amiloides/química , Péptidos beta-Amiloides/genética , Péptidos beta-Amiloides/metabolismo , Inhibidores de la Colinesterasa/uso terapéutico , Humanos , Nootrópicos/uso terapéutico , Isoformas de Proteínas/química , Isoformas de Proteínas/genética , Isoformas de Proteínas/metabolismo , Ácido Pirrolidona Carboxílico/química , Proteínas tau/química , Proteínas tau/metabolismoRESUMEN
Elastic network models (ENM) are based on the idea that the geometry of a protein structure provides enough information for computing its fluctuations around its equilibrium conformation. This geometry is represented as an elastic network (EN) that is, a network of links between residues. A spring is associated with each of these links. The normal modes of the protein are then identified with the normal modes of the corresponding network of springs. Standard approaches for generating ENs rely on a cutoff distance. There is no consensus on how to choose this cutoff. In this work, we propose instead to filter the set of all residue pairs in a protein using the concept of alpha shapes. The main alpha shape we considered is based on the Delaunay triangulation of the Cα positions; we referred to the corresponding EN as EN(∞). We have shown that heterogeneous anisotropic network models, called αHANMs, that are based on EN(∞) reproduce experimental B-factors very well, with correlation coefficients above 0.99 and root-mean-square deviations below 0.1 Å(2) for a large set of high resolution protein structures. The construction of EN(∞) is simple to implement and may be used automatically for generating ENs for all types of ENMs.
Asunto(s)
Proteínas/química , Algoritmos , Anisotropía , Simulación por Computador , Modelos Químicos , Modelos Moleculares , Conformación ProteicaRESUMEN
An iterative coarse-graining method is developed for systematically converting an atomistic force field to a model at lower resolution that is able to accurately reproduce the distribution functions defined in the coarse-grained potential. The method starts from the multiscale coarse-graining (MS-CG) approach, and it iteratively refines the distribution functions using repeated applications of the MS-CG algorithm. It is justified on the basis of the force matching normal equation, which can be considered a discrete form of the Yvon-Born-Green equation in liquid state theory. Numerical results for molecular systems involving pairwise nonbonded and three-body bonded interactions are obtained, and comparison with other approaches in literature is provided.
Asunto(s)
Algoritmos , Hexanos/química , Metanol/química , Simulación de Dinámica Molecular , Agua/químicaRESUMEN
Internal molecular forces can guide chemical reactions, yet are not straightforwardly accessible within a quantum mechanical description of the reacting molecules. Here, we present a force-matching force distribution analysis (FM-FDA) to analyze internal forces in molecules. We simulated the ring opening of trans-3,4-dimethylcyclobutene (tDCB) with on-the-fly semiempirical molecular dynamics. The self-consistent density functional tight binding (SCC-DFTB) method accurately described the force-dependent ring-opening kinetics of tDCB, showing quantitative agreement with both experimental and computational data at higher levels. Mechanical force was applied in two different ways, namely, externally by a constant pulling force and internally by embedding tDCB within a strained macrocycle-containing stiff stilbene. We analyzed the distribution of tDCB internal forces in the two different cases by FM-FDA and found that external force gave rise to a symmetric force distribution in the cyclobutene ring, which also scaled linearly with the external force, indicating that the force distribution was uniquely determined by the symmetric architecture of tDCB. In contrast, internal forces due to stiff stilbene resulted in an asymmetric force distribution within tDCB, which indicated a different geometry of force application and supported the important role of linkers in the mechanochemical reactivity of tDCB. In addition, three coordinates were identified through which the distributed forces contributed most to rate acceleration. These coordinates are mostly parallel to the coordinate connecting the two CH3 termini of tDCB. Our results confirm previous observations that the linker outside of the reactive moiety, such as a stretched polymer or a macrocycle, affects its mechanochemical reactivity. We expect FM-FDA to be of wide use to understand and quantitatively predict mechanochemical reactivity, including the challenging cases of systems within strained macrocycles.